Class eos_crust_virial (o2scl)¶
-
class
o2scl
::
eos_crust_virial
: public o2scl::eos_had_temp_pres_base¶ Virial EOS for neutrons, protons, deuterons, and alpha particles.
Virial EOS from Horowitz06 and Horowitz06b .
- Note
This class likely fails for temperatures smaller than 1 MeV and larger than 20 MeV.
- Warning
This class is implemented as a o2scl::eos_had_temp_pres_base object because it might be helpful to be able to use o2scl::eos_had_temp_pres_base::calc_temp_e(), but because of the alpha particles and deuterons, some of the other o2scl::eos_had_base methods don’t have the correct interpretation.
Interpolation objects for virial coefficients
-
std::vector<double>
Tv
¶
-
std::vector<double>
bnv
¶
-
std::vector<double>
Tbnpv
¶
-
std::vector<double>
bpnv
¶
-
std::vector<double>
Tbpnpv
¶
-
std::vector<double>
banv
¶
-
std::vector<double>
Tbanpv
¶
-
std::vector<double>
bav
¶
-
std::vector<double>
Tbapv
¶
-
interp_vec<std::vector<double>>
ibn
¶
-
interp_vec<std::vector<double>>
iTbnp
¶
-
interp_vec<std::vector<double>>
ibpn
¶
-
interp_vec<std::vector<double>>
iTbpnp
¶
-
interp_vec<std::vector<double>>
iban
¶
-
interp_vec<std::vector<double>>
iTbanp
¶
-
interp_vec<std::vector<double>>
iba
¶
-
interp_vec<std::vector<double>>
iTbap
¶
-
boson
alpha
¶ Internal alpha particle.
-
boson
deuteron
¶ Internal deuteron.
-
eos_crust_virial
()¶
-
~eos_crust_virial
()¶
-
int
calc_p
(fermion &ne, fermion &pr, thermo &th)¶ Equation of state as a function of the chemical potentials.
The virial coefficients and their temperature derivatives
These functions assume that the temperature is specified in MeV. Note that the temperature derivative functions are multiplied by the temperature after taking the derivative to make them unitless (as in Horowitz06), thus all of these functions return values without units.
-
double
bn
(double T)¶
-
double
ban
(double T)¶
-
double
ba
(double T)¶
-
double
bpn
(double T)¶
-
double
Tbn_prime
(double T)¶
-
double
Tban_prime
(double T)¶
-
double
Tba_prime
(double T)¶
-
double
Tbpn_prime
(double T)¶
-
int
calc_temp_p
(fermion &n, fermion &p, double T, thermo &th)¶ Equation of state as a function of the chemical potentials at finite temperature.
As in o2scl::eos_had_temp_base::calc_temp_p() the temperature argument should be in \( \mathrm{fm}^{-1} \), even though the virial coefficient functions take their temperature arguments in MeV.
-
int
calc_temp_p_alpha
(fermion &n, fermion &p, boson &d, boson &a, double T, thermo &th)¶ Equation of state as a function of the chemical potentials at finite temperature with alpha particles and deuterons.
-
void
fit
()¶ Fit the virial coefficients to analytical functions.