Fast neutrino flavor conversion, ejecta properties, and nucleosynthesis in newly-formed hypermassive remnants of neutron-star mergers

Manu George,1, ∗ Meng-Ru Wu,1, 2, 3, † Irene Tamborra,4, ‡ Ricard Ardevol-Pulpillo,5, § and Hans-Thomas Janka5, ¶ 1Institute of Physics, Academia Sinica, Taipei, 11529, Taiwan 2Institute of Astronomy and Astrophysics, Academia Sinica, Taipei, 10617, Taiwan 3Physics Division, National Center for Theoretical Sciences, 30013 Hsinchu, Taiwan 4Niels Bohr International Academy and DARK, Niels Bohr Institute, University of Copenhagen, Blegdamsvej 17, 2100, Copenhagen, Denmark 5Max-Planck-Institut f¨urAstrophysik, Postfach 1317, 85741 Garching, Germany (Dated: November 16, 2020) Neutrinos emitted in the coalescence of two neutron stars affect the dynamics of the outflow ejecta and the nucleosynthesis of heavy elements. In this work, we analyze the neutrino emission properties and the conditions leading to the growth of flavor instabilities in merger remnants consisting of a hypermassive neutron star and an accretion disk during the first 10 ms after the merger. The analyses are based on hydrodynamical simulations that include a modeling of neutrino emission and absorption effects via the “improved leakage-equilibration-absorption scheme” (ILEAS). We also examine the nucleosynthesis of the heavy elements via the rapid neutron-capture process (r- process) inside the material ejected during this phase. The dominant emission ofν ¯e over νe from the merger remnant leads to favorable conditions for the occurrence of fast pairwise flavor conversions of neutrinos, independent of the chosen equation of state or the mass ratio of the binary. The nucleosynthesis outcome is very robust, ranging from the first to the third r-process peaks. In 5 particular, more than 10− M of strontium are produced in these early ejecta that may account for the GW170817 kilonova observation. We find that the amount of ejecta containing free neutrons after the r-process freeze-out, which may power early-time UV emission, is reduced by roughly a factor of 10 when compared to simulations that do not include weak interactions. Finally, the potential flavor equipartition between all neutrino flavors is mainly found to affect the nucleosynthesis outcome in the polar ejecta within . 30◦, by changing the amount of the produced iron-peak and first-peak nuclei, but it does not alter the lanthanide mass fraction therein.

I. INTRODUCTION antineutrinos can be produced by the merger as matter is heated up to several tens of MeV due to the collision of two NSs. Neutrinos play an important role in deter- Compact binary systems consisting of two neutron mining the cooling of the merger remnant, changing the stars (NS) or a NS and a black hole (BH) can lose their composition of the ejecta, and altering the r-process out- angular momentum through continuous emission of grav- come and the kilonova emission properties. itational waves (GW), eventually leading to the merging of the compact objects. Such merger events have long The early-time blue color of the GW170817 kilonova been considered to be the sites producing short gamma- and the recently inferred amount of strontium produc- ray bursts (sGRB) and synthesizing heavy elements tion [11] both suggest that the merger ejecta contain via the rapid neutron-capture process (r-process) [1–3], some less neutron-rich material with electron fraction which powers electromagnetic transients in optical and per nucleon Ye & 0.3. On the other hand, numerical infrared wavelengths, the so-called kilonovae [4–7]. simulations which include weak interactions of neutrinos The first detected GW emission from a binary neu- with nucleons all suggest that neutrino emission and ab- tron star merger event by the LIGO and Virgo Col- sorption have the effect to reduce the neutron-richness laborations (GW170817) together with multi-wavelength of the outflow launched at different post-merger phases, electromagnetic observations have confirmed theoretical preferably in the direction perpendicular to the merger arXiv:2009.04046v2 [astro-ph.HE] 13 Nov 2020 predictions [8–10]. Future observations, like the one of plane [12–20]. In particular, recent work suggests that GW170817, will be able to offer further opportunities to neutrino absorption can be responsible for increasing Ye precisely determine the population of binary NS systems even for the early-time dynamical ejecta in the polar di- and the yet-uncertain rich physics involved in binary NS rection [15, 16, 19, 21–23]. Meanwhile, the νeν¯e pair anni- + − mergers, including the nuclear equation of state (EoS) hilation to e e pairs above the accretion disk, although and the properties of neutron-rich nuclei. In order to it might not be the dominant driver, can also contribute achieve these goals, solid theoretical modeling of mergers to launch the sGRB jet [2, 24–27] is needed without any doubt. Neutrinos additionally undergo flavor conversions The interaction of neutrinos with matter and their fla- above the merger remnant, altering the neutrino absorp- vor conversions in the binary NS merger environment are tion rates on nucleons as well as their pair-annhilation among the most uncertain theoretical aspects that can af- rates [28–34], thus possibly affecting the interpretation fect the observables. A copious amount of neutrinos and of the observed signals. In particular, Ref. [33] showed 2 that favorable conditions for the so-called “fast flavor II. NEUTRINO EMISSION FROM BINARY conversion” [35–43] exist nearly everywhere above the NEUTRON STAR MERGERS merger remnant because of the disk geometry and the protonization of the merger remnant in its effort to reach A. Models of binary neutron star mergers a new beta-equilibrium state for the high temperatures produced in the merger process. Fast pairwise conver- We consider models of binary NS mergers sim- sions can give rise to rapid flavor oscillations of neutrinos ulated by a three-dimensional relativistic smoothed −1 within a length scale of (G nν nν ) (1) cm. ∼ F| e − ¯e | ≈ O particle hydrodynamics code that adopts confor- Subsequently, Ref. [34] adopted time-dependent neutrino mal flatness conditions [46]. The code has re- emission characteristics from simulations of merger rem- cently been coupled to an approximate neutrino nants consisting of a central BH with an accretion disk transport scheme called “improved-leakage-equilibration- and also found favorable conditions for fast flavor con- absorption scheme (ILEAS),” which is implemented version. By assuming full flavor equilibration between on a three-dimensional (3D) Cartesian coordinate grid all neutrino flavors, it was shown that the nucleosynthesis (x, y, z) with the axis perpendicular to the merging plane outcome in the neutrino-driven outflow from the BH–disk chosen as z axis. ILEAS serves as an efficient transport remnant can be largely altered [34]. method for− multidimensional simulations; when com- pared to two-moment neutrino transport results for pro- toneutron stars and post-merger tori, it captures well neutrino energy losses from the densest regions of the sys- In this work, we focus on the newly formed hypermas- tem as well as neutrino absorption in the free-streaming sive remnants of NS mergers. In fact, in the case of a regime [23, 44]. hypermassive NS, which might transiently exist in many Our fiducial models discussed in this section are from binary NS mergers, the surrounding torus-like equatorial Ref. [23] and simulate the mergers of two non-rotating bulge (“disk”) is exposed to strong neutrino irradiation NSs with mass 1.35 M each, with the EoS of DD2 [47, from the massive core of the merger remnant. This is 48] and SFHo [49]. In addition, we consider in later sec- different from the BH-torus configuration, where there is tions cases of unequal mass binaries consisting of two no such intense central neutrino source. Moreover, the NSs with 1.25 M and 1.45 M each, based on the same torus itself produces heavy-lepton neutrinos only with numerical scheme of Refs. [23, 44]. We note here that very low luminosities, whereas the hypermassive NS ra- we have taken all simulation data averaged over the az- diates high luminosities also of the heavy-lepton neutri- imuthal angle. This is a fairly good assumption as the nos. Both aspects make it necessary to investigate the merger remnant quickly reaches an approximately axi- question of fast flavor conversions for the case where the symmetric state after 2 3 ms. hypermassive NS with a surrounding disk/torus still ex- Figure1 shows the evolution∼ − of the energy luminos- ists. To this purpose, we focus on the first 10 ms after ity of different neutrino species and their average energy the merger of two NSs, during which a central' hyper- during the first 10 ms after the merger of two NSs for massive NS surrounded by an accretion disk forms. By our symmetric mergers 1. During this time, the central relying on recent hydrodynamical simulations with two object is a hypermassive NS supported by differential ro- 53 −1 different EoS and two different mass ratios for the bi- tation. The energy luminosity ofν ¯e (Lν 10 erg s ) ¯e ' nary, developed in Refs. [23, 44] and available at [45], we is a factor of 2 or 3 larger than the one of νe through- explore the neutrino emission properties, the conditions out the whole 10 ms, leading to continuous protonization for the occurrence of fast flavor instabilities, and the ef- of the remnant. The energy luminosities of the heavy- fects of neutrino absorption and flavor conversions on the lepton neutrino flavors, denoted by νx (with νx being neutron-richness of the ejecta as well as the nucleosyn- representative of one species of heavy-lepton neutrinos), thesis outcome. are 5 1052 erg s−1 per species. The averaged energies estimated∼ × by the leakage approximation show a clear hi-

erarchy of Eνe < Eν¯e < Eνx , reflecting the ordering of the temperatureh i h at decoupling;i h i neutrinos that decou- The paper is organized as follows. In Sec.II, we in- ple in the innermost region of the remnant have higher troduce the merger remnant models and the neutrino average energy (see below and Fig.2). When comparing emission properties. In Sec.III, we outline the frame- the energy luminosity evolution of the models with DD2 work used for the linear stability analysis and present and SFHo EoS, one sees that the latter produces higher our numerical results on the occurrence of fast flavor con- luminosities in all flavors. This is related to the fact that versions. We analyze the ejecta properties, the neutrino the SFHo EoS is softer than the DD2 EoS, which re- absorption effects on the evolution of Ye, and the impact sults in a NS with smaller radius. Consequently, a more of flavor equipartition in Sec.IV. We summarize our find- ings and discuss their implications in Sec.V. We adopt natural units with ~ = c = 1 for all equations throughout 1 More precisely, the time is measured with respect to the first the paper. minimum of the lapse function (see Ref. [23]). 3

of the ν ,ν ¯ , and ν emission surfaces. The emission 14.0 DD2 e SFHo e e e x

DD2 e SFHo e surface for a given neutrino species να (νe,ν ¯e, or νx) is

) DD2 x SFHo x defined by a surface where the energy-averaged optical s /

g 10.5 depth is τ = 2/3 (see Sec.IIB for details). One can r να e

2 see from the density profiles that the remnant consists of 5

0 a rotating deformed hypermassive NS surrounded by an 1 (

y 7.0

t accretion disk. The central object retains its initial low i s

o Ye . 0.1 (see the third rows) and it is surrounded by a n i dense and neutron-rich disk, which is opaque to neutri- m 3.5 u

L nos for an extended radius out to 40 60 km. Both νe ∼ − andν ¯e decouple at locations where the matter density is approximately between 1011 1012 g cm−3 and the tem- 0.0 − perature is 5 MeV. Theν ¯e emission surfaces generally 0.0 2.5 5.0 7.5 10.0 ∼ time [ms] sit inside the ones of νe during this period, independently of the adopted EoS.

25 The size of the neutrino emission surfaces slightly ex- pands as the remnant evolves, due to the settling of the post-merger object and the redistribution of matter with 20 high angular momentum toward the equatorial plane, DD2 e SFHo e where the disk formation process proceeds. Moreover, DD2 e SFHo e 15 as the remnant keeps protonizing, Ye inside the neutrino DD2 SFHo x x surfaces gradually increases with time. The thick disk of 10 the remnant in the model with the softer SFHo EoS pro- tonizes faster as discussed above, and thus has higher Ye inside the disk when compared to the profiles with DD2

Mean energy (leakage) [MeV] 5 EoS. Note that in the polar region close to the z-axis, a 0.0 2.5 5.0 7.5 10.0 high Ye & 0.4 funnel forms at later times, as both the dif- time [ms] ferences between the νe andν ¯e luminosities and average energies become smaller. FIG. 1. Time evolution of the neutrino energy luminosity In the bottom panels of both Fig.2 and Fig.3, we and average energy for different species, νe (in red),ν ¯e (in show the ratio of the difference of the number densities green) and ν (in blue) of the 1.35 + 1.35 M NS-NS merger x ofν ¯e and νe to the sum of theν ¯e and νe densities, i.e., simulation, with DD2 (solid lines) and SFHo (dashed lines) (nν nν )/(nν + nν ). Once again, as the remnant is EoS. ¯e − e ¯e ¯e protonizing and emitting moreν ¯e than νe, nearly any location above the νe surface in both models has this ra- tio larger than zero during the entire first 10 ms. The violent collision during the merging leads to higher tem- only exception is represented by the small patches in the peratures, and correspondingly higher neutrino emission polar region at 7.5 and 10 ms for the model with DD2 rates, when the SFHo EoS was adopted. This gives rise EoS. These patches are a consequence of the neutron- to a faster protonization of the remnant. ization that takes place locally around the poles of the Comparing Fig.1 to Fig. 2 of Ref. [34], where the latter high-density core of the merger remnant in the DD2 case. displays the neutrino emission properties for a BH rem- Because the stiff EoS prevents the merger core from fur- nant, one can see that the neutrino emission properties of ther contraction, the polar regions cool quickly, and the the electron-flavor neutrinos are comparable in the two density just inside the neutrinospheres increases by grav- scenarios. However, the non-electron-flavor neutrinos are itational settling, forcing the neutrinospheres to move in- more abundant in the models investigated in this work ward to smaller radii. This explains the more pronounced and have average energies higher than the electron-flavor polar trough of the neutrinospheres in the DD2 model neutrinos. In addition, the neutrino energy luminosities compared to the SFHo merger. Striving for a new beta- in the case with the BH accretion disk quickly decrease equilibrium state, now at lower temperature, the plasma after 20 ms (see Fig. 2 of Ref. [34]). Although we only begins to neutronize again, radiating more electron neu- ' focus on the first 10 ms after the coalescence, the neu- trinos than antineutrinos in both polar directions. This trino energy luminosity reaches a plateau in the models leads to the excess of νe relative toν ¯e outside the neu- where a hypermassive NS forms. trinospheres, visible as the two blue patches in the two The first three rows in Fig.2 and Fig.3 show the bottom right panels of Fig.2. Correspondingly, νe cap- baryon density, temperature, and electron fraction pro- tures dominateν ¯e captures in this region and the polar files in the half x z plane for x > 0 at 2.5, 5, 7.5, outflow becomes more and more proton-rich (red regions − and 10 ms post merger for the models with DD2 and of Ye > 0.5 around the z-axis in the right panels of the SFHO EoS, respectively. Also shown are the locations third row of Fig.2). The same trends are visible in the 4

FIG. 2. Hydrodynamic and neutrino properties from the 1.35 + 1.35 M NS-NS merger remnant simulation with DD2 EoS in the (x, z) plane for the time snapshot taken at 2.5, 5, 7.5, and 10 ms after the coalescence (from left to right columns). The first, second and third rows show the baryon mass density ρ, temperature T , the electron number fraction Ye. The fourth row displays the ratio (nν¯e − nνe )/(nν¯e + nνe ) . A value of this ratio above (below) 0 implies larger nν¯e (nνe ) locally. Also shown with solid lines are the contours of the νe (red),ν ¯e (green), and νx (blue) emission surfaces. 5

SFHo simulation (Fig.3), though less extreme and less across the whole emission surface; it remains quite stable rapidly evolving than in the DD2 run. throughout the simulated evolution until 10 ms after the plunge. The only exception is the snapshot at 7.5 ms for the model with DD2 EoS, which shows a dip in the B. Neutrino number densities on their emission ν¯e number density around x 10 km. This is because surfaces the enhanced deleptonization above' the poles of the high- density core (described in Sec.IIA) leads to a short, tran- The inner regions of the merger remnant are dense sient excess of νe overν ¯e near the neutrinospheres in the enough to trap neutrinos. This allows us to define a neu- northern hemisphere (see Fig.2). The effect is somewhat trino emission surface for each species να above which να pathological and unusual, also because it is considerably can approximately free-stream. As we will use the prop- less strongly developed in the southern hemisphere. No- erties of νe andν ¯e on their respective surfaces to con- tably, comparing Fig.4 with the bottom panel of Fig. 1 struct their angular distributions outside the νe surface of Ref. [34], one can see that the neutrino-antineutrino in Sec.III, we discuss below the time evolution and the asymmetry is more pronounced in the present models dependence on the adopted EoS of the neutrino densities than in the BH remnant case. on their emission surfaces. For any point x inside the simulation domain, the op- tical depth along a specific path γ to another point y III. FLAVOR INSTABILITY that a neutrino with an energy E traverses is given by Z y In Sec.IIIA we first briefly introduce the theoret- x −1 x0 ical formalism, viz., the dispersion relation (DR) ap- τνα (E, , γ) = λνα (E, )ds, (1) x proach [38], widely used in the literature to investi- gate the occurrence of neutrino flavor conversions. In where x0 is a point, ds is the differential segment along Sec.IIIB we look for the conditions leading to flavor in- γ, and λ (E, x0(s)) is the corresponding mean-free-path να stabilities using the simulation data introduced in Sec.II. at x0. For each species α, we determine the position We then apply the DR formalism to investigate the oc- of the neutrino decoupling surface in the same way as currence of flavor conversions above the neutrino emis- in Ref. [23], i.e. for every point x, the minimum of the sion surfaces in our merger remnant models in Sec.IIIC. spectral-averaged τνα (x, i) is computed along the six different directionsh (i ( x,i y, z)) through the edge ∈ ± ± ± of the simulation domain. The neutrino decoupling sur- A. Dispersion relation formalism face is then defined by the location corresponding to τν (x, i) = 2/3 in the six directions. The emission sur- h α i We adopt the density matrix formalism to describe faces computed in this way for νe,ν ¯e, and νx are the ones shown in Figs.2 and3. the statistical properties of the neutrino dense gas in- corporating flavor mixing [51]. For a given density ma- Figure4 shows the number densities of νe,ν ¯e, and trix %(p, x, t), its diagonal elements in the flavor basis, νx at their respective emission surfaces as functions of x %αα, record the phase-space distributions fν of a given for different snapshots. Note that the νe andν ¯e number α densities are computed by using Fermi-Dirac distribution neutrino flavor να at the space-time location (t, x) and functions with temperature and chemical potential ex- with momentum p. The off-diagonal terms %αβ carry tracted at each location on the neutrino surfaces, consis- the information about the neutrino mixing (neutrino fla- vor correlations). In the absence of any mixing, i.e., all tently with Ref. [23]. For νx, we rescale the number den- sity according to the procedure described in AppendixA neutrinos are in their flavor eigenstates, the off-diagonal to account for the trapping effect between the energy- elements vanish. surfaces and the emission surfaces, which is known to Neglecting general-relativistic effects and collisions of neutrinos with matter, the space-time evolution of the considerably reduce the νx luminosity when compared to the analogous values directly estimated through the local density matrix %(p, x, t) is governed by a Liouville equa- emission (see, e.g., [50]). The top panels show the neu- tion trino density corresponding to the DD2 EoS, while the ∂t%(p, x, t) + vp %(p, x, t) = i[Ω(p, x, t),%(p, x, t)], bottom panels display the same quantities for the SFHo · ∇ − (2) EoS. As discussed in Sec.IIA, the softer SFHo EoS leads where Ω is the Hamiltonian that accounts for the flavor to a generally higher temperature in the merger remnant; oscillations of neutrinos. On the left hand side of Eq. (2), as a consequence, theν ¯e luminosity (Fig.1) and number the first term takes care of the explicit time dependence density on the emission surface are higher than the one of % and the second term takes into account the neutrino obtained in the model with DD2 EoS. Both nν and nν e ¯e propagation with velocity vp p/ p for ultrarelativistic are larger in a region close to the pole where the temper- neutrinos. The Hamiltonian matrix' | | Ω on the right hand ature is higher (see Figs.2 and3). Of relevance to the side can be decomposed as occurrence of fast flavor conversions is the fact that the ν¯e emission is a factor of 2–3 larger than the one of νe Ω(p, x, t) = Ωvac + ΩMSW + Ωνν , (3) 6

FIG. 3. Same as Fig.2, but for the simulation with SFHo EoS. 7

DD2 2.5 ms DD2 5 ms DD2 7.5 ms DD2 10 ms 6 6 6 6 ) ν ν ν ν

3 e e e e − ν¯e ν¯e 4 ν¯e 4 ν¯e

cm 4 4 νx νx νx νx 33

(10 2 2 2 2 α ν n 0 0 0 25 50 0 20 40 0 0 25 50 0 0 25 50 x (km) x (km) x (km) x (km) SFHo 2.5 ms SFHo 5 ms SFHo 7.5 ms SFHo 10 ms 6 6 6 ) ν 6 ν ν ν

3 e e e e − 4 ν¯e ν¯e ν¯e 4 ν¯e cm 4 νx 4 νx νx νx 33

(10 2 2 2 α 2 ν n 0 0 0 0 20 40 0 0 25 50 0 25 50 0 25 50 x (km) x (km) x (km) x (km)

FIG. 4. Number density of νe (red solid curves),ν ¯e (green dash-dotted curves), and νx (blue dash lines) on their respective emission surfaces in the northern hemisphere for the 1.35 + 1.35 M merger models with DD2 EoS (top panels) and SFHo

EoS (bottom panels) at 2.5, 5, 7.5 and 10 ms post merger. Note that the νx number density has been rescaled, see text and AppendixA for details.

∗ where the first term Ωvac takes into account flavor con- %¯ = [(fν¯e + fν¯x ) + (fν¯e fν¯x )ξ ]/2 for antineutrinos. versions in vacuum. In a simplified two-flavor scenario, The entity− ξ is a matrix defined− as Ωvac = diag(ωv/2, ωv/2) in the mass basis with ωv =   2 2 − s S (m2 m1)/2E being the vacuum oscillation frequency of ξ = , (6) − S∗ s the neutrinos with energy E. − The second term on the right hand side of Eq. (3) em- where 1 s 1 and s 2 + S 2 =1. In the absence of bodies the effects of neutrino coherent forward scattering − ≤ ≤ | | | | with electrons and nucleons. In the flavor basis, this term any flavor correlation S = 0. Furthermore, as in previ- can be expressed as ous work that studied the fast neutrino flavor conversion, we omit the vacuum oscillation term in the following dis- cussion as ωv marginally affects the linear regime [52– ΩMSW = (√2GFne)diag(1, 0), (4) 54]. With these assumptions and introducing the met- µν where ne is the net electron number density. The last ric tensor η = diag(1, 1, 1, 1) and for any µ − − ν − term in Eq. (3) is the effective Hamiltonian due to the ν– contra-variant vector A , Aµ = ηµν A , we can recast the ν interaction. For a neutrino traveling with momentum Hamiltonian defined in Eq. (3) into the following form p,Ωνν is given by Z µ σ3 0 µ 0 0 0 Ω = v λµ + dΓ v vµξ(v )g(v ), (7) Z d3q 2 Ωνν = √2GF (1 vp · vq)(%(q, x, t) %¯(q, x, t)), (2π3) − − where v = (sin θ cos φ, sin θ sin φ, cos θ) with velocity (5) µ µ v = (1, v), dΓ = sin θdθdφ and λ = √2G ne(1, v ) where% ¯ is the corresponding density matrix for antineu- F m with vm being the vector of the fluid velocity of the back- trinos. The presence of (1 vp · vq) in Eq. (5) leads − ground matter. Since the rate of pairwise conversion is to multi-angle effects, i.e., neutrinos propagating in dif- much faster than any other inverse time scale involved in ferent directions experience different Ωνν . The equation the problem, we treat the background matter as station- of motion for antineutrinos can be obtained in a similar µ ary and homogeneous: λµv = λ . fashion, by replacing ω by ω in Ω . 0 v − v vac The quantity g(v) is related to the angular distribution We focus on a simplified system that deals with two of the neutrino ELN angular distribution neutrino flavors. Under this assumption, both the den- sity matrix % and the Hamiltonian Ω are 2 2 Hermi- √ g(v) = 2GF(Φνe Φν¯e ), (8) tian matrices and hence can be expanded in× terms of − the identity matrix and three Pauli matrices. Thus, we where Φνα = dnνα /dΓ. To study the growth of S in write % = [(fν + fν ) + (fν fν )ξ]/2 for neutrinos and the linear regime, we treat the flavor correlation S as a e x e − x 8

2 perturbation and neglect all terms of O(S ) or higher. where θn is the angle with respect to the normal direc- Taking S(t, x) = Q(ω, k)e−i(ωt−k·x), the EoM becomes tion of the location on the emission surface. Ignoring the minor effect of GR bending, we can then ray-trace Z µ 0 µ 0 0 the neutrino intensities from the emission surfaces to ob- vµλ¯ Q(ω, k) = dΓ v v g(v )Q(ω, k). (9) − µ tain their angular distributions at any location above the surfaces. In the above Eq. (9) we have introduced the four vec- µ R tor λ¯ = (ω λ0 0, k ), 0 dΓg(v) and Figures5 and6 show the obtained ELN distributions  R dΓvg(v).− Inspecting− − Eq. (9), one≡ can make the as a function of the local angular variables θ (angle with ≡ µ µ µ ansatz Q(ω, k) = vµa /vµλ¯ , with a being the coeffi- respect to the z-axis) and φ (angle with respect to the cients of eigenfunction solutions. Thus, Eq. (9) becomes x-axis on the x-y plane) at selected locations above the νe surface at 2.5, 5, 7.5 and 10 ms after the coalescence, for the 1.35 + 1.35 M merger models with DD2 and µν vµΠ aν = 0, (10) SFHo EoS, respectively. Here we only show the sign of the ELN distribution to highlight the crossing. The blue where we have used the definition shade represents the region where the net ELN is positive while the red region corresponds to negative ELN. The Z vµvν Πµν ηµν + dΓg(v) . (11) top panel shows the ELN distribution at a representative ¯σ ≡ vσλ point on the z axis and the middle and bottom panels show the ELN distributions at near-center and outer re- The EoM defined in Eq. (10) holds for any vµ. Thus, gions above the ν surface. Note here that Figs.5 and we have the condition Πµν a = 0. Eigenfunctions of the e ν 6 only show the ELN distribution for φ/π 0 as the latter have non-trivial solution only if Πµν satisfies the distributions for axi-symmetric emission surfaces≥ possess condition reflection symmetry, Φν ,ν (cos θ, φ) = Φν ,ν (cos θ, φ). e ¯e e ¯e − det[Πµν (ω, k)] = 0. (12) The angular coverage of the νe andν ¯e fluxes at a point (x, z) is determined by the geometry of the νe andν ¯e Equation (12) is the DR in flavor space. The solutions emission surfaces. On the other hand, the ELN angu- of the DR have been classified into several types [39, 41]. lar distribution is determined by the combination of the k If ω is real for real values of , a perturbation in S only emission surface geometries and their respective emission propagates without growing or damping, i.e., it stays in properties, including the relative strength of νe andν ¯e the linear regime, meaning that no significant flavor con- and the angular dependence of the emission. For exam- version occurs. On the other hand, an imaginary solu- ple, the ELN angular distribution at any point on the z tion of ω with Im(ω) > 0 corresponds to exponentially axis above the emission surface is independent of φ, re- growing modes. In other words, the flavor correlation flecting the (assumed) rotational symmetry of the merger S grows exponentially with time, leading to significant | | remnants about the z axis. As we move away from the z flavor conversion. Rigorous studies have been carried axis, the ELN distribution becomes dependent on φ (sec- out to understand the characteristics of the above DR ond and third panels). For both EoS, the angular cov- with respect to the ELN angular distribution of neutri- erage in θ slightly increases with time for a given (x, z), nos [38, 39, 41]. It was shown that in the presence of a caused by the expansion of the emission surfaces. This crossing in the ELN distribution, the DR will yield com- is different from what was found in Ref. [34] where the plex ω solutions for real k, leading to temporal instabili- central remnant is a BH for which the size of the νe/ν¯e ties. In the following, we examine the ELN distributions emitting torus surfaces shrinks with time. Asν ¯e are more above the merger remnants and the corresponding flavor abundant than νe in most parts of the emission region instabilities. during the first 10 ms (see Fig.4), the overall shapes of the ELN crossings are qualitatively similar. The only ex- ception is represented by the snapshot at 7.5 ms for the B. Neutrino electron lepton number angular case with DD2 EoS, for which the inner region above the distribution merger remnant shows a double ELN crossing structure, see the left and middle panels in the third row. This is To construct the ELN distribution above the merger related to the dip of theν ¯e density on theν ¯e surface at remnant, we follow a method similar to the one adopted 10 km discussed in Sec.IIB. Also, the ELN angular ∼ in Ref. [34]. First, we assume that both νe andν ¯e freely distribution in Fig.6 is more extended toward θ = 0 com- propagate outside their respective emission surfaces de- pared to Fig.5, resulting from slightly larger radii of the fined in Sec.IIB. Second, we approximate their forward- neutrino emission surfaces in the SFHo model. For the peaked angular distributions at each point on the emis- simulations adopting 1.25 + 1.45 M binaries, we have sion surfaces as similarly checked the ELN distributions during the same time snapshots for both EoS. Unsurprisingly, ELN cross- nνe,ν¯e Φν ,ν (θn) = (1 + cos θn), (13) e ¯e 4π ings appear at all times as those shown here. 9

FIG. 5. Neutrino electron lepton number angular crossings at different positions (x, z) = (0, 25), (10, 25), (40, 25) km and times t = 2.5, 5, 7.5, 10 ms, above the remnant for the 1.35 + 1.35 M model with DD2 EoS. The red (blue) shaded region represents the region where Φνe − Φν¯e > 0 (Φνe − Φν¯e < 0).

FIG. 6. Same as Fig.5, but for the 1 .35 + 1.35 M simulation with SFHo EoS.

C. Flavor instabilities for fast pairwise conversion spect whether solutions containing non-zero Im(ω) for a

given k can be found2. As the ELN distributions above

After obtaining the ELN distributions above the neu- trino emitting surface, we numerically solve the DR 2 Note that the positive and negative Im(ω) solutions always ap- [Eq. (12)] starting from the outer neutrino surface to in- pear together. 10

DD2 Symmetry preserving DD2 Symmetry breaking

0.6 ) 1.0 2.5 ms 2.5 ms 1 5 ms 5 ms m

c 7.5 ms 0.4 7.5 ms ( | ) 0.5 10 ms 10 ms

( 0.2 m I | 0.0 0.0 30 20 10 0 10 5.0 2.5 0.0 2.5 5.0 7.5 10.0 12.5 SFHo Symmetry preserving SFHo Symmetry breaking 2

) 2.5 ms 0.75 2.5 ms 1 5 ms 5 ms m

c 7.5 ms 0.50 7.5 ms (

| 1 ) 10 ms 10 ms

( 0.25 m I | 0 0.00 50 40 30 20 10 0 10 20 5 0 5 10 15 20 k(cm 1) k(cm 1)

FIG. 7. Solutions of |Im(ω)| obtained from the DR [Eq. (12)] for k = (0, 0, k) at the location (x, z) = (10, 25) km above the merger remnant for different times. The top (bottom) panels show the results for the 1.35 + 1.35 M model with the DD2 (SFHo) EoS, while the left (right) panels show the symmetry-preserving (symmetry-breaking) solutions.

DD2 Symmetry preserving DD2 Symmetry breaking 2

) (0, 25) 0.75 (0, 25) ) 1 (10, 25) (10, 25)

m 0.50

c (40, 25) (40, 25)

( 1 | )

( 0.25 m I | 0 0.00 80 60 40 20 0 20 5 0 5 10 15 20 25 30 SFHo Symmetry preserving SFHo Symmetry breaking

1.0 ) (0, 25) (0, 25) 1 2 (10, 25) (10, 25) m

c (40, 25) (40, 25) ( |

) 0.5 1 ( m I | 0 0.0 80 60 40 20 0 5 0 5 10 15 20 25 k(cm 1) k(cm 1)

FIG. 8. Solutions of |Im(ω)| obtained from the DR [Eq. (12)] for k = (0, 0, k) at 5 ms for different locations labeled by the coordinates (x, z) expressed in km. The top (bottom) panels show the results for the 1.35 + 1.35 M model with the DD2 (SFHo) EoS, while the left (right) panels show the symmetry-preserving (symmetry-breaking) solutions. 11 the merger remnants with axial-symmetry have a reflec- parts of the ELN distribution (see e.g., Ref. [41]). Look- tion symmetry with respect to φ φ, one can obtain ing at the ELN crossing pattern shown in Figs.5 and6, two different solutions that correspond→ − to the reflection the locations above the middle part of the remnant have symmetry-preserving and symmetry-breaking cases [33]. large angular area of g(v) > 0 due to the larger separa- We show the obtained Im(ω) as a function of kz, tion of the νe andν ¯e surfaces (see Fig.2 and3). This, taking k = k = 0, at different| | times for a location in turn, enhances the corresponding values of Im(ω)/µ x y | | at (x, z) = (10, 25) km above the emitting surfaces in relative to the ones in the inner region closer to the z-axis. Fig.7 and the solutions for different locations corre- Likewise, when comparing the values of Im(ω)/µ to | | sponding to the ELN crossing shown in Fig.8 at 5 ms. the ones shown in Ref. [34], the more dominantν ¯e emis- The top (bottom) panels are for cases with DD2 (SFHo) sion relative to νe here, together with the smaller sepa- EoS while the left (right) panels show the symmetry- ration of their emission surfaces, results in smaller values preserving (symmetry-breaking) solutions. of Im(ω)/µ . | | Flavor instabilities with growth rate of (1) cm−1 ex- On the other hand, while Ref. [34] showed that the re- O ist at all locations at all times for a large range of kz. gion where the flavor instability exists shrinks on a time scale of (10) ms as the BH-disk remnant evolves, The symmetry preserving solution of the DR has two ∼ O branches while the symmetry breaking solution has only the instability region found here remains stable within one branch, similar to what was found in Ref. [33]. Com- the examined 10 ms of post-merger evolution. This can paring the results for the models with DD2 and SHFo have important consequences for the growth of the insta- EoS, the growth rate of the flavor instability is generally bilities and seems to favor the eventual development of larger in the latter, as the neutrino emission is stronger flavor conversions in the non-linear regime. The effect of with the SHFo EoS. The shape of the solutions is rather advection hindering the growth of the flavor instabilities stable over time for the case with SFHo EoS. On the for systems with non-sustained and fluctuating unstable other hand, the model with DD2 EoS shows a somewhat conditions shown by Ref. [56] may therefore not happen stronger time-dependence. In particular, the range of above the merger remnant, as shown in Ref. [55]. In the kz that leads to non-zero Im(ω) as well as the value of case of the models with unequal mass binaries, due to Im(ω) both decrease at the| snapshot| of 7.5 ms, when the the similar ELN crossing features, we find unstable re- double| | crossing shape of the ELN appears (see Fig.5). gions above the merger remnant disk and growth rates By comparing the solutions at different locations for very similar to the symmetric models (results not shown the t = 5 ms snapshot, Fig.8 shows that Im(ω) is larger here). closer to the z-axis for both the symmetry-preserving| | and We here focus on the diagnostics of flavor instabilities, symmetry-breaking solutions. This, again, is caused by but do not distinguish whether they belong to the con- vective or absolute type (see Refs. [41, 57]). This would the fact that Φνe and Φν¯e are largest in this region (see Fig.4). require solving the full dispersion relation for the com- plete (ω, k) space. As shown in Figs.7 and8, the symmetry-breaking so- lutions at all times and all locations contain non-zero Im(ω) for the mode k = 0. We further show in Fig.9 IV. MERGER EJECTA AND |the contour| plot of Im(ω) above the neutrino surface | | NUCLEOSYNTHESIS OF THE HEAVY for this mode. The plots in this figure clearly show that ELEMENTS growth rates of the instability of (1) cm−1 exist ev- erywhere above the remnant at∼ all O times. The region closer to the z-axis just above the emission surface has In this section, we first analyze the properties of the material ejected during the first 10 ms post merger, ex- the maximal growth rates. Moving away from the sur- ' face, the neutrino flux is suppressed by the geometric amine how neutrino absorption affects the evolution of Ye effects and the magnitude of the instability decreases in of the outflow material, and discuss the nucleosynthesis all cases. For the sake of comparison with the existing outcome in the absence of neutrino flavor conversion in literature, we find Im(ω)/µ (10−2) for all cases ex- Sec.IVA. In Sec.IVB, we further explore the potential amined here; this value| of the| ∼ growth O rate is similar al- effect of fast flavor conversion on the neutrino absorp- beit a bit smaller than the ones reported in Refs. [34, 55]. tion rates and the outcome of nucleosynthesis of in these The growth rates Im(ω)/µ are generally larger toward ejecta. | | the middle region above the νe surface, in agreement with the findings of Ref. [33, 34]. This is due to the relative strength of the positive vs. negative ELN strength, in the A. Ejecta properties and nucleosynthesis proximity of crossings between the νe andν ¯e angular dis- tributions [56]. For a system that is more dominated by The dynamical ejecta masses extracted at 10 ms post- νe orν ¯e, i.e., where the positive ELN distribution dom- merger for the 1.35 + 1.35 M merger simulations are −3 −3 inates the negative parts or the other way around, one 2.0 10 M and 3.3 10 M with the DD2 and expects that the value of Im(ω)/µ is smaller than in a SFHo∼ × EoS, represented by× 783 and 1263 tracer particles, | | more balanced system with similar positive and negative respectively. For the 1.25+1.45 M merger cases, the dy- 12

FIG. 9. Contour plots showing |Im(ω)| above the νe surface for k = 0 at t = 2.5, 5.0, 7.5 and 10.0 ms, from left to right, for the 1.35 + 1.35 M model. The top (bottom) panels are based on the model with the DD2 (SFHo) EoS. The red and green solid lines represent the locations ofν ¯e and νe emission surfaces, respectively.

0.4 0.2 0.2 0.16 (a) DD2 (b) (c) (d) 0.35 SFHo 0.14 0.3 0.15 0.15 0.12 0.25 0.1 ej,tot 0.2 0.1 0.1 0.08 /M

ej 0.15 0.06

M 0.1 0.05 0.05 0.04 0.05 0.02 0 0 0 0 0123456789 30 60 90 120 150 0 0.2 0.4 0.6 0.8 0 0.1 0.2 0.3 0.4 0.5 0.6 (no osc) t(r = 100km)(ms) θej(degree) (vr /c)(r = 100km) Ye (r = 100km) (no osc) Y 10 10 e 0.6 (e) DD2 0.8 (f) DD2 (g) SFHo 0.8 (h) SFHo 8 0.7 8 0.7 0.5 0.6 0.6 0.4 6 0.5 6 0.5 = 100km) = 100km) 0.3 r r 0.4 0.4 )( 4 )( 4 = 100km)(ms) = 100km)(ms) /c /c 0.3 0.3 0.2 r r r r ( ( v v t t ( ( 0.2 0.2 2 2 0.1 0.1 0.1 0 0 0 0 0 0 30 60 90 120 150 180 0 30 60 90 120 150 180 0 30 60 90 120 150 180 0 30 60 90 120 150 180 θej (degree) θej (degree) θej (degree) θej (degree)

FIG. 10. Histograms of the ejecta mass fraction when the ejecta reach r = 100 km at t(r = 100 km), [panel (a)], θej [panel (b)], radial velocity vr at r = 100 km, and Ye at r = 100 km without flavor conversions for the simulations with 1.35 + 1.35 M NS binaries. Panel (e) [(g)] and (f) [(h)] show the distributions of each tracer particle in the t(r = 100) km–θej and (vr/c)–θej planes with the DD2 (SFHo) EoS. The color coding of each tracer particle indicates its Ye value at r = 100 km without considering flavor conversions.

−3 0 0 namical ejecta masses are 3.2 10 M (DD2) and rates of νe andν ¯e on neutrons and protons, λνe and λν¯e −3 ∼ × 8.7 10 M (SFHo), represented by 1290 and 4398 (the superscript 0 here denotes the case where any neu- × tracer particles. In computing the evolution of Ye for a trino flavor conversion is omitted), as detailed in Appen- given tracer particle, we first post-process the neutrino dicesA andB. These rates are then combined with the emission data from Ref. [23] to calculate the absorption nuclear reaction network used in Refs. [58, 59] to com- 13 pute the nucleosynthesis yields in the merger ejecta. For in the SFHo model have more high-Ye ejecta than in the each tracer particle, we begin the network calculation ei- DD2 model [see panel (f)], different from the main bulk ther at a location where T = 50 GK or just outside the of the ejecta. This is because these ejecta are mainly disk with height of 25 km and radius of 55 km to make driven during the early post-merger phase. During this sure that it is outside the νe emission surface. The initial early phase, positron capture is responsible for raising Ye nuclear abundances are calculated by using the nuclear in the ejecta, rather than neutrino absorption. Thus, in statistical equilibrium (NSE) condition. For T > 10 GK, the model with SFHo EoS, the higher post-merger tem- we only compute the weak reaction rates to track the Ye perature of the remnant due to the more violent merger evolution while assuming NSE at each moment to obtain leads to relatively higher Ye material in the early fast the abundances. When T 10 GK, we instead follow ejecta. ≤ the full evolution of all nuclear species and include the We show in Fig. 11 the same quantities as in Fig. 10, feedback due to the nuclear energy release on the ejecta but for the models with 1.25 + 1.45 M mass binaries. temperature following Ref. [60]. Qualitatively, the features are similar to the ones de- In Fig. 10, we show the distributions of the ejecta scribed above, but they are quantitatively different. For mass for the 1.35 + 1.35 M models as a function of instance, the division between different episodes of mass the time when the ejecta reach the radius r = 100 km at ejection is less clear, in particular for the SHFo EoS model t(r = 100 km) in panel (a), the angle θej (relative to the [see panel (a)]. On average, most of the ejecta have z-axis) when they leave the simulation domain in panel higher velocities, peaked at 0.25c [see panel (c), (f) ∼ (b), and the radial velocity vr at r = 100 km in panel and (h)]. Moreover, a larger fraction of the ejecta has (c). We also show, in panel (d), the distributions of Ye for lower Ye . 0.2, despite the fact that similarly wide dis- the ejecta at r = 100 km, in the absence of flavor mixing. tributions of 0.01 . Ye . 0.5 are obtained. Additionally, panels (e) and (g) [(f) and (h)] show the Figure 12 shows the resulting abundance distributions distributions of each tracer particle in the plane spanned for the ejecta of the symmetric and asymmetric merger by t(r = 100 km) [vr(r = 100) km] and θej with the cor- models with the DD2 and SFHo EoS, in comparison to responding Ye at 100 km labeled in color for the DD2 and the re-scaled solar r abundance pattern [61]. Due to the SHFo EoS, respectively. These plots show that, indepen- wide-spread Ye distribution ranging from 0.02 0.5 in dently of the adopted EoS, most of the material is ejected all models, all three r-process peaks at A 90, 130,− and ◦ ' within the first 2 ms with lower Ye 0.3, within 45 195 are in relatively good agreement with the solar abun- ∼ . ∼ away from the mid-plane. The ejecta launched in the dance pattern. The difference of the high Ye component second episode [t(r = 100 km) 3 4 ms] are simi- discussed above only results in abundance variations at ' − larly distributed closer to the equatorial plane and have A 50 80. In particular, we compute the amount of ' − Ye 0.4. After that, neutrino irradiation continues to strontium present at the time of a day, which is of rele- ' power the mass ejection along the polar direction with vance to the potential identification in the spectral anal- Ye reaching 0.5. In terms of the ejecta kinematics, be- ysis of the GW170817 kilonova [11]. The total amount ∼ −5 −5 cause of the weak interactions, it is clear that Ye can be of strontium is 8.9 10 M and 3.9 10 M raised up to 0.3 0.4 mostly for the material expanding for the DD2 and' SFHo× EoS models with' equal× mass bi- − more slowly with vr/c 0.1 0.2 at r = 100 km, even naries. For the asymmetric merger models, the corre- ∼ − −4 around the equatorial plane, and reach 0.5 closer to sponding amount of strontium is 2.66 10 M and ∼ −4 ' × the polar direction. 3.25 10 M , respectively. These results are consistent × The main differences between the DD2 and the SFHo with the findings of Ref. [11]. EoS results are the following. First, for the simulation References [22, 60, 62, 63] found that some merger with SHFo EoS, the amount of ejecta is larger, particu- ejecta have low Ye and fast expansion time scale to allow larly during the first ejection episode [see panel (a)], than for a neutron-rich freeze-out during the r-process nucle- that with the DD2 EoS; this is a direct consequence of osynthesis. A potentially thin layer of “neutron-skin” at the more violent collision of the two NSs in the case with the outskirts of the ejecta may possibly power an early- softer EoS (SHFo), as discussed in Ref. [46]. Second, the time UV emission at hours post-merger due to the ∼ amount of high-velocity ejecta with vr/c & 0.5 for the radioactive heating of neutron decay [7, 62]. We find SFHo case is also higher for the same reason [see pan- that the amount of free neutrons at the end of the r- −6 els (c) and (f)]. Third, the Ye distribution [panel (d)] process, for equal-mass binaries, is 6.2 10 M and −6 ' × for both models shows that Ye(r = 100 km) is slightly 6.6 10 M with the DD2 and SFHo EoS, respectively. larger for the DD2 EoS model than for the SFHo EoS These× numbers are roughly a factor of 10 smaller than one. This reflects the higherν ¯e luminosity obtained in what was found in Ref. [60], which analyzed simulation the SHFo EoS model (see Fig.1) leading to larger ¯νe ab- trajectories without including the weak interactions. The sorption rates on protons; hence, Ye is raised less for the reduction of the amount of free neutrons at the end of bulk of the ejecta. The high Ye tail extends to values the r-process is related to the high post-merger temper- & 0.5 in the DD2 EoS model, while it remains . 0.5 in ature effect raising Ye even for the early fast ejecta. For the SFHo EoS model [see panels (c) and (e)]. Interest- the unequal mass binaries, the corresponding amount of −6 −6 ingly, it is noticeable that the high velocity components free neutrons is 9.2 10 M and 3.1 10 M for ' × × 14

0.4 0.2 0.2 0.16 (a) DD2 (b) (c) (d) 0.35 SFHo 0.14 0.3 0.15 0.15 0.12 0.25 0.1 ej,tot 0.2 0.1 0.1 0.08 /M

ej 0.15 0.06

M 0.1 0.05 0.05 0.04 0.05 0.02 0 0 0 0 0123456789 30 60 90 120 150 0 0.2 0.4 0.6 0.8 0 0.1 0.2 0.3 0.4 0.5 0.6 (no osc) t(r = 100km)(ms) θej(degree) (vr /c)(r = 100km) Ye (r = 100km) (no osc) Y 10 10 e 0.6 (e) DD2 0.8 (f) DD2 (g) SFHo 0.8 (h) SFHo 8 0.7 8 0.7 0.5 0.6 0.6 0.4 6 0.5 6 0.5 = 100km) = 100km) 0.3 r r 0.4 0.4 )( 4 )( 4 = 100km)(ms) = 100km)(ms) /c /c 0.3 0.3 0.2 r r r r ( ( v v t t ( ( 0.2 0.2 2 2 0.1 0.1 0.1 0 0 0 0 0 0 30 60 90 120 150 180 0 30 60 90 120 150 180 0 30 60 90 120 150 180 0 30 60 90 120 150 180 θej (degree) θej (degree) θej (degree) θej (degree)

FIG. 11. Same as Fig. 10, but for the 1.25 + 1.45 M models.

sumption of flavor equilibration is an extreme ansatz, es- -2 DD2, 1.35 + 1.35M pecially in the light of the findings of Ref. [55], which, SFHo, 1.35 + 1.35M⊙ however, relied on a simplified model of a relic merger DD2, 1.25 + 1.45M⊙ -3 SFHo, 1.25 + 1.45M⊙ disk. Moreover, this simple assumption does not pre-

)] ⊙ serve the total electron neutrino lepton number, which is A

( -4 a strictly conserved quantity in the case of pairwise fla- Y

[ vor conversions. Nevertheless, our extreme assumption

10 -5 for the flavor ratio is useful to explore the largest pos- sible impact that flavor conversions might have on the log -6 nucleosynthesis of the heavy elements. The correspond- ing neutrino absorption rates can be approximated by -7 osc 1 0 2 λ = λ + λν , (14) 40 60 80 100 120 140 160 180 200 220 νe 3 νe 3 x mass number, A osc 1 0 2 λ = λ + λν , (15) ν¯e 3 ν¯e 3 ¯x FIG. 12. Nucleosynthesis yields, Y (A), as a function of the nuclear mass number A at the time of 1 Gyr for all NS merger where λνx and λν¯x are the neutrino absorption rates on models investigated in this work. free nucleons assuming that all νx andν ¯x are converted to νe andν ¯e, as detailed in AppendicesA andB. We then perform the same nucleosynthesis calculations as in the DD2 EoS model and for the SFHo model. A future Sec.IVA for all tracer particles in all the merger models 0 0 osc osc (non)identification of this component may shed light on by replacing λνe and λν¯e by λνe and λν¯e . Below, we only the role of weak interactions in the post-merger environ- focus on the findings for the models with equal mass and ments. different EoS because the results obtained in the unequal mass binaries are qualitatively the same, independent of the EoS. B. Impact of flavor equipartition on Ye and We show in Fig. 13 the comparison of the ratio of nucleosynthesis λνe /λν¯e evaluated at r = 100 km for all tracer particles for the cases with and without flavor equipartition. The

Following previous work [34, 64], we assume that corresponding values of λν¯e are also shown. These figures fast flavor conversions lead to conditions close to flavor highlight that the neutrino absorption rates are orders of equipartition for neutrinos and antineutrinos. The as- magnitudes larger in the polar region than close to the 15

0 1 osc 1 (a) EoS: DD2 log10[λν¯ /(s− )] (b) EoS: DD2 log10[λν¯ /(s− )] e 4 e 4 1.4 1.4 3.5 3.5 1.2 1.2 3 3 e 1 e 1 ¯ ¯ 0 ν 2.5 osc ν 2.5 /λ 0.8 /λ 0.8

e 2 2 0 ν e osc ν λ 0.6 0.6 1.5 λ 1.5 0.4 1 0.4 1 0.2 0.5 0.2 0.5 0 0 0 0 0 20 40 60 80 100 120 140 160 180 0 20 40 60 80 100 120 140 160 180 θej (degree) θej (degree) 0 1 osc 1 (c) EoS: SFHo log10[λν¯ /(s− )] (d) EoS: SFHo log10[λν¯ /(s− )] e 4 e 4 1.4 1.4 3.5 3.5 1.2 1.2 3 3 e 1 e 1 ¯ ¯ 0 osc ν 2.5 ν 2.5 /λ 0.8 /λ 0.8

e 2 2 0 ν e osc ν λ 0.6 0.6 1.5 λ 1.5 0.4 1 0.4 1 0.2 0.5 0.2 0.5 0 0 0 0 0 20 40 60 80 100 120 140 160 180 0 20 40 60 80 100 120 140 160 180 θej (degree) θej (degree)

FIG. 13. Distribution of the ratios of the νe absorption rates on neutrons, λνe , to theν ¯e absorption rates on protons, λν¯e , as a function of θej for all tracer particles at r = 100 km in the 1.35 + 1.35 M models. Panels (a) and (c) are without flavor conversions, while panels (b) and (d) include flavor equipartition. The color of each point indicates the absolute rate of λν¯e . equator. For the case without neutrino flavor equipar- tion rates (see Fig. 13). 0 0 tition, nearly all tracer particles have λνe /λν¯e . 1, re- Figures 14(b)-(d) and 15(b)-(d) further show the com- flecting the strongerν ¯e flux emitted from its surface. parison of the Ye distribution for the cases with and with- With flavor equipartition, the nearly equal contribution out flavor conversion, for the ejecta classified into three ◦ of the converted νx andν ¯x significantly changes the ra- groups according to θej: the polar ejecta with 60 < θej osc osc ◦ ◦ ◦ ◦ | −◦ tio λνe /λν¯e for both EoS. For the case with DD2 EoS, 90 90 , the middle ejecta with 30 < θej 90 60 , ◦ ◦ | ≤ ◦ | − ◦ | ≤ ◦ most of the trajectories with θej . 60 or θej & 120 have and the equatorial ejecta with 0 θej 90 30 . osc osc osc osc ≤ | − | ≤ λνe /λν¯e & 1. On the other hand, the values of λνe /λν¯e Flavor equipartition influences the Ye distribution of the scatter around 1 for the model with SFHo EoS. polar ejecta by shifting the peak from 0.4 to 0.55 (0.5) for the DD2 (SFHo) model. The∼ larger (smaller)∼ The large change in λνe /λν¯e strongly affects the Ye distribution of the ejecta. Figures 14(a) and 15(a) show shift of the Ye peak in the DD2 (SFHo) model is re- osc osc (osc) (no osc) lated to the larger (smaller) values of λν /λν¯ shown in ∆Ye Ye Ye calculated at r = 100 km for e e ≡ − Fig. 13. For the middle ejecta, a fraction of ejecta orig- each tracer particles as a function of the corresponding (no osc) (osc) (osc) (no osc) inally with 0.3 Y 0.4 has Y 0.4 0.5 θ , where Y and Y are the Y values with and . e . e ej e e e when flavor equipartition is reached, while the' distribu-− without flavor equipartition. These results show that the tion with Y 0.3 is barely altered. As for the equatorial flavor equipartition can significantly increase Y of the e . e ejecta, the corresponding Y distribution is only affected ejecta up to 0.15 for θ < 60◦ or θ > 120◦ closer e ∼ ej ej negligibly, as expected. to the polar directions. In particular, the increase of Ye due to flavor equipartition for these ejecta is more pro- The impact of flavor equipartition in the neutrino- (no osc) driven ejecta studied in Ref. [34] is to lower Ye (see Fig. 11 nounced with Ye 0.3 0.4. For the ejecta with (no osc) (no osc)' − therein), while Ye increases in the models investigated in Ye . 0.2 or Ye 0.5, Ye is less affected. This this work, as discussed above. The main difference is '(no osc) is because the ejecta with Ye . 0.2 expand too fast that flavor equipartition leads to a larger reduction of for neutrino absorption to raise Ye either with or without λνe and λν¯e in the BH–torus case, due to the vanishingly flavor conversion. On the other hand, for the ejecta with small νx fluxes. Moreover, a significant part of the ejecta (no osc) Ye 0.5, λνe /λν¯e 1 even without flavor conver- in the BH–torus case is ejected on timescales of several sions. As' for the tracer' particles with 60◦ θ 120◦ tens of milliseconds during which the neutrino luminosi- ≤ ej ≤ closer to the disk mid-plane, Ye is barely influenced by ties decrease substantially (see Figs. 2 and 8 in Ref. [34]); flavor equipartition because of the low neutrino absorp- this leads to much smaller neutrino absorption rates even 16

(no osc) (a) EoS: DD2 Y (b) polar ejecta (60◦ < θej 90◦ 90◦) e 8 | − |≤ 0.18 0.6 (no osc) Y 0.16 7 e (osc) 0.5 Ye 0.14 ) 6 0.12 0.4 ⊙ 0.1 M 5 5 − = 100km) 0.08 0.3 4 r ( 0.06 (10 e 3 ej

Y 0.2 0.04 ∆ M 2 0.02 0.1 0 1 -0.02 0 0 0 20 40 60 80 100 120 140 160 180 0 0.1 0.2 0.3 0.4 0.5 0.6 θ (degree) ej Ye(r = 100 km) (c) middle ejecta (30◦ < θej 90◦ 60◦) (d) equatorial ejecta (0◦ < θ 90◦ 30◦) 12 | − |≤ 25 | ej − |≤ (no osc) (no osc) Y Y e (osc) e (osc) Y Y 10 e 20 e ) ) ⊙ 8 ⊙ M M 15 5 5

− 6 −

(10 (10 10

ej 4 ej M M 2 5

0 0 0 0.1 0.2 0.3 0.4 0.5 0.6 0 0.1 0.2 0.3 0.4 0.5 0.6 Ye(r = 100 km) Ye(r = 100 km)

FIG. 14. Impact of neutrino flavor equipartition on the ejecta Ye for the 1.35 + 1.35 M model with DD2 EoS. Panel (a) (osc) (no osc) shows the distribution ∆Ye = Ye − Ye , at r = 100 km, as a function of θej for all tracer particles. The color bar represents Ye without flavor conversions. Panels (b)–(d) show the histograms of the ejecta mass distributions as functions of Ye at r = 100 km, for the polar, middle, and equatorial ejecta, respectively. without assuming neutrino flavor equipartition. As a V. SUMMARY AND DISCUSSION consequence, since the ejecta start out as neutron-rich material, a largely reduced Y for the neutrino-driven e In this work, we have examined the neutrino emission outflow was found in the BH–torus model of Ref. [34]. properties and the conditions for the occurrence of fast neutrino flavor conversions during the first 10 ms after the coalescence of symmetric (1.35+1.35 M ) and asym- We show in Fig. 16 the impact of flavor equiparti- metric (1.25 + 1.45 M ) NS binaries, during which the tion on the abundance distribution for the polar ejecta remnant consists of a hypermassive NS surrounded by an ◦ ◦ ◦ accretion disk. We have also performed detailed analyses (60 θej 90 90 ) for the 1.35 + 1.35 M model with≤ both | EoSs.− | Since ≤ flavor equipartition mainly shifts regarding the properties of the material ejected during this phase and nucleosynthesis calculations for cases with the distribution of high Ye & 0.3 material, noticeable changes appear regarding the iron peak and the first and without neutrino flavor mixing. Our study is based peak elements. In particular, the amount of produced on the outputs from the general-relativistic simulations A = 56 nuclei is enhanced by a factor of 6(27) for the with approximate neutrino transport, performed with DD2(SFHo) EoS model. The amount of∼ lanthanides in two different EoS (DD2 and SHFo) based on Refs. [23, 44] the polar ejecta, relevant to the kilonova color, is affected and available at [45]. negligibly for the DD2 model and reduced by a factor Flavor instabilities that may lead to fast pairwise fla- of 2 for the SHFo model. Since the polar ejecta contribute∼ vor conversions occur throughout the whole investigated up to 10% to the total ejecta mass considered here, the post-merger evolution, independent of the adopted EoS overall∼ modifications induced by flavor equipartition on or the mass ratio of the binary. This is a direct conse- the total abundance yields are relatively small. For in- quence of theν ¯e emission dominating over the νe one due stance, the change of the produced amount of strontium to the protonization of the merger remnant, which leads at the time of a day and the amount of free neutrons at to crossings in the neutrino ELN angular distributions the end of the r-process is at the level of 10%. everywhere above the neutrino emitting surfaces. Our ∼ 17

(no osc) (a) EoS: SFHo Y (b) polar ejecta (60◦ < θej 90◦ 90◦) e 8 | − |≤ 0.18 0.6 (no osc) Y 0.16 7 e (osc) 0.5 Ye 0.14 ) 6 0.12 0.4 ⊙ 0.1 M 5 5 − = 100km) 0.08 0.3 4 r ( 0.06 (10 e 3 ej

Y 0.2 0.04 ∆ M 2 0.02 0.1 0 1 -0.02 0 0 0 20 40 60 80 100 120 140 160 180 0 0.1 0.2 0.3 0.4 0.5 0.6 θ (degree) ej Ye(r = 100 km) (c) middle ejecta (30◦ < θej 90◦ 60◦) (d) equatorial ejecta (0◦ < θ 90◦ 30◦) 18 | − |≤ 30 | ej − |≤ (no osc) (no osc) Y Y 16 e (osc) e (osc) Ye 25 Ye

) 14 ) ⊙ 12 ⊙ 20 M M

5 10 5

− − 15 8 (10 (10

ej 6 ej 10

M 4 M 5 2 0 0 0 0.1 0.2 0.3 0.4 0.5 0.6 0 0.1 0.2 0.3 0.4 0.5 0.6 Ye(r = 100 km) Ye(r = 100 km)

FIG. 15. Same as Fig. 14, but for the 1.35 + 1.35 M model with SHFo EoS.

results thus confirm the earlier conclusions of Ref. [33], powered by the decay of free neutrons [62]. which adopted a simple toy-model for the neutrino emis- By relying on the extreme ansatz that fast pairwise sion characteristics and emission surface geometry. How- conversions lead to flavor equilibration, we find that fla- ever, in contrast to the results obtained in Ref. [34], vor mixing of neutrinos mostly affects the polar ejecta which showed that the region where flavor instabilities ◦ within 30 by changing the peak Ye from 0.4 to exist shrinks on a time scale of (10) ms as the BH– ∼ ∼ ∼ O 0.5. The dominant effect is thus to reduce the first disk remnant evolves, the flavor unstable regions reported peak∼ abundances while enhancing the iron peak abun- here remain quite stable within the examined 10 ms of dances. This result quantifies the most extreme impact post-merger evolution. Since Refs. [13, 18] reported dom- of neutrino flavor conversions on the nucleosynthesis in inating emission ofν ¯e over νe on time scales longer than the early-time ejecta; most likely, the effect due to neu- (100) ms for a hypermassive NS accretion disk system, O trino flavor conversions should be between our explored we expect that the flavor instabilities found in this paper two cases with and without oscillations. Note that we may be sustained for even longer duration and affect the have only examined the flavor instability above the neu- nucleosynthesis in the disk winds. trino emitting surfaces. Beyond that, future work inves- As for the ejecta properties and the nucleosynthesis tigating the occurrence of unstable conditions inside the neutrino-trapping regime, along the lines of recent work outcome, the ejecta contain a wide Ye distribution up to 0.5 due to the effect of weak interactions including done in the context of core-collapse supernovae [65–71], neutrino absorption, allowing for the formation of heavy should be carried out. Potential effects due to the pres- elements in all three r-process peaks. In particular, a few ence of turbulence [72] may also play a role in post-merger −5 environments. times 10 M of strontium are synthesized in the ejecta in all models, consistent with the amount inferred from Multidimensional numerical simulations tracking the the GW170817 kilonova observation [11]. We also find flavor evolution in the presence of fast pairwise conver- that the amount of free neutrons left after the r-process sions (see, e.g., Refs. [54–56, 73–76]), including the col- freeze-out is roughly a factor of 10 smaller than the one lisional term in the equations of motion, are essential obtained in simulations without taking into account the to draw robust conclusions on the role of neutrino fla- effect of weak interactions. This has implications for the vor conversion for the outcome of nucleosynthesis in the prediction of the early-time UV emission that may be ejecta and the corresponding kilonova observables. 18

-2 Appendix A: Computing the neutrino number (a) EoS: DD2, polar ejecta only densities from the simulation data -3 As the simulation outputs from Refs. [23, 44] did not

)] -4 directly store the νe andν ¯e absorption rates on free nu- A

( cleons along the ejecta trajectories, we compute the rates

Y -5 [ in a post-processing fashion as follows. First, the simula-

10 tions provide the local νe andν ¯e energy luminosities and -6

log average energies at radii of 50 and 100 km as functions of the polar angle θ (with respect to the z-axis) for different -7 osc post-merger time snapshots. This allows to compute the no osc. νe andν ¯e number densities at these two radii. Then, at a -8 x 40 60 80 100 120 140 160 180 200 220 given time t, for any spatial coordinate on a trajectory with radius 50 km r 100 km (with an angle θ), we mass number, A ≤ ≤ compute the corresponding number densities of νe and -2 ν¯e by linearly interpolating the logarithmic values of the (b) EoS: SFHo, polar ejecta only densities obtained above at r = 50 and 100 km. Once we -3 have the number densities of the νe andν ¯e, the absorp- tion rates for 50 km r 100 km can be computed using ≤ ≤ )] -4 Eqs. (B1) and (B2) in Appendix.B. For r > 100 km, we

A −2 ( extrapolate the rates assuming that they scale as r :

Y -5 [ 2 10 100 km -6 λνα (r > 100 km) = λνα (r = 100 km) , log × r (A1) -7 osc For r < 50 km, we take a different form of extrapolation no osc. to partly account for the finite-size emission geometry: -8 40 60 80 100 120 140 160 180 200 220 λν (r < 50km) = λν (r = 50km) (A2) mass number, A α α × p !2 1 1 (r /r)2 − − 0 , p 2 FIG. 16. Nucleosynthesis outcome in the polar ejecta for the 1 1 (r0/(50km)) cases with and without flavor equipartition for the 1.35 + − − 1.35 M model with DD2 EoS (a) and for the SHFo EoS with r = 25 km. We have checked that even simply 0 model (b). taking the 1/r2 extrapolation for regions with r < 50 km leads to nearly identical results to those obtained by using Eq. (A2). Figure 17 compares the Ye distribution at r = 100 km, obtained by the simulation of Ref. [23] to the one ob- ACKNOWLEDGMENTS tained by using the above neutrino absorption rates in the nuclear reaction network described in Sec.IVA. It shows that the Ye distributions agree with each other MG and MRW acknowledge support from the Ministry reasonably well. of Science and Technology, Taiwan under Grant No. 107- For νx, since the simulations do not store the energy lu- 2119-M-001-038, No. 108-2112-M-001-010, No. 109-2112- minosity and average energy in the bins at 50 and 100 km, M-001-004, the Academia Sinica under project number we estimate the νx number density and average energy AS-CDA-109-M11, and the Physics Division of National along each ejecta trajectory as follows. First, the loca- Center for Theoretical Sciences. IT acknowledges sup- tions of the νx surface and the associated temperatures port from the Villum Foundation (Project No. 13164), for times between 2.5 t 10 ms are interpolated by the Danmarks Frie Forskningsfonds (Project No. 8049- using the data provided≤ at≤ 2.5, 5, 7.5 and 10 ms. Then, 00038B), the Knud Højgaard Foundation. At Garch- for a given time t, a νx number density on each point ~x, ing, funding by the European Research Council through n˜νx (~x), at the emission surface following the Fermi-Dirac grant ERC-AdG No. 341157-COCO2CASA and by the distribution with temperature T (~xi) and zero chemical Deutsche Forschungsgemeinschaft (DFG, German Re- potential can be easily calculated. We then re-normalize search Foundation) through grants SFB-1258 “Neutrinos the total neutrino number luminosity emitted from the and Dark Matter in Astro- and Particle Physics (NDM)” νx surface to the value given by simulation data, and under Germany’s Excellence Strategy through Ex- Z cellence Cluster ORIGINS (EXC 2094)—390783311 is ac- Lνx 5 ξ LN,νx = = · dSn˜νx , (A3) knowledged. Eν 12 h x i 19

nique as in the main text to compute the νx number den- 30 sities for the locations crossed by the trajectories. The simulation absorption rates of the converted ν andν ¯ on nucleons, 25 post-processing x x λ and λ , along all trajectories, are similarly com- ) νx ν¯x ⊙ 20 puted as those of νe andν ¯e given in AppendixB by re- M placing the corresponding number densities, the average 4

− 15 energies, and other higher energy moments. (10

ej 10 M 5 Appendix B: Computing the neutrino absorption rates 0 0 0.1 0.2 0.3 0.4 0.5 0.6 In order to compute the evolution of Ye for the out- Ye(r = 100km) flows for the cases with and without flavor conversions, we first compute the number densities and average ener- gies of ν ,ν ¯ and ν (without flavor conversions) along FIG. 17. Ejecta Ye histograms computed via post-processing e e x and as from the simulations of Refs. [23, 44]. the trajectories of all tracer particles by post-processing the simulation data as detailed in AppendixA. For the case without flavor conversions, we follow Ref. [77] to where Lν and Eν are the energy luminosity and mean calculate the νe andν ¯e absorption on free nucleons: x h x i energy of νx shown in Fig.1, respectively. The quantity 0 dS is the differential surface area on the ν surface. The λν = nνe σνe , (B1) x e h i factor 5/12 accounts for the forward peaked angular pro- 0 λν¯e = nν¯e σν¯e , (B2) file of νx emission consistent with Eq. (13), and ξ is the h i normalization constant. Correspondingly, the rescaled where σνe and σνe are the spectrally averaged absorp- νx number density on their emission surface is given by h i h i tion cross-sections of νe andν ¯e. By taking into account n (~x) = ξn˜ (~x). νx νx the recoil corrections and weak magnetism [78], the av- We assume the νx average energy at each location ~x erage neutrino capture cross sections are approximated on the emission surface to be by 1 " # Eνx (~x) = ( Eνx + 3.15 T (~x)) , (A4)    2 h i 2 h i ∆ ∆ σνe k Eνe ενe 1 + 2 + aνe Wνe , ± h i ' h i ενe ενe to partly account for the fact that the νx-e scatterings, (B3) which can down-scatter νx, was not included in the nu- " # merical simulations of Ref. [23]. Since the local temper-  ∆   ∆ 2 σ k E ε 1 + 2 + a W , ature on the νx emission surface within x . 20 km is ν¯e ν¯e ν¯e ν¯e ν¯e h i ' h i εν¯ εν¯ found to be & 6 MeV, the main effect of the above choice e e (B4) is to reduce the Eνx (~x) at the outer edge of the emis- sion surface. Weh havei additionally confirmed that adopt- −44 2 2 ing a location independent average energy of νx given by where k = 9.3 10 cm /MeV , ενe,ν¯e = 2 × 2 2 Eν does not qualitatively change our results shown in Eν ,ν¯ / Eνe,ν¯e , aνe,ν¯e = Eν ,ν¯ / Eνe,ν¯e , and ∆ = h x i h e e i h i h e e i h i the main text. (mn mp) = 1.293 MeV is the neutron-proton mass dif- − For t < 2.5 ms, we simply assume that the neutrino ference. The weak-magnetism and recoil correction fac- emission surface is the same as the one at t = 2.5 ms, tors Wνe,ν¯e are given by and scale the number density n (~x,t) and the average νx   energy Eν (~x,t) in the following way bνe ενe x Wν = 1 + 1.02 , (B5) h i e M   LN,ν (t)   x b¯ν εν nνx (~x,t) = nνx (~x,t = 2.5 ms) , ¯ e ¯e L (t = 2.5 ms) Wν¯e = 1 7.22 , (B6) N,νx − M (A5)   3 2 2 Eν (t) where bν ,ν = E Eν ,ν / E is the spectral E (~x,t) = E (~x,t = 2.5 ms) h x i . e ¯e νe,ν¯e e ¯e νe,ν¯e νx νx shape factor forh ν (¯ν ih) and Mi h = 940i is roughly the h i h i Eνx (t = 2.5 ms) e e h i (A6) mass of a nucleon in MeV. Note that in deriving the rates through the above equations, we have assumed zero Once we have the desired quantities on the emission chemical potentials for all neutrino species to compute surface for all times, we use the same ray-tracing tech- the i-th neutrino energy moments Ei . h να i 20

∗ [email protected] trophys. J. 848, L12 (2017), arXiv:1710.05833 [astro- † [email protected] ph.HE]. ‡ [email protected] [11] Darach Watson et al., “Identification of strontium in § [email protected] the merger of two neutron stars,” Nature 574, 497–500 ¶ [email protected] (2019), arXiv:1910.10510 [astro-ph.HE]. [1] James M. Lattimer and David N. Schramm, “Black-hole- [12] Shinya Wanajo, Yuichiro Sekiguchi, Nobuya Nishimura, neutron-star collisions,” Astrophys. J. 192, L145–L147 Kenta Kiuchi, Koutarou Kyutoku, and Masaru Shibata, (1974). “Production of all the r-process nuclides in the dynamical [2] David Eichler, Mario Livio, Tsvi Piran, and ejecta of neutron star mergers,” Astrophys. J. 789, L39 David N. Schramm, “Nucleosynthesis, Neutrino Bursts (2014), arXiv:1402.7317 [astro-ph.SR]. and Gamma-Rays from Coalescing Neutron Stars,” Na- [13] Albino Perego, Stephan Rosswog, Ruben M. Cabez´on, ture 340, 126–128 (1989). Oleg Korobkin, Roger K¨appeli, Almudena Arcones, [3] John J. Cowan, Christopher Sneden, James E. and Matthias Liebend¨orfer,“Neutrino-driven winds from Lawler, Ani Aprahamian, Michael Wiescher, Karl- neutron star merger remnants,” Mon. Not. Roy. As- heinz Langanke, Gabriel Mart´ınez-Pinedo, and tron. Soc. 443, 3134–3156 (2014), arXiv:1405.6730 [astro- Friedrich-Karl Thielemann, “Origin of the Heaviest El- ph.HE]. ements: the Rapid Neutron-Capture Process,” (2019), [14] Rodrigo Fern´andezand Brian D. Metzger, “Delayed out- arXiv:1901.01410 [astro-ph.HE]. flows from black hole accretion tori following neutron star [4] Li-Xin Li and Bohdan Paczynski, “Transient events from binary coalescence,” Mon. Not. Roy. Astron. Soc. 435, neutron star mergers,” Astrophys. J. 507, L59 (1998), 502 (2013), arXiv:1304.6720 [astro-ph.HE]. arXiv:astro-ph/9807272 [astro-ph]. [15] Yuichiro Sekiguchi, Kenta Kiuchi, Koutarou Kyutoku, [5] Shrinivas R. Kulkarni, “Modeling supernova-like explo- and Masaru Shibata, “Dynamical mass ejection from bi- sions associated with gamma-ray bursts with short dura- nary neutron star mergers: Radiation-hydrodynamics tions,” (2005), arXiv:astro-ph/0510256 [astro-ph]. study in general relativity,” Phys. Rev. D91, 064059 [6] B. D. Metzger, G. Mart´ınez-Pinedo, S. Darbha, (2015), arXiv:1502.06660 [astro-ph.HE]. E. Quataert, A. Arcones, D. Kasen, R. Thomas, P. Nu- [16] David Radice, Filippo Galeazzi, Jonas Lippuner, Luke F. gent, I. V. Panov, and N. T. Zinner, “Electromag- Roberts, Christian D. Ott, and Luciano Rezzolla, “Dy- netic Counterparts of Compact Object Mergers Powered namical Mass Ejection from Binary Neutron Star Merg- by the Radioactive Decay of R-process Nuclei,” Mon. ers,” Mon. Not. Roy. Astron. Soc. 460, 3255–3271 (2016), Not. Roy. Astron. Soc. 406, 2650 (2010), arXiv:1001.5029 arXiv:1601.02426 [astro-ph.HE]. [astro-ph.HE]. [17] Francois Foucart, Evan O’Connor, Luke Roberts, [7] Brian D. Metzger, “Kilonovae,” Living Rev. Rel. 20, 3 Lawrence E. Kidder, Harald P. Pfeiffer, and Mark A. (2017), arXiv:1610.09381 [astro-ph.HE]. Scheel, “Impact of an improved neutrino energy esti- [8] Benjamin P. Abbott et al. (Virgo and LIGO Scientific mate on outflows in neutron star merger simulations,” Collaborations), “GW170817: Observation of Gravita- Phys. Rev. D 94, 123016 (2016), arXiv:1607.07450 [astro- tional Waves from a Binary Neutron Star Inspiral,” Phys. ph.HE]. Rev. Lett. 119, 161101 (2017), arXiv:1710.05832 [gr-qc]. [18] Sho Fujibayashi, Kenta Kiuchi, Nobuya Nishimura, [9] Benjamin P. Abbott et al. (Fermi-GBM, INTEGRAL, Yuichiro Sekiguchi, and Masaru Shibata, “Mass Ejec- Virgo and LIGO Scientific Collaborations), “Gravita- tion from the Remnant of a Binary Neutron Star Merger: tional Waves and Gamma-Rays from a Binary Neutron Viscous-Radiation Hydrodynamics Study,” Astrophys. J. Star Merger: GW170817 and GRB 170817A,” Astrophys. 860, 64 (2018), arXiv:1711.02093 [astro-ph.HE]. J. 848, L13 (2017), arXiv:1710.05834 [astro-ph.HE]. [19] Jonah M. Miller, Benjamin R. Ryan, Joshua C. Do- [10] B. P. Abbott et al. (GROND, SALT Group, Oz- lence, Adam Burrows, Christopher J. Fontes, Christo- Grav, DFN, INTEGRAL, Virgo, Insight-Hxmt, MAXI pher L. Fryer, Oleg Korobkin, Jonas Lippuner, Team, Fermi-LAT, J-GEM, RATIR, IceCube, CAAS- Matthew R. Mumpower, and Ryan T. Wollaeger, TRO, LWA, ePESSTO, GRAWITA, RIMAS, SKA “Full Transport Model of GW170817-Like Disk Produces South Africa/MeerKAT, H.E.S.S., 1M2H Team, IKI-GW a Blue Kilonova,” Phys. Rev. D100, 023008 (2019), Follow-up, Fermi GBM, Pi of Sky, DWF (Deeper Wider arXiv:1905.07477 [astro-ph.HE]. Faster Program), Dark Energy Survey, MASTER, As- [20] Trevor Vincent, Francois Foucart, Matthew D. Duez, troSat Cadmium Zinc Telluride Imager Team, Swift, Roland Haas, Lawrence E. Kidder, Harald P. Pfeiffer, Pierre Auger, ASKAP, VINROUGE, JAGWAR, Chan- and Mark A. Scheel, “Unequal Mass Binary Neutron dra Team at McGill University, TTU-NRAO, GROWTH, Star Simulations with Neutrino Transport: Ejecta and AGILE Team, MWA, ATCA, AST3, TOROS, Pan- Neutrino Emission,” Phys. Rev. D 101, 044053 (2020), STARRS, NuSTAR, ATLAS Telescopes, BOOTES, Cal- arXiv:1908.00655 [gr-qc]. techNRAO, LIGO Scientific, High Time Resolution Uni- [21] S. Wanajo, “The rp-process in neutrino-driven winds,” verse Survey, Nordic Optical Telescope, Las Cumbres Astrophys. J. 647, 1323–1340 (2006), arXiv:astro- Observatory Group, TZAC Consortium, LOFAR, IPN, ph/0602488. DLT40, Texas Tech University, HAWC, ANTARES, KU, [22] S. Goriely, A. Bauswein, H.-T. Janka, S. Panebianco, J-L. Dark Energy Camera GW-EM, CALET, Euro VLBI Sida, J-F. Lemaˆıtre,S. Hilaire, and N. Dubray, “The r- Team, and ALMA Collaborations), “Multi-messenger process nucleosynthesis during the decompression of neu- Observations of a Binary Neutron Star Merger,” As- tron star crust material,” J. Phys. Conf. Ser. 665, 012052 21

(2016). [38] Ignacio Izaguirre, Georg G. Raffelt, and Irene Tamborra, [23] Ricard Ardevol-Pulpillo, H.-Thomas Janka, Oliver “Fast Pairwise Conversion of Supernova Neutrinos: A Just, and Andreas Bauswein, “Improved Leakage- Dispersion-Relation Approach,” Phys. Rev. Lett. 118, Equilibration-Absorption Scheme (ILEAS) for Neutrino 021101 (2017), arXiv:1610.01612 [hep-ph]. Physics in Compact Object Mergers,” Mon. Not. Roy. [39] Francesco Capozzi, Basudeb Dasgupta, Eligio Lisi, An- Astron. Soc. 485, 4754–4789 (2019), arXiv:1808.00006 tonio Marrone, and Alessandro Mirizzi, “Fast flavor con- [astro-ph.HE]. versions of supernova neutrinos: Classifying instabilities [24] Stan E. Woosley, “Gamma-ray bursts from stellar mass via dispersion relations,” Phys. Rev. D96, 043016 (2017), accretion disks around black holes,” Astrophys. J. 405, arXiv:1706.03360 [hep-ph]. 273 (1993). [40] Sajad Abbar and Huaiyu Duan, “Fast neutrino flavor [25] Maximillian Ruffert and H.-Thomas Janka, “Gamma- conversion: roles of dense matter and spectrum cross- ray bursts from accreting black holes in neutron star ing,” Phys. Rev. D98, 043014 (2018), arXiv:1712.07013 mergers,” Astron. Astrophys. 344, 573–606 (1999), [hep-ph]. arXiv:astro-ph/9809280 [astro-ph]. [41] Changhao Yi, Lei Ma, Joshua D. Martin, and Huaiyu [26] Ivan Zalamea and Andrei M. Beloborodov, “Neu- Duan, “Dispersion relation of the fast neutrino oscillation trino heating near hyper-accreting black holes,” Mon. wave,” Phys. Rev. D99, 063005 (2019), arXiv:1901.01546 Not. Roy. Astron. Soc. 410, 2302–2308 (2011), [hep-ph]. arXiv:1003.0710 [astro-ph.HE]. [42] Francesco Capozzi, Georg G. Raffelt, and Tobias [27] Oliver Just, Martin Obergaulinger, H.-Thomas Janka, Stirner, “Fast Neutrino Flavor Conversion: Collective Andreas Bauswein, and Nicole Schwarz, “Neutron- Motion vs. Decoherence,” JCAP 1909, 002 (2019), star merger ejecta as obstacles to neutrino-powered jets arXiv:1906.08794 [hep-ph]. of gamma-ray bursts,” Astrophys. J. 816, L30 (2016), [43] Madhurima Chakraborty and Sovan Chakraborty, arXiv:1510.04288 [astro-ph.HE]. “Three flavor neutrino conversions in supernovae: [28] A. Malkus, J. P. Kneller, G. C. McLaughlin, and R. Sur- slow & fast instabilities,” JCAP 01, 005 (2020), man, “Neutrino oscillations above black hole accretion arXiv:1909.10420 [hep-ph]. disks: disks with electron-flavor emission,” Phys. Rev. [44] Ricard Ardevol-Pulpillo, “A new scheme to treat neu- D86, 085015 (2012), arXiv:1207.6648 [hep-ph]. trino effects in neutron-star mergers: implementation, [29] A. Malkus, A. Friedland, and G. C. McLaughlin, tests and applications,” (2018), PhD Thesis, Technis- “Matter-Neutrino Resonance Above Merging Compact che Universit¨atM¨unchen. Objects,” (2014), arXiv:1403.5797 [hep-ph]. [45] MPA Data Archive, https://wwwmpa.mpa-garching. [30] Meng-Ru Wu, Huaiyu Duan, and Yong-Zhong Qian, mpg.de/ccsnarchive/data/Ardevol-Pulpillo2019/. “Physics of neutrino flavor transformation through [46] Andreas Bauswein, Stephane Goriely, and H.-Thomas matter-neutrino resonances,” Phys. Lett. B752, 89–94 Janka, “Systematics of dynamical mass ejection, nucle- (2016), arXiv:1509.08975 [hep-ph]. osynthesis, and radioactively powered electromagnetic [31] Yong-Lin Zhu, Albino Perego, and Gail C. McLaughlin, signals from neutron-star mergers,” Astrophys. J. 773, “Matter Neutrino Resonance Transitions above a Neu- 78 (2013), arXiv:1302.6530 [astro-ph.SR]. tron Star Merger Remnant,” Phys. Rev. D94, 105006 [47] S. Typel, G. R¨opke, T. Klahn, D. Blaschke, and H. H. (2016), arXiv:1607.04671 [hep-ph]. Wolter, “Composition and thermodynamics of nuclear [32] Maik Frensel, Meng-Ru Wu, Cristina Volpe, and Albino matter with light clusters,” Phys. Rev. C81, 015803 Perego, “Neutrino Flavor Evolution in Binary Neutron (2010), arXiv:0908.2344 [nucl-th]. Star Merger Remnants,” Phys. Rev. D95, 023011 (2017), [48] Matthias Hempel, Tobias Fischer, J¨urgen Schaffner- arXiv:1607.05938 [astro-ph.HE]. Bielich, and Matthias Liebend¨orfer,“New Equations [33] Meng-Ru Wu and Irene Tamborra, “Fast neu- of State in Simulations of Core-Collapse Supernovae,” trino conversions: Ubiquitous in compact binary Astrophys. J. 748, 70 (2012), arXiv:1108.0848 [astro- merger remnants,” Phys. Rev. D95, 103007 (2017), ph.HE]. arXiv:1701.06580 [astro-ph.HE]. [49] Andrew W. Steiner, Matthias Hempel, and Tobias Fis- [34] Meng-Ru Wu, Irene Tamborra, Oliver Just, and H.- cher, “Core-collapse supernova equations of state based Thomas Janka, “Imprints of neutrino-pair flavor con- on neutron star observations,” Astrophys. J. 774, 17 versions on nucleosynthesis in ejecta from neutron-star (2013), arXiv:1207.2184 [astro-ph.SR]. merger remnants,” Phys. Rev. D96, 123015 (2017), [50] H.-Thomas Janka, “When do supernova neutrinos of dif- arXiv:1711.00477 [astro-ph.HE]. ferent flavors have similar luminosities but different spec- [35] Raymond F. Sawyer, “Speed-up of neutrino transforma- tra?” Astropart. Phys. 3, 377–384 (1995), arXiv:astro- tions in a supernova environment,” Phys. Rev. D72, ph/9503068. 045003 (2005), arXiv:hep-ph/0503013 [hep-ph]. [51] G¨unter Sigl and Georg G. Raffelt, “General kinetic de- [36] Raymond F. Sawyer, “Neutrino cloud instabilities just scription of relativistic mixed neutrinos,” Nucl.Phys. above the neutrino sphere of a supernova,” Phys. B406, 423–451 (1993). Rev. Lett. 116, 081101 (2016), arXiv:1509.03323 [astro- [52] Sovan Chakraborty, Rasmus Sloth Hansen, Ignacio Iza- ph.HE]. guirre, and Georg G. Raffelt, “Self-induced neutrino fla- [37] Basudeb Dasgupta, Alessandro Mirizzi, and Manibrata vor conversion without flavor mixing,” JCAP 1603, 042 Sen, “Fast neutrino flavor conversions near the super- (2016), arXiv:1602.00698 [hep-ph]. nova core with realistic flavor-dependent angular dis- [53] Sagar Airen, Francesco Capozzi, Sovan Chakraborty, Ba- tributions,” JCAP 1702, 019 (2017), arXiv:1609.00528 sudeb Dasgupta, Georg G. Raffelt, and Tobias Stirner, [hep-ph]. “Normal-mode Analysis for Collective Neutrino Oscilla- tions,” JCAP 12, 019 (2018), arXiv:1809.09137 [hep-ph]. 22

[54] Shashank Shalgar and Irene Tamborra, “Dispelling a metric Neutrino Distributions,” Phys. Rev. D 101, myth on dense neutrino media: fast pairwise conversions 043016 (2020), arXiv:1911.01983 [astro-ph.HE]. depend on energy,” (2020), arXiv:2007.07926 [astro- [66] Milad Delfan Azari, Shoichi Yamada, Taiki Morinaga, ph.HE]. Hiroki Nagakura, Shun Furusawa, Akira Harada, Hiro- [55] Ian Padilla-Gay, Shashank Shalgar, and Irene Tamborra, tada Okawa, Wakana Iwakami, and Kohsuke Sumiyoshi, “Multi-Dimensional Solution of Fast Neutrino Conver- “Fast collective neutrino oscillations inside the neutrino sions in Binary Neutron Star Merger Remnants,” (2020), sphere in core-collapse supernovae,” Phys. Rev. D 101, arXiv:2009.01843 [astro-ph.HE]. 023018 (2020), arXiv:1910.06176 [astro-ph.HE]. [56] Shashank Shalgar, Ian Padilla-Gay, and Irene Tamborra, [67] Taiki Morinaga, Hiroki Nagakura, Chinami Kato, and “Neutrino propagation hinders fast pairwise flavor con- Shoichi Yamada, “Fast neutrino-flavor conversion in the versions,” JCAP 06, 048 (2020), arXiv:1911.09110 [astro- preshock region of core-collapse supernovae,” Phys. Rev. ph.HE]. Res. 2, 012046 (2020), arXiv:1909.13131 [astro-ph.HE]. [57] Francesco Capozzi, Basudeb Dasgupta, Alessandro Mi- [68] Hiroki Nagakura, Taiki Morinaga, Chinami Kato, and rizzi, Manibrata Sen, and G¨unter Sigl, “Collisional trig- Shoichi Yamada, “Fast-pairwise Collective Neutrino Os- gering of fast flavor conversions of supernova neutrinos,” cillations Associated with Asymmetric Neutrino Emis- Phys. Rev. Lett. 122, 091101 (2019), arXiv:1808.06618 sions in Core-collapse Supernovae,” Astrophys. J. 886, [hep-ph]. 139 (2019), arXiv:1910.04288 [astro-ph.HE]. [58] Meng-Ru Wu, Rodrigo Fern´andez, Gabriel Mart´ınez- [69] Robert Glas, H.-Thomas Janka, Francesco Capozzi, Pinedo, and Brian D. Metzger, “Production of the entire Manibrata Sen, Basudeb Dasgupta, Alessandro Mirizzi, range of r-process nuclides by black hole accretion disc and G¨unter Sigl, “Fast Neutrino Flavor Instability in outflows from neutron star mergers,” Mon. Not. Roy. the Neutron-star Convection Layer of Three-dimensional Astron. Soc. 463, 2323–2334 (2016), arXiv:1607.05290 Supernova Models,” Phys. Rev. D 101, 063001 (2020), [astro-ph.HE]. arXiv:1912.00274 [astro-ph.HE]. [59] Meng-Ru Wu, Jennifer Barnes, Gabriel Mart´ınez-Pinedo, [70] Lucas Johns, Hiroki Nagakura, George M. Fuller, and and Brian D. Metzger, “Fingerprints of heavy element Adam Burrows, “Neutrino oscillations in supernovae: an- nucleosynthesis in the late-time lightcurves of kilonovae,” gular moments and fast instabilities,” Phys. Rev. D 101, Phys. Rev. Lett. 122, 062701 (2019), arXiv:1808.10459 043009 (2020), arXiv:1910.05682 [hep-ph]. [astro-ph.HE]. [71] Sajad Abbar, “Searching for Fast Neutrino Flavor Con- [60] Joel de Jesus Mendoza-Temis, Meng-Ru Wu, Gabriel version Modes in Core-collapse Supernova Simulations,” Mart´ınez-Pinedo, Karlheinz Langanke, Andreas JCAP 05, 027 (2020), arXiv:2003.00969 [astro-ph.HE]. Bauswein, and H.-Thomas Janka, “Nuclear robustness [72] Sajad Abbar, “Turbulence Fingerprint on Collec- of the r process in neutron-star mergers,” Phys. Rev. tive Oscillations of Supernova Neutrinos,” (2020), C92, 055805 (2015), arXiv:1409.6135 [astro-ph.HE]. arXiv:2007.13655 [astro-ph.HE]. [61] John J. Cowan, B. Pfeiffer, K.L. Kratz, Friedrich K. [73] Sajad Abbar and Maria Cristina Volpe, “On Fast Neu- Thielemann, Christopher Sneden, Scott Burles, David trino Flavor Conversion Modes in the Nonlinear Regime,” Tytler, and Timothy C. Beers, “R-process abundances Phys. Lett. B790, 545–550 (2019), arXiv:1811.04215 and chronometers in metal-poor stars,” Astrophys. J. [astro-ph.HE]. 521, 194 (1999), arXiv:astro-ph/9808272. [74] Joshua D. Martin, Changhao Yi, and Huaiyu [62] Brian D. Metzger and Rodrigo Fern´andez, “Red or Duan, “Dynamic fast flavor oscillation waves in dense blue? A potential kilonova imprint of the delay until neutrino gases,” Phys. Lett. B800, 135088 (2020), black hole formation following a neutron star merger,” arXiv:1909.05225 [hep-ph]. Mon. Not. Roy. Astron. Soc. 441, 3444–3453 (2014), [75] Soumya Bhattacharyya and Basudeb Dasgupta, “Late- arXiv:1402.4803 [astro-ph.HE]. time behavior of fast neutrino oscillations,” Phys. Rev. [63] Luke Bovard, Dirk Martin, Federico Guercilena, Almu- D 102, 063018 (2020), arXiv:2005.00459 [hep-ph]. dena Arcones, Luciano Rezzolla, and Oleg Korobkin, [76] Francesco Capozzi, Madhurima Chakraborty, Sovan “r-process nucleosynthesis from matter ejected in binary Chakraborty, and Manibrata Sen, “Fast flavor con- neutron star mergers,” Phys. Rev. D 96, 124005 (2017), versions in supernovae: the rise of mu-tau neutrinos,” arXiv:1709.09630 [gr-qc]. (2020), arXiv:2005.14204 [hep-ph]. [64] Zewei Xiong, Andre Sieverding, Manibrata Sen, and [77] Else Pllumbi, Irene Tamborra, Shinya Wanajo, H.- Yong-Zhong Qian, “Potential Impact of Fast Fla- Thomas Janka, and Lorenz H¨udepohl, “Impact of neu- vor Oscillations on Neutrino-driven Winds and Their trino flavor oscillations on the neutrino-driven wind nu- Nucleosynthesis,” Astrophys. J. 900, 144 (2020), cleosynthesis of an electron-capture supernova,” Astro- arXiv:2006.11414 [astro-ph.HE]. phys. J. 808, 188 (2015), arXiv:1406.2596 [astro-ph.SR]. [65] Sajad Abbar, Huaiyu Duan, Kohsuke Sumiyoshi, To- [78] Charles J. Horowitz and Gang Li, “Charge conjuga- moya Takiwaki, and Maria Cristina Volpe, “Fast Neu- tion violating interactions in supernovae and nucleosyn- trino Flavor Conversion Modes in Multidimensional thesis,” Phys. Rev. Lett. 82, 5198 (1999), arXiv:astro- Core-collapse Supernova Models: the Role of the Asym- ph/9904171.