26 #ifndef O2SCL_SKYRME_EOS_H
27 #define O2SCL_SKYRME_EOS_H
33 #include <o2scl/constants.h>
34 #include <o2scl/mroot.h>
35 #include <o2scl/eos_had_base.h>
36 #include <o2scl/part.h>
37 #include <o2scl/fermion_nonrel.h>
39 #ifndef DOXYGEN_NO_O2NS
250 double t0,t1,t2,t3,x0,x1,x2,x3,alpha,a,b;
284 virtual double feoa(
double nb);
292 virtual double fmsom(
double nb);
303 virtual double fcomp(
double nb);
319 virtual double fesym(
double nb,
double alpha=0.0);
329 virtual double fkprime(
double nb);
363 int calpar(
double gt0=-10.0,
double gt3=70.0,
double galpha=0.2,
364 double gt1=2.0,
double gt2=-1.0);
396 double &f0,
double &g0,
double &f0p,
397 double &g0p,
double &f1,
double &g1,
398 double &f1p,
double &g1p);
410 double &f1,
double &g1);
413 virtual const char *
type() {
return "eos_had_skyrme"; }
441 (
double Crr00,
double Crr10,
double Crr0D,
double Crr1D,
442 double Crt0,
double Crt1,
double CrDr0,
double CrDr1,
443 double CrnJ0,
double CrnJ1,
double alpha2);
450 (
double &Crr00,
double &Crr10,
double &Crr0D,
double &Crr1D,
451 double &Crt0,
double &Crt1,
double &CrDr0,
double &CrDr1,
452 double &CrnJ0,
double &CrnJ1,
double &alpha2);
454 #ifndef DOXYGEN_NO_O2NS
465 double fixn0, fixeoa, fixesym, fixcomp, fixmsom;
472 #ifndef DOXYGEN_NO_O2NS
fermion_nonrel nrf
Thermodynamics of non-relativistic fermions.
std::string reference
Bibliographic reference.
virtual ~eos_had_skyrme()
Destructor.
virtual const char * type()
Return string denoting type ("eos_had_skyrme")
virtual int calc_e(fermion &ne, fermion &pr, thermo <)
Equation of state as a function of density.
void landau_neutron(double n0, double m, double &f0, double &g0, double &f1, double &g1)
Calculate the Landau parameters for neutron matter.
bool mu_at_zero_density
If true, compute the chemical potentials even at zero density.
double W0
Spin-orbit splitting (in )
A hadronic EOS at finite temperature based on a function of the densities [abstract base]...
double n0
Saturation density.
virtual double fcomp(double nb)
Calculate compressibility.
virtual double feoa(double nb)
Calculate binding energy.
virtual int calc_temp_e(fermion &ne, fermion &pr, double temper, thermo &th)
Equation of state as a function of densities.
virtual double fmsom(double nb)
Calculate effective mass.
virtual double fkprime(double nb)
skewness
int check_landau(double nb, double m)
Check the Landau parameters for instabilities.
eos_had_skyrme()
Create a blank Skyrme EOS.
virtual double fesym(double nb, double alpha=0.0)
Calculate symmetry energy.
Skyrme hadronic equation of state.
int calpar(double gt0=-10.0, double gt3=70.0, double galpha=0.2, double gt1=2.0, double gt2=-1.0)
Calculate and from the saturation properties.
bool parent_method
Use eos_had_base methods for saturation properties.
double b4
Isoscalar spin-orbit term (in )
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.
double b4p
Isovector spin-orbit term (in )
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.
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.