Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userEvolutions of Nearly Maximally Spinning Black Hole Binaries Using the Moving Puncture Approach
126
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
We demonstrate that numerical relativity codes based on the moving punctures formalism are capable of evolving nearly maximally spinning black hole binaries. We compare a new evolution of an equal-mass, aligned-spin binary with dimensionless spin chi=0.99 using puncture-based data with recent simulations of the SXS Collaboration. We find that the overlap of our new waveform with the published results of the SXS Collaboration is larger than 0.999. To generate our new waveform, we use the recently introduced HiSpID puncture data, the CCZ4 evolution system, and a modified lapse condition that helps keep the horizon radii reasonably large.
rate research
Read More
We evolve a binary black hole system bearing a mass ratio of $q=m_1/m_2=2/3$ and individual spins of $S^z_1/m_1^2=0.95$ and $S^z_2/m_2^2=-0.95$ in a configuration where the large black hole has its spin antialigned with the orbital angular momentum, $L^z$, and the small black hole has its spin aligned with $L^z$. This configuration was chosen to measure the maximum recoil of the remnant black hole for nonprecessing binaries. We find that the remnant black hole recoils at 500km/s, the largest recorded value from numerical simulations for aligned spin configurations. The remnant mass, spin, and gravitational waveform peak luminosity and frequency also provide a valuable point in parameter space for source modeling.
Numerical relativity is central to the investigation of astrophysical sources in the dynamical and strong-field gravity regime, such as binary black hole and neutron star coalescences. Current challenges set by gravitational-wave and multi-messenger astronomy call for highly performant and scalable codes on modern massively-parallel architectures. We present GR-Athena++, a general-relativistic, high-order, vertex-centered solver that extends the oct-tree, adaptive mesh refinement capabilities of the astrophysical (radiation) magnetohydrodynamics code Athena++. To simulate dynamical space-times GR-Athena++ uses the Z4c evolution scheme of numerical relativity coupled to the moving puncture gauge. We demonstrate stable and accurate binary black hole merger evolutions via extensive convergence testing, cross-code validation, and verification against state-of-the-art effective-one-body waveforms. GR-Athena++ leverages the task-based parallelism paradigm of Athena++ to achieve excellent scalability. We measure strong scaling efficiencies above $95%$ for up to $sim 1.2times10^4$ CPUs and excellent weak scaling is shown up to $sim 10^5$ CPUs in a production binary black hole setup with adaptive mesh refinement. GR-Athena++ thus allows for the robust simulation of compact binary coalescences and offers a viable path towards numerical relativity at exascale.
The standard post-Newtonian approximation to gravitational waveforms, called T-approximants, from non-spinning black hole binaries are known not to be sufficiently accurate close to the last stable orbit of the system. A new approximation, called P-approximants, is believed to improve the accuracy of the waveforms rendering them applicable up to the last stable orbit. In this study we apply P-approximants to the case of a test-particle in equatorial orbit around a Kerr black hole parameterized by a spin parameter q that takes values between -1 and 1. In order to assess the performance of the two approximants we measure their effectualness (i.e. larger overlaps with the exact signal), and faithfulness (i.e. smaller biases while measuring the parameters of the signal) with the exact (numerical) waveforms. We find that in the case of prograde orbits, that is orbits whose angular momentum is in the same sense as the spin angular momentum of the black hole, T-approximant templates obtain an effectualness of ~ 0.99 for spins q < 0.75. For 0.75 < q < 0.95, the effectualness drops to about 0.82. The P-approximants achieve effectualness of > 0.99 for all spins up to q = 0.95. The bias in the estimation of parameters is much lower in the case of P-approximants than T-approximants. We find that P-approximants are both effectual and faithful and should be more effective than T-approximants as a detection template family when q>0. For q<0 both T- and P-approximants perform equally well so that either of them could be used as a detection template family.
We construct an approximate metric that represents the spacetime of spinning binary black holes (BBH) approaching merger. We build the metric as an analytical superposition of two Kerr metrics in harmonic coordinates, where we transform each black hole term with time-dependent boosts describing an inspiral trajectory. The velocities and trajectories of the boost are obtained by solving the post-Newtonian (PN) equations of motion at 3.5 PN order. We analyze the spacetime scalars of the new metric and we show that it is an accurate approximation of Einsteins field equations in vacuum for a BBH system in the inspiral regime. Furthermore, to prove the effectiveness of our approach, we test the metric in the context of a 3D general relativistic magneto-hydrodynamical (GRMHD) simulation of accreting mini-disks around the black holes. We compare our results with a previous well-tested spacetime construction based on the asymptotic matching method. We conclude that our new spacetime is well-suited for long-term GRMHD simulations of spinning binary black holes on their way to the merger.
We explore the benefits of adapted gauges to small mass ratio binary black hole evolutions in the moving puncture formulation. We find expressions that approximate the late time behavior of the lapse and shift, $(alpha_0,beta_0)$, and use them as initial values for their evolutions. We also use a position and black hole mass dependent damping term, $eta[vec{x}_1(t),vec{x}_2(t),m_1,m_2]$, in the shift evolution, rather than a constant or conformal-factor dependent choice. We have found that this substantially reduces noise generation at the start of the numerical integration and keeps the numerical grid stable around both black holes, allowing for more accuracy with lower resolutions. We test our choices for this gauge in detail in a case study of a binary with a 7:1 mass ratio, and then use 15:1 and 32:1 binaries for a convergence study. Finally, we apply our new gauge to a 64:1 binary and a 128:1 binary to well cover the comparable and small mass ratio regimes.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
