Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userIncreasing the Accuracy of Binary Neutron Star Simulations with an improved Vacuum Treatment
89
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
Numerical relativity simulations are essential to study the last stages of the binary neutron star coalescence. Unfortunately, for stable simulations there is the need to add an artificial low-density atmosphere. Here we discuss a new framework in which we can effectively set the density surrounding the neutron stars to zero to ensure a more accurate simulation. We test our method with a number of single star test cases and for an equal mass binary neutron star simulation. While the bulk motion of the system is not influenced, and hence, there is no improvement with respect to the emitted gravitational-wave signal, we find that the new approach is superior with respect to mass conservation and it allows a much better tracking of outward moving material. This will allow a more accurate simulation of the ejected material and supports the interpretation of present and future multi-messenger observations with more accurate numerical relativity simulations.
rate research
Read More
With an increasing number of expected gravitational-wave detections of binary neutron star mergers, it is essential that gravitational-wave models employed for the analysis of observational data are able to describe generic compact binary systems. This includes systems in which the individual neutron stars are millisecond pulsars for which spin effects become essential. In this work, we perform numerical-relativity simulations of binary neutron stars with aligned and anti-aligned spins within a range of dimensionless spins of $chi sim [-0.28,0.58]$. The simulations are performed with multiple resolutions, show a clear convergence order and, consequently, can be used to test existing waveform approximants. We find that for very high spins gravitational-wave models that have been employed for the interpretation of GW170817 and GW190425 are not capable of describing our numerical-relativity dataset. We verify through a full parameter estimation study in which clear biases in the estimate of the tidal deformability and effective spin are present. We hope that in preparation of the next gravitational-wave observing run of the Advanced LIGO and Advanced Virgo detectors our new set of numerical-relativity data can be used to support future developments of new gravitational-wave models.
We perform binary neutron star merger simulations using a newly derived set of finite-temperature equations of state in the Brueckner-Hartree-Fock approach. We point out the important and opposite roles of finite temperature and rotation for stellar stability and systematically investigate the gravitational-wave properties, matter distribution, and ejecta properties in the postmerger phase for the different cases. The validity of several universal relations is also examined and the most suitable EOSs are identified.
We present the first set of numerical relativity simulations of binary neutron mergers that include spin precession effects and are evolved with multiple resolutions. Our simulations employ consistent initial data in general relativity with different spin configurations and dimensionless spin magnitudes $sim 0.1$. They start at a gravitational-wave frequency of $sim392$~Hz and cover more than $1$ precession period and about 15 orbits up to merger. We discuss the spin precession dynamics by analyzing coordinate trajectories, quasi-local spin measurements, and energetics, by comparing spin aligned, antialigned, and irrotational configurations. Gravitational waveforms from different spin configuration are compared by calculating the mismatch between pairs of waveforms in the late inspiral. We find that precession effects are not distinguishable from nonprecessing configurations with aligned spins for approximately face-on binaries, while the latter are distinguishable from a nonspinning configurations. Spin precession effects are instead clearly visible for approximately edge-on binaries. For the parameters considered here, precession does not significantly affect the characteristic postmerger gravitational-wave frequencies nor the mass ejection. Our results pave the way for the modeling of spin precession effects in the gravitational waveform from binary neutron star events.
High-accuracy numerical simulations of merging neutron stars play an important role in testing and calibrating the waveform models used by gravitational wave observatories. Obtaining high-accuracy waveforms at a reasonable computational cost, however, remains a significant challenge. One issue is that high-order convergence of the solution requires the use of smooth evolution variables, while many of the equations of state used to model the neutron star matter have discontinuities, typically in the first derivative of the pressure. Spectral formulations of the equation of state have been proposed as a potential solution to this problem. Here, we report on the numerical implementation of spectral equations of state in the Spectral Einstein Code. We show that, in our code, spectral equations of state allow for high-accuracy simulations at a lower computational cost than commonly used `piecewise polytrope equations state. We also demonstrate that not all spectral equations of state are equally useful: different choices for the low-density part of the equation of state can significantly impact the cost and accuracy of simulations. As a result, simulations of neutron star mergers present us with a trade-off between the cost of simulations and the physical realism of the chosen equation of state.
We present results from three-dimensional general relativistic simulations of binary neutron star coalescences and mergers using public codes. We considered equal mass models where the baryon mass of the two Neutron Stars (NS) is $1.4M_{odot}$, described by four different equations of state (EOS) for the cold nuclear matter (APR4, SLy, H4, and MS1; all parametrized as piecewise polytropes). We started the simulations from four different initial interbinary distances ($40, 44.3, 50$, and $60$ km), including up to the last 16 orbits before merger. That allows to show the effects on the gravitational wave phase evolution, radiated energy and angular momentum due to: the use of different EOSs, the orbital eccentricity present in the initial data and the initial separation (in the simulation) between the two stars. Our results show that eccentricity has a major role in the discrepancy between numerical and analytical waveforms until the very last few orbits, where tidal effects and missing high-order post-Newtonian coefficients also play a significant role. We test different methods for extrapolating the gravitational wave signal extracted at finite radii to null infinity. We show that an effective procedure for integrating the Newman-Penrose $psi_4$ signal to obtain the gravitational wave strain $h$ is to apply a simple high-pass digital filter to $h$ after a time domain integration, where only the two physical motivated integration constants are introduced. That should be preferred to the more common procedures of introducing additional integration constants, integrating in the frequency domain or filtering $psi_4$ before integration.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
