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 Cij. 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 3He. Predictions are also made for the thermal expectation value of the kinetic energy of solid 3He.
- Publication:
-
Physical Review B
- Pub Date:
- June 2001
- DOI:
- 10.1103/PhysRevB.63.224108
- arXiv:
- arXiv:cond-mat/0009379
- Bibcode:
- 2001PhRvB..63v4108S
- Keywords:
-
- 67.80.-s;
- 62.20.Dc;
- 02.70.Ns;
- 05.30.-d;
- Solid helium and related quantum crystals;
- Elasticity elastic constants;
- Molecular dynamics and particle methods;
- Quantum statistical mechanics;
- Condensed Matter - Materials Science;
- Condensed Matter - Statistical Mechanics
- E-Print:
- 9 pages doublecolumn, 6 figures, submitted to PRB