Determination of the neutron star Love number. More...
#include <tov_love.h>
Public Types | |
typedef boost::numeric::ublas::vector< double > | ubvector |
Public Member Functions | |
void | set_ODE (o2scl::ode_iv_solve<> &ois_new) |
Set ODE integrator. | |
void | calc_y (double &yR, double &beta, double &k2, double &lambda_km5, double &lambda_cgs, bool tabulate=false) |
Compute the love number using y. | |
void | calc_H (double &yR, double &beta, double &k2, double &lambda_km5, double &lambda_cgs) |
Compute the love number using H. | |
Public Attributes | |
o2scl::table_units | results |
A table containing the solution to the differential equation(s) | |
double | eps |
The first radial point in ![]() | |
std::shared_ptr< o2scl::table_units<> > | tab |
Pointer to the input profile. | |
Protected Member Functions | |
int | y_derivs (double r, size_t nv, const ubvector &vals, ubvector &ders) |
The derivative ![]() | |
int | H_derivs (double r, size_t nv, const ubvector &vals, ubvector &ders) |
The derivatives ![]() ![]() | |
double | eval_k2 (double beta, double yR) |
Compute ![]() | |
Protected Attributes | |
o2scl::ode_iv_solve | def_ois |
The default ODE integrator. | |
o2scl::ode_iv_solve * | oisp |
The ODE integrator. | |
double | schwarz_km |
Schwarzchild radius in km (set in constructor) | |
We use but keep factors of
, which has units
.
Following the notation in Postnikov10, define the function , which is the solution of
where (now surpressing the dependence on ),
which has units of ,
and
which has units of . The boundary conditions on
are that
and
for an arbitrary constant
(
is chosen to be equal to 1). Internally,
and
are stored in units of
.
From this we can define another (unitless) function , which obeys
with boundary condition is . Solving for
,
Define . This form for
is specified in y_derivs() .
The unitless quantity (the Love number) is defined by (this is the expression from Postnikov10 )
Hinderer10 writes the differential equation for in a slightly different (but equivalent) form,
This is the form given in H_derivs() .
The tidal deformability is then
and has units of or can be converted to
.
It is assumed that tab has been specified before-hand and has (at least) the following columns
ed
energy density in units of pr
pressure in units of cs2
sound speed squared (unitless)gm
gravitational mass in r
radius in Definition at line 140 of file tov_love.h.
|
protected |
Used in both tov_love::calc_y() and tov_love::calc_H().
Documentation generated with Doxygen. Provided under the
GNU Free Documentation License (see License Information).