Full Numerical Relativity Research Articles

Reflecting boundary conditions in numerical relativity as a model for black hole echoes

Recently, there has been much interest in black hole echoes, based on the idea that there may be some mechanism (e.g. from quantum gravity) that waves/fields falling into a black hole could partially reflect off of an interface before reaching the horizon. There does not seem to be a good understanding of how to properly model a reflecting surface in numerical relativity, as the vast majority of the literature avoids the implementation of artificial boundaries, or applies transmitting boundary conditions. Here, we present a framework for reflecting a scalar field in a fully dynamical spherically symmetric spacetime, and implement it numerically. We study the evolution of a wave packet in this situation and its numerical convergence, including when the location of a reflecting boundary is very close to the horizon of a black hole. This opens the door to model exotic near-horizon physics within full numerical relativity.

Gravitational Waveforms for Compact Binaries from Second-Order Self-Force Theory.

We produce gravitational waveforms for nonspinning compact binaries undergoing a quasicircular inspiral. Our approach is based on a two-timescale expansion of the Einstein equations in second-order self-force theory, which allows first-principles waveform production in tens of milliseconds. Although the approach is designed for extreme mass ratios, our waveforms agree remarkably well with those from full numerical relativity, even for comparable-mass systems. Our results will be invaluable in accurately modeling extreme-mass-ratio inspirals for the LISA mission and intermediate-mass-ratio systems currently being observed by the LIGO-Virgo-KAGRA Collaboration.

The gravitational afterglow of boson stars

In this work we study the long-lived post-merger gravitational wave signature of a boson-star binary coalescence. We use full numerical relativity to simulate the post-merger and track the gravitational afterglow over an extended period of time. We implement recent innovations for the binary initial data, which significantly reduce spurious initial excitations of the scalar field profiles, as well as a measure for the angular momentum that allows us to track the total momentum of the spatial volume, including the curvature contribution. Crucially, we find the afterglow to last much longer than the spin-down timescale. This prolonged gravitational wave afterglow provides a characteristic signal that may distinguish it from other astrophysical sources.

Primordial black hole formation with full numerical relativity

We study the formation of black holes from subhorizon and superhorizon perturbations in a matter dominated universe with 3+1D numerical relativity simulations. We find that there are two primary mechanisms of formation depending on the initial perturbation's mass and geometry — via direct collapse of the initial overdensity and via post-collapse accretion of the ambient dark matter. In particular, for the latter case, the initial perturbation does not have to satisfy the hoop conjecture for a black hole to form. In both cases, the duration of the formation the process is around a Hubble time, and the initial mass of the black hole is M BH ∼ 10-2 H -1 M Pl 2. Post formation, we find that the PBH undergoes rapid mass growth beyond the self-similar limit M BH α H -1, at least initially. We argue that this implies that most of the final mass of the PBH is accreted from its ambient surroundings post formation.

General Relativistic Magnetohydrodynamics Mean-Field Dynamos

Large-scale, ordered magnetic fields in several astrophysical sources are supposed to be originated, and maintained against dissipation, by the combined amplifying action of rotation and small-scale turbulence. For instance, in the solar interior, the so-called α−Ω mean-field dynamo is known to be responsible for the observed 22-years magnetic cycle. Similar mechanisms could operate in more extreme environments, like proto neutron stars and accretion disks around black holes, for which the physical modelling needs to be translated from the regime of magnetohydrodynamics (MHD) and Newtonian gravity to that of a plasma in a general relativistic curved spacetime (GRMHD). Here we review the theory behind the mean field dynamo in GRMHD, the strategies for the implementation of the relevant equations in numerical conservative schemes, and we show the most important applications to the mentioned astrophysical compact objects obtained by our group in Florence. We also present novel results, such as three-dimensional GRMHD simulations of accretion disks with dynamo and the application of our dynamo model to a super massive neutron star, remnant of a binary neutron star merger as obtained from full numerical relativity simulations.

Up-down instability of binary black holes in numerical relativity

Binary black holes with spins that are aligned with the orbital angular momentum do not precess. However, post-Newtonian calculations predict that "up-down" binaries, in which the spin of the heavier (lighter) black hole is aligned (antialigned) with the orbital angular momentum, are unstable when the spins are slightly perturbed from perfect alignment. This instability provides a possible mechanism for the formation of precessing binaries in environments where sources are preferentially formed with (anti) aligned spins. In this paper, we present the first full numerical relativity simulations capturing this instability. These simulations span $\sim 100$ orbits and $\sim 3$-$5$ precession cycles before merger, making them some of the longest numerical relativity simulations to date. Initialized with a small perturbation of $1^{\circ}$-$10^{\circ}$, the instability causes a dramatic growth of the spin misalignments, which can reach $\sim 90^{\circ}$ near merger. We show that this leaves a strong imprint on the subdominant modes of the gravitational wave signal, which can potentially be used to distinguish up-down binaries from other sources. Finally, we show that post-Newtonian and effective-one-body approximants are able to reproduce the unstable dynamics of up-down binaries extracted from numerical relativity.

Oscillon preheating in full general relativity

Oscillons are dense nonperturbative objects that may be copiously produced in the preheating period after inflation. Oscillon preheatings are usually simulated with nonlinear matter interactions but in a rigid FLRW background, without taking into account the gravitational backreactions. We investigate the oscillon preheating scenario in full general relativity with a full numerical relativity scheme, and compare the general-relativistic simulations with the traditional ones done in an FLRW background. We find that in certain parameter regions the general-relativistic corrections are significant, producing more and denser oscillons, and can be strong enough to collapse the oscillons to black holes.

Surrogate model for gravitational wave signals from comparable and large-mass-ratio black hole binaries

Gravitational wave signals from compact astrophysical sources such as those observed by LIGO and Virgo require a high-accuracy, theory-based waveform model for the analysis of the recorded signal. Current inspiral-merger-ringdown models are calibrated only up to moderate mass ratios, thereby limiting their applicability to signals from high-mass ratio binary systems. We present EMRISur1dq1e4, a reduced-order surrogate model for gravitational waveforms of 13,500M in duration and including several harmonic modes for non-spinning black hole binary systems with mass-ratios varying from 3 to 10,000 thus vastly expanding the parameter range beyond the current models. This surrogate model is trained on waveform data generated by point-particle black hole perturbation theory (ppBHPT) both for large mass-ratio and comparable mass-ratio binaries. We observe that the gravitational waveforms generated through a simple application of ppBHPT to the comparable mass-ratio cases agree remarkably (and surprisingly) well with those from full numerical relativity after a rescaling of the ppBHPT's total mass parameter. This observation and the EMRISur1dq1e4 surrogate model will enable data analysis studies in the high-mass ratio regime, including potential intermediate mass-ratio signals from LIGO/Virgo and extreme-mass ratio events of interest to the future space-based observatory LISA.

Black hole formation in relativistic Oscillaton collisions

We investigate the physics of black hole formation from the head-on collisions of boosted equal mass Oscillatons (OS) in full numerical relativity, for both the cases where the OS have equal phases or are maximally off-phase (anti-phase). While unboosted OS collisions will form a BH as long as their initial compactness \U0001d49e≡ GM/R is above a numerically determined critical value \U0001d49e>0.035, we find that imparting a small initial boost counter-intuitively prevents the formation of black holes even if \U0001d49e> 0.035. If the boost is further increased, at very high boosts γ>1/12\U0001d49e, BH formation occurs as predicted by the hoop conjecture. These two limits combine to form a “stability band” where collisions result in either the OS “passing through” (equal phase) or “bouncing back” (anti-phase), with a critical point occurring around \U0001d49e≈ 0.07. We argue that the existence of this stability band can be explained by the competition between the free fall and the interaction timescales of the collision.

Induced spins from scattering experiments of initially nonspinning black holes

When two relativistically boosted, nonspinning black holes pass by one another on a scattering trajectory, we might expect the tidal interaction to spin up each black hole. We present the first exploration of this effect, appearing at fourth post-Newtonian order, with full numerical relativity calculations. The basic setup for the calculations involves two free parameters: the initial boost of each black hole and the initial angle between the velocity vectors and a line connecting the centers of the black holes, with zero angle corresponding to a head-on trajectory. To minimize gauge effects, we measure final spins only if the black holes reach a final separation of at least $20M$. Fixing the initial boost, we find that as the initial angle decreases toward the scattering/nonscattering limit, the spin-up grows nonlinearly. In addition, as initial boosts are increased from $0.42c$ to $0.78c$, the largest observed final dimensionless spin on each black hole increases nonlinearly from $0.02$ to $0.20$. Based on these results, we conclude that much higher spin-ups may be possible with larger boosts, although achieving this will require improved numerical techniques.

Kicking gravitational wave detectors with recoiling black holes

Binary black holes emit gravitational radiation with net linear momentum leading to a retreat of the final remnant black hole that can reach up to $\sim5,000$ km/s. Full numerical relativity simulations are the only tool to accurately compute these recoils since they are largely produced when the black hole horizons are about to merge and they are strongly dependent on their spin orientations at that moment. We present eight new numerical simulations of BBH in the hangup-kick configuration family, leading to the maximum recoil. Black holes are equal mass and near maximally spinning ($|\vec{S}_{1,2}|/m_{1,2}^2=0.97$). Depending on their phase at merger, this family leads to $\sim\pm4,700$ km/s and all intermediate values of the recoil along the orbital angular momentum of the binary system. We introduce a new invariant method to evaluate the recoil dependence on the merger phase via the waveform peak amplitude used as a reference phase angle and compare it with previous definitions. We also compute the mismatch between these hangup-kick waveforms to infer their observable differentiability by gravitational wave detectors, such as advanced LIGO, finding currently reachable signal-to-noise ratios, hence allowing for the identification of highly recoiling black holes having otherwise essentially the same binary parameters.

Preheating in full general relativity

We investigate the importance of local gravity during preheating, the non-linear dynamics that may be responsible for starting the process of reheating the universe after inflation. We introduce three numerical methods that study a simple preheating scenario while relaxing gravitational assumptions, culminating in studying the process in full numerical relativity. We confirm that perturbation theory is no longer valid when one considers modes whose wavelengths are comparable to the size of the horizon at the end of inflation; however, this breakdown does not necessarily lead to a breakdown of the preheating process in non-linear gravity. For the specific model we test we find no evidence for the creation of primordial black holes from the instabilities in this model. Finally, we remark on the opportunity for future numerical study of non-linear gravitational dynamics in the early universe.

Gravitational Waves from Binary Black Hole Mergers inside Stars.

We present results from a controlled numerical experiment investigating the effect of stellar density gas on the coalescence of binary black holes (BBHs) and the resulting gravitational waves (GWs). This investigation is motivated by the proposed stellar core fragmentation scenario for BBH formation and the associated possibility of an electromagnetic counterpart to a BBH GW event. We employ full numerical relativity coupled with general-relativistic hydrodynamics and set up a 30+30 M_{⊙} BBH (motivated by GW150914) inside gas with realistic stellar densities. Our results show that at densities ρ≳10^{6}-10^{7} g cm^{-3} dynamical friction between the BHs and gas changes the coalescence dynamics and the GW signal in an unmistakable way. We show that for GW150914, LIGO observations appear to rule out BBH coalescence inside stellar gas of ρ≳10^{7} g cm^{-3}. Typical densities in the collapsing cores of massive stars are in excess of this density. This excludes the fragmentation scenario for the formation of GW150914.

Modeling the source of GW150914 with targeted numerical-relativity simulations

In fall of 2015, the two LIGO detectors measured the gravitational wave signal GW150914, which originated from a pair of merging black holes (Abbott et al Virgo, LIGO Scientific 2016 Phys. Rev. Lett. 116 061102). In the final 0.2 s (about 8 gravitational-wave cycles) before the amplitude reached its maximum, the observed signal swept up in amplitude and frequency, from 35 Hz to 150 Hz. The theoretical gravitational-wave signal for merging black holes, as predicted by general relativity, can be computed only by full numerical relativity, because analytic approximations fail near the time of merger. Moreover, the nearly-equal masses, moderate spins, and small number of orbits of GW150914 are especially straightforward and efficient to simulate with modern numerical-relativity codes. In this paper, we report the modeling of GW150914 with numerical-relativity simulations, using black-hole masses and spins consistent with those inferred from LIGO’s measurement (Abbott et al LIGO Scientific Collaboration, Virgo Collaboration 2016 Phys. Rev. Lett. 116 241102). In particular, we employ two independent numerical-relativity codes that use completely different analytical and numerical methods to model the same merging black holes and to compute the emitted gravitational waveform; we find excellent agreement between the waveforms produced by the two independent codes. These results demonstrate the validity, impact, and potential of current and future studies using rapid-response, targeted numerical-relativity simulations for better understanding gravitational-wave observations.

Quasiequilibrium sequences of binary neutron stars undergoing dynamical scalarization

We calculate quasiequilibrium sequences of equal-mass, irrotational binary neutron stars (BNSs) in a scalar-tensor (ST) theory of gravity that admits dynamical scalarization. We model neutron stars with realistic equations of state (notably through piecewise polytropic equations of state). Using these quasiequilibrium sequences we compute the binary's scalar charge and binding energy versus orbital angular frequency. We find that the absolute value of the binding energy is smaller than in general relativity (GR), differing at most by ~14% at high frequencies for the cases considered. We use the newly computed binding energy and the balance equation to estimate the number of gravitational-wave (GW) cycles during the adiabatic, quasicircular inspiral stage up to the end of the sequence, which is the last stable orbit or the mass-shedding point, depending on which comes first. We find that, depending on the ST parameters, the number of GW cycles can be substantially smaller than in GR. In particular, we obtain that when dynamical scalarization sets in around a GW frequency of ~130 Hz, the sole inclusion of the ST binding energy causes a reduction of GW cycles from ~120 Hz up to the end of the sequence (~1200 Hz) of ~11% with respect to the GR case. We estimate that when the ST energy flux is also included the reduction in GW cycles becomes of ~24%. Quite interestingly, dynamical scalarization can produce a difference in the number of GW cycles with respect to the GR point-particle case that is much larger than the effect due to tidal interactions, which is on the order of only a few GW cycles. These results further clarify and confirm recent studies that have evolved BNSs either in full numerical relativity or in post-Newtonian theory, and point out the importance of developing accurate ST-theory waveforms for systems composed of strongly self-gravitating objects, such as BNSs.

Superkicks in ultrarelativistic encounters of spinning black holes

We study ultrarelativistic encounters of two spinning, equal-mass black holes through simulations in full numerical relativity. Two initial data sequences are studied in detail: one that leads to scattering and one that leads to a grazing collision and merger. In all cases, the initial black hole spins lie in the orbital plane, a configuration that leads to the so-called superkicks. In astrophysical, quasicircular inspirals, such kicks can be as large as $\ensuremath{\sim}3000\text{ }\text{ }\mathrm{km}/\mathrm{s}$; here, we find configurations that exceed $\ensuremath{\sim}15\text{ }000\text{ }\text{ }\mathrm{km}/\mathrm{s}$. We find that the maximum recoil is to a good approximation proportional to the total amount of energy radiated in gravitational waves, but largely independent of whether a merger occurs or not. This shows that the mechanism predominantly responsible for the superkick is not related to merger dynamics. Rather, a consistent explanation is that the ``bobbing'' motion of the orbit causes an asymmetric beaming of the radiation produced by the in-plane orbital motion of the binary, and the net asymmetry is balanced by a recoil. We use our results to formulate some conjectures on the ultimate kick achievable in any black hole encounter.

Zoom-Whirl Orbits in Black Hole Binaries

Zoom-whirl behavior has the reputation of being a rare phenomenon. The concern has been that gravitational radiation would drain angular momentum so rapidly that generic orbits would circularize before zoom-whirl behavior could play out, and only rare highly tuned orbits would retain their imprint. Using full numerical relativity, we catch zoom-whirl behavior despite dissipation. The larger the mass ratio, the longer the pair can spend in orbit before merging and therefore the more zooms and whirls seen. Larger spins also enhance zoom whirliness. An important implication is that these eccentric orbits can merge during a whirl phase, before enough angular momentum has been lost to truly circularize the orbit. Waveforms will be modulated by the harmonics of zoom-whirls, showing quiet phases during zooms and louder glitches during whirls.

Numerical relativistic model of a massive particle in orbit near a Schwarzschild black hole

We present a method for computing the evolution of a spacetime containing a massive particle and a black hole. The essential idea is that the gravitational field is evolved using full numerical relativity, with the particle generating a nonzero source term in the Einstein equations. The matter fields are not evolved by hydrodynamic equations. Instead the particle is treated as a quasirigid body whose center follows a geodesic. The necessary theoretical framework is developed and then implemented in a computer code that uses the null-cone, or characteristic, formulation of numerical relativity. The performance of the code is illustrated in test runs, including a complete orbit (near $r=9M)$ of a Schwarzschild black hole.

Modeling gravitational radiation from coalescing binary black holes

With the goal of bringing theory, particularly numerical relativity, to bear on an astrophysical problem of critical interest to gravitational wave observers we introduce a model for coalescence radiation from binary black hole systems. We build our model using the Lazarus approach, a technique that bridges far and close limit approaches with full numerical relativity to solve Einstein equations applied in the truly nonlinear dynamical regime. We specifically study the post-orbital radiation from a system of equal-mass non-spinning black holes, deriving waveforms which indicate strongly circularly polarized radiation of roughly 3% of the system's total energy and 12% of its total angular momentum in just a few cycles. To support this result we first establish the reliability of the late-time part of our model, including the numerical relativity and close-limit components, with a thorough study of waveforms from a sequence of black hole configurations that varies from previously treated head-on collisions to a representative target for ``ISCO'' data corresponding to the end of the inspiral period. We then complete our model with a simple treatment for the early part of the spacetime based on a standard family of initial data for binary black holes in circular orbit. A detailed analysis shows strong robustness in the results as the initial separation of the black holes is increased from $5.0$ to $7.8M$ supporting our waveforms as a suitable basic description of the astrophysical radiation from this system. Finally, a simple fitting of the plunge waveforms is introduced as a first attempt to facilitate the task of analyzing data from gravitational wave detectors.

Gravitational waves from black hole collisions via an eclectic approach

We present the first results in a new program intended to make the best use of all available technologies to provide an effective understanding of waves from inspiralling black hole binaries in time to assist imminent observations. In particular, we address the problem of combining the close-limit approximation describing ringing black holes and full numerical relativity, required for essentially nonlinear interactions. We demonstrate the effectiveness of our approach using general methods for a model problem, the head-on collision of black holes. Our method allows a more direct physical understanding of these collisions indicating clearly when non-linear methods are important. The success of this method supports our expectation that this unified approach will be able to provide relevant results for black hole binaries in time to assist gravitational wave observations.