Arxiv:1112.0335V1 [Astro-Ph.HE] 1 Dec 2011 Ini Nitnesuc Fnurnso L Aos It ﬂavors

Home , Nuclear density

arXiv:1112.0335v1 [astro-ph.HE] 1 Dec 2011 ini nitnesuc fnurnso l aos It flavors. all of explo- neutrinos supernova of 3 source II the intense type emits an a is in sion core-collapse to sequent r nterange temperature after the and in PNS density are the matter the within number where behind core-bounce, lepton left the gradients entropy along convection and becomes and neu- by diffusion contracts, driven PNS trino are it the emission neutrino of as concomitant Cooling the seconds a and cools. on of and neutron-rich radiation tens increasingly neutrino of as scale collapse time during gained ergy togyaetn h etiosga copnigPNS accompanying signal thereby neutrino diffuse, al- the neutrinos and affecting which strongly PNS along the gradients in the number ters lepton and energy transport signal. neutrino observable den- affects the thereby high and at convection neu- influences matter) significantly dense nuclear sity between symmetric difference and energy matter the tron (which of energy measure find symmetry a We nuclear is ex- the included. of [4], is behavior convection in the of described that treatment that framework a to that basic similar cept The is evolution signal. PNS char- neutrino for temporal the the to of matter acteristics the dense to in of modifications opacities importance medium neutrino the of effect gauge un- the to and the us We convection allows with This consistent EoSs. EoS. both two derlying transport, for dif- neutrino times with fusive late along transport to convective approximate out include evolution PNS of 7–9]. the diffusion [4, affect neutrino convection either or that influencing properties by timescale microphysical mat- cooling other dense in nuclear and opacities the neutrino ter, as (EoS), such PNS, state the of of equation properties the by shaped etdt nunetenurn msina al time early at ex- emission are neutrino material the (i.e. of influence back to fall pected explosion supernova associated the and While mechanism [1–6]. respectively MeV, h o n es rt-eto tr(N)br sub- born (PNS) star proto-neutron dense and hot The ag cl ovcieoetr fmtra ildirectly will material of overturn convective scale Large models hydrostatic one-dimensional present we Here, rt-eto trCoigwt ovcin h ffc fth of Effect The Convection: with Cooling Star Proto-Neutron t . otbuc)telt ienurn inlis signal neutrino time late the bounce) post s 1 1 − eateto srnm n srpyis nvriyo Ca of University Astrophysics, and Astronomy of Department ASnmes 65.x 66.c 16.n 95.85.Ry 21.65.Mn, 26.60.-c, 26.50.+x, numbers: PACS br a a is observable. energy there clearly symmetry exten star, be nuclear spatial may proto-neutron the and that the of period within behavior ends time the convection the to find sensitive We is density. baryon high osuyissniiiyt h qaino tt,nurn o neutrino state, of equation the to sensitivity its study to 3 eatmn eFıiaAlcd,UiesttdAaat A d’Alacant, Universitat F´ısica Aplicada, de Departament 5 emdlnurn msinfo el onnurnsa sub star neutron born newly a from emission neutrino model We ..Roberts, L.F. × ρ 4 2 = 10 nttt o ula hoy nvriyo ahntn Se Washington, of University Theory, Nuclear for Institute 2 53 hoeia iiin o lmsNtoa aoaoy Los Laboratory, National Alamos Los Division, Theoretical − rso rvttoa idn en- binding gravitational of ergs 6 × 1 10 .Shen, G. 14 g/cm 3 2 and .Cirigliano, V. T 5 = − 40 2 ..Pons, J.A. γ where and epciey ovcieisaiiyst nfor matter, in dense sets in leptons instability Convective of fraction num- respectively. the baryon and baryon, density per ber entropy pressure, the gravity, iinadeeg ynurns[,1,1] al one- fingering Early neutron that 17]. 16, suggested [2, [2] compo- neutrinos work of by transport dimensional energy lateral and the to sition due instabilities sive and PNS. EoS the nuclear there within the that convection between of implies interplay patterns clearly strong the rate a growth be the will of form The h iergot ae r bandfo h standard the as from analysis obtained stability are Ledoux rates growth length. mixing linear chosen The the results to sensitive our particularly simulations, not pro- are our rea- convection throughout a since efficiently is Further, ceeds this that approximation. suggesting sonable one-dimension, mixing simple in a in scheme by simulations reproduced hydrodynamic well that were In shown two-dimensions was [2]. it theory [15] length Ref. mixing dependent time standard by signal. neutrino observable the opacities shaping modified in convection between and interplay mi- the underlying the and first to the crophysics is connection work its at this exploring 11], convection at [2, attempt of Refs. in role studied successful was the times a Although late to 11–15]. lead [2, to increase explosion could enough it luminosities that neutrino hope the the the with one dimensions, both beneath two in instability extensively and studied of been period has spheres early neutrino gra- This entropy negative outer to [10]. pas- due dients the the shock, after supernova that the soon of convection recognized sage to been unstable is long mantle has PNS It formation. n nadto,tePSmyb ujc odul diffu- doubly to subject be may PNS the addition, In emdlcneto sadffsv rcs described process diffusive a as convection model We B g = , P 3 aiis n ovcieisaiiisat instabilities convective and pacities, a ntepeitdnurn emission neutrino predicted the in eak .Cres9,000Aaat Spain Alacant, 03080 99, Correus p. ω , ∂ .Reddy, S. ∂ 2 iona at rz A904USA 95064 CA Cruz, Santa lifornia, s n bv ula est.When density. nuclear above and t ln ln , = n n vrwihcneto operates convection which over t P eun oasproaexplosion supernova a to sequent B − B tl,Wsigo 98195 Washington attle, and , γ n g lms M87544 NM Alamos, s,Y B ,4 2, L ( Y γ γ L n ..Woosley S.E. and s s ∇ = r h oa ceeaindeto due acceleration local the are ln( ymtyEnergy Symmetry e s ∂ ∂ + ) ln ln P s γ Y L n ∇ B ,Y ln( 1 L γ Y Y L L )) = , ω 2 ∂ ∂ > ln ln (1) Y 0. P L n B ,s 2

35 IAS is obtained by varying values of the poorly known three- 30 100 HIC neutron interaction [20]. 25 PDR (MeV) 20 Recent work has shown that the derivative of the sym- sym ′

E’ 15 0

n metry energy with respect to density, denoted as Esym = 10 Skin (Sn) 5 ∂Esym/∂nB, plays a crucial role both in the terrestrial 30 35

(MeV) Esym(n0 ) (MeV)

sym context where it affects the neutron density distribution E 50 in neutron-rich nuclei and in astrophysics where it affects the structure and thermal evolution of neutron stars (for IU-FSU a recent review see Ref. [22]). The pressure of neutron GM3 2 ′ QMC matter at saturation density, Pneutron(n0)= n0Esym, in- 0 fluences the radii of cold neutron stars [23]. In neutron- 0 0.1 0.2 0.3 0.4 0.5 -3 rich nuclei, the neutron-skin thickness is also sensitive n B (fm ) ′ to Esym(ρ0), so that there exists a linear correlation between the neutron-skin thickness and neutron star radius FIG. 1: The symmetry energy as function of density for the ′ [24]. IU-FSU and GM3 EoSs. Inset: n0Esym versus Esym at nuclear saturation density, for IU-FSU (circle), GM3 (square), and To study the sensitivity of PNS evolution to the nu- QMC (diamonds). The shaded regions correspond to various clear symmetry energy we employ two RMF models with ′ experimental constraints taken from Ref. [18]. different predictions for Esym(ρ0). The first EoS is the IU-FSU EoS taken from [19], which includes a non-linear coupling between the vector and iso-vector mesons that instabilities and convection enhanced the neutrino lumi- allows the symmetry energy to be tuned at high den- nosity to successfully power a neutrino driven explosion. sity. The second EoS employed is the GM3 parameter However, more recent two dimensional studies found no set, where non-linear coupling of the vector meson fields evidence of these doubly diffusive instabilities [14, 15]. is neglected [21]. The symmetry energy as a function of Because of this and the increased complexity of treating density is shown in Fig. 1 for the two EoS. The inset in the doubly diffusive instabilities, we do not include them Fig. 1 shows current theoretical estimates and experimen- ′ in our study. tal constraints on Esym and n0Esym at nuclear density. ′ The EoS and neutrino interaction rates are modeled In the rest of this letter, we demonstrate that Esym(ρ0) using a relativistic mean field (RMF) model of nuclear plays a role in stabilizing PNS convection at late times interactions. We adopt a non-linear generalization of and thereby directly affects the PNS neutrino signal. The the original Walecka model described in [19]. Here, logarithmic derivatives γs and γnB are always positive, so the nucleon-nucleon interaction energy is calculated in that negative entropy gradients always provide a destabi- the mean field approximation using effective interactions, lizing influence. For given entropy and lepton gradients, which are tuned to reproduce gross observed properties stability is then determined by the ratio γYL /γs. The of nuclei and empirical properties of symmetric nuclear sign and magnitude of γYL is strongly influenced by the matter at saturation density. Although these empiri- density dependence of the nuclear asymmetry energy, so cal constraints provide valuable guidance to constrain that negative gradients in lepton number can be either aspects of the symmetric nuclear EoS at nuclear den- stabilizing or destabilizing and the degree to which they sities, the experimental constraints on the properties of are stabilizing varies from EoS to EoS. To clarify this we neutron-rich matter are relatively weak. The difference note that at T = 0 and when the neutrino contribution between the energy of symmetric matter (equal num- to the pressure is small bers of neutrons and protons) and pure neutron matter is ∂P ≃ 4/3 1/3 − 2 ′ − called the symmetry energy, Esym(nB ), and is defined by nB Ye 4nBEsym(1 2Ye), (2) 2 ∂YL nB E(nB, xp)= E(nB, xp =1/2)+ Esym(nB)δ + ···. Here, δ = (1 − 2xp) and E(nB , xp) is the energy per particle of which is a reasonable approximation to the finite tem- uniform matter composed of neutrons and protons with perature result. The first term comes from the electron total baryon density nB and proton fraction xp. In charge contribution to the pressure, while the second term is neutral matter xp = Ye where Ye is the electron fraction. due to nucleons and is negative since both the Fermi and Various experimental probes of the nuclear symmetry en- interaction energies favor a symmetric state. For high ergy and its density dependence in nuclei and heavy-ion densities and low electron fractions, for realistic values ′ collisions are actively being pursued in terrestrial exper- of Esym, this leads to negative γYL . Therefore, a larger ′ iments, but are yet to yield strong constraints. These Esym leads to negative lepton gradients in the PNS pro- ′ constraints are shown in the inset in Fig. 1 and are dis- viding a larger stabilizing influence. Esym also partially cussed in Refs. [18, 19]. Quantum Monte Carlo (QMC) determines the equilibrium value of Ye, which can alter results are also shown in the inset in Fig. 1. The linear the value of γYL , but this is a smaller effect. In our nu- ′ correlation between Esym and Esym in the QMC results merical PNS simulations this effect is accounted for. In 3

6 0.1 0.1 outer regions to a state close to neutral buoyancy. GM3 5 has a slightly steeper entropy gradient, in part due to 4 ′ 1 1 its larger Esym, than IU-FSU. This partially accounts 3

Entropy for the slightly larger neutrino luminosity at early times 2 5 5 for GM3. As time progresses, convection steadily digs 10 1 15 10 15 30 30 deeper into the core of the PNS. For both EoSs, con- 0.4

0.35 vection proceeds all the way to the core by 15 seconds 0.1 0.1 0.3 1 1 into the simulation, but it lasts in the interior regions 0.25 5 for a much longer period of time for IU-FSU resulting L 5

Y 0.2 10 15 10 0.15 30 in more rapid lepton depletion in the core. The exact 15 0.1 30 details of how convection proceeds depend on the initial 0.05 L 0 conditions of the PNS and the behavior of γs and γY 0 0.5 1 1.5 0 0.5 1 1.5 Baryon Number (solar units) Baryon Number (Solar Units) for a given EoS. For the conditions encountered in the

PNS, the variation in γYL between the two EoSs we em- FIG. 2: Evolution of the entropy (top panels) and lepton frac- ploy is significantly larger than the variation of γs. More . M⊙ tion (bottom panels) in a 1 6 rest mass PNS for the GM3 important to the neutrino signal accompanying PNS for- EoS (left panels) and the IU-FSU EoS (right panels). The ∼ grayed regions are convectively unstable. The labels corre- mation, in GM3 convection ceases in the mantle by 5 spond to the model times in seconds. seconds, whereas convection in the mantle proceeds until ∼ 12 seconds in IU-FSU. This difference is mainly driven ′ by the difference in Esym between the two EoSs. As the contrast, the properties of zero temperature nuclear mat- mantle contracts, the second term in Eq.(2) becomes inter should have significantly less effect on the behavior creasingly dominant and is eventually able to stabilize ′ of γs. We find the variation of γs to be significantly less convection. Qualitatively, increasing Esym will shut-off than that of γYL in the two EoSs considered here, which convection at an earlier time. is consistent with expectations. The depth to which convection penetrates in the core Medium modifications of the neutrino interaction rates and how long convection proceeds in the core are depen- also influence PNS evolution. These effects have been dent upon the opacities as well as the EoS. When only investigated in earlier work where neutrino scattering mean field effects on the opacities are considered (i.e. and absorption rates on nucleons and leptons were calcu- when the neutrino mean free path is shorter), convection lated within the relativistic random phase approximation does not proceed all the way to the center of the PNS in (RPA) [25]. Here, the effects due to strong and electro- the GM3 models. When RPA effects are included, con- magnetic correlations, degeneracy, and relativistic cur- vection does proceed to the central regions of the core. rents are included in a similar way. Further, the effective An increased diffusion rate allows the core to heat up and interactions we employ in the RPA are obtained from the deleptonize more rapidly, thereby decreasing the stabi- underlying model for the EoS to ensure that the corre- lizing lepton gradients and increasing the de-stabilizing lation functions which determine the neutrino scattering entropy gradients. kernels satisfy the generalized compressibility sum-rules Of course, variations in the convective evolution of the consistent with the EoS. Because the axial portion of PNS are only interesting to the extent they are poten- the response accounts for the majority of the neutrino tially observable in the neutrino emission from a nearby scattering cross-section, the presence of an axial inter- supernova. In Fig. 3 the expected neutrino count rates action can strongly affect the total neutrino mean free for a detector similar to Super Kamiokande-III are shown path. To take this into account, we introduce an effec- for a number of PNS cooling models. We have assumed tive short-range interaction in the spin channel through a threshold energy of 7.5 MeV, a detector mass of 50 kt, the Migdal parameter, g′ [26]. The strength of this in- a detector efficiency above threshold of unity [28], and a teraction is tuned to reproduce the spin-susceptibility of distance of 10 kpc to the supernova. Equipartition has neutron matter obtained from microscopic calculations been assumed between neutrino flavors and the spectral [27]. For densities above nuclear saturation, the RPA temperatures are calculated from the average radius of causes a significant enhancement of the mean free path neutrino decoupling in the PNS. The neutrino luminosi- relative to the mean field approximation due to the re- ties and average energies for two models are included as pulsive nature of the nuclear interaction at high density. supplemental data. We now consider the evolution of the internal structure Both the GM3 and IU-FSU EoSs show enhanced lumi- of the PNS with convection and varying prescriptions for nosities at early times relative to the models not includ- the opacities. In Fig. 2, the evolutions of the entropy ing convection. There is only a small difference between and lepton fraction for the two equations of state are the two equations of state at low (sub-nuclear) density, shown. Over the first second in both models, convec- so differences prior to one second are small. The neu- tion smoothes the entropy and lepton gradients in the trino count rate is increased by about 30% relative to the 4

No Convection g’=0.6 GM3 neutrino signal is signiﬁcantly reduced and drops below Convection MF GM3 Convection g’=0.6 GM3 the detectable threshold at an earlier time. Convection g’=0.6 IU-FSU

3 10 In summary, using a self-consistent model for the PNS core physics, we ﬁnd that the late time neutrino signal ) )

∞ from a core collapse supernova is likely to contain a di- −1

0.35 rect diagnostic of the nuclear symmetry energy at high density. With current neutrino detectors, these eﬀects 2 0.3 10

Count Rate (s should be readily discernible in the neutrino light curve 0.25 of a single nearby supernova. 0.2 Support from the UCOP (09-IR-07-117968-WOOS),

Counts (3 s −> 10 s)/ (0.1 the NNSA/DOE SSGF DE-FC52-08NA28752, the 1 0.3 0.35 0.4 0.45 10 Counts (0.1 s −> 1 s)/ Counts (0.1 s −> ∞ ) US NSF grant AST-0909129, DOE grant DE-AC52-

0 1 10 10 06NA25396 (LANL), DOE TC on “Neutrinos and Nu- Time (s) cleosynthesis”, and the Spanish MEC grant AYA 2010- FIG. 3: Count rates as a function of time for a number of 21097-C03-02 is gratefully acknowledged. 1.6M⊙ PNS models with and without convection. The black line is for neutrino opacities calculated in the mean field approximation, while all the other lines are for models that use RPA opacities with g′ = 0.6. The inset plot shows the integrated number of counts from 0.1 s to 1 s divided by the total [1] A. Burrows and J. M. Lattimer, number of counts for t > 0.1 second on the horizontal axis, Astrophys. J. 307, 178 (1986). and the number of counts for t> 3 seconds on divided by the [2] J. R. Wilson and R. W. Mayle, total number of counts for t > 0.1 second. The stars corre- Phys. Rep. 163, 63 (1988). spond to the IU-FSU EoS and the circles to the GM3 EoS. [3] W. Keil and H.-T. Janka, Astron. & Astrophys. 296, 145 Symbol sizes correspond to various neutron star rest masses (1995). 513 ranging from 1.2M⊙ to 2.1M⊙. Colors correspond to different [4] J. A. Pons et al., Astrophys. J. , 780 (1999). values of the Migdal parameter, g′. [5] T. Fischer et al., Astron. & Astrophys. 517, A80+ (2010). [6] L. Hüdepohl et al., Physi. Rev. Lett. 104, 251101 (2010). [7] W. Keil, H.-T. Janka, and G. Raffelt, Phys. Rev. D 51, models that do not include convection. This is reason- 6635 (1995). [8] J. A. Pons et al., Astrophys. J. 553, 382 (2001). ably consistent with the early time enhancement seen in [9] J. A. Pons et al., Phys. Rev. Lett. 86, 5223 (2001). multi-dimensional models [15]. After a second, the count [10] R. I. Epstein, Mon. Not. Roy. Astro. Soc. 188, 305 rates between the two EoSs begin to diverge. The most (1979). obvious feature in the count rate for GM3 appears at ∼ 3 [11] A. Burrows, Astrophys. J. Lett. 318, L57 (1987). seconds, which is coincident with the end of convection [12] W. Keil, H. Janka, and E. Mueller, Astrophys. J. Lett. in the mantle. For the IU-FSU EoS, the break is also 473, L111+ (1996). 493 at the time at which mantle convection ends (∼ 10 sec- [13] A. Mezzacappa et al., Astrophys. J. , 848 (1998). et al. 645 onds), although it is hard to distinguish from the point [14] L. Dessart , Astrophys. J. , 534 (2006). [15] R. Buras et al., Astron. & Astrophys. 457, 281 (2006). at which the PNS becomes optically thin. As was ar- [16] S. W. Bruenn and T. Dineva, gued previously, the position of this break reflects the Astrophys. J. Lett. 458, L71+ (1996). density dependence of the nuclear symmetry energy at [17] J. A. Miralles, J. A. Pons, and V. A. Urpin, 543 nB > n0 and therefore provides a direct observable of Astrophys. J. , 1001 (2000). the properties of nuclear matter in the PNS neutrino sig- [18] M. B. Tsang et al., Prog. Part. Nucl. Phys. 66, 400 (2011). 82 nal. Although core convection does not seem to affect the [19] F. J. Fattoyev et al., Phys. Rev. C , 055803 (2010). [20] S. Gandolfi, J. Carlson, and S. Reddy, ArXiv e-prints break, it may impact the subsequent cooling timescale. (2011), arXiv:1101.1921 . In the inset in Fig. 3, integrated neutrino counts over [21] N. K. Glendenning and S. A. Moszkowski, two time windows are shown for a number of PNS masses. Phys. Rev. Lett. 67, 2414 (1991). There is a clear separation between the two EoSs inde- [22] A. W. Steiner et al., Phys. Rep. 411, 325 (2005). pendent of mass. The time of the convective break cre- [23] J. M. Lattimer and M. Prakash, ates this separation. This illustrates that this diagnostic Astrophys. J. 550, 426 (2001). 64 of the symmetry energy does not require an accurate de- [24] C. J. Horowitz and J. Piekarewicz, Phys. Rev. C , termination of the PNS mass. 062802 (2001). [25] S. Reddy et al., Phys. Rev. C 59, 2888 (1999). The inclusion of nucleon correlations through the RPA [26] H. Kim, J. Piekarewicz, and C. J. Horowitz, begins to significantly affect the neutrino emission after Phys. Rev. C 51, 2739 (1995). about three seconds. Initially, the luminosities are in- [27] S. Fantoni, A. Sarsa, and K. E. Schmidt, Phys. Rev. creased as energy and lepton number are able to more Lett. 87, 181101 (2001). rapidly diffuse out of the core, but at later times the [28] M. Ikeda et al., Astrophys. J. 669, 519 (2007).