Kinetics of inherent processes counteracting crystallization in supercooled monatomic liquid

• Published in: Journal of physics. Condensed matter Vol. 34; no. 45; pp. 454002 - 454009
• Main Authors: Galimzyanov, B N, Yarullin, D T, Mokshin, A V
• Format: Journal Article
• Language: English
• Published: IOP Publishing 09.11.2022
• Subjects: crystal growth; crystal nucleation; crystallization kinetics; molecular dynamics; phase transitions
• ISSN: 0953-8984; 1361-648X
• DOI: 10.1088/1361-648X/ac8fd1

Abstract Crystallization of supercooled liquids is mainly determined by two competing processes associated with the transition of particles (atoms) from liquid phase to crystalline one and, vice versa, with the return of particles from crystalline phase to liquid one. The quantitative characteristics of these processes are the so-called attachment rate g + and the detachment rate g − , which determine how particles change their belonging from one phase to another. In the present study, a correspondence rule between the rates g + and g − as functions of the size N of growing crystalline nuclei is defined for the first time. In contrast to the well-known detailed balance condition, which relates g + ( N ) and g − ( N ) at N = n c (where n c is the critical nucleus size) and is satisfied only at the beginning of the nucleation regime, the found correspondence rule is fulfilled at all the main stages of crystallization kinetics (crystal nucleation, growth and coalescence). On the example of crystallizing supercooled Lennard–Jones liquid, the rate g − was calculated for the first time at different supercooling levels and for the wide range of nucleus sizes N ∈ [ n c ; 40 n c ] . It was found that for the whole range of nucleus sizes, the detachment rate g − is only ≈ 2 % less than the attachment rate g + . This is direct evidence that the role of the processes that counteract crystallization remains significant at all the stages of crystallization. Based on the obtained results, a kinetic equation was formulated for the time-dependent distribution function of the nucleus sizes, that is an alternative to the well-known kinetic Becker–Döring–Zeldovich–Frenkel equation.