Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Ejecta”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

The multiwavelength view of shocks in the fastest nova V1674 Her

Classical novae are shock-powered multiwavelength transients triggered by a thermonuclear runaway on an accreting white dwarf. V1674 Her is the fastest nova ever recorded (time to declined by two magnitudes is t 2 = 1.1 d) that challenges our understanding of shock formation in novae. We investigate the physical mechanisms behind nova emission from GeV γ-rays to cm-band radio using coordinated Fermi-LAT, NuSTAR, Swift, and VLA observations supported by optical photometry. Fermi-LAT detected short-lived (18 h) 0.1–100 GeV emission from V1674 Her that appeared 6 h after the eruption began; this was at a level of (1.6 ± 0.4) × 10 −6 photons cm −2 s −1 . Eleven days later, simultaneous NuSTAR and Swift X-ray observations revealed optically thin thermal plasma shock-heated to k $T$ shock = 4 keV. The lack of a detectable 6.7 keV Fe Kα emission suggests super-solar CNO abundances. The radio emission from V1674 Her was consistent with thermal emission at early times and synchrotron at late times. The radio spectrum steeply rising with frequency may be a result of either free-free absorption of synchrotron and thermal emission by unshocked outer regions of the nova shell or the Razin–Tsytovich effect attenuating synchrotron emission in dense plasma. The development of the shock inside the ejecta is unaffected by the extraordinarily rapid evolution and the intermediate polar host of this nova.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

How negative feedback and the ambient environment limit the influence of recombination in common envelope evolution

ABSTRACT We perform 3D hydrodynamical simulations to study recombination and ionization during the common envelope (CE) phase of binary evolution, and develop techniques to track the ionic transitions in time and space. We simulate the interaction of a $2\, \mathrm{M_\odot }$ red giant branch primary and a $1\, \mathrm{M_\odot }$ companion modelled as a particle. We compare a run employing a tabulated equation of state (EOS) that accounts for ionization and recombination, with a run employing an ideal gas EOS. During the first half of the simulations, ∼15 per cent more mass is unbound in the tabulated EOS run due to the release of recombination energy, but by simulation end the difference has become negligible. We explain this as being a consequence of (i) the tabulated EOS run experiences a shallower inspiral and hence smaller orbital energy release at late times because recombination energy release expands the envelope and reduces drag, and (ii) collision and mixing between expanding envelope gas, ejecta and circumstellar ambient gas assists in unbinding the envelope, but does so less efficiently in the tabulated EOS run where some of the energy transferred to bound envelope gas is used for ionization. The rate of mass unbinding is approximately constant in the last half of the simulations and the orbital separation steadily decreases at late times. A simple linear extrapolation predicts a CE phase duration of ${\sim}2\, {\rm yr}$, after which the envelope would be unbound.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

What if GW190425 did not produce a black hole promptly?

ABSTRACT It is widely believed that the binary neutron star merger GW190425 produced a black hole promptly upon merger. Motivated by the potential association with the fast radio burst FRB 20190425A, which took place 2.5 h after the merger, we revisit the question of the outcome of GW190425 by means of numerical relativity simulations. We show that current laboratory and astrophysical constraints on the equation of state of dense matter do not rule out the formation of a long-lived remnant. However, the formation of a stable remnant would have produced a bright kilonova, in tension with upper limits by ZTF at the location and time of FRB 20190425A. Moreover, the ejecta would have been optically thick to radio emission for days to months, preventing a putative FRB from propagating out. The predicted dispersion measure is also several orders of magnitude larger than that observed for FRB 20190425A. Our results indicate that FRB 20190425A and GW190425 are not associated. However, we cannot completely rule out the formation of a long-lived remnant, due to the incomplete coverage of the relevant sky regions. More observations of GW190425-like events, including potential upper limit, have the potential to constrain nuclear physics. To this aim, it is important that follow-up observational campaigns of gravitational wave events are informed by the properties of the source, such as their chirp mass, and we urge the LIGO-Virgo-KAGRA collaboration to promptly release them publicly.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

The luminous, slow-rising orphan afterglow AT2019pim as a candidate moderately relativistic outflow

ABSTRACT Classical gamma-ray bursts (GRBs) have two distinct emission episodes: prompt emission from ultrarelativistic ejecta and afterglow from shocked circumstellar material. While both components are extremely luminous in known GRBs, a variety of scenarios predict the existence of luminous afterglow emission with little or no associated high-energy prompt emission. We present AT 2019pim, the first spectroscopically confirmed afterglow with no observed high-energy emission to be identified. Serendipitously discovered during follow-up observations of a gravitational-wave trigger and located in a contemporaneous TESS sector, it is hallmarked by a fast-rising ($t \approx 2$ h), luminous ($M_{\rm UV,peak} \approx -24.4$ mag) optical transient with accompanying luminous X-ray and radio emission. No gamma-ray emission consistent with the time and location of the transient was detected by Fermi-GBM or by Konus, placing constraining limits on an accompanying GRB. We investigate several independent observational aspects of the afterglow in the context of constraints on relativistic motion and find all of them are consistent with an initial Lorentz factor of $\Gamma _0 \approx$ 10–30 for the on-axis material, significantly lower than in any well-observed GRB and consistent with the theoretically predicted ‘dirty fireball’ scenario in which the high-energy prompt emission is stifled by pair production. However, we cannot rule out a structured jet model in which only the line-of-sight material was ejected at low-$\Gamma$, off-axis from a classical high-$\Gamma$ jet core, and an on-axis GRB with below-average gamma-ray efficiency also remains a possibility. This event represents a milestone in orphan afterglow searches, demonstrating that luminous optical afterglows lacking detected GRB counterparts can be identified and spectroscopically confirmed in real time.

Perley, Daniel A. (ORCID:0000000184721996)↗
DOE OSTI2025

Optical evolution of AT 2024wpp: the high-velocity outflows in Cow-like transients are consistent with high spherical symmetry

ABSTRACT We present the analysis of optical/near-infrared (NIR) data and host galaxy properties of a bright, extremely rapidly evolving transient, AT 2024wpp, which resembles the enigmatic AT 2018cow. AT 2024wpp rose to a peak brightness of $c=-21.9$ mag in 4.3 d and remained above the half-maximum brightness for only 6.7 d. The blackbody fits to the photometry show that the event remained persistently hot ($T\gtrsim 20\, 000$ K) with a rapidly receding photosphere ($v\sim 11\, 500$ km s$^{-1}$), similarly to AT 2018cow albeit with a several times larger photosphere. $JH$ photometry reveals an NIR excess over the thermal emission at $\sim +20$ d, indicating a presence of an additional component. The spectra are consistent with blackbody emission throughout our spectral sequence ending at $+21.9$ d, showing a tentative, very broad emission feature at $\sim 5500$ Å – implying that the optical photosphere is likely within a near-relativistic outflow. Furthermore, reports of strong X-ray and radio emission cement the nature of AT 2024wpp as a likely Cow-like transient. AT 2024wpp is the second event of the class with optical polarimetry. Our $BVRI$ observations obtained from $+6.1$ to $+14.4$ d show a low polarization of $P\lesssim 0.5$ per cent across all bands, similar to AT 2018cow that was consistent with $P\sim 0$ per cent during the same outflow-driven phase. In the absence of evidence for a preferential viewing angle, it is unlikely that both events would have shown low polarization in the case that their photospheres were aspherical. As such, we conclude that the near-relativistic outflows launched in these events are likely highly spherical, but polarimetric observations of further events are crucial to constrain their ejecta geometry and stratification in detail.

Pursiainen, M.↗
DOE OSTI2026

Collapse of magnetized white dwarfs as site of heavy-element formation and kilonova signal

We present the first end-to-end calculation connecting the accretion-induced collapse (AIC) of a magnetized, rapidly rotating white dwarf to observable kilonova signatures, combining two-dimensional (2D) general-relativistic neutrino-magnetohydrodynamic simulations, followed by radiation hydrodynamics with in-situ nuclear network and 2D Monte Carlo radiative transfer with spatially resolved heating rates. Unlike all previous unmagnetized AIC models – which predicted proton-rich, $^{56}$Ni-dominated ejecta – strong magnetic fields eject ${\approx }\, 0.2\, \mathrm{ M}_\odot$ of neutron-rich material ($\langle Y_e \rangle \sim 0.24$) on dynamical time-scales, before neutrino irradiation can raise the electron fraction, enabling strong r-process nucleosynthesis up to and beyond the third peak. The resulting kilonova is lanthanide-rich ($X_{\rm lan} \approx 8~{{\ \rm per\ cent}}$) and dominated by near-infrared emission. We compute synthetic light curves in the Large Synoptic Survey Telescope and J ames Webb Space Telescope bands and find striking agreement, without parameter tuning, between the observations of AT 2023vfi/GRB 230307A and our broadband light curves for polar viewing angles. These results establish magnetized AIC as a viable channel for heavy r-process element production and a compelling progenitor candidate for long-duration gamma-ray bursts with kilonova signatures.

MHD↗
DOE OSTI2020

Type Ia supernovae deflagration-to-detonation transition explosions powered by the Zel’dovich reactivity gradient mechanism

ABSTRACT Our aim in this work is to identify and explain the necessary conditions required for an energetic explosion of a Chandrasekhar-mass white dwarf. We construct and analyse weakly compressible turbulence models with nuclear burning effects for carbon/oxygen plasma at a density expected for the deflagration-to-detonation transition (DDT) to occur. We observe the formation of carbon deflagrations and transient carbon detonations at early times. As turbulence becomes increasingly inhomogeneous, sustained carbon detonations are initiated by the Zel’dovich reactivity gradient mechanism. The fuel is suitably preconditioned by the action of compressive turbulent modes with wavelength comparable to the size of resolved turbulent eddies; no acoustic wave is involved in this process. Oxygen detonations are initiated, aided either by reactivity gradients or by collisions of carbon detonations. The observed evolutionary time-scales are found to be sufficiently short for the above process to occur in the expanding, centrally ignited massive white dwarf. The inhomogeneous conditions produced prior to the DDT might be of consequence for the chemical composition of the outer ejecta regions of Type Ia supernovae from the single degenerate channel, and offer the potential for validation of the proposed model.

Brooker, E.↗
DOE OSTI2025

Distinguishing prompt-collapse binary neutron star mergers from binary black Holes: Tidal effects and remnant properties

We study the properties of remnants formed in prompt-collapse binary neutron star mergers. We consider nonspinning neutron star binaries over a range of total masses and mass ratios across a set of 22 equations of state, totaling 107 numerical relativity simulations. We report the final mass and spin of the systems (including the accretion disk and ejecta) to be constrained in a narrow range—0.98 ≲ 𝑀 𝑓 /𝑀 ≲ 0.99 for the mass and 0.85 ≲ 𝑎 𝑓 ≲ 0.95 for the dimensionless spin—regardless of the binary configuration and matter effects. This sets them apart from binary black hole merger remnants. We assess the detectability of the postmerger signal in a future 40 km Cosmic Explorer observatory and find that the signal-to-noise ratio in the postmerger of an optimally located and oriented binary at a distance of 100 Mpc can range from <1 to 8, depending on the binary configuration and equation of state, with a majority of them greater than 4 in the set of simulations that we consider. We also consider the distinguishability between prompt-collapse binary neutron star and binary black hole mergers with the same masses and spins. We find that Cosmic Explorer will be able to distinguish such systems primarily via the measurement of tidal effects in the late inspiral. Neutron star binaries with reduced tidal deformability $\tilde{Λ}$ as small as ∼ 3.5 can be identified up to a distance of 100 Mpc, while neutron star binaries with $\tilde{Λ}$ ∼ 22 can be identified to distances greater than 250 Mpc. This is larger than the distance up to which the postmerger will be visible. Finally, we discuss the possible implications of our findings for the equation of state of neutron stars from the gravitational wave event GW230529.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Neutrino Flavor Transformation in Neutron Star Mergers

Here, we present the first numerical relativity simulations including neutrino flavor transformations that could result from flavor instabilities, quantum many-body effects, or potential beyond standard model physics in neutron star mergers. We find that neutrino flavor transformations impact the composition and structure of the remnant, potentially leaving an imprint on the postmerger gravitational-wave signal. They also have a significant impact on the composition and nucleosynthesis yields of the ejecta.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Efficient method for estimation of fission fragment yields of $\textit{r}$-process nuclei

Background: More than half of all the elements heavier than iron are made by the rapid neutron capture process (or $\textit{r}$ process). For very-neutron-rich astrophysical conditions, such at those found in the tidal ejecta of neutron stars, nuclear fission determines the $\textit{r}$-process endpoint, and the fission-fragment yields shape the final abundances of 110 ≤ $\textit{A}$ ≤ 170 nuclei. The knowledge of fission-fragment yields of hundreds of nuclei inhabiting very-neutron-rich regions of the nuclear landscape is thus crucial for the modeling of heavy-element nucleosynthesis. Purpose: In this study, we propose a model for the fast calculation of fission-fragment yields based on the concept of shell-stabilized prefragments defined with help of the nucleonic localization functions. Methods: To generate realistic potential-energy surfaces and nucleonic localizations, we apply Skyrme density-functional theory. In this work, the distribution of the neck nucleons among the two prefragments is obtained by means of a statistical model. Results: We benchmark the method by studying the fission yields of 178 Pt, 240 Pu, 254 Cf, and 254,256,258 Fm and show that it satisfactorily explains the experimental data. We then make predictions for 254 Pu and 290 Fm as two representative cases of fissioning nuclei that are expected to significantly contribute during the $\textit{r}$-process nucleosynthesis occurring in neutron-star mergers. Conclusions: The proposed framework provides an efficient alternative to microscopic approaches based on the evolution of the system in a space of collective coordinates all the way to scission. It can be used to carry out global calculations of fission-fragment distributions across the $\textit{r}$-process region.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Shell-model studies of the astrophysical r p -process reactions S 34 ( p , γ ) Cl 35 and Cl 34 g , m ( p , γ ) Ar 35

Background: Dust grains condensed in the outflows of presolar classical novae should have been present in the protosolar nebula. Candidates for such presolar nova grains have been found in primitive meteorites and can in principle be identified by their isotopic ratios, but the ratios predicted by state-of-the-art one-dimensional hydrodynamic models are uncertain due to nuclear-physics uncertainties. Purpose: To theoretically calculate the thermonuclear rates and uncertainties of the S 34 ( p , γ ) Cl 35 and Cl 34 g , m ( p , γ ) Ar 35 reactions and investigate their impacts on the predicted 34 S/ 32 S isotopic ratio for presolar nova grains. Method: A shell-model approach in a (0 +1 ) ℏ ω model space was used to calculate the properties of resonances in the 34 S(p ,γ ) 35 Cl and 34 g,m Cl (p ,γ ) 35 Ar reactions and their thermonuclear rates. Uncertainties were estimated using a Monte Carlo method. The implications of these rates and their uncertainties on sulfur isotopic nova yields were investigated using a postprocessing nucleosynthesis code. The rates for transitions from the ground state of 34 Cl as well as from the isomeric first excited state of 34 Cl were explicitly calculated. Results: At energies in the resonance region near the proton-emission threshold, many negative-parity states appear. Energies, spectroscopic factors, and proton-decay widths are reported. The resulting thermonuclear rates are compared with previous determinations. Conclusions: The shell-model calculations alone are sufficient to constrain the variation of the 34 S/ 32 S ratios to within about 30%. Uncertainties associated with other reactions must also be considered, but in general we find that the 34 S/ 32 S ratios are not a robust diagnostic to clearly identify presolar grains made from nova ejecta.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Developing the S 32 ( p , d ) S * 31 ( p ) ( γ ) reaction to probe the P 30 ( p , γ ) S 31 reaction rate in classical novae

Background: The 30 P(p,γ) 31 S reaction rate is one of the largest remaining sources of uncertainty in the final abundances of nuclei created in a classical nova involving a ONe white dwarf. The reaction rate directly influences silicon isotopic ratios, which are used as identifiers of presolar grains with nova origins. Additionally, the uncertainty in the 30 P(p,γ) 31 S reaction rate has been found to limit the use of nova nuclear thermometers based on observations of elemental ratios in nova ejecta. Purpose: Reduce uncertainties in the nuclear data for proton-unbound states in 31 S, which act as resonances for the 30 P(p,γ) 31 S reaction at classical nova temperatures, and develop a technique for high efficiency, high-resolution reaction-decay coincidence measurements. Methods: Here, the 32 S(p,d) 31 S reaction was used to populate the states of interest in 31 S. The experiment was performed at the Texas A&M Cyclotron Institute using the LLNL Hyperion array for the detection of charged particles and γ rays. A downstream silicon telescope was used to select reaction deuterons, and a single upstream silicon detector was used to measure protons emitted in the decay of unbound 31 S levels. Results: Several states in 31 S above the proton separation energy were observed to have been populated. Decay protons from the resonant states in 31 S were identified as events in the upstream silicon detectors that came in coincidence with deuterons in the downstream telescope. Protons emitted from these states were measured and branching ratios extracted. Conclusions: While no new reaction rate is derived, spin-parity assignments for several higher-lying proton unbound states have been confirmed. Measured p 0 branching ratios for these levels have been compared to previous measurements with good agreement, and in some cases provided a reduction in uncertainty. The previously identified T = 3/2 state may have been incorrectly assigned a large p 0 branching ratio in a previous measurement. The technique of measuring reaction-decay coincidences with a particle-gamma setup appears promising.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

First inverse kinematics measurement of resonances in Be 7 ( α , γ ) C 11 relevant to neutrino-driven wind nucleosynthesis using DRAGON

A possible mechanism to explain the origin of the light $\textit{p}$ nuclei in the Galaxy is the nucleosynthesis in the proton-rich neutrino-driven wind ejecta of core-collapse supernovas via the $\textit{νp}$ process. However, this production scenario is very sensitive to the underlying supernova dynamics and the nuclear physics input. As far as the nuclear uncertainties are concerned, the breakout from the $\textit{pp}$ chains via the 7 Be(α,γ) 11 C reaction has been identified as an important link which can influence the nuclear flow and, therefore, the efficiency of the $\textit{νp}$ process. However, its reaction rate is poorly known over the relevant temperature range, $\textit{T}$ = 1.5–3 GK. We report on the first direct measurement of two resonances of the 7 Be(α,γ) 11 C reaction with previously unknown strengths using an intense radioactive 7 Be beam from the Isotope Separator and Accelerator (ISAC-I) Center facility and the DRAGON recoil separator in inverse kinematics. In this work, we have decreased the 7 Be(α,γ) 11 C reaction rate uncertainty to ≈9.4–10.7% over the relevant temperature region.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Shake-up and shake-off effects in neutrinoless double- β decay

Here, we revisit the role of shake-up and shake-off effects in neutrinoless double-β decay, following earlier work by Drukarev et al. [Phys. Rev. C 94, 035504 (2016)]. We find in agreement with Drukarev et al. that the Q value of the decay is reduced by any atomic excitation in the final-state but not by the difference in atomic binding energy of the initial and final state atoms, as was recently suggested by Mei et al. [Mod. Phys. Lett. A 37, 2250058 (2022)] and Mei et al. [Nucl. Phys. A. 1032, 122623 (2023)]. We discuss how the absorption of subsequent atomic de-excitation ejecta shifts the neutrinoless double-β decay peak back to its nominal value. We propose an in situ experimental verification of these ideas, and discuss the impact of shake-up and shake-off on two-neutrino double-β decay spectral shapes.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Charm contribution to ultrahigh-energy neutrinos from newborn magnetars

Newborn, strongly magnetized neutron stars (so-called magnetars) surrounded by their stellar or merger ejecta are expected to be sources of ultrahigh-energy neutrinos via decay of mesons produced in hadronic interactions of protons which are accelerated to ultrahigh energies by magnetic dissipation of the spindown energy. In this work, we show that not only pions and kaons but also charm hadrons, which are typically neglected due to their small production cross sections, can represent dominant contributions to neutrino fluence at ultrahigh energies, because of their short lifetimes, while the ultrahigh-energy neutrino fluence from pion and kaon production is suppressed at early times due to their significant cooling before their decay. We show that the next-generation detectors such as Probe Of Extreme Multi-Messenger Astrophysics (POEMMA), Giant Radio Array for Neutrino Detection (GRAND) and IceCube-Gen2 have a good chance of observing neutrinos, primarily originating from charm hadrons, from nearby magnetars. We also show that neutrinos from nearby magnetar-driven merger novae could be observed in the time interval between 10 2 s and 10 3 s, where the charm hadron contribution is dominant for neutrino energies above 10 8 GeV, of relevance to next generation detectors. We also comment on potential impacts of the charm hadron contribution to the diffuse neutrino flux.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

New method for detecting fast neutrino flavor conversions in core-collapse supernova models with two-moment neutrino transport

Fast-pairwise neutrino oscillations potentially affect many aspects of core-collapse supernovae (CCSNe): the explosion mechanism, neutrino signals, and nucleosynthesis in the ejecta. This particular mode of collective neutrino oscillations has a deep connection to the angular structure of neutrinos in momentum space; for instance, the appearance of electron-neutrino lepton number (ELN) angular crossings in momentum space is a good indicator of the occurrence of flavor conversions. However, many multidimensional CCSN simulations are carried out with approximate neutrino transport (such as two-moment methods), which limits access to the angular distributions of neutrinos, i.e., inhibits ELNcrossing searches. In this paper, we develop a new method of searching for ELN crossing in these CCSN simulations. The required data is the zeroth and first angular moments of neutrinos and the matter profile, all of which are available in CCSN models with the two-moment method. One of the novelties of our new method is the use of a ray-tracing neutrino transport to determine ELNs in the direction of the stellar center. It is designed to compensate for shortcomings of the crossing searches with only two angular moments. We assess the capability of the method by carrying out a detailed comparison with results of full Boltzmann neutrino transport in 1D and 2D CCSN models. In this work, we find that the ray-tracing neutrino transport improves the accuracy of crossing searches; indeed, the appearance/disappearance of the crossings is accurately detected even in the region of forward-peaked angular distributions. The new method is computationally cheap and has the benefit of efficient parallelization; hence, it will be useful for ELN-crossing searches in any CCSN models that employ two-moment neutrino transport.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Estimating outflow masses and velocities in merger simulations: Impact of r -process heating and neutrino cooling

Here, the determination of the mass, composition, and geometry of matter outflows in black hole-neutron star and neutron star-neutron star binaries is crucial to current efforts to model kilonovae and to understand the role of neutron star merger in r-process nucleosynthesis. In this manuscript, we review the simple criteria currently used in merger simulations to determine whether matter is unbound and what the asymptotic velocity of ejected material will be. We then show that properly accounting for both heating and cooling during r-process nucleosynthesis is important to accurately predict the mass and kinetic energy of the outflows. These processes are also likely to be crucial to predict the fall-back timescale of any bound ejecta. We derive a model for the asymptotic velocity of unbound matter and binding energy of bound matter that accounts for both of these effects and that can easily be implemented in merger simulations. We show, however, that the detailed velocity distribution and geometry of the outflows can currently only be captured by full three-dimensional fluid simulations of the outflows, as nonlocal effect ignored by the simple criteria used in merger simulations cannot be safely neglected when modeling these effects. Finally, we propose the introduction of simple source terms in the fluid equations to approximately account for heating/cooling from r-process nucleosynthesis in future seconds-long three-dimensional simulations of merger remnants, without the explicit inclusion of out-of-nuclear statistical equilibrium reactions in the simulations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Fast flavor instability in hypermassive neutron star disk outflows

Here, we examine the effect of neutrino flavor transformation by the fast flavor instability (FFI) on long-term mass ejection from accretion disks formed after neutron star mergers. Neutrino emission and absorption in the disk set the composition of the disk ejecta, which subsequently undergoes r-process nucleosynthesis upon expansion and cooling. Here we perform 28 time-dependent, axisymmetric, viscous-hydrodynamic simulations of accretion disks around hypermassive neutron stars (HMNSs) of variable lifetime, using a 3-species neutrino leakage scheme for emission and an annular-lightbulb scheme for absorption. We include neutrino flavor transformation due the FFI in a parametric way, by modifying the absorbed neutrino fluxes and temperatures, allowing for flavor mixing at various levels of flavor equilibration, and also in a way that aims to respect the lepton-number preserving symmetry of the neutrino self-interaction Hamiltonian. We find that for a promptly-formed black hole (BH), the FFI lowers the average electron fraction of the disk outflow due to a decrease in neutrino absorption, driven primarily by a drop in electron neutrino/antineutrino flux upon flavor mixing. For a long-lived HMNS, the disk emits more heavy lepton neutrinos and reabsorbs more electron neutrinos than for a BH, with a smaller drop in flux compensated by a higher neutrino temperature upon flavor mixing. The resulting outflow has a broader electron fraction distribution, a more proton-rich peak, and undergoes stronger radiative driving. Disks with intermediate HMNS lifetimes show results that fall in between these two limits. In most cases, the impact of the FFI on the outflow is moderate, with changes in mass ejection, average velocity, and average electron fraction of order ~10%, and changes in the lanthanide/actinide mass fraction of up to a factor ~2.

79 ASTRONOMY AND ASTROPHYSICS↗