Class eos_had_rmf_hyp (o2scl)¶
-
class eos_had_rmf_hyp : public o2scl::eos_had_rmf¶
Relativistic mean field theory EOS with hyperons.
Based on [Glendenning91ro], but generalized for higher-order couplings as in ref eos_had_rmf .
Todo
In class eos_had_rmf_hyp:
The couplings in the test code match the table but the maximum masses appear much smaller than GM91. I need to check that muons are added correctly, and it might be good to compare with a different reference. This also might be due to a different crust EOS.
The interpretation of the calc_e() function is a bit unclear, so I need to more clearly figure out what that function ought to do. I don’t think it’s really used at the moment.
Hyperon objects
-
fermion def_lambda¶
The default Lambda hyperon.
-
fermion def_sigma_p¶
The default Sigma plus hyperon.
-
fermion def_sigma_z¶
The default Sigma zero hyperon.
-
fermion def_sigma_m¶
The default Sigma minus hyperon.
-
fermion def_cascade_z¶
The default Xi zero hyperon.
-
fermion def_cascade_m¶
The default Xi minus hyperon.
Hyperon-meson couplings
-
double xs¶
-
double xw¶
-
double xr¶
-
bool inc_cascade¶
If true, include cascade hyperons (default true)
-
virtual int calc_eq_hyp_p(fermion &ne, fermion &pr, fermion &lam, fermion &sigp, fermion &sigz, fermion &sigm, fermion &casz, fermion &casm, double sig, double ome, double lrho, double &f1, double &f2, double &f3, thermo <h)¶
Equation of state and meson field equations as a function of chemical potentials.
-
virtual int calc_eq_temp_hyp_p(fermion &ne, fermion &pr, fermion &lam, fermion &sigp, fermion &sigz, fermion &sigm, fermion &casz, fermion &casm, double T, double sig, double ome, double lrho, double &f1, double &f2, double &f3, thermo <h)¶
Equation of state and meson field equations as a function of chemical potentials at finite temperature.
-
void calc_xs(double lam_be)¶
Compute xs assuming a fixed value of the \( \Lambda \) binding energy in nuclear matter in \( \mathrm{fm}^{-1} \).
-
void calc_xw(double lam_be)¶
Compute xs assuming a fixed value of the \( \Lambda \) binding energy in nuclear matter in \( \mathrm{fm}^{-1} \).
-
virtual int calc_hyp_e(double n_baryon, double n_charge, fermion &ne, fermion &pr, fermion &lam, fermion &sigp, fermion &sigz, fermion &sigm, fermion &casz, fermion &casm, thermo <h)¶
Equation of state as a function of density at fixed baryon and charge density presuming the hyperons are in beta-equilibrium with the nucleons.
Initial guesses for the chemical potentials are taken from the user-given values. Initial guesses for the fields can be set by set_fields(), or default values will be used. After the call to calc_e(), the final values of the fields can be accessed through get_fields().
-
virtual int calc_temp_hyp_e(double n_baryon, double n_charge, fermion &ne, fermion &pr, fermion &lam, fermion &sigp, fermion &sigz, fermion &sigm, fermion &casz, fermion &casm, double T, thermo <h)¶
Equation of state as a function of density at fixed baryon and charge density presuming the hyperons are in beta-equilibrium with the nucleons.
Initial guesses for the chemical potentials are taken from the user-given values. Initial guesses for the fields can be set by set_fields(), or default values will be used. After the call to calc_e(), the final values of the fields can be accessed through get_fields().
-
virtual int calc_hyp_e_nobeta(double nB, double Ye, double Ys, fermion &ne, fermion &pr, fermion &lam, fermion &sigp, fermion &sigz, fermion &sigm, fermion &casz, fermion &casm, thermo <h)¶
Equation of state as a function of density out of beta equilibrium.
-
virtual int calc_hyp_e_nobeta_np(double Ys, fermion &ne, fermion &pr, fermion &lam, fermion &sigp, fermion &sigz, fermion &sigm, fermion &casz, fermion &casm, thermo <h)¶
Equation of state as a function of density out of beta equilibrium at fixed neutron and proton densities.
-
virtual void set_hyp(fermion &lam, fermion &sigp, fermion &sigz, fermion &sigm, fermion &casz, fermion &casm)¶
Set the hyperon objects.
Compute the EOS in beta-equilibrium at zero temperature.
Public Functions
-
eos_had_rmf_hyp()¶
Protected Functions
-
virtual int calc_e_solve_fun(size_t nv, const ubvector &ex, ubvector &ey)¶
The function for calc_e()
-
int calc_e_nobeta_fun(size_t nv, const ubvector &ex, ubvector &ey, double nB, double Ye, double Ys)¶
The function for calc_e_hyp_nobeta()
-
int calc_e_nobeta_np_fun(size_t nv, const ubvector &ex, ubvector &ey, double nn, double np, double Ys)¶
The function for calc_e_hyp_nobeta_np()
-
virtual int solve_beta_eq_T0(size_t nv, const ubvector &x, ubvector &y, const double &nB, eos_leptons &elep)¶
Equation for solving for beta-equilibrium at T=0.