Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userGRAthena++: puncture evolutions on vertex-centered oct-tree AMR
56
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
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.
rate research
Read More
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.
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.
Information about the last stages of a binary neutron star inspiral and the final merger can be extracted from quasi-equilibrium configurations and dynamical evolutions. In this article, we construct quasi-equilibrium configurations for different spins, eccentricities, mass ratios, compactnesses, and equations of state. For this purpose we employ the SGRID code, which allows us to construct such data in previously inaccessible regions of the parameter space. In particular, we consider spinning neutron stars in isolation and in binary systems; we incorporate new methods to produce highly eccentric and eccentricity reduced data; we present the possibility of computing data for significantly unequal-mass binaries; and we create equal-mass binaries with individual compactness up to 0.23. As a proof of principle, we explore the dynamical evolution of three new configurations. First, we simulate a $q=2.06$ mass ratio which is the highest mass ratio for a binary neutron star evolved in numerical relativity to date. We find that mass transfer from the companion star sets in a few revolutions before merger and a rest mass of $sim10^{-2}M_odot$ is transferred between the two stars. This configuration also ejects a large amount of material during merger, imparting a substantial kick to the remnant. Second, we simulate the first merger of a precessing binary neutron star. We present the dominant modes of the gravitational waves for the precessing simulation, where a clear imprint of the precession is visible in the (2,1) mode. Finally, we quantify the effect of an eccentricity reduction procedure on the gravitational waveform. The procedure improves the waveform quality and should be employed in future precision studies, but also other errors, notably truncation errors, need to be reduced in order for the improvement due to eccentricity reduction to be effective. [abridged]
A popular approach in numerical simulations of black hole binaries is to model black holes as punctures in the fabric of spacetime. The location and the properties of the black hole punctures are tracked with apparent horizons, namely outermost marginally outer trapped surfaces (MOTSs). As the holes approach each other, a common apparent horizon suddenly appears, engulfing the two black holes and signaling the merger. The evolution of common apparent horizons and their connection with gravitational wave emission have been studied in detail with the framework of dynamical horizons. We present a study of the dynamics of the MOTSs and their punctures in the interior of the final black hole. The study focuses on head-on mergers for various initial separations and mass ratios. We find that MOTSs intersect for most of the parameter space. We show that for those situations in which they do not, it is because of the singularity avoidance property of the moving puncture gauge condition used in the study. Although we are unable to carry out evolutions that last long enough to show the ultimate fate of the punctures, our results suggest that MOTSs always intersect and that at late times their overlap is only partial. As a consequence, the punctures inside the MOTSs, although close enough to each other to act effectively as a single puncture, do not merge.
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.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
