Computation of elastic constants of solids using molecular simulation: comparison of constant volume and constant pressure ensemble methods

We compute the elastic stiffness tensor of fcc argon at 60 K and 1 bar using molecular simulation tools. Three different methods are investigated: explicit deformations of the simulation box, strain fluctuations at constant pressure and stress fluctuations at constant volume. Statistical ensemble sampling is done using molecular dynamics and Monte Carlo simulations. We observe a good agreement between the different methods and sampling algorithms excepted with molecular dynamics simulations in the (NpT) ensemble. There, we notice a strong dependence of the computed elastic constants with the barostat parameter, whereas molecular dynamics simulations in the (NVT) ensemble are not affected by the thermostat parameter.