Black holes in binaries with other compact objects can provide natural venues for indirect detection of axions or other ultralight fields. The superradiant instability associated with a rapidly spinning black hole leads to the creation of an axion cloud which carries energy and angular momentum from the black hole. This cloud will then decay via gravitational wave emission. We show that the energy lost as a result of this process tends toward an outspiraling of the binary orbit. A given binary system is sensitive to a narrow range of axion masses, determined by the mass of the black hole. Pulsar-black hole binaries, once detected in the electromagnetic band, will allow high-precision measurements of black hole mass loss via timing measurements of the companion pulsar. This avenue of investigation is particularly promising in light of the recent preliminary announcements of two candidate black hole-neutron star mergers by LIGO/VIRGO (#S190814bv and #S190426c). We demonstrate that for such a binary system with a typical millisecond pulsar and a 3-solar-mass black hole, axions with masses between $2.7 \times 10^{-12}$ eV and $3.2 \times 10^{-12}$ eV are detectable. Recent gravitational wave observations by LIGO/VIRGO of binary black hole mergers imply that, for these binaries, gravitational radiation from the rotating quadrupole moment is a dominant effect, causing an inspiraling orbit. With some reasonable assumptions about the period of the binary when it formed and the spins of the black holes, these observations rule out possible axion masses between $3 \times 10^{-13}$ eV and $6 \times 10^{-13}$ eV. Future binary black hole observations, for example by LISA, are expected to provide more robust bounds. In some circumstances, neutron stars may also undergo superradiant instabilities, and binary pulsars could be used to exclude axions with certain masses and matter couplings.

The conventional view of stellar-mass binary black hole (sBBH) mergers is that there should not be enough matter present to produce a detectable electromagnetic transient. However, there ARE a number of mechanisms for producing such a counterpart through accretion of matter from various reserves, charged black holes, or interactions with magnetic or exotic fields. After reviewing these mechanisms, we describe what we can learn from multi-messenger observations of sBBH mergers in the areas of stellar evolution and compact binary formation, cosmological measurements, searches for charged black holes, and tests of general relativity and fundamental physics. We end with some remarks about the need to support both observing capabilities and modeling in the next decade and beyond.

The strong version of the nonviolent nonlocality proposal of Giddings predicts "strong but soft" quantum metric fluctuations near black hole horizons in an attempt to resolve the information paradox. To study observable signatures of this proposal, we numerically solve Einstein's equations modified by these fluctuations and analyze the gravitational wave signal from the inspiral and merger of two black holes. In a model of evolution for such fluctuations, we show that they lead to significant deviations in the observed waveform, even when the black holes are still well separated, and could potentially be observed by aLIGO.

Observations of gravitational radiation from compact binary systems provide an unprecedented opportunity to test General Relativity in the strong field dynamical regime. In this paper, we investigate how future observations of gravitational radiation from binary neutron star mergers might provide constraints on finite-range forces from a universally coupled massive scalar field. Such scalar degrees of freedom are a characteristic feature of many extensions of General Relativity. For concreteness, we work in the context of metric $f(R)$ gravity, which is equivalent to General Relativity and a universally coupled scalar field with a non-linear potential whose form is fixed by the choice of $f(R)$. In theories where neutron stars (or other compact objects) obtain a significant scalar charge, the resulting attractive finite-range scalar force has implications for both the inspiral and merger phases of binary systems. We first present an analysis of the inspiral dynamics in Newtonian limit, and forecast the constraints on the mass of the scalar and charge of the compact objects for the Advanced LIGO gravitational wave observatory. We then perform a comparative study of binary neutron star mergers in General Relativity with those of a one-parameter model of $f(R)$ gravity using fully relativistic hydrodynamical simulations. These simulations elucidate the effects of the scalar on the merger and post-merger dynamics. We comment on the utility of the full waveform (inspiral, merger, post-merger) to probe different regions of parameter space for both the particular model of $f(R)$ gravity studied here and for finite-range scalar forces more generally.

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 {\em 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 {\em 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.

The observation of a possible electromagnetic counterpart by the Fermi GBM group to the aLIGO detection of the merger of a black hole binary has spawned a number of ideas about its source. Furthermore, observations of fast radio bursts (FRBs) have similarly resulted in a range of new models that might endow black hole binaries with electromagnetic signatures. In this context, even the unlikely idea that astrophysical black holes may have significant charge is worth exploring, and here we present results from the simulation of weakly charged black holes as they orbit and merge. Our simulations suggest that a black hole binary with mass comparable to that observed in GW150914 could produce the level of electromagnetic luminosity observed by Fermi GBM ($10^{49}$ ergs/s) with a non-dimensional charge of $q \equiv Q/M = 10^{-4}$ assuming good radiative efficiency. However even a charge such as this is difficult to imagine avoiding neutralization long enough for the binary to produce its electromagnetic counterpart, and so this value would likely serve simply as an upper bound. On the other hand, one can equivalently consider the black holes as having acquired a magnetic monopole charge that would be easy to maintain and would generate an identical electromagnetic signature as the electric charges. The observation of such a binary would have significant cosmological implications, not the least of which would be an explanation for the quantization of charge itself. We also study such a magnetically charged binary in the force-free regime and find it much more radiative, reducing even further the requirements to produce the counterpart.

We examine the development and detectability of the $m=1$ instability in the remnant of binary neutron star mergers. The detection of the gravitational mode associated with the $m=1$ degree of freedom could potentially reveal details of the equation of state. We analyze the post-merger epoch of simulations of both equal and non-equal mass neutron star mergers using three realistic, microphysical equations of state and neutrino cooling. Our studies show such an instability develops generically and within a short dynamical time to strengths that are comparable or stronger than the $m=2$ mode which is the strongest during the early post-merger stage. We estimate the signal to noise ratio that might be obtained for the $m=1$ mode and discuss the prospects for observing this signal with available Earth-based detectors. Because the $m=1$ occurs at roughly half the frequency of the more powerful $m=2$ signal and because it can potentially be long-lived, targeted searches could be devised to observe it. We estimate that with constant amplitude direct detection of the mode could occur up to a distance of roughly $14\,\mathrm{Mpc}$ whereas a search triggered by the inspiral signal could extend this distance to roughly $100\,\mathrm{Mpc}$.

We study the merger of binary neutron stars with different mass ratios adopting three different realistic, microphysical nuclear equations of state, as well as incorporating neutrino cooling effects. In particular, we concentrate on the influence of the equation of state on the gravitational wave signature and also on its role, in combination with neutrino cooling, in determining the properties of the resulting hypermassive neutron star, of the neutrinos produced, and of the ejected material. The ejecta we find are consistent with other recent studies that find that small mass ratios produce more ejecta than equal mass cases (up to some limit) and this ejecta is more neutron rich. This trend indicates the importance with future kilonovae observations of measuring the individual masses of an associated binary neutron star system, presumably from concurrent gravitational wave observations, in order to be able to extract information about the nuclear equation of state

We study the merger of binary neutron stars using different realistic, microphysical nuclear equations of state, as well as incorporating magnetic field and neutrino cooling effects. In particular, we concentrate on the influence of the equation of state on the gravitational wave signature and also on its role, in combination with cooling and electromagnetic effects, in determining the properties of the hypermassive neutron star resulting from the merger, the production of neutrinos, and the characteristics of ejecta from the system. The ejecta we find are consistent with other recent studies that find soft equations of state produce more ejecta than stiffer equations of state. Moreover, the degree of neutron richness increases for softer equations of state. In light of reported kilonova observations (associated to GRB~130603B and GRB~060614) and the discovery of relatively low abundances of heavy, radioactive elements in deep sea deposits (with respect to possible production via supernovae), we speculate that a soft EoS might be preferred---because of its significant production of sufficiently neutron rich ejecta---if such events are driven by binary neutron star mergers. We also find that realistic magnetic field strengths, obtained with a sub-grid model tuned to capture magnetic amplification via the Kelvin-Helmholtz instability at merger, are generally too weak to affect the gravitational wave signature post-merger within a time scale of $\approx 10$~ms but can have subtle effects on the post-merger dynamics.

We study the dependence of the electromagnetic luminosity --produced by interactions of force-free magnetospheres-- on dipole inclinations in binary neutron star systems. We show that this interaction extracts kinetic energy from the system and powers a Poynting flux with a strong dependence on the dipole orientations. This dependence can be linked to the reconnection and redistribution of magnetic field as the stars interact. Although the details of the Poynting luminosity are very much dependent on the orientation, all the cases considered here nevertheless radiate a large Poynting flux. This robust emission suggests that the pre-merger stage of binary neutron star systems can yield interesting electromagnetic counterparts to gravitational wave events.

The theory of eternal inflation in an inflaton potential with multiple vacua predicts that our universe is one of many bubble universes nucleating and growing inside an ever-expanding false vacuum. The collision of our bubble with another could provide an important observational signature to test this scenario. We develop and implement an algorithm for accurately computing the cosmological observables arising from bubble collisions directly from the Lagrangian of a single scalar field. We first simulate the collision spacetime by solving Einstein's equations, starting from nucleation and ending at reheating. Taking advantage of the collision's hyperbolic symmetry, simulations are performed with a 1+1-dimensional fully relativistic code that uses adaptive mesh refinement. We then calculate the comoving curvature perturbation in an open Friedmann-Robertson-Walker universe, which is used to determine the temperature anisotropies of the cosmic microwave background radiation. For a fiducial Lagrangian, the anisotropies are well described by a power law in the cosine of the angular distance from the center of the collision signature. For a given form of the Lagrangian, the resulting observational predictions are inherently statistical due to stochastic elements of the bubble nucleation process. Further uncertainties arise due to our imperfect knowledge about inflationary and pre-recombination physics. We characterize observational predictions by computing the probability distributions over four phenomenological parameters which capture these intrinsic and model uncertainties. This represents the first fully-relativistic set of predictions from an ensemble of scalar field models giving rise to eternal inflation, yielding significant differences from previous non-relativistic approximations. Thus, our results provide a basis for a rigorous confrontation of these theories with cosmological data.

We incorporate realistic, tabulated equations of state into fully relativistic simulations of magnetized neutron stars along with a neutrino leakage scheme which accounts for cooling via neutrino emission. Both these improvements utilize open-source code (GR1D) and tables from http://stellarcollapse.org. Our implementation makes use of a novel method for the calculation of the optical depth which simplifies its use with distributed adaptive mesh refinement. We present various tests with and without magnetization and preliminary results both from single stars and from the merger of a binary system.

We expand on our study of the gravitational and electromagnetic emissions from the late stage of an inspiraling neutron star binary as presented in Ref. \cite{Palenzuela:2013hu}. Interactions between the stellar magnetospheres, driven by the extreme dynamics of the merger, can yield considerable outflows. We study the gravitational and electromagnetic waves produced during the inspiral and merger of a binary neutron star system using a full relativistic, resistive MHD evolution code. We show that the interaction between the stellar magnetospheres extracts kinetic energy from the system and powers radiative Poynting flux and heat dissipation. These features depend strongly on the configuration of the initial stellar magnetic moments. Our results indicate that this power can strongly outshine pulsars in binaries and have a distinctive angular and time-dependent pattern. Our discussion provides more detail than Ref. \cite{Palenzuela:2013hu}, showing clear evidence of the different effects taking place during the inspiral. Our simulations include a few milliseconds after the actual merger and study the dynamics of the magnetic fields during the formation of the hypermassive neutron star. We also briefly discuss the possibility of observing such emissions.

The late stage of an inspiraling neutron star binary gives rise to strong gravitational wave emission due to its highly dynamic, strong gravity. Moreover, interactions between the stellar magnetospheres can produce considerable electromagnetic radiation. We study this scenario using fully general relativistic, resistive magneto-hydrodynamics simulations. We show that these interactions extract kinetic energy from the system, dissipate heat, and power radiative Poynting flux, as well as develop current sheets. Our results indicate that this power can: (i) outshine pulsars in binaries, (ii) display a distinctive angular- and time-dependent pattern, and (iii) radiate within large opening angles. These properties suggest that some binary neutron star mergers are ideal candidates for multimessenger astronomy.

Interferometric detectors will very soon give us an unprecedented view of the gravitational-wave sky, and in particular of the explosive and transient Universe. Now is the time to challenge our theoretical understanding of short-duration gravitational-wave signatures from cataclysmic events, their connection to more traditional electromagnetic and particle astrophysics, and the data analysis techniques that will make the observations a reality. This paper summarizes the state of the art, future science opportunities, and current challenges in understanding gravitational-wave transients.

We study the gravitational collapse of a magnetized neutron star using a novel numerical approach able to capture both the dynamics of the star and the behavior of the surrounding plasma. In this approach, a fully general relativistic magnetohydrodynamics implementation models the collapse of the star and provides appropriate boundary conditions to a force-free model which describes the stellar exterior. We validate this strategy by comparing with known results for the rotating monopole and aligned rotator solutions and then apply it to study both rotating and non-rotating stellar collapse scenarios, and contrast the behavior with what is obtained when employing the electrovacuum approximation outside the star. The non-rotating electrovacuum collapse is shown to agree qualitatively with a Newtonian model of the electromagnetic field outside a collapsing star. We illustrate and discuss a fundamental difference between the force-free and electrovacuum solutions, involving the appearance of large zones of electric-dominated field in the vacuum case. This provides a clear demonstration of how dissipative singularities appear generically in the non-linear time-evolution of force-free fluids. In both the rotating and non-rotating cases, our simulations indicate that the collapse induces a strong electromagnetic transient. In the case of sub-millisecond rotation, the magnetic field experiences strong winding and the transient carries much more energy. This result has important implications for models of gamma-ray bursts.

The extraction of rotational energy from a spinning black hole via the Blandford-Znajek mechanism has long been understood as an important component in models to explain energetic jets from compact astrophysical sources. Here we show more generally that the kinetic energy of the black hole, both rotational and translational, can be tapped, thereby producing even more luminous jets powered by the interaction of the black hole with its surrounding plasma. We study the resulting Poynting jet that arises from single boosted black holes and binary black hole systems. In the latter case, we find that increasing the orbital angular momenta of the system and/or the spins of the individual black holes results in an enhanced Poynting flux.

We investigate the merger of a neutron star (of compaction ratio $0.1$) in orbit about a spinning black hole in full general relativity with a mass ratio of $5:1$, allowing for the star to have an initial magnetization of $10^{12} {\rm Gauss}$. We present the resulting gravitational waveform and analyze the fallback accretion as the star is disrupted. The evolutions suggest no significant effects from the initial magnetization. We find that only a negligible amount of matter becomes unbound; $99%$ of the neutron star material has a fallback time of 10 seconds or shorter to reach the region of the central engine and that $99.99%$ of the star will interact with the central disk and black hole within 3 hours.

The interaction of black holes with ambient magnetic fields is important for a variety of highly energetic astrophysical phenomena. We study this interaction within the force-free approximation in which a tenuous plasma is assumed to have zero inertia. Blandford and Znajek (BZ) used this approach to demonstrate the conversion of some of the black hole's energy into electromagnetic Poynting flux in stationary and axisymmetric single black hole systems. We adopt this approach and extend it to examine asymmetric and, most importantly, dynamical systems by implementing the fully nonlinear field equations of general relativity coupled to Maxwell's equations. For single black holes, we study in particular the dependence of the Poynting flux and show that, even for misalignments between the black hole spin and the direction of the asymptotic magnetic field, a Poynting flux is generated with a luminosity dependent on such misalignment. For binary black hole systems, we show both in the head-on and orbiting cases that the moving black holes generate a Poynting flux.

We study the evolution of magnetized and rigidly rotating neutron stars within a fully general relativistic implementation of ideal magnetohydrodynamics with no assumed symmetries in three spatial dimensions. The stars are modeled as rotating, magnetized polytropic stars and we examine diverse scenarios to study their dynamics and stability properties. In particular we concentrate on the stability of the stars and possible critical behavior. In addition to their intrinsic physical significance, we use these evolutions as further tests of our implementation which incorporates new developments to handle magnetized systems.

Supermassive black holes are found at the centers of most galaxies and their inspiral is a natural outcome when galaxies merge. The inspiral of these systems is of utmost astrophysical importance as prodigious producers of gravitational waves and in their possible role in energetic electromagnetic events. We study such binary black hole coalescence under the influence of an external magnetic field produced by the expected circumbinary disk surrounding them. Solving the Einstein equations to describe the spacetime and using the force-free approach for the electromagnetic fields and the tenuous plasma, we present numerical evidence for possible jets driven by these systems. Extending the process described by Blandford and Znajek for a single spinning black hole, the picture that emerges suggests the electromagnetic field extracts energy from the orbiting black holes, which ultimately merge and settle into the standard Blandford-Znajek scenario. Emissions along dual and single jets would be expected that could be observable to large distances.

In addition to producing gravitational waves (GW), the dynamics of a binary black hole system could induce emission of electromagnetic (EM) radiation by affecting the behavior of plasmas and electromagnetic fields in their vicinity. We here study how the electromagnetic fields are affected by a pair of orbiting black holes through the merger. In particular, we show how the binary's dynamics induce a variability in possible electromagnetically induced emissions as well as a possible enhancement of electromagnetic fields during the late-merge and merger epochs. These time dependent features will likely leave their imprint in processes generating detectable emissions and can be exploited in the detection of electromagnetic counterparts of gravitational waves.

The merger process of a binary black hole system can have a strong impact on a circumbinary disk. In the present work we study the effect of both central mass reduction (due to the energy loss through gravitational waves) and a possible black hole recoil (due to asymmetric emission of gravitational radiation). For the mass reduction case and recoil directed along the disk's angular momentum, oscillations are induced in the disk which then modulate the internal energy and bremsstrahlung luminosities. On the other hand, when the recoil direction has a component orthogonal to the disk's angular momentum, the disk's dynamics are strongly impacted, giving rise to relativistic shocks. The shock heating leaves its signature in our proxies for radiation, the total internal energy and bremsstrahlung luminosity. Interestingly, for cases where the kick velocity is below the smallest orbital velocity in the disk (a likely scenario in real AGN), we observe a common, characteristic pattern in the internal energy of the disk. Variations in kick velocity simply provide a phase offset in the characteristic pattern implying that observations of such a signature could yield a measure of the kick velocity through electromagnetic signals alone.

We investigate the influence of magnetic fields upon the dynamics of and resulting gravitational waves from a binary neutron star merger in full general relativity coupled to ideal magnetohydrodynamics (MHD). We consider two merger scenarios, one where the stars begin with initially aligned poloidal magnetic fields and one with no magnetic field. Both mergers result in a strongly differentially rotating object. In comparison to the non-magnetized scenario, the aligned magnetic fields delay the final merger of the two stars. During and after merger we observe phenomena driven by the magnetic field, including Kelvin-Helmholtz instabilities in shear layers, winding of the field lines, and transition from poloidal to toroidal fields. These effects not only produce electromagnetic radiation, but also can have a strong influence on the gravitational waves. Thus, there are promising prospects for studying such systems with both types of waves.

We model two mergers of orbiting binary neutron stars, the first forming a black hole and the second a differentially rotating neutron star. We extract gravitational waveforms in the wave zone. Comparisons to a post-Newtonian analysis allow us to compute the orbital kinematics, including trajectories and orbital eccentricities. We verify our code by evolving single stars and extracting radial perturbative modes, which compare very well to results from perturbation theory. The Einstein equations are solved in a first order reduction of the generalized harmonic formulation, and the fluid equations are solved using a modified convex essentially non-oscillatory method. All calculations are done in three spatial dimensions without symmetry assumptions. We use the \had computational infrastructure for distributed adaptive mesh refinement.

We extend our previous studies of head-on collisions of boson stars by considering orbiting binary boson stars. We concentrate on equal mass binaries and study the dynamical behavior of boson/boson and boson/antiboson pairs. We examine the gravitational wave output of these binaries and compare with other compact binaries. Such a comparison lets us probe the apparent simplicity observed in gravitational waves produced by black hole binary systems. In our system of interest however, there is an additional internal freedom which plays a significant role in the system's dynamics, namely the phase of each star. Our evolutions show rather simple behavior at early times, but large differences occur at late times for the various initial configurations.

We study head-on collisions of boson stars in three dimensions. We consider evolutions of two boson stars which may differ in their phase or have opposite frequencies but are otherwise identical. Our studies show that these phase differences result in different late time behavior and gravitational wave output.