Search

We study the dynamical formation of scalar monopole and dipole hair in scalar Gauss-Bonnet theory and dynamical Chern-Simons theory. We prove that the spherically-symmetric mode of the dipole hair is completely determined by the product of the mass of the spacetime and the value of the monopole hair. We then show that the dynamics of the $\ell=1$ mode of the dipole hair is intimately tied to the appearance of the event horizon during axisymmetric collapse, which results in the radiation of certain modes that could have been divergent in the future of the collapse. We confirm these analytical predictions by simulating the gravitational collapse of a rapidly rotating neutron star in the decoupling limit, both in scalar Gauss-Bonnet and dynamical Chern-Simons theory. Our results, combined with those of Ref.~\cite{R:2022cwe}, provide a clear physical picture of the dynamics of scalar monopole and dipole radiation in axisymmetric and spherical gravitational collapse in these theories., v2-matches published version in PRD

K. G. Arun et al., The Laser Interferometer Space Antenna (LISA) has the potential to reveal wonders about the fundamental theory of nature at play in the extreme gravity regime, where the gravitational interaction is both strong and dynamical. In this white paper, the Fundamental Physics Working Group of the LISA Consortium summarizes the current topics in fundamental physics where LISA observations of gravitational waves can be expected to provide key input. We provide the briefest of reviews to then delineate avenues for future research directions and to discuss connections between this working group, other working groups and the consortium work package teams. These connections must be developed for LISA to live up to its science potential in these areas., E. Berti is supported by NSF Grants No. PHY-1912550 and AST-2006538, NASA ATP Grants No. 17-ATP17-0225 and 19-ATP19-0051, NSF-XSEDE Grant No. PHY-090003, and NSF Grant PHY-20043. D. Blas is supported by a ‘Ayuda Beatriz Galindo Senior’ from the Spanish ‘Ministerio de Universidades’, grant BG20/00228. IFAE is partially funded by the CERCA program of the Generalitat de Catalunya. The research leading of to these results has received funding from the Spanish Ministry of Science and Innovation (PID2020-115845GB-I00/AEI/10.13039/501100011033). K. Clough is supported by funding from the European Research Council (ERC) under the European Union’s Horizon 2020 research and innovation programme (grant agreement No 693024). A. Cárdenas-Avendaño acknowledges funding from the Fundación Universitaria Konrad Lorenz (Project 5INV1) and from Will and Kacie Snellings. M. Crisostomi and E. Barausse are supported by the European Union’s H2020 ERC Consolidator Grant “GRavity from Astrophysical to Microscopic Scales” (Grant No. GRAMS-815673). P. Fleury received the support of a fellowship from “la Caixa” Foundation (ID 100010434). The fellowship code is LCF/BQ/PI19/11690018. C. Herdeiro thanks the support of the Center for Research and Development in Mathematics and Applications (CIDMA) through the Portuguese Foundation for Science and Technology (FCT - Fundação para a Ciência e a Tecnologia), references UIDB/04106/2020 and UIDP/04106/2020, the projects PTDC/FIS-OUT/28407/2017, CERN/FIS-PAR/0027/2019, PTDC/FIS-AST/3041/2020 and the European Union’s Horizon 2020 research and innovation (RISE) programme H2020-MSCA-RISE-2017 Grant No. FunFiCO-777740. P. Pani and E. Maggio acknowledge financial support provided under the European Union’s H2020 ERC, Starting Grant agreement no. DarkGRA–757480, and under the MIUR PRIN and FARE programmes (GW-NEXT, CUP: B84I20000100001), and support from the Amaldi Research Center funded by the MIUR program “Dipartimento di Eccellenza” (CUP: B81I18001170001). N.Frusciante was supported by Fundação para a Ciência e a Tecnologia (FCT) through the research grants UIDB/04434/2020, UIDP/04434/2020, PTDC/FIS-OUT/29048/2017, CERN/FIS-PAR/0037/2019, the FCT project “CosmoTests—Cosmological tests of gravity theories beyond General Relativity” with ref. number CEECIND/00017/2018 and the FCT project “BEYLA –BEYond LAmbda” with ref. number PTDC/FIS-AST/0054/2021. L.Lombriser was supported by a Swiss National Science Foundation Professorship grant (No. 170547). S.N. acknowledges support from the Alexander von Humboldt Foundation. D.N. acknowledges support from the NSF Grant No. PHY-2011784. R.B. acknowledges financial support from FCT – Fundação para a Ciência e a Tecnologia, I.P., under the Scientific Employment Stimulus - Individual Call - 2020.00470.CEECIND. V. Paschalidis acknowledges support from NSF Grant PHY-1912619 and NASA Grant 80NSSC20K1542 to the University of Arizona. B.S.S. is supported by NSF grants No. AST-2006384 and PHY-2012083. C.F.S. is supported by contracts ESP2017-90084-P and PID2019-106515GB-I00/AEI/10.13039/501100011033 (Spanish Ministry of Science and Innovation) and 2017-SGR-1469 (AGAUR, Generalitat de Catalunya). T. P. S. acknowledges partial support from the STFC Consolidated Grant No. ST/P000703/1. M. Ruiz acknowledges support from NASA Grant 80NSSC17K0070 to the University of Illinois at Urbana-Champaign. I.D. Saltas is supported by the Czech Science Foundation GAČR, Grant No. 21-16583M. N. Stergioulas is supported by the ESA Prodex grant PEA:4000132310 “LISA Stochastic Signals Analysis Pipeline”. F.V. acknowledges partial support from CNES. K.Y. acknowledges support from NSF Grant PHY-1806776, NASA Grant 80NSSC20K0523, a Sloan Foundation Research Fellowship and the Owens Family Foundation. K.Y. would like to also acknowledge support by the COST Action GWverse CA16104 and JSPS KAKENHI Grants No. JP17H06358. N. Yunes acknowledges support from NASA Grants No. NNX16AB98G, 80NSSC17M0041 and 80NSSC18K1352, NSF Award No. 1759615, and the Simons Foundation through MPS Award Number 896696. D.D. acknowledge financial support via an Emmy Noether Research Group funded by the German Research Foundation (DFG) under grant no. DO 1771/1-1.

Particular couplings between a scalar field and the Gauss-Bonnet invariant lead to spontaneous scalarization of black holes. Here we continue our work on simulating this phenomenon in the context of binary black hole systems. We consider a negative coupling for which the black-hole spin plays a major role in the scalarization process. We find two main phenomena: (i) dynamical descalarization, in which initially scalarized black holes form an unscalarized remnant, and (ii) dynamical scalarization, whereby the late merger of initially unscalarized black holes can cause scalar hair to grow. An important consequence of the latter case is that modifications to the gravitational waveform due to the scalar field may only occur post-merger, as its presence is hidden during the entirety of the inspiral. However, with a sufficiently strong coupling, we find that scalarization can occur before the remnant has even formed. We close with a discussion of observational implications for gravitational-wave tests of general relativity., 16 pages, 12 figures

In this work, we estimate how much bulk viscosity driven by Urca processes is likely to affect the gravitational wave signal of a neutron star coalescence. In the late inspiral, we show that bulk viscosity affects the binding energy at fourth post-Newtonian (PN) order. Even though this effect is enhanced by the square of the gravitational compactness, the coefficient of bulk viscosity is likely too small to lead to observable effects in the waveform during the late inspiral, when only considering the orbital motion itself. In the post-merger, however, the characteristic time-scales and spatial scales are different, potentially leading to the opposite conclusion. We post-process data from a state-of-the-art equal-mass binary neutron star merger simulation to estimate the effects of bulk viscosity (which was not included in the simulation itself). In that scenario, we find that bulk viscosity can reach high values in regions of the merger. We compute several estimates of how much it might directly affect the global dynamics of the considered merger scenario, and find that it could become significant. Even larger effects could arise in different merger scenarios or in simulations that include non-linear effects. This assessment is reinforced by a quantitative comparison with relativistic heavy-ion collisions where such effects have been explored extensively., Comment: 12 pages, 6 figures

Black holes in certain modified gravity theories that contain a scalar field coupled to curvature invariants are known to possess (monopole) scalar hair while non-black-hole spacetimes (like neutron stars) do not. Therefore, as a neutron star collapses to a black hole, scalar hair must grow until it settles to the stationary black hole solution with (monopole) hair. In this paper, we study this process in detail and show that the growth of scalar hair is tied to the appearance and growth of the event horizon (before an apparent horizon forms), which forces scalar modes that would otherwise (in the future) become divergent to be radiated away. We prove this result rigorously in general first for a large class of modified theories, and then we exemplify the results by studying the temporal evolution of the scalar field in scalar Gauss-Bonnet gravity in two backgrounds: (i) a collapsing Oppenheimer-Snyder background, and (ii) a collapsing neutron star background. In case (i), we find an exact scalar field solution analytically, while in case (ii) we solve for the temporal evolution of the scalar field numerically, with both cases supporting the conclusion presented above. Our results suggest that the emission of a burst of scalar field radiation is a necessary condition for black hole formation in a large class of modified theories of gravity., Comment: v2 : 23 pages, minor corrections in text, matches published version

Scalar fields coupled to the Gauss-Bonnet invariant can undergo a tachyonic instability, leading to spontaneous scalarization of black holes. Studies of this effect have so far been restricted to single black hole spacetimes. We present the first results on dynamical scalarization in head-on collisions and quasicircular inspirals of black hole binaries with numerical relativity simulations. We show that black hole binaries can either form a scalarized remnant or dynamically descalarize by shedding off its initial scalar hair. The observational implications of these findings are discussed., 9 pages, 6 figures. (v2) new references added; (v3) matches published version. For movies of the simulations see https://bhscalarization.bitbucket.io

The ability to test general relativity in extreme gravity regimes using gravitational wave observations from current ground-based or future space-based detectors motivates the mathematical study of the symmetries of black holes in modified theories of gravity. In this paper we focus on spinning black hole solutions in two quadratic gravity theories: dynamical Chern-Simons and scalar Gauss-Bonnet gravity. We compute the principal null directions, Weyl scalars, and complex null tetrad in the small-coupling, slow rotation approximation for both theories, confirming that both spacetimes are Petrov type I. Additionally, we solve the Killing equation through rank 6 in dynamical Chern-Simons gravity and rank 2 in scalar Gauss-Bonnet gravity, showing that there is no nontrivial Killing tensor through those ranks for each theory. We therefore conjecture that the still-unknown, exact, quadratic-gravity, black-hole solutions do not possess a fourth constant of motion., Comment: 23 pages, 2 figures

Gravitational waves emitted by black hole binary inspiral and mergers enable unprecedented strong-field tests of gravity, requiring accurate theoretical modelling of the expected signals in extensions of General Relativity. In this paper we model the gravitational wave emission of inspiraling binaries in scalar Gauss-Bonnet gravity theories. Going beyond the weak-coupling approximation, we derive the gravitational waveform to first post-Newtonian order beyond the quadrupole approximation and calculate new contributions from nonlinear curvature terms. We quantify the effect of these terms and provide ready-to-implement gravitational wave and scalar waveforms as well as the Fourier domain phase for quasi-circular binaries. We also perform a parameter space study, which indicates that the values of black hole scalar charges play a crucial role in the detectability of deviation from General Relativity. We also compare the scalar waveforms to numerical relativity simulations to assess the impact of the relativistic corrections to the scalar radiation. Our results provide important foundations for future precision tests of gravity., 45 pages, 5 figures, comments welcome

The metric of a spacetime can be greatly simplified if the spacetime is circular. We prove that in generic effective theories of gravity, the spacetime of a stationary, axisymmetric and asymptotically flat solution must be circular if the solution can be obtained perturbatively from a solution in the general relativity limit. This result applies to a broad class of gravitational theories that include arbitrary scalars and vectors in their light sector, so long as their nonstandard kinetic terms and nonmininal couplings to gravity are treated perturbatively., Comment: 7+5 pages, 1 figure, text and references updated, v2 matches version published in PRL

Since 2015 the gravitational-wave observations of LIGO and Virgo have transformed our understanding of compact-object binaries. In the years to come, ground-based gravitational-wave observatories such as LIGO, Virgo, and their successors will increase in sensitivity, discovering thousands of stellar-mass binaries. In the 2030s, the space-based LISA will provide gravitational-wave observations of massive black holes binaries. Between the $\sim 10$-$10^3~\mathrm{Hz}$ band of ground-based observatories and the $\sim10^{-4}$-$10^{-1}~\mathrm{Hz}$ band of LISA lies the uncharted decihertz gravitational-wave band. We propose a Decihertz Observatory to study this frequency range, and to complement observations made by other detectors. Decihertz observatories are well suited to observation of intermediate-mass ($\sim10^2$-$10^4 M_\odot$) black holes; they will be able to detect stellar-mass binaries days to years before they merge, providing early warning of nearby binary neutron star mergers and measurements of the eccentricity of binary black holes, and they will enable new tests of general relativity and the Standard Model of particle physics. Here we summarise how a Decihertz Observatory could provide unique insights into how black holes form and evolve across cosmic time, improve prospects for both multimessenger astronomy and multiband gravitational-wave astronomy, and enable new probes of gravity, particle physics and cosmology., Comment: 13 pages, 1 figure. Published in Experimental Astronomy. Summarising white paper arXiv:1908.11375

The gravitational-wave astronomical revolution began in 2015 with LIGO's observation of the coalescence of two stellar-mass black holes. Over the coming decades, ground-based detectors like LIGO will extend their reach, discovering thousands of stellar-mass binaries. In the 2030s, the space-based LISA will enable gravitational-wave observations of the massive black holes in galactic centres. Between LISA and ground-based observatories lies the unexplored decihertz gravitational-wave frequency band. Here, we propose a Decihertz Observatory to cover this band, and complement observations made by other gravitational-wave observatories. The decihertz band is uniquely suited to observation of intermediate-mass ($\sim 10^2$-$10^4 M_\odot$) black holes, which may form the missing link between stellar-mass and massive black holes, offering a unique opportunity to measure their properties. Decihertz observations will be able to detect stellar-mass binaries days to years before they merge and are observed by ground-based detectors, providing early warning of nearby binary neutron star mergers, and enabling measurements of the eccentricity of binary black holes, providing revealing insights into their formation. Observing decihertz gravitational-waves also opens the possibility of testing fundamental physics in a new laboratory, permitting unique tests of general relativity and the Standard Model of particle physics. Overall, a Decihertz Observatory will answer key questions about how black holes form and evolve across cosmic time, open new avenues for multimessenger astronomy, and advance our understanding of gravitation, particle physics and cosmology., 52 pages, 5 figures, 4 tables. Submitted to Classical & Quantum Gravity. Based upon a white paper for ESA's Voyage 2050 on behalf of the LISA Consortium 2050 Task Force

Scalar Gauss-Bonnet gravity is the only theory with quadratic curvature corrections to general relativity whose field equations are of second differential order. This theory allows for nonperturbative dynamical corrections and is therefore one of the most compelling case studies for beyond-general relativity effects in the strong-curvature regime. However, having second-order field equations is not a guarantee for a healthy time evolution in generic configurations. As a first step towards evolving black-hole binaries in this theory, we here derive the 3+1 decomposition of the field equations for any (not necessarily small) coupling constant and we discuss potential challenges of its implementation., 8 pages

We study the nonlinear dynamics of black holes that carry scalar hair and binaries composed of such black holes. The scalar hair is due to a linear or exponential coupling between the scalar and the Gauss--Bonnet invariant. We work perturbatively in the coupling constant of that interaction but nonperturbatively in the fields. We first consider the dynamical formation of hair for isolated black holes of arbitrary spin and determine the final state. This also allows us to compute for the first time the scalar quasinormal modes of rotating black holes in the presence of this coupling. We then study the evolution of nonspinning black-hole binaries with various mass ratios and produce the first scalar waveform for a coalescence. An estimate of the energy loss in scalar radiation and the effect this has on orbital dynamics and the phase of the GWs (entering at quadratic order in the coupling) show that GW detections can set the most stringent constraint to date on theories that exhibit a coupling between a scalar field and the Gauss--Bonnet invariant., Comment: 24 pages; 12 figures; Numerical code based on the Einstein Toolkit with new thorns available at https://bitbucket.org/canuda/. v2: updated references and discussion

Ultralight bosonic fields in the mass range $\sim (10^{-20}-10^{-11})\,{\rm eV}$ can trigger a superradiant instability that extracts energy and angular momentum from an astrophysical black hole with mass $M\sim(5,10^{10})M_\odot$, forming a nonspherical, rotating condensate around it. So far, most studies of the evolution and end-state of the instability have been limited to initial data containing only the fastest growing superradiant mode. By studying the evolution of multimode data in a quasi-adiabatic approximation, we show that the dynamics is much richer and depend strongly on the energy of the seed, on the relative amplitude between modes, and on the gravitational coupling. If the seed energy is a few percent of the black-hole mass, a black hole surrounded by a mixture of superradiant and nonsuperradiant modes with comparable amplitudes might not undergo a superradiant unstable phase, depending on the value of the boson mass. If the seed energy is smaller, as in the case of an instability triggered by quantum fluctuations, the effect of nonsuperradiant modes is negligible. We discuss the implications of these findings for current constraints on ultralight fields with electromagnetic and gravitational-wave observations., 21 pages, 12 figures; matches version accepted in PRD

A linear coupling between a scalar field and the Gauss-Bonnet invariant is the only known interaction term between a scalar and the metric that: respects shift symmetry; does not lead to higher order equations; inevitably introduces black hole hair in asymptotically flat, 4-dimensional spacetimes. Here we focus on the simplest theory that includes such a term and we explore the dynamical formation of scalar hair. In particular, we work in the decoupling limit that neglects the backreaction of the scalar onto the metric and evolve the scalar configuration numerically in the background of a Schwarzschild black hole or a collapsing dust star described by the Oppenheimer-Snyder solution. For all types of initial data that we consider, the scalar relaxes at late times to the known, static, analytic configuration that is associated with a hairy, spherically symmetric black hole. This suggests that the corresponding black hole solutions are indeed endpoints of collapse., Comment: 27 pages, 7 figures, contribution to the Classical and Quantum Gravity Focus Issue "Hairy black holes", v3 updated acknowledgements

Scalar fields coupled to the Gauss-Bonnet invariant evade the known no-hair theorems and have nontrivial configurations around black holes. We focus on a scalar field that couples linearly to the Gauss-Bonnet invariant and hence exhibits shift symmetry. We study its dynamical evolution and the formation of scalar hair in a Schwarzschild background. We show that the evolution eventually settles to the known static hairy solutions in the appropriate limit., Comment: 5 pages, 5 figures. v2 updated acknowledgements

Physics beyond the Standard Model is an important candidate for dark matter, and an interesting testing ground for strong-field gravity: the equivalence principle "forces" all forms of matter to fall in the same way, and it is therefore natural to look for imprints of these fields in regions with strong gravitational fields, such as compact stars or black holes. Here we study General Relativity minimally coupled to a massive vector field, and how black holes in this theory lose "hair". Our results indicate that black holes can sustain Proca field condensates for extremely long time-scales., Comment: 23 pages, 5 figures. Invited contribution to the Focus Issue on "Black holes and fundamental fields" to appear in Classical and Quantum Gravity. Matches published version

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., 188 pages, 46 figures, 6 tables, 903 references. Matches version published in Classical and Quantum Gravity. Supplementary data files available at http://www.phy.olemiss.edu/~berti/research/ and http://centra.tecnico.ulisboa.pt/network/grit/files/

We derive an initial value formulation for dynamical Chern-Simons gravity, a modification of general relativity involving parity-violating higher derivative terms. We investigate the structure of the resulting system of partial differential equations thinking about linearization around arbitrary backgrounds. This type of consideration is necessary if we are to establish well-posedness of the Cauchy problem. Treating the field equations as an effective field theory we find that weak necessary conditions for hyperbolicity are satisfied. For the full field equations we find that there are states from which subsequent evolution is not determined. Generically the evolution system closes, but the full field equations are in no sense hyperbolic. In a cursory mode analysis we find that the equations of motion contain terms that may cause ill-posedness of the initial value problem., 16 pages; matches published version

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, 51 pages, 10 figures; published version

Fundamental fields are a natural outcome in cosmology and particle physics and might therefore serve as a proxy for more complex interactions. The equivalence principle implies that all forms of matter gravitate, and one therefore expects relevant, universal imprints of new physics in strong field gravity, such as that encountered close to black holes. Fundamental fields in the vicinities of supermassive black holes give rise to extremely long-lived, or even unstable, configurations which slowly extract angular momentum from the black hole or simply evolve non-linearly over long timescales, with important implications for particle physics and gravitational-wave physics. Here, we perform a fully non-linear study of scalar-field condensates around rotating black holes. We provide novel ways to specify initial data for the Einstein-Klein-Gordon system, with potential applications in a variety of scenarios. Our numerical results confirm the existence of long-lived bar-modes which act as lighthouses for gravitational wave emission: the scalar field condenses outside the black hole geometry and acts as a constant frequency gravitational-wave source for very long timescales. This effect could turn out to be a potential signature of beyond standard model physics and also a promising source of gravitational waves for future gravitational wave detectors., 24 pages, 15 figures; v2. Overall revision including substantial improvements in the discussion of the initial data, refs added

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., 7 pages, RevTex4

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., 29 pages, 13 figures, matched to published version

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., 6 pages, 5 figures

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, nonrotating black holes of mass ratios up to 1∶100 and making contact with a classical calculation in general relativity: the infall of a pointlike particle into a massive black hole. Our results agree well with the predictions coming from linearized calculations of the infall of pointlike particles into nonrotating 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 nonlinear 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.

15 páginas., The puncture method for dealing with black holes in the numerical simulation of vacuum spacetimes is remarkably successful when combined with the Baumgarte-Shapiro-Shibata-Nakamura (BSSN) formulation of the Einstein equations. We examine a generalized class of formulations modeled along the lines of the Laguna-Shoemaker system and 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. Surprisingly, a significant portion of the parameter space yields (long-term) stable simulations, including the standard BSSN formulation as a special case. Furthermore, nonstandard parameter choices typically result in smoother behavior of the evolution variables close to the puncture.

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., Comment: 11 pages, 3 figures; corrected typos, added one Table and other minor changes. matches published version in PRD

We present unequal mass head-on collisions of black holes in D = 5 dimensional space-times. We have simulated BH systems with mass ratios q = 1,1/2, 1/3, 1/4. We extract the total energy radiated throughout the collision and compute the linear momentum flux and the recoil velocity of the final black hole. The numerical results show very good agreement with point particle calculations when extrapolated to this limit.

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 nonglobally hyperbolic space-times and the Cauchy problem may only be well defined if the initial data are supplemented by boundary conditions at the timelike 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 mirrorlike 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 observed only in the first interaction. We discuss the possible role of superradiant scattering for this result. Unlike the studies with outgoing boundary conditions, both of the Newman-Penrose scalars Psi(4) and Psi(0) are nontrivial in our setup, and we show that the numerical data verifies the expected relations between them.

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., 31 pages, 12 figures, RevTex4. v2: Published version. Further details can be found at http://blackholes.ist.utl.pt/

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., Comment: 31 pages, 6 figures; v2 Minor changes and added two references. Matches the published version in PRD.

Geodesic motion determines important features of spacetimes. Null unstable geodesics are closely related to the appearance of compact objects to external observers and have been associated with the characteristic modes of black holes. By computing the Lyapunov exponent, which is the inverse of the instability timescale associated with this geodesic motion, we show that, in the eikonal limit, quasinormal modes of black holes in any dimensions are determined by the parameters of the circular null geodesics. This result is independent of the field equations and only assumes a stationary, spherically symmetric and asymptotically flat line element, but it does not seem to be easily extendable to anti-de Sitter spacetimes. We further show that (i) in spacetime dimensions greater than four, equatorial circular timelike geodesics in a Myers-Perry black hole background are unstable, and (ii) the instability timescale of equatorial null geodesics in Myers-Perry spacetimes has a local minimum for spacetimes of dimension d > 5., 13 pages, 2 Figs, RevTex4. v2: Minor corrections. v3: more minor corrections

Black holes are among the most exciting phenomena predicted by General Relativity and play a key role in fundamental physics. Many interesting phenomena involve dynamical black hole configurations in the high curvature regime of gravity. In these lecture notes I will summarize the main numerical relativity techniques to explore highly dynamical phenomena, such as black hole collisions, in generic D-dimensional space–times. The present notes are based on my lectures given at the NR/HEP2 spring school at IST/Lisbon (Portugal) from March 11–14, 2013.

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., Comment: 68 pages, 4 Figures

The evolution of BHs in "confining boxes" is interesting for a number of reasons, particularly because it mimics some aspects of anti-de Sitter spacetimes. These admit no Cauchy surface and are a simple example of a non-globally hyperbolic spacetime. We are here interested in the potential role that boundary conditions play in the evolution of a BH system. For that, we imprison a binary BH in a box, at which boundary we set mirror-like boundary conditions.

The grand challenges of contemporary fundamental physics-dark matter, dark energy, vacuum energy, inflation and early universe cosmology, singularities and the hierarchy problem-all involve gravity as a key component. And of all gravitational phenomena, black holes stand out in their elegant simplicity, while harbouring some of the most remarkable predictions of General Relativity: event horizons, singularities and ergoregions. The hitherto invisible landscape of the gravitational Universe is being unveiled before our eyes: the historical direct detection of gravitational waves by the LIGO-Virgo collaboration marks the dawn of a new era of scientific exploration. Gravitational-wave astronomy will allow us to test models of black hole formation, growth and evolution, as well as models of gravitational-wave generation and propagation. It will provide evidence for event horizons and ergoregions, test the theory of General Relativity itself, and may reveal the existence of new fundamental fields. The synthesis of these results has the potential to radically reshape our understanding of the cosmos and of the laws of Nature. The purpose of this work is to present a concise, yet comprehensive overview of the state of the art in the relevant fields of research, summarize important open problems, and lay out a roadmap for future progress. This write-up is an initiative taken within the framework of the European Action on 'Black holes, Gravitational waves and Fundamental Physics'., SCOPUS: re.j, info:eu-repo/semantics/published