We present a detailed, three-dimensional hydrodynamic study of the neutrino-driven winds that emerge from the remnant of a neutron star merger. Our simulations are performed with the Newtonian, Eulerian code FISH, augmented by a detailed, spectral neutrino leakage scheme that accounts for heating due to neutrino absorption in optically thin conditions. Consistent with the earlier, two-dimensional study of Dessart et al. (2009), we find that a strong baryonic wind is blown out along the original binary rotation axis within ≈100 milliseconds after the merger. We compute a lower limit on the expelled mass of 3.5 × 10 −3 M , large enough to be relevant for heavy element nucleosynthesis. The physical properties vary significantly between different wind regions. For example, due to stronger neutrino irradiation, the polar regions show substantially larger electron fractions than those at lower latitudes. This has its bearings on the nucleosynthesis: the polar ejecta produce interesting r-process contributions from A ≈ 80 to about 130, while the more neutron-rich, lower-latitude parts produce in addition also elements up to the third r-process peak near A ≈ 195. We also calculate the properties of electromagnetic transients that are powered by the radioactivity in the wind, in addition to the "macronova" transient that stems from the dynamic ejecta. The high-latitude (polar) regions produce UV/optical transients reaching luminosities up to 10 41 erg s −1 , which peak around 1 day in optical and 0.3 days in bolometric luminosity. The lower-latitude regions, due to their contamination with high-opacity heavy elements, produce dimmer and more red signals, peaking after ∼ 2 days in optical and infrared. Our numerical experiments indicate that it will be difficult to infer the collapse time-scale of the hypermassive neutron star to a black hole based on the wind electromagnetic transient, at least for collapse time-scales larger than the wind production time-scale.

Massive stars end their lives in explosions with kinetic energies on the order of 10 51 erg. Immediately after the explosion has been launched, a region of low density and high entropy forms behind the ejecta, which is continuously subject to neutrino heating. The neutrinos emitted from the remnant at the center, the protoneutron star (PNS), heat the material above the PNS surface. This heat is partly converted into kinetic energy, and the material accelerates to an outflow that is known as the neutrino-driven wind. For the first time we simulate the collapse, bounce, explosion, and the neutrino-driven wind phases consistently over more than 20 s. Our numerical model is based on spherically symmetric general relativistic radiation hydrodynamics using spectral three-flavor Boltzmann neutrino transport. In simulations where no explosions are obtained naturally, we model neutrino-driven explosions for low-and intermediatemass Fe-core progenitor stars by enhancing the charged current reaction rates. In the case of a special progenitor star, the 8.8 M O-Ne-Mg-core, the explosion in spherical symmetry was obtained without enhanced opacities. The post-explosion evolution is in qualitative agreement with static steady-state and parametrized dynamic models of the neutrino-driven wind. On the other hand, we generally find lower neutrino luminosities and mean neutrino energies, as well as a different evolutionary behavior of the neutrino luminosities and mean neutrino energies. The neutrino-driven wind is proton-rich for more than 10 s and the contraction of the PNS differs from the assumptions made for the conditions at the inner boundary in previous neutrino-driven wind studies. Despite the moderately high entropies of about 100 k B /baryon and the fast expansion timescales, the conditions found in our models are unlikely to favor r-process nucleosynthesis. The simulations are carried out until the neutrino-driven wind settles down to a quasi-stationary state. About 5 s after the bounce, the peak temperature inside the PNS already starts to decrease because of the continued deleptonization. This moment determines the beginning of a cooling phase dominated by the emission of neutrinos. We discuss the physical conditions of the quasi-static PNS evolution and take the effects of deleptonization and mass accretion from early fallback into account.

We examine magnetorotationally driven supernovae as sources of r-process elements in the early Galaxy. On the basis of thermodynamic histories of tracer particles from a three-dimensional magnetohydrodynamical core-collapse supernova model with approximated neutrino transport, we perform nucleosynthesis calculations with and without considering the effects of neutrino absorption reactions on the electron fraction (Y e ) during post-processing. We find that the peak distribution of Y e in the ejecta is shifted from ∼ 0.15 to ∼ 0.17 and broadened toward higher Y e due to neutrino absorption. Nevertheless, in both cases the second and third peaks of the solar r-process element distribution can be well reproduced. The rare progenitor configuration that was used here, characterized by a high rotation rate and a large magnetic field necessary for the formation of bipolar jets, could naturally provide a site for the strong r-process in agreement with observations of the early galactic chemical evolution.

Accurate neutrino transport has been built into spherically symmetric simulations of stellar core collapse and postbounce evolution. The results of such simulations agree that spherically symmetric models with standard microphysical input fail to explode by the delayed, neutrino-driven mechanism. Independent groups implemented fundamentally different numerical methods to tackle the Boltzmann neutrino transport equation. Here we present a direct and detailed comparison of such neutrino radiation-hydrodynamical simulations for two codes, agile-boltztran of the Oak Ridge-Basel group and vertex of the Garching group. The former solves the Boltzmann equation directly by an implicit, general relativistic discrete angle method on the adaptive grid of a conservative implicit hydrodynamics code with second-order TVD advection. In contrast, the latter couples a variable Eddington factor technique with an explicit, moving-grid, conservative high-order Riemann solver with important relativistic effects treated by an effective gravitational potential. The presented study is meant to test both neutrino radiation-hydrodynamics implementations and to provide a data basis for comparisons and verifications of supernova codes to be developed in the future. Results are discussed for simulations of the core collapse and post-bounce evolution of a 13 M ⊙ star with Newtonian gravity and a 15 M ⊙ star with relativistic gravity.