26 #ifndef O2SCL_BPS_EOS_H 27 #define O2SCL_BPS_EOS_H 29 #include <o2scl/part.h> 30 #include <o2scl/mroot.h> 31 #include <o2scl/eos_base.h> 32 #include <o2scl/mroot_hybrids.h> 33 #include <o2scl/fermion.h> 34 #include <o2scl/nucmass.h> 36 #ifndef DOXYGEN_NO_O2NS 132 virtual const char *
type() {
return "eos_crust"; }
155 #ifndef DOXYGEN_INTERNAL 161 virtual int eq274(
size_t nv,
const ubvector &nx, ubvector &ny,
165 double gibbs(
int Z,
int A);
168 double energy(
double barn,
int Z,
int A);
180 #ifndef DOXYGEN_NO_O2NS Equation of state base class.
Baym-Pethick-Sutherland equation of state.
virtual int calc_density(double barn, thermo &th, int &Z, int &A)
Calculate the equation of state as a function of the baryon number density barn.
nucmass * nmp
The nuclear mass formula.
int calc_density_fixedA(double barn, thermo &th, int &Z, int A)
Compute the ground state assuming a fixed atomic number.
virtual int calc_pressure(thermo &th, double &barn, int &Z, int &A)
Calculate the equation of state as a function of the pressure.
nucmass_semi_empirical def_mass
Default mass formula.
virtual const fermion & get_electron()
Get a pointer to the electron.
mroot_hybrids< mm_funct11 > gs
A solver to solve Eq. 2.7.4.
int set_mass_formula(nucmass &nm)
Set the nuclear mass formula to be used.
virtual double mass_formula(int Z, int A)
The mass formula.
double gibbs(int Z, int A)
The Gibbs free energy.
virtual const char * type()
Return string denoting type ("eos_crust")
virtual double lattice_energy(int Z)
The electron lattice energy.
double energy(double barn, int Z, int A)
The energy density.
virtual int eq274(size_t nv, const ubvector &nx, ubvector &ny, int &Zt)
Solve Equation 2.7.4 for a given pressure.
fermion_zerot fzt
Zero-temperature thermodynamics for the electrons.
fermion e
The electron thermodynamics.