Fortran codes
Last upgrade:
June 15, 2022.
All updates since the first distribution are listed as comments
within the program files.
|
The Fortran codes implement the formulae for the EOS of nonideal Coulomb plasmas, described by Potekhin and Chabrier (2010,2013) (see References). The subroutines are collected in two files (right-click on their names and Save Link to download):
- eos22.f contains the complete set of subroutines for nonmagnetized plasmas;
- eosmag22.f
is a supplement to eos22.f, which is needed when either a quantizing magnetic field or positrons are present. In those cases the two files must be linked together.
The first file (eos22.f) contains the following subroutines:
- MELANGE9 - for arbitrary ionic mixture, renders total (ion+electron)
pressure, internal energy, entropy, heat capacity (all
normalized to the ionic ideal-gas values), logarithmic
derivatives of pressure over temperature and density.
- EOSFI22 - nonideal (ion-ion + ion-electron + electron-electron)
contributions to the free and internal energies, pressure,
entropy, heat capacity, derivatives of pressure over
logarithm of temperature and over logarithm of density (all
normalized to the ionic ideal-gas values) for one ionic
component in a mixture.
- FITION9 - ion-ion interaction contributions to the free and internal
energies, pressure, entropy, heat capacity, derivatives of
pressure over logarithms of temperature and density.
- LIQBC - quantum additions to classical OCP liquid.
- FSCRliq8 - ion-electron (screening) contributions to the free and
internal energies, pressure, entropy, heat capacity,
derivatives of pressure over logarithms of temperature and
density in the liquid phase for one ionic component in a
mixture.
- FSCRsol8 - ion-electron (screening) contributions to the free and
internal energies, pressure, entropy, heat capacity,
derivatives of pressure over logarithms of temperature and
density for monoionic solid.
- FHARM12 - harmonic (including static-lattice and zero-point)
contributions to the free and internal energies, pressure,
entropy, heat capacity, derivatives of pressure over
logarithms of temperature and density for solid OCP.
- HLfit12 - the same as FHARM12, but only for thermal contributions.
- ANHBC - anharmonic corrections to the previous two.
- CORMIX - correction to the linear mixing rule for the Coulomb
contributions to the thermodynamic functions.
- ELECT11 - for an ideal electron gas of arbitrary
degeneracy and relativity at given temperature and
electron chemical potential, renders number density
(in atomic units), free energy, pressure, internal
energy, entropy, heat capacity (normalized to the
electron ideal-gas values), logarithmic derivatives
of pressure over temperature and density. This subroutine is actually split into ELECT11a, which realizes general fitting expressions via the Fermi-Dirac integrals, and ELECT11b, which employs the Sommerfeld expansions in the case of a very strong electron degeneracy.
- SOMMERF provides the Sommerfeld expansion up to order T2 for the relativistic Fermi-Dirac integrals of orders 1/2, 3/2, 5/2, together with the corresponding approximations for their first, second, and some third derivatives.
- SUBFERMJ - auxiliary to SOMMERF.
- EXCOR7 - electron-electron (exchange-correlation) contributions to
the free and internal energies, pressure, entropy, heat
capacity, derivatives of pressure over logarithm of
temperature and over logarithm of density (all normalized
to the classical electron ideal-gas values).
- FERINV7 - inverse non-relativistic Fermi integrals of orders -1/2,
1/2, 3/2, 5/2, and their first and second derivatives.
- BLIN9 - relativistic Fermi-Dirac integrals of orders 1/2, 3/2, 5/2, together with their first, second, and some third derivatives; uses BLIN9a, BLIN9b, and BLIN9c - three segments of the fit by Blinnikov et al. (1996).
- CHEMFIT7 - electron chemical potential at given density and
temperature, and its first derivatives over density and
temperature and the second derivative over temperature.
- CHEMFIT - redirection to CHEMFIT7 using another list of arguments (auxiliary).
- FERMI10 - Fermi distribution function (auxiliary).
The second file (eosmag22.f) contains the following subroutines:
- EOSM20 - EOS of a homogeneous plasma of electrons and one ion species.
- Input: ion charge and mass numbers (Z, A), mass density ρ in CGS, temperature T and magnetic field strength B in atomic units.
-
Output:
- electron number density ne [a.u.];
- positron number density [a.u.];
- ion Coulomb coupling constant Γ;
- electron chemical potential divided by kBT;
- ion quantum parameter: ion plasma temperature Tp divided by thermodynamic temperature T;
- indicator of the phase state (liquid or solid);
- pressure P divided by nikBT, where ni=ne/Z is the ion number density;
- internal energy per ion per kBT;
- entropy per 1 ion (assuming spin=1/2) in units of kB;
- specific heat CV;
- logarithmic derivatives
of P over ρ and T.
- CHEM20 - EOS of magnetized electron-positron ideal gas (input are magnetic field strength, temperature, and the excess of electron over positron number density (assuming overall charge neutrality) ne* = ne - npos = Znion.
- CHEMAG8 - EOS of magnetized electron ideal gas (without positrons; input are density, temperature, and magnetic field strength).
- ELECTMAG - EOS of magnetized electron ideal gas (input are density, temperature, and the chemical potential of the ideal electron gas).
- MAGNION - EOS of nonrelativistic nondegenerate ion gas in magnetic field.
- SPINION - contributions to the EOS of nonrelativistic
nondegenerate ion gas due to the spin degeneracy
and magnetic moments.
- EOSFIM14 - non-ideal parts of thermodynamic functions.
- EXCORM - exchange-correlation contribution for the electron gas in magnetic field.
- HLMAG - thermal phonon contributions to the internal energy, heat capacity, entropy, pressure, logarithmic pressure derivatives, and several auxiliary quantities used in FHARMAG.
- FHARMAG - harmonic (including static-lattice and
zero-point) contributions to the free and internal
energies, pressure, entropy, heat capacity,
derivatives of pressure over logarithms of
temperature and density for solid OCP, and the energy difference
for two characteristic orientations of the Coulomb crystal in a quantizing magnetic field.
- QDMAG - quantum-diffraction correction to the EOS of the magnetized liquid EIP.
- function darcth = inverse hyperbolic tangent (inverse dtanh)
Detailed descriptions of all subroutines are given by comment lines in the corresponding codes.
The second file (magnetic EOS), unlike the first (nonmagnetic) one, is designed for an EIP containing only one kind of the ions, but not a mixture.
[TOP]
[MAIN]
Page created on November 28, 2008,
updated on June 15, 2022
by Alexander Potekhin.
Last revision of the programs: June 15, 2022