Effect of NOM on copper sulfide nanoparticle growth, stability, and oxidative dissolution †

Metal sulfide nanoparticles have recently been discovered in natural environments subject to fluctuating redox conditions such as wetlands and river floodplains, which are often very rich in natural organic matter (NOM). Strong binding of chalcophile metals like copper by NOM is expected to interfere with sulfide precipitation and is likely to affect the growth of these nanoparticles. We therefore conducted experiments on the formation of copper sulfide nanoparticles (Cu x S NPs) in anoxic solutions with varying Cu (50, 500 μ mol L − 1 ) and sulfide (100, 1000 μ mol L − 1 ) concentrations in the absence and presence of Suwannee River fulvic acid (SRFA, 0, 5, 50 mg C L − 1 ). The size development of the Cu x S nanoparticles and their stability were tracked over 4 weeks using transmission electron microscopy (TEM) and size-exclusion chromatography coupled to ICP-MS (SEC-ICP-MS). Additional dissolution experiments were performed in closed containers in the absence and presence of O 2 and fulvic acid over several months providing insights into the oxidative dissolution behavior of Cu x S. Our results highlight the high colloidal stability of Cu x S nanoparticles in anoxic environments irrespective of the NOM concentration. Median particle diameters ranged between a few and a few tens of nanometers with larger particles forming in more concentrated suspensions. At low Cu and S concentrations, fulvic acid restricted particle growth by up to 25% compared to SRFA-free suspensions and metal sulfide clusters even smaller than 1 nm were detected. The Cu x S nanoparticles exhibited a remarkable chemical stability against oxidative dissolution and were only dissolved when both fulvic acid and O 2 were present. and river floodplains, known as typical hotspots for natural organic matter (NOM) production, have recently been recognized as sources of natural metal sulfide nanoparticles. These sulfide particles are expected to have a profound impact on trace metal cycling in these environments. We show that NOM hampers the growth of Cu x S nanoparticles and slows down their morphological development. Further, we demonstrate the importance of NOM in facilitating the oxidative dissolution of Cu x S nanoparticles. Our findings suggest an increased mobility of the Cu x S nanoparticles in NOM-rich environments, but also indicate enhanced oxidation of Cu x S and release of Cu 2+ ions in the presence of NOM upon aeration of anoxic waters.


Introduction
The release of trace metals from enriched parent rock material by weathering as well as from mining, ore processing, industrial activities and agri-or vinicultural practices can lead to trace metal accumulation in natural environments. Many floodplain and wetland areas around the world are contaminated with metals due to upland mine waste discharge or tailings dam failures. [1][2][3][4][5][6][7] Among those metals, copper (Cu) is an essential nutrient for all living organisms since it is required for the synthesis of numerous enzymes and proteins involved in redox-related biological functions. But its uptake is merely required in trace amounts. 8 Exposure to elevated concentrations of Cu is highly toxic to a wide range of organisms with fish, insects, algae, and microorganisms being among the most susceptible. 8,9 Therefore, excessive releases of Cu into natural environments pose a threat to ecosystem health, particularly when high biological diversity areas such as wetlands and marshes are affected. In this respect, a detailed understanding of trace metal behavior in these redox-variable systems is crucial for a risk assessment. Unlike oxic upland soils or permanently anoxic sediments, soils in these areas are subject to periodical flooding, leading to intense redoxdynamics. 10 In periods of prolonged waterlogging, anoxic conditions can establish within a few days as a consequence of microbial respiration quickly consuming molecular dioxygen. The absence of molecular dioxygen favors the use of alternative electron acceptors (e.g., Fe III , Mn III/IV , and S VI ) in biotic and abiotic reactions. [10][11][12] Bisulfides (HS − ) resulting from sulfate reduction have been shown to form sparingly soluble nanoscale precipitates with chalcophile trace metals (e.g., Cu x S, CdS) in the pore water of floodplain soils. [13][14][15] Likewise, other micro-and nanosized metal sulfides (mostly FeS, ZnS and HgS) have been observed in various natural environments, such as treated anaerobic wetlands, 16 anoxic lakes, 17 estuarine and marine waters and sediments, [18][19][20][21] biofilms, [22][23][24] and even in oxic surface waters, [25][26][27][28] thus proving their persistence in oxic environments over several months. Considering the general stability of metal sulfide nanoparticles under various redox conditions, natural copper sulfide nanoparticles (Cu x S NPs) may contribute to longrange transport of Cu in aquatic environments. However, indepth knowledge about key factors controlling their distribution over various environmental compartments is still limited.
Particle characteristics of Cu x S such as size, morphology, surface charge, and aggregation behavior are decisive for their colloidal transport, [29][30][31][32] control their reactivity towards the surrounding medium, 31,33,34 the availability to cells and organisms 35 and determine whether they can pass through a narrow soil pore. 32,[36][37][38] Key environmental factors such as ionic strength, pH, as well as copper, sulfide, and natural organic matter (NOM) concentrations control particle properties and may vary over a wide range depending on soil pore water composition. Large variations in the pore water concentrations of CuĲII) and SĲ-II) in floodplain soils are due to legacy contamination, regional parent rock composition, inflow of sulfur-rich waters, the flooding regime of the catchment or agricultural practices among others. In this context, the effect of elevated CuĲII) and SĲ-II) concentrations on the nucleation and growth of Cu x S nanoparticles has not been addressed yet.
The poor aeration of soils, being typical for wetlands, was shown to slow down organic matter decomposition and favor the accumulation of humic substances. 11,39 Since Cu exhibits a remarkable binding affinity towards NOM, [40][41][42] this affects Cu cycling in these systems. Furthermore, Cu can take part in redox interactions with NOM, especially under anoxic conditions, which are known to induce high affinity proton binding groups in humic substances. 39 However, inorganic HS − formed in anoxic pore waters is a much stronger ligand and competes for Cu binding so that Cu availability is limited by the formation of Cu-bearing sulfide (nano)phases. 13,43 The relevance of NOM as key modulator for particle growth and stability has been shown for various engineered nanoparticles [44][45][46][47][48][49][50] and other metal sulfides, e.g., ZnS and HgS. [51][52][53] Moreover, several studies have demonstrated the ability of NOM coatings to protect various nanoparticles from aggregation via electrostatic or electrosteric stabilization. [44][45][46][47][48][49][50][51][52][53][54] Up to now, knowledge on NOM effects on Cu x S nanoparticle formation and growth is, however, still lacking.
The fate of nano-Cu x S in variable redox environments can only be predicted reasonably if dissolution is taken into account. Even though bulk Cu x S solubility is extremely low at circumneutral pH and under anoxic conditions, 43 dissolution has to be considered as it can be triggered upon re-aeration of wetland soils in case of poorly crystalline nanoscale precipitates. When present in their nanoparticulate forms, the solubility of mineral phases can even be increased by several orders of magnitude. [55][56][57] It has been suggested that partial oxidative dissolution of metal sulfides only occurs in the presence of NOM, e.g. HgS dissolution by various DOM isolates. 58,59 Contrarily, dissolved CuĲII) has been detected after exposing sulfidized CuO nanoparticles to oxic solutions in the absence of NOM. 60 Based on these, to some extent contradictory findings, together with the observations of metal sulfides in oxic natural waters, warrants further studies on the role of NOM in nano-Cu x S dissolution.
Given the limited knowledge on Cu x S NP properties and behavior under relevant environmental conditions, the aim of our study was to track particle growth of Cu x S and to investigate the effects of absolute metal and sulfide concentration as well as NOM concentration on particle size, particle morphology, and colloidal stability of Cu x S over four weeks. For this purpose, we used transmission electron microscopy (TEM) and adapted a hyphenated size-exclusion chromatography (SEC)-ICP-MS technique, which allowed the detection and sizing of Cu x S nanoparticles between 2-50 nm in diameter. Additionally, oxidative dissolution experiments were conducted to study the influence of fulvic acid and dioxygen on Cu x S NP solubility.

Materials and methods
Materials Filtered (0.22 μm, sterile Millipak® 40, Millipore) ultrapure deionized water (Milli-Q®, 18 MΩ cm, Millipore) was used to prepare all solutions. All chemicals used were of at least analytical grade. For the preparation of anoxic solutions, ultrapure water was purged with N 2 for at least 2 h under vigorous stirring. Anoxic work was performed in a glovebox (Labmaster 130, MBraun, Germany) under N 2 atmosphere (O 2 < 1 ppm). All glassware was acid-washed by soaking overnight in 1 mol L −1 HCl/1.5 mol L −1 HNO 3 followed by three rinses with ultrapure water. All stock solutions were prepared in the glovebox and filtered (0.22 μm mixed cellulose ester membrane filters, BGB Analytik AG, Switzerland).
Stock solutions of dissolved CuĲII) were prepared with CuCl 2 ·2H 2 O (Merck). Using 0.1 mol L −1 HCl and 0.1 mol L −1 NaOH, the stock of 3-(N-morpholino)propanesulfonic acid (MOPS, Amresco) buffer was adjusted to pH 7.5. Likewise, Environmental Science: Nano Paper the stock of Suwannee River fulvic acid II (SRFA, Cat. No. 2S101F, IHSS), serving as the representative for NOM, was adjusted to pH 6.1. Closed air-tight, the SRFA stock was stored outside the glovebox in the dark at 4°C. Before use, the concentrations of CuĲII) and SRFA stocks were determined by inductively coupled plasma optical emission spectroscopy (5100 ICP-OES, Agilent) and by TOC analysis (DIMATOC® 2000, DIMATEC, Germany), respectively. Prior to the preparation of sulfide stock solutions, crystals of Na 2 S·9H 2 O (Sigma-Aldrich) were thoroughly rinsed with ultrapure water in order to remove potentially oxidized or polysulfidic surface coatings 61 and dried for at least 24 h within the glovebox. The concentration of the sulfide stock was determined via iodometric titration, which was standardized against a thiosulfate solution (Titrisol®, Merck). Stability of the sulfide stock over 4 weeks within the glovebox atmosphere was verified using the same method. Sulfide working solutions were diluted from this stock. Ammonium nitrate (NH 4 NO 3 , Merck) and FL-70 (Fisher Scientific, USA) were used to prepare the mobile phase for SEC-ICP-MS analyses. Citratecapped gold nanoparticles (Au NPs) with nominal sizes of 2, 5, 10, 20, 30, 40, and 50 nm as given by the manufacturer (BBI solutions, UK) served as the size calibration analogs for the Cu x S NPs during the growth experiment. A quality assurance report of these well-characterized Au NPs is provided by the manufacturer. 62 For the quantification of sulfur contents in the nanoparticles, sodium sulfate (Na 2 SO 4 , Merck) was used as the calibration standard. Calcium chloride (CaCl 2 ·2H 2 O) required for the coagulation procedure in the dissolution experiments was obtained from Merck (Germany).

Preparation of Cu x S suspensions for growth experiments
Transmission electron microscopy (TEM) and size-exclusion chromatography (SEC) coupled with inductively coupled plasma mass spectrometry (ICP-MS) were used to follow the size of the Cu x S nanoparticles over time. For this purpose, we conducted experiments on Cu x S nanoparticle formation from anoxic solutions with low copper (50 μmol L −1 ) and sulfide (100 μmol L −1 ) concentrations at a 1 : 2 ratio (henceforth denoted 'dilute suspensions'), in absence and presence of SRFA (0, 5, 50 mg C L −1 ). Additionally, the influence of higher copper (500 μmol L −1 ) and sulfide (1000 μmol L −1 ) concentrations (henceforth called 'concentrated suspensions') on particle characteristics was investigated. All experiments were conducted at pH 7.5 (1 mmol L −1 MOPS) in a 10 mmol L −1 NaCl electrolyte at 25 ± 2°C. Starting solutions were prepared in 100 mL serum glass bottles by mixing NaCl, MOPS buffer, SRFA, and CuĲII) at given concentrations. Subsequently, SĲ-II) was added in one spike while stirring vigorously (with a magnetic stirrer) to induce particle formation. After that, serum bottles were closed with a butyl rubber stopper (Rubberbv, Netherlands) and an aluminium cap (Chromacol 20-ACB, Thermo Scientific, Germany). The bottle headspace volume after addition of all aqueous reagents was 16 mL N 2 . To a large part of this study, the analyses (scanning transmission electron microscopy (STEM), energy dispersive X-ray spectroscopy (EDX)) were conducted with samples equilibrated for 24 h. Apart from these, growth experiments for selected treatments were conducted to investigate nanoparticle growth over time. For this, aliquots were withdrawn at selected time points (after 45 and 100 min, then 3, 5, 8, 24 and 48 h, then 5, 7, 14 and 28 d from physical triplicate samples for SEC-ICP-MS analysis and after 100 min, 24 h, 7 d and 28 d from one of the same triplicate samples for TEM analyses) for subsequent analyses using a plastic syringe (Amefa, Braun, Germany) equipped with a stainless steel needle (Sterican, Braun, Germany). For SEC-ICP-MS analysis, aliquots were diluted to ∼100 μg L −1 Cu x S (assuming a Cu : S ratio around 1.0) in 10 mmol L −1 NaCl and 1% FL-70 into 2 mL HPLC glass vials and closed with a plastic screw cap within the glovebox before they were transferred to the instrument. The chemical speciation was calculated for all prepared solutions with Visual MINTEQ 3.1 63 to compare experimental conditions with thermodynamic equilibria (Tables S6a-e and S7a-e †).

SEC-ICP-MS analysis
For size determination and quantification of the Cu x S NPs, we adapted a SEC-ICP-MS method from Zhou et al. (2017), 64 using Nucleosil-SiOH columns (250 × 4.6 mm) with pore sizes of 1000 Å or 4000 Å (Macherey-Nagel). The mobile phase was composed of 10 mmol L −1 NH 4 NO 3 and 0.5% FL-70 at a flow rate of 0.5 mL min −1 and the injection volume was set to 50 μL. The setup consisted of a quaternary solvent pump (1200 series, Agilent Technologies), housing one of the above mentioned columns, coupled via a regular PEEK tubing to a triple quadrupole ICP-MS/MS instrument (ICP-MS 8800, Agilent Technologies) operated in tandem mode with either H 2 or a mixture of O 2 and H 2 as reaction gas. The instrument was tuned daily before analysis with a multi-elemental solution and a solution containing Sc and Y (1000 μg L −1 and 50 μg L −1 , respectively) was continuously admitted postcolumn through a T-piece to account for sensitivity variations during each session. The same solution was used over 28 d to correct for Cu x S quantification. A summary of instrumental parameters is given in Table S1. † A series of 10 μg L −1 Au NPs was used for size calibration (R 2 > 0.95, Fig.  S7 †), while CuĲII) and SO 4 2− served as standards for quantification (range 0-500 μg L −1 , R 2 > 0.99). Matrices of all standards were matched with the sample matrix and prepared in 10 mmol L −1 NaCl and 1% FL-70, accordingly. The calibration series for quantification was measured once before every sample batch, whereas Au NP size calibration series was always measured at least once before and after each sample batch in order to check for minor retention time variations. In order to achieve conversion of peak retention times (t R ) into particle diameters by means of established calibrations using Au nanoparticles, chromatographic data were first extracted from the Agilent Mass Hunter software.

View Article Online
Then, Cu and S signals were corrected with the internal standards. After that, peak fitting of all Au NP and Cu x S NP peaks was carried out with Origin 2018 (version b9.5.0.193) to determine peak properties such as peak maxima, center of gravity (COG) of the peaks and their full width at half maximum (FWHM). Finally, retention times at peak maxima for the internal standard corrected Cu signals were used to determine SEC diameters.

TEM analysis
Cu x S particles were investigated by scanning transmission electron microscopy (STEM, HD2700-Cs operated at 200 kV, Hitachi, Japan) using a high-angle annular dark-field (HAADF) detector for image formation. Elemental composition of individual particles was determined by energy dispersive X-ray spectroscopy (EDX, EDAX, USA). Samples for TEM analyses were prepared in the glovebox. Nickel grids with a holey carbon support film (S147N4, 400 mesh, Plano) were functionalized with 0.1% (w/v) poly-L-lysine (Sigma-Aldrich) to establish a positively charged collector surface. A drop of 25 μL of Cu x S suspension collected from the growth experiments was placed on the grid, either undiluted (for the lower concentrated samples, final concentration in the droplet: ∼50 μmol L −1 Cu x S) or up to 10 times diluted (for the higher concentrated samples, final concentration in the droplet: ∼50-250 μmol L −1 Cu x S). Particles were allowed to attach to the grid surface for 10 min before the grid was washed with one drop of ultrapure water. Excess liquid from the grid was wicked away from underneath with a lint-free tissue after each step. Prepared grids were placed into plastic grid racks (Plano, Germany) wrapped in aluminum foil to prevent light exposure and stored under N 2 atmosphere in the glovebox until analysis. Recorded TEM images were processed using different image analysis software packages (Digital Micrograph v3.01, ImageJ Fiji v1.51n, Adobe Illustrator CC v2015.0.0, Adobe Photoshop CC v2005.1.2). Particle sizing was performed either by using the Nanodefine Particle Sizer plugin 65 incorporated into the ImageJ Fiji software or by manual size determination and counting using the Digital Micrograph software. A forward Fourier transformation was applied to high-resolution TEM images that showed crystal lattice fringes of Cu x S particles to obtain corresponding selected area diffraction (SAED) patterns. Calculations on particle size distributions and statistical significance tests were performed with Origin 2018. To test for significant differences in particle sizes among treatments, pairwise Wilcoxon-Mann-Whitney-U-tests were performed with the TEM particle size distributions acquired during the growth experiments (Table S8-S10 †).

XRD analysis
Identification of mineralogical phases was accomplished by X-ray diffraction (XRD, D8 Advance, Bruker) analysis. For this, 10 mL of the concentrated sample suspensions were filled into polycarbonate tubes and ultracentrifuged (40 000 rpm, 193 000g, 3 h, 15°C, fixed-angle rotor, Centrikon T-1080, Kontron). The supernatant was carefully discarded and the fragile pellet at the bottom of the tubes was resuspended in 1 mL of ultrapure water. The dispersed suspension was then centrifuged again for 30 min at otherwise same conditions. The supernatant was removed and the new pellet was dispersed in 0.2 mL ethanol. As a final step, 10 μL of the ethanolic suspension was repetitively (5-10 times) deposited and dried onto a zero background silicon wafer (orientation (510), Siltronix, France) until a thin Cu x S layer was visible. For the measurement, the loaded silicon wafers were analyzed in Bragg-Brentano geometry employing Cu Kα radiation (λ = 1.5418 Å, 40 kV, 40 mA) and a high-resolution energy dispersive 1D detector (LynxEye XE). X-ray diffractograms were recorded from 4-90°2θ in 0.02°steps with 10 s acquisition time per step. Qualitative analysis of the measured diffractograms was done via peak matching of the observed reflections with the ICDD PDF2 database used with the Bruker EVA software.

Dissolution experiments
To study the stability of Cu x S against oxidative dissolution, a set of Cu x S nanoparticle suspensions was prepared in triplicates in the absence and presence of fulvic acid (0, 5, 50 mg C L −1 ) and at lower (50 μmol L −1 CuĲII), 100 μmol L −1 SĲ-II)) and higher concentrations (500 μmol L −1 CuĲII), 1000 μmol L −1 SĲ-II)) in the same way as already described for growth experiments. The only two deviations from this procedure were the removal of the magnetic stirrer from the serum bottles (in order to prevent any potential adsorption of Cu x S on the stirrer surface) together with the 1 : 10 dilution of the concentrated samples to 50 μmol L −1 Cu x S, both after one hour to ensure comparable conditions. After 90 min from particle formation point, the N 2 atmosphere in the headspace of the closed serum sample bottles was changed to 21% (v/v) O 2 (synthetic air) by replacing equivalent volumes of N 2 with O 2 gas using a gas-tight syringe. The amounts of spiked O 2 exceeded the amounts required to theoretically oxidize all Cu x S in the samples 7.5-fold. Additionally, controls without oxygen gas spikes were run. Subsequently, suspensions were incubated in the dark under constant shaking at 25°C. Aliquots of suspensions were collected with a gas-tight glass syringe at different time steps over several weeks (after 1, 4, 14, 22, 28, 42, 63, 87 d) and immediately processed by (i) an induced coagulation of the precipitates in 1 mmol L −1 CaCl 2 solution followed by (ii) a filtration using nylon membrane filters (0.22 μm), thus effectively separating the initial solid from the liquid phase. Subsequently, Cu concentrations in the filtrate were determined with ICP-OES (method detection limit: 0.2 μmol L −1 Cu or 0.4% of total Cu). Preliminary tests with chemically stable, anoxic Cu x S suspensions verified the suitability of this coagulation and filtration method in retaining all nanoparticulate Cu x S (>1 nm) and even clusters between 0.5-1.0 nm as no Cu was found in the filtrates.

Effects of reactant and NOM concentrations on Cu x S nanoparticles
Influence of reactant concentrations. The morphology of Cu x S NPs formed in the absence of fulvic acid after 24 h equilibration was compared between dilute (50 μmol L −1 Cu) and concentrated (500 μmol L −1 Cu) suspensions (Cu/S = 0.5). In dilute suspensions, the shapes of most Cu x S nanoparticles remained irregular or roundish after 24 h and only a few larger particles exhibited triangular, hexatriangular or hexagonal shapes. Similar triangular and hexagonal shapes of CuS nanoparticles have been described before when formed in non-aqueous solvents. 66,67 In contrast, most nanoparticles in the concentrated suspensions (500 μmol L −1 Cu) had already developed clear hexagonal shapes after the same reaction time (Fig. 1). In general, Cu x S particles were plate-like as can be seen for example in Fig. 1 and S18 † where Cu x S NPs are observed standing upright on the grid. In the dilute suspensions, Cu x S particles had a median minimum Feret diameter (from now on denoted median TEM diameter) of d TEM = 3.9 nm with very few larger particles between 14 and 20 nm as compared to d TEM = 15.8 nm in concentrated suspensions with about 5% of the particles being larger than 45 nm ( Fig. S1b and S4b †). In both cases, the Cu x S nanoparticles appeared welldispersed and no evidence for particle aggregation was observed at near-neutral pH. The comparison of the TEM images recorded for dilute and concentrated suspensions revealed clear differences in the speed at which particle morphology evolved, and also evidenced pronounced differences in particle sizes.
Influence of NOM. The effect of fulvic acid concentration (0, 5 and 50 mg C L −1 ) on Cu x S morphology and size was evaluated for both dilute and concentrated suspensions after 24 h equilibration. In the dilute SRFA-free Cu x S suspensions, Cu x S particles exhibited both irregular or roundish and to a lesser extent triangular or hexagonal particle shapes. On the contrary, almost exclusively irregular or roundish particle shapes were observed when fulvic acid was present at 5 mg C L −1 (Fig. 1a and b). Likewise, Cu x S particles were notably smaller in dilute suspensions in the presence of fulvic acid (d TEM = 3.9 nm without vs. d TEM = 2.9 nm with 5 mg C L −1 SRFA). Similar to these dilute SRFA-free suspensions, a 10fold increase in reactant and SRFA concentrations of the SRFA-containing suspensions led to the formation of hexagonal platelets and larger particle sizes ( Fig. 1b and d).
Median TEM diameters were much larger in concentrated (d TEM = 15.8 nm without SRFA, d TEM = 20.8 nm with 50 mg C L −1 SRFA) than in dilute (d TEM = 3.9 nm without SRFA, d TEM = 2.9 nm with 5 mg C L −1 SRFA) suspensions (Fig. 1c, d, a and b). In all 24 h experiments with dilute suspensions, the fulvic acid hampered the development of the Cu x S in particle size and morphology. This effect was much less pronounced in concentrated suspensions. However, in these concentrated suspensions, the fulvic acid hindered the growth of large particles exceeding sizes of 50 nm, which were only found in SRFA-free suspensions ( Fig. S4b and S5b †).
A common observation for all investigated conditions was the extraordinary colloidal stability of the Cu x S precipitates. Even without NOM, the suspensions remained stable under anoxic conditions for several months, likely due to the strongly negatively charged surfaces of the Cu x S particles (zeta potentials around −50 mV, Table S5 †) under the experimental conditions (pH 7.5, 10 mmol L −1 NaCl). This is consistent with a reported point of zero charge (PZC) for Cu x S ranging between pH 1 and 3. 68-72 Even a decrease in pH to 6 did not cause notable aggregation and only occasional aggregates ∼100 nm were observed in an otherwise welldispersed suspension after 24 h (Fig. S10 †).

View Article Online
Structure and composition of the Cu x S nanoparticles Microscopic and spectroscopic characterization. The Cu x S nanoparticles were characterized by STEM and XRD to obtain information on their elemental composition, stoichiometry and crystal structure. Recorded EDX spectra confirmed that the nanoparticles consisted of copper and sulfur, but could not give definitive information on stoichiometry due to the relatively low counts (Fig. S6 †). Therefore, molar Cu : S ratios were determined via SEC-ICP-MS measurements in a later phase (Fig. S9, † see also section below) and were found to be close to 1.0 (std. dev. = 0.1, min = 0.8, max = 1.4). Highresolution STEM images of Cu x S nanoparticles (∼10 nm, 24 h) featured clear lattice fringes (Fig. 2). The identified lattice spacing of 0.19 nm corresponds to the crystallographic (110) plane d-spacing in the crystal structure of covellite (CuS) as well as to the (111) plane d-spacing in the crystal structure of yarrowite (Cu 1.12 S). Metastable Cu x S mineral phases such as djurleite (Cu 1.96 S), digenite (Cu 1.8 S), anilite (Cu 1.75 S), geerite (Cu 1.6 S), spionkopite (Cu 1.4 S) and yarrowite (Cu 1.12 S) have been reported in earlier studies. 73,74 Peak matching of the X-ray diffraction patterns of our samples suggested that the particles mainly consisted of covellite (CuS) but it cannot be excluded that metastable mineral phases coexisted during the first 24 h, since the theoretical positions of their XRD reflections partly overlapped with the observed broad diffraction peaks. The X-ray diffractograms depicted in Fig. S19 † show the structural evolution of the Cu x S particles over time and suggest a growth of crystalline domains (e.g., sharpening of (110) reflection peak with aging). This result together with STEM images confirmed that the larger nanoparticles developed crystalline domains.
Temporal change in color. In all Cu x S suspensions, a change in color from initially brownish to pale and then dark greenish color was observed. This change was most pronounced during the first 24 h and colors intensified upon aging for 4 weeks. Furthermore, the rate of this change in color increased with increasing Cu-S concentrations and was completed in less than 1 h in the concentrated, SRFA-free Cu x S suspensions. As optical properties of particles are related to their structure, information of particle development can be inferred from the observed change in color. 73 In the case of Cu x S, it is believed that soluble CuĲII)-aquo (e.g. [CuĲH 2 O) 6 ] 2+ ) complexes or organically complexed CuĲII) are the most common metal precursors in natural systems to form covellite, eventually. 26,75 As such, the formation of metal sulfide nanoparticles is a multi-step process by which (in the case of Cu x S) rapidly formed soluble copper bisulfide complexes react under the loss of water and protons to Cu 3 S 3 ring structures that condense into tetrameric clusters (e.g. Cu 4 S 5 , Cu 4 S 6 ). 76 Subsequently, reduction of CuĲII) to CuĲI) occurs within these clusters and induces a structural conversion of the initial CuĲII) from 5-or 6-fold coordination to CuĲI) with 3-and 4-fold coordination. 76 The structural unit is then formed by S 2 -layer incorporation leading to the final Cu 3 S-CuS 3 -S 2 -Cu 3 S-CuS 3 layer structure, e.g. for covellite. 73 The structural reordering of these layers causes visible color changes in concentrated Cu x S suspensions 73 as it was also observed in our experiments. In the mentioned study of Pattrick et al. (1997), 73 such poorly crystalline greenish precipitates of covellite-like structure developed from brownish, metastable primitive-type Cu x S precipitates that were lacking in tetrahedrally-coordinated Cu: the sites that link the S 2 -layers with the Cu 3 S-CuS 3 -layers in covellite. Upon aging, the structure of this primitive phase was reordered and the linking of these layers evolved. 73 The linkage of these layers has been shown to allow an electron transfer that induces a considerable shift in the visible View Article Online spectrum towards the green fraction, explaining the green color development. 77 Apparently, such a structural transition was also occurring in our study since the initial brownish poorly crystalline suspension evolved during 6 weeks into a slightly more crystalline one with covellite character and a greenish color (Fig. S19 †). The rapid primary color transition in our experiments suggests a fast structural reordering already within the first 24 h of growth.
Sizing of the Cu x S nanoparticles by SEC-ICP-MS and comparison with TEM The results presented above covered Cu x S particle properties at 24 h aging. Particle properties also varied over time during particle growth and aging as suggested by the change in colors as discussed above. Therefore, we used SEC-ICP-MS to investigate the particle growth and composition (molar Cu : S ratios, Fig. S9 †) over time in more detail. Typical chromatograms of the Cu and S signals recorded for a growth experiment at 100 min and 24 h are shown in Fig. 3a. Notice that the very large S peak starting at 5.75 min (Fig. 3a, inset) originates from the sulfonic acid group of the MOPS buffer. The S signal related to the Cu x S nanoparticles could always be unambiguously distinguished from these large S background peaks given the excellent separation efficiency of the SEC columns with respect to dissolved species vs. nanoparticles. However, the Cu signal overall and especially the background signal was less noisy than the S signal, which made the Cu signal best suited for further use in data analysis of SEC diameters. The recoveries of Au nanoparticles ranged from 82 to 110% (ionic Au: 93%) and were even better for Cu x S nanoparticles, i.e. 95-104% (ionic Cu: 109%; Table  S2 †), demonstrating the preservation and absence of loss during their transport through the size-exclusion column. Fig. 3b presents the Cu signals as corrected with the 45 Sc signals (using the same Sc solution over 28 d) at various time points along a growth experiment of Cu x S particles formed in a concentrated, fulvic acid rich suspension. Overall, differences between peak retention time among the growing particles were very distinct, gradually decreasing from 5.45 min after 45 minutes to 4.95 min after 28 d, which corresponded to an increase in average diameter (d SEC ± std. dev.) from 3.0 ± 0.2 nm to 16.4 ± 0.1 nm when applying the size calibration (Fig. S7 †) using the Au nanoparticles. Interestingly, the peak shape markedly changed between 3 and 8 h with a clear bimodal distribution at 5 h, highlighting an important growing stage over that time period. Retention time differences assessed by either the peak maxima (PM) or peak center of gravity (COG) varied by less than 1.3% (n = 42) between same-sized Au nanoparticles, thus allowing a precise estimation of the median diameter of the nanoparticle population. Even small retention time shifts of the peaks ≥0.06 min were detectable over the whole course of the experiment, suggesting that changes in SEC diameters in the sub-nanometer range were detected. The applied technique exhibited a very high repeatability. Based on the retention  time of the dissolved Au and Cu analyzed at the same concentration than the Au and Cu x S nanoparticles, the size detection limit (lower size cutoff) was estimated to be 1 nm for Au and 1.6 nm for Cu x S nanoparticles. The upper size cutoff of the method was not determined precisely but 50 nm Au NPs were still well detectable, even though peaks started to broaden above 40 nm (Fig. S8 †), and recoveries were still excellent (Table S2 †). We found a good agreement between the diameters derived from TEM and SEC analyses (Fig. 3c) for small Cu x S nanoparticles (2.5-10 nm) that showed a rather narrow distribution (Fig. S1-S3 †). In contrast, SEC diameters of samples containing few very large particles in an otherwise small particle size population led to mismatches with TEM diameters (cf. Fig. 3c purple dataset and Fig. S4 †). The reason is that the TEM particle size corresponds to the number based average particle size whereas SEC-ICP-MS gives a mass based average particle size. 31 As the mass of a particle scales with the third power of its size (diameter), the mass based average size of a particle population is shifted towards larger sizes compared to the number based average size. Hence, the large particles measured in SEC-ICP-MS increased the mass based average nanoparticle diameter compared to the number based average diameter. This effect becomes more pronounced with increasing polydispersity of the particle populations. Consequently, a larger discrepancy in diameters derived from SEC and TEM was found in Cu x S samples that were formed in concentrated suspensions. In an attempt to explain this discrepancy in the SEC and TEM central tendency values in more detail, we converted the number based TEM particle size distributions of the Cu x S suspensions into mass based TEM PSDs by assuming a spherical shape for all particles in the conversion (Fig. S11-S15 †). While this procedure reduced the discrepancy between SEC and TEM average diameters in the concentrated SRFAfree suspensions, it led to considerable disagreement among the calculated mass based TEM average diameters and the SEC diameters in all other suspensions (Fig. S11-S16 †). As particle shapes were different among the suspensions and changed during particle development (cf. Fig. 1, and S1-S5 † addressed in later sections), the applied conversion procedure assuming a spherical particle model was prone to errors in the calculation of mass based diameters, particularly with suspensions that contained particles whose shape deviated the most from spheres. This is reflected by a larger shift from the 1 : 1 line in the comparison of converted mass TEM diameters and SEC diameters (Fig. S16 †). In addition, recorded particle population statistics need to be considered with a few dozens to hundreds of particles in TEM datasets used for the conversion compared to millions of particles in SEC datasets. In this regard, much weight is given to a few large particles when applying this conversion, especially in low particle number TEM datasets, which then leads to a distortion in the calculated mass based TEM diameter. This clearly shows that such rather simple conversion procedures for particle size distributions should be applied and interpreted with great caution as respecting the correct particle geometry in x-, y-and z-direction as well as counting statistics are crucial.
The discrepancies were also slightly larger when peak maxima were used for the identification of peak retention times and the subsequent conversion into diameters. If instead the center of gravity (COG) obtained from peak fitting was used for the assignment of retention times and subsequent size conversion, average diameters of more polydisperse samples were in closer agreement with median TEM diameters (Fig. S17, Tables S3 and S4 †). For all other samples, the SEC diameters derived from peak maxima agreed with the median TEM diameters. Therefore, peak maxima were preferably used to convert SEC retention times to equivalent diameters for the tracking of the nanoparticle size development in the growth experiments. Future studies could make use of such a diameter determination from sizeexclusion chromatograms based on a central parameter of a Gaussian curve (used to fit the elution peak), in order to characterize different nanoparticle types with different coatings.
Nucleation and growth of Cu x S nanoparticles in the absence of NOM Dilute suspensions. The growth of Cu x S nanoparticles in the absence of SRFA was followed at different reactant concentrations over the course of 4 weeks (Fig. 4, gray symbols and lines). The average SEC diameter of Cu x S nanoparticles formed in dilute SRFA-free suspensions (Fig. 4a) measured after 100 min was d SEC = 3.0 nm. The corresponding median TEM minimum Feret diameter was only slightly larger (d TEM = 3.4 nm). Thereafter, the growth of the Cu x S particles proceeded slowly, with an average SEC diameter increase of less than 1 nm over the next 24 h. The major growth phase occurred between day 1 and day 7 with nearly a doubling in the average SEC diameter (from d SEC = 3.7 nm to d SEC = 7.1 nm). After 7 d, nanoparticle growth leveled off in all dilute suspensions, indicating that their size development was completed. Median TEM diameters for the selected dilute samples were very close to the SEC ones (max. difference <2 nm, Fig. 4a). TEM particle size distributions of all dilute suspensions widened in the course of the experiment without developing a clear bimodal size distribution (Fig. S1 †).
Thermodynamic calculations for the precursor solution (before sulfide addition) suggested the possible precipitation of 98.4% of the Cu as atacamite (Table S6 †), a copperĲII) chloride hydroxide mineral, and we actually found and isolated CuĲII)-hydroxide-like phases in the precursor solution after 30 min of equilibration time (Fig. S20 †). Due to kinetic constraints on (par)atacamite precipitation, 78 it is more likely that a simple CuĲII) hydroxide phase would form under these conditions. For nucleation, this means that the initial situation for metal sulfide formation is different than direct nucleation from solution as (par)atacamite seeds were already View Article Online present in the system. Since there were no visible precipitates in the precursor solution and sodium sulfide was added within a maximum timespan of one minute, the growth of the potential nucleation seeds for Cu x S was interrupted at a very early stage. We, therefore, assume that the formation of (par)atacamite was limited to nanoclusters or very small nanoparticles. In this case, growth initiation of Cu x S occurs via CuĲII) that is released from the (par)atacamite structure by dissolution, before the sulfide forms Cu x S complexes. Consequently, Cu x S growth would also depend on the actual size and number of the (par)atacamite seeds. In this scenario, however, the supply of Cu might be slower compared to the situation where Cu is present as an aqueous CuĲII) complex, thereby resulting in a reduced nucleation rate. Generally, if the number of nucleation sites is low, particle size at the end of the main growth phase will be larger, and vice versa. 79 Since reaction rates for metal sulfide cluster formation can be faster than 1 s −1 , 75 this primary growth phase was not resolved by our analytical methodology. According to Luther III (2016), 75 cluster size definition ranges from 0.4-1.0 nm. In the SRFA-free suspensions, the smallest Cu x S nanoparticles in the first TEM sample (after 100 min) showed minimum Feret diameters between 1.1-1.5 nm, indicating that cluster formation was already finished at this stage. With increasing reaction time the particle size distributions of Cu x S in SRFA-free suspensions considerably broadened, which we believe is due to an Ostwald ripening process that favors the growth of larger particles at the expense of smaller particles, which dissolve, reform clusters and contribute to the growth of the larger particles. [80][81][82][83][84] The growth phase eventually ends due to the termination of the Ostwald ripening process, [84][85][86] when the supply of small particles is exhausted and the dissolution of the bigger particles becomes negligible (Gibbs-Thomson effect 86 ).
Concentrated suspensions. With a 10-fold increase in copper and sulfide concentrations in the absence of SRFA, the first measured median TEM diameter (after 100 min) was 7.1 nm and it rapidly increased within the first 8 to 24 h to 15.8 nm (Fig. 4b). At this time point a plateau was reached and the median TEM diameter remained constant until the end of the experiment. Unlike in dilute suspensions, TEM particle size distributions of the concentrated Cu x S suspensions showed a high degree of polydispersity (Fig.  S4 †). In those, both very small particles and relatively large particles coexisted in these suspensions, especially in the initial phase. In this phase, TEM images of the suspensions showed a formation of small aggregates only consisting of the smallest particles of the sample. The contribution of small particles continuously diminished over time as they most likely dissolved due to Ostwald ripening and the number of large hexagonal Cu x S particles increased. Overall, both the initial and final size of the Cu x S nanoparticles in concentrated suspensions were higher than in the dilute suspensions and the growth was much faster.
The relative proportions of precipitated (par)atacamite predicted by thermodynamic calculations were identical to dilute suspensions (Table S6 †) and can therefore not explain the observed differences in particle sizes. However, the supersaturation with respect to (par)atacamite was higher in concentrated suspensions. Consequently, more (par) atacamite seeds may have been formed at the beginning in concentrated suspensions. 87 When compared to the dilute suspensions, the (par)atacamite seeds may take longer to dissolve when the sulfide is added to the system to form Cu x S particles. Such a retarded nucleus supply at initial stage in higher concentrated systems has occasionally been observed to induce the formation of condensed nucleus clusters at the cost of 'true' nucleus abundance. 87,88 With less nuclei at growth initiation, the terminal particle size is expected to be larger, as was observed in our experiments with concentrated SRFA-free suspensions. 79 With an increase in concentrations of the main reactants, a much larger reservoir for the View Article Online reaction of copper and sulfide is available, resulting in a larger and prolonged supply of reactants during growth. On this account, some effects already observed in the SRFA-free experiments in dilute suspensions are magnified due to the extreme affinity of CuĲII) for SĲ-II), which drives crystal development and more rapidly leads to higher crystalline structures than in dilute Cu x S variants (cf. hexagonal platelets already after 100 min in Fig. S4 † vs. the still more spherical particles after 100 min in Fig. S1 †). In concentrated suspensions, the color change from brownish to greenish was much more rapid (<1 h in concentrated suspensions vs. ≥24 h in dilute suspensions) indicating a strongly accelerated restructuring process from soluble copper bisulfide complexes over tetrameric clusters to the Cu 3 S-CuS 3 -S 2 -Cu 3 S-CuS 3 layer structure of covellite. The increased supply of Cu assures a more continuous and sustained delivery of new reactants onto the surfaces of the precipitating clusters and nanocrystals promoting the fast structural and sizerelated development as observed in our experiments. 89 Hence, the aforementioned plateau and thus also the final particle size is reached much faster for all concentrated suspensions.

Effect of NOM concentrations on Cu x S nucleation and growth
Dilute suspensions. The influence of different fulvic acid concentrations on the growth of Cu x S nanoparticles from dilute suspensions is depicted in Fig. 4a (blue and red symbols/lines). Average SEC diameters of Cu x S nanoparticles formed in the presence of 5 or 50 mg C L −1 of SRFA (measured after 100 min) were 2.5 and 2.2 nm, respectively. Median minimum Feret TEM diameters of the same samples were slightly larger (3.1 and 2.4 nm, respectively). Significant differences in average SEC diameters among samples with and without fulvic acid were already observable during the first 8 h, with increasing SRFA concentrations resulting in smaller Cu x S particles. From day 7 to the end of the experiment, average SEC diameters of Cu x S in samples containing SRFA were always notably smaller than in SRFA-free suspensions. After 4 weeks, Cu x S particles formed in dilute suspensions with the highest SRFA concentration were up to 25% smaller than the ones formed in the absence of SRFA. Similar to their SRFA-free counterparts, the TEM particle size distributions of the dilute, SRFA-containing suspensions broadened during the 4 weeks ( Fig. S2 and S3 †). However, this broadening was less pronounced when higher fulvic acid concentrations were present in the Cu x S suspensions.
The relative proportion of Cu initially precipitated as metal hydroxide depends on the concentration of SRFA. According to our thermodynamic calculations for the precursor solution before adding sulfide, 81% of the CuĲII) may have been bound in (par)atacamite nano-seeds in the presence of 5 mg C L −1 SRFA whereas the remaining fraction (19%) may have been bound by SRFA (Table S6 †). CuĲII) exhibits a high affinity to major carboxylic and phenolic functional groups, but also to minor N-containing and reduced-sulfur (thiol) functional groups in fulvic acids. 42,90,91 Accordingly, the precipitation of (par)atacamite was almost completely suppressed at the high SRFA concentrations (50 mg C L −1 ) as 99.8% of the Cu was predicted to be present as CuĲII)-SRFA-complexes (Table S6 †). Therefore, the formation of Cu x S nanoparticles from (par)atacamite nano-templates was only expected in the low SRFA case (5 mg C L −1 ) while Cu x S particle growth at high SRFA concentration (50 mg C L −1 ) was primarily initiated from SRFA-CuĲII)-SĲ-II) complexes. An excess in SRFA was available in both cases for surface capping of the nanoparticles as only a small fraction of the fulvic acid was needed to bind the Cu in each case at initial formation stages (29% for low SRFA, 17% for high SRFA). The growth of Cu x S in dilute suspensions was effectively slowed down by the fulvic acid as the first measured diameters were significantly lower compared to the ones found in the corresponding SRFA-free suspensions. Cu x S nanoparticles with minimum diameters between 0.5-1.0 nm were occasionally observed in the first TEM sample (after 100 min), only in the presence of SRFA. This indicates that the cluster formation was not yet finished and that the length of the nucleation phase was longer when SRFA was present, as also discussed by Luther III (2016). 75 If we consider nucleation from the different initial CuĲII) species, a nucleation of Cu x S from SRFA-CuĲII) complexes would not only proceed faster than from (par)atacamite templates, it would also yield more nuclei in the same period of time. As SRFA concentration influenced the concentration level of (par)atacamite template seeds in the suspensions at initial stage, only 19% of the CuĲII) was readily available as SRFA-CuĲII) for cluster formation in suspensions containing 5 mg C L −1 SRFA, whereas most CuĲII) (99.8%) was quickly available as SRFA-CuĲII) for cluster formation in the high SRFA variants (50 mg C L −1 SRFA). Consequently, the number of nuclei at t = 0 was greater under the latter conditions, thus leading to smaller particle sizes after 4 weeks. This is consistent with the results of our experiments, where the final particle sizes of Cu x S formed in dilute, SRFA-free suspensions were larger than in all NOM spiked suspensions that exhibited lower supersaturation levels.
SRFA may have interacted with early Cu x S polynuclear complexes and clusters through chelation with (surface) copper cations or through sorption of hydrophobic fulvic acid sites to the newly formed clusters. 92 Even though termination of Cu x S surfaces by sulfur is expected due to the sulfide excess in the reactors, these early SRFA-CuĲII)-SĲ-II) associations may delay the process of further cluster formation especially in the beginning by steric hindrance as also observed by Poulin et al. (2017) in experiments involving DOM-HgS associations. 93 If so, electrostatic repulsion of new clusters (negatively charged) and already existing SRFA-CuĲII)-SĲ-II) associations (negatively charged) additionally retard the reaction of CuĲII) and SĲ-II). In this way, SRFA inherently caps newly formed Cu x S surfaces and protects the clusters and nanoparticles from extensive growth. Furthermore, it protects the nanoparticles from attachment to other Cu x S particles. As a consequence, the presence of fulvic acid limited the broadening of the particle size distribution and hampered the development of large particles. Since Cu x S particles were significantly smaller in dilute suspensions containing SRFA, a more pronounced Gibbs-Thomson effect may have facilitated Ostwald ripening. Additionally, SRFA functional groups, known to be strong Cu ligands, could have stimulated dissolution-reprecipitation reactions. Moreover, we observed a distinct deceleration of structural and morphological development of Cu x S particles with the increase in fulvic acid concentration. A similar effect has been recently discovered in aged HgS-DOM suspensions. 93 In accordance with Lee et al. (2005), 83 the lack of hexagonally shaped particles with a higher structural order at the highest NOM concentration may be a hint for a strong contribution of Ostwald ripening in these suspensions. Surface ligands such as fulvic acid can also selectively bind to certain facets of growing nanocrystals to minimize their surface site energy. 94 Like this, reactant or cluster delivery to these surfaces can be diminished or even be completely blocked, leading to a slower and restricted structural development with altered particle shapes compared to ligand-free suspensions 89 as seen in our dilute Cu x S suspensions. Concentrated suspensions. The growth of Cu x S in concentrated suspensions containing 50 mg C L −1 SRFA is shown in Fig. 4b (orange symbols/lines). After 100 min, measured diameters were between d SEC = 3.3 nm and d TEM = 3.5 nm and increased to d SEC = 15.1 nm and d TEM = 20.8 nm within the first 8 to 24 h. At this time point, nanoparticle growth leveled off and the average SEC diameter remained nearly constant. Only the median TEM diameter of these suspensions decreased to 17.9 nm after 4 weeks, a value close to the corresponding SEC diameter (d SEC = 16.4 nm, 28 d). A temporary bimodal phase in the development of the particles was detected with SEC as the 5 h chromatogram (Fig. 3b) constituted the rapid shift from a particle population that was mostly consisting of small particles <10 nm to a population that was dominated by larger particles ≥20 nm (cf. Fig. S4a and b †). This considerable change in the particle population coincided with median TEM and average SEC diameters recorded in the abovementioned intensive growth phase (Fig. 4b).
A comparison of the concentrated Cu x S suspensions showed that after 4 weeks both the SRFA-free and SRFAspiked samples were in a similar size range (d TEM = 15.5 nm, no SRFA vs. d TEM = 17.9 nm and d SEC = 16.4 nm, with SRFA). When comparing the dilute SRFA-spiked suspensions to their concentrated counterparts, a 10-fold increase in reactant concentrations increased the particle size of Cu x S in the most pronounced case by a factor of about 2.5 after 4 weeks irrespective of which measure was used (d TEM = 6.9 nm and d SEC = 6.5 nm in dilute, 5 mg C L −1 SRFA vs. d TEM = 17.9 nm and d SEC = 16.4 nm in concentrated, 50 mg C L −1 SRFA). With regard to the rapid transition in the particle populations from the initial to the plateau phase, we suggest that our SEC chromatograms show the ongoing process of Ostwald ripening as we observed the decline in contribution of smaller particles <5 nm which almost disappeared in later chromatograms and TEM particle size distributions.
Based on the capability of the fulvic acid to restrict particle growth in dilute suspensions, it is surprising that this does not seem to apply to the concentrated suspensions in the same way. In high concentrated suspensions, the fulvic acid may stabilize the Cu x S particles only after crystal development into larger particles (d TEM_high_28d = 17.9 nm vs. d TEM_low_28d = 6.9 nm). In this case, the comparably large reactional driving force in the concentrated Cu x S suspensions would favor the direct reaction of CuĲII) and SĲ-II) at the expense of early stage SRFA-CuĲII) complexation to functional groups that are responsible for effectively decreasing the size and structural order of the crystals. At CuĲII) concentrations this high, these high affinity CuĲII) binding sites can be saturated so that their size-and structure-restricting activity is considerably impeded. 93 Poulin and coworkers have found similar decoupling of such size limiting effects of NOM when reactant concentrations exceeded a critical level at which saturation of strong binding sites prevailed. 93 Additionally, the fast reaction kinetics in concentrated SRFA-containing suspensions provoke collisions that are more effective in removing surface bound fulvic acid. This also provides a reasonable explanation for the indifference in final Cu x S particle size among concentrated suspensions in absence and presence of SRFA. Interestingly, even though the effectiveness in decreasing particle size is largely impaired, the occurrence of extraordinary large particles was still suppressed in the presence of SRFA, suggesting that other than high affinity functional groups are responsible for the limitation of

Dissolution of Cu x S nanoparticles
The oxidative dissolution of Cu x S nanoparticles in the absence and presence of fulvic acid and with Cu x S particles of different morphologies is illustrated in Fig. 5. Under anoxic conditions, no dissolved Cu was detected (<0.2 μmol L −1 or <0.4% of total Cu) over several months at circumneutral pH, irrespective of the presence of fulvic acid. Cu x S particles remained stable when exposed to O 2 in the absence of fulvic acid and only the concurrent presence of O 2 and SRFA led to their dissolution. In this case, the dissolution proceeded quickly (20 days) and levelled off between 22 μmol L −1 and 26 μmol L −1 , values close to the Cu-SRFA complexation limit of the system (22.6 μmol L −1 Cu). The remaining solid fraction was determined by XRD and consisted exclusively of covellite (Fig. S21 †). No significant difference was found among differently treated SRFA-containing samples. Accordingly, it was irrelevant whether fulvic acid was added before or 1 hour after particle formation. Likewise, the reaction was independent of the particle aging since the original dilute suspensions (containing more roundish-platy or nearly spherical shapes) behaved similar to samples from 10-fold diluted high concentrated suspensions (containing larger amounts of particles with hexagonal symmetries).
In the oxic, SRFA-free media, the particles maintained their integrity owing to the generally very low solubility of Cu x S phases, that prevents dissolved O 2 to successfully disintegrate the particle surface. The excess of reduced dissolved sulfur species in the reaction medium was most likely oxidized first, before S-terminating surface sites of the Cu x S particles underwent oxidation, thus hindering dissolution. Our finding is opposed to the experiments of Sukola et al. (2005), who observed partial oxidative dissolution of CuS in synthetic suspensions that did not contain DOM. 95 However, CuS samples were not kept strictly in the dark in their experiments, which could have provoked photochemical reactions facilitating dissolution and would be an explanation for the observed differences. An oxidation of the particle surface would lead to the transformation into a copper hydroxide phase on the surface, which would possess a much higher solubility than any Cu x S phase. This may facilitate the dissolution promoted by strong ligands of the fulvic acid that help disintegrate the particle by liberating Cu at defect sites via complexation. Samples that did not contain fulvic acid were lacking this trigger for a ligandpromoted dissolution of the oxidized particle surfaces. The dissolution process continues until the maximum binding capacity of the fulvic acid for CuĲII) is reached, which corresponds to the Cu level of the plateau in the dissolution data. Depending on the batch variant taken for calculation, the half-life of Cu x S was determined to be 27-32 d, which agreed well with the half-life of 22 d found for CuS in natural river water containing NOM, too. 25 In natural waters, dissolution of Cu x S could be influenced by its spatial residence as UV-light absorption by abundant aquatic species (DOM, NO 3 − and FeĲIII)-complexes) leads to the formation of reactive oxygen species (ROS) in the upper layers of surface waters, which can take part in the oxidative dissolution process. 96 In addition, several studies have identified oxicanoxic interfaces in natural environments as hotspots for dark HO˙formation when DOM and metal-based reductants (such as CuS) are present. 97,98 A recent study even observed the formation of ROS in the dark in oxic solutions containing Cl − ions, SRFA and CuĲII), thus interfering in the redox cycle of such systems. 99 Furthermore, the presence of different NOM types that exhibit increased binding affinities towards CuĲII), such as humic acids, 40,42,100 could enhance Cu x S dissolution.

Environmental implications
In this study, we demonstrated that fulvic acid as a NOM representative influences the formation of Cu x S colloids by altering their nanoparticle size and structural development. Fulvic acid hampered the growth of the Cu x S particles so that their sizes remained smaller than in fulvic acid free suspensions. Herein, the absolute concentrations of Cu and S had a substantial effect on the temporal evolution of Cu x S particle development, resulting in morphologically more mature particles at higher reactant concentrations. Additionally, we observed a remarkable colloidal stability of Cu x S suspensions at near neutral pH values relevant to natural or farmed wetland systems. 10 Concerning oxic conditions, NOM acted as a key driver in the oxidative dissolution of Cu x S nanoparticles. The growth-suppressing properties of NOM impact the potential mobility of small (<10 nm) Cu x S nanoparticles forming at concentrations typical to soil pore waters, and also restrict size of particles precipitating from heavily contaminated sites with elevated Cu concentrations. This results in a high potential mobility of Cu x S nanoparticles. Our findings suggest that Cu x S will not grow into aggregated microparticles in natural environments, even at elevated Cu and bisulfide concentrations. Due to the strongly negatively charged surfaces resulting in an extraordinary long-term colloidal stability, Cu x S is expected to contribute to the longrange transport of Cu, especially in NOM-rich, oxygendepleted waters. In addition, our results provide strong evidence that CuĲII) is released from Cu x S nanoparticles in NOM-rich oxic aquatic environments or upon aeration of (fertilized) wetland soils.

Conflicts of interest
There are no conflicts of interest to declare.