amotornenko / tovsolver Goto Github PK

1. I have compiled your codes using these lines

m_arr = []
R_arr = []

for dens_c in np.logspace(-0,3.7,350):
R, M, prof = tov_s.solve(dens_c, rmax=50e5, dr=100)
m_arr.append(M)
R_arr.append(R)
print(R, M, dens_c) # extra I have written to check which central density gives the maximum mass of the neutron star.

I found that the radius and central density corresponded to the maximum mass 2.1498973233892715 M_sun are

12.004 km 1076.692506643068 MeV/fm^3

2. Now, I have fixed only this central density to check whether i am getting the same or not. But it gives me wrong values.

R, M, prof = tov_s.solve(1076.692506643068, rmax=50e5, dr=100) # 2.1498973233892715 M_sun, 12.004 km

12.004 km 4.274931084233547 M_sun

If you compare the masses there is huge differences 4.274931084233547-2.1498973233892715 = 2.125033761 M_sun. But the radius is same.

3. For arbitrary star say mass 1.3820322317287081 M_sun, the radius and central density are 13.144 km and 386.20641587213856 Mev/fm^3

R, M, prof = tov_s.solve(386.20641587213856, rmax=50e5, dr=100) # 1.3820322317287081, 13.144

print(max(prof[0]/1e5), max(prof[3]/1e33)) gives

13.144 km 2.748081260697474 M_sun

Please help me in this regard.

Hello,
When I used my own EOS for calculation, I found that there was an exception in it, as shown in the figure, but I did not quite understand this exception and was not clear about its impact on the calculation results. I hope I can get your help.
Thank you.

Hello,
I have installed the TOVsolver code as instructed.

Now I want to use my EoS to solve the TOV equation and calculate Mass, radius, etc.
Where do I need to save my EoS and how do I then run tov.py

Please respond back.
I will be waiting

For example, why can't the left side of the neutron star mass-radius relation curve continue to be drawn?

Hi ,
Can u tell me why this issue is coming