We present the new code NADA-FLD to solve multi-dimensional neutrino-hydrodynamics in full general relativity (GR) in spherical polar coordinates. The energy-dependent neutrino transport assumes the flux-limited diffusion (FLD) approximation and evolves the neutrino energy densities measured in the frame comoving with the fluid. Operator splitting is used to avoid multi-dimensional coupling of grid cells in implicit integration steps involving matrix inversions. Terms describing lateral diffusion and advection are integrated explicitly using the Allen-Cheng or the Runge-Kutta-Legendre method, which remain stable even in the optically thin regime. We discuss several toy-model problems in one and two dimensions to test the basic functionality and individual components of the transport scheme. We also perform fully dynamic core-collapse supernova (CCSN) simulations in spherical symmetry. For a Newtonian model we find good agreement with the M1 code ALCAR, and for a GR model we reproduce the main effects of GR in CCSNe already found by previous works.

A set of eight self-consistent, time-dependent supernova (SN) simulations in three spatial dimensions (3D) for 9 solar-mass and 20 solar-mass progenitors is evaluated for the presence of dipolar asymmetries of the electron lepton-number emission as discovered by Tamborra et al. and termed lepton-number emission self-sustained asymmetry (LESA). The simulations were performed with the Aenus-Alcar neutrino/hydrodynamics code, which treats the energy- and velocity-dependent transport of neutrinos of all flavors by a two-moment scheme with algebraic M1 closure. For each of the progenitors, results with fully multi-dimensional (FMD) neutrino transport and with ray-by-ray-plus (RbR+) approximation are considered for two different grid resolutions. While the 9 solar-mass models develop explosions, the 20 solar-mass progenitor does not explode with the employed version of simplified neutrino opacities. In all 3D models we observe the growth of substantial dipole amplitudes of the lepton-number (electron neutrino minus antineutrino) flux with stable or slowly time-evolving direction and overall properties fully consistent with the LESA phenomenon. Models with RbR+ transport develop LESA dipoles somewhat faster and with temporarily higher amplitudes, but the FMD calculations exhibit cleaner hemispheric asymmetries with a far more dominant dipole. In contrast, the RbR+ results display much wider multipole spectra of the neutrino-emission anisotropies with significant power also in the quadrupole and higher-order modes. Our results disprove speculations that LESA is a numerical artifact of RbR+ transport. We also discuss LESA as consequence of a dipolar convection flow inside of the nascent neutron star and establish, tentatively, a connection to Chandrasekhar's linear theory of thermal instability in spherical shells.

We describe a multi-messenger interpretation of GW170817, which yields a robust lower limit on NS radii. This excludes NSs with radii smaller than about 10.7 km and thus rules out very soft nuclear matter. We stress the potential of this type of constraints when future detections become available. A very similar argumentation may yield an upper bound on the maximum mass of nonrotating NSs. We also discuss simulations of NS mergers, which undergo a first-order phase transition to quark matter. We point out a different dynamical behavior. Considering the gravitational-wave signal, we identify an unambiguous signature of the QCD phase transition in NS mergers. The occurrence of quark matter through a strong first-order phase transition during merging leads to a characteristic shift of the dominant postmerger frequency. The frequency shift is indicative for a phase transition if it is compared to the postmerger frequency which is expected for purely hadronic EoS models. A very strong deviation of several 100 Hz is observed for hybrid EoSs in an otherwise tight relation between the tidal deformability and the postmerger frequency. We address the potential impact of a first-order phase transition on the electromagnetic counterpart of NS mergers. Our simulations suggest that there would be no significant qualitative differences between a system undergoing a phase transition to quark matter and purely hadronic mergers. The quantitative differences are within the spread which is found between different hadronic EoS models. This implies on the one hand that GW170817 is compatible with a possible transition to quark matter. On the other hand these considerations show that it may not be easy to identify quantitative differences between purely hadronic mergers and events in which quark matter occurs considering solely their electromagnetic counterpart or their nucleosynthesis products. (abridged)

We present a new, computationally efficient, energy-integrated approximation for neutrino effects in hot and dense astrophysical environments such as supernova cores and compact binary mergers and their remnants. Our new method, termed ILEAS for Improved Leakage-Equilibration-Absorption Scheme, improves the lepton-number and energy losses of traditional leakage descriptions by a novel prescription of the diffusion time-scale based on a detailed energy integral of the flux-limited diffusion equation. The leakage module is supplemented by a neutrino-equilibration treatment that ensures the proper evolution of the total lepton number and medium plus neutrino energies as well as neutrino-pressure effects in the neutrino-trapping domain. Moreover, we employ a simple and straightforwardly applicable ray-tracing algorithm for including re-absorption of escaping neutrinos especially in the decoupling layer and during the transition to semi-transparent conditions. ILEAS is implemented on a three-dimensional (3D) Cartesian grid with a minimum of free and potentially case-dependent parameters and exploits the basic physics constraints that should be fulfilled in the neutrino-opaque and free-streaming limits. We discuss a suite of tests for stationary and time-dependent proto-neutron star models and post-merger black-hole-torus configurations, for which 3D ILEAS results are demonstrated to agree with energy-dependent 1D and 2D two-moment (M1) neutrino transport on the level of 10--15 percent in basic neutrino properties. This also holds for the radial profiles of the neutrino luminosities and of the electron fraction. Even neutrino absorption maps around torus-like neutrino sources are qualitatively similar without any fine-tuning, confirming that ILEAS can satisfactorily reproduce local losses and re-absorption of neutrinos as found in sophisticated transport calculations.

Self-consistent, time-dependent supernova (SN) simulations in three spatial dimensions (3D) are conducted with the Aenus-Alcar code, comparing, for the first time, calculations with fully multi-dimensional (FMD) neutrino transport and the ray-by-ray-plus (RbR+) approximation, both based on a two-moment solver with algebraic M1 closure. We find good agreement between 3D results with FMD and RbR+ transport for both tested grid resolutions in the cases of a 20 solar-mass progenitor, which does not explode with the employed simplified set of neutrino opacities, and of an exploding 9 solar-mass model. This is in stark contrast to corresponding axisymmetric (2D) simulations, which confirm previous claims that the RbR+ approximation can foster explosions in 2D in particular in models with powerful axial sloshing of the stalled shock due to the standing accretion shock instability (SASI). However, while local and instantaneous variations of neutrino fluxes and heating rates can still be considerably higher with RbR+ transport in 3D, the time-averaged quantities are very similar to FMD results because of the absence of a fixed, artificial symmetry axis that channels the flow. Therefore, except for stochastic fluctuations, the neutrino signals and the post-bounce evolution of 3D simulations with FMD and RbR+ transport are also very similar, in particular for our calculations with the better grid resolution. Higher spatial resolution has clearly a more important impact than the differences by the two transport treatments. Our results back up the use of the RbR+ approximation for neutrino transport in 3D SN modeling.

We present spherically symmetric (1D) and axisymmetric (2D) supernova simulations for a convection-dominated 9 Msun and a 20 Msun progenitor that develops violent activity by the standing-accretion-shock instability (SASI). We compare in detail the Aenus-Alcar code, which uses fully multidimensional two-moment neutrino transport with an M1 closure, with a ray-by-ray-plus (RbR+) version of this code and with the Prometheus-Vertex code that employs RbR+ two-moment transport with a Boltzmann closure. Besides testing consequences of ignored non-radial neutrino-flux components in the RbR+ approximation, we also discuss the influence of various transport ingredients applied or not applied in recent literature, namely simplified neutrino-pair processes, neutrino-electron scattering, velocity-dependent and gravitational-redshift terms, and strangeness and many-body corrections for neutrino-nucleon scattering. Alcar and Vertex show excellent agreement in 1D and 2D despite a slightly but systematically smaller radius (~1km) and stronger convection of the proto-neutron star with Alcar. As found previously, the RbR+ approximation is conducive to explosions, but much less severely in the convection-dominated 9 Msun case than in the marginally exploding 20 Msun model, where the onset time of explosion also exhibits big stochastic variations, and the RbR+ approximation has no distinctly stronger supportive effect than simplified pair processes or strangeness and many-body corrections. Neglecting neutrino-electron scattering has clearly unfavorable effects for explosions, while ignoring velocity and gravitational-redshift effects can both promote or delay the explosion. The ratio of advection timescale to neutrino-heating timescale in 1D simulations is a sensitive indicator of the influence of physics ingredients on explosions also in multidimensional simulations.

We study the effects of magnetic fields and rotation on the core collapse of a star of an initial mass of M = 20 solar masses using axisymmetric simulations coupling special relativistic magnetohydrodynamics, an approximately relativistic gravitational potential, and spectral neutrino transport. We compare models of the same core with different, artificially added profiles of rotation and magnetic field. A model with weak field and slow rotation does not produce an explosion, while stronger fields and fast rotation open the possibility of explosions. Whereas the neutrino luminosities of the exploding models are the same as or even less than those of the non-exploding model, magnetic fields locally in equipartition with the gas pressure provide a strong contribution to the shock revival and the acceleration of bipolar outflows. Among the amplification processes generating such strong fields, we find the magneto-rotational instability. However, our limited grid resolution allows us to find it only in limited regions of the model with the strongest pre-collapse field ($10^{11}$ G) and fastest rotation.

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.

The remnant of neutron star mergers is dense in neutrinos. By employing inputs from one hydrodynamical simulation of a binary neutron star merger remnant with a black hole of $3\ M_\odot$ in the center, dimensionless spin parameter $0.8$ and an accretion torus of $0.3\ M_\odot$, the neutrino emission properties are investigated as the merger remnant evolves. Initially, the local number density of $\bar{\nu}_e$ is larger than that of $\nu_e$ everywhere above the remnant. Then, as the torus approaches self-regulated equilibrium, the local abundance of neutrinos overcomes that of antineutrinos in a funnel around the polar region. The region where the fast pairwise flavor conversions can occur shrinks accordingly as time evolves. Still, we find that fast flavor conversions do affect most of the neutrino-driven ejecta. Assuming that fast flavor conversions lead to flavor equilibration, a significant enhancement of nuclei with mass numbers $A>130$ is found as well as a change of the lanthanide mass fraction by more than a factor of a thousand. Our findings hint towards a potentially relevant role of neutrino flavor oscillations for the prediction of the kilonova (macronova) lightcurves and motivate further work in this direction.

We introduce a new, powerful method to constrain properties of neutron stars (NSs). We show that the total mass of GW170817 provides a reliable constraint on the stellar radius if the merger did not result in a prompt collapse as suggested by the interpretation of associated electromagnetic emission. The radius R_1.6 of nonrotating NSs with a mass of 1.6 M_sun can be constrained to be larger than 10.68_{-0.04}^{+0.15} km, and the radius R_max of the nonrotating maximum mass configuration must be larger than 9.60_{-0.03}^{+0.14} km. We point out that detections of future events will further improve these constraints. Moreover, we show that a future event with a signature of a prompt collapse of the merger remnant will establish even stronger constraints on the NS radius from above and the maximum mass M_max of NSs from above. These constraints are particularly robust because they only require a measurement of the chirp mass and a distinction between prompt and delayed collapse of the merger remnant, which may be inferred from the electromagnetic signal or even from the presence/absence of a ringdown gravitational-wave (GW) signal. This prospect strengthens the case of our novel method of constraining NS properties, which is directly applicable to future GW events with accompanying electromagnetic counterpart observations. We emphasize that this procedure is a new way of constraining NS radii from GW detections independent of existing efforts to infer radius information from the late inspiral phase or postmerger oscillations, and it does not require particularly loud GW events.

The merger of two neutron stars may give birth to a long-lived hypermassive neutron star. If it harbours a strong magnetic field of magnetar strength, its spin-down could explain several features of short gamma-ray burst afterglows. The magnetorotational instability (MRI) has been proposed as a mechanism to amplify the magnetic field to the required strength. Previous studies have, however, neglected neutrinos, which may have an important impact on the MRI by inducing a viscosity and drag. We investigate the impact of these neutrino effects on the linear growth of the MRI by applying a local stability analysis to snapshots of a neutron star merger simulation. We find that neutrinos have a significant impact inside the hypermassive neutron star, but have at most a marginal effect in the torus surrounding it. Inside the hypermassive neutron star, the MRI grows in different regimes depending on the radius and on the initial magnetic field strength. For magnetic fields weaker than $10^{13}-10^{14}\,{\rm G}$, the growth rate of the MRI is significantly reduced due to the presence of neutrinos. We conclude that neutrinos should be taken into account when studying the growth of the MRI from realistic initial magnetic fields. Current numerical simulations, which neglect neutrino viscosity, are only consistent, i.e. in the adopted ideal regime, if they start from artificially strong initial magnetic fields above $\sim10^{14}\,{\rm G}$. One should be careful when extrapolating these results to lower initial magnetic fields, where the MRI growth is strongly affected by neutrino viscosity or drag.

We present the first special relativistic, axisymmetric hydrodynamic simulations of black hole-torus systems (approximating general relativistic gravity) as remnants of binary-neutron star (NS-NS) and neutron star-black hole (NS-BH) mergers, in which the viscously driven evolution of the accretion torus is followed with self-consistent energy-dependent neutrino transport and the interaction with the cloud of dynamical ejecta expelled during the NS-NS merging is taken into account. The modeled torus masses, BH masses and spins, and the ejecta masses, velocities, and spatial distributions are adopted from relativistic merger simulations. We find that energy deposition by neutrino annihilation can accelerate outflows with initially high Lorentz factors along polar low-density funnels, but only in mergers with extremely low baryon pollution in the polar regions. NS-BH mergers, where polar mass ejection during the merging phase is absent, provide sufficiently baryon-poor environments to enable neutrino-powered, ultrarelativistic jets with terminal Lorentz factors above 100 and considerable dynamical collimation, favoring short gamma-ray bursts (sGRBs), although their typical energies and durations might be too small to explain the majority of events. In the case of NS-NS mergers, however, neutrino emission of the accreting and viscously spreading torus is too short and too weak to yield enough energy for the outflows to break out from the surrounding ejecta shell as highly relativistic jets. We conclude that neutrino annihilation alone cannot power sGRBs from NS-NS mergers.

We present a comprehensive study of r-process element nucleosynthesis in the ejecta of compact binary mergers (CBMs) and their relic black-hole (BH)-torus systems. The evolution of the BH-accretion tori is simulated for seconds with a Newtonian hydrodynamics code including viscosity effects, pseudo-Newtonian gravity for rotating BHs, and an energy-dependent two-moment closure scheme for the transport of electron neutrinos and antineutrinos. The investigated cases are guided by relativistic double neutron star (NS-NS) and NS-BH merger models, producing ~3-6 Msun BHs with rotation parameters of A~0.8 and tori of 0.03-0.3 Msun. Our nucleosynthesis analysis includes the dynamical (prompt) ejecta expelled during the CBM phase and the neutrino and viscously driven outflows of the relic BH-torus systems. While typically ~20-25% of the initial accretion-torus mass are lost by viscously driven outflows, neutrino-powered winds contribute at most another ~1%, but neutrino heating enhances the viscous ejecta significantly. Since BH-torus ejecta possess a wide distribution of electron fractions (0.1-0.6) and entropies, they produce heavy elements from A~80 up to the actinides, with relative contributions of A>130 nuclei being subdominant and sensitively dependent on BH and torus masses and the exact treatment of shear viscosity. The combined ejecta of CBM and BH-torus phases can reproduce the solar abundances amazingly well for A>90. Varying contributions of the torus ejecta might account for observed variations of lighter elements with 40<Z<56 relative to heavier ones, and a considerable reduction of the prompt ejecta compared to the torus ejecta, e.g. in highly asymmetric NS-BH mergers, might explain the composition of heavy-element deficient stars.

We present the new code ALCAR developed to model multidimensional, multi energy-group neutrino transport in the context of supernovae and neutron-star mergers. The algorithm solves the evolution equations of the 0th- and 1st-order angular moments of the specific intensity, supplemented by an algebraic relation for the 2nd-moment tensor to close the system. The scheme takes into account frame-dependent effects of order O(v/c) as well as the most important types of neutrino interactions. The transport scheme is significantly more efficient than a multidimensional solver of the Boltzmann equation, while it is more accurate and consistent than the flux-limited diffusion method. The finite-volume discretization of the essentially hyperbolic system of moment equations employs methods well-known from hydrodynamics. For the time integration of the potentially stiff moment equations we employ a scheme in which only the local source terms are treated implicitly, while the advection terms are kept explicit, thereby allowing for an efficient computational parallelization of the algorithm. We investigate various problem setups in one and two dimensions to verify the implementation and to test the quality of the algebraic closure scheme. In our most detailed test, we compare a fully dynamic, one-dimensional core-collapse simulation with two published calculations performed with well-known Boltzmann-type neutrino-hydrodynamics codes and we find very satisfactory agreement.

We investigate beta-interactions of free nucleons and their impact on the electron fraction (Y_e) and r-process nucleosynthesis in ejecta characteristic of binary neutron star mergers (BNSMs). For that we employ trajectories from a relativistic BNSM model to represent the density-temperature evolutions in our parametric study. In the high-density environment, positron captures decrease the neutron richness at the high temperatures predicted by the hydrodynamic simulation. Circumventing the complexities of modelling three-dimensional neutrino transport, (anti)neutrino captures are parameterized in terms of prescribed neutrino luminosities and mean energies, guided by published results and assumed as constant in time. Depending sensitively on the adopted neutrino-antineutrino luminosity ratio, neutrino processes increase Y_e to values between 0.25 and 0.40, still allowing for a successful r-process compatible with the observed solar abundance distribution and a significant fraction of the ejecta consisting of r-process nuclei. If the electron neutrino luminosities and mean energies are relatively large compared to the antineutrino properties, the mean Y_e might reach values >0.40 so that neutrino captures seriously compromise the success of the r-process. In this case, the r-abundances remain compatible with the solar distribution, but the total amount of ejected r-material is reduced to a few percent, because the production of iron-peak elements is favored. Proper neutrino physics, in particular also neutrino absorption, have to be included in BNSM simulations before final conclusions can be drawn concerning r-processing in this environment and concerning observational consequences like kilonovae, whose peak brightness and color temperature are sensitive to the composition-dependent opacity of the ejecta.

We present the first comprehensive study of r-process element nucleosynthesis in the ejecta of compact binary mergers (CBMs) and their relic black-hole (BH)-torus systems. The evolution of the BH-accretion tori is simulated for seconds with a Newtonian hydrodynamics code including viscosity effects, pseudo-Newtonian gravity for rotating BHs, and an energy-dependent two-moment closure scheme for the transport of electron neutrinos and antineutrinos. The investigated cases are guided by relativistic double neutron star (NS-NS) and NS-BH merger models, producing ~3-6 Msun BHs with rotation parameters of A~0.8 and tori of 0.03-0.3 Msun. Our nucleosynthesis analysis includes the dynamical (prompt) ejecta expelled during the CBM phase and the neutrino and viscously driven outflows of the relic BH-torus systems. While typically ~20-25% of the initial accretion-torus mass are lost by viscously driven outflows, neutrino-powered winds contribute at most another ~1%, but neutrino heating enhances the viscous ejecta significantly. Since BH-torus ejecta possess a wide distribution of electron fractions (0.1-0.6) and entropies, they produce heavy elements from A~80 up to the actinides, with relative contributions of A>130 nuclei being subdominant and sensitively dependent on BH and torus masses and the exact treatment of shear viscosity. The combined ejecta of CBM and BH-torus phases can reproduce the solar abundances amazingly well for A>90. Varying contributions of the torus ejecta might account for observed variations of lighter elements with 40<Z<56 relative to heavier ones, and a considerable reduction of the prompt ejecta compared to the torus ejecta, e.g. in highly asymmetric NS-BH mergers, might explain the composition of heavy-element deficient stars.

The likely outcome of a compact object merger event is a central black hole surrounded by a rapidly accreting torus of debris. This disk of debris is a rich source of element synthesis, the outcome of which is needed to predict electromagnetic counterparts of individual events and to understand the contribution of mergers to galactic chemical evolution. Here we study disk outflow nucleosynthesis in the context of a two-dimensional, time-dependent black hole-neutron star merger accretion disk model. We use two time snapshots from this model to examine the impact of the evolution of the neutrino fluxes from the disk on the element synthesis. While the neutrino fluxes from the early-time disk snapshot appear to favor neutron-rich outflows, by the late-time snapshot the situation is reversed. As a result we find copious production of Nickel-56 in the outflows.