Miguel Zilhão

Testing general relativity with present and future astrophysical observations

One century after its formulation, Einsteins general relativity (GR) 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 GR 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 GR. We present a (necessarily incomplete) catalog of modified theories of gravity for which strong-field predictions have been computed and contrasted to Einsteins 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.

Thermodynamics, transport and relaxation in non-conformal theories

A bstractWe study the equilibrium and near-equilibrium properties of a holographic five-dimensional model consisting of Einstein gravity coupled to a scalar field with a non-trivial potential. The dual four-dimensional gauge theory is not conformal and, at zero temperature, exhibits a renormalisation group flow between two different fixed points. We quantify the deviations from conformality both in terms of thermodynamic observables and in terms of the bulk viscosity of the theory. The ratio of bulk over shear viscosity violates Buchel’s bound. We study relaxation of small-amplitude, homogeneous perturbations by computing the quasi-normal modes of the system at zero spatial momentum. In this approximation we identify two different relaxation channels. At high temperatures, the different pressures first become approximately equal to one another, and subsequently this average pressure evolves towards the equilibrium value dictated by the equation of state. At low temperatures, the average pressure first evolves towards the equilibrium pressure, and only later the different pressures become approximately equal to one another.

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≤a/M≤0.994 and 0≤Q/M≤0.996, corresponding to an extremality parameter a/a_max (amax≡ √(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_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.

Holographic collisions in non-conformal theories

A bstractWe numerically simulate gravitational shock wave collisions in a holographic model dual to a non-conformal four-dimensional gauge theory. We find two novel effects associated to the non-zero bulk viscosity of the resulting plasma. First, the hydrodynamization time increases. Second, if the bulk viscosity is large enough then the plasma becomes well described by hydrodynamics before the energy density and the average pressure begin to obey the equilibrium equation of state. We discuss implications for the quark-gluon plasma created in heavy ion collision experiments.

An Introduction to the Einstein Toolkit

We give an introduction to the Einstein Toolkit, a mature, open-source computational infrastructure for numerical relativity based on the Cactus Framework, for the target group of new users. This toolkit is composed of several different modules, is developed by researchers from different institutions throughout the world and is in active continuous development. Documentation for the toolkit and its several modules is often scattered across different locations, a difficulty new users may at times have to struggle with. Scientific papers exist describing the toolkit and its methods in detail, but they might be overwhelming at first. With these lecture notes we hope to provide an initial overview for new users. We cover how to obtain, compile and run the toolkit, and give an overview of some of the tools and modules provided with it.

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_rad/M=(9.0±0.5)×10^(−4) in five dimensions and E_rad/M=(8.1±0.4)×10^(−4) in six dimensions.

Phase Transitions, Inhomogeneous Horizons and Second-Order Hydrodynamics

A bstractWe use holography to study the spinodal instability of a four-dimensional, strongly-coupled gauge theory with a first-order thermal phase transition. We place the theory on a cylinder in a set of homogeneous, unstable initial states. The dual gravity configurations are black branes afflicted by a Gregory-Laflamme instability. We numerically evolve Einstein’s equations to follow the instability until the system settles down to a stationary, inhomogeneous black brane. The dual gauge theory states have constant temperature but non-constant energy density. We show that the time evolution of the instability and the final states are accurately described by second-order hydrodynamics. In the static limit, the latter reduces to a single, second-order, non-linear differential equation from which the inhomogeneous final states can be derived.

Inspiraling black-hole binary spacetimes: Challenges in transitioning from analytical to numerical techniques

We explore how a recently developed analytical black-hole binary spacetime can be extended using numerical simulations to go beyond the slow-inspiral phase. The analytic spacetime solves the Einstein field equations approximately, with the approximation error becoming progressively smaller the more separated the binary. To continue the spacetime beyond the slow-inspiral phase, we need to transition. Such a transition was previously explored at smaller separations. Here, we perform this transition at a separation of D=20M (large enough that the analytical metric is expected to be accurate), and evolve for six orbits. We find that small constraint violations can have large dynamical effects, but these can be removed by using a constraint-damping system like the conformal covariant formulation of the Z4 system. We find agreement between the subsequent numerical spacetime and the predictions of post-Newtonian theory for the waveform and inspiral rate that is within the post-Newtonian truncation error.

Black Hole Collisions in Asymptotically de Sitter Spacetimes

We report on the first dynamical evolutions of black holes in asymptotically de Sitter spacetimes. We focus on the head-on collision of equal mass binaries and 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.