Einstein originally proposed a nonsymmetric tensor field, with its symmetric part associated with the spacetime metric and its antisymmetric part associated with the electromagnetic field, as an approach to a unified field theory. Here we interpret it more modestly as an alternative to Einstein-Maxwell theory, approximating the coupling between the electromagnetic field and spacetime curvature in the macroscopic classical regime. Previously it was shown that the Lorentz force can be derived from this theory, albeit with deviation on the scale of a universal length constant $\ell$. Here we assume that $\ell$ is of galactic scale and show that the modified coupling of the electromagnetic field with charged particles allows a non-Maxwellian equilibrium of non-neutral plasma. The resulting electromagnetic field is "dark" in the sense that its modified Lorentz force on the plasma vanishes, yet through its modified coupling to the gravitational field it engenders a nonvanishing, effective mass density. We obtain a solution for which this mass density asymptotes approximately to that of the pseudo-isothermal model of dark matter. The resulting gravitational field produces radial acceleration, in the context of a post-Minkowskian approximation, which is negligible at small radius but yields a flat rotation curve at large radius. We further exhibit a family of such solutions which, like the pseudo-isothermal model, has a free parameter to set the mass scale (in this case related to the charge density) and a free parameter to set the length scale (in this case an integer multiple of $\ell$). Moreover, these solutions are members of a larger family with more general angular and radial dependence. They thus show promise as approximations of generalized pseudo-isothermal models, which in turn are known to fit a wide range of mass density profiles for galaxies and clusters.

Coalescing supermassive black hole binaries are produced by the mergers of galaxies and are the most powerful sources of gravitational waves accessible to space-based gravitational observatories. Some such mergers may occur in the presence of matter and magnetic fields and hence generate an electromagnetic counterpart. In this Letter, we present the first general relativistic simulations of magnetized plasma around merging supermassive black holes using the general relativistic magnetohydrodynamic code Whisky. By considering different magnetic field strengths, going from non-magnetically dominated to magnetically dominated regimes, we explore how magnetic fields affect the dynamics of the plasma and the possible emission of electromagnetic signals. In particular we observe a total amplification of the magnetic field of ~2 orders of magnitude which is driven by the accretion onto the binary and that leads to much stronger electromagnetic signals, more than a factor of 10^4 larger than comparable calculations done in the force-free regime where such amplifications are not possible.

Recent breakthroughs in the field of numerical relativity have led to dramatic progress in understanding the predictions of General Relativity for the dynamical interactions of two black holes in the regime of very strong gravitational fields. Such black-hole binaries are important astrophysical systems and are a key target of current and developing gravitational-wave detectors. The waveform signature of strong gravitational radiation emitted as the black holes fall together and merge provides a clear observable record of the process. After decades of slow progress, these mergers and the gravitational-wave signals they generate can now be routinely calculated using the methods of numerical relativity. We review recent advances in understanding the predicted physics of events and the consequent radiation, and discuss some of the impacts this new knowledge is having in various areas of astrophysics.

Understanding the predictions of general relativity for the dynamical interactions of two black holes has been a long-standing unsolved problem in theoretical physics. Black-hole mergers are monumental astrophysical events, releasing tremendous amounts of energy in the form of gravitational radiation, and are key sources for both ground- and space-based gravitational-wave detectors. The black-hole merger dynamics and the resulting gravitational waveforms can only be calculated through numerical simulations of Einstein's equations of general relativity. For many years, numerical relativists attempting to model these mergers encountered a host of problems, causing their codes to crash after just a fraction of a binary orbit could be simulated. Recently, however, a series of dramatic advances in numerical relativity has allowed stable, robust black-hole merger simulations. This remarkable progress in the rapidly maturing field of numerical relativity, and the new understanding of black-hole binary dynamics that is emerging is chronicled. Important applications of these fundamental physics results to astrophysics, to gravitational-wave astronomy, and in other areas are also discussed.

Although the gravitational wave kick velocity in the orbital plane of coalescing black holes has been understood for some time, apparently conflicting formulae have been proposed for the dominant out-of-plane kick, each a good fit to different data sets. This is important to resolve because it is only the out-of-plane kicks that can reach more than 500 km/s and can thus eject merged remnants from galaxies. Using a different ansatz for the out-of-plane kick, we show that we can fit almost all existing data to better than 5 %. This is good enough for any astrophysical calculation, and shows that the previous apparent conflict was only because the two data sets explored different aspects of the kick parameter space.

Coalescing massive black hole binaries are produced by the mergers of galaxies. The final stages of the black hole coalescence produce strong gravitational radiation that can be detected by the space-borne LISA. In cases where the black hole merger takes place in the presence of gas and magnetic fields, various types of electromagnetic signals may also be produced. Modeling such electromagnetic counterparts of the final merger requires evolving the behavior of both gas and fields in the strong-field regions around the black holes. We have taken a step towards solving this problem by mapping the flow of pressureless matter in the dynamic, 3-D general relativistic spacetime around the merging black holes. We find qualitative differences in collision and outflow speeds, including a signature of the merger when the net angular momentum of the matter is low, between the results from single and binary black holes, and between nonrotating and rotating holes in binaries. If future magnetohydrodynamic results confirm these differences, it may allow assessment of the properties of the binaries as well as yielding an identifiable electromagnetic counterpart to the attendant gravitational wave signal.

Recent numerical relativistic results demonstrate that the merger of comparable-mass spinning black holes has a maximum ``recoil kick'' of up to $\sim 4000 \kms$. However the scaling of these recoil velocities with mass ratio is poorly understood. We present new runs showing that the maximum possible kick perpendicular to the orbital plane does not scale as $\sim\eta^2$ (where $\eta$ is the symmetric mass ratio), as previously proposed, but is more consistent with $\sim\eta^3$, at least for systems with low orbital precession. We discuss the effect of this dependence on galactic ejection scenarios and retention of intermediate-mass black holes in globular clusters.

We present a multipolar analysis of the gravitational recoil computed in recent numerical simulations of binary black hole (BH) coalescence, for both unequal masses and non-zero, non-precessing spins. We show that multipole moments up to and including l=4 are sufficient to accurately reproduce the final recoil velocity (within ~2%) and that only a few dominant modes contribute significantly to it (within ~5%). We describe how the relative amplitudes, and more importantly, the relative phases, of these few modes control the way in which the recoil builds up throughout the inspiral, merger, and ringdown phases. We also find that the numerical results can be reproduced by an ``effective Newtonian'' formula for the multipole moments obtained by replacing the radial separation in the Newtonian formulae with an effective radius computed from the numerical data. Beyond the merger, the numerical results are reproduced by a superposition of three Kerr quasi-normal modes (QNMs). Analytic formulae, obtained by expressing the multipole moments in terms of the fundamental QNMs of a Kerr BH, are able to explain the onset and amount of ``anti-kick'' for each of the simulations. Lastly, we apply this multipolar analysis to help explain the remarkable difference between the amplitudes of planar and non-planar kicks for equal-mass spinning black holes.

Several groups have recently computed the gravitational radiation recoil produced by the merger of two spinning black holes. The results suggest that spin can be the dominant contributor to the kick, with reported recoil speeds of hundreds to even thousands of kilometers per second. The parameter space of spin kicks is large, however, and it is ultimately desirable to have a simple formula that gives the approximate magnitude of the kick given a mass ratio, spin magnitudes, and spin orientations. As a step toward this goal, we perform a systematic study of the recoil speeds from mergers of black holes with mass ratio $q\equiv m_1/m_2=2/3$ and dimensionless spin parameters of $a_1/m_1$ and $a_2/m_2$ equal to 0 or 0.2, with directions aligned or anti-aligned with the orbital angular momentum. We also run an equal-mass $a_1/m_1=-a_2/m_2=0.2$ case, and find good agreement with previous results. We find that, for currently reported kicks from aligned or anti-aligned spins, a simple kick formula inspired by post-Newtonian analyses can reproduce the numerical results to better than $\sim$10%.

Coalescing black-hole binaries are expected to be the strongest sources of gravitational waves for ground-based interferometers as well as the space-based interferometer LISA. Recent progress in numerical relativity now makes it possible to calculate the waveforms from the strong-field dynamical merger and is revolutionizing our understanding of these systems. We review these dramatic developments, emphasizing applications to issues in gravitational wave observations. These new capabilities also make possible accurate calculations of the recoil or kick imparted to the final remnant black hole when the merging components have unequal masses, or unequal or unaligned spins. We highlight recent work in this area, focusing on results of interest to astrophysics.

Recent developments in numerical relativity have made it possible to follow reliably the coalescence of two black holes from near the innermost stable circular orbit to final ringdown. This opens up a wide variety of exciting astrophysical applications of these simulations. Chief among these is the net kick received when two unequal mass or spinning black holes merge. The magnitude of this kick has bearing on the production and growth of supermassive black holes during the epoch of structure formation, and on the retention of black holes in stellar clusters. Here we report the first accurate numerical calculation of this kick, for two nonspinning black holes in a 1.5:1 mass ratio, which is expected based on analytic considerations to give a significant fraction of the maximum possible recoil. We have performed multiple runs with different initial separations, orbital angular momenta, resolutions, extraction radii, and gauges. The full range of our kick speeds is 86--116 km s$^{-1}$, and the most reliable runs give kicks between 86 and 97 km s$^{-1}$. This is intermediate between the estimates from two recent post-Newtonian analyses and suggests that at redshifts $z\gtrsim 10$, halos with masses $\lesssim 10^9 M_\odot$ will have difficulty retaining coalesced black holes after major mergers.

We study dynamics and radiation generation in the last few orbits and merger of a binary black hole system, applying recently developed techniques for simulations of moving black holes. Our analysis of the gravitational radiation waveforms and dynamical black hole trajectories produces a consistent picture for a set of simulations with black holes beginning on circular-orbit trajectories at a variety of initial separations. We find profound agreement at the level of one percent among the simulations for the last orbit, merger and ringdown. We are confident that this part of our waveform result accurately represents the predictions from Einstein's General Relativity for the final burst of gravitational radiation resulting from the merger of an astrophysical system of equal-mass non-spinning black holes. The simulations result in a final black hole with spin parameter a/m=0.69. We also find good agreement at a level of roughly 10 percent for the radiation generated in the preceding few orbits.

We present new techniqes for evolving binary black hole systems which allow the accurate determination of gravitational waveforms directly from the wave zone region of the numerical simulations. Rather than excising the black hole interiors, our approach follows the "puncture" treatment of black holes, but utilzing a new gauge condition which allows the black holes to move successfully through the computational domain. We apply these techniques to an inspiraling binary, modeling the radiation generated during the final plunge and ringdown. We demonstrate convergence of the waveforms and good conservation of mass-energy, with just over 3% of the system's mass converted to gravitional radiation.