26 #ifndef O2SCL_SCHEMATIC_EOS_H 27 #define O2SCL_SCHEMATIC_EOS_H 31 #include <o2scl/eos_had_base.h> 33 #ifndef DOXYGEN_NO_O2NS 137 double alpha=n*(1.0-2.0*x);
138 double ret=(kpp*n*(5.0-2.0*
n0)*(n-
n0)*(n-
n0)+
140 kprime*n*(2.0*n*n-3.0*n*
n0+
n0*
n0)+3.0*pow(n0,3.0)*alpha*
141 (-1.0*a*pow(n/n0,2.0/3.0)*(-6.0+alpha)+
142 9*b*pow(n/n0,gamma)*(2.0+alpha*(-1.0+
gamma))))/
143 (486.0*n*pow(n0,4.0));
148 virtual const char *
type() {
return "eos_had_schematic"; }
152 #ifndef DOXYGEN_NO_O2NS virtual int set_kprime_zeroden()
Set kprime so that the energy per baryon of zero-density matter is zero.
double comp
Compression modulus in .
double a
The kinetic energy symmetry coefficient in inverse fm (default )
virtual int calc_e(fermion &ln, fermion &lp, thermo <h)
Equation of state as a function of density.
double n0
Saturation density in .
virtual int set_a_from_mstar(double u_msom, double mnuc)
Fix the kinetic energy symmetry coefficient from the reduced nucleon effective mass and the saturatio...
double gamma
The exponent of the high-density symmetry energy (unitless, default 1.0)
A hadronic EOS based on a function of the densities [abstract base].
double b
The potential energy symmetry coefficient in inverse fm (default )
Schematic hadronic equation of state.
virtual int set_kpp_zeroden()
Set kpp so that the energy per baryon of zero-density matter is zero.
double kprime
Skewness in .
virtual double baryon_suscep(double n, double x)
Return the baryon number susceptibility, in .
virtual double eoa_zeroden()
Return the energy per baryon of matter at zero density.
double kpp
The coefficient of a density to the fourth term in inverse fm (default 0)
double eoa
Binding energy (without the rest mass) in .
virtual const char * type()
Return string denoting type ("eos_had_schematic")