Efficient simulation of second-order phase transitions in quantum anharmonic materials
Andrea Baldanza, Lorenzo Monacelli
Abstract
When a crystal undergoes a second-order structural phase transition, such as in ferroelectrics, Peierls, and charge-density waves, the diverging fluctuations of the order parameter lead to the break- down of the standard phonon quasiparticle picture. Simulating these highly anharmonic regimes is notoriously challenging, as methods such as molecular dynamics suffer from a critical slowdown near the transition point, while the harmonic approximation fails dramatically at saddle points of the energy landscape stabilized by quantum or thermal ionic fluctuations. This work introduces a new approach, based on the variational free-energy principle, to predict critical long-range behavior and dynamical spectra in strongly anharmonic systems, even when quantum ionic fluctuations dominate. The proposed framework builds upon the stochastic self- consistent harmonic approximation but reduces its computational scaling with the number of atoms, N, from O(N6) to O(N2) and the memory requirement from O(N4) to O(N). We benchmark the method on the prototypical lead-free metal-halide perovskite CsSnI3, a promising candidate for photovoltaic engineering, simulating its phase stability and Raman spectrum near the phase transition, where the breakdown of the quasiparticle picture becomes evident. We demonstrate the effectiveness of the method by computing the full free-energy Hessian and the critical temperature in a supercell with 1080 atoms. Such simulations would have required tens of thousands of years with the legacy approach; it is now feasible in a few hours on consumer hardware.
Create a lesson
Related papers
Divergence between long- and short-wavelength magnon damping in spinel ferrites
Christopher T. Parzyck, Octave Duros, Hari Paudyal et al.
An Atlas and Design Rules for Single- and Dual-Atom Alloys
Fabian Berger, Yicheng Wang, E. Charles H. Sykes et al.
Epitaxial inversion of spontaneous polarization in ε-Ga2O3
Yan Wang, Zhigao Xie, Weihua Tang et al.
Gauge-including neural-network quantum Monte Carlo for molecules in magnetic fields
Chengye Lü, Weizhong Fu, Xin-gao Gong et al.
Photoresponse properties of single-crystalline thick film based on high-entropy topological insulator (Bi3/4Sb1/4)2(Te2/5Se2/5S1/5)3
Alexei Vasilev, Marina Zhezhu, Oleg Ivanov
Adaptive Substrate Support Based on Thin-Film Piezoelectric Actuators
Ertuğ Şimşek, Bas Jansen, Marcelo Ackermann et al.