Gibbs2: A new version of the quasiharmonic model code. II. Models for solid-state thermodynamics, features and implementation

Abstract

In the second article of the series, we present the Gibbs2 code, a Fortran90 reimplementation of the original Gibbs program [Comput. Phys. Commun. 158 (2004) 57] for the calculation of pressure–temperature dependent thermodynamic properties of solids under the quasiharmonic approximation. We have taken advantage of the detailed analysis carried out in the first paper to implement robust fitting techniques. In addition, new models to introduce temperature effects have been incorporated, from the simple Debye model contained in the original article to a full quasiharmonic model that requires the phonon density of states at each calculated volume. Other interesting novel features include the empirical energy corrections, that rectify systematic errors in the calculation of equilibrium volumes caused by the choice of the exchange-correlation functional, the electronic contributions to the free energy and the automatic computation of phase diagrams. Full documentation in the form of a userʼs guide and a complete set of tests and sample data are provided along with the source code. Program summaryProgram title:Gibbs2Catalogue identifier: AEJI_v1_0Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AEJI_v1_0.htmlProgram obtainable from: CPC Program Library, Queenʼs University, Belfast, N. IrelandLicensing provisions: GNU General Public License, v3No. of lines in distributed program, including test data, etc.: 936 087No. of bytes in distributed program, including test data, etc.: 8 596 671Distribution format: tar.gzProgramming language: Fortran90Computer: Any running Unix/LinuxOperating system: Unix, GNU/LinuxClassification: 7.8External routines: Part of the minpack, pppack and slatec libraries (downloaded from www.netlib.org) are distributed along with the program.Nature of problem: Given the static E(V) curve, and possibly vibrational information such as the phonon density of states, calculate the equilibrium volume and thermodynamic properties of a solid at arbitrary temperatures and pressures in the framework of the quasiharmonic approximation.Additional comments: A detailed analysis concerning the fitting of equations of state has been carried out in the first part of this article, and implemented in the code presented here.Running time: The tests provided only take a few seconds to run.