Temperature effect on the phonon dispersion stability of zirconium by machine learning driven atomistic simulations

It is well known that conventional harmonic lattice dynamics cannot be applied to dynamically unstable crystals at 0 K, such as the high-temperature body centered cubic (bcc) phase of crystalline Zr. Predicting phonon spectra at finite temperature requires the calculation of force constants to the third, fourth, and even higher orders; however, it remains challenging to determine at which order the Taylor expansion of the potential-energy surface for different materials should be cut off. Molecular dynamics, on the other hand, intrinsically includes arbitrary orders of phonon anharmonicity; however, its accuracy is severely limited by the empirical potential field used. Recently, machine learning algorithms have emerged as promising tools to build accurate potentials for molecular dynamics simulation. In this paper, we approach the problem of predicting phonon dispersion at finite temperature by performing molecular dynamics simulations with machine learning driven potential fields. We developed Gaussian approximation potential models for both the hexagonal close packed (hcp) phase and the bcc phase of zirconium crystals. The developed potential field is first validated with static properties including energy-volume relationship, elastic constants, and phonon dispersions at 0 K. Molecular dynamics simulations are then performed to stochastically sample the potential-energy surface and to calculate the phonon dispersion at elevated temperatures. The phonon renormalization in bcc Zr is successfully captured by the molecular dynamics simulation at 1188 K. The instability of the bcc structure is found to originate from the double-well shape of the potential-energy surface where the local maximum is located in an unstable equilibrium position. The stabilization of the bcc phase at high temperature is due to the dynamical average of the low-symmetry minima of the double well due to atomic vibrations.