Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userAdapted gauge to small mass ratio binary black hole evolutions
88
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
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.
rate research
Read More
We perform a sequence of binary black hole simulations with increasingly small mass ratios, reaching to a 128:1 binary that displays 13 orbits before merger. Based on a detailed convergence study of the $q=m_1/m_2=1/15$ nonspinning case, we apply additional mesh refinements levels around the smaller hole horizon to reach successively the $q=1/32$, $q=1/64$, and $q=1/128$ cases. Roughly a linear computational resources scaling with $1/q$ is observed on 8-nodes simulations. We compute the remnant properties of the merger: final mass, spin, and recoil velocity, finding precise consistency between horizon and radiation measures. We also compute the gravitational waveforms: its peak frequency, amplitude, and luminosity. We compare those values with predictions of the corresponding phenomenological formulas, reproducing the particle limit within 2%, and we then use the new results to improve their fitting coefficients.
We explore different gauge choices in the moving puncture formulation in order to improve the accuracy of a linear momentum measure evaluated on the horizon of the remnant black hole produced by the merger of a binary. In particular, motivated by constant values studies, we design a gauge via a variable shift parameter $meta(vec{r}(t))$ such that it takes a low asymptotic (and at the orbiting punctures) value, while about the standard value of 2 at the final hole horizon. This choice then follows the remnant black hole as it moves due to its net recoil velocity. We find that this choice keeps the accuracy of the binary evolution and, once the asymptotic value of the parameter $meta$ is chosen about or below 1.0, it produces more accurate results for the recoil velocity than the corresponding evaluation of the radiated linear momentum at infinity, for typical numerical resolutions. We also find that the choice of the $partial_t$-gauge (at our working resolutions) is more accurate in this regard of computing recoil velocities than the $partial_0$-gauge. Detailed studies of an unequal mass $q=m_1/m_2=1/3$ nonspinning binary are provided and then verified for other mass ratios $(q=1/2,1/5)$ and spinning $(q=1)$ binary black hole mergers.
We perform a detailed parameter estimation study of binary black hole merger events reported in Zackay et al. and Venumadhav et al.. These are some of the faintest signals reported so far, and hence, relative to the loud events in the GWTC-1 catalog, the data should have lesser constraining power on their intrinsic parameters. Hence we examine the robustness of parameter inference to choices made in the analysis, as well as any potential systematics. We check the impact of different methods of estimating the noise power spectral density, different waveform models, and different priors for the compact object spins. For most of the events, the resulting differences in the inferred values of the parameters are much smaller than their statistical uncertainties. The estimation of the effective spin parameter $chi_{mathrm{eff}}$, i.e. the projection of the mass-weighted total spin along the angular momentum, can be sensitive to analysis choices for two of the sources with the largest effective spin magnitudes, GW151216 and GW170403. The primary differences arise from using a 3D isotropic spin prior: the tails of the posterior distributions should be interpreted with care and due consideration of the other data analysis choices.
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.
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]
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
