Neutrinos play a critical role of transporting energy and changing the lepton density within core-collapse supernovae and neutron star mergers. The quantum kinetic equations (QKEs) combine the effects of neutrino-matter interactions treated in classical Boltzmann transport with the neutrino flavor-changing effects treated in neutrino oscillation calculations. We present a method for extending existing neutrino interaction rates to full QKE source terms for use in numerical calculations. We demonstrate the effects of absorption and emission by nucleons and nuclei, electron scattering, electron-positron pair annihilation, nucleon-nucleon bremsstrahlung, neutrino-neutrino scattering. For the first time, we include all these collision terms self-consistently in a simulation of the full isotropic QKEs in conditions relevant to core-collapse supernovae and neutron star mergers. For our choice of parameters, the long-term evolution of the neutrino distribution function proceeds similarly with and without the oscillation term, though with measurable differences. We demonstrate that electron scattering, nucleon-nucleon bremsstrahlung processes, and four-neutrino processes dominate flavor decoherence in the protoneutron star (PNS), absorption dominates near the shock, and all of the considered processes except elastic nucleon scattering are relevant in the decoupling region. Finally, we propose an effective decoherence opacity that at most energies predicts decoherence rates to within a factor of 10 in our model PNS and within 20% outside of the PNS.

We perform an extensive study of the influence of nuclear weak interactions on core-collapse supernovae (CCSNe), paying particular attention to consistency between nuclear abundances in the equation of state (EOS) and nuclear weak interactions. We compute properties of uniform matter based on the variational method. For inhomogeneous nuclear matter, we take a full ensemble of nuclei into account with various finite-density and thermal effects and directly use the nuclear abundances to compute nuclear weak interaction rates. To quantify the impact of a consistent treatment of nuclear abundances on CCSN dynamics, we carry out spherically symmetric CCSN simulations with full Boltzmann neutrino transport, systematically changing the treatment of weak interactions, EOSs, and progenitor models. We find that the inconsistent treatment of nuclear abundances between the EOS and weak interaction rates weakens the EOS dependence of both the dynamics and neutrino signals. We also test the validity of two artificial prescriptions for weak interactions of light nuclei and find that both prescriptions affect the dynamics. Furthermore, there are differences in neutrino luminosities by ~10% and in average neutrino energies by 0.25-1 MeV from those of the fiducial model. We also find that the neutronization burst neutrino signal depends on the progenitor more strongly than on the EOS, preventing a detection of this signal from constraining the EOS.

In July 2018 an FRIB Theory Alliance program was held on the implications of GW170817 and its associated kilonova for r-process nucleosynthesis. Topics of discussion included the astrophysical and nuclear physics uncertainties in the interpretation of the GW170817 kilonova, what we can learn about the astrophysical site or sites of the r process from this event, and the advances in nuclear experiment and theory most crucial to pursue in light of the new data. Here we compile a selection of scientific contributions to the workshop, broadly representative of progress in r-process studies since the GW170817 event.

The mechanism driving core-collapse supernovae is sensitive to the interplay between matter and neutrino radiation. However, neutrino radiation transport is very difficult to simulate, and several radiation transport methods of varying levels of approximation are available. We carefully compare for the first time in multiple spatial dimensions the discrete ordinates (DO) code of Nagakura, Yamada, and Sumiyoshi and the Monte Carlo (MC) code Sedonu, under the assumptions of a static fluid background, flat spacetime, elastic scattering, and full special relativity. We find remarkably good agreement in all spectral, angular, and fluid interaction quantities, lending confidence to both methods. The DO method excels in determining the heating and cooling rates in the optically thick region. The MC method predicts sharper angular features due to the effectively infinite angular resolution, but struggles to drive down noise in quantities where subtractive cancellation is prevalent, such as the net gain in the protoneutron star and off-diagonal components of the Eddington tensor. We also find that errors in the angular moments of the distribution functions induced by neglecting velocity dependence are sub-dominant to those from limited momentum-space resolution. We briefly compare directly computed second angular moments to those predicted by popular algebraic two-moment closures, and find that the errors from the approximate closures are comparable to the difference between the DO and MC methods. Included in this work is an improved Sedonu code, which now implements a fully special relativistic, time-independent version of the grid-agnostic Monte Carlo random walk approximation.

Gravitational waves (GWs) generated by axisymmetric rotating collapse, bounce, and early postbounce phases of a galactic core-collapse supernova will be detectable by current-generation gravitational wave observatories. Since these GWs are emitted from the quadrupole-deformed nuclear-density core, they may encode information on the uncertain nuclear equation of state (EOS). We examine the effects of the nuclear EOS on GWs from rotating core collapse and carry out 1824 axisymmetric general-relativistic hydrodynamic simulations that cover a parameter space of 98 different rotation profiles and 18 different EOS. We show that the bounce GW signal is largely independent of the EOS and sensitive primarily to the ratio of rotational to gravitational energy, and at high rotation rates, to the degree of differential rotation. The GW frequency of postbounce core oscillations shows stronger EOS dependence that can be parameterized by the core's EOS-dependent dynamical frequency $\sqrt{G\bar{\rho}_c}$. We find that the ratio of the peak frequency to the dynamical frequency follows a universal trend that is obeyed by all EOS and rotation profiles and that indicates that the nature of the core oscillations changes when the rotation rate exceeds the dynamical frequency. We find that differences in the treatments of low-density nonuniform nuclear matter, of the transition from nonuniform to uniform nuclear matter, and in the description of nuclear matter up to around twice saturation density can mildly affect the GW signal. We find that approximations and uncertainties in electron capture rates can lead to variations in the GW signal that are of comparable magnitude to those due to different nuclear EOS. This emphasizes the need for reliable nuclear electron capture rates and for self-consistent multi-dimensional neutrino radiation-hydrodynamic simulations of rotating core collapse.

We present Sedonu, a new open source, steady-state, special relativistic Monte Carlo (MC) neutrino transport code, available at bitbucket.org/srichers/sedonu. The code calculates the energy- and angle-dependent neutrino distribution function on fluid backgrounds of any number of spatial dimensions, calculates the rates of change of fluid internal energy and electron fraction, and solves for the equilibrium fluid temperature and electron fraction. We apply this method to snapshots from two-dimensional simulations of accretion disks left behind by binary neutron star mergers, varying the input physics and comparing to the results obtained with a leakage scheme for the case of a central black hole and a central hypermassive neutron star. Neutrinos are guided away from the densest regions of the disk and escape preferentially around 45 degrees from the equatorial plane. Neutrino heating is strengthened by MC transport a few scale heights above the disk midplane near the innermost stable circular orbit, potentially leading to a stronger neutrino-driven wind. Neutrino cooling in the dense midplane of the disk is stronger when using MC transport, leading to a globally higher cooling rate by a factor of a few and a larger leptonization rate by an order of magnitude. We calculate neutrino pair annihilation rates and estimate that an energy of 2.8e46 erg is deposited within 45 degrees of the symmetry axis over 300 ms when a central BH is present. Similarly, 1.9e48 erg is deposited over 3 s when an HMNS sits at the center, but neither estimate is likely to be sufficient to drive a GRB jet.

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.

Global disk simulations provide a powerful tool for investigating accretion and the underlying magnetohydrodynamic turbulence driven by the magneto-rotational instability (MRI). Using them to predict accurately quantities such as stress, accretion rate, and surface brightness profile requires that purely numerical effects, arising from both resolution and algorithm, be understood and controlled. We use the flux-conservative Athena code to conduct a series of experiments on disks having a variety of magnetic topologies to determine what constitutes adequate resolution. We develop and apply several resolution metrics: Qz and Qphi, the ratio of the grid zone size to the characteristic MRI wavelength, alpha_mag, the ratio of the Maxwell stress to the magnetic pressure, and the ratio of radial to toroidal magnetic field energy. For the initial conditions considered here, adequate resolution is characterized by Qz > 15, Qphi > 20, alpha_mag = 0.45, and a field energy ratio of 0.2. These values are associated with > 35 zones per scaleheight, a result consistent with shearing box simulations. Numerical algorithm is also important. Use of the HLLE flux solver or second-order interpolation can significantly degrade the effective resolution compared to the HLLD flux solver and third-order interpolation. Resolution at this standard can be achieved only with large numbers of grid zones, arranged in a fashion that matches the symmetries of the problem and the scientific goals of the simulation.

Centaurus A is the closest active galactic nucleus. High resolution imaging using Very Long Baseline Interferometry (VLBI) enables us to study the spectral and kinematic behavior of the radio jet-counterjet system on sub-parsec scales, providing essential information for jet emission and formation models. Our aim is to study the structure and spectral shape of the emission from the central-parsec region of Cen A. As a target of the Southern Hemisphere VLBI monitoring program TANAMI (Tracking Active Galactic Nuclei with Milliarcsecond Interferometry), VLBI observations of Cen A are made regularly at 8.4 and 22.3 GHz with the Australian Long Baseline Array (LBA) and associated telescopes in Antarctica, Chile, and South Africa. The first dual-frequency images of this source are presented along with the resulting spectral index map. An angular resolution of 0.4 mas x 0.7 mas is achieved at 8.4 GHz, corresponding to a linear scale of less than 0.013 pc. Hence, we obtain the highest resolution VLBI image of Cen A, comparable to previous space-VLBI observations. By combining with the 22.3 GHz image, which has been taken without contributing transoceanic baselines at somewhat lower resolution, we present the corresponding dual-frequency spectral index distribution along the sub-parsec scale jet revealing the putative emission regions for recently detected gamma-rays from the core region by Fermi/LAT. We resolve the innermost structure of the milliarcsecond scale jet and counterjet system of Cen A into discrete components. The simultaneous observations at two frequencies provide the highest resolved spectral index map of an AGN jet allowing us to identify multiple possible sites as the origin of the high energy emission.