Program title: Phasego2
Catalogue identifier: AEVQ_v2_0
Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AEVQ_v2_0.html
Program obtainable from: CPC Program Library, Queen’s University, Belfast, N. Ireland
Licensing provisions: GNU General Public License, version 3
No. of lines in distributed program, including test data, etc.: 5867093
No. of bytes in distributed program, including test data, etc.: 56696690
Distribution format: tar.gz
Programming language: Python (versions 2.4 and later).
Computer: Any computer that can run Python (versions 2.4 and later).
Operating system: Any operating system that can run Python.
RAM: 50 M bytes
Classification: 7.8.
Catalogue identifier of previous version: AEVQ_v1_0
Journal reference of previous version: Comput. Phys. Comm. 191 (2015) 150
External routines: Numpy [1], Scipy [2], Matplotlib [3]
Does the new version supersede the previous version?: Yes
Nature of problem: The anharmonic lattice vibrations have more and more contributions to Gibbs free energy when temperature is beyond Debye temperature and goes up to melting temperature. The quasi-harmonic approximation (QHA) only includes part of the anharmonic effects due to the volume variation of frequency. At high temperature, the neglect of phonon–phonon interactions in the QHA will introduce errors and even arrive at wrong conclusions. However, the calculation of the anharmonic free energy from the phonon–phonon interactions is very complicated to perform.
Solution method: The contribution of phonon–phonon interactions to free energy can be naturally extracted from the high-temperature phonon density of states (DOS). The extraction and inclusion of anharmonic free energy are simplified and automated in the updated version. The thermal properties of materials are then corrected by including full anharmonic effects with respective to the QHA results.
Reasons for new version: We have improved the package considerably to include full anharmonic contributions of phonon–phonon interactions. The calculations of thermal properties of materials are automated and easy to implement after anharmonic corrections.
Summary of revisions:
The anharmonic free energy extraction functionality is added. Now the Gibbs free energy can include full anharmonic effects.
The thermal expansion coefficients are deduced by taking into account full anharmonic contributions.
The bulk moduli, the heat capacities, the thermal pressures can include and reflect the contributions of phonon–phonon interactions.
The Hugoniot pressure–volume–temperature relations, the Grüneisen parameters, and the Debye temperatures are corrected by anharmonic effects.
The Polynomial curve fitting of Helmholtz free energy is added.
The unit of Gibbs free energy in the output files is changed.
A number of bugs have been corrected.
Unusual features: The Gibbs free energy of phonon–phonon interactions are automatically extracted from the high-temperature phonon DOS. The anharmonic effects of all the thermal properties are automatically corrected.
Additional comments: The new version of this package can treat the high-temperature phonon density of states data from many methods, including the molecular dynamics (MD) simulations [4], the self-consistent ab initio lattice dynamics (SCAILD) calculations [5], or other simulation methods.
The distribution file for this program is over 56 Mbytes and therefore is not delivered directly when download or Email is requested. Instead a html file giving details of how the program can be obtained is sent.
Running time: The examples provided in the distribution take less than 5 minute to run.
References:
[1] www.numpy.org.
[2] www.scipy.org.
[3] www.matplotlib.org.
[4] O. Hellman, I. A. Abrikosov, and S. I. Simak, Phys. Rev. B, 84 (2011) 180301.
[5] P. Souvatzis, O. Eriksson, M. I. Katsnelson, S. P. Rudin, Phys. Rev. Lett., 100 (2008) 095901.