phq: A Fortran code to compute phonon quasiparticle properties and dispersions

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 summary Program Title: phq Program Files doi: http://dx.doi.org/10.17632/sk4jsjc6p9.1 Licensing provisions: GPLv3 Programming language: Fortran 90 Nature 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.

[1]  Hoover,et al.  Canonical dynamics: Equilibrium phase-space distributions. , 1985, Physical review. A, General physics.

[2]  Andrew G. Glen,et al.  APPL , 2001 .

[3]  M. Born,et al.  Dynamical Theory of Crystal Lattices , 1954 .

[4]  R. Wentzcovitch,et al.  Spin transition in (Mg,Fe)SiO3 perovskite under pressure , 2008 .

[5]  M. Balkanski,et al.  Anharmonic effects in light scattering due to optical phonons in silicon , 1983 .

[6]  D. Wallace,et al.  Thermodynamics of Crystals , 1972 .

[7]  Kresse,et al.  Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set. , 1996, Physical review. B, Condensed matter.

[8]  A. Hofmeister,et al.  Mantle values of thermal conductivity and the geotherm from phonon lifetimes , 1999, Science.

[9]  Burke,et al.  Generalized Gradient Approximation Made Simple. , 1996, Physical review letters.

[10]  Bartomeu Monserrat,et al.  Lattice dynamics and electron-phonon coupling calculations using nondiagonal supercells , 2015, 1510.04418.

[11]  L. Paulatto,et al.  Phonon Collapse and Second-Order Phase Transition in Thermoelectric SnSe. , 2018, Physical review letters.

[12]  A. McGaughey,et al.  Predicting phonon properties and thermal conductivity from anharmonic lattice dynamics calculations and molecular dynamics simulations , 2009 .

[13]  S. Jain,et al.  Earth as a Planet , 2014 .

[14]  R. Wentzcovitch,et al.  Dynamic stabilization of cubic Ca Si O 3 perovskite at high temperatures and pressures from ab initio molecular dynamics , 2014 .

[15]  Dario Alfè,et al.  PHON: A program to calculate phonons using the small displacement method , 2009, Comput. Phys. Commun..

[16]  William H. Press,et al.  Numerical Recipes in Fortran 77 , 1992 .

[17]  E. Wang,et al.  Vibrational recognition of hydrogen-bonded water networks on a metal surface. , 2002, Physical review letters.

[18]  Abel Carreras,et al.  DynaPhoPy: A code for extracting phonon quasiparticles from molecular dynamics simulations , 2017, Comput. Phys. Commun..

[19]  R. Wentzcovitch,et al.  Premelting hcp to bcc Transition in Beryllium. , 2017, Physical review letters.

[20]  Hafner,et al.  Ab initio molecular-dynamics simulation of the liquid-metal-amorphous-semiconductor transition in germanium. , 1994, Physical review. B, Condensed matter.

[21]  Tianli Feng,et al.  Phonon anharmonic frequency shift induced by four-phonon scattering calculated from first principles , 2018, Journal of Applied Physics.

[22]  Stefano de Gironcoli,et al.  QUANTUM ESPRESSO: a modular and open-source software project for quantum simulations of materials , 2009, Journal of physics. Condensed matter : an Institute of Physics journal.

[23]  O. Delaire,et al.  Phonon self-energy and origin of anomalous neutron scattering spectra in SnTe and PbTe thermoelectrics. , 2013, Physical review letters.

[25]  William H. Press,et al.  Numerical Recipes: FORTRAN , 1988 .

[26]  Wang,et al.  Tight-binding molecular-dynamics study of phonon anharmonic effects in silicon and diamond. , 1990, Physical review. B, Condensed matter.

[27]  Raphael Tsu,et al.  Temperature dependence of silicon Raman lines , 1982 .

[28]  Blöchl,et al.  Projector augmented-wave method. , 1994, Physical review. B, Condensed matter.

[29]  Stefano de Gironcoli,et al.  Phonons and related crystal properties from density-functional perturbation theory , 2000, cond-mat/0012092.

[30]  Stefan Goedecker,et al.  ABINIT: First-principles approach to material and nanosystem properties , 2009, Comput. Phys. Commun..

[31]  Dong-Bo Zhang,et al.  Phonon quasiparticles and anharmonic free energy in complex systems. , 2013, Physical review letters.

[32]  Ladd,et al.  Lattice thermal conductivity: A comparison of molecular dynamics and anharmonic lattice dynamics. , 1986, Physical review. B, Condensed matter.