Joint multi-messenger observations with gravitational waves and electromagnetic data offer new insights into the astrophysical studies of compact objects. The third Advanced LIGO and Advanced Virgo observing run began on April 1, 2019; during the eleven months of observation, there have been 14 compact binary systems candidates for which at least one component is potentially a neutron star. Although intensive follow-up campaigns involving tens of ground and space-based observatories searched for counterparts, no electromagnetic counterpart has been detected. Following on a previous study of the first six months of the campaign, we present in this paper the next five months of the campaign from October 2019 to March 2020. We highlight two neutron star - black hole candidates (S191205ah, S200105ae), two binary neutron star candidates (S191213g and S200213t) and a binary merger with a possible neutron star and a "MassGap" component, S200115j. Assuming that the gravitational-wave candidates are of astrophysical origin and their location was covered by optical telescopes, we derive possible constraints on the matter ejected during the events based on the non-detection of counterparts. We find that the follow-up observations during the second half of the third observing run did not meet the necessary sensitivity to constrain the source properties of the potential gravitational-wave candidate. Consequently, we suggest that different strategies have to be used to allow a better usage of the available telescope time. We examine different choices for follow-up surveys to optimize sky localization coverage vs.\ observational depth to understand the likelihood of counterpart detection.

Of the three main types of binaries detectable through ground-based gravitational wave observations, black hole-neutron star (BHNS) mergers remain the most elusive. While candidates BHNS exist in the triggers released during the third observing run of the Advanced LIGO/Virgo collaboration, no detection has been confirmed so far. As for binary neutron star systems, BHNS binaries allow us to explore a wide range of physical processes, including the neutron star equation of state, nucleosynthesis, stellar evolution, high-energy astrophysics, and the expansion of the Universe. Here, we review some of the main features of BHNS systems: the distinction between disrupting and non-disrupting binaries, the types of outflows that BHNS mergers can produce, and the information that can be extracted from the observation of their gravitational wave and electromagnetic signals. We also emphasize that for the most likely binary parameters, BHNS mergers seem less likely to power electromagnetic signals than binary neutron star systems. Finally, we discuss some of the issues that still limit our ability to model and interpret electromagnetic signals from BHNS binaries.

We investigate mass ejection from accretion disks formed in mergers of black holes (BHs) and neutron stars (NSs). The third observing run of the LIGO/Virgo interferometers provided BH-NS candidate events that yielded no electromagnetic (EM) counterparts. The broad range of disk configurations expected from BH-NS mergers motivates a thorough exploration of parameter space to improve EM signal predictions. Here we conduct 27 high-resolution, axisymmetric, long-term hydrodynamic simulations of the viscous evolution of black hole accretion disks that include neutrino emission/absorption effects and post-processing with a nuclear reaction network. In the absence of magnetic fields, these simulations provide a lower-limit to the fraction of the initial disk mass ejected. We find a nearly linear inverse dependence of this fraction on disk compactness (BH mass over initial disk radius). The dependence is related to the fraction of the disk mass accreted at the time when weak interactions freeze out. We also characterize a trend of decreasing ejected fraction and decreasing lanthanide/actinide content with increasing disk mass at fixed BH mass. This trend results from a longer time to reach weak freezout and an increasingly dominant role of neutrino absorption at higher disk masses. We provide estimates of the radioactive luminosity from the disk outflow alone available to power kilonovae over the range of BH and NS combinations studied, finding a spread of two orders of magnitude. For most of the BH-NS parameter space, the disk outflow contribution is well below the kilonova mass upper limits for S190814bv.

There is irresistible observational evidence that binary systems of compact objects with at least one neutron star are progenitors of short gamma-ray bursts, as well as a production site for r-process elements, at least when some matter is ejected by the merger and an accretion disk is formed. The recent observations of gravitational waves in conjunction with electromagnetic counterparts fuel the need for models predicting the outcome of a given merger and the properties of the associated matter outflows as a function of the initial parameters of the binary. In this manuscript, we provide updated fitting formulae that estimate the disk mass for double neutron star binaries and ejecta masses for black hole-neutron star and double neutron star binaries, fitted to the results of numerical simulations. Our proposed fitting formulae improve on existing models by aiming for analytical simplicity, by covering a larger region of parameter space, and by accounting for regions of parameter space not covered by numerical simulations but with physically manifest merger outcomes.

The discovery of GW170817 with gravitational waves (GWs) and electromagnetic (EM) radiation is prompting new questions in strong-gravity astrophysics. Importantly, it remains unknown whether the progenitor of the merger comprised two neutron stars (NSs), or a NS and a black hole (BH). Using new numerical-relativity simulations and incorporating modeling uncertainties we produce novel GW and EM observables for NS-BH mergers with similar masses. A joint analysis of GW and EM measurements reveals that if GW170817 is a NS-BH merger, <40% of the binary parameters consistent with the GW data are compatible with EM observations.

On 2019 August 14, the Advanced LIGO and Virgo interferometers detected the high-significance gravitational wave (GW) signal S190814bv. The GW data indicated that the event resulted from a neutron star--black hole (NSBH) merger, or potentially a low-mass binary black hole merger. Due to the low false alarm rate and the precise localization (23 deg$^2$ at 90\%), S190814bv presented the community with the best opportunity yet to directly observe an optical/near-infrared counterpart to a NSBH merger. To search for potential counterparts, the GROWTH collaboration performed real-time image subtraction on 6 nights of public Dark Energy Camera (DECam) images acquired in the three weeks following the merger, covering $>$98\% of the localization probability. Using a worldwide network of follow-up facilities, we systematically undertook spectroscopy and imaging of optical counterpart candidates. Combining these data with a photometric redshift catalog, we ruled out each candidate as the counterpart to S190814bv and we placed deep, uniform limits on the optical emission associated with S190814bv. For the nearest consistent GW distance, radiative transfer simulations of NSBH mergers constrain the ejecta mass of S190814bv to be $M_\mathrm{ej} < 0.04$~$M_{\odot}$ at polar viewing angles, or $M_\mathrm{ej} < 0.03$~$M_{\odot}$ if the opacity is $\kappa < 2$~cm$^2$g$^{-1}$. Assuming a tidal deformability for the neutron star at the high end of the range compatible with GW170817 results, our limits would constrain the BH spin component aligned with the orbital momentum to be $ \chi < 0.7$ for mass ratios $Q < 6$, with weaker constraints for more compact neutron stars. We publicly release the photometry from this campaign at http://www.astro.caltech.edu/~danny/static/s190814bv.

GW170817 showed that neutron star mergers not only emit gravitational waves but also can release electromagnetic signatures in multiple wavelengths. Within the first half of the third observing run of the Advanced LIGO and Virgo detectors, there have been a number of gravitational wave candidates of compact binary systems for which at least one component is potentially a neutron star. In this article, we look at the candidates S190425z, S190426c, S190510g, S190901ap, and S190910h, predicted to have potentially a non-zero remnant mass, in more detail. All these triggers have been followed up with extensive campaigns by the astronomical community doing electromagnetic searches for their optical counterparts; however, according to the released classification, there is a high probability that some of these events might not be of extraterrestrial origin. Assuming that the triggers are caused by a compact binary coalescence and that the individual source locations have been covered during the EM follow-up campaigns, we employ three different kilonova models and apply them to derive possible constraints on the matter ejection consistent with the publicly available gravitational-wave trigger information and the lack of a kilonova detection. These upper bounds on the ejecta mass can be related to limits on the maximum mass of the binary neutron star candidate S190425z and to constraints on the mass-ratio, spin, and NS compactness for the potential black hole-neutron star candidate S190426c. Our results show that deeper electromagnetic observations for future gravitational wave events near the horizon limit of the advanced detectors are essential.

High-accuracy numerical simulations of merging neutron stars play an important role in testing and calibrating the waveform models used by gravitational wave observatories. Obtaining high-accuracy waveforms at a reasonable computational cost, however, remains a significant challenge. One issue is that high-order convergence of the solution requires the use of smooth evolution variables, while many of the equations of state used to model the neutron star matter have discontinuities, typically in the first derivative of the pressure. Spectral formulations of the equation of state have been proposed as a potential solution to this problem. Here, we report on the numerical implementation of spectral equations of state in the Spectral Einstein Code. We show that, in our code, spectral equations of state allow for high-accuracy simulations at a lower computational cost than commonly used `piecewise polytrope' equations state. We also demonstrate that not all spectral equations of state are equally useful: different choices for the low-density part of the equation of state can significantly impact the cost and accuracy of simulations. As a result, simulations of neutron star mergers present us with a trade-off between the cost of simulations and the physical realism of the chosen equation of state.

Neutron star mergers are unique laboratories of accretion, ejection, and r-process nucleosynthesis. We used 3D general relativistic magnetohydrodynamic simulations to study the role of the post-merger magnetic geometry in the evolution of merger remnant discs around stationary Kerr black holes. Our simulations fully capture mass accretion, ejection, and jet production, owing to their exceptionally long duration exceeding $4$ s. Poloidal post-merger magnetic field configurations produce jets with energies $E_\mathrm{jet} \sim (4{-}30)\times10^{50}$ erg, isotropic equivalent energies $E_\mathrm{iso}\sim(4{-}20)\times10^{52}$ erg, opening angles $\theta_\mathrm{jet}\sim6{-}13^\circ$, and durations $t_j\lesssim1$ s. Accompanying the production of jets is the ejection of $f_\mathrm{ej}\sim30{-}40\%$ of the post-merger disc mass, continuing out to times $> 1$ s. We discover that a more natural, purely toroidal post-merger magnetic field geometry generates large-scale poloidal magnetic flux of alternating polarity and striped jets. The first stripe, of $E_\mathrm{jet}\simeq2\times10^{48}\,\mathrm{erg}$, $E_\mathrm{iso}\sim10^{51}$ erg, $\theta_\mathrm{jet}\sim3.5{-}5^\circ$, and $t_j\sim0.1$ s, is followed by $\gtrsim4$ s of striped jet activity with $f_\mathrm{ej}\simeq27\%$. The dissipation of such stripes could power the short gamma-ray burst (sGRB) prompt emission. Our simulated jet energies and durations span the range of sGRBs. We find that although the blue kilonova component is initially hidden from view by the red component, it expands faster, outruns the red component, and becomes visible to off-axis observers. In comparison to GW 170817/GRB 170817A, our simulations under-predict the mass of the blue relative to red component by a factor of few. Including the dynamical ejecta and neutrino absorption may reduce this tension.

Simulations of neutron star-black hole (NSBH) binaries generally consider black holes with masses in the range $(5-10)M_\odot$, where we expect to find most stellar mass black holes. The existence of lower mass black holes, however, cannot be theoretically ruled out. Low-mass black holes in binary systems with a neutron star companion could mimic neutron star-neutron (NSNS) binaries, as they power similar gravitational wave (GW) and electromagnetic (EM) signals. To understand the differences and similarities between NSNS mergers and low-mass NSBH mergers, numerical simulations are required. Here, we perform a set of simulations of low-mass NSBH mergers, including systems compatible with GW170817. Our simulations use a composition and temperature dependent equation of state (DD2) and approximate neutrino transport, but no magnetic fields. We find that low-mass NSBH mergers produce remnant disks significantly less massive than previously expected, and consistent with the post-merger outflow mass inferred from GW170817 for moderately asymmetric mass ratio. The dynamical ejecta produced by systems compatible with GW170817 is negligible except if the mass ratio and black hole spin are at the edge of the allowed parameter space. That dynamical ejecta is cold, neutron-rich, and surprisingly slow for ejecta produced during the tidal disruption of a neutron star : $v\sim (0.1-0.15)c$. We also find that the final mass of the remnant black hole is consistent with existing analytical predictions, while the final spin of that black hole is noticeably larger than expected -- up to $\chi_{\rm BH}=0.84$ for our equal mass case.

Gamma-ray astronomy explores the most energetic photons in nature to address some of the most pressing puzzles in contemporary astrophysics. It encompasses a wide range of objects and phenomena: stars, supernovae, novae, neutron stars, stellar-mass black holes, nucleosynthesis, the interstellar medium, cosmic rays and relativistic-particle acceleration, and the evolution of galaxies. MeV gamma-rays provide a unique probe of nuclear processes in astronomy, directly measuring radioactive decay, nuclear de-excitation, and positron annihilation. The substantial information carried by gamma-ray photons allows us to see deeper into these objects, the bulk of the power is often emitted at gamma-ray energies, and radioactivity provides a natural physical clock that adds unique information. New science will be driven by time-domain population studies at gamma-ray energies. This science is enabled by next-generation gamma-ray instruments with one to two orders of magnitude better sensitivity, larger sky coverage, and faster cadence than all previous gamma-ray instruments. This transformative capability permits: (a) the accurate identification of the gamma-ray emitting objects and correlations with observations taken at other wavelengths and with other messengers; (b) construction of new gamma-ray maps of the Milky Way and other nearby galaxies where extended regions are distinguished from point sources; and (c) considerable serendipitous science of scarce events -- nearby neutron star mergers, for example. Advances in technology push the performance of new gamma-ray instruments to address a wide set of astrophysical questions.

Gravitational waveforms from numerical simulations are a critical tool to test and analytically calibrate the waveform models used to study the properties of merging compact objects. In this paper, we present a series of high-accuracy waveforms produced with the SpEC code for systems involving at least one neutron star. We provide for the first time waveforms with sub-radian accuracy over more than twenty cycles for low-mass black hole-neutron star binaries, including binaries with non-spinning objects, and binaries with rapidly spinning neutron stars that maximize the impact on the gravitational wave signal of the near-resonant growth of the fundamental excitation mode of the neutron star (f-mode). We also provide for the first time with SpEC a high-accuracy neutron star-neutron star waveform. These waveforms are made publicly available as part of the SxS catalogue. We compare our results to analytical waveform models currently implemented in data analysis pipelines. For most simulations, the models lie outside of the predicted numerical errors in the last few orbits before merger, but do not show systematic deviations from the numerical results: comparing different models appears to provide reasonable estimates of the modeling errors. The sole exception is the equal-mass simulation using a rapidly counter-rotating neutron star to maximize the impact of the excitation of the f-mode, for which all models perform poorly. This is however expected, as even the single model that takes f-mode excitation into account ignores the significant impact of the neutron star spin on the f-mode excitation frequency.

Mergers of compact binaries containing two neutron stars (NS-NS), or a neutron star and a stellar-mass black hole (NS-BH), are likely progenitors of short-duration gamma ray bursts (SGRBs). A fraction >20% of SGRBs are followed by temporally-extended (>~ minute-long), variable X-ray emission, attributed to ongoing activity of the central engine. One source of late-time engine activity is fall-back accretion of bound tidal ejecta; however, observed extended emission light curves do not track the naively-anticipated, uninterrupted ~t^{-5/3} power-law decay, instead showing a lull or gap in emission typically lasting tens of seconds after the burst. Here, we re-examine the impact of heating due to rapid neutron capture (r-process) nucleosynthesis on rate of the fall-back accretion, using ejecta properties obtained from numerical relativity simulations of NS-BH mergers in a toy model for the dynamical influence of nuclear heating. Depending on the mass of the remnant black hole, r-process heating can imprint a variety of fall-back curve shapes, ranging from temporal lulls of up to tens of seconds to complete late-time cut-off in the fall-back rate. This behavior is robust to realistic variations in the nuclear heating experienced by different parts of the ejecta. Central black holes with masses <~ 3 Msun typically experience absolute cut-offs in the fall-back rate, while more massive >6-8 Msun black holes instead show temporal gaps. We thus propose that SGRBs with extended emission arise from NS-BH, rather than NS-NS, mergers. Our model implies a NS-BH merger detection rate by LIGO which, in steady-state, is comparable to or greater than that of NS-NS mergers.

We investigate the long-term evolution of black hole accretion disks formed in neutron star mergers. These disks expel matter that contributes to an $r$-process kilonova, and can produce relativistic jets powering short gamma-ray bursts. Here we report the results of a three-dimensional, general-relativistic magnetohydrodynamic (GRMHD) simulation of such a disk which is evolved for long enough ($\sim 9$s, or $\sim 6\times 10^5 r_{\rm g}/c$) to achieve completion of mass ejection far from the disk. Our model starts with a poloidal field, and fully resolves the most unstable mode of the magnetorotational instability. We parameterize the dominant microphysics and neutrino cooling effects, and compare with axisymmetric hydrodynamic models with shear viscosity. The GRMHD model ejects mass in two ways: a prompt MHD-mediated outflow and a late-time, thermally-driven wind once the disk becomes advective. The total amount of unbound mass ejected ($0.013M_\odot$, or $\simeq 40\%$ of the initial torus mass) is twice as much as in hydrodynamic models, with higher average velocity ($0.1c$) and a broad electron fraction distribution with a lower average value ($0.16$). Scaling the ejected fractions to a disk mass of $\sim 0.1M_\odot$ can account for the red kilonova from GW170817 but underpredicts the blue component. About $\sim 10^{-3}M_\odot$ of material should undergo neutron freezout and could produce a bright kilonova precursor in the first few hours after the merger. With our idealized initial magnetic field configuration, we obtain a robust jet and sufficient ejecta with Lorentz factor $\sim 1-10$ to (over)produce the non-thermal emission from GW1708107.

We identify various contributors of systematic effects in the measurement of the neutron star (NS) tidal deformability and quantify their magnitude for several types of neutron star - black hole (NSBH) binaries. Gravitational waves from NSBH mergers contain information about the components' masses and spins as well as the NS equation of state. Extracting this information requires comparison of the signal in noisy detector data with theoretical templates derived from some combination of post-Newtonian (PN) approximants, effective one-body (EOB) models and %analytic fits to numerical relativity (NR) simulations. The accuracy of these templates is limited by errors in the NR simulations, by the approximate nature of the PN/EOB waveforms, and by the hybridization procedure used to combine them. In this paper, we estimate the impact of these errors by constructing and comparing a set of PN-NR hybrid waveforms, for the first time with NR waveforms from two different codes, namely, SpEC and SACRA, for such systems. We then attempt to recover the parameters of the binary using two non-precessing template approximants. We find that systematic errors are too large for tidal effects to be accurately characterized for any realistic NS equation of state model. We conclude that NSBH waveform models must be significantly improved if they are to be useful for the extraction of NS equation of state information or even for distinguishing NSBH systems from binary black holes.

Gravitational-wave (GW) and electromagnetic (EM) signals from the merger of a Neutron Star (NS) and a Black Hole (BH) are a highly anticipated discovery in extreme gravity, nuclear-, and astrophysics. We develop a simple formula that distinguishes between merger outcomes and predicts the post-merger remnant mass, validated with 75 simulations. Our formula improves on existing results by describing critical unexplored regimes: comparable masses and higher BH spins. These are important to differentiate NSNS from NSBH mergers, and to infer source physics from EM signals.

We present a covariant ray tracing algorithm for computing high-resolution neutrino distributions in general relativistic numerical spacetimes with hydrodynamical sources. Our formulation treats the very important effect of elastic scattering of neutrinos off of nuclei and nucleons (changing the neutrino's direction but not energy) by incorporating estimates of the background neutrino fields. Background fields provide information about the spectra and intensities of the neutrinos scattered into each ray. These background fields may be taken from a low-order moment simulation or be ignored, in which case the method reduces to a standard state-of-the-art ray tracing formulation. The method handles radiation in regimes spanning optically thick to optically thin. We test the new code, highlight its strengths and weaknesses, and apply it to a simulation of a neutron star merger to compute neutrino fluxes and spectra, and to demonstrate a neutrino flavor oscillation calculation. In that environment, we find qualitatively different fluxes, spectra, and oscillation behaviors when elastic scattering is included.

Neutrino-matter interactions play an important role in the post-merger evolution of neutron star-neutron star and black hole-neutron star mergers. Most notably, they determine the properties of the bright optical/infrared transients observable after a merger. Unfortunately, Boltzmann's equations of radiation transport remain too costly to be evolved directly in merger simulations. Simulations rely instead on approximate transport algorithms with unquantified modeling errors. In this paper, we use for the first time a time-dependent general relativistic Monte-Carlo (MC) algorithm to solve Boltzmann's equations and estimate important properties of the neutrino distribution function ~10ms after a neutron star merger. We do not fully couple the MC algorithm to the fluid evolution, but use a short evolution of the merger remnant to critically assess errors in our approximate gray two-moment transport scheme. We demonstrate that the analytical closure used by the moment scheme is highly inaccurate in the polar regions, but performs well elsewhere. While the average energy of polar neutrinos is reasonably well captured by the two-moment scheme, estimates for the neutrino energy become less accurate at lower latitudes. The two-moment formalism also overestimates the density of neutrinos in the polar regions by ~50%, and underestimates the neutrino pair-annihilation rate at the poles by factors of 2-3. Although the latter is significantly more accurate than one might have expected before this study, our results indicate that predictions for the properties of polar outflows and for the creation of a baryon-free region at the poles are likely to be affected by errors in the two-moment scheme, thus limiting our ability to reliably model kilonovae and gamma-ray bursts.

This is an exciting time for the study of r-process nucleosynthesis. Recently, a neutron star merger GW170817 was observed in extraordinary detail with gravitational waves and electromagnetic radiation from radio to gamma rays. The very red color of the associated kilonova suggests that neutron star mergers are an important r-process site. Astrophysical simulations of neutron star mergers and core collapse supernovae are making rapid progress. Detection of both, electron neutrinos and antineutrinos from the next galactic supernova will constrain the composition of neutrino-driven winds and provide unique nucleosynthesis information. Finally FRIB and other rare-isotope beam facilities will soon have dramatic new capabilities to synthesize many neutron-rich nuclei that are involved in the r-process. The new capabilities can significantly improve our understanding of the r-process and likely resolve one of the main outstanding problems in classical nuclear astrophysics. However, to make best use of the new experimental capabilities and to fully interpret the results, a great deal of infrastructure is needed in many related areas of astrophysics, astronomy, and nuclear theory. We will place these experiments in context by discussing astrophysical simulations and observations of r-process sites, observations of stellar abundances, galactic chemical evolution, and nuclear theory for the structure and reactions of very neutron-rich nuclei. This review paper was initiated at a three-week International Collaborations in Nuclear Theory program in June 2016 where we explored promising r-process experiments and discussed their likely impact, and their astrophysical, astronomical, and nuclear theory context.

Each of the potential signals from a black hole-neutron star merger should contain an imprint of the neutron star equation of state: gravitational waves via its effect on tidal disruption, the kilonova via its effect on the ejecta, and the gamma ray burst via its effect on the remnant disk. These effects have been studied by numerical simulations and quantified by semi-analytic formulae. However, most of the simulations on which these formulae are based use equations of state without finite temperature and composition-dependent nuclear physics. In this paper, we simulate black hole-neutron star mergers varying both the neutron star mass and the equation of state, using three finite-temperature nuclear models of varying stiffness. Our simulations largely vindicate formulae for ejecta properties but do not find the expected dependence of disk mass on neutron star compaction. We track the early evolution of the accretion disk, largely driven by shocking and fallback inflow, and do find notable equation of state effects on the structure of this early-time, neutrino-bright disk.

Turbulence in an accretion disk launches Alfv\'en waves (AWs) that propagate away from the disk along magnetic field lines. Because the Alfv\'en speed varies with distance from the disk, the AWs undergo partial non-WKB reflection, and counter-propagating AWs subsequently interact, causing AW energy to cascade to small scales and dissipate. To investigate this process, we introduce an Elsasser-like formulation of general relativistic MHD (GRMHD) and develop the theory of general relativistic reduced MHD in an inhomogeneous medium. We then derive a set of equations for the mean-square AW amplitude $M_+$ and turbulent heating rate $Q$ under the assumption that, in the plasma rest frame, AWs propagating away from the disk are much more energetic than AWs propagating toward the disk. For the case in which the background flow is axisymmetric and time-independent, we solve these equations analytically to determine $M_+$ and $Q$ as functions of position. We find that, for an idealized thin disk threaded by a large-scale poloidal magnetic field, the AW energy flux is $\sim (\rho_{\rm b}/\rho_{\rm d})^{1/2} \beta_{\rm net,d}^{-1/2}$ times the disk's radiative flux, where $\rho_{\rm b}$ and $\rho_{\rm d}$ are the mass densities at the coronal base and disk midplane, respectively, and $\beta_{\rm net,d}$ is the ratio (evaluated at the disk midplane) of plasma-plus-radiation pressure to the pressure of the average vertical magnetic field. This energy flux could have a significant impact on disk coronae and outflows. To lay the groundwork for future global simulations of turbulent disk coronae and jets, we derive a set of averaged GRMHD equations that account for reflection-driven AW turbulence using a sub-grid model.

General relativistic radiation hydrodynamics simulations are necessary to accurately model a number of astrophysical systems involving black holes and neutron stars. Photon transport plays a crucial role in radiatively dominated accretion disks, while neutrino transport is critical to core-collapse supernovae and to the modeling of electromagnetic transients and nucleosynthesis in neutron star mergers. However, evolving the full Boltzmann equations of radiative transport is extremely expensive. Here, we describe the implementation in the general relativistic SpEC code of a cheaper radiation hydrodynamics method which theoretically converges to a solution of Boltzmann's equation in the limit of infinite numerical resources. The algorithm is based on a gray two-moment scheme, in which we evolve the energy density and momentum density of the radiation. Two-moment schemes require a closure which fills in missing information about the energy spectrum and higher-order moments of the radiation. Instead of the approximate analytical closure currently used in core-collapse and merger simulations, we complement the two-moment scheme with a low-accuracy Monte-Carlo evolution. The Monte-Carlo results can provide any or all of the missing information in the evolution of the moments, as desired by the user. As a first test of our methods, we study a set of idealized problems demonstrating that our algorithm performs significantly better than existing analytical closures. We also discuss the current limitations of our method, in particular open questions regarding the stability of the fully coupled scheme.

Black hole-torus systems from compact binary mergers are possible engines for gamma-ray bursts (GRBs). During the early evolution of the post-merger remnant, the state of the torus is determined by a combination of neutrino cooling and magnetically-driven heating processes, so realistic models must include both effects. In this paper, we study the post-merger evolution of a magnetized black hole-neutron star binary system using the Spectral Einstein Code (SpEC) from an initial post-merger state provided by previous numerical relativity simulations. We use a finite-temperature nuclear equation of state and incorporate neutrino effects in a leakage approximation. To achieve the needed accuracy, we introduce improvements to SpEC's implementation of general-relativistic magnetohydrodynamics (MHD), including the use of cubed-sphere multipatch grids and an improved method for dealing with supersonic accretion flows where primitive variable recovery is difficult. We find that a seed magnetic field triggers a sustained source of heating, but its thermal effects are largely cancelled by the accretion and spreading of the torus from MHD-related angular momentum transport. The neutrino luminosity peaks at the start of the simulation, and then drops significantly over the first 20\,ms but in roughly the same way for magnetized and nonmagnetized disks. The heating rate and disk's luminosity decrease much more slowly thereafter. These features of the evolution are insensitive to grid structure and resolution, formulation of the MHD equations, and seed field strength, although turbulent effects are not fully converged

We investigate the nucleosynthesis of heavy elements in the winds ejected by accretion disks formed in neutron star mergers. We compute the element formation in disk outflows from hypermassive neutron star (HMNS) remnants of variable lifetime, including the effect of angular momentum transport in the disk evolution. We employ long-term axisymmetric hydrodynamic disk simulations to model the ejecta, and compute r-process nucleosynthesis with tracer particles using a nuclear reaction network containing $\sim 8000$ species. We find that the previously known strong correlation between HMNS lifetime, ejected mass, and average electron fraction in the outflow is directly related to the amount of neutrino irradiation on the disk, which dominates mass ejection at early times in the form of a neutrino-driven wind. Production of lanthanides and actinides saturates at short HMNS lifetimes ($\lesssim 10$ ms), with additional ejecta contributing to a blue optical kilonova component for longer-lived HMNSs. We find good agreement between the abundances from the disk outflow alone and the solar r-process distribution only for short HMNS lifetimes ($\lesssim 10$ ms). For longer lifetimes, the rare-earth and third r-process peaks are significantly under-produced compared to the solar pattern, requiring additional contributions from the dynamical ejecta. The nucleosynthesis signature from a spinning black hole (BH) can only overlap with that from a HMNS of moderate lifetime ($\lesssim 60$ ms). Finally, we show that angular momentum transport not only contributes with a late-time outflow component, but that it also enhances the neutrino-driven component by moving material to shallower regions of the gravitational potential, in addition to providing additional heating.

We introduce a new relativistic astrophysics code, SpECTRE, that combines a discontinuous Galerkin method with a task-based parallelism model. SpECTRE's goal is to achieve more accurate solutions for challenging relativistic astrophysics problems such as core-collapse supernovae and binary neutron star mergers. The robustness of the discontinuous Galerkin method allows for the use of high-resolution shock capturing methods in regions where (relativistic) shocks are found, while exploiting high-order accuracy in smooth regions. A task-based parallelism model allows efficient use of the largest supercomputers for problems with a heterogeneous workload over disparate spatial and temporal scales. We argue that the locality and algorithmic structure of discontinuous Galerkin methods will exhibit good scalability within a task-based parallelism framework. We demonstrate the code on a wide variety of challenging benchmark problems in (non)-relativistic (magneto)-hydrodynamics. We demonstrate the code's scalability including its strong scaling on the NCSA Blue Waters supercomputer up to the machine's full capacity of 22,380 nodes using 671,400 threads.

We investigate the ejecta from black hole - neutron star mergers by modeling the formation and interaction of mass ejected in a tidal tail and a disk wind. The outflows are neutron-rich, giving rise to optical/infrared emission powered by the radioactive decay of $r$-process elements (a kilonova). Here we perform an end-to-end study of this phenomenon, where we start from the output of a fully-relativistic merger simulation, calculate the post-merger hydrodynamical evolution of the ejecta and disk winds including neutrino physics, determine the final nucleosynthetic yields using post-processing nuclear reaction network calculations, and compute the kilonova emission with a radiative transfer code. We study the effects of the tail-to-disk mass ratio by scaling the tail density. A larger initial tail mass results in fallback matter becoming mixed into the disk and ejected in the subsequent disk wind. Relative to the case of a disk without dynamical ejecta, the combined outflow has lower mean electron fraction, faster speed, larger total mass, and larger absolute mass free of high-opacity Lanthanides or Actinides. In most cases, the nucleosynthetic yield is dominated by the heavy $r$-process contribution from the unbound part of the tidal tail. A Solar-like abundance distribution can however be obtained when the total mass of the dynamical ejecta is comparable to the mass of the disk outflows. The kilonova has a characteristic duration of 1 week and a luminosity of ~$10^{41}$ erg/s, with orientation effects leading to variations of a factor ~2 in brightness. At early times (< 1 day) the emission includes an optical component from the (hot) Lanthanide-rich material, but the spectrum evolves quickly to the infrared thereafter.

Black holes with accretion rates well below the Eddington rate are expected to be surrounded by low-density, hot, geometrically thick accretion disks. This includes the two black holes being imaged at sub-horizon resolution by the Event Horizon Telescope. In these disks, the mean free path for Coulomb interactions between charged particles is large, and the accreting matter is a nearly collisionless plasma. Despite this, numerical simulations have so far modeled these accretion flows using ideal magnetohydrodynamics. Here, we present the first global, general relativistic, 3D simulations of accretion flows onto a Kerr black hole including the non-ideal effects most likely to affect the dynamics of the disk: the anisotropy between the pressure parallel and perpendicular to the magnetic field, and the heat flux along magnetic field lines. We show that for both standard and magnetically arrested disks, the pressure anisotropy is comparable to the magnetic pressure, while the heat flux remains dynamically unimportant. Despite this large pressure anisotropy, however, the time-averaged structure of the accretion flow is strikingly similar to that found in simulations treating the plasma as an ideal fluid. We argue that these similarities are largely due to the interchangeability of the viscous and magnetic shear stresses as long as the magnetic pressure is small compared to the gas pressure, and to the sub-dominant role of pressure/viscous effects in magnetically arrested disks. We conclude by highlighting outstanding questions in modeling the dynamics of low collisionality accretion flows.

During the merger of a black hole and a neutron star, baryonic mass can become unbound from the system. Because the ejected material is extremely neutron-rich, the r-process rapidly synthesizes heavy nuclides as the material expands and cools. In this work, we map general relativistic models of black hole-neutron star (BHNS) mergers into a Newtonian smoothed particle hydrodynamics (SPH) code and follow the evolution of the thermodynamics and morphology of the ejecta until the outflows become homologous. We investigate how the subsequent evolution depends on our mapping procedure and find that the results are robust. Using thermodynamic histories from the SPH particles, we then calculate the expected nucleosynthesis in these outflows while varying the level of neutrino irradiation coming from the postmerger accretion disk. We find that the ejected material robustly produces r-process nucleosynthesis even for unrealistically high neutrino luminosities, due to the rapid velocities of the outflow. Nonetheless, we find that neutrinos can have an impact on the detailed pat- tern of the r-process nucleosynthesis. Electron neutrinos are captured by neutrons to produce protons while neutron capture is occurring. The produced protons rapidly form low mass seed nuclei for the r-process. These low mass seeds are eventually incorporated into the first r-process peak at A~78, producing mainly Ge and Se. We consider the mechanism of this process in detail and discuss if it can impact galactic chemical evolution of the first peak r-process nuclei.

Hot, diffuse, relativistic plasmas such as sub-Eddington black hole accretion flows are expected to be collisionless, yet are commonly modeled as a fluid using ideal general relativistic magnetohydrodynamics (GRMHD). Dissipative effects such as heat conduction and viscosity can be important in a collisionless plasma and will potentially alter the dynamics and radiative properties of the flow from that in ideal fluid models; we refer to models that include these processes as Extended GRMHD. Here we describe a new conservative code, grim, that enables all the above and additional physics to be efficiently incorporated. grim combines time evolution and primitive variable inversion needed for conservative schemes into a single step using an algorithm that only requires the residuals of the governing equations as inputs. This algorithm enables the code to be physics agnostic as well as flexibility regarding time-stepping schemes. grim runs on CPUs, as well as on GPUs, using the same code. We formulate a performance model, and use it to show that our implementation runs optimally on both architectures. grim correctly captures classical GRMHD test problems as well as a new suite of linear and nonlinear test problems with anisotropic conduction and viscosity in special and general relativity. As tests and example applications, we resolve the shock substructure due to the presence of dissipation, and report on relativistic versions of the magneto-thermal instability and heat flux driven buoyancy instability, which arise due to anisotropic heat conduction, and of the firehose instability, which occurs due to anisotropic pressure (i.e. viscosity). Finally, we show an example integration of an accretion flow around a Kerr black hole, using Extended GRMHD.

Binary neutron star mergers are promising sources of gravitational waves for ground-based detectors such as Advanced LIGO. Neutron-rich material ejected by these mergers may also be the main source of r-process elements in the Universe, while radioactive decays in the ejecta can power bright electromagnetic post-merger signals. Neutrino-matter interactions play a critical role in the evolution of the composition of the ejected material, which significantly impacts the outcome of nucleosynthesis and the properties of the associated electromagnetic signal. In this work, we present a simulation of a binary neutron star merger using an improved method for estimating the average neutrino energies in our energy-integrated neutrino transport scheme. These energy estimates are obtained by evolving the neutrino number density in addition to the neutrino energy and flux densities. We show that significant changes are observed in the composition of the polar ejecta when comparing our new results with earlier simulations in which the neutrino spectrum was assumed to be the same everywhere in optically thin regions. In particular, we find that material ejected in the polar regions is less neutron rich than previously estimated. Our new estimates of the composition of the polar ejecta make it more likely that the color and timescale of the electromagnetic signal depend on the orientation of the binary with respect to an observer's line-of-sight. These results also indicate that important observable properties of neutron star mergers are sensitive to the neutrino energy spectrum, and may need to be studied through simulations including a more accurate, energy-dependent neutrino transport scheme.

Neutron star-black hole binaries are among the strongest sources of gravitational waves detectable by current observatories. They can also power bright electromagnetic signals (gamma-ray bursts, kilonovae), and may be a significant source of production of r-process nuclei. A misalignment of the black hole spin with respect to the orbital angular momentum leads to precession of that spin and of the orbital plane, and has a significant effect on the properties of the post-merger remnant and of the material ejected by the merger. We present a first set of simulations of precessing neutron star-black hole mergers using a hot, composition dependent, nuclear-theory based equation of state (DD2). We show that the mass of the remnant and of the dynamical ejecta are broadly consistent with the result of simulations using simpler equations of state, while differences arise when considering the dynamics of the merger and the velocity of the ejecta. We show that the latter can easily be understood from assumptions about the composition of low-density, cold material in the different equations of state, and propose an updated estimate for the ejecta velocity which takes those effects into account. We also present an updated mesh-refinement algorithm which allows us to improve the numerical resolution used to evolve neutron star-black hole mergers.

Extracting the unique information on ultradense nuclear matter from the gravitational waves emitted by merging, neutron-star binaries requires robust theoretical models of the signal. We develop a novel effective-one-body waveform model that includes, for the first time, dynamic (instead of only adiabatic) tides of the neutron star, as well as the merger signal for neutron-star--black-hole binaries. We demonstrate the importance of the dynamic tides by comparing our model against new numerical-relativity simulations of nonspinning neutron-star--black-hole binaries spanning more than 24 gravitational-wave cycles, and to other existing numerical simulations for double neutron-star systems. Furthermore, we derive an effective description that makes explicit the dependence of matter effects on two key parameters: tidal deformability and fundamental oscillation frequency.

Neutrino emission significantly affects the evolution of the accretion tori formed in black hole-neutron star mergers. It removes energy from the disk, alters its composition, and provides a potential power source for a gamma-ray burst. To study these effects, simulations in general relativity with a hot microphysical equation of state and neutrino feedback are needed. We present the first such simulation, using a neutrino leakage scheme for cooling to capture the most essential effects and considering a moderate mass (1.4 M_{\odot} neutron star, 5.6 M_{\odot} black hole), high spin (black hole J/M^2=0.9) system with the K_0=220 MeV Lattimer-Swesty equation of state. We find that about 0.08 M_{\odot} of nuclear matter is ejected from the system, while another 0.3 M_{\odot} forms a hot, compact accretion disk. The primary effects of the escaping neutrinos are (i) to make the disk much denser and more compact, (ii) to cause the average electron fraction Y_e of the disk to rise to about 0.2 and then gradually decrease again, and (iii) to gradually cool the disk. The disk is initially hot (T~6 MeV) and luminous in neutrinos (L_{\nu}~10^{54} erg s^{-1}), but the neutrino luminosity decreases by an order of magnitude over 50 ms of post-merger evolution.

Neutron star mergers are among the most promising sources of gravitational waves for advanced ground-based detectors. These mergers are also expected to power bright electromagnetic signals, in the form of short gamma-ray bursts, infrared/optical transients, and radio emission. Simulations of these mergers with fully general relativistic codes are critical to understand the merger and post-merger gravitational wave signals and their neutrinos and electromagnetic counterparts. In this paper, we employ the SpEC code to simulate the merger of low-mass neutron star binaries (two $1.2M_\odot$ neutron stars) for a set of three nuclear-theory based, finite temperature equations of state. We show that the frequency peaks of the post-merger gravitational wave signal are in good agreement with predictions obtained from simulations using a simpler treatment of gravity. We find, however, that only the fundamental mode of the remnant is excited for long periods of time: emission at the secondary peaks is damped on a millisecond timescale in the simulated binaries. For such low-mass systems, the remnant is a massive neutron star which, depending on the equation of state, is either permanently stable or long-lived. We observe strong excitations of l=2, m=2 modes, both in the massive neutron star and in the form of hot, shocked tidal arms in the surrounding accretion torus. We estimate the neutrino emission of the remnant using a neutrino leakage scheme and, in one case, compare these results with a gray two-moment neutrino transport scheme. We confirm the complex geometry of the neutrino emission, also observed in previous simulations with neutrino leakage, and show explicitly the presence of important differences in the neutrino luminosity, disk composition, and outflow properties between the neutrino leakage and transport schemes.

Black holes accreting well below the Eddington rate are believed to have geometrically thick, optically thin, rotationally supported accretion discs in which the Coulomb mean free path is large compared to $GM/c^2$. In such an environment, the disc evolution may differ significantly from ideal magnetohydrodynamic predictions. We present non-ideal global axisymmetric simulations of geometrically thick discs around a rotating black hole. The simulations are carried out using a new code ${\rm\it grim}$, which evolves a covariant extended magnetohydrodynamics model derived by treating non-ideal effects as a perturbation of ideal magnetohydrodynamics. Non-ideal effects are modeled through heat conduction along magnetic field lines, and a difference between the pressure parallel and perpendicular to the field lines. The model relies on an effective collisionality in the disc from wave-particle scattering and velocity-space (mirror and firehose) instabilities. We find that the pressure anisotropy grows to match the magnetic pressure, at which point it saturates due to the mirror instability. The pressure anisotropy produces outward angular momentum transport with a magnitude comparable to that of MHD turbulence in the disc, and a significant increase in the temperature in the wall of the jet. We also find that, at least in our axisymmetric simulations, conduction has a small effect on the disc evolution because (1) the heat flux is constrained to be parallel to the field and the field is close to perpendicular to temperature gradients, and (2) the heat flux is choked by an increase in effective collisionality associated with the mirror instability.

Black holes that accrete far below the Eddington limit are believed to accrete through a geometrically thick, optically thin, rotationally supported plasma that we will refer to as a radiatively inefficient accretion flow (RIAF). RIAFs are typically collisionless in the sense that the Coulomb mean free path is large compared to $GM/c^2$, and relativistically hot near the event horizon. In this paper we develop a phenomenological model for the plasma in RIAFs, motivated by the application to sources such as Sgr A* and M87. The model is derived using Israel-Stewart theory, which considers deviations up to second order from thermal equilibrium, but modified for a magnetized plasma. This leads to thermal conduction along magnetic field lines and a difference in pressure, parallel and perpendicular to the field lines (which is equivalent to anisotrotropic viscosity). In the non-relativistic limit, our model reduces to the widely used Braginskii theory of magnetized, weakly collisional plasmas. We compare our model to the existing literature on dissipative relativistic fluids, describe the linear theory of the plasma, and elucidate the physical meaning of the free parameters in the model. We also describe limits of the model when the conduction is saturated and when the viscosity implies a large pressure anisotropy. In future work, the formalism developed in this paper will be used in numerical models of RIAFs to assess the importance of non-ideal processes for the dynamics and radiative properties of slowly accreting black holes.

We present a first simulation of the post-merger evolution of a black hole-neutron star binary in full general relativity using an energy-integrated general relativistic truncated moment formalism for neutrino transport. We describe our implementation of the moment formalism and important tests of our code, before studying the formation phase of a disk after a black hole-neutron star merger. We use as initial data an existing general relativistic simulation of the merger of a neutron star of 1.4 solar mass with a black hole of 7 solar mass and dimensionless spin a/M=0.8. Comparing with a simpler leakage scheme for the treatment of the neutrinos, we find noticeable differences in the neutron to proton ratio in and around the disk, and in the neutrino luminosity. We find that the electron neutrino luminosity is much lower in the transport simulations, and that the remnant is less neutron-rich. The spatial distribution of the neutrinos is significantly affected by relativistic effects. Over the short timescale evolved, we do not observe purely neutrino-driven outflows. However, a small amount of material (3e-4Msun) is ejected in the polar region during the circularization of the disk. Most of that material is ejected early in the formation of the disk, and is fairly neutron rich. Through r-process nucleosynthesis, that material should produce high-opacity lanthanides in the polar region, and could thus affect the lightcurve of radioactively powered electromagnetic transients. We also show that by the end of the simulation, while the bulk of the disk is neutron-rich, its outer layers have a higher electron fraction. As that material would be the first to be unbound by disk outflows on longer timescales, the changes in Ye experienced during the formation of the disk could have an impact on the nucleosynthesis outputs from neutrino-driven and viscously-driven outflows. [Abridged]

For highly compact neutron stars, constructing numerical initial data for black hole-neutron star binary evolutions is very difficult. We describe improvements to an earlier method that enable it to handle these more challenging cases. We examine the case of a 6:1 mass ratio system in inspiral close to merger, where the star is governed by a polytropic $\Gamma=2$, an SLy, or an LS220 equation of state. In particular, we are able to obtain a solution with a realistic LS220 equation of state for a star with compactness 0.26 and mass 1.98 $M_\odot$, which is representative of the highest reliably determined neutron star masses. For the SLy equation of state, we can obtain solutions with a comparable compactness of 0.25, while for a family of polytropic equations of state, we obtain solutions with compactness up to 0.21, the largest compactness that is stable in this family. These compactness values are significantly higher than any previously published results. We find that improvements in adapting the computational domain to the neutron star surface and in accounting for the center of mass drift of the system are the key ingredients allowing us to obtain these solutions.

Warped accretion discs are expected in many protostellar binary systems. In this paper, we study the long-term evolution of disc warp and precession for discs with dimensionless thickness $H/r$ larger than their viscosity parameter $\alpha$, such that bending waves can propagate and dominate the warp evolution. For small warps, these discs undergo approximately rigid-body precession. We derive analytical expressions for the warp/twist profiles of the disc and the alignment timescale for a variety of models. Applying our results to circumbinary discs, we find that these discs align with the orbital plane of the binary on a timescale comparable to the global precession time of the disc, and typically much smaller than its viscous timescale. We discuss the implications of our finding for the observations of misaligned circumbinary discs (such as KH 15D) and circumbinary planetary systems (such as Kepler-413); these observed misalignments provide useful constraints on the uncertain aspects of the disc warp theory. On the other hand, we find that circumstellar discs can maintain large misalignments with respect to the plane of the binary companion over their entire lifetime. We estimate that inclination angles larger than $\sim 20^\circ$ can be maintained for typical disc parameters. Overall, our results suggest that while highly misaligned circumstellar discs in binaries are expected to be common, such misalignments should be rare for circumbinary discs. These expectations are consistent with current observations of protoplanetary discs and exoplanets in binaries, and can be tested with future observations.

We present a first exploration of the results of neutron star-black hole mergers using black hole masses in the most likely range of $7M_\odot-10M_\odot$, a neutrino leakage scheme, and a modeling of the neutron star material through a finite-temperature nuclear-theory based equation of state. In the range of black hole spins in which the neutron star is tidally disrupted ($\chi_{\rm BH}\gtrsim 0.7$), we show that the merger consistently produces large amounts of cool ($T\lesssim 1\,{\rm MeV}$), unbound, neutron-rich material ($M_{\rm ej}\sim 0.05M_\odot-0.20M_\odot$). A comparable amount of bound matter is initially divided between a hot disk ($T_{\rm max}\sim 15\,{\rm MeV}$) with typical neutrino luminosity $L_\nu\sim 10^{53}\,{\rm erg/s}$, and a cooler tidal tail. After a short period of rapid protonization of the disk lasting $\sim 10\,{\rm ms}$, the accretion disk cools down under the combined effects of the fall-back of cool material from the tail, continued accretion of the hottest material onto the black hole, and neutrino emission. As the temperature decreases, the disk progressively becomes more neutron-rich, with dimmer neutrino emission. This cooling process should stop once the viscous heating in the disk (not included in our simulations) balances the cooling. These mergers of neutron star-black hole binaries with black hole masses $M_{\rm BH}\sim 7M_\odot-10M_\odot$ and black hole spins high enough for the neutron star to disrupt provide promising candidates for the production of short gamma-ray bursts, of bright infrared post-merger signals due to the radioactive decay of unbound material, and of large amounts of r-process nuclei.

Dynamical instabilities in protoneutron stars may produce gravitational waves whose observation could shed light on the physics of core-collapse supernovae. When born with sufficient differential rotation, these stars are susceptible to a shear instability (the "low-T/|W| instability"), but such rotation can also amplify magnetic fields to strengths where they have a considerable impact on the dynamics of the stellar matter. Using a new magnetohydrodynamics module for the Spectral Einstein Code, we have simulated a differentially-rotating neutron star in full 3D to study the effects of magnetic fields on this instability. Though strong toroidal fields were predicted to suppress the low-T/|W| instability, we find that they do so only in a small range of field strengths. Below 4e13 G, poloidal seed fields do not wind up fast enough to have an effect before the instability saturates, while above 5e14 G, magnetic instabilities can actually amplify a global quadrupole mode (this threshold may be even lower in reality, as small-scale magnetic instabilities remain difficult to resolve numerically). Thus, the prospects for observing gravitational waves from such systems are not in fact diminished over most of the magnetic parameter space. Additionally, we report that the detailed development of the low-T/|W| instability, including its growth rate, depends strongly on the particular numerical methods used. The high-order methods we employ suggest that growth might be considerably slower than found in some previous simulations.

We present the first direct comparison of numerical simulations of neutron star-black hole and black hole-black hole mergers in full general relativity. We focus on a configuration with non spinning objects and within the most likely range of mass ratio for neutron star-black hole systems (q=6). In this region of the parameter space, the neutron star is not tidally disrupted prior to merger, and we show that the two types of mergers appear remarkably similar. The effect of the presence of a neutron star on the gravitational wave signal is not only undetectable by the next generation of gravitational wave detectors, but also too small to be measured in the numerical simulations: even the plunge, merger and ringdown signals appear in perfect agreement for both types of binaries. The characteristics of the post-merger remnants are equally similar, with the masses of the final black holes agreeing within dM< 5 10^{-4}M_BH and their spins within da< 10^{-3}M_BH. The rate of periastron advance in the mixed binary agrees with previously published binary black hole results, and we use the inspiral waveforms to place constraints on the accuracy of our numerical simulations independent of algorithmic choices made for each type of binary. Overall, our results indicate that non-disrupting neutron star-black hole mergers are exceptionally well modeled by black hole-black hole mergers, and that given the absence of mass ejection, accretion disk formation, or differences in the gravitational wave signals, only electromagnetic precursors could prove the presence of a neutron star in low-spin systems of total mass ~10Msun, at least until the advent of gravitational wave detectors with a sensitivity comparable to that of the proposed Einstein Telescope.

Black hole-neutron star (BHNS) binaries are important sources of gravitational waves for second-generation interferometers, and BHNS mergers are also a proposed engine for short, hard gamma-ray bursts. The behavior of both the spacetime (and thus the emitted gravitational waves) and the neutron star matter in a BHNS merger depend strongly and nonlinearly on the black hole's spin. While there is a significant possibility that astrophysical black holes could have spins that are nearly extremal (i.e. near the theoretical maximum), to date fully relativistic simulations of BHNS binaries have included black-hole spins only up to $S/M^2$=0.9, which corresponds to the black hole having approximately half as much rotational energy as possible, given the black hole's mass. In this paper, we present a new simulation of a BHNS binary with a mass ratio $q=3$ and black-hole spin $S/M^2$=0.97, the highest simulated to date. We find that the black hole's large spin leads to the most massive accretion disk and the largest tidal tail outflow of any fully relativistic BHNS simulations to date, even exceeding the results implied by extrapolating results from simulations with lower black-hole spin. The disk appears to be remarkably stable. We also find that the high black-hole spin persists until shortly before the time of merger; afterwards, both merger and accretion spin down the black hole.

The Kepler satellite has discovered a number of transiting planets around close binary stars. These circumbinary systems have highly aligned planetary and binary orbits. In this paper, we explore how the mutual inclination between the planetary and binary orbits may reflect the physical conditions of the assembly of protoplanetary disks and the interaction between protostellar binaries and circumbinary disks. Given the turbulent nature of star-forming molecular clouds, it is possible that the gas falling onto the outer region of a circumbinary disk and the central protostellar binary have different axes of rotation. Thus, the newly assembled circumbinary disk can be misaligned with respect to the binary. However, the gravitational torque from the binary produces a warp and twist in the disk, and the back-reaction torque tends to align the disk and the binary orbital plane. We present a new, analytic calculation of this alignment torque, and show that the binary-disk inclination angle can be reduced appreciably after the binary accretes a few percent of its mass from the disk. Our calculation suggests that in the absence of other disturbances, circumbinary disks and planets around close (sub-AU) stellar binaries, for which mass accretion onto the proto-binary is very likely to have occurred, are expected to be highly aligned with the binary orbits, while disks and planets around wide binaries can be misaligned. Measurements of the mutual inclinations of circumbinary planetary systems can provide a clue to the birth environments of such systems.

Black hole-neutron star mergers resulting in the disruption of the neutron star and the formation of an accretion disk and/or the ejection of unbound material are prime candidates for the joint detection of gravitational-wave and electromagnetic signals when the next generation of gravitational-wave detectors comes online. However, the disruption of the neutron star and the properties of the post-merger remnant are very sensitive to the parameters of the binary. In this paper, we study the impact of the radius of the neutron star and the alignment of the black hole spin for systems within the range of mass ratio currently deemed most likely for field binaries (M_BH ~ 7 M_NS) and for black hole spins large enough for the neutron star to disrupt (J/M^2=0.9). We find that: (i) In this regime, the merger is particularly sensitive to the radius of the neutron star, with remnant masses varying from 0.3M_NS to 0.1M_NS for changes of only 2 km in the NS radius; (ii) 0.01-0.05M_sun of unbound material can be ejected with kinetic energy >10^51 ergs, a significant increase compared to low mass ratio, low spin binaries. This ejecta could power detectable optical and radio afterglows. (iii) Only a small fraction (<3%) of the Advanced LIGO events in this parameter range have gravitational-wave signals which could offer constraints on the equation of state of the neutron star. (iv) A misaligned black hole spin works against disk formation, with less neutron star material remaining outside of the black hole after merger, and a larger fraction of that material remaining in the tidal tail instead of the forming accretion disk. (v) Large kicks (v>300 km/s) can be given to the final black hole as a result of a precessing BHNS merger, when the disruption of the neutron star occurs just outside or within the innermost stable spherical orbit.

Determining the final result of black hole-neutron star mergers, and in particular the amount of matter remaining outside the black hole at late times and its properties, has been one of the main motivations behind the numerical simulation of these systems. Black hole-neutron star binaries are amongst the most likely progenitors of short gamma-ray bursts --- as long as massive (probably a few percents of a solar mass), hot accretion disks are formed around the black hole. Whether this actually happens strongly depends on the physical characteristics of the system, and in particular on the mass ratio, the spin of the black hole, and the radius of the neutron star. We present here a simple two-parameter model, fitted to existing numerical results, for the determination of the mass remaining outside the black hole a few milliseconds after a black hole-neutron star merger (i.e. the combined mass of the accretion disk, the tidal tail, and the potential ejecta). This model predicts the remnant mass within a few percents of the mass of the neutron star, at least for remnant masses up to 20% of the neutron star mass. Results across the range of parameters deemed to be the most likely astrophysically are presented here. We find that, for 10 solar mass black holes, massive disks are only possible for large neutron stars (R>12km), or quasi-extremal black hole spins (a/M>0.9). We also use our model to discuss how the equation of state of the neutron star affects the final remnant, and the strong influence that this can have on the rate of short gamma-ray bursts produced by black hole-neutron star mergers.

General relativistic simulations of black hole-neutron star mergers have currently been limited to low-mass black holes (less than 7 solar mass), even though population synthesis models indicate that a majority of mergers might involve more massive black holes (10 solar mass or more). We present the first general relativistic simulations of black hole-neutron star mergers with 10 solar mass black holes. For massive black holes, the tidal forces acting on the neutron star are usually too weak to disrupt the star before it reaches the innermost stable circular orbit of the black hole. Varying the spin of the black hole in the range a/M = 0.5-0.9, we find that mergers result in the disruption of the star and the formation of a massive accretion disk only for large spins a/M>0.7-0.9. From these results, we obtain updated constraints on the ability of BHNS mergers to be the progenitors of short gamma-ray bursts as a function of the mass and spin of the black hole. We also discuss the dependence of the gravitational wave signal on the black hole parameters, and provide waveforms and spectra from simulations beginning 7-8 orbits before merger.

The spin of black holes in black hole-neutron star (BHNS) binaries can have a strong influence on the merger dynamics and the postmerger state; a wide variety of spin magnitudes and orientations are expected to occur in nature. In this paper, we report the first simulations in full general relativity of BHNS mergers with misaligned black hole spin. We vary the spin magnitude from a/m=0 to a/m=0.9 for aligned cases, and we vary the misalignment angle from 0 to 80 degrees for a/m=0.5. We restrict our study to 3:1 mass ratio systems and use a simple Gamma-law equation of state. We find that the misalignment angle has a strong effect on the mass of the postmerger accretion disk, but only for angles greater than ~ 40 degrees. Although the disk mass varies significantly with spin magnitude and misalignment angle, we find that all disks have very similar lifetimes ~ 100ms. Their thermal and rotational profiles are also very similar. For a misaligned merger, the disk is tilted with respect to the final black hole's spin axis. This will cause the disk to precess, but on a timescale longer than the accretion time. In all cases, we find promising setups for gamma-ray burst production: the disks are hot, thick, and hyperaccreting, and a baryon-clear region exists above the black hole.

Magnetic interactions between a protostar and its accretion disc tend to induce warping in the disc and produce secular changes in the stellar spin direction, so that the spin axis may not always be perpendicular to the disc. This may help explain the recently observed spin-orbit misalignment in a number of exoplanetary systems. We study the dynamics of warped protoplanetary discs under the combined effects of magnetic warping/precession torques and internal stresses in the disc, including viscous damping of warps and propagation of bending waves. We show that when the outer disc axis is misaligned with the stellar spin axis, the disc evolves towards a warped steady-state on a timescale that depends on the disc viscosity or the bending wave propagation speed, but in all cases is much shorter than the timescale for the spin evolution (of order of a million years). Moreover, for the most likely physical parameters characterizing magnetic protostars, circumstellar discs and their interactions, the steady-state disc has a rather small warp, such that the whole disc lies approximately in a single plane determined by the outer disc boundary conditions, although more extreme parameters may give rise to larger disc warps. In agreement with our recent analysis (Lai et al. 2010) based on flat discs, we find that the back-reaction magnetic torques of the slightly warped disc on the star can either align the stellar spin axis with the disc axis or push it towards misalignment, depending on the parameters of the star-disc system. This implies that newly formed planetary systems may have a range of inclination angles between the stellar spin axis and the symmetry axis of the planetary orbits.

Recent observations have shown that in many exoplanetary systems the spin axis of the parent star is misaligned with the planet's orbital axis. These have been used to argue against the scenario that short-period planets migrated to their present-day locations due to tidal interactions with their natal discs. However, this interpretation is based on the assumption that the spins of young stars are parallel to the rotation axes of protostellar discs around them. We show that the interaction between a magnetic star and its circumstellar disc can (but not always) have the effect of pushing the stellar spin axis away from the disc angular momentum axis toward the perpendicular state and even the retrograde state. Planets formed in the disc may therefore have their orbital axes misaligned with the stellar spin axis, even before any additional planet-planet scatterings or Kozai interactions take place. In general, magnetosphere--disc interactions lead to a broad distribution of the spin--orbit angles, with some systems aligned and other systems misaligned.

The merger dynamics of a black hole-neutron star (BHNS) binary is influenced by the neutron star equation of state (EoS) through the latter's effect on the neutron star's radius and on the character of the mass transfer onto the black hole. We study these effects by simulating a number of BHNS binaries in full general relativity using a mixed pseudospectral/finite difference code. We consider several models of the neutron star matter EoS, including Gamma=2 and Gamma=2.75 polytropes and the nuclear-theory based Shen EoS. For models using the Shen EoS, we consider two limits for the evolution of the composition: source-free advection and instantaneous beta-equilibrium. To focus on EoS effects, we fix the mass ratio to 3:1 and the initial aligned black hole spin to a/m=0.5 for all models. We confirm earlier studies which found that more compact stars create a stronger gravitational wave signal but a smaller postmerger accretion disk. We also vary the EoS while holding the compaction fixed. All mergers are qualitatively similar, but we find signatures of the EoS in the waveform and in the tail and disk structures.

We present a code for solving the coupled Einstein-hydrodynamics equations to evolve relativistic, self-gravitating fluids. The Einstein field equations are solved in generalized harmonic coordinates on one grid using pseudospectral methods, while the fluids are evolved on another grid using shock-capturing finite difference or finite volume techniques. We show that the code accurately evolves equilibrium stars and accretion flows. Then we simulate an equal-mass nonspinning black hole-neutron star binary, evolving through the final four orbits of inspiral, through the merger, to the final stationary black hole. The gravitational waveform can be reliably extracted from the simulation.

We present a new numerical scheme to solve the initial value problem for black hole-neutron star binaries. This method takes advantage of the flexibility and fast convergence of a multidomain spectral representation of the initial data to construct high-accuracy solutions at a relatively low computational cost. We provide convergence tests of the method for both isolated neutron stars and irrotational binaries. In the second case, we show that we can resolve the small inconsistencies that are part of the quasi-equilibrium formulation, and that these inconsistencies are significantly smaller than observed in previous works. The possibility of generating a wide variety of initial data is also demonstrated through two new configurations inspired by results from binary black holes. First, we show that choosing a modified Kerr-Schild conformal metric instead of a flat conformal metric allows for the construction of quasi-equilibrium binaries with a spinning black hole. Second, we construct binaries in low-eccentricity orbits, which are a better approximation to astrophysical binaries than quasi-equilibrium systems.