The 2017 detection of the inspiral and merger of two neutron stars in gravitational waves and gamma rays was accompanied by a quickly-reddening transient. Such a transient was predicted to occur following a rapid neutron capture (r-process) nucleosynthesis event, which synthesizes neutronrich, radioactive nuclei and can take place in both dynamical ejecta and in the wind driven off the accretion torus formed after a neutron star merger. We present the first three-dimensional general relativistic, full transport neutrino radiation magnetohydrodynamics (GRRMHD) simulations of the black hole-accretion disk-wind system produced by the GW170817 merger. We show that the small but non-negligible optical depths lead to neutrino transport globally coupling the disk electron fraction, which we capture by solving the transport equation with a Monte Carlo method. The resulting absorption drives up the electron fraction in a structured, continuous outflow, with electron fraction as high as Ye∼0.4 in the extreme polar region. We show via nuclear reaction network and radiative transfer calculations that nucleosynthesis in the disk wind will produce a blue kilonova.

GRS 1915+105 harbors one of the most massive known stellar black holes in the Galaxy. In 2007 May, we observed GRS 1915+105 for ∼117 ks in the low/hard state using Suzaku. We collected and analyzed the data with the Hard X-ray Detector/Positive Intrinsic Negative and X-ray Spectrometer cameras spanning the energy range from 2.3 to 55 keV. Fits to the spectra with simple models reveal strong disk reflection through an Fe K emission line and a Compton backscattering hump. We report constraints on the spin parameter of the black hole in GRS 1915 + 105 using relativistic disk reflection models. The model for the soft X-ray spectrum (i.e., < 10 keV) suggestŝ a = 0.56 +0.02 −0.02 and excludes zero spin at the 4σ level of confidence. The model for the full broadband spectrum suggests that the spin may be higher,â = 0.98 +0.01 −0.01 (1σ confidence), and again excludes zero spin at the 2σ level of confidence. We discuss these results in the context of other spin constraints and inner disk studies in GRS 1915 + 105.

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 taskbased 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. variety of astrophysics codes (e.g., Refs. [6,[9][10][11][12][13][14]) have been designed based on these fundamental building blocks.These strategies work well when the computations are reasonably homogeneous or when one seeks good parallelization to only a few thousand cores. As the number of MPI processes increases, so does the cost of communication which, together with non-uniform workload typical of astrophysics problems, limits the maximum number of useful cores that codes can run on. Efficient core utilization becomes non-trivial, often requiring careful optimization by hand to achieve good scalability [15]. Standard finite-volume and finitedifference methods achieve higher order accuracy with increasingly large (overlapping) stencil sizes, and may require additional effort to achieve scalability on massively parallel machines.As one looks ahead to the arrival of exascale computing, it will become increasingly important to focus on developing algorithms that can take full advantage of these very large machines.Discontinuous Galerkin (DG) methods [16-21], together with a task-based parallelization strategy, have the potential to tackle many of these problems. DG methods offer high-order accuracy in smooth regions (although, for stability, increasing the scheme's order requires decreasing the timestep, which restricts the largest usable order in practice), robustness for shocks and other discontinuities, and grid flexibility including a formulation that allows for comparatively straightforward hp-adaptivity and local timestepping. DG methods can be combined with positivity preserving strategies [22][23][24] or "atmosphere treatments" [25] which seek to maintain non-negative values of the pressure and density in challenging regions such as those containing high-speed astrophysical flow. DG methods are also well suited for parallelization: Their formulation in terms of l...

The detailed observations of GW170817 proved for the first time directly that neutron star mergers are a major production site of heavy elements. The observations could be fit by a number of simulations that qualitatively agree, but can quantitatively differ (e.g., in total r-process mass) by an order of magnitude. We categorize kilonova ejecta into several typical morphologies motivated by numerical simulations, and apply a radiative transfer Monte Carlo code to study how the geometric distribution of the ejecta shapes the emitted radiation. We find major impacts on both spectra and light curves. The peak bolometric luminosity can vary by two orders of magnitude and the timing of its peak by a factor of five. These findings provide the crucial implication that the ejecta masses inferred from observations around the peak brightness are uncertain by at least an order of magnitude. Mixed two-component models with lanthanide-rich ejecta are particularly sensitive to geometric distribution. A subset of mixed models shows very strong viewing angle dependence due to lanthanide “curtaining,” which persists even if the relative mass of lanthanide-rich component is small. The angular dependence is weak in the rest of our models, but different geometric combinations of the two components lead to a highly diverse set of light curves. We identify geometry-dependent P Cygni features in late spectra that directly map out strong lines in the simulated opacity of neodymium, which can help to constrain the ejecta geometry and to directly probe the r-process abundances.

We model a compact black hole-accretion disk system in the collapsar scenario with full transport, frequency dependent, general relativistic radiation magnetohydrodynamics. We examine whether or not winds from a collapsar disk can undergo rapid neutron capture (r-process) nucleosynthesis and significantly contribute to solar r-process abundances. We find the inclusion of accurate transport has significant effects on outflows, raising the electron fraction above and preventing third-peak r-process material from being synthesized. We analyze the time evolution of neutrino processes and electron fraction in the disk and present a simple one-dimensional model for the vertical structure that emerges. We compare our simulation to semi-analytic expectations and argue that accurate neutrino transport and realistic initial and boundary conditions are required to capture the dynamics and nucleosynthetic outcome of a collapsar.

We study a lattice regularization of the gravitational path integral-causal dynamical triangulationsfor (2 + 1)-dimensional Einstein gravity with positive cosmological constant in the presence of past and future spacelike boundaries of fixed intrinsic geometries. For spatial topology of a 2-sphere, we determine the form of the Einstein-Hilbert action supplemented by the Gibbons-Hawking-York boundary terms within the Regge calculus of causal triangulations. Employing this action we numerically simulate a variety of transition amplitudes from the past boundary to the future boundary. To the extent that we have so far investigated them, these transition amplitudes appear consistent with the gravitational effective action previously found to characterize the ground state of quantum spacetime geometry within the Euclidean de Sitter-like phase. Certain of these transition amplitudes convincingly demonstrate that the so-called stalks present in this phase are numerical artifacts of the lattice regularization, seemingly indicate that the quantization technique of causal dynamical triangulations differs in detail from that of the no-boundary proposal of Hartle and Hawking, and possibly represent the first numerical simulations of portions of temporally unbounded quantum spacetime geometry within the causal dynamical triangulations approach. We also uncover tantalizing evidence suggesting that Lorentzian not Euclidean de Sitter spacetime dominates the ground state on sufficiently large scales.

The merger of neutron star binaries is believed to eject a wide range of heavy elements into the universe. By observing the emission from this ejecta, scientists can probe the ejecta properties (mass, velocity, and composition distributions). The emission (a.k.a. kilonova) is powered by the radioactive decay of the heavy isotopes produced in the merger and this emission is reprocessed by atomic opacities to optical and infrared wavelengths. Understanding the ejecta properties requires calculating the dependence of this emission on these opacities. The strong lines in the optical and infrared in lanthanide opacities have been shown to significantly alter the light curves and spectra in these wavelength bands, arguing that the emission in these wavelengths can probe the composition of this ejecta. Here we study variations in the kilonova emission by varying individual lanthanide (and the actinide uranium) concentrations in the ejecta. The broad forest of lanthanide lines makes it difficult to determine the exact fraction of individual lanthanides. Nd is an exception. Its opacities above 1 μm are higher than other lanthanides and observations of kilonovae can potentially probe increased abundances of Nd. Similarly, at early times when the ejecta is still hot (first day), the U opacity is strong in the 0.2–1 μm wavelength range and kilonova observations may also be able to constrain these abundances.

Multi-messenger astrophysics is a fast-growing, interdisciplinary field that combines data, which vary in volume and speed of data processing, from many different instruments that probe the Universe using different cosmic messengers: electromagnetic waves, cosmic rays, gravitational waves and neutrinos. In this Expert Recommendation, we review the key challenges of real-time observations of gravitational wave sources and their electromagnetic and astroparticle counterparts, and make a number of recommendations to maximize their potential for scientific discovery. These recommendations refer to the design of scalable and computationally efficient machine learning algorithms; the cyber-infrastructure to numerically simulate astrophysical sources, and to process and interpret multi-messenger astrophysics data; the management of gravitational wave detections to trigger real-time alerts for electromagnetic and astroparticle follow-ups; a vision to harness future developments of machine learning and cyber-infrastructure resources to cope with the big-data requirements; and the need to build a community of experts to realize the goals of multi-messenger astrophysics.