Class eos_sn_ls (o2scl)

O2scl : Class List

class eos_sn_ls : public o2scl::eos_sn_base

The Lattimer-Swesty supernova EOS.

This class is experimental.

Note that the tables on this website are different than what is generated from the LS Fortran code. See eos_sn_oo to read O’Connor and Ott’s tables generated from the LS Fortran code.

The four models are

  • LS (K=370, Sv=31)

  • SKI’ (K=371, Sv=30.4)

  • SKa (K=263, Sv=34.5)

  • SKM* (K=217, Sv=31.4)

See also the documentation at ref eos_sn_base and the general description in the Finite-temperature Equation of State Tables section of the User’s guide.

See also [Lattimer91] and [Lattimer85].

Todo

  • In class eos_sn_ls: There are still a few points for which the electron/photon EOS seems to be off, but this may be the result of small inaccuracies from finite-differencing the LS table.

Note

Does not contain the Lattimer-Swesty EOS, only provides some code to manipulate it. This class is designed to be used with the files ls.dat, sk1.dat, ska.dat and skm.dat as provided on Jim Lattimer’s website, http://www.astro.sunysb.edu/lattimer/EOS/main.html .

Note

In the original table, the full internal energy per baryon (data section 4 of 26) is apparently based on a rest mass of \( Y_e m_p + (1-Y_e) m_n \), while the baryon part of the internal energy per baryon (data section 13 of 26) is based on a rest mass of \( m_n \). This means that

\[ E - E_{\mathrm{int}} = E_{\mathrm{eg}} - Y_e (m_n - m_p) \]
where \( E_{\mathrm{eg}} \) is the energy per baryon of electrons and photons. In order to keep things consistent with the other EOS tables, when the EOS table is loaded, eos_sn_base::Eint is rescaled to a rest mass of \( Y_e m_p + (1-Y_e) m_n \) .

Additional data included in this EOS

tensor_grid3 &fill

Filling factor for nuclei.

tensor_grid3 &nb_in

Baryon number density inside nuclei in \( \mathrm{fm}^{-3} \).

tensor_grid3 &dPdn

Derivative of pressure with respect to baryon density.

tensor_grid3 &dPdT

Derivative of pressure with respect to temperature.

tensor_grid3 &dPdY

Derivative of pressure with respect to electron fraction.

tensor_grid3 &dsdT

Derivative of entropy with respect to temperature.

tensor_grid3 &dsdY

Derivative of entropy with respect to electron fraction.

tensor_grid3 &Nskin

Number of neutrons in skin.

tensor_grid3 &nb_out

Baryon density outside nuclei in \( \mathrm{fm}^{-3} \).

tensor_grid3 &x_out

Proton fraction outside nuclei.

tensor_grid3 &mu

Out of whackness parameter, \( \mu_n-\mu_p-\mu_e+1.293~\mathrm{MeV} \), in MeV.

inline eos_sn_ls()
virtual void load(std::string fname, size_t mode)

Load table from filename fname.

virtual double check_eg()

Check electrons and photons.

This checks that the electron and photon thermodynamics generated by is consistent with the data in E, Eint, F, Fint, P, Pint, S, and Sint.