Source author record

Carlos Herdeiro

Carlos Herdeiro appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

56works
7topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

56 published item(s)

preprint2026arXiv

Stability and collisions of excited spherical boson stars: glimpses of chains and rings

Scalar, spherically symmetric, radially excited boson stars were previously shown to be stabilized, against spherical dynamics, by sufficiently strong self-interactions. Here, we further test their stability now in a full 3+1D evolution. We show that the stable stars in the former case become afflicted by a non-spherical instability. Then, we perform head-on collisions of both (stable) fundamental and (sufficiently long-lived) excited boson stars. Depending on the stars chosen, either a black hole or a bosonic remnant are possible. In particular, collisions of excited stars result in a bosonic bound state which resembles a dynamical superposition of chains and rings, akin to the ones found as equilibrium solutions in Liang:2025myf. These evolutions emphasize a key difference concerning the dynamical robustness of fundamental vs. excited spherical boson stars, when generic (beyond spherical) dynamics is considered.

preprint2022arXiv

Head-on collisions of $\ell$-boson stars

Fully non-linear numerical evolutions of the Einstein-(multi)--Klein-Gordon equations are performed to study head-on collisions of $\ell$-boson stars. Despite being spherically symmetric, $\ell$-boson stars have a (hidden) frame of reference, used in defining their individual multipolar fields. To assess the impact of their relative orientation, we perform simulations with different angles between the axes of the two colliding stars. Additionally, two scenarios are considered for the colliding stars: that they are composites of either the same or different scalar fields. Despite some model-specific behaviours, the simulations generically indicate that: 1) the collision of two sufficiently (and equally) massive stars leads to black hole formation; 2) below a certain mass threshold the end result of the evolution is a bound state of the composite scalar fields, that neither disperses nor collapses into a black hole within the simulation time; 3) this end product (generically) deviates from spherical symmetry and the equipartition of the number of bosonic particles between the different scalar fields composing the initial boson stars is lost, albeit not dramatically. This last observation indicates, albeit without being conclusive, that the end result of these collisions belongs to the previously reported larger family of equilibrium multi-field boson stars, generically non-spherical, and of which $\ell$-boson stars are a symmetry enhanced point. We also extract and discuss the waveforms from the collisions studied.

preprint2022arXiv

New Horizons for Fundamental Physics with LISA

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 GWs 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.

preprint2021arXiv

GW190521 as a merger of Proca stars: a potential new vector boson of $8.7 \times 10^{-13}$ eV

Advanced LIGO-Virgo reported a short gravitational-wave signal (GW190521) interpreted as a quasi-circular merger of black holes, one populating the pair-instability supernova gap, forming a remnant black hole of $M_f\sim 142 M_\odot$ at a luminosity distance of $d_L \sim 5.3$ Gpc. With barely visible pre-merger emission, however, GW190521 merits further investigation of the pre-merger dynamics and even of the very nature of the colliding objects. We show that GW190521 is consistent with numerically simulated signals from head-on collisions of two (equal mass and spin) horizonless vector boson stars (aka Proca stars), forming a final black hole with $M_f = 231^{+13}_{-17}\,M_\odot$, located at a distance of $d_L = 571^{+348}_{-181}$ Mpc. The favoured mass for the ultra-light vector boson constituent of the Proca stars is $μ_{\rm V}= 8.72^{+0.73}_{-0.82}\times10^{-13}$ eV. This provides the first demonstration of close degeneracy between these two theoretical models, for a real gravitational-wave event. Confirmation of the Proca star interpretation, which we find statistically slightly preferred, would provide the first evidence for a long sought dark matter particle.

preprint2021arXiv

Self-interactions can stabilize excited boson stars

We study the time evolution of spherical, excited ($i.e.$ nodeful) boson star models. We consider a model including quartic self-interactions, controlled by a coupling $Λ$. Performing non-linear simulations of the Einstein-(complex)-Klein-Gordon system, using as initial data equilibrium boson stars solutions of that system, we assess the impact of $Λ$ in the stability properties of the boson stars. In the absence of self-interactions ($Λ=0$), we observe the known behaviour that the excited stars in the (candidate) stable branch decay to a non-excited star without a node; however, we show that for large enough values of the self-interactions coupling, these excited stars do not decay (up to timescales of about $t\sim 10^4$). The stabilization of the excited states for large enough self-interactions is further supported by evidence that the nodeful states dynamically form through the gravitational cooling mechanism, starting from dilute initial data. Our results support the healing power (against dynamical instabilities) of self-interactions, recently unveiled in the context of the non-axisymmetric instabilities of spinning boson stars.

preprint2021arXiv

Spinning gauged boson and Dirac stars: a comparative study

Scalar boson stars and Dirac stars are solitonic solutions of the Einstein{Klein-Gordon and Einstein- Dirac classical equations, respectively. Despite the different bosonic vs: fermionic nature of the matter field, these solutions to the classical field equations have been shown to have qualitatively similar features [1]. In particular, for spinning stars the most fundamental configurations can be in both cases toroidal, unlike spinning Proca stars that are spheroidal [2]. In this paper we gauge the scalar and Dirac fields, by minimally coupling them to standard electromagnetism. We explore the impact of the gauge coupling on the resulting solutions. One of the most relevant difference concerns the gyromagnetic ratio, which for the scalar stars takes values around 1, whereas for Dirac stars takes values around 2.

preprint2020arXiv

Dynamical $\ell$-boson stars: generic stability and evidence for non-spherical solutions

$\ell$-boson stars are static, spherical, multi-field self-gravitating solitons. They are asymptotically flat, finite energy solutions of Einstein's gravity minimally coupled to an odd number of massive, complex scalar fields. A previous study assessed the stability of $\ell$-boson stars under spherical perturbations, finding that there are both stable and unstable branches of solutions, as for single-field boson stars ($\ell=0$). In this work we probe the stability of $\ell$-boson stars against non-spherical perturbations by performing numerical evolutions of the Einstein-Klein-Gordon system, with a 3D code. For the timescales explored, the $\ell$-boson stars belonging to the spherical stable branch do not exhibit measurable growing modes. We find, however, evidence of zero modes; that is, non-spherical perturbations that neither grow nor decay. This suggests the branching off towards a larger family of equilibrium solutions: we conjecture that $\ell$-boson stars are the enhanced isometry point of a larger family of static (and possibly stationary), non-spherical multi-field self-gravitating solitons.

preprint2020arXiv

Static black holes without spatial isometries: from AdS multipoles to electrovacuum scalarization

We review recent results on the existence of static black holes without spatial isometries in four spacetime dimensions and propose a general framework for their study. These configurations are regular on and outside a horizon of spherical topology. Two different mechanisms allowing for their existence are identified. The first one relies on the presence of a solitonic limit of the black holes; when the solitons have no spatial isometries, the black holes, being a non-linear bound state between the solitons and a horizon, inherit this property. The second one is related to black hole scalarization, and the existence of zero modes of the scalar field without isometries around a spherical horizon. When the zero modes have no spatial isometries, the backreaction of their non-linear continuation makes the scalarized black holes inherit the absence of spatial continuous symmetries. A number of general features of the solutions are discussed together with possible generalizations.

preprint2019arXiv

Black hole spontaneous scalarisation with a positive cosmological constant

A scalar field non-minimally coupled to certain geometric [or matter] invariants which are sourced by [electro]vacuum black holes (BHs) may spontaneously grow around the latter, due to a tachyonic instability. This process is expected to lead to a new, dynamically preferred, equilibrium state: a scalarised BH. The most studied geometric [matter] source term for such spontaneous BH scalarisation is the Gauss-Bonnet quadratic curvature [Maxwell invariant]. This phenomenon has been mostly analysed for asymptotically flat spacetimes. Here we consider the impact of a positive cosmological constant, which introduces a cosmological horizon. The cosmological constant does not change the local conditions on the scalar coupling for a tachyonic instability of the scalar-free BHs to emerge. But it leaves a significant imprint on the possible new scalarised BHs. It is shown that no scalarised BH solutions exist, under a smoothness assumption, if the scalar field is confined between the BH and cosmological horizons. Admitting the scalar field can extend beyond the cosmological horizon, we construct new scalarised BHs. These are asymptotically de Sitter in the (matter) Einstein-Maxwell-scalar model, with only mild difference with respect to their asymptotically flat counterparts. But in the (geometric) extended-scalar-tensor-Gauss-Bonnet-scalar model, they have necessarily non-standard asymptotics, as the tachyonic instability dominates in the far field. This interpretation is supported by the analysis of a test tachyon on a de Sitter background.

preprint2019arXiv

Charged Dirac perturbations on Reissner-Nordström-Anti-de Sitter spacetimes: quasinormal modes with Robin boundary conditions

We study charged Dirac quasinormal modes (QNMs) on Reissner-Nordström-Anti-de Sitter (RN-AdS) black holes with generic Robin boundary conditions, by extending our earlier work of neutral Dirac QNMs on Schwarzschild-AdS black holes. We first derive the equations of motion for charged Dirac fields on a RN-AdS background. To solve these equations we impose a requirement on the Dirac field: that its energy flux should vanish at asymptotic infinity. A set of two Robin boundary conditions compatible with QNMs is consequently found. By employing both analytic and numeric methods, we then obtain the quasinormal spectrum for charged Dirac fields, and analyse the impact of various parameters, in particular of electric charge. An analytic calculation shows explicitly that the charge coupling between the black hole and the Dirac field does not trigger superradiant instabilities. Numeric calculations, on the other hand, show quantiatively that Dirac QNMs may change substantially due to the electric charge. Our results illustrate how vanishing energy flux boundary conditions, as a generic principle, are applicable not only to neutral but also to electrically charged fields.

preprint2019arXiv

Dynamically and thermodynamically stable black holes in Einstein-Maxwell-dilaton gravity

We consider Einstein-Maxwell-dilaton gravity with the non-minimal exponential coupling between the dilaton and the Maxwell field emerging from low energy heterotic string theory. The dilaton is endowed with a potential that originates from an electromagnetic Fayet-Iliopoulos (FI) term in $\mathcal{N}=2$ extended supergravity in four spacetime dimensions. For the case we are interested in, this potential introduces a single parameter $α$. When $α\rightarrow 0$, the static black holes (BHs) of the model are the Gibbons-Maeda-Garfinkle-Horowitz-Strominger (GMGHS) solutions. When $α\rightarrow \infty$, the BHs become the standard Reissner-Nordström (RN) solutions of electrovacuum General Relativity. The BH solutions for finite non-zero $α$ interpolate between these two families. In this case, the dilaton potential regularizes the extremal limit of the GMGHS solution yielding a set of zero temperature BHs with a near horizon $AdS_2\times S^2$ geometry. We show that, in the neighborhood of these extremal solutions, there is a subset of BHs that are dynamically and thermodynamically stable, all of which have charge to mass ratio larger than unity. By dynamical stability we mean that no growing quasi-normal modes are found; thus they are stable against linear perturbations (spherical and non-spherical). Moreover, non-linear numerical evolutions lend support to their non-linear stability. By thermodynamical stability we mean the BHs are stable both in the canonical and grand-canonical ensemble. In particular, both the specific heat at constant charge and the isothermal permittivity are positive. This is not possible for RN and GMGHS BHs. We discuss the different thermodynamical phases for the BHs in this model and comment on what may allow the existence of both dynamically and thermodynamically stable BHs.

preprint2016arXiv

Dynamical formation of a Reissner-Nordström black hole with scalar hair in a cavity

In a recent letter, we presented numerical relativity simulations, solving the full Einstein--Maxwell--Klein-Gordon equations, of superradiantly unstable Reissner-Nordström black holes (BHs), enclosed in a cavity. Low frequency, spherical perturbations of a charged scalar field, trigger this instability. The system's evolution was followed into the non-linear regime, until it relaxed into an equilibrium configuration, found to be a $\textit{hairy}$ BH: a charged horizon in equilibrium with a scalar field condensate, whose phase is oscillating at the (final) critical frequency. Here, we investigate the impact of adding self-interactions to the scalar field. In particular, we find sufficiently large self-interactions suppress the exponential growth phase, known from linear theory, and promote a non-monotonic behaviour of the scalar field energy. Furthermore, we discuss in detail the influence of the various parameters in this model: the initial BH charge, the initial scalar perturbation, the scalar field charge, mass, and the position of the cavity's boundary (mirror). We also investigate the "explosive" non-linear regime previously reported to be akin to a bosenova. A mode analysis shows that the "explosions" can be interpreted as the decay into the BH of modes that exit the superradiant regime.

preprint2016arXiv

Einstein-Maxwell-Anti-de-Sitter spinning solitons

Electrostatics on global Anti-de-Sitter (AdS) spacetime is sharply different from that on global Minkowski spacetime. It admits a multipolar expansion with everywhere regular, finite energy solutions, for every multipole moment except the monopole (arXiv:1507.04370). A similar statement holds for global AdS magnetostatics. We show that everywhere regular, finite energy, electric plus magnetic fields exist on AdS in three distinct classes: $(I)$ with non-vanishing total angular momentum $J$; $(II)$ with vanishing $J$ but non-zero angular momentum density, $T^t_φ$; $(III)$ with vanishing $J$ and $T^t_φ$. Considering backreaction, these configurations remain everywhere smooth and finite energy, and we find, for example, Einstein-Maxwell-AdS solitons that are globally - Type I - or locally (but not globally) - Type II - spinning. This backreaction is considered first perturbatively, using analytical methods and then non-perturbatively, by constructing numerical solutions of the fully non-linear Einstein-Maxwell-AdS system. The variation of the energy and total angular momentum with the boundary data is explicitly exhibited for one example of a spinning soliton.

preprint2016arXiv

Explosion and final state of an unstable Reissner-Nordstrom black hole

A Reissner-Nordström black hole (BH) is superradiantly unstable against spherical perturbations of a charged scalar field, enclosed in a cavity, with frequency lower than a critical value. We use numerical relativity techniques to follow the development of this unstable system -- dubbed a charged BH bomb -- into the non-linear regime, solving the full Einstein--Maxwell--Klein-Gordon equations, in spherical symmetry. We show that: $i)$ the process stops before all the charge is extracted from the BH; $ii)$ the system settles down into a hairy BH: a charged horizon in equilibrium with a scalar field condensate, whose phase is oscillating at the (final) critical frequency. For low scalar field charge, $q$, the final state is approached smoothly and monotonically. For large $q$, however, the energy extraction overshoots and an explosive phenomenon, akin to a $bosenova$, pushes some energy back into the BH. The charge extraction, by contrast, does not reverse.

preprint2016arXiv

Inside black holes with synchronized hair

Recently, various examples of asymptotically flat, rotating black holes (BHs) with synchronized hair have been explicitly constructed, including Kerr BHs with scalar or Proca hair, and Myers-Perry BHs with scalar hair and a mass gap, showing there is a general mechanism at work. All these solutions have been found numerically, integrating the fully non-linear field equations of motion from the event horizon outwards. Here, we address the spacetime geometry of these solutions inside the event horizon. Firstly, we provide arguments, within linear theory, that there is no regular inner horizon for these solutions. Then, we address this question fully non-linearly, using as a tractable model five dimensional, equal spinning, Myers-Perry hairy BHs. We find that, for non-extremal solutions: $(1)$ the inside spacetime geometry in the vicinity of the event horizon is smooth and the equations of motion can be integrated inwards; $(2)$ before an inner horizon is reached, the spacetime curvature grows (apparently) without bound. In all cases, our results suggest the absence of a smooth Cauchy horizon, beyond which the metric can be extended, for hairy BHs with synchronized hair.

preprint2016arXiv

Kerr black holes with Proca hair

Bekenstein proved that in Einstein's gravity minimally coupled to one (or many) real, Abelian, Proca field, stationary black holes (BHs) cannot have Proca hair. Dropping Bekenstein's assumption that matter inherits spacetime symmetries, we show this model admits asymptotically flat, stationary, axi-symmetric, regular on and outside an event horizon BHs with Proca hair, for an even number of real (or an arbitrary number of complex) Proca fields. To establish it, we start by showing that a test, complex Proca field can form bound states, with real frequency, around Kerr BHs: stationary Proca clouds. These states exist at the threshold of superradiance. It was conjectured in arXiv:1403.2757, that the existence of such clouds at the linear level implies the existence of a new family of BH solutions at the non-linear level. We confirm this expectation and explicitly construct examples of such Kerr black holes with Proca hair (KBHsPH). For a single complex Proca field, these BHs form a countable number of families with three continuous parameters (ADM mass, ADM angular momentum and Noether charge). They branch off from the Kerr solutions that can support stationary Proca clouds and reduce to Proca stars when the horizon size vanishes. We present the domain of existence of one family of KBHsPH, as well as its phase space in terms of ADM quantities. Some physical properties of the solutions are discussed; in particular, and in contrast with Kerr BHs with scalar hair, some spacetime regions can be counter-rotating with respect to the horizon. We further establish a no-Proca-hair theorem for static, spherically symmetric BHs but allowing the complex Proca field to have a harmonic time dependence, which shows BHs with Proca hair in this model require rotation and have no static limit. KBHsPH are also disconnected from Kerr-Newman BHs with a real, massless vector field.

preprint2016arXiv

Numerical Relativity and High Energy Physics: Recent Developments

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.

preprint2015arXiv

Acoustic clouds: standing sound waves around a black hole analogue

Under certain conditions sound waves in fluids experience an acoustic horizon with analogue properties to those of a black hole event horizon. In particular, a draining bathtub-like model can give rise to a rotating acoustic horizon and hence a rotating black hole (acoustic) analogue. We show that sound waves, when enclosed in a cylindrical cavity, can form stationary waves around such rotating acoustic black holes. These acoustic perturbations display similar properties to the scalar clouds that have been studied around Kerr and Kerr-Newman black holes; thus they are dubbed acoustic clouds. We make the comparison between scalar clouds around Kerr black holes and acoustic clouds around the draining bathtub explicit by studying also the properties of scalar clouds around Kerr black holes enclosed in a cavity. Acoustic clouds suggest the possibility of testing, experimentally, the existence and properties of black hole clouds, using analog models.

preprint2015arXiv

Anti-de-Sitter regular electric multipoles: towards Einstein-Maxwell-AdS solitons

We discuss electrostatics in Anti-de-Sitter ($AdS$) spacetime, in global coordinates. We observe that the multipolar expansion has two crucial differences to that in Minkowski spacetime. First, there are everywhere regular solutions, with finite energy, for every multipole moment except for the monopole. Second, all multipole moments decay with the same inverse power of the areal radius, $1/r$, as spatial infinity is approached. The first observation suggests there may be regular, self-gravitating, Einstein-Maxwell solitons in $AdS$ spacetime. The second observation, renders a Lichnerowicz-type no-soliton theorem inapplicable. Consequently, we suggest Einstein-Maxwell solitons exist in $AdS$, and we support this claim by computing the first order metric perturbations sourced by test electric field multipoles, which are obtained analytically in closed form.

preprint2015arXiv

Construction and physical properties of Kerr black holes with scalar hair

Kerr black holes with scalar hair are solutions of the Einstein-Klein-Gordon field equations describing regular (on and outside an event horizon), asymptotically flat black holes with scalar hair (arXiv:1403.2757). These black holes interpolate continuously between the Kerr solution and rotating boson stars in D=4 spacetime dimensions. Here we provide details on their construction, discussing properties of the ansatz, the field equations, the boundary conditions and the numerical strategy. Then, we present an overview of the parameter space of the solutions, and describe in detail the space-time structure of the black holes exterior geometry and of the scalar field for a sample of reference solutions. Phenomenological properties of potential astrophysical interest are also discussed, and the stability properties and possible generalizations are commented on. As supplementary material to this paper we make available numerical data files for the sample of reference solutions discussed, for public use.

preprint2015arXiv

Exploring New Physics Frontiers Through Numerical Relativity

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.

preprint2015arXiv

Maxwell perturbations on asymptotically anti-de Sitter spacetimes: Generic boundary conditions and a new branch of quasinormal modes

Perturbations of asymptotically Anti-de-Sitter (AdS) spacetimes are often considered by imposing field vanishing boundary conditions (BCs) at the AdS boundary. Such BCs, of Dirichlet-type, imply a vanishing energy flux at the boundary, but the converse is, generically, not true. Regarding AdS as a gravitational box, we consider vanishing energy flux (VEF) BCs as a more fundamental physical requirement and we show that these BCs can lead to a new branch of modes. As a concrete example, we consider Maxwell perturbations on Kerr-AdS black holes in the Teukolsky formalism, but our formulation applies also for other spin fields. Imposing VEF BCs, we find a set of two Robin BCs, even for Schwarzschild-AdS black holes. The Robin BCs on the Teukolsky variables can be used to study quasinormal modes, superradiant instabilities and vector clouds. As a first application, we consider here the quasinormal modes of Schwarzschild-AdS black holes. We find that one of the Robin BCs yields the quasinormal spectrum reported in the literature, while the other one unveils a new branch for the quasinormal spectrum.

preprint2015arXiv

Maxwell perturbations on Kerr-anti-de Sitter: quasinormal modes, superradiant instabilities and vector clouds

Scalar and gravitational perturbations on Kerr-anti-de Sitter (Kerr-AdS) black holes have been addressed in the literature and have been shown to exhibit a rich phenomenology. In this paper we complete the analysis of bosonic fields on this background by studying Maxwell perturbations, focusing on superradiant instabilities and vector clouds. For this purpose, we solve the Teukolsky equations numerically, imposing the boundary conditions we have proposed in\cite{Wang:2015goa} for the radial Teukolsky equation. As found therein, two Robin boundary conditions can be imposed for Maxwell fields on Kerr-AdS black holes, one of which produces a new set of quasinormal modes even for Schwarzschild-AdS black holes. Here, we show these different boundary conditions produce two different sets of superradiant modes. Interestingly the "new modes" may be unstable in a larger parameter space. We then study stationary Maxwell clouds, that exist at the threshold of the superradiant instability, with the two Robin boundary conditions. These clouds, obtained at the linear level, indicate the existence of a new family of black hole solutions at the nonlinear level, within the Einstein-Maxwell-AdS system, branching off from the Kerr-Newman-AdS family. As a comparison with the Maxwell clouds, scalar clouds on Kerr-AdS black holes are also studied, and it is shown there are Kerr-AdS black holes that are stable against scalar, but not vector modes, with the same "quantum numbers".

preprint2015arXiv

Myers-Perry black holes with scalar hair and a mass gap: unequal spins

We construct rotating boson stars and Myers-Perry black holes with scalar hair (MPBHsSH) as fully non-linear solutions of five dimensional Einstein gravity minimally coupled to a complex, massive scalar field. The MPBHsSH are, in general, regular on and outside the horizon, asymptotically flat, and possess angular momentum in a single rotation plane. They are supported by rotation and have no static limit. Such hairy BHs may be thought of as bound states of boson stars and singly spinning, vacuum MPBHs and inherit properties of both these building blocks. When the horizon area shrinks to zero, the solutions reduce to (in a single plane) rotating boson stars; but the extremal limit also yields a zero area horizon, as for singly spinning MPBHs. Similarly to the case of equal angular momenta, and in contrast to Kerr black holes with scalar hair, singly spinning MPBHsSH are disconnected from the vacuum black holes, due to a mass gap. We observe that for the general case, with two unequal angular momenta, the equilibrium condition for the existence of MPBHsSH is $w=m_1Ω_1+m_2Ω_2$, where $Ω_i$ are the horizon angular velocities in the two independent rotation planes and $w,m_i$, $i=1,2$, are the scalar field's frequency and azimuthal harmonic indices.

preprint2015arXiv

Testing General Relativity with Present and Future Astrophysical Observations

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.

preprint2014arXiv

Collisions of oppositely charged black holes

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.

preprint2014arXiv

Ergo-spheres, ergo-tori and ergo-Saturns for Kerr black holes with scalar hair

We have recently reported (arXiv:1403.2757) the existence of Kerr black holes with scalar hair in General Relativity minimally coupled to a massive, complex scalar field. These solutions interpolate between boson stars and Kerr black holes. The latter have a well known topologically S^2 ergo-surface (ergo-sphere) whereas the former develop a S^1*S^1 ergo-surface (ergo-torus) in a region of parameter space. We show that hairy black holes always have an ergo-region, and that this region is delimited by either an ergo-sphere or an ergo-Saturn -- i.e. a S^2\oplus (S^1*S^1) ergo-surface. In the phase space of solutions, the ergo-torus can either appear disconnected from the ergo-sphere or pinch off from it. We provide a heuristic argument, based on a measure of the size of the ergo-region, that superradiant instabilities - which are likely to be present - are weaker for hairy black holes than for Kerr black holes with the same global charges. We observe that Saturn-like, and even more remarkable ergo-surfaces, should also arise for other rotating `hairy' black holes.

preprint2014arXiv

Higher dimensional Numerical Relativity: code comparison

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.

preprint2014arXiv

Kerr-Newman scalar clouds

Massive complex scalar fields can form bound states around Kerr black holes. These bound states -- dubbed scalar clouds -- are generically non-zero and finite on and outside the horizon; they decay exponentially at spatial infinity, have a real frequency and are specified by a set of integer "quantum" numbers (n,l,m). For a specific set of these numbers, the clouds are only possible along a 1-dimensional subset of the 2-dimensional parameter space of Kerr black holes, called an existence line. In this paper we make a thorough investigation of the scalar clouds due to neutral (charged) scalar fields around Kerr(-Newman) black holes. We present the location of the existence lines for a variety of quantum numbers, their spatial representation and compare analytic approximation formulas in the literature with our exact numerical results, exhibiting a sometimes remarkable agreement.

preprint2014arXiv

Marginal scalar and Proca clouds around Reissner-Nordström black holes

Massive scalar test fields around Kerr black holes can form quasi-bound states with complex frequencies. Some of these states decay in time, but some others -- in the superradiant regime -- grow, causing a superradiant instability. Precisely at the threshold between decaying and growing modes, there exist bound states with real frequency, known as scalar clouds. Fully non-linear counterparts of these clouds have been shown to yield Kerr black holes with scalar hair. Here, we consider massive, (electrically) charged, test scalar and Proca fields on the background of Reissner-Nordström black holes. By computing the quasi-bound state frequencies, we establish that no such states exist in the superradiant regime for the Proca field -- a similar behavior to that known for scalar fields. But decaying quasi-bound states with an arbitrarily small imaginary part of the frequency exist and thus which are arbitrarily long lived. In the limit of vanishing imaginary part of the frequency, the field does not trivialize and we dub the corresponding configurations as marginal (charged) scalar or Proca clouds, since they are only marginally bound. We comment on the possible existence of non-linear counterparts to these marginal clouds.

preprint2014arXiv

Myers-Perry black holes with scalar hair and a mass gap

We construct a family of asymptotically flat, rotating black holes with scalar hair and a regular horizon, within five dimensional Einstein's gravity minimally coupled to a complex, massive scalar field doublet. These solutions are supported by rotation and have no static limit. They are described by their mass $M$, two equal angular momenta $J_1=J_2\equiv J$ and a conserved Noether charge $Q$, measuring the scalar hair. For vanishing horizon size the solutions reduce to five dimensional boson stars. In the limit of vanishing Noether charge density, the scalar field becomes point-wise arbitrarily small and the geometry becomes, locally, arbitrarily close to that of a specific set of Myers-Perry black holes (MPBHs); but there remains a global difference with respect to the latter, manifest in a finite mass gap. Thus, the scalar hair never becomes a linear perturbation of MPBHs. This is a qualitative difference when compared to Kerr black holes with scalar hair~\cite{Herdeiro:2014goa}. Whereas the existence of the latter can be anticipated in linear theory, from the existence of scalar bound states on the Kerr geometry (i.e. scalar clouds), the hair of these MPBHs is intrinsically non-linear.

preprint2014arXiv

Non-linear Q-clouds around Kerr black holes

Q-balls are regular extended `objects' that exist for some non-gravitating, self-interacting, scalar field theories with a global, continuous, internal symmetry, on Minkowski spacetime. Here, analogous objects are also shown to exist around rotating (Kerr) black holes, as non-linear bound states of a test scalar field. We dub such configurations Q-clouds. We focus on a complex massive scalar field with quartic plus hexic self-interactions. Without the self-interactions, linear clouds have been shown to exist, in synchronous rotation with the black hole horizon, along 1-dimensional subspaces - existence lines - of the Kerr 2-dimensional parameter space. They are zero modes of the superradiant instability. Non-linear Q-clouds, on the other hand, are also in synchronous rotation with the black hole horizon; but they exist on a 2-dimensional subspace, delimited by a minimal horizon angular velocity and by an appropriate existence line, wherein the non-linear terms become irrelevant and the Q-cloud reduces to a linear cloud. Thus, Q-clouds provide an example of scalar bound states around Kerr black holes which, generically, are not zero modes of the superradiant instability. We describe some physical properties of Q-clouds, whose backreaction leads to a new family of hairy black holes, continuously connected to the Kerr family.

preprint2014arXiv

PPN expansion and FRW scalar perturbations in n-DBI gravity

n-DBI gravity explicitly breaks Lorentz invariance by the introduction of a unit time-like vector field, thereby giving rise to an extra (scalar) degree of freedom. We look for observational consequences of this mode in two setups. Firstly, we compute the parametrized post-Newtonian (PPN) expansion of the metric to first post-Newtonian order. Surprisingly, we find that the PPN parameters are exactly the same as in General Relativity (GR), and no preferred-frame effects are produced. In particular this means that n-DBI gravity is consistent with all GR solar system experimental tests. We discuss the origin of such degeneracy between n-DBI gravity and GR, and suggest it may also hold in higher post-Newtonian order. Secondly, we study gravitational scalar perturbations of a Friedmann-Robertson-Walker space-time with a cosmological constant $Λ\geq 0$. In the case of de Sitter space, we show that the scalar mode grows as the universe expands and, in contrast with a canonical scalar field coupled to GR, it does not freeze on super horizon scales.

preprint2014arXiv

Radiation from a $D$-dimensional collision of shock waves: proof of first order formula and angular factorisation at all orders

In two previous papers we have computed the inelasticity $ε$ in a head-on collision of two $D$-dimensional Aichelburg-Sexl shock waves, using perturbation theory to calculate the geometry in the future light-cone of the collision. The first order result, obtained as an accurate numerical fit, yielded the remarkably simple formula $ε_{\rm 1st\, order} = 1/2 - 1/D$. Here we show, analytically, that this result is exact in first order perturbation theory. Moreover, we clarify the relation between perturbation theory and an angular series of the inelasticity's angular power around the symmetry axis of the collision $(θ= 0,π)$. To establish these results, firstly, we show that at null infinity the angular dependence factorises order by order in perturbation theory, as a result of a hidden symmetry. Secondly, we show that a consistent truncation of the angular series in powers of $\sin^2 θ$ at some order $O(n)$ requires knowledge of the metric perturbations up to $O(n+1)$. In particular, this justifies the isotropy assumption used in first order perturbation theory. We then compute, analytically, all terms that contribute to the inelasticity and depend linearly on the initial conditions (surface terms), including second order contributions.

preprint2014arXiv

Superradiant instabilities in a $D$-dimensional small Reissner-Nordström-Anti-de Sitter black hole

We investigate the superradiant instability for a charged scalar field in a $D$-dimensional small Reissner-Nordström-Anti-de Sitter (RN-AdS) black hole. Firstly, we solve the charged Klein-Gordon equation analytically by a matching method. We show that the general $D$-dimensional quasinormal frequencies depend on the relation between the angular momentum quantum number, $\ell$, and $D$. When $\ell$ is a (non-negative) integer multiple of $D-3$, i.e $\ell=p(D-3)$, we find an analytical quasinormal frequency formula adding a purely imaginary correction to the AdS normal frequencies. This is the case for all $\ell$ modes in $D=4$. For general $D$ there are two more cases: i) when $\ell$ obeys $\ell=(p+\frac{1}{2})(D-3)$, which can occur for odd $D$, we observe that the matching method fails, since the near and far region solutions have different functional dependences; ii) for all other cases, the analytical quasinormal frequency formula gives a complex correction to the AdS normal frequencies. Secondly, we perform a numerical calculation which confirms the analytical formulas obtained with the matching method and allows us to explore the case where that method failed. In the latter case, as in the former, we verify that all $\ell$ modes for all $D$ may become superradiant, which contradicts a previous claim in the literature.

preprint2014arXiv

Testing the nonlinear stability of Kerr-Newman black holes

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.

preprint2013arXiv

Hawking radiation for a Proca field in D dimensions II: charged field in a brane charged black hole

We generalise our first analysis of the wave equation for a massive vector boson in the background of a D-dimensional Schwarzschild black hole, by adding charge both to the field and the black hole, on the 3+1 dimensional Standard Model brane. A detailed numerical study is performed to obtain the transmission factor for the coupled (as well as decoupled) system of equations describing the Proca field modes, varying the angular momentum number, mass, charge and space-time dimensions. A qualitatively new feature arising from the introduction of charge is the appearance of superradiant modes, which we investigate. We then compute the Hawking fluxes, and analyse the effect of the charge. In particular we observe an inverted charge splitting effect for small energies and for two or more extra dimensions. For neutral fields, we compare the emission of massive particles with spins up to one and also compare the Proca bulk-to-brane ratio of energy emission, showing that, as for a scalar field, most of the energy is emitted on the brane.

preprint2013arXiv

n-DBI Gravity in a nutshell

We present a new model of gravity which explicitly breaks Lorentz-invariance by the introduction of a unit time-like vector field, thereby giving rise to an extra (scalar) degree of freedom. We discuss its cosmology, exact solutions and the dynamics of the scalar mode. We show that it predicts inflation without an inflaton and admits the black hole solutions of General Relativity (GR). We argue that the scalar mode is well behaved and contains none of the pathologies previously found in similar models.

preprint2013arXiv

n-DBI gravity, maximal slicing and the Kerr geometry

Recently, in arXiv:1110.0832, we have established that solutions of Einstein's gravity admitting foliations with a certain geometric condition are also solutions of n-DBI gravity, arXiv:1109.1468. Here we observe that, in vacuum, the required geometric condition is fulfilled by the well known maximal slicing, often used in numerical relativity. As a corollary, we establish that the Kerr geometry is a solution of n-DBI gravity in the foliation adapted to Boyer-Lindquist coordinates.

preprint2013arXiv

Radiation from a D-dimensional collision of shock waves: higher order set up and perturbation theory validity

The collision of two D-dimensional, ultra-relativistic particles, described in General Relativity as Aichelberg-Sexl shock waves, is inelastic. In first order perturbation theory, the fraction of the initial centre of mass energy radiated away was recently shown to be 1/2 - 1/D. Here, we extend the formalism to higher orders in perturbation theory, and derive a general expression to extract the inelasticity, valid non-perturbatively, based on the Bondi mass loss formula. Then, to clarify why perturbation theory captures relevant physics of a strong field process in this problem, we provide one variation of the problem where the perturbative framework breaks down: the collision of ultra-relativistic charged particles. The addition of charge, and the associated repulsive nature of the source, originates an extra radiation burst, which we argue to be an artifact of the perturbative framework, veiling the relevant physics.

preprint2012arXiv

Collisions of charged black holes

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.

preprint2012arXiv

Dynamics of black holes in de Sitter spacetimes

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.

preprint2012arXiv

NR/HEP: roadmap for the future

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.

preprint2012arXiv

On the scalar graviton in n-DBI gravity

n-DBI gravity is a gravitational theory which yields near de Sitter inflation spontaneously at the cost of breaking Lorentz invariance by a preferred choice of foliation. We show that this breakdown endows n-DBI gravity with one extra physical gravitational degree of freedom: a scalar graviton. Its existence is established by Dirac's theory of constrained systems. Firstly, studying scalar perturbations around Minkowski space-time, we show that there exists one scalar degree of freedom and identify it in terms of the metric perturbations. Then, a general analysis is made in the canonical formalism, using ADM variables. It is useful to introduce an auxiliary scalar field, which allows recasting n-DBI gravity in an Einstein-Hilbert form but in a Jordan frame. Identifying the constraints and their classes we confirm the existence of an extra degree of freedom in the full theory, besides the two usual tensorial modes of the graviton. We then argue that, unlike the case of (the original proposal for) Horava-Lifschitz gravity, there is no evidence that the extra degree of freedom originates pathologies, such as vanishing lapse, instabilities and strong self-coupling at low energy scales.

preprint2012arXiv

Radiation from a D-dimensional collision of shock waves: a remarkably simple fit formula

Recently, in arXiv:1105.2298 [hep-th], we have estimated the energy radiated in the head-on collision of two equal D-dimensional Aichelburg-Sexl shock waves, for even D, by solving perturbatively, to first order, the Einstein equations in the future of the collision. Here, we report on the solution for the odd D case. After finding the wave forms, we extract the estimated radiated energy for D=5,7,9 and 11 and unveil a remarkably simple pattern, given the complexity of the framework: (for all D) the estimated fraction of radiated energy matches the analytic expression 1/2-1/D, within the numerical error (less than 0.1%). Both this fit and the apparent horizon bound converge to 1/2 as D goes to infinity.

preprint2012arXiv

Scale invariance and a gravitational model with non-eternal inflation

We propose a 3 + 1 dimensional model of gravity which results in inflation at early times, followed by radiation- and matter-dominated epochs and a subsequent acceleration at late times. Both the inflation and late time acceleration are nearly de Sitter with a large hierarchy between the effective cosmological constants. There is no scalar field agent of inflation, and the transition from the inflation to the radiation-dominated period is smooth. This model is designed so that it yields, at the cost of giving up on Lorentz invariance in the gravitational sector, the Dirac-Born-Infeld type conformal scalar theory when the universe is conformally flat. It, however, resembles Einstein's gravity with the Gibbons-Hawking-York boundary term in weakly curved space-times.

preprint2011arXiv

Hawking radiation for a Proca field in D-dimensions

We study the wave equation of a massive vector boson in the background of a D-dimensional Schwarzschild black hole. The mass term introduces a coupling between two physical degrees of freedom of the field, and we solve the resulting system of ODEs numerically, without decoupling. We show how to define decoupled transmission factors from an S-matrix and compute them for various modes, masses and space-time dimensions. The mass term lifts the degeneracy between transverse modes, in D=4, and excites the longitudinal modes, in particular the s-wave. Moreover, it increases the contribution of waves with larger angular momentum, which can be dominant at intermediate energies. The transmission factors are then used to obtain the Hawking fluxes in this channel. Our results alert for the importance of modelling the longitudinal modes correctly, instead of treating them as decoupled scalars as in current black hole event generators; thus they can be used to improve such generators for phenomenological studies of TeV gravity scenarios.

preprint2011arXiv

Higher-dimensional puncture initial data

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.

preprint2011arXiv

n-DBI gravity

n-DBI gravity is a gravitational theory introduced in arXiv:1109.1468 [hep-th], motivated by Dirac-Born-Infeld type conformal scalar theory and designed to yield non-eternal inflation spontaneously. It contains a foliation structure provided by an everywhere time-like vector field n, which couples to the gravitational sector of the theory, but decouples in the small curvature limit. We show that any solution of Einstein gravity with a particular curvature property is a solution of n-DBI gravity. Amongst them is a class of geometries isometric to Reissner-Nordstrom-(Anti) de Sitter black hole, which is obtained within the spherically symmetric solutions of n-DBI gravity minimally coupled to the Maxwell field. These solutions have, however, two distinct features from their Einstein gravity counterparts: 1) the cosmological constant appears as an integration constant and can be positive, negative or vanishing, making it a variable quantity of the theory; 2) there is a non-uniqueness of solutions with the same total mass, charge and effective cosmological constant. Such inequivalent solutions cannot be mapped to each other by a foliation preserving diffeomorphism. Physically they are distinguished by the expansion and shear of the congruence tangent to n, which define scalar invariants on each leave of the foliation.

preprint2011arXiv

Radiation from a D-dimensional collision of shock waves: first order perturbation theory

We study the spacetime obtained by superimposing two equal Aichelburg-Sexl shock waves in D dimensions traveling, head-on, in opposite directions. Considering the collision in a boosted frame, one shock becomes stronger than the other, and a perturbative framework to compute the metric in the future of the collision is setup. The geometry is given, in first order perturbation theory, as an integral solution, in terms of initial data on the null surface where the strong shock has support. We then extract the radiation emitted in the collision by using a D-dimensional generalisation of the Landau-Lifschitz pseudo-tensor and compute the percentage of the initial centre of mass energy epsilon emitted as gravitational waves. In D=4 we find epsilon=25.0%, in agreement with the result of D'Eath and Payne. As D increases, this percentage increases monotonically, reaching 40.0% in D=10. Our result is always within the bound obtained from apparent horizons by Penrose, in D=4, yielding 29.3%, and Eardley and Giddings, in D> 4, which also increases monotonically with dimension, reaching 41.2% in D=10. We also present the wave forms and provide a physical interpretation for the observed peaks, in terms of the null generators of the shocks.

preprint2010arXiv

Black holes in a box: towards the numerical evolution of black holes in AdS

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.

preprint2010arXiv

Numerical relativity for D dimensional axially symmetric space-times: formalism and code tests

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.

preprint2010arXiv

Numerical relativity for D dimensional space-times: head-on collisions of black holes and gravitational wave extraction

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.

preprint2010arXiv

On the Bekenstein-Hawking area law for black objects with conical singularities

We argue that, when working with the appropriate set of thermodynamical variables, the Bekenstein-Hawking law still holds for asymptotically flat black objects with conical singularities. The mass-energy which enters the first law of thermodynamics does not, however, coincide with the ADM mass; it differs from the latter by the energy associated with the conical singularity, as seen by an asymptotic, static observer. These statements are supported by a number of examples: the Bach-Weyl (double-Schwarzschild) solution, its dihole generalisation in Einstein-Maxwell theory and the five dimensional static black ring.

preprint2010arXiv

Thermodynamical description of stationary, asymptotically flat solutions with conical singularities

We examine the thermodynamical properties of a number of asymptotically flat, stationary (but not static) solutions having conical singularities, with both connected and non-connected event horizons, using the thermodynamical description recently proposed in arXiv:0912.3386 [gr-qc]. The examples considered are the double-Kerr solution, the black ring rotating in either S^2 or S^1 and the black Saturn, where the balance condition is not imposed for the latter two solutions. We show that not only the Bekenstein-Hawking area law is recovered from the thermodynamical description but also the thermodynamical angular momentum is the ADM angular momentum. We also analyse the thermodynamical stability and show that, for all these solutions, either the isothermal moment of inertia or the specific heat at constant angular momentum is negative, at any point in parameter space. Therefore, all these solutions are thermodynamically unstable in the grand canonical ensemble.

preprint2001arXiv

String Theory and Hybrid Inflation/Acceleration

We find a description of hybrid inflation in (3+1)-dimensions using brane dynamics of Hanany-Witten type. P-term inflation/acceleration of the universe with the hybrid potential has a slow-roll de Sitter stage and a waterfall stage which leads towards an N=2 supersymmetric ground state. We identify the slow-roll stage of inflation with a non-supersymmetric `Coulomb phase' with Fayet-Iliopoulos term. This stage ends when the mass squared of one of the scalars in the hypermultiplet becomes negative. At that moment the brane system starts undergoing a phase transition via tachyon condensation to a fully Higgsed supersymmetric vacuum which is the absolute ground state of P-term inflation. A string theory/cosmology dictionary is provided, which leads to constraints on parameters of the brane construction from cosmological experiments. We display a splitting of mass levels reminiscent of the Zeeman effect due to spontaneous supersymmetry breaking.