- #1

Stella.Physics

- 63

- 13

I am working on this paper by Lattimer on binary systems:

https://arxiv.org/abs/astro-ph/0002203

I am trying to reproduce Figure 11 but something seems to go wrong.

Here's what I have done so far:

we set ##q=\frac{m_{ns}}{M_{bh}}## the ratio of the neutron star mass and the black hole mass.

then we set the radius

##\frac{R_{l}}{a} = 0.46\bigg(\frac{q}{1+q}\bigg)^{1/3}## with a being the semi major axis that follows

## a<2.8 \bigg[ \frac{m_{ns} m_{bh} (m_{ns}+m_{bh})}{〖M^3_ {\odot}〗} \bigg]^{\frac{1}{4}} R_{\odot}##

The first thing to plot is the neutron star radius versus the neutron star mass.

Then we must plot ##\frac{d \ln{R}}{d \ln{m_{ns}}} ## versus ## m_{ns}## which derivative simply leads to ##\frac{m_{ns}}{R} \frac{dR}{dm_{ns}}##

I've been trying to plot these functions for days but what I get is far from the desired results...

Any help would be very appreciated :)