The Buchdahl EOS for the TOV solver. More...
#include <eos_tov.h>
The Buchdahl EOS is
which can be inverted to give
Physical solutions are obtained only for (ensuring that the argument to the square root is positive) and
(ensuring that the EOS is not acausal).
The baryon chemical potential is
The baryon density is
In the case that one assumes and
, the baryon density can be simplified to
c.f. Eq. 10 in Lattimer01.
The mass-radius curve is the solution of the equation
where and
(which has units of inverse km) is defined by
The central pressure and energy density are
To obtain energy density and pressure profiles can be obtained by solving
for the new coordinate where
is defined by
Using these, the profiles are
and
Based on Lattimer01 .
Public Member Functions | |
void | set_baryon_density (double nb, double ed) |
Set the baryon density. | |
virtual double | ed_from_pr (double pr) |
From the pressure, return the energy density. | |
virtual double | pr_from_ed (double ed) |
From the energy density, return the pressure. | |
virtual double | nb_from_ed (double ed) |
From the energy density, return the baryon density. | |
virtual double | nb_from_pr (double pr) |
From the pressure, return the baryon density. | |
virtual double | ed_from_nb (double nb) |
From the baryon density, return the energy density. | |
virtual double | pr_from_nb (double nb) |
From the baryon density, return the pressure. | |
virtual void | ed_nb_from_pr (double pr, double &ed, double &nb) |
Given the pressure, produce the energy and number densities. More... | |
virtual double | rad_from_gm (double gm) |
Given the gravitational mass, compute the radius. | |
virtual double | ed_from_r_gm (double r, double gm) |
Compute the energy density at radius r for a fixed gravitational mass. More... | |
virtual double | pr_from_r_gm (double r, double gm) |
Compute the pressure at radius r for a fixed gravitational mass. More... | |
![]() | |
bool | has_baryons () |
Return true if a baryon density is available. | |
void | check_nb (double &avg_abs_dev, double &max_abs_dev) |
Check that the baryon density is consistent with the ![]() | |
Public Attributes | |
double | Pstar |
The parameter with units of pressure in units of solar masses per km cubed (default value ![]() | |
![]() | |
int | verbose |
Control for output (default 1) | |
Protected Member Functions | |
virtual double | solve_rad (double rad, double gm) |
Solve for the radius at fixed gravitational mass. | |
virtual double | solve_rp (double rp, double r, double gm, double rad) |
Solve for ![]() ![]() | |
Protected Attributes | |
double | nb1 |
The baryon density at ed1 . | |
double | ed1 |
The energy density for which the baryon density is known. | |
double | pr1 |
The pressure at ed1 . | |
root_brent_gsl | rbg |
Solver. | |
![]() | |
bool | baryon_column |
Set to true if the baryon density is provided in the EOS (default false) | |
|
virtual |
|
virtual |
If the baryon density is not specified, it should be set to zero or baryon_column should be set to false
Implements o2scl::eos_tov.
|
virtual |
Documentation generated with Doxygen. Provided under the
GNU Free Documentation License (see License Information).