We search for an isotropic stochastic gravitational-wave background (GWB) in the newly released 11-year dataset from the North American Nanohertz Observatory for Gravitational Waves (NANOGrav). While we find no evidence for a GWB, we place constraints on a population of inspiraling supermassive black hole (SMBH) binaries, a network of decaying cosmic strings, and a primordial GWB. For the first time, we find that the GWB constraints are sensitive to the Solar System ephemeris (SSE) model used, and that SSE errors arXiv:1801.02617v2 [astro-ph.HE] 7 Jun 2018 2 THE NANOGRAV COLLABORATION can mimic a GWB signal. We developed an approach that bridges systematic SSE differences, producing the first PTA constraints that are robust against SSE errors. We thus place a 95% upper limit on the GW strain amplitude of A GWB < 1.45 × 10 −15 at a frequency of f = 1-yr −1 for a fiducial f −2/3 power-law spectrum, and with inter-pulsar correlations modeled. This is a factor of ∼ 2 improvement over the NANOGrav 9-year limit, calculated using the same procedure. Previous PTA upper limits on the GWB (as well as their astrophysical and cosmological interpretations) will need revision in light of SSE systematic errors. We use our constraints to characterize the combined influence on the GWB of the stellar mass-density in galactic cores, the eccentricity of SMBH binaries, and SMBH-galactic-bulge scaling relationships. We constrain cosmic-string tension using recent simulations, yielding an SSE-marginalized 95% upper limit of Gµ < 5.3 × 10 −11 -a factor of ∼ 2 better than the published NANOGrav 9-year constraints. Our SSE-marginalized 95% upper limit on the energy density of a primordial GWB (for a radiation-dominated post-inflation Universe) is Ω GWB ( f )h 2 < 3.4 × 10 −10 .

We describe the Einstein Toolkit, a community-driven, freely accessible computational infrastructure intended for use in numerical relativity, relativistic astrophysics, and other applications. The Toolkit, developed by a collaboration involving researchers from multiple institutions around the world, combines a core set of components needed to simulate astrophysical objects such as black holes, compact objects, and collapsing stars, as well as a full suite of analysis tools. The Einstein Toolkit is currently based on the Cactus Framework for high-performance computing and the Carpet adaptive mesh refinement driver. It implements spacetime evolution via the BSSN evolution system and general-relativistic hydrodynamics in a finite-volume discretization. The toolkit is under continuous development and contains many new code components that have been publicly released for the first time and are described in this article. We discuss the motivation behind the release of the toolkit, the philosophy underlying its development, and the goals of the project. A summary of the implemented numerical techniques is included, as are results of numerical test covering a variety of sample astrophysical problems.

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 an accretion 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 mass 1.4M with a black hole of mass 7M and dimensionless spin χBH = 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 both the disk and the disk outflows are less neutron-rich. The spatial distribution of the neutrinos is significantly affected by relativistic effects, due to large velocities and curvature in the regions of strongest emission. Over the short timescale evolved, we do not observe purely neutrino-driven outflows. However, a small amount of material (3 × 10 −4 M ) 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 (electron fraction Ye ∼ 0.15 − 0.25). 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 remains neutron-rich (Ye ∼ 0.15 − 0.2 and decreasing), its outer layers have a higher electron fraction: 10% of the remaining mass has Ye > 0.3. As that material would be the first to be unbound by disk outflows on longer timescales, and as composition evolution is slower at later times, the changes in Ye experienced during the formation phase of the disk could have an impact on nucleosynthesis outputs from neutrino-driven and viscously-driven outflows. Finally, we find that the effective viscosity due to momentum transport by neutrinos is unlikely to have a strong effect on the growth of the magnetorotational instability in the post-merger accretion disk.

We present results of new three-dimensional (3D) general-relativistic magnetohydrodynamic simulations of rapidly rotating strongly magnetized core collapse. These simulations are the first of their kind and include a microphysical finite-temperature equation of state and a leakage scheme that captures the overall energetics and lepton number exchange due to postbounce neutrino emission. Our results show that the 3D dynamics of magnetorotational core-collapse supernovae are fundamentally different from what was anticipated on the basis of previous simulations in axisymmetry (2D). A strong bipolar jet that develops in a simulation constrained to 2D is crippled by a spiral instability and fizzles in full 3D. While multiple (magneto-)hydrodynamic instabilities may be present, our analysis suggests that the jet is disrupted by an m = 1 kink instability of the ultra-strong toroidal field near the rotation axis. Instead of an axially symmetric jet, a completely new, previously unreported flow structure develops. Highly magnetized spiral plasma funnels expelled from the core push out the shock in polar regions, creating wide secularly expanding lobes. We observe no runaway explosion by the end of the full 3D simulation at 185 ms after bounce. At this time, the lobes have reached maximum radii of ∼900 km.

Neutron star mergers are among the most promising sources of gravitational waves for advanced groundbased detectors. These mergers are also expected to power bright electromagnetic signals, in the form of short gamma-ray bursts, infrared/optical transients powered by r-process nucleosynthesis in neutron-rich material ejected by the merger, and radio emission from the interaction of that ejecta with the interstellar medium. 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 Spectral Einstein Code (SpEC) to simulate the merger of low-mass neutron star binaries (two 1.2M 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 recent 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 (i.e. rapid uniform rotation is sufficient to prevent its collapse). 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.

Magnetohydrodynamic turbulence is important in many high-energy astrophysical systems, where instabilities can amplify the local magnetic field over very short timescales. Specifically, the magnetorotational instability and dynamo action have been suggested as a mechanism for the growth of magnetar-strength magnetic fields (of 10(15) gauss and above) and for powering the explosion of a rotating massive star. Such stars are candidate progenitors of type Ic-bl hypernovae, which make up all supernovae that are connected to long γ-ray bursts. The magnetorotational instability has been studied with local high-resolution shearing-box simulations in three dimensions, and with global two-dimensional simulations, but it is not known whether turbulence driven by this instability can result in the creation of a large-scale, ordered and dynamically relevant field. Here we report results from global, three-dimensional, general-relativistic magnetohydrodynamic turbulence simulations. We show that hydromagnetic turbulence in rapidly rotating protoneutron stars produces an inverse cascade of energy. We find a large-scale, ordered toroidal field that is consistent with the formation of bipolar magnetorotationally driven outflows. Our results demonstrate that rapidly rotating massive stars are plausible progenitors for both type Ic-bl supernovae and long γ-ray bursts, and provide a viable mechanism for the formation of magnetars. Moreover, our findings suggest that rapidly rotating massive stars might lie behind potentially magnetar-powered superluminous supernovae.

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 ⊙ -10M ⊙ , 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 (χ BH ≳ 0.7), we show that the merger consistently produces large amounts of cool (T ≲ 1 MeV), unbound, neutron-rich material (M ej ∼ 0.05M ⊙ -0.20M ⊙ ). A comparable amount of bound matter is initially divided between a hot disk (T max ∼ 15 MeV) with typical neutrino luminosity of L ν ∼ 10 53 erg=s, and a cooler tidal tail. After a short period of rapid protonization of the disk lasting ∼10 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 of M BH ∼ 7M ⊙ -10M ⊙ , 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 postmerger signals due to the radioactive decay of unbound material, and of large amounts of r-process nuclei.

We study the three-dimensional (3D) hydrodynamics of the post-core-bounce phase of the collapse of a 27 M star and pay special attention to the development of the standing accretion shock instability (SASI) and neutrino-driven convection. To this end, we perform 3D general-relativistic simulations with a three-species neutrino leakage scheme. The leakage scheme captures the essential aspects of neutrino cooling, heating, and lepton number exchange as predicted by radiation-hydrodynamics simulations. The 27 M progenitor was studied in 2D by Müller et al., who observed strong growth of the SASI while neutrino-driven convection was suppressed. In our 3D simulations, neutrino-driven convection grows from numerical perturbations imposed by our Cartesian grid. It becomes the dominant instability and leads to large-scale non-oscillatory deformations of the shock front. These will result in strongly aspherical explosions without the need for large-scale SASI shock oscillations. Low--mode SASI oscillations are present in our models, but saturate at small amplitudes that decrease with increasing neutrino heating and vigor of convection. Our results, in agreement with simpler 3D Newtonian simulations, suggest that once neutrino-driven convection is started, it is likely to become the dominant instability in 3D. Whether it is the primary instability after bounce will ultimately depend on the physical seed perturbations present in the cores of massive stars. The gravitational wave signal, which we extract and analyze for the first time from 3D general-relativistic models, will serve as an observational probe of the postbounce dynamics and, in combination with neutrinos, may allow us to determine the primary hydrodynamic instability.