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.

The generalized quasiharmonic approximation via space group irreducible derivatives

Approximations in first-principles volumetric thermal expansion determination

Thermodynamics of Solids Including Anharmonicity Through Quasiparticle Theory

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

Inherent Anharmonicity of Harmonic Solids

The Phonon Quasiparticle Approach for Anharmonic Properties of Solids

Quasiharmonic elastic constants corrected for deviatoric thermal stresses

Group-theoretical Method for Physical Property Tensors of Quasicrystals

Ab initio elasticity at finite temperature and stress in ferroelectrics

Quasiharmonic Thermal Elasticity of Crystals: an Analytical Approach

Time-Dependent Self Consistent Harmonic Approximation: Anharmonic nuclear quantum dynamics and time correlation functions

Anharmonic lattice dynamics via the special displacement method

An efficient method for calculating thermoelastic properties

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

Certain problems in constrained cubic quasicrystals: General solutions and infinite-space Green’s functions

Adding anisotropy to the standard quasi-harmonic approximation still fails in several ways to capture organic crystal thermodynamics

Generalized dynamics of moving dislocations in quasicrystals

Anharmonic Vibrational States of Double-Well Potentials in the Solid State from DFT Calculations

Renormalized Phonon Microstructures At High Temperatures From First-Principles Calculations: Methodologies And Applications In Studying Strong Anharmonic Vibrations Of Solids

Beyond Gaussian fluctuations of quantum anharmonic nuclei

Primitive to conventional geometry projection for efficient phonon transport calculations