Critical Phenomena in Gravitational Collapse

by: Carsten Gundlach and José M. Martín-García

As first discovered by Choptuik, the black hole threshold in the space of initial data for general relativity shows both surprising structure and surprising simplicity. Universality, power-law scaling of the black hole mass, and scale echoing have given rise to the term “critical phenomena”. They are explained by the existence of exact solutions which are attractors within the black hole threshold, that is, attractors of codimension one in phase space, and which are typically self-similar. Critical phenomena give a natural route from smooth initial data to arbitrarily large curvatures visible from infinity, and are therefore likely to be relevant for cosmic censorship, quantum gravity, astrophysics, and our general understanding of the dynamics of general relativity.

The Hubble Constant

by: Neal Jackson

I review the current state of determinations of the Hubble constant, which gives the length scale of the Universe by relating the expansion velocity of objects to their distance. In the last 20 years, much progress has been made and estimates now range between 60 and 75 km s^-1 Mpc^-1, with most now between 70 and 75 km s^-1 Mpc^-1, a huge improvement over the factor-of-2 uncertainty which used to prevail. Further improvements which gave a generally agreed margin of error of a few percent rather than the current 10% would be vital input to much other interesting cosmology. There are several programmes which are likely to lead us to this point in the next 10 years.

Event and Apparent Horizon Finders for 3+1 Numerical Relativity

by: Jonathan Thornburg

Event and apparent horizons are key diagnostics for the presence and properties of black holes. In this article I review numerical algorithms and codes for finding event and apparent horizons in numerically-computed spacetimes, focusing on calculations done using the 3+1 ADM formalism. The event horizon of an asymptotically-flat spacetime is the boundary between those events from which a future-pointing null geodesic can reach future null infinity and those events from which no such geodesic exists. The event horizon is a (continuous) null surface in spacetime. The event horizon is defined nonlocally in time: it is a global property of the entire spacetime and must be found in a separate post-processing phase after all (or at least the nonstationary part) of spacetime has been numerically computed.

There are three basic algorithms for finding event horizons, based on integrating null geodesics forwards in time, integrating null geodesics backwards in time, and integrating null surfaces backwards in time. The last of these is generally the most efficient and accurate.

In contrast to an event horizon, an apparent horizon is defined locally in time in a spacelike slice and depends only on data in that slice, so it can be (and usually is) found during the numerical computation of a spacetime. A marginally outer trapped surface (MOTS) in a slice is a smooth closed 2-surface whose future-pointing outgoing null geodesics have zero expansion Theta. An apparent horizon is then defined as a MOTS not contained in any other MOTS. The MOTS condition is a nonlinear elliptic partial differential equation (PDE) for the surface shape, containing the ADM 3-metric, its spatial derivatives, and the extrinsic curvature as coefficients. Most “apparent horizon” finders actually find MOTSs.

There are a large number of apparent horizon finding algorithms, with differing trade-offs between speed, robustness, accuracy, and ease of programming. In axisymmetry, shooting algorithms work well and are fairly easy to program. In slices with no continuous symmetries, spectral integral-iteration algorithms and elliptic-PDE algorithms are fast and accurate, but require good initial guesses to converge. In many cases, Schnetter’s “pretracking” algorithm can greatly improve an elliptic-PDE algorithm’s robustness. Flow algorithms are generally quite slow but can be very robust in their convergence. Minimization methods are slow and relatively inaccurate in the context of a finite differencing simulation, but in a spectral code they can be relatively faster and more robust.

The Post-Newtonian Approximation for Relativistic Compact Binaries

by: Toshifumi Futamase and Yousuke Itoh

We discuss various aspects of the post-Newtonian approximation in general relativity. After presenting the foundation based on the Newtonian limit, we show a method to derive post-Newtonian equations of motion for relativistic compact binaries based on a surface integral approach and the strong field point particle limit. As an application we derive third post-Newtonian equations of motion for relativistic compact binaries which respect the Lorentz invariance in the post-Newtonian perturbative sense, admit a conserved energy, and are free from any ambiguity.

Relativistic Fluid Dynamics: Physics for Many Different Scales

by: Nils Andersson and Gregory L. Comer

The relativistic fluid is a highly successful model used to describe the dynamics of many-particle, relativistic systems. It takes as input basic physics from microscopic scales and yields as output predictions of bulk, macroscopic motion. By inverting the process, an understanding of bulk features can lead to insight into physics on the microscopic scale. Relativistic fluids have been used to model systems as “small” as heavy ions in collisions, and as large as the Universe itself, with “intermediate” sized objects like neutron stars being considered along the way. The purpose of this review is to discuss the mathematical and theoretical physics underpinnings of the relativistic (multiple) fluid model. We focus on the variational principle approach championed by Brandon Carter and his collaborators, in which a crucial element is to distinguish the momenta that are conjugate to the particle number density currents. This approach differs from the “standard” text-book derivation of the equations of motion from the divergence of the stress-energy tensor in that one explicitly obtains the relativistic Euler equation as an “integrability” condition on the relativistic vorticity. We discuss the conservation laws and the equations of motion in detail, and provide a number of (in our opinion) interesting and relevant applications of the general theory.