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 (default 0.02)  
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 and .  
double  eval_k2 (double beta, double yR) 
Compute using the analytic expression. More...  
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 beforehand 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).