Abstract:We here introduce a Fortran code that computes anharmonic free energy of solids from first-principles based on our phonon quasiparticle approach. In this code, phonon quasiparticle properties, i.e., renormalized phonon frequencies and lifetimes, are extracted from mode-projected velocity auto-correlation functions (VAF) of modes sampled by molecular dynamics (MD) trajectories. Using renormalized frequencies as input, the code next constructs an effective harmonic force constant matrix to calculate anharmonic phonon dispersions over the whole Brillouin zone and thus the anharmonic free energy in the thermodynamic limit (N→∞). A detailed description of all the input parameters and the subroutines is provided as well. We illustrate the use of this code to compute ab initio temperature-dependent anharmonic phonons of Si in the diamond structure.Program summaryProgram Title: phqProgram Files doi: http://dx.doi.org/10.17632/sk4jsjc6p9.1Licensing provisions: GPLv3Programming language: Fortran 90Nature of problem: Accurate free energy calculations deliver predictive thermodynamic properties of solids. Although the quasi-harmonic approximation (QHA) has been widely used for various materials, its validity at very high temperatures is not guaranteed because of intrinsic anharmonic effects. The QHA cannot be used even at low temperatures for strongly anharmonic systems, e.g., those that are stabilized by anharmonic fluctuations. When anharmonicity is non-negligible, phonon frequencies display a pronounced dependence on temperature, which impacts thermodynamic properties. Therefore, the calculation of anharmonic phonon dispersions is critical for an accurate estimation of the free energy.Solution method: We calculate phonon quasiparticle properties, i.e., renormalized phonon frequencies and lifetimes, by combining molecular dynamics (MD) trajectories and harmonic normal modes (HNM). MD velocities are projected into HNMs and the mode-projected velocity auto-correlation functions (VAF) are then calculated. Quasiparticle properties are then obtained by analyzing the VAFs, if they are well defined. Using the renormalized phonon frequencies and HNMs we build an effective (temperature-dependent) harmonic force constant matrix from which the anharmonic phonon dispersions can be obtained.

Approximations in first-principles volumetric thermal expansion determination

The generalized quasiharmonic approximation via space group irreducible derivatives

Efficient Calculation of the Lattice Thermal Conductivity by Atomistic Simulations with Ab Initio Accuracy

Efficient first-principles approach to Gibbs free energy with thermal expansion

Thermodynamics of Solids Including Anharmonicity Through Quasiparticle Theory

An Assessment of Approximate Methods for Anharmonic Free Energies

Thermal expansion in insulating solids from first principles

Phq: A Fortran Code to Compute Phonon Quasiparticle Properties and Dispersions.

Inherent Anharmonicity of Harmonic Solids

Descriptors for thermal expansion in solids

How Accurate for Phonon Models to Predict the Thermodynamics Properties of Crystals.

Revisiting first-principles thermodynamics by quasiharmonic approach: Application to study thermal expansion of additively-manufactured Inconel 625

Elastic moduli and thermal conductivity of quantum materials at finite temperature

The Phonon Quasiparticle Approach for Anharmonic Properties of Solids

Thermal conductivity prediction by atomistic simulation methods: Recent advances and detailed comparison

Primitive to conventional geometry projection for efficient phonon transport calculations

Ab initio study of the volume dependence of dynamical and thermodynamical properties of silicon

Lattice thermal expansion of complex oxides with intrinsic anharmonicity

Calculation of the Anisotropic Coefficients of Thermal Expansion: A First-Principles Approach

Understanding mechanisms of thermal expansion in PbTiO$_3$ thin-films from first principles: role of high-order phonon-strain anharmonicity