Молекулярно-динамическое моделирование стеклования переохлажденного расплава алюминия



Vitrification of aluminum melt. A molecular dynamics study

Molecular dynamics study of transition from liquid aluminum into amorphous solid state is carried out. Aluminum is described by a form of embedded atom method potential. The system was first equilibrated at T = 1500 K, which is above the melting temperature of aluminum. Then the velocities of atoms were gradually rescaled until the system reached T = 300 K. Influence of cooling rate on the final state of the system is studied. MD simulations were held for cooling rates from 1010 K/s to 1013 K/s. At cooling rates below 1010 K/s the melt eventually crystallizes. At higher cooling rates the melt freezes into an amorphous structure. Structural analysis of amorphous aluminum based on radial distribution function of atoms is performed. Splitting of the second peak is used as a criterion for the amorphous structure. The results show the vitrifying temperature and final pressure in glass depends on cooling rate. At cooling rate higher than 1013 K/s the final amorphous state is nonequilibrium and relaxation time is approximately 100 ps. Equation of state hysteresis takes place during cooling and heating. Phase transition from amorphous to solid state can take place at higher temperatures (800 − 900 K). The results are compared with the common data.

molecular modeling, material science, glass transition, crystallization, metglassy, metastable states

Лада Николаевна Колотова, Генри Эдгарович Норман, Василий Вячеславович Писарев

Том 14, выпуск 2, 2013 год



В данной работе с помощью метода молекулярной динамики (МД) рассматривается влияние скорости охлаждения на положение точки стеклования расплава алюминия. Для МД моделирования алюминия используется потенциал погруженного атома. Рассматривается режим охлаждения при постоянном объеме. При этом используется изменение парной корреляционной функции, как критерий, различающий стекло и переохлажденный расплав. Определены диапазоны скоростей охлаждения, при которых наблюдается стеклование расплава, критическая скорость охлаждения ~ 1010 K/с. Показано, что температура стеклования повышается с увеличением скорости охлаждения, но при скоростях охлаждения ниже 1012 K/с остается практически постоянной. Аморфное состояние является неравновесным и время релаксации давления при 300 K составляет 100 пс. Аморфное состояние может переходить в кристаллическое при температурах вблизи точки стеклования. При низких температурах кристаллизация не происходит в связи с низкой подвижностью атомов.

молекулярное моделирование, материаловедение, стеклование, кристаллизация, метглассы, метастабильные состояния

Лада Николаевна Колотова, Генри Эдгарович Норман, Василий Вячеславович Писарев

Том 14, выпуск 2, 2013 год



1. Хоник В. А. // СОЖ. 2001. № 3. С.95-103.
2. Черноуцан А. И. // СОЖ. 2001. №3. С.103-110.
3. Золотухин И. В. // СОЖ. 1997. №4. С.73-79.
4. Fedorchenko A. I., Chernov A.A. // Russ. J. Eng. Termophys.
2000. V. 10. N. 3. Pp. 201-206.
5. Fedorchenko A. I., Chernov A.A. // Int. J. Heat Mass Transf.
2003. V. 46. Pp. 921-929.
6. Khonik V.A., Kitagawa K., Morii H. // J. Appl. Phys. 2000. V.
87. No. 12. Pp. 8440-8443.
7. Mendelev M. I., Schmalian J., Wang C., Morris J., Ho K. //
Phys.Rev.B. 2006. V. 74. No. 10. P. 104206.
8. Щульц М. М. // СОЖ. 1996. №3. С.49-56.
9. Полухин В. А., Ватолин Н. А. Моделирование аморфных
металлов. 1985. М.:Наука.
10. Bengtzelius U., Gotze W., Sjolander A. // J. Phys. C. 1984. V.
17. Pp. 5915–5934.
11. V. N. Ryzhov, E. E. Tareyeva, and Yu. D. Fomin // J. TMP.
2011. V. 167(2). Pp. 645–653.
12. Daw M. S., Baskes M. I. // Phys. Rev. Lett. 1983. V. 50. P.
1285.
13. Liu X.Y., Wei X., Foiles S.M., Adams J.B. // Appl. Phys. Lett.
2006. V. 72. Pp. 1578-1580.
14. Fomin Yu. D., Gribova N. V., Ryzhov V. N.// J. Chem. Phys.
2008. V. 129. P. 064512
15. Белащенко Д. К. // ЖФХ. 2008. Т. 82. № 3. С. 443-455.
16. Plimpton S.G. // J. Comp. Phys. 1995. V. 117. Pp. 1-19.