Public Member Functions | List of all members
o2scl::eos_had_sym4_skyrme Class Reference

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

#include <eos_had_sym4.h>

Inheritance diagram for o2scl::eos_had_sym4_skyrme:
o2scl::eos_had_skyrme o2scl::eos_had_sym4_base o2scl::eos_had_temp_eden_base o2scl::eos_had_temp_base o2scl::eos_had_base o2scl::eos_base

Detailed Description


References:

Created for Steiner06.

Definition at line 144 of file eos_had_sym4.h.

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_skyrme
 eos_had_skyrme ()
 Create a blank Skyrme EOS.
 
virtual ~eos_had_skyrme ()
 Destructor.
 
virtual int calc_temp_e (fermion &ne, fermion &pr, double temper, thermo &th)
 Equation of state as a function of densities at finite temperature.
 
virtual int calc_deriv_temp_e (fermion_deriv &ne, fermion_deriv &pr, double temper, thermo &th, thermo_np_deriv_helm &thd)
 Equation of state as a function of the densities at finite temperature and including second derivatives.
 
virtual int calc_e (fermion &ne, fermion &pr, thermo &lt)
 Equation of state as a function of densities at zero temperature.
 
virtual const char * type ()
 Return string denoting type ("eos_had_skyrme")
 
virtual double feoa (double nb)
 Calculate binding energy. More...
 
virtual double fmsom (double nb)
 Calculate effective mass. More...
 
virtual double fcomp_nuc (double nb)
 Calculate compressibility in nuclear (isospin-symmetric matter) More...
 
virtual double fesym (double nb, double alpha=0.0)
 Calculate symmetry energy. More...
 
virtual double fkprime_nuc (double nb)
 Skewness in nuclear (isospin-symetric) matter. More...
 
int check_landau (double nb, double m)
 Check the Landau parameters for instabilities. More...
 
void landau_nuclear (double n0, double m, double &f0, double &g0, double &f0p, double &g0p, double &f1, double &g1, double &f1p, double &g1p)
 Calculate the Landau parameters for nuclear matter. More...
 
void landau_neutron (double n0, double m, double &f0, double &g0, double &f1, double &g1)
 Calculate the Landau parameters for neutron matter. More...
 
template<class fermion_t >
void eff_mass (fermion_t &ne, fermion_t &pr, double &term, double &term2)
 Evaluate the effective masses for neutrons and protons.
 
int calpar (double gt0=-10.0, double gt3=70.0, double galpha=0.2, double gt1=2.0, double gt2=-1.0)
 Calculate $ t_0,t_1,t_2,t_3 $ and $ \alpha $ from the saturation properties. More...
 
void alt_params_set (double Crr00, double Crr10, double Crr0D, double Crr1D, double Crt0, double Crt1, double CrDr0, double CrDr1, double CrnJ0, double CrnJ1, double alpha2)
 Set using alternate parameterization. More...
 
void alt_params_get (double &Crr00, double &Crr10, double &Crr0D, double &Crr1D, double &Crt0, double &Crt1, double &CrDr0, double &CrDr1, double &CrnJ0, double &CrnJ1, double &alpha2)
 Get alternate parameterization. More...
 
void alt_params_saturation (double n0, double EoA, double K, double Ms_star, double a, double L, double Mv_star, double CrDr0, double CrDr1, double CrnJ0, double CrnJ1)
 Use the specified saturation properties and couplings and the function alt_params_set() to set the Skyrme coefficients. More...
 
- Public Member Functions inherited from o2scl::eos_had_temp_eden_base
virtual int calc_p (fermion &n, fermion &p, thermo &th)
 Equation of state as a function of the chemical potentials.
 
virtual 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.
 
- 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 double calc_temp_mun_e (double nn, double np, double T)
 Neutron chemical potential as a function of the densities.
 
virtual double calc_temp_mup_e (double nn, double np, double T)
 Proton chemical potential as a function of the densities.
 
virtual double calc_temp_nn_p (double mun, double mup, double T)
 Neutron density as a function of the chemical potentials.
 
virtual double calc_temp_np_p (double mun, double mup, double T)
 Proton density as a function of the chemical potentials.
 
double calc_fr (double nn, double np, double T)
 Compute the free energy as a function of the temperature and the densities.
 
virtual void f_number_suscept_T (double mun, double mup, double T, double &dPdnn, double &dPdnp, double &dPdpp)
 Compute the number susceptibilities as a function of the chemical potentials, $ \partial^2 P / \partial \mu_i \mu_j $ at a fixed temperature.
 
virtual void f_inv_number_suscept_T (double mun, double mup, double T, double &dednn, double &dednp, double &dedpp)
 Compute the 'inverse' number susceptibilities as a function of the densities, $ \partial^2 \varepsilon / \partial n_i n_j $ at a fixed temperature.
 
void check_en (fermion &n, fermion &p, double T, thermo &th, double &en_deriv, double &en_err)
 Check the entropy by computing the derivative numerically.
 
void check_mu_T (fermion &n, fermion &p, double T, thermo &th, double &mun_deriv, double &mup_deriv, double &mun_err, double &mup_err)
 Check the chemical potentials at finite temperature by computing the derivative numerically.
 
virtual void set_fermion_thermo (fermion_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_err (double nb, double delta, double &unc)
 Calculate symmetry energy of matter and its uncertainty in $ \mathrm{fm}^{-1} $. More...
 
virtual double fesym_slope (double nb, double delta=0.0)
 The symmetry energy slope parameter in $ \mathrm{fm}^{-1} $. More...
 
virtual double fesym_curve (double nb, double delta=0.0)
 The curvature of the symmetry energy in $ \mathrm{fm}^{-1} $.
 
virtual double fesym_skew (double nb, double delta=0.0)
 The skewness of the symmetry energy in $ \mathrm{fm}^{-1} $.
 
virtual double fesym_diff (double nb)
 Calculate symmetry energy of matter as energy of neutron matter minus the energy of nuclear matter in $ \mathrm{fm}^{-1} $. 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 in $ \mathrm{fm}^{-1} $ 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 (reduced) effective mass.
 
virtual double f_effm_prot (double nb, double delta=0.0)
 Proton (reduced) 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=1.0)
 Vector effective mass. More...
 
virtual double fn0 (double delta, double &leoa)
 Calculate saturation density using calc_e() More...
 
virtual void f_number_suscept (double mun, double mup, double &dPdnn, double &dPdnp, double &dPdpp)
 Compute the number susceptibilities as a function of the chemical potentials, $ \partial^2 P / \partial \mu_i \mu_j $.
 
virtual void f_inv_number_suscept (double mun, double mup, double &dednn, double &dednp, double &dedpp)
 Compute the 'inverse' number susceptibilities as a function of the densities, $ \partial^2 \varepsilon / \partial n_i n_j $.
 
virtual int saturation ()
 Calculates some of the EOS properties at the saturation density. More...
 
double calc_mun_e (double nn, double np)
 Compute the neutron chemical potential at fixed density. More...
 
double calc_ed (double nn, double np)
 Compute the energy density as a function of the nucleon densities.
 
double calc_pr (double nn, double np)
 Compute the pressure as a function of the nucleon chemical potentials.
 
double calc_mup_e (double nn, double np)
 Compute the proton chemical potential at fixed density. More...
 
double calc_nn_p (double mun, double mup)
 Compute the neutron density at fixed chemical potential. More...
 
double calc_np_p (double mun, double mup)
 Compute the proton density at fixed chemical potential. 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...
 
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...
 
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 check_mu (fermion &n, fermion &p, thermo &th, double &mun_deriv, double &mup_deriv, double &mun_err, double &mup_err)
 Check the chemical potentials by computing the derivatives numerically.
 
void check_den (fermion &n, fermion &p, thermo &th, double &nn_deriv, double &np_deriv, double &nn_err, double &np_err)
 Check the densities by computing the derivatives numerically.
 
- 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
 
- Public Attributes inherited from o2scl::eos_had_skyrme
double t0
 
double t1
 
double t2
 
double t3
 
double x0
 
double x1
 
double x2
 
double x3
 
double alpha
 
double a
 
double b
 
double W0
 Spin-orbit splitting (in $ \mathrm{fm}^{-1} $) More...
 
double b4
 Isoscalar spin-orbit term (in $ \mathrm{fm}^{-1} $)
 
double b4p
 Isovector spin-orbit term (in $ \mathrm{fm}^{-1} $)
 
std::string reference
 Bibliographic reference.
 
bool parent_method
 Use eos_had_base methods for saturation properties. More...
 
fermion_nonrel nrf
 Thermodynamics of non-relativistic fermions.
 
fermion_deriv_nr nrfd
 Thermodynamics of non-relativistic fermions with derivatives.
 
- Public Attributes inherited from o2scl::eos_had_temp_base
fermion_eff def_fet
 Default fermion thermodynamics object.
 
- Public Attributes inherited from o2scl::eos_had_base
double eoa
 Binding energy (without the rest mass) in $ \mathrm{fm}^{-1} $. More...
 
double comp
 Compression modulus in $ \mathrm{fm}^{-1} $.
 
double esym
 Symmetry energy in $ \mathrm{fm}^{-1} $.
 
double n0
 Saturation density in $ \mathrm{fm}^{-3} $.
 
double msom
 Effective mass (neutron)
 
double kprime
 Skewness in $ \mathrm{fm}^{-1} $.
 
bool err_nonconv
 If true, call the error handler if msolve() or msolve_de() does not converge (default true)
 
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...
 
fermion def_neutron
 The defaut neutron. More...
 
fermion def_proton
 The defaut proton. More...
 
- Public Attributes inherited from o2scl::eos_base
thermo def_thermo
 The default thermo object.
 
- Protected Member Functions inherited from o2scl::eos_had_skyrme
void hamiltonian_coeffs (double &ham1, double &ham2, double &ham3, double &ham4, double &ham5, double &ham6)
 Compute the coefficients of the potential energy which have unique dependence on the densities. More...
 
template<class fermion_t >
void base_thermo (fermion_t &ne, fermion_t &pr, double ltemper, thermo &locth, double term, double term2, double ham1, double ham2, double ham3, double ham4, double ham5, double ham6)
 Compute the base thermodynamic quantities. More...
 
template<class fermion_t >
void second_deriv (fermion_t &ne, fermion_t &pr, double ltemper, thermo &locth, thermo_np_deriv_helm &locthd, double term, double term2, double ham1, double ham2, double ham3, double ham4, double ham5, double ham6)
 Compute second derivatives of the free energy.
 
template<class fermion_t >
void check_input (fermion_t &ne, fermion_t &pr, double T)
 
int calparfun (size_t nv, const ubvector &x, ubvector &y)
 
int calparfun2 (size_t nv, const ubvector &x, ubvector &y)
 
- 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. More...
 
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. More...
 
- 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_skyrme
double fixn0
 
double fixeoa
 
double fixesym
 
double fixcomp
 
double fixmsom
 
- Protected Attributes inherited from o2scl::eos_had_temp_base
fermion_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.
 
fermion * neutron
 The neutron object.
 
fermion * proton
 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).