All Data Structures Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Pages
Public Member Functions
o2scl::eos_had_sym4_rmf Class Reference

A version of eos_had_rmf to separate potential and kinetic contributions. More...

Detailed Description


References:

Created for Steiner06.

Definition at line 104 of file eos_had_sym4.h.

#include <eos_had_sym4.h>

Inheritance diagram for o2scl::eos_had_sym4_rmf:
o2scl::eos_had_rmf o2scl::eos_had_sym4_base o2scl::eos_had_temp_pres_base o2scl::eos_had_temp_base o2scl::eos_had_base o2scl::eos_base

Public Member Functions

virtual int calc_e_sep (fermion &ne, fermion &pr, double &ed_kin, double &ed_pot, double &mu_n_kin, double &mu_p_kin, double &mu_n_pot, double &mu_p_pot)
 Compute the potential and kinetic parts separately.
 
- Public Member Functions inherited from o2scl::eos_had_rmf
virtual const char * type ()
 Return string denoting type ("eos_had_rmf")
 
virtual int calc_e (fermion &ne, fermion &pr, thermo &lth)
 Equation of state as a function of density. More...
 
virtual int calc_p (fermion &ne, fermion &pr, thermo &lth)
 Equation of state as a function of chemical potential. More...
 
virtual int calc_eq_p (fermion &neu, fermion &p, double sig, double ome, double rho, double &f1, double &f2, double &f3, thermo &th)
 Equation of state and meson field equations as a function of chemical potentials. More...
 
virtual int calc_eq_temp_p (fermion &ne, fermion &pr, double temper, double sig, double ome, double rho, double &f1, double &f2, double &f3, thermo &th)
 Equation of state and meson field equations as a function of chemical potentials at finite temperature. More...
 
virtual int calc_temp_p (fermion &ne, fermion &pr, double T, thermo &lth)
 Equation of state as a function of chemical potential. More...
 
int calc_temp_e (fermion &ne, fermion &pr, double T, thermo &lth)
 Equation of state as a function of densities at finite temperature.
 
int fix_saturation (double guess_cs=4.0, double guess_cw=3.0, double guess_b=0.001, double guess_c=-0.001)
 Calculate cs, cw, cr, b, and c from the saturation properties. More...
 
virtual void saturation ()
 Calculate properties of nuclear matter at the saturation density. More...
 
double fesym_fields (double sig, double ome, double nb)
 Calculate symmetry energy assuming the field equations have already been solved. More...
 
double fcomp_fields (double sig, double ome, double nb)
 Calculate the compressibility assuming the field equations have already been solved. More...
 
void fkprime_fields (double sig, double ome, double nb, double &k, double &kprime)
 Calculate compressibilty and kprime assuming the field equations have already been solved. More...
 
int field_eqs (size_t nv, const ubvector &x, ubvector &y)
 A function for solving the field equations. More...
 
int field_eqsT (size_t nv, const ubvector &x, ubvector &y)
 A function for solving the field equations at finite temperature. More...
 
virtual int set_fields (double sig, double ome, double lrho)
 Set a guess for the fields for the next call to calc_e(), calc_p(), or saturation()
 
int get_fields (double &sig, double &ome, double &lrho)
 Return the most recent values of the meson fields. More...
 
void check_naturalness (eos_had_rmf &re)
 Set the coefficients of a eos_had_rmf object to their limits from naturalness. More...
 
void naturalness_limits (double value, eos_had_rmf &re)
 Provide the maximum values of the couplings assuming a limit on naturalness. More...
 
virtual int set_sat_mroot (mroot< mm_funct11, ubvector, jac_funct11 > &mrx)
 Set class mroot object for use calculating saturation density.
 
- Public Member Functions inherited from o2scl::eos_had_temp_base
virtual int calc_liqgas_dens_temp_e (fermion &n1, fermion &p1, fermion &n2, fermion &p2, double T, thermo &th1, thermo &th2)
 Compute liquid-gas phase transition densities using eos_had_temp_base::calc_temp_e() . More...
 
virtual int calc_liqgas_temp_e (fermion &n1, fermion &p1, fermion &n2, fermion &p2, double nB, double Ye, double T, thermo &th1, thermo &th2, double &chi)
 Compute the liquid-gas phase transition using eos_had_temp_base::calc_temp_e() . More...
 
virtual int calc_liqgas_beta_temp_e (fermion &n1, fermion &p1, fermion &n2, fermion &p2, double nB, double T, thermo &th1, thermo &th2, double &Ye, double &chi)
 Compute the liquid-gas phase transition in beta-equilibrium using eos_had_temp_base::calc_temp_e() . More...
 
virtual double fesym_T (double nb, double T, double delta=0.0)
 Compute the symmetry energy at finite temperature.
 
virtual double fsyment_T (double nb, double T, double delta=0.0)
 Compute the symmetry entropy at finite temperature.
 
virtual void set_fermion_eval_thermo (fermion_eval_thermo &f)
 Computing finite-temperature integrals. More...
 
- Public Member Functions inherited from o2scl::eos_had_base
virtual double fcomp (double nb, double delta=0.0)
 Calculate the incompressibility in $ \mathrm{fm}^{-1} $ using calc_e() More...
 
virtual double fcomp_err (double nb, double delta, double &unc)
 Compute the incompressibility and its uncertainty. More...
 
virtual double feoa (double nb, double delta=0.0)
 Calculate the energy per baryon in $ \mathrm{fm}^{-1} $ using calc_e() More...
 
virtual double fesym (double nb, double delta=0.0)
 Calculate symmetry energy of matter in $ \mathrm{fm}^{-1} $ using calc_dmu_delta() . More...
 
virtual double fesym_err (double nb, double delta, double &unc)
 Calculate symmetry energy of matter and its uncertainty. More...
 
virtual double fesym_slope (double nb, double delta=0.0)
 The symmetry energy slope parameter. More...
 
virtual double fesym_curve (double nb, double delta=0.0)
 The curvature of the symmetry energy.
 
virtual double fesym_skew (double nb, double delta=0.0)
 The skewness of the symmetry energy.
 
virtual double fesym_diff (double nb)
 Calculate symmetry energy of matter as energy of neutron matter minus the energy of nuclear matter. More...
 
virtual double feta (double nb)
 The strength parameter for quartic terms in the symmetry energy.
 
virtual double feta_prime (double nb)
 The derivative of the strength parameter for quartic terms in the symmetry energy.
 
virtual double fkprime (double nb, double delta=0.0)
 Calculate skewness of nuclear matter using calc_e() More...
 
virtual double fmsom (double nb, double delta=0.0)
 Calculate reduced neutron effective mass using calc_e() More...
 
virtual double f_effm_neut (double nb, double delta=0.0)
 Neutron effective mass.
 
virtual double f_effm_prot (double nb, double delta=0.0)
 Proton effective mass.
 
virtual double f_effm_scalar (double nb, double delta=0.0)
 Scalar effective mass. More...
 
virtual double f_effm_vector (double nb, double delta=0.0)
 Vector effective mass. More...
 
virtual double fn0 (double delta, double &leoa)
 Calculate saturation density using calc_e() More...
 
double calc_dmu_delta (double delta, double nb)
 Compute the difference between neutron and proton chemical potentials as a function of the isospin asymmetry. More...
 
double calc_musum_delta (double delta, double nb)
 Compute the sum of the neutron and proton chemical potentials as a function of the isospin asymmetry. More...
 
double calc_pressure_nb (double nb, double delta=0.0)
 Compute the pressure as a function of baryon density at fixed isospin asymmetry. More...
 
double calc_edensity_nb (double nb, double delta=0.0)
 Compute the energy density as a function of baryon density at fixed isospin asymmetry. More...
 
void const_pf_derivs (double nb, double pf, double &dednb_pf, double &dPdnb_pf)
 Compute derivatives at constant proton fraction.
 
double calc_press_over_den2 (double nb, double delta=0.0)
 Calculate pressure / baryon density squared in nuclear matter as a function of baryon density at fixed isospin asymmetry. More...
 
double calc_edensity_delta (double delta, double nb)
 Calculate energy density as a function of the isospin asymmetry at fixed baryon density. More...
 
int nuc_matter_p (size_t nv, const ubvector &x, ubvector &y, double nn0, double np0)
 Solve for the chemical potentials given the densities. More...
 
int nuc_matter_e (size_t nv, const ubvector &x, ubvector &y, double mun0, double mup0)
 Solve for the densities given the chemical potentials. More...
 
virtual void set_mroot (mroot<> &mr)
 Set class mroot object for use in calculating chemical potentials from densities. More...
 
virtual void set_sat_root (root<> &mr)
 Set class mroot object for use calculating saturation density. More...
 
virtual void set_sat_deriv (deriv_base<> &de)
 Set deriv_base object to use to find saturation properties.
 
virtual void set_sat_deriv2 (deriv_base<> &de)
 Set the second deriv_base object to use to find saturation properties. More...
 
virtual void set_n_and_p (fermion &n, fermion &p)
 Set neutron and proton.
 
void gradient_qij (fermion &n, fermion &p, thermo &th, double &qnn, double &qnp, double &qpp, double &dqnndnn, double &dqnndnp, double &dqnpdnn, double &dqnpdnp, double &dqppdnn, double &dqppdnp)
 Calculate coefficients for gradient part of Hamiltonian. More...
 
- Public Member Functions inherited from o2scl::eos_base
virtual void set_thermo (thermo &th)
 Set class thermo object.
 
virtual const thermoget_thermo ()
 Get class thermo object.
 
- Public Member Functions inherited from o2scl::eos_had_sym4_base
virtual int calc_e_alpha (fermion &ne, fermion &pr, thermo &lth, double &alphak, double &alphap, double &alphat, double &diff_kin, double &diff_pot, double &ed_kin_nuc, double &ed_pot_nuc)
 Compute alpha at the specified density.
 
virtual double calc_muhat (fermion &ne, fermion &pr)
 Compute $ \hat{\mu} $, the out-of-whack parameter.
 

Additional Inherited Members

- Public Types inherited from o2scl::eos_had_base
typedef
boost::numeric::ublas::vector
< double > 
ubvector
 
- Data Fields inherited from o2scl::eos_had_rmf
size_t calc_e_steps
 The number of separate calls to the solver that the calc_e functions take (default 20) More...
 
bool calc_e_relative
 If true, solve for relative densities rather than absolute densities (default false) More...
 
bool zm_mode
 Modifies method of calculating effective masses (default false)
 
int verbose
 Verbosity parameter. More...
 
bool err_nonconv
 If true, throw exceptions when the function calc_e() does not converge (default true)
 
double mnuc
 The scale $ M $. More...
 
double ms
 $ \sigma $ mass (in $ \mathrm{fm}^{-1} $ )
 
double mw
 $ \omega $ mass (in $ \mathrm{fm}^{-1} $ )
 
double mr
 $ \rho $ mass (in $ \mathrm{fm}^{-1} $ )
 
double cs
 
double cw
 
double cr
 
double b
 
double c
 
double zeta
 
double xi
 
double a1
 
double a2
 
double a3
 
double a4
 
double a5
 
double a6
 
double b1
 
double b2
 
double b3
 
mroot_hybrids def_sat_mroot
 The default solver for calculating the saturation density. More...
 
- Data Fields inherited from o2scl::eos_had_temp_base
fermion_eff def_fet
 Default fermion thermodynamics object.
 
- Data Fields inherited from o2scl::eos_had_base
double eoa
 Binding energy.
 
double comp
 Compression modulus.
 
double esym
 Symmetry energy.
 
double n0
 Saturation density.
 
double msom
 Effective mass (neutron)
 
double kprime
 Skewness.
 
fermion def_neutron
 The defaut neutron. More...
 
fermion def_proton
 The defaut proton. More...
 
deriv_gsl def_deriv
 The default object for derivatives. More...
 
deriv_gsl def_deriv2
 The second default object for derivatives. More...
 
mroot_hybrids def_mroot
 The default solver. More...
 
root_cern def_sat_root
 The default solver for calculating the saturation density. More...
 
- Data Fields inherited from o2scl::eos_base
thermo def_thermo
 The default thermo object.
 
- Protected Member Functions inherited from o2scl::eos_had_rmf
int fix_saturation_fun (size_t nv, const ubvector &x, ubvector &y)
 The function for fix_saturation()
 
virtual int zero_pressure (size_t nv, const ubvector &ex, ubvector &ey)
 Compute matter at zero pressure (for saturation())
 
virtual int calc_e_solve_fun (size_t nv, const ubvector &ex, ubvector &ey)
 The function for calc_e()
 
virtual int calc_temp_e_solve_fun (size_t nv, const ubvector &ex, ubvector &ey)
 The function for calc_temp_e()
 
int calc_cr (double sig, double ome, double nb)
 Calculate the cr coupling given sig and ome at the density 'nb'. More...
 
- Protected Member Functions inherited from o2scl::eos_had_temp_base
int nuc_matter_temp_e (size_t nv, const ubvector &x, ubvector &y, double nn0, double np0, double T)
 Solve for nuclear matter at finite temperature given density.
 
int nuc_matter_temp_p (size_t nv, const ubvector &x, ubvector &y, double mun0, double mup0, double T)
 Solve for nuclear matter at finite temperature given mu.
 
int liqgas_dens_solve (size_t nv, const ubvector &x, ubvector &y, fermion &n1, fermion &p1, fermion &n2, fermion &p2, double T, thermo &th1, thermo &th2)
 Solve for the liquid gas phase transition as a function of the densities.
 
int liqgas_solve (size_t nv, const ubvector &x, ubvector &y, fermion &n1, fermion &p1, fermion &n2, fermion &p2, double nB0, double Ye0, double T, thermo &th1, thermo &th2)
 Solve for the liquid-gas phase transition at fixed baryon density and electron fraction.
 
int liqgas_beta_solve (size_t nv, const ubvector &x, ubvector &y, fermion &n1, fermion &p1, fermion &n2, fermion &p2, double nB0, double T, thermo &th1, thermo &th2, fermion &e)
 Solve for the liquid-gas phase transition in beta-equilibrium.
 
double calc_entropy_delta (double delta, double nb, double T)
 Compute the entropy.
 
double calc_dmu_delta_T (double delta, double nb, double T)
 Compute the difference between the neutron and proton chemical potentials.
 
- Protected Member Functions inherited from o2scl::eos_had_base
double t1_fun (double barn)
 Compute t1 for gradient_qij().
 
double t2_fun (double barn)
 Compute t2 for gradient_qij().
 
- Protected Attributes inherited from o2scl::eos_had_rmf
double n_baryon
 Temporary baryon density.
 
double n_charge
 Temporary charge density. More...
 
double fe_temp
 Temperature for solving field equations at finite temperature.
 
bool ce_neut_matter
 For calc_e(), if true, then solve for neutron matter.
 
bool ce_prot_matter
 For calc_e(), if true, then solve for proton matter.
 
bool guess_set
 True if a guess for the fields has been given.
 
mroot< mm_funct11, ubvector,
jac_funct11 > * 
sat_mroot
 The solver to compute saturation properties.
 
double ce_temp
 Temperature storage for calc_temp_e()
 
double sigma
 
double omega
 
double rho
 
- Protected Attributes inherited from o2scl::eos_had_temp_base
fermion_eval_thermofet
 Fermion thermodynamics (default is def_fet)
 
- Protected Attributes inherited from o2scl::eos_had_base
mrooteos_mroot
 The EOS solver.
 
rootsat_root
 The solver to compute saturation properties.
 
deriv_basesat_deriv
 The derivative object for saturation properties.
 
deriv_basesat_deriv2
 The second derivative object for saturation properties.
 
fermionneutron
 The neutron object.
 
fermionproton
 The proton object.
 
- Protected Attributes inherited from o2scl::eos_base
thermoeos_thermo
 A pointer to the thermo object.
 
- Protected Attributes inherited from o2scl::eos_had_sym4_base
fermion e
 An electron for the computation of the $ \hat{\mu}$.
 
fermion_zerot fzt2
 Zero-temperature fermion thermodynamics.
 

The documentation for this class was generated from the following file:

Documentation generated with Doxygen. Provided under the GNU Free Documentation License (see License Information).
Hosted at Get Object-oriented Scientific Computing
Lib at SourceForge.net. Fast, secure and Free Open Source software
downloads..