- #1
Stella.Physics
- 63
- 13
Hello to all. I've been trying to plot some graphs on Mathematica but I have faced some troubles so far.
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 :)
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 :)