[Users] mass estimation of neutron star

[林家暉]

Dear Roland,
First thanks for the explanation of notions and I would check the paper of GRHydro.

To my understanding, "sum" output multiply dxdydz is the total mass of the neutron star(s) in the simulation.Is this correct ? In my case which is gallery example,  initial sum = 0.00012605 and final sum is 0.00012701 so that the mass increases from the initial binary neutron star to the final merged neutron star is (0.00012701- 0.00012605)*18*18*18*4=0.0224 M.
This is somehow different to what I thought. I originally think that during the inspiral and merging . some mass would outflow due to tidal force , so the total mass of neutron star would decrease. Contrast to this, the mass gained. IS this means accretion mass is larger than the summation of outflow mass and energy radiated away by the gravitational wave ?

Thank you.

Best regards,

Chia-Hui

________________________________
寄件者: Roland Haas <rhaas at illinois.edu>
寄件日期: 2018年9月25日 上午 01:22:16
收件者: 林家暉
副本: Einstein Toolkit Users
主旨: Re: [Users] mass estimation of neutron star

Hello Chia-Hui,

> Thanks for your information. However I have further questions:
> 1.In your expression :\int \sqrt{\gamma} \rho W d^3x ,is the gamma Lorentz factor ? And what is W?
Uhm, gamma is the determinant of the three metric, W is the Lorentz
factor. The equation was for illustration only, I would definitely not
take it very seriously. See eg the GRHydro paper:
https://arxiv.org/pdf/1304.5544 eq 10 for the definitions of the quantities.

> 2. For baryon mass case.What does the cell mass( "sum"*dx*dy*dz ) mean? In my case of binary neutron star, dx=dy=dz=18 and initial sum = 0.000126, so the cell mass=0.735 which is still much less than the initial mass of each neutron star(1.45 Ms).What is the physical meaning of the cell mass ?

What you call "cell mass" is the "cell volume". You also have to take
any symmetries into account. Eg if you used both rotation180 symmetry
and z symmetry then you are only evolving 1/4 of the physical volume
and thus have to multiply the number by 4. This then gives you the
expected result. This multiplication is needed because (for technical
reason) Carpet computes the "sum" reduction assuming a cell volume of 1.


> What I want to get is the mass lost during merge, which is 2*initial neutron star mass - final merged neutron star mass (for the equal mass binary). My first thought is integrating the density in rho.xy.h5 and assume spherical symmetry to calculate each mass.

The density in rho.xy.h5 is not sufficient since it is only the
quantity "rho" above so you are missing eg the Lorentz factor W and the
metric determinant "gamma". Also rho.xy.h5 contains only data in the xy
plane and not the full 3 volume (which would be in rho.xyz.file_*.h5).
Finally you would have to manually take care of not double counting in
regions where finer grid overlaps a coarser grid. It is much easier to
rely on Cactus to do this for you via the "sum" scalar output.

Yours,
Roland

--
My email is as private as my paper mail. I therefore support encrypting
and signing email messages. Get my PGP key from http://pgp.mit.edu .
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://lists.einsteintoolkit.org/pipermail/users/attachments/20180925/cf7c927f/attachment-0001.html