BZPEERPapers, people, evidence and research pathsDiscover
Research tools
Network
Opportunities
Account
Source author record
Ulrich Sperhake appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.
Catalog footprint
Research graph
Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.
BZPEER is loading the nearby papers, people, topics and institutions for this page.
Published work
GRChombo is an open-source code for performing Numerical Relativity time evolutions, built on top of the publicly available Chombo software for the solution of PDEs. Whilst GRChombo uses standard techniques in NR, it focusses on applications in theoretical physics where adaptability, both in terms of grid structure, and in terms of code modification, are key drivers.
We demonstrate the flexibility and utility of the Berger-Rigoutsos Adaptive Mesh Refinement (AMR) algorithm used in the open-source numerical relativity code GRChombo for generating gravitational waveforms from binary black-hole inspirals, and for studying other problems involving non-trivial matter configurations. We show that GRChombo can produce high quality binary black-hole waveforms through a code comparison with the established numerical relativity code Lean. We also discuss some of the technical challenges involved in making use of full AMR (as opposed to, e.g. moving box mesh refinement), including the numerical effects caused by using various refinement criteria when regridding. We suggest several "rules of thumb" for when to use different tagging criteria for simulating a variety of physical phenomena. We demonstrate the use of these different criteria through example evolutions of a scalar field theory. Finally, we also review the current status and general capabilities of GRChombo.
The mass of a Kerr black hole can be separated into irreducible and rotational components -the former is a lower limit to the energy that cannot be possibly extracted from the event horizon and is related to its area. Here we compute the irreducible masses of the stellar-mass black holes observed by gravitational-wave interferometers LIGO and Virgo. Using single-event data, we present a re-parametrization of the posterior distribution that explicitly highlights the irreducible and rotational contributions to the total energy. We exploit the area law to rank the black-hole mergers observed to date according to their irreversibility, thus providing a guide to selecting events for targeted tests of General Relativity. Using population fits, we compute the rate by which the total area of black-hole horizons increases due to the observable mergers.
Through numerical simulations of boson-star head-on collisions, we explore the quality of binary initial data obtained from the superposition of single-star spacetimes. Our results demonstrate that evolutions starting from a plain superposition of individual boosted boson-star spacetimes are vulnerable to significant unphysical artefacts. These difficulties can be overcome with a simple modification of the initial data suggested in [PRD 99 (2018) 044046] for collisions of oscillatons. While we specifically consider massive complex scalar field boson star models up to a 6th-order-polynomial potential, we argue that this vulnerability is universal and present in other kinds of exotic compact systems and hence needs to be addressed.
This paper provides an extended exploration of the inverse-chirp gravitational-wave signals from stellar collapse in massive scalar-tensor gravity reported in [Phys. Rev. Lett. {\bf 119}, 201103]. We systematically explore the parameter space that characterizes the progenitor stars, the equation of state and the scalar-tensor theory of the core collapse events. We identify a remarkably simple and straightforward classification scheme of the resulting collapse events. For any given set of parameters, the collapse leads to one of three end states, a weakly scalarized neutron star, a strongly scalarized neutron star or a black hole, possibly formed in multiple stages. The latter two end states can lead to strong gravitational-wave signals that may be detectable in present continuous-wave searches with ground-based detectors. We identify a very sharp boundary in the parameter space that separates events with strong gravitational-wave emission from those with negligible radiation.
We present in detail the Einstein equations in the Baumgarte-Shapiro-Shibata-Nakamura formulation for the case of $D$ dimensional spacetimes with $SO(D-d)$ isometry based on a method originally introduced in Ref.1. Regularized expressions are given for a numerical implementation of this method on a vertex centered grid including the origin of the quasi-radial coordinate that covers the extra dimensions with rotational symmetry. Axisymmetry, corresponding to the value $d=D-2$, represents a special case with fewer constraints on the vanishing of tensor components and is conveniently implemented in a variation of the general method. The robustness of the scheme is demonstrated for the case of a black-hole head-on collision in $D=7$ spacetime dimensions with $SO(4)$ symmetry.
Gravitational waves from coalescing binary black holes encode the evolution of their spins prior to merger. In the post-Newtonian regime and on the precession timescale, this evolution has one of three morphologies, with the spins either librating around one of two fixed points ("resonances") or circulating freely. In this work we perform full parameter estimation on resonant binaries with fixed masses and spin magnitudes, changing three parameters: a conserved "projected effective spin" $ξ$ and resonant family $ΔΦ=0,π$ (which uniquely label the source), the inclination $θ_{JN}$ of the binary's total angular momentum with respect to the line of sight (which determines the strength of precessional effects in the waveform), and the signal amplitude. We demonstrate that resonances can be distinguished for a wide range of binaries, except for highly symmetric configurations where precessional effects are suppressed. Motivated by new insight into double-spin evolution, we introduce new variables to characterize precessing black hole binaries which naturally reflects the timescale separation of the system and therefore better encode the dynamical information carried by gravitational waves.
We review recent progress in the application of numerical relativity techniques to astrophysics and high-energy physics. We focus on some developments that took place within the "Numerical Relativity and High Energy Physics" network, a Marie Curie IRSES action that we coordinated, namely: spin evolution in black hole binaries, high-energy black hole collisions, compact object solutions in scalar-tensor gravity, superradiant instabilities and hairy black hole solutions in Einstein's gravity coupled to fundamental fields, and the possibility to gain insight into these phenomena using analog gravity models.
We present numerical-relativity simulations of spherically symmetric core collapse and compact-object formation in scalar-tensor theories of gravity. The additional scalar degree of freedom introduces a propagating monopole gravitational-wave mode. Detection of monopole scalar waves with current and future gravitational-wave experiments may constitute smoking gun evidence for strong-field modifications of General Relativity. We collapse both polytropic and more realistic pre-supernova profiles using a high-resolution shock-capturing scheme and an approximate prescription for the nuclear equation of state. The most promising sources of scalar radiation are protoneutron stars collapsing to black holes. In case of a Galactic core collapse event forming a black hole, Advanced LIGO may be able to place independent constraints on the parameters of the theory at a level comparable to current Solar-System and binary-pulsar measurements. In the region of the parameter space admitting spontaneously scalarised stars, transition to configurations with prominent scalar hair before black-hole formation further enhances the emitted signal. Although a more realistic treatment of the microphysics is necessary to fully investigate the occurrence of spontaneous scalarisation of neutron star remnants, we speculate that formation of such objects could constrain the parameters of the theory beyond the current bounds obtained with Solar-System and binary-pulsar experiments.
We derive an effective potential for binary black-hole (BBH) spin precession at second post-Newtonian order. This effective potential allows us to solve the orbit-averaged spin-precession equations analytically for arbitrary mass ratios and spins. These solutions are quasiperiodic functions of time: after a fixed period the BBH spins return to their initial relative orientations and jointly precess about the total angular momentum by a fixed angle. Using these solutions, we classify BBH spin precession into three distinct morphologies between which BBHs can transition during their inspiral. We also derive a precession-averaged evolution equation for the total angular momentum that can be integrated on the radiation-reaction time and identify a new class of spin-orbit resonances that can tilt the direction of the total angular momentum during the inspiral. Our new results will help efforts to model and interpret gravitational waves from generic BBH mergers and predict the distributions of final spins and gravitational recoils.
The demand to obtain answers to highly complex problems within strong-field gravity has been met with significant progress in the numerical solution of Einstein's equations - along with some spectacular results - in various setups. We review techniques for solving Einstein's equations in generic spacetimes, focusing on fully nonlinear evolutions but also on how to benchmark those results with perturbative approaches. The results address problems in high-energy physics, holography, mathematical physics, fundamental physics, astrophysics and cosmology.
We continue our series of studies of high-energy collisions of black holes investigating unequal-mass, boosted head-on collisions in four dimensions. We show that the fraction of the center-of-mass energy radiated as gravitational waves becomes independent of mass ratio and approximately equal to $13\%$ at large energies. We support this conclusion with calculations using black hole perturbation theory and Smarr's zero-frequency limit approximation. These results lend strong support to the conjecture that the detailed structure of the colliding objects is irrelevant at high energies.
The dynamics of precessing binary black holes (BBHs) in the post-Newtonian regime has a strong timescale hierarchy: the orbital timescale is very short compared to the spin-precession timescale which, in turn, is much shorter than the radiation-reaction timescale on which the orbit is shrinking due to gravitational-wave emission. We exploit this timescale hierarchy to develop a multi-scale analysis of BBH dynamics elaborating on the analysis of Kesden et al. (2015). We solve the spin-precession equations analytically on the precession time and then implement a quasi-adiabatic approach to evolve these solutions on the longer radiation-reaction time. This procedure leads to an innovative "precession-averaged" post-Newtonian approach to studying precessing BBHs. We use our new solutions to classify BBH spin precession into three distinct morphologies, then investigate phase transitions between these morphologies as BBHs inspiral. These precession-averaged post-Newtonian inspirals can be efficiently calculated from arbitrarily large separations, thus making progress towards bridging the gap between astrophysics and numerical relativity.
Binary black holes on quasicircular orbits with spins aligned with their orbital angular momentum have been testbeds for analytic and numerical relativity for decades, not least because symmetry ensures that such configurations are equilibrium solutions to the spin-precession equations. In this work, we show that these solutions can be unstable when the spin of the higher-mass black hole is aligned with the orbital angular momentum and the spin of the lower-mass black hole is anti-aligned. Spins in these configurations are unstable to precession to large misalignment when the binary separation $r$ is between the values $r_{\rm ud\pm}= (\sqrt{χ_1} \pm \sqrt{q χ_2})^4 (1-q)^{-2} M$, where $M$ is the total mass, $q \equiv m_2/m_1$ is the mass ratio, and $χ_1$ ($χ_2$) is the dimensionless spin of the more (less) massive black hole. This instability exists for a wide range of spin magnitudes and mass ratios and can occur in the strong-field regime near merger. We describe the origin and nature of the instability using recently developed analytical techniques to characterize fully generic spin precession. This instability provides a channel to circumvent astrophysical spin alignment at large binary separations, allowing significant spin precession prior to merger affecting both gravitational-wave and electromagnetic signatures of stellar-mass and supermassive binary black holes.
Gravitational theories with multiple scalar fields coupled to the metric and each other --- a natural extension of the well studied single-scalar-tensor theories --- are interesting phenomenological frameworks to describe deviations from general relativity in the strong-field regime. In these theories, the $N$-tuple of scalar fields takes values in a coordinate patch of an $N$-dimensional Riemannian target-space manifold whose properties are poorly constrained by weak-field observations. Here we introduce for simplicity a non-trivial model with two scalar fields and a maximally symmetric target-space manifold. Within this model we present a preliminary investigation of spontaneous scalarization for relativistic, perfect fluid stellar models in spherical symmetry. We find that the scalarization threshold is determined by the eigenvalues of a symmetric scalar-matter coupling matrix, and that the properties of strongly scalarized stellar configurations additionally depend on the target-space curvature radius. In preparation for numerical relativity simulations, we also write down the $3+1$ decomposition of the field equations for generic tensor-multi-scalar theories.
One century after its formulation, Einstein's general relativity has made remarkable predictions and turned out to be compatible with all experimental tests. Most of these tests probe the theory in the weak-field regime, and there are theoretical and experimental reasons to believe that general relativity should be modified when gravitational fields are strong and spacetime curvature is large. The best astrophysical laboratories to probe strong-field gravity are black holes and neutron stars, whether isolated or in binary systems. We review the motivations to consider extensions of general relativity. We present a (necessarily incomplete) catalog of modified theories of gravity for which strong-field predictions have been computed and contrasted to Einstein's theory, and we summarize our current understanding of the structure and dynamics of compact objects in these theories. We discuss current bounds on modified gravity from binary pulsar and cosmological observations, and we highlight the potential of future gravitational wave measurements to inform us on the behavior of gravity in the strong-field regime.
The first fully non-linear numerical simulations of colliding charged black holes in D=4 Einstein-Maxwell theory were recently reported arXiv:1205.1063. These collisions were performed for black holes with equal charge-to-mass ratio, for which initial data can be found in closed analytic form. Here we generalize the study of collisions of charged black holes to the case of unequal charge-to-mass ratios. We focus on oppositely charged black holes, as to maximize acceleration-dependent effects. As |Q|/M increases from 0 to 0.99, we observe that the gravitational radiation emitted increases by a factor of ~ 2.7; the electromagnetic radiation emission becomes dominant for |Q|/M >~ 0.37 and at |Q|/M=0.99 is larger, by a factor of ~ 5.8, than its gravitational counterpart. We observe that these numerical results exhibit a precise and simple scaling with the charge. Furthermore, we show that the results from the numerical simulations are qualitatively captured by a simple analytic model that computes the electromagnetic dipolar radiation and the gravitational quadrupolar radiation of two non-relativistic interacting particles in Minkowski spacetime.
If binary black holes form following the successive core collapses of sufficiently massive binary stars, precessional dynamics may align their spins $\mathbf S_1$ and $\mathbf S_2$ and the orbital angular momentum $\mathbf L$ into a plane in which they jointly precess about the total angular momentum $\mathbf J$. These spin orientations are known as spin-orbit resonances since $\mathbf S_1$, $\mathbf S_2$, and $\mathbf L$ all precess at the same frequency to maintain their planar configuration. Two families of such spin-orbit resonances exist, alike in dignity but differentiated by whether the components of the two spins in the orbital plane are either aligned or antialigned. The fraction of binary black holes in each family is determined by the stellar evolution of their progenitors, so if gravitational-wave detectors could measure this fraction they could provide important insights into astrophysical formation scenarios for binary black holes. In this paper, we show that even under the conservative assumption that binary black holes are observed along the direction of $\mathbf J$ (where precession-induced modulations to the gravitational waveforms are minimized), the waveforms of many members of each resonant family can be distinguished from all members of the opposing family in events with signal-to-noise ratios $ρ\simeq 10$, typical of those expected for the first detections with Advanced LIGO/Virgo. We hope that our preliminary findings inspire a greater appreciation of the capability of gravitational-wave detectors to constrain stellar astrophysics and stimulate further studies of the distinguishability of spin-orbit resonant families in more expanded regions of binary black-hole parameter space.
The nonlinear behavior of higher dimensional black hole spacetimes is of interest in several contexts, ranging from an understanding of cosmic censorship to black hole production in high-energy collisions. However, nonlinear numerical evolutions of higher dimensional black hole spacetimes are tremendously complex, involving different diagnostic tools and "dimensional reduction methods". In this work we compare two different successful codes to evolve Einstein's equations in higher dimensions, and show that the results of such different procedures agree to numerical precision, when applied to the collision from rest of two equal-mass black holes. We calculate the total radiated energy to be E/M=9x10^{-4} in five dimensions and E/M=8.1x10^{-4} in six dimensions.
The nonlinear stability of Kerr-Newman black holes (KNBHs) is investigated by performing numerical simulations within the full Einstein-Maxwell theory. We take as initial data a KNBH with mass $M$, angular momentum to mass ratio $a$ and charge $Q$. Evolutions are performed to scan this parameter space within the intervals $0\le a/M\le 0.994$ and $0\le Q/M\le 0.996$, corresponding to an extremality parameter $a/a_{\rm max}$ ($a_{\rm max} \equiv \sqrt{M^2-Q^2}$) ranging from $0$ to $0.995$. These KNBHs are evolved, together with a small bar-mode perturbation, up to a time of order $120M$. Our results suggest that for small $Q/a$, the quadrupolar oscillation modes depend solely on $a/a_{\rm max}$, a universality also apparent in previous perturbative studies in the regime of small rotation. Using as a stability criterion the absence of significant relative variations in the horizon areal radius and BH spin, we find no evidence for any developing instability.
The evolution of black-hole binaries in vacuum spacetimes constitutes the two-body problem in general relativity. The solution of this problem in the framework of the Einstein field equations is a substantially more complex exercise than that of the dynamics of two point masses in Newtonian gravity, but it also presents us with a wealth of new exciting physics. Numerical methods are likely the only method to compute the dynamics of black-hole systems in the fully non-linear regime and have been pursued since the 1960s, culminating in dramatic breakthroughs in 2005. Here we review the methodology and the developments that finally gave us a solution of this fundamental problem of Einstein's theory and discuss the breakthrough's implication for the wide range of contemporary black-hole physics.
The Numerical-Relativity-Analytical-Relativity (NRAR) collaboration is a joint effort between members of the numerical relativity, analytical relativity and gravitational-wave data analysis communities. The goal of the NRAR collaboration is to produce numerical-relativity simulations of compact binaries and use them to develop accurate analytical templates for the LIGO/Virgo Collaboration to use in detecting gravitational-wave signals and extracting astrophysical information from them. We describe the results of the first stage of the NRAR project, which focused on producing an initial set of numerical waveforms from binary black holes with moderate mass ratios and spins, as well as one non-spinning binary configuration which has a mass ratio of 10. All of the numerical waveforms are analysed in a uniform and consistent manner, with numerical errors evaluated using an analysis code created by members of the NRAR collaboration. We compare previously-calibrated, non-precessing analytical waveforms, notably the effective-one-body (EOB) and phenomenological template families, to the newly-produced numerical waveforms. We find that when the binary's total mass is ~100-200 solar masses, current EOB and phenomenological models of spinning, non-precessing binary waveforms have overlaps above 99% (for advanced LIGO) with all of the non-precessing-binary numerical waveforms with mass ratios <= 4, when maximizing over binary parameters. This implies that the loss of event rate due to modelling error is below 3%. Moreover, the non-spinning EOB waveforms previously calibrated to five non-spinning waveforms with mass ratio smaller than 6 have overlaps above 99.7% with the numerical waveform with a mass ratio of 10, without even maximizing on the binary parameters.
We present an overview of recent developments in numerical relativity studies of higher dimensional spacetimes with a focus on time evolutions of black-hole systems. After a brief review of the numerical techniques employed for these studies, we summarize results grouped into the following three areas: (i) Numerical studies of fundamental properties of black holes, (ii) Applications of black-hole collisions to the modeling of Trans-Planckian scattering, (iii) Numerical studies of asymptotically anti-de Sitter spacetimes in the context of the gauge-gravity duality.
Scalar-tensor theories are a compelling alternative to general relativity and one of the most accepted extensions of Einstein's theory. Black holes in these theories have no hair, but could grow "wigs" supported by time-dependent boundary conditions or spatial gradients. Time-dependent or spatially varying fields lead in general to nontrivial black hole dynamics, with potentially interesting experimental consequences. We carry out a numerical investigation of the dynamics of single and binary black holes in the presence of scalar fields. In particular we study gravitational and scalar radiation from black-hole binaries in a constant scalar-field gradient, and we compare our numerical findings to analytical models. In the single black hole case we find that, after a short transient, the scalar field relaxes to static configurations, in agreement with perturbative calculations. Furthermore we predict analytically (and verify numerically) that accelerated black holes in a scalar-field gradient emit scalar radiation. For a quasicircular black-hole binary, our analytical and numerical calculations show that the dominant component of the scalar radiation is emitted at twice the binary's orbital frequency.
We study the influence of astrophysical formation scenarios on the precessional dynamics of spinning black-hole binaries by the time they enter the observational window of second- and third-generation gravitational-wave detectors, such as Advanced LIGO/Virgo, LIGO-India, KAGRA and the Einstein Telescope. Under the plausible assumption that tidal interactions are efficient at aligning the spins of few-solar mass black-hole progenitors with the orbital angular momentum, we find that black-hole spins should be expected to preferentially lie in a plane when they become detectable by gravitational-wave interferometers. This "resonant plane" is identified by the conditions ΔΦ=0° or ΔΦ=+/-180°, where ΔΦis the angle between the components of the black-hole spins in the plane orthogonal to the orbital angular momentum. If the angles ΔΦcan be accurately measured for a large sample of gravitational-wave detections, their distribution will constrain models of compact binary formation. In particular, it will tell us whether tidal interactions are efficient and whether a mechanism such as mass transfer, stellar winds, or supernovae can induce a mass-ratio reversal (so that the heavier black hole is produced by the initially lighter stellar progenitor). Therefore our model offers a concrete observational link between gravitational-wave measurements and astrophysics. We also hope that it will stimulate further studies of precessional dynamics, gravitational-wave template placement and parameter estimation for binaries locked in the resonant plane.
Light bosonic degrees of freedom have become a serious candidate for dark matter. The evolution of these fields around curved spacetimes is poorly understood but is expected to display interesting effects. In particular, the interaction of light bosonic fields with supermassive black holes, key players in most galaxies, could provide colourful examples of superradiance and nonlinear bosenova-like collapse. In turn, the observation of spinning black holes is expected to impose stringent bounds on the mass of putative massive bosonic fields in our universe. Our purpose here is to present a comprehensive study of the evolution of linearized massive scalar and vector fields in the vicinities of rotating black holes. For a certain boson field mass range, the field can become trapped in a potential barrier outside the horizon and transition to a bound state. Because there are a number of such quasi-bound states, the generic outcome is an amplitude modulated sinusoidal, or beating, signal. We believe that the appearance of such beatings has gone unnoticed in the past, and in fact mistaken for exponential growth. The amplitude modulation of the signal depends strongly on the relative excitation of the overtones, which in turn is strongly tied to the bound-state geography. For the first time we explore massive vector fields in generic BH background which are hard, if not impossible, to separate in the Kerr background. Our results show that spinning BHs are generically strongly unstable against massive vector fields.
Interferometric detectors will very soon give us an unprecedented view of the gravitational-wave sky, and in particular of the explosive and transient Universe. Now is the time to challenge our theoretical understanding of short-duration gravitational-wave signatures from cataclysmic events, their connection to more traditional electromagnetic and particle astrophysics, and the data analysis techniques that will make the observations a reality. This paper summarizes the state of the art, future science opportunities, and current challenges in understanding gravitational-wave transients.
We explore the impact of black hole spins on the dynamics of high-energy black hole collisions. We report results from numerical simulations with gamma-factors up to 2.49 and dimensionless spin parameter +0.85, +0.6, 0, -0.6, -0.85. We find that the scattering threshold becomes independent of spin at large center-of-mass energies, confirming previous conjectures that structure does not matter in ultrarelativistic collisions. It has further been argued that in this limit all of the kinetic energy of the system may be radiated by fine tuning the impact parameter to threshold. On the contrary, we find that only about 60% of the kinetic energy is radiated for gamma=2.49. By monitoring apparent horizons before and after scattering events we show that the "missing energy" is absorbed by the individual black holes in the encounter, and moreover the individual black-hole spins change significantly. We support this conclusion with perturbative calculations. An extrapolation of our results to the limit gamma->infinity suggests that about half of the center-of-mass energy of the system can be emitted in gravitational radiation, while the rest must be converted into rest-mass and spin energy.
We perform fully non-linear numerical simulations of charged-black-hole collisions, described by the Einstein-Maxwell equations, and contrast the results against analytic expectations. We focus on head-on collisions of non-spinning black holes, starting from rest and with the same charge to mass ratio, Q/M. The addition of charge to black holes introduces a new interesting channel of radiation and dynamics, most of which seem to be captured by Newtonian dynamics and flat-space intuition. The waveforms can be qualitatively described in terms of three stages; (i) an infall phase prior to the formation of a common apparent horizon; (ii) a nonlinear merger phase which corresponds to a peak in gravitational and electromagnetic energy; (iii) the ringdown marked by an oscillatory pattern with exponentially decaying amplitude and characteristic frequencies that are in good agreement with perturbative predictions. We observe that the amount of gravitational-wave energy generated throughout the collision decreases by about three orders of magnitude as the charge-to-mass ratio Q/M is increased from 0 to 0.98. We interpret this decrease as a consequence of the smaller accelerations present for larger values of the charge. In contrast, the ratio of energy carried by electromagnetic to gravitational radiation increases, reaching about 22% for the maximum Q/M ratio explored, which is in good agreement with analytic predictions.
Nonlinear dynamics in cosmological backgrounds has the potential to teach us immensely about our universe, and also to serve as prototype for nonlinear processes in generic curved spacetimes. Here we report on dynamical evolutions of black holes in asymptotically de Sitter spacetimes. We focus on the head-on collision of equal mass binaries and for the first time compare analytical and perturbative methods with full blown nonlinear simulations. Our results include an accurate determination of the merger/scatter transition (consequence of an expanding background) for small mass binaries and a test of the Cosmic Censorship conjecture, for large mass binaries. We observe that, even starting from small separations, black holes in large mass binaries eventually lose causal contact, in agreement with the conjecture.
Recent numerical relativity simulations have shown that the final black hole produced in a binary merger can recoil with a velocity as large as 5,000 km/s. Because of enhanced gravitational-wave emission in the so-called "hang-up" configurations, this maximum recoil occurs when the black-hole spins are partially aligned with the orbital angular momentum. We revisit our previous statistical analysis of post-Newtonian evolutions of black-hole binaries in the light of these new findings. We demonstrate that despite these new configurations with enhanced recoil velocities, spin alignment during the post-Newtonian stage of the inspiral will still significantly suppress (or enhance) kick magnitudes when the initial spin of the more massive black hole is more (or less) closely aligned with the orbital angular momentum than that of the smaller hole. We present a preliminary study of how this post-Newtonian spin alignment affects the ejection probabilities of supermassive black holes from their host galaxies with astrophysically motivated mass ratio and initial spin distributions. We find that spin alignment suppresses (enhances) ejection probabilities by ~ 40% (20%) for an observationally motivated mass-dependent galactic escape velocity, and by an even greater amount for a constant escape velocity of 1,000 km/s. Kick suppression is thus at least a factor two more efficient than enhancement.
Physics in curved spacetime describes a multitude of phenomena, ranging from astrophysics to high energy physics. The last few years have witnessed further progress on several fronts, including the accurate numerical evolution of the gravitational field equations, which now allows highly nonlinear phenomena to be tamed. Numerical relativity simulations, originally developed to understand strong field astrophysical processes, could prove extremely useful to understand high-energy physics processes like trans-Planckian scattering and gauge-gravity dualities. We present a concise and comprehensive overview of the state-of-the-art and important open problems in the field(s), along with guidelines for the next years. This writeup is a summary of the "NR/HEP Workshop" held in Madeira, Portugal from August 31st to September 3rd 2011.
The Numerical INJection Analysis (NINJA) project is a collaborative effort between members of the numerical relativity and gravitational wave data analysis communities. The purpose of NINJA is to study the sensitivity of existing gravitational-wave search and parameter-estimation algorithms using numerically generated waveforms, and to foster closer collaboration between the numerical relativity and data analysis communities. The first NINJA project used only a small number of injections of short numerical-relativity waveforms, which limited its ability to draw quantitative conclusions. The goal of the NINJA-2 project is to overcome these limitations with long post-Newtonian - numerical relativity hybrid waveforms, large numbers of injections, and the use of real detector data. We report on the submission requirements for the NINJA-2 project and the construction of the waveform catalog. Eight numerical relativity groups have contributed 63 hybrid waveforms consisting of a numerical portion modelling the late inspiral, merger, and ringdown stitched to a post-Newtonian portion modelling the early inspiral. We summarize the techniques used by each group in constructing their submissions. We also report on the procedures used to validate these submissions, including examination in the time and frequency domains and comparisons of waveforms from different groups against each other. These procedures have so far considered only the $(\ell,m)=(2,2)$ mode. Based on these studies we judge that the hybrid waveforms are suitable for NINJA-2 studies. We note some of the plans for these investigations.
We analyse an eleven-orbit inspiral of a non-spinning black-hole binary with mass ratio q=M1/M2=4. The numerically obtained gravitational waveforms are compared with post-Newtonian (PN) predictions including several sub-dominant multipoles up to multipolar indices (l=5,m=5). We find that (i) numerical and post-Newtonian predictions of the phase of the (2,2) mode accumulate a phase difference of about 0.35 rad at the PN cut off frequency 0.1 for the Taylor T1 approximant when numerical and PN waveforms are matched over a window in the early inspiral phase; (ii) in contrast to previous studies of equal-mass and specific spinning binaries, we find the Taylor T4 approximant to agree less well with numerical results, provided the latter are extrapolated to infinite extraction radius; (iii) extrapolation of gravitational waveforms to infinite extraction radius is particularly important for subdominant multipoles with l unequal m; (iv) 3PN terms in post-Newtonian multipole expansions significantly improve the agreement with numerical predictions for sub-dominant multipoles.
Numerical relativity has seen incredible progress in the last years, and is being applied with success to a variety of physical phenomena, from gravitational-wave research and relativistic astrophysics to cosmology and high-energy physics. Here we probe the limits of current numerical setups, by studying collisions of unequal mass, non-rotating black holes of mass-ratios up to 1:100 and making contact with a classical calculation in General Relativity: the infall of a point-like particle into a massive black hole. Our results agree well with the predictions coming from linearized calculations of the infall of point-like particles into non-rotating black holes. In particular, in the limit that one hole is much smaller than the other, and the infall starts from an infinite initial separation, we recover the point-particle limit. Thus, numerical relativity is able to bridge the gap between fully non-linear dynamics and linearized approximations, which may have important applications. Finally, we also comment on the "spurious" radiation content in the initial data and the linearized predictions.
We perform simulations of general relativistic rotating stellar core collapse and compute the gravitational waves (GWs) emitted in the core bounce phase of three representative models via multiple techniques. The simplest technique, the quadrupole formula (QF), estimates the GW content in the spacetime from the mass quadrupole tensor. It is strictly valid only in the weak-field and slow-motion approximation. For the first time, we apply GW extraction methods in core collapse that are fully curvature-based and valid for strongly radiating and highly relativistic sources. We employ three extraction methods computing (i) the Newman-Penrose (NP) scalar Psi_4, (ii) Regge-Wheeler-Zerilli-Moncrief (RWZM) master functions, and (iii) Cauchy-Characteristic Extraction (CCE) allowing for the extraction of GWs at future null infinity, where the spacetime is asymptotically flat and the GW content is unambiguously defined. The latter technique is the only one not suffering from residual gauge and finite-radius effects. All curvature-based methods suffer from strong non-linear drifts. We employ the fixed-frequency integration technique as a high-pass waveform filter. Using the CCE results as a benchmark, we find that finite-radius NP extraction yields results that agree nearly perfectly in phase, but differ in amplitude by ~1-7% at core bounce, depending on the model. RWZM waveforms, while in general agreeing in phase, contain spurious high-frequency noise of comparable amplitudes to those of the relatively weak GWs emitted in core collapse. We also find remarkably good agreement of the waveforms obtained from the QF with those obtained from CCE. They agree very well in phase but systematically underpredict peak amplitudes by ~5-11% which is comparable to the NP results and is within the uncertainties associated with core collapse physics. (abridged)
We calculate puncture initial data, corresponding to single and binary black holes with linear momenta, which solve the constraint equations of D dimensional vacuum gravity. The data are generated by a modification of the pseudo-spectral code presented in arXiv:gr-qc/0404056 and made available as the TwoPunctures thorn inside the Cactus computational toolkit. As examples, we exhibit convergence plots, the violation of the Hamiltonian constraint as well as the initial data for D=4,5,6,7. These initial data are the starting point to perform high energy collisions of black holes in D dimensions.
The puncture method for dealing with black holes in the numerical simulation of vacuum spacetimes is remarkably successful when combined with the BSSN formulation of the Einstein equations. We examine a generalized class of formulations modeled along the lines of the Laguna-Shoemaker system, including BSSN as a special case. The formulation is a two parameter generalization of the choice of variables used in standard BSSN evolutions. Numerical stability of the standard finite difference methods is proven for the formulation in the linear regime around flat space, a special case of which is the numerical stability of BSSN. Numerical evolutions are presented and compared with a standard BSSN implementation. We find that a significant portion of the parameter space leads to stable evolutions and that standard BSSN is located near the edge of the stability region. Non-standard parameter choices typically result in smoother behaviour of the evolution variables close to the puncture and thus hold promise for improved accuracy in, e.g., long-term BH binary inspirals, and for overcoming (numerical) stability problems still encountered in some types of black-hole simulations, e.g., in $D \ge 6$ dimensions.
The evolution of black holes in "confining boxes" is interesting for a number of reasons, particularly because it mimics the global structure of Anti-de Sitter geometries. These are non-globally hyperbolic space-times and the Cauchy problem may only be well defined if the initial data is supplemented by boundary conditions at the time-like conformal boundary. Here, we explore the active role that boundary conditions play in the evolution of a bulk black hole system, by imprisoning a black hole binary in a box with mirror-like boundary conditions. We are able to follow the post-merger dynamics for up to two reflections off the boundary of the gravitational radiation produced in the merger. We estimate that about 15% of the radiation energy is absorbed by the black hole per interaction, whereas transfer of angular momentum from the radiation to the black hole is only observed in the first interaction. We discuss the possible role of superradiant scattering for this result. Unlike the studies with outgoing boundary conditions, both the Newman-Penrose scalars Ψ_4 and Ψ_0 are non-trivial in our setup, and we show that the numerical data verifies the expected relations between them.
The inspiral of binary black holes is governed by gravitational radiation reaction at binary separations r < 1000 M, yet it is too computationally expensive to begin numerical-relativity simulations with initial separations r > 10 M. Fortunately, binary evolution between these separations is well described by post-Newtonian equations of motion. We examine how this post-Newtonian evolution affects the distribution of spin orientations at separations r ~ 10 M where numerical-relativity simulations typically begin. Although isotropic spin distributions at r ~ 1000 M remain isotropic at r ~ 10 M, distributions that are initially partially aligned with the orbital angular momentum can be significantly distorted during the post-Newtonian inspiral. Spin precession tends to align (anti-align) the binary black hole spins with each other if the spin of the more massive black hole is initially partially aligned (anti-aligned) with the orbital angular momentum, thus increasing (decreasing) the average final spin. Spin precession is stronger for comparable-mass binaries, and could produce significant spin alignment before merger for both supermassive and stellar-mass black hole binaries. We also point out that precession induces an intrinsic accuracy limitation (< 0.03 in the dimensionless spin magnitude, < 20 degrees in the direction) in predicting the final spin resulting from the merger of widely separated binaries.
In Ref. [1, 2] a formalism to deal with high-order perturbations of a general spherical background was developed. In this article, we apply it to the particular case of a perfect fluid background. We have expressed the perturbations of the energy-momentum tensor at any order in terms of the perturbed fluid's pressure, density and velocity. In general, these expressions are not linear and have sources depending on lower order perturbations. For the second-order case we make the explicit decomposition of these sources in tensor spherical harmonics. Then, a general procedure is given to evolve the perturbative equations of motions of the perfect fluid for any value of the harmonic label. Finally, with the problem of a spherical collapsing star in mind, we discuss the high-order perturbative matching conditions across a timelike surface, in particular the surface separating the perfect fluid interior from the exterior vacuum.
The numerical evolution of Einstein's field equations in a generic background has the potential to answer a variety of important questions in physics: from applications to the gauge-gravity duality, to modelling black hole production in TeV gravity scenarios, analysis of the stability of exact solutions and tests of Cosmic Censorship. In order to investigate these questions, we extend numerical relativity to more general space-times than those investigated hitherto, by developing a framework to study the numerical evolution of D dimensional vacuum space-times with an SO(D-2) isometry group for D\ge 5, or SO(D-3) for D\ge 6. Performing a dimensional reduction on a (D-4)-sphere, the D dimensional vacuum Einstein equations are rewritten as a 3+1 dimensional system with source terms, and presented in the Baumgarte, Shapiro, Shibata and Nakamura (BSSN) formulation. This allows the use of existing 3+1 dimensional numerical codes with small adaptations. Brill-Lindquist initial data are constructed in D dimensions and a procedure to match them to our 3+1 dimensional evolution equations is given. We have implemented our framework by adapting the LEAN code and perform a variety of simulations of non-spinning black hole space-times. Specifically, we present a modified moving puncture gauge which facilitates long term stable simulations in D=5. We further demonstrate the internal consistency of the code by studying convergence and comparing numerical versus analytic results in the case of geodesic slicing for D=5,6.
Black objects in higher dimensional space-times have a remarkably richer structure than their four dimensional counterparts. They appear in a variety of configurations (e.g. black holes, black branes, black rings, black Saturns), and display complex stability phase diagrams. They might also play a key role in high energy physics: for energies above the fundamental Planck scale, gravity is the dominant interaction which, together with the hoop-conjecture, implies that the trans-Planckian scattering of point particles should be well described by black hole scattering. Higher dimensional scenarios with a fundamental Planck scale of the order of TeV predict, therefore, black hole production at the LHC, as well as in future colliders with yet higher energies. In this setting, accurate predictions for the production cross-section and energy loss (through gravitational radiation) in the formation of black holes in parton-parton collisions is crucial for accurate phenomenological modelling in Monte Carlo event generators. In this paper, we use the formalism and numerical code reported in arXiv:1001.2302 to study the head-on collision of two black holes. For this purpose we provide a detailed treatment of gravitational wave extraction in generic D-dimensional space-times, which uses the Kodama-Ishibashi formalism. For the first time, we present the results of numerical simulations of the head-on collision in five space-time dimensions, together with the relevant physical quantities. We show that the total radiated energy, when two black holes collide from rest at infinity, is approximately (0.089\pm 0.006)% of the centre of mass energy, slightly larger than the 0.055% obtained in the four dimensional case, and that the ringdown signal at late time is in very good agreement with perturbative calculations.
Numerical-relativity simulations indicate that the black hole produced in a binary merger can recoil with a velocity up to v_max ~ 4,000 km/s with respect to the center of mass of the initial binary. This challenges the paradigm that most galaxies form through hierarchical mergers, yet retain supermassive black holes at their centers despite having escape velocities much less than v_max. Interaction with a circumbinary disk can align the binary black hole spins with their orbital angular momentum, reducing the recoil velocity of the final black hole produced in the subsequent merger. However, the effectiveness of this alignment depends on highly uncertain accretion flows near the binary black holes. In this Letter, we show that if the spin S_1 of the more massive binary black hole is even partially aligned with the orbital angular momentum L, relativistic spin precession on sub-parsec scales can align the binary black hole spins with each other. This alignment significantly reduces the recoil velocity even in the absence of gas. For example, if the angle between S_1 and L at large separations is 10 degrees while the second spin S_2 is isotropically distributed, the spin alignment discussed in this paper reduces the median recoil from 864 km/s to 273 km/s for maximally spinning black holes with a mass ratio of 9/11. This reduction will greatly increase the fraction of galaxies retaining their supermassive black holes.
Ultrarelativistic collisions of black holes are ideal gedanken experiments to study the nonlinearities of general relativity. In this paper we use semianalytical tools to better understand the nature of these collisions and the emitted gravitational radiation. We explain many features of the energy spectra extracted from numerical relativity simulations using two complementary semianalytical calculations. In the first calculation we estimate the radiation by a "zero-frequency limit" analysis of the collision of two point particles with finite impact parameter. In the second calculation we replace one of the black holes by a point particle plunging with arbitrary energy and impact parameter into a Schwarzschild black hole, and we explore the multipolar structure of the radiation paying particular attention to the near-critical regime. We also use a geodesic analogy to provide qualitative estimates of the dependence of the scattering threshold on the black hole spin and on the dimensionality of the spacetime.
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 ~3,000 km/s; here, we find configurations that exceed ~15,000 km/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.
We present the first numerical simulations of an initially non-spinning black-hole binary with mass ratio as large as 10:1 in full general relativity. The binary completes approximately 3 orbits prior to merger and radiates about 0.415% of the total energy and 12.48% of the initial angular momentum in the form of gravitational waves. The single black hole resulting from the merger acquires a kick of about 66.7 km/s relative to the original center of mass frame. The resulting gravitational waveforms are used to validate existing formulas for the recoil, final spin and radiated energy over a wider range of the symmetric mass ratio parameter eta=M1*M2/(M1+M2)^2 than previously possible. The contributions of l > 2 multipoles are found to visibly influence the gravitational wave signal obtained at fixed inclination angles.
Research on extracting science from binary-black-hole (BBH) simulations has often adopted a "scattering matrix" perspective: given the binary's initial parameters, what are the final hole's parameters and the emitted gravitational waveform? In contrast, we are using BBH simulations to explore the nonlinear dynamics of curved spacetime. Focusing on the head-on plunge, merger, and ringdown of a BBH with transverse, antiparallel spins, we explore numerically the momentum flow between the holes and the surrounding spacetime. We use the Landau-Lifshitz field-theory-in-flat-spacetime formulation of general relativity to define and compute the density of field energy and field momentum outside horizons and the energy and momentum contained within horizons, and we define the effective velocity of each apparent and event horizon as the ratio of its enclosed momentum to its enclosed mass-energy. We find surprisingly good agreement between the horizons' effective and coordinate velocities. To investigate the gauge dependence of our results, we compare pseudospectral and moving-puncture evolutions of physically similar initial data; although spectral and puncture simulations use different gauge conditions, we find remarkably good agreement for our results in these two cases. We also compare our simulations with the post-Newtonian trajectories and near-field energy-momentum. [Abstract abbreviated; full abstract also mentions additional results.]
The 2008 NRDA conference introduced the Numerical INJection Analysis project (NINJA), a new collaborative effort between the numerical relativity community and the data analysis community. NINJA focuses on modeling and searching for gravitational wave signatures from the coalescence of binary system of compact objects. We review the scope of this collaboration and the components of the first NINJA project, where numerical relativity groups shared waveforms and data analysis teams applied various techniques to detect them when embedded in colored Gaussian noise.
The Numerical INJection Analysis (NINJA) project is a collaborative effort between members of the numerical relativity and gravitational-wave data analysis communities. The purpose of NINJA is to study the sensitivity of existing gravitational-wave search algorithms using numerically generated waveforms and to foster closer collaboration between the numerical relativity and data analysis communities. We describe the results of the first NINJA analysis which focused on gravitational waveforms from binary black hole coalescence. Ten numerical relativity groups contributed numerical data which were used to generate a set of gravitational-wave signals. These signals were injected into a simulated data set, designed to mimic the response of the Initial LIGO and Virgo gravitational-wave detectors. Nine groups analysed this data using search and parameter-estimation pipelines. Matched filter algorithms, un-modelled-burst searches and Bayesian parameter-estimation and model-selection algorithms were applied to the data. We report the efficiency of these search methods in detecting the numerical waveforms and measuring their parameters. We describe preliminary comparisons between the different search methods and suggest improvements for future NINJA analyses.
We discuss possible origins for discrepancies observed in the radiated energies in head-on collisions of non-spinning binaries starting from Brill-Lindquist and superposed Kerr-Schild data. For this purpose, we discuss the impact of different choices of gauge parameters and a small initial boost of the black holes.
Currently the most popular method to evolve black-hole binaries is the ``moving puncture'' method. It has recently been shown that when puncture initial data for a Schwarzschild black hole are evolved using this method, the numerical slices quickly lose contact with the second asymptotically flat end, and end instead on a cylinder of finite Schwarzschild coordinate radius. These slices are stationary, meaning that their geometry does not evolve further. We will describe these results in the context of maximal slices, and present time-independent puncture-like data for the Schwarzschild spacetime.