Abstract:The quasiharmonic approximation (QHA) in combination with density-functional theory is the main computational method used to calculate thermodynamic properties under arbitrary temperature and pressure conditions. QHA can predict thermodynamic phase diagrams, elastic properties, and thermal conductivities, all of which are important in various fields of knowledge. The main drawback of QHA is that it makes spurious predictions for the volume and other properties in the high temperature limit due to its approximate treatment of anharmonicity. In this work, we propose a simple extension to QHA that fixes this problem. Our approach is based on four ingredients: i) the calculation of the n-th order force constants using randomly displaced configurations and regularized regression, ii) the calculation of temperature-dependent effective harmonic frequencies with self-consistent harmonic approximation (SCHA), iii) Allen's quasiparticle (QP) theory, which allows the calculation of the anharmonic entropy from the effective frequencies, and iv) a simple Debye-like numerical model that enables the calculation of all other thermodynamic properties from the QP entropies. The proposed method is conceptually simple, with a computational complexity similar to QHA, and allows incorporating anharmonic effects to any order. The predictions of the new method coincide with QHA in the low-temperature limit and eliminate the QHA blowout at high temperature, recovering the experimentally observed behavior of all thermodynamic properties tested. The performance of our new method is demonstrated by calculating the thermodynamic properties of geologically relevant minerals MgO and CaO. We expect this new method to be an important tool in geochemistry and materials discovery.

Path-integral simulation of solids

Path-integral approaches to strongly-coupled quantum many-body systems

Quantum Monte Carlo applied to solids

Capturing the Nuclear Quantum Effects in Molecular Dynamics for Lattice Thermal Conductivity Calculations: Using Ice As Example

Path-integral methodology and simulations of quantum thermal transport: Full counting statistics approach

Thermodynamics of Solids Including Anharmonicity Through Quasiparticle Theory

Ab Initio Path Integral Monte Carlo Simulations of the Uniform Electron Gas on Large Length Scales

Path Integral Monte Carlo Simulations for Fermion Systems: Pairing in the Electron-Hole Plasma

Introduction to Configuration Path Integral Monte Carlo

Efficient Quantum Simulation Algorithms in the Path Integral Formulation

Quantum effects on the phase diagram of nuclear-like systems

Phaseless auxiliary field quantum Monte Carlo with projector-augmented wave method for solids

Path Integral Monte Carlo Simulations of liquid $^3$He without Fixed Nodes: Structural Properties and Collective Excitations

The solid phase of He-4: A Monte Carlo simulation study

Simulation of surface processes

Solvable Hydrodynamics of Quantum Integrable Systems.

Bridging lattice-scale physics and continuum field theory with quantum Monte Carlo simulations

Development of a new quantum trajectory molecular dynamics framework

Path integral spin dynamics for quantum paramagnets

Accurate and Fast Estimation of the Continuum Limit in Path Integral Simulations of Quantum Oscillators and Crystals