Study of the thermodynamic characteristics of gold nanoparticles in the processes of melting and crystallization by the molecular dynamics method

The article presents the results of computer modeling of melting and crystallization processes of spherical gold nanoparticles. The molecular dynamics method was used to analyze the thermodynamic characteristics (temperature, heat, entropy of melting, and crystallization) of nanoparticles for different heating and cooling rates. Such studies allow choosing the most suitable temperature ranges for the formation of nanocrystalline structures and predicting their sizes. Reducing the size of gold nanoparticles (less than 100 nm) leads to a significant increase in the ratio of the surface area to the volume of particles, as a result of which the physical and chemical characteristics of the material change significantly compared to materials of the same material in bulk form. Interest in gold nanoparticles is due to enhanced photoemission, high electrical and thermal conductivity, and increased catalytic activity of the surface. Gold nanoparticles have strong optical absorption and scattering properties in the visible region of the spectrum due to surface plasmon oscillations of free electrons. In the known studies from the literature, it was found that with an increase in the size of nanoparticles, the hysteresis between the melting and crystallization temperatures increases, while in theory the macroscopic melting and crystallization temperatures should be the same. The novelty of the study presented in this paper is to reveal a previously unobserved tendency for macroscopic temperatures, heat and entropy of melting and crystallization to converge with a decrease in the heating and cooling rates. The classical molecular dynamics method was used to study the thermodynamic properties of gold nanoparticles. The subject of modeling was gold nanoparticles of various sizes, spherical in shape with a face-centered cubic lattice. In the modeling process, the interatomic potential was used corresponding to the embedded atom method which was developed for gold, using an improved force matching methodology. Heating and cooling of nanoparticles were modeled at temperature change rates of 0.1 TK/s, 1 TK/s, 3 TK/s. By analyzing the relationship between the potential energy of gold nanoparticles and temperature, the dependences of the melting and crystallization temperatures of nanoparticles on their size were revealed. A relationship was established between the size of nanoparticles, heat, entropy of melting and crystallization at different heating and cooling rates. It was shown that with a decrease in the heating and cooling rates from 3 TK/s to 0.1 TK/s, there is a convergence of the macroscopic values of the melting and crystallization temperatures (a decrease in the difference from 467 K to 158 K), macroscopic values of the heat of melting and crystallization (from 4.24 kJ/mol to 0.67 kJ/mol), entropy of melting and crystallization (from 1.99 J/(mol·K) to 0.16 J/(mol·K)). It was assumed that this is due to a decrease in the proportion of the formed nanostructures other than face-centered cubic ones. Prediction of temperature regimes of melting and crystallization of gold nanoparticles allows controlling phase transitions in the production of nanocrystals with specified properties. This phenomenon can find application in microelectronics for the formation of thin films with a high degree of homogeneity and in catalysis for the formation of nanoparticles with the required structures and properties. The obtained data can be used to verify the results of real experiments on phase transitions and to adjust molecular dynamics models.

1. Skripov V.P., Koverda V.P. Spontaneous Crystallization of Supercooled Liquids: Nucleation of Crystals in Liquids and Amorphous Solids. Moscow, Nauka Publ.,1984. 232 p. (in Russian)

2. Sdobnyakov N.Yu., Komarov P.V., Sokolov D.N., Samsonov V.M. Study of the thermodynamic characteristics of gold nanoclusters using a Gupta many-body potential, The Physics of Metals and Metallography, 2011, vol. 111, no. 1, pp. 13–20. https://doi.org/10.1134/s0031918X11010121

3. Samsonov V.M., Kharechkin S.S., Gafner S.L., Redel’ L.V., Gafner Yu.Ya. Molecular dynamics study of the melting and crystallization of nanoparticles, Crystallography Reports, 2009, vol. 54, no. 3, pp. 526–531. https://doi.org/10.1134/S1063774509030250

4. Sdobnyakov N.Y., Sokolov D.N., Bazulev A.N., Samsonov V.M., Zykov T.Y., Antonov A.S. Relation between the size dependences of the melting and crystallization temperatures of metallic nanoparticles. Russian Metallurgy (Metally), 2013, vol. 2013, no. 2, pp. 100–105. https://doi.org/10.1134/S0036029513020110

5. Gupta R.P. Lattice relaxation at a metal surface. Physical Review B, 1981, vol. 23, no. 12, pp. 6265–6270. https://doi.org/10.1103/PhysRevB.23.6265

6. Cleri F., Rosato V. Tight-binding potentials for transition metals and alloys. Physical Review B, 1993, vol. 48, no. 1, pp. 22–33. https://doi.org/10.1103/PhysRevB.48.22

7. Samsonov V.M., Vasilyev S.A., Talyzin I.V., Ryzhkov Y.A. On reasons for the hysteresis of melting and crystallization of nanoparticles. JETP Letters, 2016, vol. 103, no. 2, pp. 97–99. https://doi.org/10.1134/S0021364016020119

8. Grochola G., Russo S.P., Snook I.K. On fitting a gold embedded atom method potential using the force matching method. The Journal of Chemical Physics, 2005, vol. 123, no. 20, pp. 204719. https://doi.org/10.1063/1.2124667

9. Tsydypov D.G., Nomoev A.V. Calculation of the dependence of the melting point temperature of a silver nanoparticle on the nanoparticle size using different embedded atom potentials. BSU bulletin. Chemistry. Physics, 2022, no. 2-3, pp. 22–30. (in Russian). https://doi.org/10.1063/1.2124667

10. Tsydypov D.G., Nomoev A.V., Garmaev B.Z. Determination of thermophysical properties and the radial distribution function of a silver nanoparticle using numerical methods. Proceedings of the Russian Higher School Academy of Sciences, 2023, no. 3 (60), pp. 20–31. (in Russian). https://doi.org/10.17212/1727-2769-2023-3-20-31

11. Thomson W. LX. On the equilibrium of vapour at a curved surface of liquid. The London, Edinburgh, and Dublin Philosophical Magazine and Journal of Science, 1871, vol. 42, no. 282, pp. 448–452. https://doi.org/10.1080/14786447108640606

12. Babichev A.P. Physical Quantities. Handbook. Moscow, Energoatomizdat Publ., 1991. 1232 p. (in Russian)

13. Sokolov D.N., Sdobnyakov N.Yu., Komarov P.V. The calculation of the size dependencies of the heat of fusion for metal nanoparticles. Physical and Chemical Aspects of the Study of Custers, Nanostructures and Nanomaterials, 2011, no. 3, pp. 229–138. (in Russian)

14. Sdobnyakov N.Yu., Komarov P.V., Kolosov A.Yu., Novozhilov N.V., Sokolov D.N., Kulpin D.A. The calculation of the size dependencies of the heat of fusion for metal nanoparticles. Condensed Matter and Interphases, 2013, vol. 15, no. 3, pp. 334–344. (in Russian)

15. Samsonov V.M., Sdobnyakov N.Yu., Vasilyev S.A., Sokolov D.N. On the size dependence of the heats of melting of metal nanoclusters. Bulletin of the Russian Academy of Sciences: Physics, 2016, vol. 80, no. 5, pp. 494–496. https://doi.org/10.3103/S1062873816050166

16. Samsonov V.M., Vasilyev S.A., Bembel A.G., Samsonov T.E., Skopich V.L. Molecular dynamics investigation of the size dependence of the heat of melting of metal nanoclusters. Physics of the Solid State, 2014, vol. 56, no. 12, pp. 2369–2373. https://doi.org/10.1134/S1063783414120270

17. Sdobnyakov N.Yu., Sokolov D.N., Myasnichenko V.S., Bazulev A.N. Cаlculation of the heat of fusion and the heat of freezing size dependencies for metal nanoparticles. Physical and Chemical Aspects of the Study of Custers, Nanostructures and Nanomaterials, 2014, no. 6, pp. 342–348. (in Russian)

18. Sdobnyakov N.Yu., Veselov A.D., Ershov P.M., Sokolov D.N., Samsonov V.M., Vasilyev S.A., Myasnichenko V.S. Size dependence of the entropies of melting and crystallisation of metal nanoparticles. Computational Materials Science, 2018, vol. 153, pp. 153–158. https://doi.org/10.1016/j.commatsci.2018.06.037