Jakse Noel, Sandberg Johannes, Granz Leon F, Saliou Anthony, Jarry Philippe, Devijver Emilie, Voigtmann Thomas, Horbach Jürgen, Meyer Andreas
Université Grenoble Alpes, CNRS, Grenoble INP, SIMaP, F-38000 Grenoble, France.
Institut für Materialphysik im Weltraum, Deutsches Zentrum für Luft- und Raumfahrt (DLR), 51170 Köln, Germany.
J Phys Condens Matter. 2022 Nov 15;51(3). doi: 10.1088/1361-648X/ac9d7d.
In studying solidification process by simulations on the atomic scale, the modeling of crystal nucleation or amorphization requires the construction of interatomic interactions that are able to reproduce the properties of both the solid and the liquid states. Taking into account rare nucleation events or structural relaxation under deep undercooling conditions requires much larger length scales and longer time scales than those achievable bymolecular dynamics (AIMD). This problem is addressed by means of classical molecular dynamics simulations using a well established high dimensional neural network potential trained on a set of configurations generated by AIMD relevant for solidification phenomena. Our dataset contains various crystalline structures and liquid states at different pressures, including their time fluctuations in a wide range of temperatures. Applied to elemental aluminium, the resulting potential is shown to be efficient to reproduce the basic structural, dynamics and thermodynamic quantities in the liquid and undercooled states. Early stages of crystallization are further investigated on a much larger scale with one million atoms, allowing us to unravel features of the homogeneous nucleation mechanisms in the fcc phase at ambient pressure as well as in the bcc phase at high pressure with unprecedented accuracy close to theone. In both cases, a single step nucleation process is observed.
在通过原子尺度模拟研究凝固过程时,晶体成核或非晶化的建模需要构建能够再现固态和液态性质的原子间相互作用。考虑到深过冷条件下的罕见成核事件或结构弛豫,所需的长度尺度和时间尺度比分子动力学(AIMD)所能达到的要大得多。通过经典分子动力学模拟解决了这个问题,该模拟使用了一个成熟的高维神经网络势,该势是在一组由与凝固现象相关的AIMD生成的构型上训练得到的。我们的数据集包含不同压力下的各种晶体结构和液态,包括它们在很宽温度范围内的时间波动。应用于元素铝时,结果表明所得势能够有效地再现液态和过冷态的基本结构、动力学和热力学量。在更大的规模上,用一百万个原子进一步研究了结晶的早期阶段,这使我们能够以前所未有的精度揭示常压下fcc相以及高压下bcc相的均匀成核机制的特征。在这两种情况下,都观察到了单步成核过程。