Abstract
This work presents modelling for the aggregation process of metal nanoparticles following a theory proposed in literature. In this theory, metal atoms aggregate into particles of bigger size due to Van der Waal’s forces of attraction. Then, owing to the electrostatic forces of repulsion, the particles eventually stop aggregating and become stabilized. Based on this mechanistic description, we developed a model for the aggregation process. Because this process often occurs along with other processes, such as growth, we employed as case study the synthesis of gold nanoparticles by the citrate synthesis method for conditions where the aggregation process is decoupled from the growth process. Using this model, we calculated the seed particle sizes and compared them with the values previously reported in literature. Furthermore, we calculated the final particle sizes and compared them with experimental data. The results show excellent agreement.
Similar content being viewed by others
Abbreviations
- Symbol :
-
Meaning Units
- C :
-
Concentration mol/m3
- C C :
-
Citrate mol/m3
- CT, \( {C}_{Au{Cl_4}^{-}} \) :
-
Concentration of tetrachloroauric ion mol/m3
- C Ct :
-
Concentration of all citrate species mol/m3
- \( {C}_{H^{+}} \) :
-
Concentration of H+ ions mol/m3
- \( {C}_{OH^{-}} \) :
-
Concentration of OH− ions mol/m3
- C Au :
-
Concentration of gold in the particle phase mol/m3
- C Pr1 :
-
Concentration of all other products from the reduction step, lumped together mol/m3
- C Pr2 :
-
Concentration of all other products from the growth step, lumped together mol/m3
- D s :
-
Seed correlation parameter −
- D z :
-
Stability gradient correlation parameter −
- e :
-
The charge on an electron, whose value C
- E A :
-
Energy due to Van der Waal’s force of attraction J
- E R :
-
Energy due to the charge repulsion J
- E T :
-
the sum of the particles’ interaction energy due to Van der Waal’s force of attraction and that due the charge repulsion J
- E agg :
-
the energy barrier to particle aggregation J
- E max :
-
The maximum particles’ interactive energy attainable J
- f(s):
-
number of particles per particle-length per total volume of fluid-particle mixture
1/(m3. m)
- F z :
-
Stability gradient correlation parameter −
- G s :
-
linear growth rate m/s
- k n :
-
Rate constant for the nucleation step (m3)4/(mol5. s)
- k h :
-
Rate constant for the growth step m/s
- k r :
-
Rate constant for the reduction step [m3/mol]1.851/s
- k p :
-
Rate constant for the passivation step m3/(mol. s)
- k g :
-
Rate constant for the growth step m4/(mol. s)
- K B :
-
Boltzmann constant J/K
- m a :
-
Particles area shape factor −
- m v :
-
Particles volume shape factor −
- N av :
-
Avogadro ’ s number = 6.02e23−
- p i :
-
The number concentration of an ion in the bulk of the solution 1/m3
- R :
-
Universal gas constant 8.31J/(mol. K)
- s :
-
Size m
- s 0 :
-
Nucleus diameter (size of gold atom, i.e. 0.272 nm in the aggregation model)
m
- s s :
-
Seed diameter m
- s m :
-
Mean diameter with time m
- s f :
-
Final mean diameter m
- S :
-
Selectivity of the reduction step over the passivation step −
- t :
-
Time s
- t s :
-
Synthesis time s
- T :
-
Temperature K
- V :
-
volume of synthesis solution m3
- W :
-
Stability factor −
- \( {\overset{\sim }{\omega}}_A\left(\overline{s},\hat{s},t\right) \) :
-
Aggregation kernel m3/s
- x :
-
Separation of particles m
- y x :
-
relative mole fraction of CtH2− at the quasi-equilibrium pH −
- y y :
-
relative mole fraction of CtH2− at the quasi-equilibrium pH −
- Y i :
-
Molar mass of species i kg/mol
- z i :
-
The charge on the ion C
- Z :
-
The numerator of the stability gradient J
- ϵ 0 :
-
The permittivity of free space F/m
- ϵ c :
-
The dielectric constant of the solution −
- κ :
-
the Debye-Huckel parameter 1/m
- μ :
-
Fluid viscosity kg/(m. s)
- ρ :
-
Molar density of gold mol/m3
- ψ x :
-
Electrostatic charge potential V
- τ :
-
Characteristic time s
- τ p :
-
Reaction time for the passivation s
- τ g :
-
Time for the growth step only in the citrate method s
References
Agunloye E, Gavriilidis A, Mazzei L (2017) A mathematical investigation of the Turkevich organizer theory in the citrate method for the synthesis of gold nanoparticles. Chem Eng Sci. https://doi.org/10.1016/j.ces.2017.07.032
Agunloye E, Panariello L, Gavriilidis A, Mazzei L (2018) A model for the formation of gold nanoparticles in the citrate method, Chem Eng Sci, 173, 275–286. https://doi.org/10.1016/j.ces.2018.06.046
Biggs S, Chow MK, Zukoski CF, Grieser F (1993) The role of colloidal stability in the formation of gold sols. J Colloid Interface Sci 160:511–513
Bogush GH, Zukoski CF IV (1991) Uniform silica particle precipitation: an aggregative growth model. J Colloid Interface Sci 142:19–34. https://doi.org/10.1016/0021-9797(91)90030-C
Cordero B, Gómez V, Ana EP-P, Marc R, Jorge E, Eduard C, Flavia B, Santiago A (2008) Covalent radii revisited. Dalton Trans 21:2832–2838. https://doi.org/10.1039/b801115j
Derjaguin B, Landau L (1941) Theory of the stability of strongly charged lyophobic sols and of the adhesion of strongly charged particles in solutions of electrolytes. Prog Surf Sci 43:30–59. https://doi.org/10.1016/0079-6816(93)90013-L
Elimelech M, Gregory J, Jia X, Williams R (1995) Particle deposition and aggregation - measurement, modelling and simulation. Elsevier
Finke RG, Watzky MA (1997) Nanocluster size-control and “magic number investigations”. Chem Mater 9:3083–3095. https://doi.org/10.1021/cm9704387
Fuchs NA (1964) The mechanics of aerosols, Pergamon
Israelachvili JN (2011) Intermolecular and surface forces, 3rd edn. Academic Press, Burlinton
Ji XH, Song XN, Li J, Bai YB, Yang WS, Peng XG (2007) Size control of gold nanocrystals in citrate reduction: the third role of citrate. J Am Chem Soc 129:13939–13948
Kumar S, Kumar R, Gandhi KS (2007) Modeling of formation of gold nanoparticles by citrate method. Ind Eng Chem Res 46:3128–3136
LaMer VK, Dinegar RH (1950) Theory, production and mechanism of formation of monodispersed hydrosols. J Am Chem Soc 72:4847–4854. https://doi.org/10.1021/ja01167a001
Liveri VT (2006) Controlled synthesis of nanoparticles in microheterogeneous systems, Nanostructure Science and Technology, Springer Science
Marchisio DL, Fox RO (2013) Computational models for polydisperse particulate and multiphase systems. Cambridge University Press
Polte J (2015) Fundamental growth principles of colloidal metal nanoparticles—a new perspective. CrystEngComm 17(36):6809–6830
Ramkrishna D (2000) Population balances. Academic Press
Reerink H, Overbeek JTG (1954) The rate of coagulation as a measure of the stability of silver iodide sols. Discuss Faraday Soc 18:74–84. https://doi.org/10.1039/DF9541800074
Thünemann AF, Kraehnert R 2010. Mechanism of Gold Nanoparticle Formation in the Classical Citrate Synthesis Method Derived from Coupled In Situ XANES and SAXS Evaluation. J Am Chem Soc 132:1296–1301, https://doi.org/10.1021/ja906506j
Turkevich J, Stevenson P, Hillier J (1951) A study of the nucleation and growth process in the synthesis of colloidal gold. Discuss Faraday Soc 11, 55
Verwey EJW, Overbeek JTG (1948) Theory of the stability of lyophobic colloids. Elsevier, Amsterdam
Wuithschick M, Witte S, Kettemann F, Rademann K, Polte J (2015) Turkevich in new robes: key questions answered for the most common gold nanoparticle synthesis. Phys Chem Chem Phys 17:19895–19900. https://doi.org/10.1021/acsnano.5b01579
Zabetakis K, Ghann WE, Kumar S, Daniel MC (2012) Effect of high gold salt concentration on the size and polydispersity of gold nanoparticles prepared by and extended Turkevich-Frens method. Gold Bull 45:203–211. https://doi.org/10.1007/s13404-012-0069-2
Acknowledgements
The authors are grateful to Dr. Luca Mazzei, Associate Professor in the Department of Chemical Engineering, University College London, and Dr. Michael Wulkow, managing director of Computing in Technology GmbH (CiT), the company that developed the numerical code Parsival, for their support and stimulating discussions. Emmanuel holds both Luca and Michael in very high esteem. Emmanuel Agunloye would also like to thank the Nigerian government for funding his PhD programme via the Petroleum Technology Development Fund and the National University Commission.
Author information
Authors and Affiliations
Corresponding author
Ethics declarations
Conflict of interest
The authors declare that they have no conflict of interest.
Additional information
Research highlights
• New mathematical modelling for the aggregation of metal nanoparticles is presented.
• The modelling is based on the theory proposed by Polte (2015).
• This theory accounts for the interplay of attractive and repulsive forces among metal nanoparticles.
• The evolution of predicted nanoparticle size distributions closely agrees with those measured experimentally.
Publisher’s note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Electronic supplementary material
ESM 1
(DOCX 139 kb)
Rights and permissions
About this article
Cite this article
Agunloye, E., Usman, M. Modelling of the aggregation process using the citrate synthesis of gold as case study. J Nanopart Res 22, 183 (2020). https://doi.org/10.1007/s11051-020-04871-1
Received:
Accepted:
Published:
DOI: https://doi.org/10.1007/s11051-020-04871-1