Equation of State Sub-Library: Version 0.910
Public Member Functions | Data Fields
schematic_eos Class Reference

Schematic hadronic equation of state. More...

#include <schematic_eos.h>

Inheritance diagram for schematic_eos:
hadronic_eos_eden hadronic_eos eos

Detailed Description

A schematic equation of state defined by the energy density:

\[ \epsilon = n_n m_n + n_p m_p + n \left\{ eoa+\frac{comp}{18}(n/n0-1)^2+ \frac{kprime}{162}(n/n0-1)^3+ \frac{kpp}{1944}(n/n0-1)^4+(1- 2 x)^2 \left[a \left(\frac{n}{n0}\right)^{2/3}+ b \left(\frac{n}{n0}\right)^{\gamma} \right] \right\} \]

Symmetry energy at nuclear matter density is a+b.

Note that it doesn't really matter what kind of particle object is used, since the calc_e() function doesn't use any of the particle thermodynamics functions.

Definition at line 52 of file schematic_eos.h.

Public Member Functions

virtual int calc_e (fermion &ln, fermion &lp, thermo &lth)
 Equation of state as a function of density.
virtual int set_kprime_zeroden ()
 Set kprime so that the energy per baryon of zero-density matter is zero.
virtual int set_kpp_zeroden ()
 Set kpp so that the energy per baryon of zero-density matter is zero.
virtual int set_a_from_mstar (double u_msom, double mnuc)
 Fix the kinetic energy symmetry coefficient from the nucleon effective mass and the saturation density.
virtual double eoa_zeroden ()
 Return the energy per baryon of matter at zero density.
virtual const char * type ()
 Return string denoting type ("schematic_eos")

Data Fields

double a
 The kinetic energy symmetry coefficient in inverse fm (default 17/hc)
double b
 The potential energy symmetry coefficient in inverse fm (default 13/hc)
double kpp
 The coefficient of a density to the fourth term (default 0)
double gamma
 The exponent of the high-density symmetry energy (default 1.0)

Member Function Documentation

virtual int schematic_eos::set_a_from_mstar ( double  u_msom,
double  mnuc 
) [inline, virtual]

This assumes the nucleons are non-relativistic and that the neutrons and protons have equal mass. The relativistic corrections are around 1 part in $ 10^{6} $.

Todo:
This was computed in schematic_sym.nb, which might be added to the documentation?

Definition at line 107 of file schematic_eos.h.

virtual double schematic_eos::eoa_zeroden ( ) [inline, virtual]

This is inaccessible from calc_e() so is available separately here. Using set_kprime_zeroden() or set_kpp_zeroden() will fix kprime or kpp (respectively) to ensure that this is zero.

The result provided here does not include the nucleon mass and is given in $ \mathrm{fm}^{-1} $.

Definition at line 122 of file schematic_eos.h.


Field Documentation

The default value corresponds to an effective mass of about 0.7.

Definition at line 62 of file schematic_eos.h.


The documentation for this class was generated from the following file:
 All Data Structures Namespaces Files Functions Variables Typedefs Enumerations Friends

Documentation generated with Doxygen. Provided under the GNU Free Documentation License (see License Information).

Get Object-oriented Scientific Computing
Lib at SourceForge.net. Fast, secure and Free Open Source software
downloads.