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

被引:814
|
作者
Otero-de-la-Roza, A. [1 ]
Abbasi-Perez, David [1 ]
Luana, Victor [1 ]
机构
[1] Univ Oviedo, Fac Quim, Dept Quim Fis & Analit, E-33006 Oviedo, Spain
关键词
Equations of state in solids; Equilibrium properties of solids; Data analysis; Treatment of noisy data; Discontinuous data; Thermal effects; MGO; COMPRESSION; PERICLASE; MINERALS; GPA;
D O I
10.1016/j.cpc.2011.05.009
中图分类号
TP39 [计算机的应用];
学科分类号
081203 ; 0835 ;
摘要
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 summary Program title: GIBBS2 Catalogue identifier: AEJI_v1_0 Program summary URL: http://cpc.cs.qub.ac.uk/summaries/AEJI_v1_0.html Program obtainable from: CPC Program Library, Queen's University, Belfast, N. Ireland Licensing provisions: GNU General Public License, v3 No. of lines in distributed program, including test data, etc.: 936 087 No. of bytes in distributed program, including test data, etc.: 8 596 671 Distribution format: tar.gz Programming language: Fortran90 Computer: Any running Unix/Linux Operating system: Unix, GNU/Linux Classification: 7.8 External 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. (C) 2011 Elsevier B.V. All rights reserved.
引用
收藏
页码:2232 / 2248
页数:17
相关论文
共 2 条