We perform general-relativistic simulations of charged black holes targeting GW150914. We show that the early inspiral is most efficient for detecting black hole charge through gravitational waves and that GW150914 places only weak constraints on the charge-to-mass ratio, ${Q}/{M}< 0.4$. Our work applies to electric and magnetic charge, and to theories with black holes endowed with U(1) (hidden or dark) charges. With our results, we place an upper bound on the deviation of the so-called theory of MOdified Gravity (MOG) from general relativity in the dynamical, strong-field regime.

Spin can have significant effects on the electromagnetic transients accompanying binary neutron star mergers. The measurement of spin can provide important information about binary formation channels. In the absence of a strong neutron star spin prior, the degeneracy of spin with other parameters leads to significant uncertainties in their estimation, in particular limiting the power of gravitational waves to place tight constraints on the nuclear equation of state. Thus detailed studies of highly spinning neutron star mergers are essential to understand all aspects of multimessenger observation of such events. We perform a systematic investigation of the impact of neutron star spin---considering dimensionless spin values up to $a_{\rm NS}=0.33$---on the merger of equal mass, quasicircular binary neutron stars using fully general-relativistic simulations. We find that the peak frequency of the post-merger gravitational wave signal is only weakly influenced by the neutron star spin, with cases where the spin is aligned (antialigned) with the orbital angular momentum giving slightly lower (higher) values compared to the irrotational case. We find that the one-arm instability arises in a number of cases, with some dependence on spin. Spin has a pronounced impact on the mass, velocity, and angular distribution of the dynamical ejecta, and the mass of the disk that remains outside the merger remnant. We discuss the implications of these findings on anticipated electromagnetic signals, and on constraints that have been placed on the equation of state based on multimessenger observations of GW170817.

Inspiraling and merging binary neutron stars (BNSs) are important sources of both gravitational waves and coincident electromagnetic counterparts. If the BNS total mass is larger than a threshold value, a black hole ensues promptly after merger. Through a statistical study in conjunction with recent LIGO/Virgo constraints on the nuclear equation of state, we estimate that up to $\sim 25\%$ of BNS mergers may result in prompt collapse. Moreover, we find that most models of the BNS mass function we study here predict that the majority of prompt-collapse BNS mergers have $q\gtrsim 0.8$. Prompt-collapse BNS mergers with mass ratio $q \gtrsim 0.8$ may not be accompanied by detectable kilonovae or short gamma-ray bursts, because they unbind a negligible amount of mass and form negligibly small accretion disks onto the remnant black hole. We call such BNS mergers "orphan". However, recent studies have found that ${10^{41-43}(B_p/10^{12}\rm G)^2 erg\, s^{-1}}$ electromagnetic signals can be powered by magnetospheric interactions several milliseconds prior to merger. Moreover, the energy stored in the magnetosphere of an orphan BNS merger remnant will be radiated away in ${\mathcal O}(1\ \rm ms)$. Through simulations in full general relativity of BNSs endowed with an initial dipole magnetosphere, we find that the energy in the magnetosphere following black hole formation is $E_B \sim 10^{39-41} (B_p/10^{12}\rm G)^2$ erg. Radiating $\sim 1\%$ of $E_B$ in 1 ms, as has been found in previous studies, matches the premerger magnetospheric luminosity. These magnetospheric signals are not beamed, and their duration and power agrees with those of non-repeating fast radio bursts (FRBs). These results combined with our statistical study suggest that a non-repeating FRB may be the most likely electromagnetic counterpart of prompt-collapse BNSs.

We present Searches After Gravitational-waves Using ARizona Observatories (SAGUARO), a comprehensive effort dedicated to the discovery and characterization of optical counterparts to gravitational wave (GW) events. SAGUARO utilizes ground-based facilities ranging from 1.5m to 10m in diameter, located primarily in the Northern Hemisphere. We provide an overview of SAGUARO's telescopic resources, pipeline for transient detection, and database for candidate visualization. We describe SAGUARO's discovery component, which utilizes the $5$~deg$^2$ field-of-view optical imager on the Mt. Lemmon 1.5m telescope, reaching limits of $\approx 21.3$~AB mag while rapidly tiling large areas. We also describe the follow-up component of SAGUARO, used for rapid vetting and monitoring of optical candidates. With the onset of Advanced LIGO/Virgo's third observing run, we present results from the first three SAGUARO searches following the GW events S190408an, S190425z and S190426c, which serve as a valuable proof-of-concept of SAGUARO. We triggered and searched 15, 60 and 60 deg$^{2}$ respectively, 17.6, 1.4 and 41.8 hrs after the initial GW alerts. We covered 7.8, 3.0 and 5.1\% of the total probability within the GW event localizations, reaching 3$\sigma$ limits of 19.8, 21.3 and 20.8 AB mag, respectively. Although no viable counterparts associated with these events were found, we recovered 6 known transients and ruled out 5 potential candidates. We also present Large Binocular Telescope spectroscopy of PS19eq/SN2019ebq, a promising kilonova candidate that was later determined to be a supernova. With the ability to tile large areas and conduct detailed follow-up, SAGUARO represents a significant addition to GW counterpart searches.

We investigate the dynamical stability of relativistic, differentially rotating, quasi-toroidal models of neutron stars through hydrodynamical simulations in full general relativity. We find that all quasi-toroidal configurations studied in this work are dynamically unstable against the growth of non-axisymmetric modes. Both one-arm and bar mode instabilities grow during their evolution. We find that very high rest mass configurations collapse to form black holes. Our calculations suggest that configurations whose rest mass is less than the binary neutron star threshold mass for prompt collapse to black hole transition dynamically to spheroidal, differentially rotating stars that are dynamically stable, but secularly unstable. Our study shows that the existence of extreme quasi-toroidal neutron star equilibrium solutions does not imply that long-lived binary neutron star merger remnants can be much more massive than previously found. Finally, we find models that are initially supra-Kerr ($J/M^2>1$) and undergo catastrophic collapse on a dynamical timescale, in contrast to what was found in earlier works. However, cosmic censorship is respected in all of our cases. Our work explicitly demonstrates that exceeding the Kerr bound in rotating neutron star models does not imply dynamical stability.

A general relativistic, stationary and axisymmetric black hole in a four-dimensional asymptotically-flat spacetime is fully determined by its mass, angular momentum and electric charge. The expectation that astrophysically relevant black holes do not posses charge has resulted in a limited number of investigations of moving and charged black holes in the dynamical, strong-field gravitational (and electromagnetic) regime, where numerical studies are necessary. Apart from having a theoretical interest, the advent of multimessenger astronomy with gravitational waves offers new ways to think about charged black holes. In this work, we initiate an exploration of charged binary black holes by generating valid initial data for general relativistic simulations of black hole systems that have generic electric charge, linear and angular momenta. We develop our initial data formalism within the framework of the conformal transverse-traceless (Bowen-York) technique using the puncture approach, and apply the theory of isolated horizons to attribute physical parameters (mass, charge and angular momentum) to each hole. We implemented our formalism in the case of a binary system by modifying the publicly available TwoPunctures and QuasiLocalMeasures codes. We demonstrate that our code can recover existing solutions and that it has excellent self-convergence properties for a generic configuration of two black holes.

We perform long-term simulations of spinning binary neutron stars, with our highest dimensionless spin being $\chi \sim 0.32$. To assess the importance of spin during the inspiral we vary the spin, and also use two equations of state, one that consists of plain nuclear matter and produces compact stars (SLy), and a hybrid one that contains both nuclear and quark matter and leads to larger stars (ALF2). Using high resolution that has grid spacing $\Delta x\sim 98$ m on the finest refinement level, we find that the effects of spin in the phase evolution of a binary system can be larger than the one that comes from tidal forces. Our calculations demonstrate explicitly that although tidal effects are dominant for small spins ($\lesssim 0.1$), this is no longer true when the spins are larger, but still much smaller than the Keplerian limit.

We construct equilibrium models of uniformly and differentially rotating hybrid hadron-quark stars using equations of state (EOSs) with a first-order phase transition that gives rise to a third family of compact objects. We find that the ratio of the maximum possible mass of uniformly rotating configurations - the supramassive limit - to the Tolman-Oppenheimer-Volkoff (TOV) limit mass is not EOS-independent, and is between 1.15 and 1.31,in contrast with the value of 1.20 previously found for hadronic EOSs. Therefore, some of the constraints placed on the EOS from the observation of the gravitational wave event GW170817 do not apply to hadron-quark EOSs. However, the supramassive limit mass for the family of EOSs we treat is consistent with limits set by GW170817, strengthening the possibility of interpreting GW170817 with a hybrid hadron-quark EOSs. We also find that along constant angular momentum sequences of uniformly rotating stars, the third family maximum and minimum mass models satisfy approximate EOS-independent relations, and the supramassive limit of the third family is approximately 16.5 % larger than the third family TOV limit. For differentially rotating spheroidal stars, we find that a lower-limit on the maximum supportable rest mass is 123 % more than the TOV limit rest mass. Finally, we verify that the recently discovered universal relations relating angular momentum, rest mass and gravitational mass for turning-point models hold for hybrid hadron-quark EOSs when uniform rotation is considered, but have a clear dependence on the degree of differential rotation.

Events GW170817 and GRB 170817A provide the best confirmation so far that compact binary mergers where at least one of the companions is a neutron star (NS) can be the progenitors of short gamma-ray bursts (sGRBs). An open question for GW170817 remains the values and impact of the initial NS spins. The initial spins could possibly affect the remnant black hole (BH) mass and spin, the remnant disk and the formation and lifetime of a jet and its luminosity. Here we summarize our general relativistic magnetohydrodynamic simulations of spinning, NS binaries undergoing merger and delayed collapse to a BH. The binaries consist of two identical NSs, modeled as $\Gamma=2$ polytropes, in quasicircular orbit, each with spins $\chi_{\rm{NS}}=-0.053,\,0,\,0.24$, or $0.36$. The stars are endowed initially with a dipolar magnetic field extending from the interior into the exterior, as in a radio pulsar. Following merger, the redistribution of angular momentum by magnetic braking and magnetic turbulent viscosity in the hypermassive neutron star (HMNS) remnant, along with the loss of angular momentum due to gravitational radiation, induce the formation of a massive, nearly uniformly rotating inner core surrounded by a magnetized Keplerian disk-like envelope. The HMNS eventually collapses to a BH, with spin $a/M_{\rm BH} \simeq 0.78$ independent of the initial spin of the NSs, surrounded by a magnetized accretion disk. The larger the initial NS spin the heavier the disk. After $\Delta t\sim 3000-4000 M \sim 45-60(M_{\rm NS}/1.625M_\odot)\rm ms$ following merger, a mildly relativistic jet is launched. The lifetime of the jet [$\Delta t\sim 100-140(M_{\rm NS}/1.625M_\odot)\rm ms$] and its outgoing Poynting luminosity [$L_{\rm EM}\sim 10^{51.5\pm 1}\rm erg/s$] are consistent with typical sGRBs, as well as with the Blandford--Znajek mechanism for launching jets and their associated Poynting luminosities.

Pulsar timing arrays (PTAs) are on the verge of detecting low-frequency gravitational waves (GWs) from supermassive black hole binaries (SMBHBs). With continued observations of a large sample of millisecond pulsars, PTAs will reach this major milestone within the next decade. Already, SMBHB candidates are being identified by electromagnetic surveys in ever-increasing numbers; upcoming surveys will enhance our ability to detect and verify candidates, and will be instrumental in identifying the host galaxies of GW sources. Multi-messenger (GW and electromagnetic) observations of SMBHBs will revolutionize our understanding of the co-evolution of SMBHs with their host galaxies, the dynamical interactions between binaries and their galactic environments, and the fundamental physics of accretion. Multi-messenger observations can also make SMBHBs 'standard sirens' for cosmological distance measurements out to $z\simeq0.5$. LIGO has already ushered in breakthrough insights in our knowledge of black holes. The multi-messenger detection of SMBHBs with PTAs will be a breakthrough in the years $2020-2030$ and beyond, and prepare us for LISA to help complete our views of black hole demographics and evolution at higher redshifts.

We study the solution space of general relativistic, axisymmetric, equilibria of differentially rotating neutron stars with realistic, nuclear equations of state. We find that different types of stars, which were identified by earlier works for polytropic equations of state, arise for realistic equations of state, too. Scanning the solution space for the sample of realistic equations of state we treat, we find lower limits on the maximum rest masses supported by cold, differentially rotating stars for each type of stars. We often discover equilibrium configurations that can support more than 2 times the mass of a static star. We call these equilibria "\"ubermassive", and in our survey we find \"ubermassive stars that can support up to 2.5 times the maximum rest mass that can be supported by a cold, non-rotating star with the same equation of state. This is nearly two times larger than what previous studies employing realistic equations of state had found, and which did not uncover \"ubermassive neutron stars. Moreover, we find that the increase in the maximum rest mass with respect to the non-spinning stellar counterpart is larger for moderately stiff equations of state. These results may have implications for the lifetime and the gravitational wave and electromagnetic counterparts of hypermassive neutron stars formed following binary neutron star mergers.

We report the analysis of the gamma-ray variability of NGC 1275--the radio galaxy at the center of the Perseus cluster. NGC 1275 has been observed continuously with the Fermi Large Area Telescope over the last nine years. We applied different time-domain analysis methods including Fourier, wavelets and Bayesian methods, in order to search for quasi-periodic oscillations in the gamma-ray emission. We found no evidence for periodicities of astrophysical origin.

This work is the first in a series of studies aimed at understanding the dynamics of highly eccentric binary neutron stars, and constructing an appropriate gravitational-waveform model for detection. Such binaries are possible sources for ground-based gravitational wave detectors, and are expected to form through dynamical scattering and multi-body interactions in globular clusters and galactic nuclei. In contrast to black holes, oscillations of neutron stars are generically excited by tidal effects after close pericenter passage. Depending on the equation of state, this can enhance the loss of orbital energy by up to tens of percent over that radiated away by gravitational waves during an orbit. Under the same interaction mechanism, part of the orbital angular momentum is also transferred to the star. We calculate the impact of the neutron star oscillations on the orbital evolution of such systems, and compare these results to full numerical simulations. Utilizing a Post-Newtonian flux description we propose a preliminary model to predict the timing of different pericenter passages. A refined version of this model (taking into account Post-Newtonian corrections to the tidal coupling and the oscillations of the stars) may serve as a waveform model for such highly eccentric systems.

Gravitational wave observations of GW170817 placed bounds on the tidal deformabilities of compact stars allowing one to probe equations of state for matter at supranuclear densities. Here we design new parametrizations for hybrid hadron-quark equations of state, that give rise to low-mass twin stars, and test them against GW170817. We find that GW170817 is consistent with the coalescence of a binary hybrid star--neutron star. We also test and find that the I-Love-Q relations for hybrid stars in the third family agree with those for purely hadronic and quark stars within $\sim 3\%$ for both slowly and rapidly rotating configurations, implying that these relations can be used to perform equation-of-state independent tests of general relativity and to break degeneracies in gravitational waveforms for hybrid stars in the third family as well.

We perform magnetohydrodynamic simulations in full general relativity of disk accretion onto nonspinning black hole binaries with mass ratio 36:29. We survey different disk models which differ in their scale height, total size and magnetic field to quantify the robustness of previous simulations on the initial disk model. Scaling our simulations to LIGO GW150914 we find that such systems could explain possible gravitational wave and electromagnetic counterparts such as the Fermi GBM hard X-ray signal reported 0.4s after GW150915 ended. Scaling our simulations to supermassive binary black holes, we find that observable flow properties such as accretion rate periodicities, the emergence of jets throughout inspiral, merger and post-merger, disk temperatures, thermal frequencies, and the time-delay between merger and the boost in jet outflows that we reported in earlier studies display only modest dependence on the initial disk model we consider here.

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.

We perform MHD simulations in full GR of uniformly rotating stars that are marginally unstable to collapse. Our simulations model the direct collapse of supermassive stars (SMSs) to seed black holes (BHs) that can grow to become the supermassive BHs at the centers of quasars and AGNs. They also crudely model the collapse of massive Pop III stars to BHs, which could power a fraction of distant, long gamma-ray bursts (GRBs). The initial stellar models we adopt are $\Gamma = 4/3$ polytropes seeded with a dynamically unimportant dipole magnetic field (B field). We treat initial B-field configurations either confined to the stellar interior or extending out from the interior into the stellar exterior. The BH formed following collapse has mass $M_{BH} \simeq 0.9M$ (where $M$ is the mass of the initial star) and spin $a_{BH}/M_{BH}\simeq 0.7$. A massive, hot, magnetized torus surrounds the remnant BH. At $\Delta t\sim 400-550M\approx 2000-2700(M/10^6M_\odot)$s following the gravitational wave (GW) peak amplitude, an incipient jet is launched. The disk lifetime is $\Delta t\sim 10^5(M/10^6M_\odot)$s, and the jet luminosity is $L_{EM}\sim 10^{51-52}$ ergs/s. If $\gtrsim 1-10\%$ of this power is converted into gamma rays, SWIFT and FERMI could potentially detect these events out to large redshifts $z\sim 20$. Thus, SMSs could be sources of ultra-long GRBs and massive Pop III stars could be the progenitors that power a fraction of the long GRBs observed at redshift $z \sim 5-8$. GWs are copiously emitted during the collapse, and peak at $\sim 15(10^6 M_{\odot}/M)\rm mHz$ ($\sim 0.15(10^4 M_{\odot}/M)\rm Hz$), i.e., in the LISA (DECIGO/BBO) band; optimally oriented SMSs could be detectable by LISA (DECIGO/BBO) at $z \lesssim 3$ ($z \lesssim 11$). Hence $10^4 M_{\odot}$ SMSs collapsing at $z\sim 10$ are promising multimessenger sources of coincident gravitational and electromagnetic waves.

Targets for ground-based gravitational wave interferometers include continuous, quasiperiodic sources of gravitational radiation, such as isolated, spinning neutron stars. In this work we perform evolution simulations of uniformly rotating, triaxially deformed stars, the compressible analogues in general relativity of incompressible, Newtonian Jacobi ellipsoids. We investigate their stability and gravitational wave emission. We employ five models, both normal and supramassive, and track their evolution with different grid setups and resolutions, as well as with two different evolution codes. We find that all models are dynamically stable and produce a strain that is approximately one-tenth the average value of a merging binary system. We track their secular evolution and find that all our stars evolve towards axisymmetry, maintaining their uniform rotation, kinetic energy, and angular momentum profiles while losing their triaxiality.

Black hole - neutron star (BHNS) and neutron star - neutron star (NSNS) binaries are among the favored candidates for the progenitors of the black hole - disk systems that may be the engines powering short-hard gamma ray bursts. After almost two decades of simulations of binary NSNSs and BHNSs in full general relativity we are now beginning to understand the ingredients that may be necessary for these systems to launch incipient jets. Here, we review our current understanding, and summarize the surprises and lessons learned from state-of-the-art (magnetohydrodynamic) simulations in full general relativity of BHNS and NSNS mergers as jet engines for short-hard gamma-ray bursts.

We continue our investigations of the development and importance of the one-arm spiral instability in long-lived hypermassive neutron stars (HMNSs) formed in dynamical capture binary neutron star mergers. Employing hydrodynamic simulations in full general relativity, we find that the one-arm instability is generic in that it can develop in HMNSs within a few tens of milliseconds after merger for all equations of state in our survey. We find that mergers with stiffer equations of state tend to produce HMNSs with stronger $m=2$ azimuthal mode density deformations, and weaker $m=1$ components, relative to softer equations of state. We also find that for equations of state that can give rise to double-core HMNSs, large $m=1$ density modes can already be present due to asymmetries in the two cores. This results in the generation of $l=2$, $m=1$ gravitational wave modes even before the dominance of a one-arm mode that ultimately arises following merger of the two cores. Our results further suggest that stiffer equations of state give rise to HMNSs generating lower $m=1$ gravitational wave frequencies. Thus, if gravitational waves from the one-arm instability are detected, they could in principle constrain the neutron star equation of state. We estimate that, depending on the equation of state, the one-arm mode could potentially be detectable by second generation gravitational wave detectors at $\sim 10$ Mpc and by third generation ones at $\sim 100$ Mpc. Finally, we provide estimates of the properties of dynamical ejecta, as well as the accompanying kilonovae signatures.

We perform magnetohydrodynamic simulations in full general relativity (GRMHD) of quasi-circular, equal-mass, binary neutron stars that undergo merger. The initial stars are irrotational, $n=1$ polytropes and are magnetized. We explore two types of magnetic-field geometries: one where each star is endowed with a dipole magnetic field extending from the interior into the exterior, as in a pulsar, and the other where the dipole field is initially confined to the interior. In both cases the adopted magnetic fields are initially dynamically unimportant. The merger outcome is a hypermassive neutron star that undergoes delayed collapse to a black hole (spin parameter $a/M_{\rm BH} \sim 0.74$) immersed in a magnetized accretion disk. About $4000M \sim 60(M_{\rm NS}/1.625M_\odot)$ ms following merger, the region above the black hole poles becomes strongly magnetized, and a collimated, mildly relativistic outflow --- an incipient jet --- is launched. The lifetime of the accretion disk, which likely equals the lifetime of the jet, is $\Delta t \sim 0.1 (M_{\rm NS}/1.625M_\odot)$ s. In contrast to black hole--neutron star mergers, we find that incipient jets are launched even when the initial magnetic field is confined to the interior of the stars.

We perform general-relativistic hydrodynamical simulations of dynamical capture binary neutron star mergers, emphasizing the role played by the neutron star spin. Dynamical capture mergers may take place in globular clusters, as well as other dense stellar systems, where most neutron stars have large spins. We find significant variability in the merger outcome as a function of initial neutron star spin. For cases where the spin is aligned with the orbital angular momentum, the additional centrifugal support in the remnant hypermassive neutron star can prevent the prompt collapse to a black hole, while for antialigned cases the decreased total angular momentum can facilitate the collapse to a black hole. We show that even moderate spins can significantly increase the amount of ejected material, including the amount unbound with velocities greater than half the speed of light, leading to brighter electromagnetic signatures associated with kilonovae and interaction of the ejecta with the interstellar medium. Furthermore, we find that the initial neutron star spin can strongly affect the already rich phenomenology in the postmerger gravitational wave signatures that arise from the oscillation modes of the hypermassive neutron star. In several of our simulations, the resulting hypermassive neutron star develops the one-arm ($m=1$) spiral instability, the most pronounced cases being those with small but non-negligible neutron star spins. For long-lived hypermassive neutron stars, the presence of this instability leads to improved prospects for detecting these events through gravitational waves, and thus may give information about the neutron star equation of state.

Using general-relativistic hydrodynamical simulations, we show that merging binary neutron stars can form hypermassive neutrons stars that undergo the one-arm spiral instability. We study the particular case of a dynamical capture merger where the stars have a small spin, as may arise in globular clusters, and focus on an equal-mass scenario where the spins are aligned with the orbital angular momentum. We find that this instability develops when post-merger fluid vortices lead to the generation of a toroidal remnant - a configuration whose maximum density occurs in a ring around the center-of-mass - with high vorticity along its rotation axis. The instability quickly saturates on a timescale of $\sim 10$ ms, with the $m=1$ azimuthal density multipole mode dominating over higher modes. The instability also leaves a characteristic imprint on the post-merger gravitational wave signal that could be detectable if the instability persists in long-lived remnants.

In the extreme violence of merger and mass accretion, compact objects like black holes and neutron stars are thought to launch some of the most luminous outbursts of electromagnetic and gravitational wave energy in the Universe. Modeling these systems realistically is a central problem in theoretical astrophysics, but has proven extremely challenging, requiring the development of numerical relativity codes that solve Einstein's equations for the spacetime, coupled to the equations of general relativistic (ideal) magnetohydrodynamics (GRMHD) for the magnetized fluids. Over the past decade, the Illinois Numerical Relativity (ILNR) Group's dynamical spacetime GRMHD code has proven itself as a robust and reliable tool for theoretical modeling of such GRMHD phenomena. However, the code was written "by experts and for experts" of the code, with a steep learning curve that would severely hinder community adoption if it were open-sourced. Here we present IllinoisGRMHD, which is an open-source, highly-extensible rewrite of the original closed-source GRMHD code of the ILNR Group. Reducing the learning curve was the primary focus of this rewrite, with the goal of facilitating community involvement in the code's use and development, as well as the minimization of human effort in generating new science. IllinoisGRMHD also saves computer time, generating roundoff-precision identical output to the original code on adaptive-mesh grids, but nearly twice as fast at scales of hundreds to thousands of cores.

We study dynamical capture binary black hole-neutron star (BH-NS) mergers focusing on the effects of the neutron star spin. These events may arise in dense stellar regions, such as globular clusters, where the majority of neutron stars are expected to be rapidly rotating. We initialize the BH-NS systems with positions and velocities corresponding to marginally unbound Newtonian orbits, and evolve them using general-relativistic hydrodynamical simulations. We find that even moderate spins can significantly increase the amount of mass in unbound material. In some of the more extreme cases, there can be up to a third of a solar mass in unbound matter. Similarly, large amounts of tidally stripped material can remain bound and eventually accrete onto the BH---as much as a tenth of a solar mass in some cases. These simulations demonstrate that it is important to treat neutron star spin in order to make reliable predictions of the gravitational wave and electromagnetic transient signals accompanying these sources.

We perform magnetohydrodynamic simulations in full general relativity (GRMHD) of a binary black hole-neutron star on a quasicircular orbit that undergoes merger. The binary mass ratio is 3:1, the black hole initial spin parameter $a/m=0.75$ ($m$ is the black hole Christodoulou mass) aligned with the orbital angular momentum, and the neutron star is an irrotational $\Gamma=2$ polytrope. About two orbits prior to merger (at time $t=t_B$), we seed the neutron star with a dynamically weak interior dipole magnetic field that extends into the stellar exterior. At $t=t_B$ the exterior has a low-density atmosphere with constant plasma parameter $\beta\equiv P_{\rm gas}/P_{\rm mag}$. Varying $\beta$ at $t_B$ in the exterior from $0.1$ to $0.01$, we find that at a time $\sim 4000{\rm M} \sim 100(M_{\rm NS}/1.4M_\odot)$ms [M is the total (ADM) mass] following the onset of accretion of tidally disrupted debris, magnetic winding above the remnant black hole poles builds up the magnetic field sufficiently to launch a mildly relativistic, collimated outflow - an incipient jet. The duration of the accretion and the lifetime of the jet is $\Delta t\sim 0.5(M_{\rm NS}/1.4M_\odot)$s. Our simulations furnish the first explicit examples in GRMHD which show that a jet can emerge following a black hole - neutron star merger.

We report results from simulations in general relativity of magnetized disks accreting onto merging black hole binaries, starting from relaxed disk initial data. The simulations feature an effective, rapid radiative cooling scheme as a limiting case of future treatments with radiative transfer. Here we evolve the systems after binary-disk decoupling through inspiral and merger, and analyze the dependence on the binary mass ratio with $q\equiv m_{\rm bh}/M_{\rm BH}=1,1/2,$ and $1/4$. We find that the luminosity associated with local cooling is larger than the luminosity associated with matter kinetic outflows, while the electromagnetic (Poynting) luminosity associated with bulk transport of magnetic field energy is the smallest. The cooling luminosity around merger is only marginally smaller than that of a single, non-spinning black hole. Incipient jets are launched independently of the mass ratio, while the same initial disk accreting on a single non-spinning black hole does not lead to a jet, as expected. For all mass ratios we see a transient behavior in the collimated, magnetized outflows lasting $2-5 ( M/10^8M_\odot ) \rm days$ after merger: the outflows become increasingly magnetically dominated and accelerated to higher velocities, boosting the Poynting luminosity. These sudden changes can alter the electromagnetic emission across the jet and potentially help distinguish mergers of black holes in AGNs from single accreting black holes based on jet morphology alone.

Robust gauge conditions are critically important to the stability and accuracy of numerical relativity (NR) simulations involving compact objects. Most of the NR community use the highly robust---though decade-old---moving-puncture (MP) gauge conditions for such simulations. It has been argued that in binary black hole (BBH) evolutions adopting this gauge, noise generated near adaptive-mesh-refinement (AMR) boundaries does not converge away cleanly with increasing resolution, severely limiting gravitational waveform accuracy at computationally feasible resolutions. We link this noise to a sharp (short-wavelength), initial outgoing gauge wave crossing into progressively lower resolution AMR grids, and present improvements to the standard MP gauge conditions that focus on stretching, smoothing, and more rapidly settling this outgoing wave. Our best gauge choice greatly reduces gravitational waveform noise during inspiral, yielding less fluctuation in convergence order and $\sim 40%$ lower waveform phase and amplitude errors at typical resolutions. Noise in other physical quantities of interest is also reduced, and constraint violations drop by more than an order of magnitude. We expect these improvements will carry over to simulations of all types of compact binary systems, as well as other $N$+1 formulations of gravity for which MP-like gauge conditions can be chosen.

Adopting our new method for matching general relativistic, ideal magnetohydrodynamics to its force-free limit, we perform the first systematic simulations of force-free pulsar magnetospheres in general relativity. We endow the neutron star with a general relativistic dipole magnetic field, model the interior with ideal magnetohydrodynamics, and adopt force-free electrodynamics in the exterior. Comparing the spin-down luminosity to its corresponding Minkowski value, we find that general relativistic effects give rise to a modest enhancement: the maximum enhancement for $n=1$ polytropes is $\sim 23\%$. Evolving a rapidly rotating $n=0.5$ polytrope we find an even greater enhancement of $\sim 35\%$. Using our simulation data, we derive fitting formulas for the pulsar spin-down luminosity as a function of the neutron star compaction, angular speed, and dipole magnetic moment. We expect stiffer equations of state and more rapidly spinning neutron stars to lead to even larger enhancements in the spin-down luminosity.

We adopt the conduction fluid approximation to model the steady-state distribution of matter around a massive black hole at the center of a weakly collisional cluster of particles. By "`weakly collisional" we mean a cluster in which the mean free time between particle collisions is much longer than the characteristic particle crossing (dynamical) time scale, but shorter than the cluster lifetime. When applied to a star cluster, we reproduce the familiar Bahcall-Wolf power-law cusp solution for the stars bound to the black hole. Here the star density scales with radius as $r^{-7/4}$ and the velocity dispersion as $r^{-1/2}$ throughout most of the gravitational well of the black hole. When applied to a relaxed, self-interacting dark matter (SIDM) halo with a velocity-dependent cross section $\sigma \sim v^{-a}$, the gas again forms a power-law cusp, but now the SIDM density scales as $r^{-\beta}$, where $\beta = (a+3)/4$, while its velocity dispersion again varies as $r^{-1/2}$. Results are obtained first in Newtonian theory and then in full general relativity. Although the conduction fluid model is a simplification, it provides a reasonable first approximation to the matter profiles and is much easier to implement than a full Fokker-Planck treatment or an $N$-body simulation of the Boltzmann equation with collisional perturbations.

The Numerical INJection Analysis (NINJA) project is a collaborative effort between members of the numerical relativity and gravitational-wave astrophysics communities. The purpose of NINJA is to study the ability to detect gravitational waves emitted from merging binary black holes and recover their parameters with next-generation gravitational-wave observatories. We report here on the results of the second NINJA project, NINJA-2, which employs 60 complete binary black hole hybrid waveforms consisting of a numerical portion modelling the late inspiral, merger, and ringdown stitched to a post-Newtonian portion modelling the early inspiral. In a "blind injection challenge" similar to that conducted in recent LIGO and Virgo science runs, we added 7 hybrid waveforms to two months of data recolored to predictions of Advanced LIGO and Advanced Virgo sensitivity curves during their first observing runs. The resulting data was analyzed by gravitational-wave detection algorithms and 6 of the waveforms were recovered with false alarm rates smaller than 1 in a thousand years. Parameter estimation algorithms were run on each of these waveforms to explore the ability to constrain the masses, component angular momenta and sky position of these waveforms. We also perform a large-scale monte-carlo study to assess the ability to recover each of the 60 hybrid waveforms with early Advanced LIGO and Advanced Virgo sensitivity curves. Our results predict that early Advanced LIGO and Advanced Virgo will have a volume-weighted average sensitive distance of 300Mpc (1Gpc) for $10M_{\odot}+10M_{\odot}$ ($50M_{\odot}+50M_{\odot}$) binary black hole coalescences. We demonstrate that neglecting the component angular momenta in the waveform models used in matched-filtering will result in a reduction in sensitivity for systems with large component angular momenta. [Abstract abridged for ArXiv, full version in PDF]

We report on simulations in general relativity of magnetized disks onto black hole binaries. We vary the binary mass ratio from 1:1 to 1:10 and evolve the systems when they orbit near the binary-disk decoupling radius. We compare (surface) density profiles, accretion rates (relative to a single, non-spinning black hole), variability, effective $\alpha$-stress levels and luminosities as functions of the mass ratio. We treat the disks in two limiting regimes: rapid radiative cooling and no radiative cooling. The magnetic field lines clearly reveal jets emerging from both black hole horizons and merging into one common jet at large distances. The magnetic fields give rise to much stronger shock heating than the pure hydrodynamic flows, completely alter the disk structure, and boost accretion rates and luminosities. Accretion streams near the horizons are among the densest structures; in fact, the 1:10 no-cooling evolution results in a refilling of the cavity. The typical effective temperature in the bulk of the disk is $\sim 10^5 (M/10^8 M_\odot)^{-1/4} (L/L_{\rm edd})^{1/4} {\rm K}$ yielding characteristic thermal frequencies $\sim 10^{15} (M/10^8 M_\odot)^{-1/4} (L/L_{\rm edd})^{1/4}(1+z)^{-1}{\rm Hz} $. These systems are thus promising targets for many extragalactic optical surveys, such as LSST, WFIRST, and PanSTARRS.

We present a new computational method for smoothly matching general relativistic ideal magnetohydrodynamics (MHD) to its force-free limit. The method is based on a flux-conservative formalism for MHD and its force-free limit, and a vector potential formulation for the induction equation to maintain the zero divergence constraint for the magnetic field. The force-free formulation we adopt evolves the magnetic field and the Poynting vector, instead of the magnetic and electric fields. We show that our force-free code passes a robust suite of tests, performed both in 1D flat spacetime and in 3D black hole spacetimes. We also demonstrate that our matching technique successfully reproduces the aligned rotator force-free solution. Our new techniques are suitable for studying electromagnetic effects and predicting electromagnetic signals arising in many different curved spacetime scenarios. For example, we can treat spinning neutron stars, either in isolation or in compact binaries, that have MHD interiors and force-free magnetospheres.

We perform the first general relativistic force-free simulations of neutron star (NS) magnetospheres in orbit about spinning and non-spinning black holes. We find promising precursor electromagnetic emission: typical Poynting luminosities at, e.g., an orbital separation of 6.6 times the NS radius are L ~ 6 x 10^{42} erg/s for a 1.4 solar-mass NS with a 10^{13}G polar magnetic field. The Poynting flux peaks within a broad beam of ~40 degrees in the azimuthal direction and within ~60 degrees from the orbital plane, establishing a possible lighthouse effect. Our calculations, though preliminary, preview more detailed simulations of these systems that we plan to perform in the future.

TwoPunctures is perhaps the most widely-adopted code for generating binary black hole "puncture" initial data and interpolating these (spectral) data onto evolution grids. In typical usage, the bulk of this code's run time is spent in its spectral interpolation routine. We announce a new publicly-available spectral interpolation routine that improves the performance of the original interpolation routine by a factor of ~100, yielding results consistent with the original spectral interpolation routine to roundoff precision. This note serves as a guide for installing this routine both in the original standalone TwoPunctures code and the Einstein Toolkit supported version of this code.

We present new numerical techniques we developed for launching the first parameter study of {\it magnetized} black hole--neutron star (BHNS) mergers, varying the magnetic fields seeded in the initial neutron star. We found that magnetic fields have a negligible impact on the gravitational waveforms and bulk dynamics of the system during merger, regardless of magnetic field strength or BH spin. In a recent simulation, we seeded the remnant disk from an unmagnetized BHNS merger simulation with large-scale, purely poloidal magnetic fields, which are otherwise absent in the full simulation. The outcome appears to be a viable sGRB central engine.

We present results from the first fully general relativistic, magnetohydrodynamic (GRMHD) simulations of an equal-mass black hole binary (BHBH) in a magnetized, circumbinary accretion disk. We simulate both the pre and post-decoupling phases of a BHBH-disk system and both "cooling" and "no-cooling" gas flows. Prior to decoupling, the competition between the binary tidal torques and the effective viscous torques due to MHD turbulence depletes the disk interior to the binary orbit. However, it also induces a two-stream accretion flow and mildly relativistic polar outflows from the BHs. Following decoupling, but before gas fills the low-density "hollow" surrounding the remnant, the accretion rate is reduced, while there is a prompt electromagnetic (EM) luminosity enhancement following merger due to shock heating and accretion onto the spinning BH remnant. This investigation, though preliminary, previews more detailed GRMHD simulations we plan to perform in anticipation of future, simultaneous detections of gravitational and EM radiation from a merging BHBH-disk system.

Black hole--neutron star (BHNS) binary mergers can form disks in which magnetorotational instability (MRI)-induced turbulence may drive accretion onto the remnant BH, supporting relativistic jets and providing the engine for a short-hard gamma-ray burst (SGRB). Our earlier study of magnetized BHNSs showed that NS tidal disruption winds the magnetic field into a toroidal configuration, with poloidal fields so weak that capturing MRI with full-disk simulations would require $\sim 10^8$ CPU-hours. In that study we imposed equatorial symmetry, suppressing poloidal magnetic fields that might be generated from plasma crossing the orbital plane. Here we show that initial conditions that break this symmetry (i.e., {\it tilted} poloidal magnetic fields in the NS) generate much stronger poloidal fields in the disk, indicating that asymmetric initial conditions may be necessary for establishing BHNS mergers as SGRB progenitors via fully general relativistic MHD simulations. We demonstrate that BHNS mergers may form an SGRB engine under the right conditions by seeding the remnant disk from an unmagnetized BHNS simulation with purely poloidal fields dynamically unimportant initially, but strong enough to resolve MRI. Magnetic turbulence occurs in the disk, driving accretion and supporting Poynting-dominated jet outflows sufficient to power an SGRB.

The inspiral and merger of a binary neutron star (NSNS) can lead to the formation of a hypermassive neutron star (HMNS). As the HMNS loses thermal pressure due to neutrino cooling and/or centrifugal support due to gravitational wave (GW) emission, and/or magnetic breaking of differential rotation it will collapse to a black hole. To assess the importance of shock-induced thermal pressure and cooling, we adopt an idealized equation of state and perform NSNS simulations in full GR through late inspiral, merger, and HMNS formation, accounting for cooling. We show that thermal pressure contributes significantly to the support of the HMNS against collapse and that thermal cooling accelerates its "delayed" collapse. Our simulations demonstrate explicitly that cooling can induce the catastrophic collapse of a hot hypermassive neutron star formed following the merger of binary neutron stars. Thus, cooling physics is important to include in NSNS merger calculations to accurately determine the lifetime of the HMNS remnant and to extract information about the NS equation of state, cooling mechanisms, bar instabilities and B-fields from the GWs emitted during the transient phase prior to BH formation.

As a neutron star (NS) is tidally disrupted by a black hole (BH) companion at the end of a BH-NS binary inspiral, its magnetic fields will be stretched and amplified. If sufficiently strong, these magnetic fields may impact the gravitational waveforms, merger evolution and mass of the remnant disk. Formation of highly-collimated magnetic field lines in the disk+spinning BH remnant may launch relativistic jets, providing the engine for a short-hard GRB. We analyze this scenario through fully general relativistic, magnetohydrodynamic (GRMHD) BHNS simulations from inspiral through merger and disk formation. Different initial magnetic field configurations and strengths are chosen for the NS interior for both nonspinning and moderately spinning (a/M=0.75) BHs aligned with the orbital angular momentum. Only strong interior (Bmax~10^17 G) initial magnetic fields in the NS significantly influence merger dynamics, enhancing the remnant disk mass by 100% and 40% in the nonspinning and spinning BH cases, respectively. However, detecting the imprint of even a strong magnetic field may be challenging for Advanced LIGO. Though there is no evidence of mass outflows or magnetic field collimation during the preliminary simulations we have performed, higher resolution, coupled with longer disk evolutions and different initial magnetic field configurations, may be required to definitively assess the possibility of BHNS binaries as short-hard GRB progenitors.

We present fully general relativistic (GR) simulations of binary white dwarf-neutron star (WDNS) inspiral and merger. The initial binary is in a circular orbit at the Roche critical separation. The goal is to determine the ultimate fate of such systems. We focus on binaries whose total mass exceeds the maximum mass (Mmax) a cold, degenerate EOS can support against gravitational collapse. The time and length scales span many orders of magnitude, making fully general relativistic hydrodynamic (GRHD) simulations computationally prohibitive. For this reason, we model the WD as a "pseudo-white dwarf" (pWD) as in our binary WDNS head-on collisions study [PRD83:064002,2011]. Our GRHD simulations of a pWDNS system with a 0.98-solar-mass WD and a 1.4-solar-mass NS show that the merger remnant is a spinning Thorne-Zytkow-like Object (TZlO) surrounded by a massive disk. The final total rest mass exceeds Mmax, but the remnant does not collapse promptly. To assess whether the object will ultimately collapse after cooling, we introduce radiative thermal cooling. We first apply our cooling algorithm to TZlOs formed in WDNS head-on collisions, and show that these objects collapse and form black holes on the cooling time scale, as expected. However, when we cool the spinning TZlO formed in the merger of a circular-orbit WDNS binary, the remnant does not collapse, demonstrating that differential rotational support is sufficient to prevent collapse. Given that the final total mass exceeds Mmax, magnetic fields and/or viscosity may redistribute angular momentum and ultimately lead to delayed collapse to a BH. We infer that the merger of realistic massive WDNS binaries likely will lead to the formation of spinning TZlOs that undergo delayed collapse.

We recently developed a new general relativistic magnetohydrodynamic code with adaptive mesh refinement that evolves the electromagnetic (EM) vector potential (A) instead of the magnetic fields directly. Evolving A enables one to use any interpolation scheme on refinement level boundaries and still guarantee that the magnetic field remains divergenceless. As in classical EM, a gauge choice must be made when evolving A, and we chose a straightforward "algebraic" gauge condition to simplify the A evolution equation. However, magnetized black hole-neutron star (BHNS) simulations in this gauge exhibit unphysical behavior, including the spurious appearance of strong magnetic fields on refinement level boundaries. This spurious behavior is exacerbated when matter crosses refinement boundaries during tidal disruption of the NS. Applying Kreiss-Oliger dissipation to the evolution of the magnetic vector potential A slightly weakens this spurious magnetic effect, but with undesired consequences. We demonstrate via an eigenvalue analysis and a numerical study that zero-speed modes in the algebraic gauge, coupled with the frequency filtering that occurs on refinement level boundaries, are responsible for the creation of spurious magnetic fields. We show that the EM Lorenz gauge exhibits no zero-speed modes, and as a consequence, spurious magnetic effects are quickly propagated away, allowing for long-term, stable magnetized BHNS evolutions. Our study demonstrates how the EM gauge degree of freedom can be chosen to one's advantage, and that for magnetized BHNS simulations the Lorenz gauge constitutes a major improvement over the algebraic gauge.

Theories of gravity other than general relativity (GR) can explain the observed cosmic acceleration without a cosmological constant. One such class of theories of gravity is f(R). Metric f(R) theories have been proven to be equivalent to Brans-Dicke (BD) scalar-tensor gravity without a kinetic term. Using this equivalence and a 3+1 decomposition of the theory it has been shown that metric f(R) gravity admits a well-posed initial value problem. However, it has not been proven that the 3+1 evolution equations of metric f(R) gravity preserve the (hamiltonian and momentum) constraints. In this paper we show that this is indeed the case. In addition, we show that the mathematical form of the constraint propagation equations in BD-equilavent f(R) gravity and in f(R) gravity in both the Jordan and Einstein frames, is exactly the same as in the standard ADM 3+1 decomposition of GR. Finally, we point out that current numerical relativity codes can incorporate the 3+1 evolution equations of metric f(R) gravity by modifying the stress-energy tensor and adding an additional scalar field evolution equation. We hope that this work will serve as a starting point for relativists to develop fully dynamical codes for valid f(R) models.

We simulate head-on collisions from rest at large separation of binary white dwarf -- neutron stars (WDNSs) in full general relativity. Our study serves as a prelude to our analysis of the circular binary WDNS problem. We focus on compact binaries whose total mass exceeds the maximum mass that a cold degenerate star can support, and our goal is to determine the fate of such systems. A fully general relativistic hydrodynamic computation of a realistic WDNS head-on collision is prohibitive due to the large range of dynamical time scales and length scales involved. For this reason, we construct an equation of state (EOS) which captures the main physical features of NSs while, at the same time, scales down the size of WDs. We call these scaled-down WD models "pseudo-WDs (pWDs)". Using pWDs, we can study these systems via a sequence of simulations where the size of the pWD gradually increases toward the realistic case. We perform two sets of simulations; One set studies the effects of the NS mass on the final outcome, when the pWD is kept fixed. The other set studies the effect of the pWD compaction on the final outcome, when the pWD mass and the NS are kept fixed. All simulations show that 14%-18% of the initial total rest mass escapes to infinity. All remnant masses still exceed the maximum rest mass that our cold EOS can support (1.92 solar masses), but no case leads to prompt collapse to a black hole. This outcome arises because the final configurations are hot. All cases settle into spherical, quasiequilibrium configurations consisting of a cold NS core surrounded by a hot mantle, resembling Thorne-Zytkow objects. Extrapolating our results to realistic WD compactions, we predict that the likely outcome of a head-on collision of a realistic, massive WDNS system will be the formation of a quasiequilibrium Thorne-Zytkow-like object.

White dwarf-neutron star binaries generate detectable gravitational radiation. We construct Newtonian equilibrium models of corotational white dwarf-neutron star (WDNS) binaries in circular orbit and find that these models terminate at the Roche limit. At this point the binary will undergo either stable mass transfer (SMT) and evolve on a secular time scale, or unstable mass transfer (UMT), which results in the tidal disruption of the WD. The path a given binary will follow depends primarily on its mass ratio. We analyze the fate of known WDNS binaries and use population synthesis results to estimate the number of LISA-resolved galactic binaries that will undergo either SMT or UMT. We model the quasistationary SMT epoch by solving a set of simple ordinary differential equations and compute the corresponding gravitational waveforms. Finally, we discuss in general terms the possible fate of binaries that undergo UMT and construct approximate Newtonian equilibrium configurations of merged WDNS remnants. We use these configurations to assess plausible outcomes of our future, fully relativistic simulations of these systems. If sufficient WD debris lands on the NS, the remnant may collapse, whereby the gravitational waves from the inspiral, merger, and collapse phases will sweep from LISA through LIGO frequency bands. If the debris forms a disk about the NS, it may fragment and form planets.

We present an analysis of well-posedness of constrained evolution of 3+1 formulations of GR. In this analysis we explicitly take into account the energy and momentum constraints as well as possible algebraic constraints on the evolution of high-frequency perturbations of solutions of Einstein's equations. In this respect, our approach is principally different from standard analyses of well-posedness of free evolution in general relativity. Our study reveals the existence of subsets of the linearized Einstein's equations that control the well-posedness of constrained evolution. It is demonstrated that the well-posedness of ADM, BSSN and other 3+1 formulations derived from ADM by adding combinations of constraints to the right-hand-side of ADM and/or by linear transformation of the dynamical ADM variables depends entirely on the properties of the gauge. For certain classes of gauges we formulate conditions for well-posedness of constrained evolution. This provides a new basis for constructing stable numerical integration schemes for a classical Arnowitt--Deser--Misner (ADM) and many other 3+1 formulations of general relativity.