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.

We model the large scale late time universe as a Lambda-CDM cosmology driven by cosmological constant and perfect dust fluid. Our aim is to find new solutions in the matter and Lambda epoch consistent with inflationary initial conditions, namely that to the far past in the matter era the cosmology tends to a flat FLRW solution. We identify the moduli degrees of freedom that parametrize the flat Lambda-dust FLRW solution and then promote these moduli to slowly varying functions of the spatial coordinates and show how to solve the Einstein equations in a comoving gradient expansion, controlled by the cosmological constant length scale. Our initial conditions ensure that the approximation remains under control to the far past of the matter era, and to the far future of Lambda domination. The solution is fully non-perturbative in the amplitude of the metric deformation, and we explicitly construct it to fourth order in derivatives. A general Lambda-dust universe dominated by Lambda in the future is characterized by a 3-metric and a stress tensor (with positive trace) defined on the future conformal boundary. The new cosmologies with inflationary initial conditions are characterized only by the boundary 3-metric, the stress tensor being locally determined entirely in terms of that metric.

Assuming homogeneous isotropic Lambda-CDM cosmology allows Lambda, spatial curvature and dark matter density to be inferred from large scale structure observations such as supernovae. The purpose of this paper is to extend this to allow observations to measure or constrain inhomogeneity and anisotropy. We obtain the general inhomogeneous anisotropic Lambda-CDM solution which is locally asymptotic to an expanding de Sitter solution as a late time expansion using Starobinsky's method (analogous to the `holographic renormalization' technique in AdS/CFT) together with a resummation of the series. The dark matter is modeled as perfect dust fluid. The terms in the expansion systematically describe inhomogeneous and anisotropic deformations of an expanding FLRW solution, and are given as a spatial derivative expansion in terms of data characterizing the solution - a 3-metric and a perturbation of that 3-metric. Leading terms describe inhomogeneity and anisotropy on the scale set by the cosmological constant, approximately the horizon scale today. Higher terms in the expansion describe shorter scale variations. We compute the luminosity distance-redshift relation and argue that comparison with current and future observation would allow a partial reconstruction of the characterizing data. We also comment on smoothing these solutions noting that geometric flows (such as Ricci flow) applied to the characterizing data provide a canonical averaging method.

TeVeS uses a dynamical vector field with timelike unit norm constraint to specify a preferred local frame. When matter moves slowly in this frame - the so-called quasi-static regime - Modified Newtonian Dynamics (MoND) results. Theories with such vectors (such as Einstein-aether) are prone to the vector dynamics forming singularities which render their classical evolution problematic. Here we analyse the dynamics of the vector in TeVeS in various situations. We find that, quite generically, the vector field develops caustic singularities on time scales of order the gravitational in-fall time. Having shown singularity formation is generic with or without matter, Bekenstein's original formulation of TeVeS appears dynamically problematic. We argue that by modifying the vector field kinetic terms to the more general form used by Einstein-Aether this problem may be avoided.

We investigate the universal low-energy dynamics of the simplest Higgs phase for gravity, `ghost condensation.' We show that the nonlinear dynamics of the `ghostone' field dominate for all interesting gravitational sources. Away from caustic singularities, the dynamics is equivalent to the irrotational flow of a perfect fluid with equation of state p \propto \rho^2, where the fluid particles can have negative mass. We argue that this theory is free from catastrophic instabilities due to growing modes, even though the null energy condition is violated. Numerical simulations show that solutions generally have singularities in which negative energy regions shrink to zero size. We exhibit partial UV completions of the theory in which these singularities are smoothly resolved, so this does not signal any inconsistency in the effective theory. We also consider the bounds on the symmetry breaking scale M in this theory. We argue that the nonlinear dynamics cuts off the Jeans instability of the linear theory, and allows M \lsim 100MeV.

The non-linear behaviour of Randall-Sundrum gravity with one brane is examined. Due to the non-compact extra dimension, the perturbation spectrum has no mass gap, and the long wavelength effective theory is only understood perturbatively. The full 5-dimensional Einstein equations are solved numerically for static, spherically symmetric matter localized on the brane, yielding regular geometries in the bulk with axial symmetry. An elliptic relaxation method is used, allowing both the brane and asymptotic radiation boundary conditions to be simultaneously imposed. The same data that specifies stars in 4-dimensional gravity, uniquely constructs a 5-dimensional solution. The algorithm performs best for small stars (radius less than the AdS length) yielding highly non-linear solutions. An upper mass limit is observed for these small stars, and the geometry shows no global pathologies. The geometric perturbation is shown to remain localized near the brane at high densities, the confinement interestingly increasing for both small and large stars as the upper mass limit is approached. Furthermore, the static spatial sections are found to be approximately conformal to those of AdS. We show that the intrinsic geometry of large stars, with radius several times the AdS length, is described by 4-dimensional General Relativity far past the perturbative regime. This indicates that the non-linear long wavelength effective action remains local, even though the perturbation spectrum has no mass gap. The implication is that Randall-Sundrum gravity, with localized brane matter, reproduces relativistic astrophysical solutions, such as neutron stars and massive black holes, consistent with observation.

For the 2-brane Randall-Sundrum model, we calculate the bulk geometry for strong gravity, in the low matter density regime, for slowly varying matter sources. This is relevant for astrophysical or cosmological applications. The warped compactification means the radion can not be written as a homogeneous mode in the orbifold coordinate, and we introduce it by extending the coordinate patch approach of the linear theory to the non-linear case. The negative tension brane is taken to be in vacuum. For conformally invariant matter on the positive tension brane, we solve the bulk geometry as a derivative expansion, formally summing the `Kaluza-Klein' contributions to all orders. For general matter we compute the Einstein equations to leading order, finding a scalar-tensor theory with $\omega(\Psi) \propto \Psi / (1 - \Psi)$, and geometrically interpret the radion. We comment that this radion scalar may become large in the context of strong gravity with low density matter. Equations of state allowing $(\rho - 3 P)$ to be negative, can exhibit behavior where the matter decreases the distance between the 2 branes, which we illustrate numerically for static star solutions using an incompressible fluid. For increasing stellar density, the branes become close before the upper mass limit, but after violation of the dominant energy condition. This raises the interesting question of whether astrophysically reasonable matter, and initial data, could cause branes to collide at low energy, such as in dynamical collapse.