Elastic Constants of Quantum Solids by Path Integral Simulations

Abstract

Two methods are proposed to evaluate the second-order elastic constants of quantum mechanically treated solids. One method is based on path-integral simulations in the (NVT) ensemble using an estimator for elastic constants. The other method is based on simulations in the (NpT) ensemble exploiting the relationship between strain fluctuations and elastic constants. The strengths and weaknesses of the methods are discussed thoroughly. We show how one can reduce statistical and systematic errors associated with so-called primitive estimators. The methods are then applied to solid argon at atmospheric pressures and solid helium 3 (hcp, fcc, and bcc) under varying pressures. Good agreement with available experimental data on elastic constants is found for helium 3. Predictions are made for the thermal expectation value of the kinetic energy of solid helium 3.

9 pages doublecolumn, 6 figures, submitted to PRB

References in corpus (2)

• The Debye-Waller Factor in solid 3He and 4He
• Simulation of Material Properties Below the Debye Temperature: A Path-Integral Molecular Dynamics Case Study of Quartz

Cited by in corpus (7)

• Molecular Dynamics Simulations
• Simulation and understanding of quantum crystals
• Thermal conductivity of polymers: A simple matter where complexity matters
• On quantum effects near the liquid-vapor transition in helium
• Elastic constants of hcp ${}^4$He: Path-integral Monte Carlo results versus experiment
• Analytic elastic constants in molecular calculations: Finite strain, non-affine displacements, and many-body interatomic potentials
• Computing the thermal transport coefficient of neutral amorphous polymers using exact vibrational density of states: Comparison with experiments