First-Principles Modeling of Quantum Nuclear Effects and Atomic Interactions in Solid 4He at High Pressure
arXiv:1409.4945 · doi:10.1103/PhysRevB.91.024103
Abstract
We present a first-principles computational study of solid 4He at T=0K and pressures up to 160GPa. Our computational strategy consists in using van der Waals density functional theory (DFT-vdW) to describe the electronic degrees of freedom in this material, and the diffusion Monte Carlo (DMC) method to solve the Schrodinger equation describing the behavior of the quantum nuclei. For this, we construct an analytical interaction function based on the pairwise Aziz potential that closely matches the volume variation of the cohesive energy calculated with DFT-vdW in dense helium. Interestingly, we find that the kinetic energy of solid 4He does not increase appreciably with compression for P > 85GPa. Also, we show that the Lindemann ratio in dense solid 4He amounts to 0.10 almost independently of pressure. The reliability of customary quasi-harmonic DFT (QH DFT) approaches in the description of quantum nuclear effects in solids is also studied. We find that QH DFT simulations, although provide a reasonable equation of state in agreement with experiments, are not able to reproduce correctly these critical effects in compressed 4He. In particular, we disclose huge discrepancies of at least ~50% in the calculated 4He kinetic energies using both the QH DFT and present DFT-DMC methods.
10 pages, 9 figures
References in corpus (8)
- Ab initio melting curve of molybdenum by the phase coexistence method
- Population size bias in Diffusion Monte Carlo
- Bose-Einstein Quantum Statistics and the Ground State of Solid 4He
- Quantum Monte Carlo, Density Functional Theory, and Pair Potential Studies of Solid Neon
- Comment on 'Molybdenum at High Pressure and Temperature: Melting from Another Solid Phase'
- Zero-temperature generalized phase diagram of the 4d transition metals under pressure
- Atomic kinetic energy, momentum distribution and structure of solid neon at zero-temperature
- Two-dimensional molecular para-hydrogen and ortho-deuterium at zero temperature