-
Realistic binary neutron star initial data with Elliptica
Authors:
Alireza Rashti,
Andrew Noe
Abstract:
This work introduces the Elliptica pseudo-spectral code for generating initial data of binary neutron star systems. Building upon the recent Elliptica code update, we can now construct initial data using not only piecewise polytropic equations of state, but also tabulated equations of state for these binary systems. Furthermore, the code allows us to endow neutron stars within the binary system wi…
▽ More
This work introduces the Elliptica pseudo-spectral code for generating initial data of binary neutron star systems. Building upon the recent Elliptica code update, we can now construct initial data using not only piecewise polytropic equations of state, but also tabulated equations of state for these binary systems. Furthermore, the code allows us to endow neutron stars within the binary system with spins. These spins can have a magnitude close to the mass shedding limit and can point in any direction.
△ Less
Submitted 1 July, 2024;
originally announced July 2024.
-
Scattering and dynamical capture of two black holes: synergies between numerical and analytical methods
Authors:
Simone Albanesi,
Alireza Rashti,
Francesco Zappa,
Rossella Gamba,
William Cook,
Boris Daszuta,
Sebastiano Bernuzzi,
Alessandro Nagar,
David Radice
Abstract:
We study initially unbound systems of two black holes using numerical relativity (NR) simulations performed with GR-Athena++. We focus on regions of the parameter space close to the transition from scatterings to dynamical captures, considering equal mass and spin-aligned configurations, as well as unequal mass and nonspinning ones. The numerical results are then used to validate the effective-one…
▽ More
We study initially unbound systems of two black holes using numerical relativity (NR) simulations performed with GR-Athena++. We focus on regions of the parameter space close to the transition from scatterings to dynamical captures, considering equal mass and spin-aligned configurations, as well as unequal mass and nonspinning ones. The numerical results are then used to validate the effective-one-body (EOB) model TEOBResumS-Dalí for dynamical captures and scatterings. We find good agreement for the waveform phenomenologies, scattering angles, mismatches, and energetics in the low energy regime ($E_0\lesssim 1.02\,M$). In particular, mismatches are typically below or around the $1\%$ level, with only a few cases -- corresponding to spinning binaries -- slightly above the $3\%$ threshold. We also discuss dynamical captures in the test-mass limit by solving numerically the Zerilli equation with the time domain code RWZHyp. The latter analysis provides valuable insights into both the analytical noncircular corrections of the EOB waveform and the integration of NR Weyl scalars.
△ Less
Submitted 30 May, 2024;
originally announced May 2024.
-
Adaptive mesh refinement in binary black holes simulations
Authors:
Alireza Rashti,
Maitraya Bhattacharyya,
David Radice,
Boris Daszuta,
William Cook,
Sebastiano Bernuzzi
Abstract:
We discuss refinement criteria for the Berger-Rigoutsos (block-based) refinement algorithm in our numerical relativity code GR-Athena++ in the context of binary black hole merger simulations. We compare three different strategies: the "box-in-box" approach, the "sphere-in-sphere" approach and a local criterion for refinement based on the estimation of truncation error of the finite difference sche…
▽ More
We discuss refinement criteria for the Berger-Rigoutsos (block-based) refinement algorithm in our numerical relativity code GR-Athena++ in the context of binary black hole merger simulations. We compare three different strategies: the "box-in-box" approach, the "sphere-in-sphere" approach and a local criterion for refinement based on the estimation of truncation error of the finite difference scheme. We extract and compare gravitational waveforms using the three different mesh refinement methods and compare their accuracy against a calibration waveform and demonstrate that the sphere-in-sphere approach provides the best strategy overall when considering computational cost and the waveform accuracy. Ultimately, we demonstrate the capability of each mesh refinement method in accurately simulating gravitational waves from binary black hole systems -- a crucial aspect for their application in next-generation detectors. We quantify the mismatch achievable with the different strategies by extrapolating the gravitational wave mismatch to higher resolution.
△ Less
Submitted 25 March, 2024; v1 submitted 8 December, 2023;
originally announced December 2023.
-
Analytically improved and numerical-relativity informed effective-one-body model for coalescing binary neutron stars
Authors:
Rossella Gamba,
Matteo Breschi,
Sebastiano Bernuzzi,
Alessandro Nagar,
William Cook,
Georgios Doulis,
Francesco Fabbri,
Néstor Ortiz,
Amit Poudel,
Alireza Rashti,
Wolfgang Tichy,
Maximiliano Ujevic
Abstract:
Gravitational wave astronomy pipelines rely on template waveform models for searches and parameter estimation purposes. For coalescing binary neutron stars (BNS), such models need to accurately reproduce numerical relativity (NR) up to merger, in order to provide robust estimate of the stars' equation of state - dependent parameters. In this work we present an improved version of the Effective One…
▽ More
Gravitational wave astronomy pipelines rely on template waveform models for searches and parameter estimation purposes. For coalescing binary neutron stars (BNS), such models need to accurately reproduce numerical relativity (NR) up to merger, in order to provide robust estimate of the stars' equation of state - dependent parameters. In this work we present an improved version of the Effective One Body (EOB) model $\tt TEOBResumS$ for gravitational waves from BNS systems. Building upon recent post-Newtonian calculations, we include subleading order tidal terms in the waveform multipoles and EOB metric potentials, as well as add up to 5.5PN terms in the gyro-gravitomagnetic functions entering the spin-orbit sector of the model. In order to further improve the EOB-NR agreement in the last few orbital cycles before merger, we introduce next-to-quasicircular corrections in the waveform -- informed by a large number of BNS NR simulations -- and introduce a new NR-informed parameter entering the tidal sector of our conservative dynamics. The performance of our model is then validated against 14 new eccentricity reduced simulations of unequal mass, spinning binaries with varying equation of state. A time-domain phasing analysis and mismatch computations demonstrate that the new model overall improves over the previous version of $\tt TEOBResumS$. Finally, we present a closed-form frequency domain representation of the (tidal) amplitude and phase of the new model. This representation accounts for mass-ratio, aligned spin and (resummed) spin-quadrupole effects in the tidal phase and -- within the calibration region -- it is faithful to the original model.
△ Less
Submitted 27 July, 2023;
originally announced July 2023.
-
The new discontinuous Galerkin methods based numerical relativity program Nmesh
Authors:
Wolfgang Tichy,
Liwei Ji,
Ananya Adhikari,
Alireza Rashti,
Michal Pirog
Abstract:
Interpreting gravitational wave observations and understanding the physics of astrophysical compact objects such as black holes or neutron stars requires accurate theoretical models. Here, we present a new numerical relativity computer program, called Nmesh, that has the design goal to become a next generation program for the simulation of challenging relativistic astrophysics problems such as bin…
▽ More
Interpreting gravitational wave observations and understanding the physics of astrophysical compact objects such as black holes or neutron stars requires accurate theoretical models. Here, we present a new numerical relativity computer program, called Nmesh, that has the design goal to become a next generation program for the simulation of challenging relativistic astrophysics problems such as binary black hole or neutron star mergers. In order to efficiently run on large supercomputers, Nmesh uses a discontinuous Galerkin method together with a domain decomposition and mesh refinement that parallelizes and scales well. In this work, we discuss the various numerical methods we use. We also present results of test problems such as the evolution of scalar waves, single black holes and neutron stars, as well as shock tubes. In addition, we introduce a new positivity limiter that allows us to stably evolve single neutron stars without an additional artificial atmosphere, or other more traditional limiters.
△ Less
Submitted 12 December, 2022;
originally announced December 2022.
-
Second release of the CoRe database of binary neutron star merger waveforms
Authors:
Alejandra Gonzalez,
Francesco Zappa,
Matteo Breschi,
Sebastiano Bernuzzi,
David Radice,
Ananya Adhikari,
Alessandro Camilletti,
Swami Vivekanandji Chaurasia,
Georgios Doulis,
Surendra Padamata,
Alireza Rashti,
Maximiliano Ujevic,
Bernd Brügmann,
William Cook,
Tim Dietrich,
Albino Perego,
Amit Poudel,
Wolfgang Tichy
Abstract:
We present the second data release of gravitational waveforms from binary neutron star merger simulations performed by the Computational Relativity (CoRe) collaboration. The current database consists of 254 different binary neutron star configurations and a total of 590 individual numerical-relativity simulations using various grid resolutions. The released waveform data contain the strain and the…
▽ More
We present the second data release of gravitational waveforms from binary neutron star merger simulations performed by the Computational Relativity (CoRe) collaboration. The current database consists of 254 different binary neutron star configurations and a total of 590 individual numerical-relativity simulations using various grid resolutions. The released waveform data contain the strain and the Weyl curvature multipoles up to $\ell=m=4$. They span a significant portion of the mass, mass-ratio,spin and eccentricity parameter space and include targeted configurations to the events GW170817 and GW190425. CoRe simulations are performed with 18 different equations of state, seven of which are finite temperature models, and three of which account for non-hadronic degrees of freedom. About half of the released data are computed with high-order hydrodynamics schemes for tens of orbits to merger; the other half is computed with advanced microphysics. We showcase a standard waveform error analysis and discuss the accuracy of the database in terms of faithfulness. We present ready-to-use fitting formulas for equation of state-insensitive relations at merger (e.g. merger frequency), luminosity peak, and post-merger spectrum.
△ Less
Submitted 28 March, 2023; v1 submitted 28 October, 2022;
originally announced October 2022.
-
High-accuracy high-mass ratio simulations for binary neutron stars and their comparison to existing waveform models
Authors:
Maximiliano Ujevic,
Alireza Rashti,
Henrique Gieg,
Wolfgang Tichy,
Tim Dietrich
Abstract:
The subsequent observing runs of the advanced gravitational-wave detector network will likely provide us with various gravitational-wave observations of binary neutron star systems. For an accurate interpretation of these detections, we need reliable gravitational-wave models. To test and to point out how existing models could be improved, we perform a set of high-resolution numerical-relativity s…
▽ More
The subsequent observing runs of the advanced gravitational-wave detector network will likely provide us with various gravitational-wave observations of binary neutron star systems. For an accurate interpretation of these detections, we need reliable gravitational-wave models. To test and to point out how existing models could be improved, we perform a set of high-resolution numerical-relativity simulations for four different physical setups with mass ratios $q$ = $1.25$, $1.50$, $1.75$, $2.00$, and total gravitational mass $M = 2.7M_\odot$ . Each configuration is simulated with five different resolutions to allow a proper error assessment. Overall, we find approximately 2nd order converging results for the dominant $(2,2)$, but also subdominant $(2,1)$, $(3,3)$, $(4,4)$ modes, while, generally, the convergence order reduces slightly for an increasing mass ratio. Our simulations allow us to validate waveform models, where we find generally good agreement between state-of-the-art models and our data, and to prove that scaling relations for higher modes currently employed for binary black hole waveform modeling also apply for the tidal contribution. Finally, we also test if the current NRTidal model to describe tidal effects is a valid description for high-mass ratio systems. We hope that our simulation results can be used to further improve and test waveform models in preparation for the next observing runs.
△ Less
Submitted 18 February, 2022;
originally announced February 2022.
-
Elliptica: a new pseudo-spectral code for the construction of initial data
Authors:
Alireza Rashti,
Francesco Maria Fabbri,
Bernd Brügmann,
Swami Vivekanandji Chaurasia,
Tim Dietrich,
Maximiliano Ujevic,
Wolfgang Tichy
Abstract:
Numerical studies of the dynamics of gravitational systems, e.g., black hole-neutron star systems, require physical and constraint-satisfying initial data. In this article, we present the newly developed pseudo-spectral code Elliptica, an infrastructure for construction of initial data for various binary and single gravitational systems of all kinds. The elliptic equations under consideration are…
▽ More
Numerical studies of the dynamics of gravitational systems, e.g., black hole-neutron star systems, require physical and constraint-satisfying initial data. In this article, we present the newly developed pseudo-spectral code Elliptica, an infrastructure for construction of initial data for various binary and single gravitational systems of all kinds. The elliptic equations under consideration are solved on a single spatial hypersurface of the spacetime manifold. Using coordinate maps, the hypersurface is covered by patches whose boundaries can adapt to the surface of the compact objects. To solve elliptic equations with arbitrary boundary condition, Elliptica deploys a Schur complement domain decomposition method with a direct solver. In this version, we use cubed sphere coordinate maps and the fields are expanded using Chebyshev polynomials of the first kind. Here, we explain the building blocks of Elliptica and the initial data construction algorithm for a black hole-neutron star binary system. We perform convergence tests and evolve the data to validate our results. Within our framework, the neutron star can reach spin values close to breakup with arbitrary direction, while the black hole can have arbitrary spin with dimensionless spin magnitude $\sim 0.8$.
△ Less
Submitted 7 May, 2022; v1 submitted 29 September, 2021;
originally announced September 2021.
-
High-accuracy simulations of highly spinning binary neutron star systems
Authors:
Reetika Dudi,
Tim Dietrich,
Alireza Rashti,
Bernd Bruegmann,
Jan Steinhoff,
Wolfgang Tichy
Abstract:
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 nu…
▽ 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 $χ\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.
△ Less
Submitted 23 August, 2021;
originally announced August 2021.
-
Constructing Binary Neutron Star Initial Data with High Spins, High Compactness, and High Mass-Ratios
Authors:
Wolfgang Tichy,
Alireza Rashti,
Tim Dietrich,
Reetika Dudi,
Bernd Brügmann
Abstract:
The construction of accurate and consistent initial data for various binary parameters is a critical ingredient for numerical relativity simulations of the compact binary coalescence. In this article, we present an upgrade of the pseudospectral SGRID code, which enables us to access even larger regions of the binary neutron star parameter space. As a proof of principle, we present a selected set o…
▽ More
The construction of accurate and consistent initial data for various binary parameters is a critical ingredient for numerical relativity simulations of the compact binary coalescence. In this article, we present an upgrade of the pseudospectral SGRID code, which enables us to access even larger regions of the binary neutron star parameter space. As a proof of principle, we present a selected set of first simulations based on initial configurations computed with the new code version. In particular, we simulate two millisecond pulsars close to their breakup spin, highly compact neutron stars with masses at about $98\%$ of the maximum supported mass of the employed equation of state, and an unequal mass systems with mass ratios even outside the range predicted by population synthesis models ($q = 2.03$). The discussed code extension will help us to simulate previously unexplored binary configurations. This is a necessary step to construct and test new gravitational wave approximants and to interpret upcoming binary neutron star merger observations. When we construct initial data, one has to specify various parameters, such as a rotation parameter for each star. Some of these parameters do not have direct physical meaning, which makes comparisons with other methods or models difficult. To facilitate this, we introduce simple estimates for the initial spin, momentum, mass, and center of mass of each individual star.
△ Less
Submitted 21 October, 2019;
originally announced October 2019.