Increased Mobility of Metal Oxide Nanoparticles Due to Photo and Thermal Induced Disagglomeration

Abstract

Significant advances have been made on our understanding of the fate and transport of engineered nanomaterials. One unexplored aspect of nanoparticle aggregation is how environmental stimuli such as light exposure and temperature variations affect the mobility of engineered nanoparticles. In this study, TiO2, ZnO, and CeO2 were chosen as model materials for investigating the mobility of nanoparticles under three external stimuli: heat, light and sonication. Sunlight and high power sonication were able to partially disagglomerate metal oxide clusters, but primary particles bonded by solid state necks were left intact. A cycle of temperature increase from 25°C to 65°C and then decrease back was found to disagglomerate the compact clusters in the heating phase and reagglomerate them as more open fractal structures during the cooling phase. A fractal model summing the pair-wise DLVO interactions between primary particles within two fractal agglomerates predicts weak attractions on the order of a few kT. Our study shows that common environmental stimuli such as light exposure or temperature variation can disagglomerate nanoparticle clusters and enhance their mobility in open waters. This phenomenon warrants attention since it is likely that metal oxide nanoparticles will experience these natural stimuli during their transport in the environment.

About the Authors:

Dongxu Zhou

Samuel W. Bennett

Arturo A. Keller

* E-mail: [email protected]

Affiliation: Bren School of Environmental Science and Management and University of California Center for Environmental Implications of Nanotechnology, University of California Santa Barbara, Santa Barbara, California, United States of America

Introduction

Recently a great amount of effort has been devoted to evaluating the potential environmental impacts of nanotechnology [1]–[4] due to the increasing applications of nanomaterials. Understanding nanoparticle aggregation in the aqueous environment is critical for assessing the fate, transport and toxicity of nanomaterials [5]. Studies conducted in different aqueous media, including synthetic matrices [6]–[16], natural waters [17], [18], and culture media [19], [20] have indicate that pH, ionic strength [7], [9]–[11], nanoparticle concentration [10], and natural organic matter [6], [8] affect nanoparticle aggregation. However, questions such as how the sintered structure and the presence of clay minerals affect the aggregation process remain to be answered.

Nanoparticles may form aggregate structures during the synthesis process. The bonding attractions between primary particles range from strong chemical bonds, established during the cooling phase of their synthesis, that form essentially unbreakable aggregates to weak physical interactions such as van der Waals forces which form agglomerates that are readily disrupted by sample preparation methods [21]–[23]. For metal oxide nanoparticles that are produced via high temperature processes, chemical bonding form during the synthesis stage when the temperature is high enough for sintering but not high enough for fast coalescence; agglomerates form at a later stage when the temperature decreases below sintering [21]. Following the suggestion of Mandzy et al. [23], we refer to particle clusters bound by irreversible chemical bonds as aggregates, and those held together by weak physical interactions as agglomerates. Once released in the environment, nanoparticles will very likely exist as agglomerated aggregates, i.e. aggregate clusters that have weaker bonds between them. This structural conformation and its effects on the stability and mobility of nanoparticles have not been addressed to any significant extent in the environmental context.

In this study, widely used metal oxide nanoparticles (TiO2, ZnO, and CeO2) were selected to study how the agglomerated aggregate structure controls nanoparticle mobility. We first probed the metal oxide agglomerate/aggregate structure via sonication and light exposure. Then the sedimentation behavior of the fractal agglomerates and aggregates was examined, followed by investigating the response to temperature variation. We report for the first time that the temperature variations can cause either agglomeration or disagglomeration of agglomerated aggregate structure depending on the heating and cooling paths. This finding is very relevant in evaluating the transport of nanoparticle transport, since it indicates that ambient temperature change, constantly occurring in open waters, can alter nanoparticle mobility. Finally a fractal aggregate model was developed to better understand the interparticle interaction between two clusters and provide an explanation for the observed phenomena, i.e. thermally-induced disagglomeration.

Materials and Methods

Materials

P25 TiO2 dry powder was obtained from Evonik Degussa (USA). CeO2 and ZnO dry powders were obtained from Melorium Technologies (USA). The characterization of the three metal oxides is presented in Table S2. The stock suspension was prepared by weighing the nanoparticle dry powder and suspending it in NanoPure water (NanoPure Diamond, Barnstead, MA) to achieve a 1.0 g/L concentration. The dispersion was then sonicated in a sonication bath for 30 min (Branson 2510, total power output 100 W, Danbury, CT). Chowdhury and coworkers found that 30 min in a sonication bath produced metal oxide nanoparticle suspensions with minimal hydrodynamic diameter, and further sonication led to reagglomeration of the nanoparticles [13]. The supernatant was then taken out as stock. Fresh stock suspensions were prepared daily. Samples were prepared by directly diluting the stock suspension with NanoPure water to achieve a 100 mg/L concentration. Because of the simple composition of the samples, the pH and conductivity of the samples were consistent and stable with no further pH or ionic strength adjustments (Table S2).

Aggregate size and zeta-potential characterization

Dynamic light scattering (DLS) (Zetasizer Nano ZS-90, Malvern Instruments) was used to determine the hydrodynamic size of the metal oxide nanoparticles. A 633 nm laser source and a detection angle of 90° were used. The intensity weighted mean hydrodynamic diameter and polydispersity index (PDI) were obtained. The PDI values of all of the measurements presented in this study were below 0.4. Data was collected for 30 s for each sizing measurement. Laser Doppler Velocimetry was used to measure the electrophoretic mobility (EPM) of the nanoparticles using the Zetasizer. EPMs were converted to zeta-potential using the Smoluchowski equation [24]. Most studies were conducted at 25°C in the thermo-regulated chamber of the Zetasizer, except the set of experiments investigating the effect of heat, as described below.

Sonication experiments

In an effort to break the agglomerates more, a 1 mL sample (100 mg/L metal oxide) was subjected to high energy sonication using a sonication probe fitted with a microtip (S-4000, Misonix Ultrasonic, USA). Samples were placed in the Zetasizer sampling chamber for sizing immediately after 2 s of sonication at 40% amplitude (input power 7 W). The input power was chosen based on preliminary results indicating that higher sonication power could not decrease the hydrodynamic size further. The 2 s sonication duration was chosen to avoid overheating, which could lead to fast aggregation and false DLS readings. Several cycles of sonication/sizing measurements were completed until no further decrease in sample hydrodynamic size was observed.

Sedimentation experiments

Metal oxide stock suspensions (1 g/L) were observed to settle slowly over the time course of days. To quantify the sedimentation process of the agglomerated aggregates, the hydrodynamic size of 100 mg/L metal oxide samples diluted from freshly prepared stock was monitored by DLS for a week.

Light exposure experiments

Cuvettes with 1 mL sample (100 mg/L metal oxide) were placed in sunlight for 30 min. Preliminary results showed that the bulk temperature change was less than 5°C after light exposure. Samples were subsequently placed in the Zetasizer's thermo-regulated chamber for 5 min to allow for temperature equilibration. Preliminary experiments showed that 5 min was sufficient for samples to reach the 25°C measurement temperature. Thirty 30 s sizing measurements were then carried out and the results were averaged. The experiments were carried out between 10 am to 2 pm with no cloud cover to ensure sufficient sunlight intensity.

Temperature experiment

A cuvette with 1 mL sample (100 mg/L metal oxide) was placed in the Zetasizer chamber, and the temperature of the chamber was programmed to increase from 25°C to 65°C and then decrease from 65°C to 25°C at 1°C intervals. At each temperature, a 120 s equilibration period was allowed before size and zeta-potential data were collected in triplicate.

TEM imaging

Transmission electron microscopy (TEM) samples were prepared by placing a drop of nanoparticle suspension on a copper grid (Ted Pella, CA) and letting it air dry overnight. Imaging was performed on a JEOL 123 microscope operated at 80 kV (JEOL, USA).

Fractal model

Aggregates are generally fractal in nature [25], where the structure can be characterized by an exponential relationship [26],(1)where N is the number of primary particles in an aggregate, kg is the fractal prefactor, Rg is the radius of gyration, a is the radius of a primary particle, and dF is the non-integer dimensionality. For an aggregate with the coordinates of each primary particle known, Rg can be determined by [27],(2)where dij is the distance between primary particle i and j.

A number of studies have reported the successful application of DLVO theory [28] to evaluate colloidal stability of nanoparticle systems [14], [16], [29]–[32] by considering van der Waals attraction and electrostatic repulsion. Studies to date use DLVO calculation results between primary particles to predict the aggregation behavior of nanoparticle suspensions. Most studies have assumed that the energy-separation distance profile between primary particles can represent the interaction between aggregates [29]–[32]. However, the validity of this assumption is questionable. To evaluate this assumption we developed a 3-D MATLAB simulation code to investigate the DLVO interactions between aggregates as the agglomeration process evolves, considering a fractal configuration following Brasil, et al. [26].

The logic in the MATLAB code is as follows. Primary particles are placed in 3D space one at a time. The first primary particle is placed at the origin (0, 0, 0). For the second and subsequent primary particles, random numbers are generated for the coordinates with the following two conditions: (a) touch at least one of the previously placed particles; and (b) not overlap with any of the existing particles. Then Rg is calculated using Equation (2) and dF is calculated with Equation (1). If the calculated dF is within 5% difference of the predefined dF, the new coordinates are accepted and the code adds the next primary particle; otherwise the new coordinates are rejected, another set of coordinates is generated, and the dF evaluation is repeated. Figure 1A shows a typical fractal aggregate generated by the code, and Figure 1B a scanning electron micrograph of a typical TiO2 aggregate in an aqueous matrix, for comparison. After generating two fractal aggregates, the DLVO energy-separation distance profiles between them is evaluated by calculating pair-wise DLVO interactions between individual primary particles within the two fractal aggregates. The DLVO parameters for each system are presented in Table S1. The solution chemistry in the simulation can be adjusted by changing ionic strength and zeta-potential. At each condition the simulation was run 1000 times and the results were averaged. Rg was converted to hydrodynamic radius (RH) using a Rg/RH conversion ratio determined by a previous study [33] to compare simulation results with experimental data.

[Figure omitted. See PDF.]

Figure 1. A: A representative aggregate generated by the Matlab code; B: SEM image of a TiO2 aggregate (Image adapted from Thio, et al. [33], with permission).

https://doi.org/10.1371/journal.pone.0037363.g001

Results and Discussion

Probing the agglomerate/aggregate structure by sonication

The nanoparticle suspension procedure using just the sonication bath for 30 min generated TiO2 agglomerates with hydrodynamic diameters between 250–350 nm. Several cycles of sonication/size measurement using the Misonix sonication probe were conducted until no further size decrease could be observed (Figure 2). Cavitation via sonication provides mechanical and thermal energy that promotes disagglomeration. The Misonix sonication probe (input power 7 W/mL sample volume) was unable to completely break up the agglomerates to primary particles, with the smallest observed agglomerate size as ca. 240 nm.Similar results were observed for ZnO and CeO2 (Figure S1).

[Figure omitted. See PDF.]

Figure 2. Sonication effect on particle size for a fresh stock of TiO2.

https://doi.org/10.1371/journal.pone.0037363.g002

Since DLS measurements are temperature sensitive, an increase in bulk temperature can lead to an incorrect size measurement, resulting in an artificially decreased particle size due to increased particle diffusion. Preliminary experiments with longer sonication times (5 s and 1 min) using the Misonix sonication probe heated the sample considerably leading to incorrect readings (Figure S2).

Given the high power input of the sonication probe, it is reasonable to assume that clusters without the sonication cycles were soft agglomerates of the aggregates (hydrodynamic diameter denoted as ds) while clusters after the sonication cycles were hard aggregates (hydrodynamic diameter denoted as dh). The ratio h = ds/dh represents the degree of agglomeration, which is 1.29 for the sample in Figure 2.

Sedimentation

A significant decrease in ds, from 370 nm to 260 nm, was observed during the first 6 hr of the sedimentation experiments (Figure 3A). PDI decreased from 0.3 to 0.2. Subsequently, ds decreased at a much slower rate. Over the next 100 hr ds decreased from 250 nm to 200 nm, while PDI decreased from 0.2 to 0.1. In DLS, the hydrodynamic diameter (dH) is calculated assuming a single size population following a Gaussian distribution, while PDI = σ2/dH2, where σ is the standard deviation of the Gaussian distribution. Based on the sedimentation equation that Kajihara developed for porous particles [34], it takes roughly 10 days for a 200 nm diameter TiO2 agglomerate (density 4.23 g/cm3) to settle 1 cm (the sample depth in this study), while a 1000 nm diameter agglomerate will settle the same distance in only 20 hr (Figure S3A). A freshly prepared TiO2 suspension contains many large clusters, which settled down relatively quickly. As the sedimentation process continues the percentage of larger clusters decreases, therefore ds and PDI decrease. The sedimentation of a sonicated TiO2 sample (sonicated by sonication probe for a few cycles to break agglomerates) was also monitored by DLS for 10 hr (Figure 3B). Distinctively different from sedimentation of the freshly prepared sample, no noticeable sedimentation was observed for the sonicated sample.

[Figure omitted. See PDF.]

Figure 3. Sedimentation monitored by DLS for (A) TiO2 agglomerates, prepared by diluting a TiO2 stock suspension; (B) TiO2 aggregates, prepared by diluting a TiO2 stock suspension and sonicating with a sonication probe several times until no further reduction in hydrodynamic size was observed. (C) Hydrodynamic diameter of TiO2 suspension supernatant monitored for 6 days.

https://doi.org/10.1371/journal.pone.0037363.g003

Daily monitoring of the supernatant of a 1 g/L TiO2 stock suspension for 6 days, measuring dz and PDI for sonicated and unsonicated samples, indicated that ds decreased from above 400 nm to around 200 nm, while dh decreased from 250 nm to 160 nm (Figure 3C). The degree of agglomeration, h, decreased from 1.68 to 1.17, indicating percentage-wise more agglomerates than aggregates settled during the course of 6 days. TEM images revealed that agglomerates form compact structures while aggregates have rather open structures (Figure 4). The sedimentation rate of fractal (porous) clusters is strongly affected by the agglomerate porosity [34]. Figure S3B shows the calculated sedimentation rate as a function of fractal dimension for clusters of equivalent hydrodynamic size (dH = 200 nm). A cluster with dF = 2.1 settles 3 times faster than a cluster with dF = 1.5. This explains the decreased h during the 6 days.

[Figure omitted. See PDF.]

Figure 4. TEM micrographs of TiO2 agglomerate and aggregate.

https://doi.org/10.1371/journal.pone.0037363.g004

Photo-induced disaggregation

Light has been observed to disaggregate TiO2 nanoparticles [35]. Bennett and coworkers have shown that the hydrodynamic diameter of TiO2 particles was reduced significantly after the samples were subjected to either sunlight or a Xenon lamp exposure and that the particle diameters returned to their equilibrium size in the dark. Consistent behavior is observed for CeO2 and ZnO nanoparticles, where light can disagglomerate the dispersions (Figure 5). Similar to the sonication experiments, there appears to be an agglomerate fraction that can be disagglomerated by light, but a hard aggregate core that cannot be disaggregated remains. Localized heating of the nanoparticle agglomerates, due to exposure to natural or artificial light, provided sufficient thermokinetic energy for the clusters to disagglomerate. In this case, h is 1.05 for TiO2, 1.14 for ZnO, and 1.09 for CeO2.

[Figure omitted. See PDF.]

Figure 5. Light induced disagglomeration of metal oxide nanoparticles.

https://doi.org/10.1371/journal.pone.0037363.g005

Effect of temperature

The hydrodynamic size remained unchanged when the temperature was increased from 25°C to 65°C; however, the size increased slowly when the dispersion was cooled back to 25°C (Figure 6A). The PDI increased from about 0.1 to 0.25 when the temperature was increased, and it decreased back to 0.1 when the sample was cooled back to 25°C (Figure 6B). The change in PDI indicates a broadened size distribution at elevated temperature, which is supported by the intensity distribution data as temperature increases (Figure 6C). At 25°C, the sample was monodisperse with a peak around 255 nm. A secondary peak emerged below 100 nm once the temperature reached 55°C, and this peak became more pronounced when the sample was heated up to 65°C. During the cooling phase, the secondary peak gradually shifted towards the main peak and disappeared after the sample temperature was below 45°C (Figure 6D). Since intensity-weighted distribution is strongly biased towards larger particles (intensity is proportional to the sixth power of particle size [9]), it is very likely that TiO2 clusters smaller than 100 nm have emerged before the temperature reached 45–55°C, but the light scattered by the smaller TiO2 clusters was too weak to make a distinguishable peak. Similar disagglomeration/reagglomeration behavior was also observed for both ZnO and CeO2 (Figure S4).

[Figure omitted. See PDF.]

Figure 6. Temperature effect on: (A) hydrodynamic diameter; (B) PDI; (C) intensity weighted size distribution with increasing temperature; (D) intensity weighted size distribution with decreasing temperature.

https://doi.org/10.1371/journal.pone.0037363.g006

The zeta-potential of the TiO2 nanoparticles decreased when the sample was heated, and increased when the temperature was reduced back to room temperature (Figure S5). A similar trend in zeta-potential with temperature was reported by another study [36]. The reduction in zeta-potential with increased temperature is not the cause of the disagglomeration, since a smaller zeta-potential would lead to less electrostatic repulsion and thus agglomeration. The higher thermal kinetic energy at elevated temperatures disrupted weak physical attractions inside the agglomerates but was unable to break the strongly bonded aggregates. The process of the disagglomeration from compact agglomerates to loose aggregates did not reduce the measured hydrodynamic size, since looser structures with the same mass yield greater hydrodynamic size. When the temperature was reduced, the weak attractions became sufficient again to hold aggregates together as agglomerates. However, the new agglomerates were not as compact as before since reagglomeration happened too fast for optimal configurations. This led to increased hydrodynamic size (Figure 6A). The schematic in Figure 7 illustrates this process.

[Figure omitted. See PDF.]

Figure 7. Schematic of temperature induced disagglomeration and reagglomeration.

https://doi.org/10.1371/journal.pone.0037363.g007

Interpretation via DLVO Modeling

Given the complex fractal nature of the agglomerate/aggregate mixture, it is likely that the van der Waals and electrostatic interactions between TiO2 agglomerates/aggregate mixture deviate from those of nearly spherical particles. Figure 8A presents the interaction-separation distance profiles predicted for single particles and fractal agglomerates. The DLVO theory predicts a rather deep primary minimum (∼1000 kT), a large energy barrier (∼30 kT), and a shallow secondary minimum (∼1.8 kT) for two particles interacting. In contrast, a comparable energy barrier (∼35 kT), a shallow secondary minimum (∼3.7 kT) and essentially no primary minimum are predicted for two fractal agglomerates. This explains the soft interactions between agglomerates that can be broken by heating. It is worth noting that the simulation result is not sensitive to the dF value.

[Figure omitted. See PDF.]

Figure 8. Simulation results considering (A) ionic strength = 0.1 mM; and (B) increasing ionic strength for fractal agglomerates. For particle-particle calculation, d = 200 nm. For fractal agglomerates calculation, dprimary = 54 nm, 20 primary particles in each aggregate, dF = 1.82.

https://doi.org/10.1371/journal.pone.0037363.g008

With increasing ionic strength, the electrostatic repulsion is screened and the energy barrier shrinks (Figure 8B). At the critical coagulation concentration (CCC), DLVO theory predicts irreversible aggregation for the particles at the primary minimum, while for fractal agglomerates additional reversible agglomeration will occur at the slightly increased secondary minimum. To evaluate this, diffusion-limited coagulation of TiO2 was induced at 15 mM NaCl (CCC of TiO2 [37]), with the hydrodynamic diameter growing over 1000 nm. The sample was then sonicated for 2 s with 7 W using the sonication probe. Sonication completely disagglomerated the clusters back to the initial size (Figure S6), suggesting the coagulation was reversible, supporting the prediction of the fractal agglomerate model.

Conclusion

Metal oxide nanoparticle synthesis may produce permanently sintered aggregate structures, which can further form reversible agglomerates when dispersed in aqueous media. This study shows that in open water these soft (weakly bonded) agglomerates can be disagglomerated by common environmental stimuli, such as exposure to sunlight or an increase in temperature from diurnal variations. Although not evaluated, it is likely that mechanical shocks may also result in temporary disagglomeration. The released aggregates can be much more mobile and bioavailable while the stimuli is present. Although in our experimental setting we observe reagglomeration once the stimuli are removed, in the environment it may be that the probability of interacting with another nanoparticle aggregate is much lower. Most toxicological studies are conducted at constant temperature and subdued light or under dark conditions. The effect of disagglomeration on toxicity has not been considered, or systematically evaluated. This phenomenon warrants attention since it is likely that these metal oxide nanoparticles will experience these natural stimuli during their transport in the environment.

Supporting Information

[Figure omitted. See PDF.]

Figure S1.

Sonication effect on ZnO (A) and CeO2 (B) nanoparticle dispersion. Each sonication cycle lasted for 2 s with power input 7 W.

https://doi.org/10.1371/journal.pone.0037363.s001

(TIF)

Figure S2.

Incorrect size measurement of TiO2 samples due to sonication induced temperature variation. Sonication duration = 1 min. Shaded area indicates incorrect measurements.

https://doi.org/10.1371/journal.pone.0037363.s002

(TIF)

Figure S3.

(A) Calculated sedimentation rate of TiO2 as a function of hydrodynamic size, dF = 2.1; (B) calculated sedimentation rate of TiO2 as a function of fractal dimension, hydrodynamic diameter = 200 nm. Equation adapted from Kajihara, 1971.1

https://doi.org/10.1371/journal.pone.0037363.s003

(TIF)

Figure S4.

Temperature effect on ZnO (A) and CeO2(B) hydrodynamic size.

https://doi.org/10.1371/journal.pone.0037363.s004

(TIF)

Figure S5.

Effect of temperature on the zeta-potential of TiO2.

https://doi.org/10.1371/journal.pone.0037363.s005

(TIF)

Figure S6.

Reversibility of the diffusion-limited coagulation of TiO2.agglomerates. At every point where size decreased back to below 400 nm, sample had been sonicated by sonication probe at 7 W for 2 s. Size grows again due to the high ionic strength.

https://doi.org/10.1371/journal.pone.0037363.s006

(TIF)

Table S1.

Parameters used in DLVO simulations.

https://doi.org/10.1371/journal.pone.0037363.s007

(DOC)

Table S2.

Characterization of metal oxide Nanoparticles.

https://doi.org/10.1371/journal.pone.0037363.s008

(DOC)

Acknowledgments

The authors would like to thank the lab assistance of Stephen Tjan, Adeel Lakhani, Danielle Vallone, and Ryan Tjan.

Author Contributions

Conceived and designed the experiments: DZ SWB AAK. Performed the experiments: DZ SWB. Analyzed the data: DZ AAK. Contributed reagents/materials/analysis tools: SWB. Wrote the paper: DZ AAK.

Citation: Zhou D, Bennett SW, Keller AA (2012) Increased Mobility of Metal Oxide Nanoparticles Due to Photo and Thermal Induced Disagglomeration. PLoS ONE7(5): e37363. https://doi.org/10.1371/journal.pone.0037363

References

1. Wiesner MR, Lowry GV, Alvarez P, Dionysiou D, Biswas P (2006) Assessing the risks of manufactured nanomaterials. Environ Sci Technol 40: 4336–4345.MR WiesnerGV LowryP. AlvarezD. DionysiouP. Biswas2006Assessing the risks of manufactured nanomaterials.Environ Sci Technol4043364345

2. Peralta-Videa JR, Zhao L, Lopez-Moreno ML, la Rosa de G, Hong J, et al. (2011) Nanomaterials and the environment: A review for the biennium 2008–2010. J Hazard Mater 186: 1–15.JR Peralta-VideaL. ZhaoML Lopez-MorenoG. la Rosa deJ. Hong2011Nanomaterials and the environment: A review for the biennium 2008–2010.J Hazard Mater186115

3. Godwin HA, Chopra K, Bradley KA, Cohen Y, Harthorn BH, et al. (2009) The University of California Center for the Environmental Implications of Nanotechnology. Environ Sci Technol 43: 6453–6457.HA GodwinK. ChopraKA BradleyY. CohenBH Harthorn2009The University of California Center for the Environmental Implications of Nanotechnology.Environ Sci Technol4364536457

4. Hotze EM, Phenrat T, Lowry GV (2010) Nanoparticle aggregation: challenges to understanding transport and reactivity in the environment. J Environ Qual 39: 1909–1924.EM HotzeT. PhenratGV Lowry2010Nanoparticle aggregation: challenges to understanding transport and reactivity in the environment.J Environ Qual3919091924

5. Petosa AR, Jaisi DP, Quevedo IR, Elimelech M, Tufenkji N (2010) Aggregation and deposition of engineered nanomaterials in aquatic environments: role of physicochemical interactions. Environ Sci Technol 44: 6532–6549.AR PetosaDP JaisiIR QuevedoM. ElimelechN. Tufenkji2010Aggregation and deposition of engineered nanomaterials in aquatic environments: role of physicochemical interactions.Environ Sci Technol4465326549

6. Chen K, Mylon S, Elimelech M (2006) Aggregation kinetics of alginate-coated hematite nanoparticles in monovalent and divalent electrolytes. Environ Sci Technol 40: 1516–1523.K. ChenS. MylonM. Elimelech2006Aggregation kinetics of alginate-coated hematite nanoparticles in monovalent and divalent electrolytes.Environ Sci Technol4015161523

7. Saleh NB, Pfefferle LD, Elimelech M (2008) Aggregation kinetics of multiwalled carbon nanotubes in aquatic systems: measurements and environmental implications. Environ Sci Technol 42: 7963–7969.NB SalehLD PfefferleM. Elimelech2008Aggregation kinetics of multiwalled carbon nanotubes in aquatic systems: measurements and environmental implications.Environ Sci Technol4279637969

8. Domingos RF, Tufenkji N, Wilkinson KJ (2009) Aggregation of titanium dioxide nanoparticles: role of a fulvic acid. Environ Sci Technol 43: 1282–1286.RF DomingosN. TufenkjiKJ Wilkinson2009Aggregation of titanium dioxide nanoparticles: role of a fulvic acid.Environ Sci Technol4312821286

9. French RA, Jacobson AR, Kim B, Isley SL, Penn RL, Baveye PC (2009) Influence of ionic strength, pH, and cation valence on aggregation kinetics of titanium dioxide nanoparticles. Environ Sci Technol 43: 1354–1359.RA FrenchAR JacobsonB. KimSL IsleyRL PennPC Baveye2009Influence of ionic strength, pH, and cation valence on aggregation kinetics of titanium dioxide nanoparticles.Environ Sci Technol4313541359

10. Zhou D, Keller AA (2010) Role of morphology in the aggregation kinetics of ZnO nanoparticles. Water Res 44: 2948–2956.D. ZhouAA Keller2010Role of morphology in the aggregation kinetics of ZnO nanoparticles.Water Res4429482956

11. Buettner KM, Rinciog CI, Mylon SE (2010) Aggregation kinetics of cerium oxide nanoparticles in monovalent and divalent electrolytes. Colloids Surf A 366: 74–79.KM BuettnerCI RinciogSE Mylon2010Aggregation kinetics of cerium oxide nanoparticles in monovalent and divalent electrolytes.Colloids Surf A3667479

12. Li X, Lenhart JJ, Walker HW (2010) Dissolution-accompanied aggregation kinetics of silver nanoparticles. Langmuir 26: 16690–16698.X. LiJJ LenhartHW Walker2010Dissolution-accompanied aggregation kinetics of silver nanoparticles.Langmuir261669016698

13. Chowdhury I, Hong Y, Walker SL (2010) Container to characterization: impacts of metal oxide handling, preparation, and solution chemistry on particle stability. Colloids Surf A 368: 91–95.I. ChowdhuryY. HongSL Walker2010Container to characterization: impacts of metal oxide handling, preparation, and solution chemistry on particle stability.Colloids Surf A3689195

14. Liu HH, Surawanvijit S, Rallo R, Orkoulas G, Cohen Y (2011) Analysis of nanoparticle agglomeration in aqueous suspensions via constant-number monte carlo simulation. Environ Sci Technol 45: 9284–9292.HH LiuS. SurawanvijitR. RalloG. OrkoulasY. Cohen2011Analysis of nanoparticle agglomeration in aqueous suspensions via constant-number monte carlo simulation.Environ Sci Technol4592849292

15. Bian SW, Mudunkotuwa IA, Rupasinghe T, Grassian VH (2011) Aggregation and dissolution of 4 nm ZnO nanoparticles in aqueous environments: influence of pH, ionic strength, size, and adsorption of humic acid. Langmuir 27: 6059–6068.SW BianIA MudunkotuwaT. RupasingheVH Grassian2011Aggregation and dissolution of 4 nm ZnO nanoparticles in aqueous environments: influence of pH, ionic strength, size, and adsorption of humic acid.Langmuir2760596068

16. Huynh KA, Chen KL (2011) Aggregation kinetics of citrate and polyvinylpyrrolidone coated silver nanoparticles in monovalent and divalent electrolyte solutions. Environ Sci Technol 45: 5564–5571.KA HuynhKL Chen2011Aggregation kinetics of citrate and polyvinylpyrrolidone coated silver nanoparticles in monovalent and divalent electrolyte solutions.Environ Sci Technol4555645571

17. Findlay A, Thompson D, Tipping E (1996) The aggregation of silica and haematite particles dispersed in natural water samples. Colloids Surf A 118: 97–105.A. FindlayD. ThompsonE. Tipping1996The aggregation of silica and haematite particles dispersed in natural water samples.Colloids Surf A11897105

18. Keller AA, Wang H, Zhou D, Lenihan HS, Cherr G, et al. (2010) Stability and aggregation of metal oxide nanoparticles in natural aqueous matrices. Environ Sci Technol 44: 1962–1967.AA KellerH. WangD. ZhouHS LenihanG. Cherr2010Stability and aggregation of metal oxide nanoparticles in natural aqueous matrices.Environ Sci Technol4419621967

19. Ji Z, Jin X, George S, Xia T, Meng H, et al. (2010) Dispersion and stability optimization of TiO2 nanoparticles in cell culture media. Environ Sci Technol 44: 7309–7314.Z. JiX. JinS. GeorgeT. XiaH. Meng2010Dispersion and stability optimization of TiO2 nanoparticles in cell culture media.Environ Sci Technol4473097314

20. Horst AM, Neal AC, Mielke RE, Sislian PR, Suh WH, et al. (2010) Dispersion of TiO2 nanoparticle agglomerates by pseudomonas aeruginosa. Appl Environ Microbiol 76: 7292–7298.AM HorstAC NealRE MielkePR SislianWH Suh2010Dispersion of TiO2 nanoparticle agglomerates by pseudomonas aeruginosa.Appl Environ Microbiol7672927298

21. Grass RN, Tsantilis S, Pratsinis SE (2006) Design of high-temperature, gas-phase synthesis of hard or soft TiO2 agglomerates. AIChE J 52: 1318–1325.RN GrassS. TsantilisSE Pratsinis2006Design of high-temperature, gas-phase synthesis of hard or soft TiO2 agglomerates.AIChE J5213181325

22. Teleki A, Wengeler R, Wengeler L, Nirschl H, Pratsinis SE (2008) Distinguishing between aggregates and agglomerates of flame-made TiO2 by high-pressure dispersion. Powder Technol 181: 292–300.A. TelekiR. WengelerL. WengelerH. NirschlSE Pratsinis2008Distinguishing between aggregates and agglomerates of flame-made TiO2 by high-pressure dispersion.Powder Technol181292300

23. Mandzy N, Grulke E, Druffel T (2005) Breakage of TiO2 agglomerates in electrostatically stabilized aqueous dispersions. Powder Technol 160: 121–126.N. MandzyE. GrulkeT. Druffel2005Breakage of TiO2 agglomerates in electrostatically stabilized aqueous dispersions.Powder Technol160121126

24. Hunter RJ (1981) Zeta potential in colloid science: principles and applications. RJ Hunter1981Zeta potential in colloid science: principles and applications.Academic Press, London. Academic Press, London.

25. Kim A, Berg J (2000) Fractal aggregation: Scaling of fractal dimension with stability ratio. Langmuir 16: 2101–2104.A. KimJ. Berg2000Fractal aggregation: Scaling of fractal dimension with stability ratio.Langmuir1621012104

26. Brasil A, Farias T, Carvalho M, Koylu U (2001) Numerical characterization of the morphology of aggregated particles. J Aerosol Sci 32: 489–508.A. BrasilT. FariasM. CarvalhoU. Koylu2001Numerical characterization of the morphology of aggregated particles.J Aerosol Sci32489508

27. Warren PB (1994) Hydrodynamics of fractal aggregates. IL Nuovo Cimento 16: 1231–1236.PB Warren1994Hydrodynamics of fractal aggregates.IL Nuovo Cimento1612311236

28. Elimelech M, Gregory J, Jia X, Williams R (1995) Particle deposition and aggregation: measurement, modelling and simulation. M. ElimelechJ. GregoryX. JiaR. Williams1995Particle deposition and aggregation: measurement, modelling and simulation.Elsevier. Elsevier.

29. Yotsumoto H, Yoon R (1993) Application of extended DLVO theory. 1. Stability of rutile suspensions. J Colloid Interface Sci 157: 426–433.H. YotsumotoR. Yoon1993Application of extended DLVO theory. 1. Stability of rutile suspensions.J Colloid Interface Sci157426433

30. Croll S (2002) DLVO theory applied to TiO2 pigments and other materials in latex paints. Prog Org Coat 44: 131–146.S. Croll2002DLVO theory applied to TiO2 pigments and other materials in latex paints.Prog Org Coat44131146

31. Kim T, Lee CH, Joo SW, Lee K (2008) Kinetics of gold nanoparticle aggregation: Experiments and modeling. J Colloid Interface Sci 318: 238–243.T. KimCH LeeSW JooK. Lee2008Kinetics of gold nanoparticle aggregation: Experiments and modeling.J Colloid Interface Sci318238243

32. Zhang Y, Chen Y, Westerhoff P, Crittenden J (2009) Impact of natural organic matter and divalent cations on the stability of aqueous nanoparticles. Water Res 43: 4249–4257.Y. ZhangY. ChenP. WesterhoffJ. Crittenden2009Impact of natural organic matter and divalent cations on the stability of aqueous nanoparticles.Water Res4342494257

33. Lattuada M, Wu H, Morbidelli M (2001) Estimation of fractal dimension of colloidal gels in the presence of multiple scattering. Phys Rev E 64: 061404.M. LattuadaH. WuM. Morbidelli2001Estimation of fractal dimension of colloidal gels in the presence of multiple scattering.Phys Rev E64061404

34. Kajihara M (1971) Settling velocity and porosity of large suspended particle. J Oceanogr 27: 158–162.M. Kajihara1971Settling velocity and porosity of large suspended particle.J Oceanogr27158162

35. Bennett SW, Zhou D, Montes MO, Mielke R, Keller AA (2012) Photoinduced disaggregation of nanoparticles enables transdermal penetration. PLoS ONE. SW BennettD. ZhouMO MontesR. MielkeAA Keller2012Photoinduced disaggregation of nanoparticles enables transdermal penetration.PLoS ONEIn press. In press.

36. Akratopulu KC, Kordulis C, Lycourghiotis A (1990) Effect of temperature on the point of zero charge and surface charge of TiO2. Faraday Trans 86: 3437–3440.KC AkratopuluC. KordulisA. Lycourghiotis1990Effect of temperature on the point of zero charge and surface charge of TiO2.Faraday Trans8634373440

37. Thio BJR, Zhou D, Keller AA (2011) Influence of natural organic matter on the aggregation and deposition of titanium dioxide nanoparticles. J Hazard Mater 189: 556–563.BJR ThioD. ZhouAA Keller2011Influence of natural organic matter on the aggregation and deposition of titanium dioxide nanoparticles.J Hazard Mater189556563