neost.eos package

Submodules

neost.eos.base module

class neost.eos.base.BaseEoS(crust='ceft-Hebeler', rho_t=200000000000000.0)[source]

Bases: object

Base class representing an equation of state object.

Parameters

crust: str

The name of the EoS crust model to use. Can be either ‘ceft-Hebeler’, ‘ceft-Drischler’, ‘ceft-Lynn’, ‘ceft-Tews’, ‘ceft-old’, ‘BPS’, or None if a tabulated EoS with a crust model already included is used.

rho_t: float

The transition density between the crust EOS and the high density parameterization in cgs.

Methods

update(eos_params, max_edsc=True)

Update the EoS object with a given set of parameters

get_eos_crust()

Construct the crust of the equation of state, with or without cEFT.

plot()

Plot the equation of state.

plot_massradius()

Plot the mass-radius curve of the equation of state.

CofE(E, epsgrid, presgrid)[source]
EofP(P, epsgrid, presgrid)[source]
EofRho(rho, rhos, eps)[source]
Mass_Radius(epscent, epscent_dm)[source]
PofE(E, epsgrid, presgrid)[source]
SLYfit(rho)[source]
ceft_band_func(rho, norm, min_norm, max_norm, min_index, max_index)[source]
edens(eps, rho, eos)[source]
f_chi_calc(epscent, epscent_dm)[source]

Method to calculate the ADM mass-fraction given the baryonic and ADM central densities, respectively.

Parameters

epscentfloat

Baryonic central energy density in cgs units for mass-density, i.e., divided by the speed of light squared.

epscent_dm: float

ADM central energy density in cgs units for mass-density, i.e., divided by the speed of light squared.

find_epsdm_cent(ADM_fraction, epscent)[source]

Method to calculate the ADM central energy density given the baryonic central energy density and ADM mass-fraction. Uses a wide array of different intervals of central energy densities to determine the ADM central energy density as a root finding problem.

Parameters

ADM_fraction: float

Given ADM mass-fraction as a percentage [%]

epscentfloat

Baryonic central energy density in cgs units for mass-density, i.e., divided by the speed of light squared.

find_max_edsc()[source]
get_BPS()[source]
get_eos_crust()[source]
get_minmax_edsc_chirp(chirp)[source]
plot(dm='None')[source]

Plot the EoS. If dm is not ‘None’ will include ADM contribution.

plot_massradius()[source]

Plot the mass-radius curve.

polytropic_func(rho, K, index)[source]
rhodens(rho, eps)[source]
update(eos_params, max_edsc=True)[source]

Method to update a given EoS object with specified parameters.

Parameters

eos_paramsdict

A dictionary containing the parameter values of the EoS model.

max_edsc: bool

If True, compute the maximum central energy density allowed by this set of parameters (default is True).

neost.eos.polytropes module

class neost.eos.polytropes.PolytropicEoS(crust, rho_t, adm_type='None', dm_halo=False, two_fluid_tidal=False)[source]

Bases: BaseEoS

Class representing a polytropic equation of state object.

Parameters

rho_t: float

The transition density between the crust EOS and the high density parameterization in cgs.

ceft: bool

If True a low-density cEFT parameterization is used.

ceft_method: str

The name of the cEFT calculations used at low density. Can be one of ‘Hebeler’, ‘Drischler’, ‘Lynn’, ‘Keller-N2LO’, ‘Keller-N3L0’, or ‘Tews’.

adm_type: str

The name of the ADM particle type. Can be ‘None’, ‘Bosonic’, or ‘Fermionic’

dm_halo: bool

If True, ADM halos will be allowed.

two_fluid_tidal: bool

If True, the two-fluid tidal deformability solver will be used. See TidalDef.py file

Methods

update(eos_params, max_edsc=True)

Update the EoS object with a given set of parameters

get_eos_crust()

Construct the crust of the equation of state, with or without cEFT.

get_eos()

Construct the high-density parameterization of the equation of state.

eos_core_pp()

Function to compute the polytropic equation of state parameterization.

add_adm_eos()

Function to compute the ADM EoS, whether it be bosonic or fermionic in nature.

plot()

Plot the equation of state.

plot_massradius()

Plot the mass-radius curve of the equation of state.

add_adm_eos()[source]
check_constraints()[source]
eos_core_pp(rho, P_t)[source]
get_eos()[source]
polytropic_func(rho, K, index)[source]

neost.eos.speedofsound module

class neost.eos.speedofsound.SpeedofSoundEoS(crust, rho_t, adm_type='None', dm_halo=False, two_fluid_tidal=False)[source]

Bases: BaseEoS

Class representing an equation of state object.

Parameters

eos_namestr

The name of the EOS parameterization. Should be one of ‘polytropes’, ‘speedofsound’ or ‘spectral’.

rho_t: float

The transition density between the crust EOS and the high density parameterization in cgs.

ceft: bool

If True a low-density cEFT parameterization is used.

ceft_method: str

The name of the cEFT calculations used at low density. Can be one of ‘Hebeler’, ‘Drischler’, ‘Lynn’, ‘Keller-N2LO’, ‘Keller-N3L0’, or ‘Tews’.

adm_type: str

The name of the ADM particle type. Can be ‘None’, ‘Bosonic’, or ‘Fermionic’

dm_halo: bool

If True, ADM halos will be allowed.

two_fluid_tidal: bool

If True, the two-fluid tidal deformability solver will be used. See TidalDef.py file

Methods

update(eos_params, max_edsc=True)

Update the EoS object with a given set of parameters

get_eos_crust()

Construct the crust of the equation of state, with or without cEFT.

get_eos()

Construct the high-density parameterization of the equation of state.

add_adm_eos()

Function to compute the ADM EoS, whether it be bosonic or fermionic in nature.

plot()

Plot the equation of state.

plot_massradius()

Plot the mass-radius curve of the equation of state.

Cs_model_core(x, norm, negative=-1e+100)[source]
Cs_model_total(x, norm)[source]
add_adm_eos()[source]
check_constraints()[source]
get_eos()[source]
logistic(x, norm=0.0, xt=5, s=1.0)[source]
match_func(x, cscrust)[source]