Nonaffine lattice dynamics with the Ewald method reveals strongly nonaffine elasticity of α-quartz

J Chem Phys. 2019 Dec 14;151(22):224509. doi: 10.1063/1.5129025.

Abstract

A lattice dynamical formalism based on nonaffine response theory is derived for noncentrosymmetric crystals, accounting for long-range interatomic interactions using the Ewald method. The framework takes equilibrated static configurations as input to compute the elastic constants in excellent agreement with both experimental data and calculations under strain. Besides this methodological improvement, which enables faster evaluation of elastic constants without the need of explicitly simulating the deformation process, the framework provides insights into the nonaffine contribution to the elastic constants of α-quartz. It turns out that, due to the noncentrosymmetric lattice structure, the nonaffine (softening) correction to the elastic constants is very large, such that the overall elastic constants are at least 3-4 times smaller than the affine Born-Huang estimate.