Numerical Hydrodynamics and Magnetohydrodynamics in General Relativity

Neutron Star Mergers and the Quark Matter Equation of State

As neutron stars merge they can approach very high nuclear density. Here, we summarized recent results for the evolution and gravitational wave emission from binary-neutron star mergers using a a variety of nuclear equations of state with and without a crossover transition to quark matter. We discuss how the late time gravitational wave emission from binary neutron star mergers may possibly reveal the existence of a crossover transition to quark matter.

Machine Learning for Conservative-to-Primitive in Relativistic Hydrodynamics

The numerical solution of relativistic hydrodynamics equations in conservative form requires root-finding algorithms that invert the conservative-to-primitive variables map. These algorithms employ the equation of state of the fluid and can be computationally demanding for applications involving sophisticated microphysics models, such as those required to calculate accurate gravitational wave signals in numerical relativity simulations of binary neutron stars. This work explores the use of machine learning methods to speed up the recovery of primitives in relativistic hydrodynamics. Artificial neural networks are trained to replace either the interpolations of a tabulated equation of state or directly the conservative-to-primitive map. The application of these neural networks to simple benchmark problems shows that both approaches improve over traditional root finders with tabular equation-of-state and multi-dimensional interpolations. In particular, the neural networks for the conservative-to-primitive map accelerate the variable recovery by more than an order of magnitude over standard methods while maintaining accuracy. Neural networks are thus an interesting option to improve the speed and robustness of relativistic hydrodynamics algorithms.

Neutron star merger remnants

Binary neutron star mergers observations are a unique way to constrain fundamental physics and astrophysics at the extreme. The interpretation of gravitational-wave events and their electromagnetic counterparts crucially relies on general-relativistic models of the merger remnants. Quantitative models can be obtained only by means of numerical relativity simulations in 3 + 1 dimensions including detailed input physics for the nuclear matter, electromagnetic and weak interactions. This review summarizes the current understanding of merger remnants focusing on some of the aspects that are relevant for multimessenger observations.

Magnetohydrodynamics as Superfluidity

We show that relativistic magnetohydrodynamics (MHD) can be recast as a novel theory of superfluidity. This new theory formulates MHD just in terms of conservation equations, including dissipative effects, by introducing appropriate variables such as a magnetic scalar potential, and providing necessary and sufficient conditions to obtain equilibrium configurations. We show that this scalar potential can be interpreted as a Goldstone mode originating from the spontaneous breaking of a one-form symmetry, and present the most generic constitutive relations at one derivative order for a parity-preserving plasma in this new superfluid formulation.

Numerical Simulations of Jets from Active Galactic Nuclei

Numerical simulations have been playing a crucial role in the understanding of jets from active galactic nuclei (AGN) since the advent of the first theoretical models for the inflation of giant double radio galaxies by continuous injection in the late 1970s. In the almost four decades of numerical jet research, the complexity and physical detail of simulations, based mainly on a hydrodynamical/magneto-hydrodynamical description of the jet plasma, have been increasing with the pace of the advance in theoretical models, computational tools and numerical methods. The present review summarizes the status of the numerical simulations of jets from AGNs, from the formation region in the neighborhood of the supermassive central black hole up to the impact point well beyond the galactic scales. Special attention is paid to discuss the achievements of present simulations in interpreting the phenomenology of jets as well as their current limitations and challenges.

Mapping nonlinear gravity into General Relativity with nonlinear electrodynamics

We show that families of nonlinear gravity theories formulated in a metric-affine approach and coupled to a nonlinear theory of electrodynamics can be mapped into general relativity (GR) coupled to another nonlinear theory of electrodynamics. This allows to generate solutions of the former from those of the latter using purely algebraic transformations. This correspondence is explicitly illustrated with the Eddington-inspired Born-Infeld theory of gravity, for which we consider a family of nonlinear electrodynamics and show that, under the map, preserve their algebraic structure. For the particular case of Maxwell electrodynamics coupled to Born-Infeld gravity we find, via this correspondence, a Born-Infeld-type nonlinear electrodynamics on the GR side. Solving the spherically symmetric electrovacuum case for the latter, we show how the map provides directly the right solutions for the former. This procedure opens a new door to explore astrophysical and cosmological scenarios in nonlinear gravity theories by exploiting the full power of the analytical and numerical methods developed within the framework of GR.

Rotating stars in relativity

Rotating relativistic stars have been studied extensively in recent years, both theoretically and observationally, because of the information they might yield about the equation of state of matter at extremely high densities and because they are considered to be promising sources of gravitational waves. The latest theoretical understanding of rotating stars in relativity is reviewed in this updated article. The sections on equilibrium properties and on nonaxisymmetric oscillations and instabilities in f-modes and r-modes have been updated. Several new sections have been added on equilibria in modified theories of gravity, approximate universal relationships, the one-arm spiral instability, on analytic solutions for the exterior spacetime, rotating stars in LMXBs, rotating strange stars, and on rotating stars in numerical relativity including both hydrodynamic and magnetohydrodynamic studies of these objects.

Dynamical boson stars

The idea of stable, localized bundles of energy has strong appeal as a model for particles. In the 1950s, John Wheeler envisioned such bundles as smooth configurations of electromagnetic energy that he called geons, but none were found. Instead, particle-like solutions were found in the late 1960s with the addition of a scalar field, and these were given the name boson stars. Since then, boson stars find use in a wide variety of models as sources of dark matter, as black hole mimickers, in simple models of binary systems, and as a tool in finding black holes in higher dimensions with only a single Killing vector. We discuss important varieties of boson stars, their dynamic properties, and some of their uses, concentrating on recent efforts.

Effects of Nonlinear Inhomogeneity on the Cosmic Expansion with Numerical Relativity

We construct a three-dimensional, fully relativistic numerical model of a universe filled with an inhomogeneous pressureless fluid, starting from initial data that represent a perturbation of the Einstein-de Sitter model. We then measure the departure of the average expansion rate with respect to this homogeneous and isotropic reference model, comparing local quantities to the predictions of linear perturbation theory. We find that collapsing perturbations reach the turnaround point much earlier than expected from the reference spherical top-hat collapse model and that the local deviation of the expansion rate from the homogeneous one can be as high as 28% at an underdensity, for an initial density contrast of 10^{-2}. We then study, for the first time, the exact behavior of the backreaction term Q_{D}. We find that, for small values of the initial perturbations, this term exhibits a 1/a scaling, and that it is negative with a linearly growing absolute value for larger perturbation amplitudes, thereby contributing to an overall deceleration of the expansion. Its magnitude, on the other hand, remains very small even for relatively large perturbations.

Modeling the Complete Gravitational Wave Spectrum of Neutron Star Mergers

In the context of neutron star mergers, we study the gravitational wave spectrum of the merger remnant using numerical relativity simulations. Postmerger spectra are characterized by a main peak frequency f2 related to the particular structure and dynamics of the remnant hot hypermassive neutron star. We show that f(2) is correlated with the tidal coupling constant κ(2)^T that characterizes the binary tidal interactions during the late-inspiral merger. The relation f(2)(κ(2)^T) depends very weakly on the binary total mass, mass ratio, equation of state, and thermal effects. This observation opens up the possibility of developing a model of the gravitational spectrum of every merger unifying the late-inspiral and postmerger descriptions.

Testing General Relativity with Low-Frequency, Space-Based Gravitational-Wave Detectors

We review the tests of general relativity that will become possible with space-based gravitational-wave detectors operating in the ∼ 10-5 - 1 Hz low-frequency band. The fundamental aspects of gravitation that can be tested include the presence of additional gravitational fields other than the metric; the number and tensorial nature of gravitational-wave polarization states; the velocity of propagation of gravitational waves; the binding energy and gravitational-wave radiation of binaries, and therefore the time evolution of binary inspirals; the strength and shape of the waves emitted from binary mergers and ringdowns; the true nature of astrophysical black holes; and much more. The strength of this science alone calls for the swift implementation of a space-based detector; the remarkable richness of astrophysics, astronomy, and cosmology in the low-frequency gravitational-wave band make the case even stronger.

Foundations of Black Hole Accretion Disk Theory

This review covers the main aspects of black hole accretion disk theory. We begin with the view that one of the main goals of the theory is to better understand the nature of black holes themselves. In this light we discuss how accretion disks might reveal some of the unique signatures of strong gravity: the event horizon, the innermost stable circular orbit, and the ergosphere. We then review, from a first-principles perspective, the physical processes at play in accretion disks. This leads us to the four primary accretion disk models that we review: Polish doughnuts (thick disks), Shakura-Sunyaev (thin) disks, slim disks, and advection-dominated accretion flows (ADAFs). After presenting the models we discuss issues of stability, oscillations, and jets. Following our review of the analytic work, we take a parallel approach in reviewing numerical studies of black hole accretion disks. We finish with a few select applications that highlight particular astrophysical applications: measurements of black hole mass and spin, black hole vs. neutron star accretion disks, black hole accretion disk spectral states, and quasi-periodic oscillations (QPOs).

Continuum and Discrete Initial-Boundary Value Problems and Einstein's Field Equations

Many evolution problems in physics are described by partial differential equations on an infinite domain; therefore, one is interested in the solutions to such problems for a given initial dataset. A prominent example is the binary black-hole problem within Einstein's theory of gravitation, in which one computes the gravitational radiation emitted from the inspiral of the two black holes, merger and ringdown. Powerful mathematical tools can be used to establish qualitative statements about the solutions, such as their existence, uniqueness, continuous dependence on the initial data, or their asymptotic behavior over large time scales. However, one is often interested in computing the solution itself, and unless the partial differential equation is very simple, or the initial data possesses a high degree of symmetry, this computation requires approximation by numerical discretization. When solving such discrete problems on a machine, one is faced with a finite limit to computational resources, which leads to the replacement of the infinite continuum domain with a finite computer grid. This, in turn, leads to a discrete initial-boundary value problem. The hope is to recover, with high accuracy, the exact solution in the limit where the grid spacing converges to zero with the boundary being pushed to infinity. The goal of this article is to review some of the theory necessary to understand the continuum and discrete initial boundary-value problems arising from hyperbolic partial differential equations and to discuss its applications to numerical relativity; in particular, we present well-posed initial and initial-boundary value formulations of Einstein's equations, and we discuss multi-domain high-order finite difference and spectral methods to solve them.

Binary Neutron Star Mergers

We review the current status of studies of the coalescence of binary neutron star systems. We begin with a discussion of the formation channels of merging binaries and we discuss the most recent theoretical predictions for merger rates. Next, we turn to the quasi-equilibrium formalisms that are used to study binaries prior to the merger phase and to generate initial data for fully dynamical simulations. The quasi-equilibrium approximation has played a key role in developing our understanding of the physics of binary coalescence and, in particular, of the orbital instability processes that can drive binaries to merger at the end of their lifetimes. We then turn to the numerical techniques used in dynamical simulations, including relativistic formalisms, (magneto-)hydrodynamics, gravitational-wave extraction techniques, and nuclear microphysics treatments. This is followed by a summary of the simulations performed across the field to date, including the most recent results from both fully relativistic and microphysically detailed simulations. Finally, we discuss the likely directions for the field as we transition from the first to the second generation of gravitational-wave interferometers and while supercomputers reach the petascale frontier.

Dynamical Boson Stars

The idea of stable, localized bundles of energy has strong appeal as a model for particles. In the 1950s, John Wheeler envisioned such bundles as smooth configurations of electromagnetic energy that he called geons, but none were found. Instead, particle-like solutions were found in the late 1960s with the addition of a scalar field, and these were given the name boson stars. Since then, boson stars find use in a wide variety of models as sources of dark matter, as black hole mimickers, in simple models of binary systems, and as a tool in finding black holes in higher dimensions with only a single Killing vector. We discuss important varieties of boson stars, their dynamic properties, and some of their uses, concentrating on recent efforts.

Characteristic Evolution and Matching

I review the development of numerical evolution codes for general relativity based upon the characteristic initial-value problem. Progress in characteristic evolution is traced from the early stage of 1D feasibility studies to 2D-axisymmetric codes that accurately simulate the oscillations and gravitational collapse of relativistic stars and to current 3D codes that provide pieces of a binary black-hole spacetime. Cauchy codes have now been successful at simulating all aspects of the binary black-hole problem inside an artificially constructed outer boundary. A prime application of characteristic evolution is to extend such simulations to null infinity where the waveform from the binary inspiral and merger can be unambiguously computed. This has now been accomplished by Cauchy-characteristic extraction, where data for the characteristic evolution is supplied by Cauchy data on an extraction worldtube inside the artificial outer boundary. The ultimate application of characteristic evolution is to eliminate the role of this outer boundary by constructing a global solution via Cauchy-characteristic matching. Progress in this direction is discussed.

Coalescence of Black Hole-Neutron Star Binaries

We review the current status of general relativistic studies for the coalescence of black hole-neutron star (BH-NS) binaries. First, procedures for a solution of BH-NS binaries in quasi-equilibrium circular orbits and the numerical results, such as quasi-equilibrium sequence and mass-shedding limit, of the high-precision computation, are summarized. Then, the current status of numerical-relativity simulations for the merger of BH-NS binaries is described. We summarize our understanding for the merger and/or tidal disruption processes, the criterion for tidal disruption, the properties of the remnant formed after the tidal disruption, gravitational waveform, and gravitational-wave spectrum.