An official website of the United States government
Abstract We explore the effects of rapid rotation on the properties of neutrino-heated winds from proto-neutron stars (PNS) formed in core-collapse supernovae or neutron-star mergers by means of three-dimensional general-relativistic hydrodynamical simulations with M0 neutrino transport. We focus on conditions characteristic of a few seconds into the PNS cooling evolution when the neutrino luminosities obey erg s−1, and over which most of the wind mass loss will occur. After an initial transient phase, all of our models reach approximately steady-state outflow solutions with positive energies and sonic surfaces captured on the computational grid. Our nonrotating and slower rotating models (angular velocity relative to Keplerian Ω/ΩK≲ 0.4; spin periodP≳ 2 ms) generate approximately spherically symmetric outflows with properties in good agreement with previous PNS wind studies. By contrast, our most rapidly spinning PNS solutions (Ω/ΩK≳ 0.75;P≈ 1 ms) generate outflows focused in the rotational equatorial plane with much higher mass-loss rates (by over an order of magnitude), lower velocities, lower entropy, and lower asymptotic electron fractions, than otherwise similar nonrotating wind solutions. Although such rapidly spinning PNS are likely rare in nature, their atypical nucleosynthetic composition and outsized mass yields could render them important contributors of light neutron-rich nuclei compared to more common slowly rotating PNS birth. Our calculations pave the way to including the combined effects of rotation and a dynamically important large-scale magnetic field on the wind properties within a three-dimensional GRMHD framework.
more »
« less
Three-dimensional General-relativistic Simulations of Neutrino-driven Winds from Magnetized Proto–Neutron Stars
Abstract Formed in the aftermath of a core-collapse supernova or neutron star merger, a hot proto–neutron star (PNS) launches an outflow driven by neutrino heating lasting for up to tens of seconds. Though such winds are considered potential sites for the nucleosynthesis of heavy elements via the rapid neutron capture process (r-process), previous work has shown that unmagnetized PNS winds fail to achieve the necessary combination of high entropy and/or short dynamical timescale in the seed nucleus formation region. We present three-dimensional general-relativistic magnetohydrodynamical simulations of PNS winds which include the effects of a dynamically strong (B≳ 1015G) dipole magnetic field. After initializing the magnetic field, the wind quickly develops a helmet-streamer configuration, characterized by outflows along open polar magnetic field lines and a “closed” zone of trapped plasma at lower latitudes. Neutrino heating within the closed zone causes the thermal pressure of the trapped material to rise in time compared to the polar outflow regions, ultimately leading to the expulsion of this matter from the closed zone on a timescale of ∼60 ms, consistent with the predictions of Thompson. The high entropies of these transient ejecta are still growing at the end of our simulations and are sufficient to enable a successful second-peakr-process in at least a modest ≳1% of the equatorial wind ejecta.
more »
« less
- Award ID(s):
- 2002577
- PAR ID:
- 10458159
- Publisher / Repository:
- DOI PREFIX: 10.3847
- Date Published:
- Journal Name:
- The Astrophysical Journal
- Volume:
- 954
- Issue:
- 2
- ISSN:
- 0004-637X
- Format(s):
- Medium: X Size: Article No. 192
- Size(s):
- Article No. 192
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract The remnant black hole–accretion disk system resulting from binary neutron star mergers has proven to be a promising site for synthesizing the heaviest elements via rapid neutron capture (r-process). A critical factor in determining the fullr-process pattern in these environments is the neutron richness of the ejecta, which is strongly influenced by neutrino interactions. One key ingredient shaping these interactions is fast neutrino flavor conversions (FFCs), which arise due to angular crossings in neutrino distributions and occur on nanosecond timescales. We present the first three-dimensional in situ angle-dependent modeling of FFCs in postmerger disks, implemented within general relativistic magnetohydrodynamics with Monte Carlo neutrino transport. Our results reveal that, by suppressing electron neutrinos, FFCs more efficiently cool the disk and weaken the early thermally driven wind. Less releptonization due to electron neutrino absorption makes this cooler wind more neutron rich, producing a more robustr-process at higher latitudes of the outflow. This study underscores the necessity of incorporating FFCs in realistic simulations.more » « less
-
Abstract We present a 3D general-relativistic magnetohydrodynamic simulation of a short-lived neutron star remnant formed in the aftermath of a binary neutron star merger. The simulation uses an M1 neutrino transport scheme to track neutrino–matter interactions and is well suited to studying the resulting nucleosynthesis and kilonova emission. A magnetized wind is driven from the remnant and ejects neutron-rich material at a quasi-steady-state rate of 0.8 × 10−1M⊙s−1. We find that the ejecta in our simulations underproducer-process abundances beyond the secondr-process peak. For sufficiently long-lived remnants, these outflowsalonecan produce blue kilonovae, including the blue kilonova component observed for AT2017gfo.more » « less
-
ABSTRACT Neutron star merger accretion discs can launch neutron-rich winds of >10−2M⊙. This ejecta is a prime site for r-process nucleosynthesis, which will produce a range of radioactive heavy nuclei. The decay of these nuclei releases enough energy to accelerate portions of the wind by ∼0.1c. Here, we investigate the effect of r-process heating on the dynamical evolution of disc winds. We extract the wind from a 3D general relativistic magnetohydrodynamic simulation of a disc from a post-merger system. This is used to create inner boundary conditions for 2D hydrodynamic simulations that continue the original 3D simulation. We perform two such simulations: one that includes the r-process heating, and another one that does not. We follow the hydrodynamic simulations until the winds reach homology (60 s). Using time-dependent multifrequency multidimensional Monte Carlo radiation transport simulations, we then calculate the kilonova light curves from the winds with and without dynamical r-process heating. We find that the r-process heating can substantially alter the velocity distribution of the wind, shifting the mass-weighted median velocity from 0.06c to 0.12c. The inclusion of the dynamical r-process heating makes the light curve brighter and bluer at $$\sim 1\, \mathrm{d}$$ post-merger. However, the high-velocity tail of the ejecta distribution and the early ($$\lesssim 1\, \mathrm{d}$$) light curves are largely unaffected.more » « less
-
ABSTRACT We investigate r-process nucleosynthesis and kilonova emission resulting from binary neutron star (BNS) mergers based on a three-dimensional (3D) general-relativistic magnetohydrodynamic (GRMHD) simulation of a hypermassive neutron star (HMNS) remnant. The simulation includes a microphysical finite-temperature equation of state (EOS) and neutrino emission and absorption effects via a leakage scheme. We track the thermodynamic properties of the ejecta using Lagrangian tracer particles and determine its composition using the nuclear reaction network SkyNet. We investigate the impact of neutrinos on the nucleosynthetic yields by varying the neutrino luminosities during post-processing. The ejecta show a broad distribution with respect to their electron fraction Ye, peaking between ∼0.25–0.4 depending on the neutrino luminosity employed. We find that the resulting r-process abundance patterns differ from solar, with no significant production of material beyond the second r-process peak when using luminosities recorded by the tracer particles. We also map the HMNS outflows to the radiation hydrodynamics code SNEC and predict the evolution of the bolometric luminosity as well as broadband light curves of the kilonova. The bolometric light curve peaks on the timescale of a day and the brightest emission is seen in the infrared bands. This is the first direct calculation of the r-process yields and kilonova signal expected from HMNS winds based on 3D GRMHD simulations. For longer-lived remnants, these winds may be the dominant ejecta component producing the kilonova emission.more » « less
