Arxiv:1112.0335V1 [Astro-Ph.HE] 1 Dec 2011 Ini Nitnesuc Fnurnso L Aos It flavors
Total Page:16
File Type:pdf, Size:1020Kb
Proto-Neutron Star Cooling with Convection: The Effect of the Symmetry Energy L.F. Roberts,1 G. Shen,2 V. Cirigliano,2 J.A. Pons,3 S. Reddy,2, 4 and S.E. Woosley1 1Department of Astronomy and Astrophysics, University of California, Santa Cruz, CA 95064 USA 2Theoretical Division, Los Alamos National Laboratory, Los Alamos, NM 87544 3Departament de F´ısica Aplicada, Universitat d’Alacant, Ap. Correus 99, 03080 Alacant, Spain 4Institute for Nuclear Theory, University of Washington, Seattle, Washington 98195 We model neutrino emission from a newly born neutron star subsequent to a supernova explosion to study its sensitivity to the equation of state, neutrino opacities, and convective instabilities at high baryon density. We find the time period and spatial extent over which convection operates is sensitive to the behavior of the nuclear symmetry energy at and above nuclear density. When convection ends within the proto-neutron star, there is a break in the predicted neutrino emission that may be clearly observable. PACS numbers: 26.50.+x, 26.60.-c, 21.65.Mn, 95.85.Ry The hot and dense proto-neutron star (PNS) born sub- formation. It has long been recognized that the outer sequent to core-collapse in a type II supernova explo- PNS mantle is unstable to convection soon after the pas- sion is an intense source of neutrinos of all flavors. It sage of the supernova shock, due to negative entropy gra- emits the 3 − 5 × 1053 ergs of gravitational binding en- dients [10]. This early period of instability beneath the ergy gained during collapse as neutrino radiation on a neutrino spheres has been studied extensively in both one time scale of tens of seconds as it contracts, becomes and two dimensions, with the hope that it could increase increasingly neutron-rich and cools. Cooling of the PNS the neutrino luminosities enough to lead to a successful and the concomitant neutrino emission are driven by neu- explosion [2, 11–15]. Although the role of convection at trino diffusion and convection along the lepton number late times was studied in Refs. [2, 11], this work is the first and entropy gradients left behind within the PNS after attempt at exploring its connection to the underlying mi- core-bounce, where the matter density and temperature crophysics and the interplay between modified opacities are in the range ρ =2 − 6 × 1014 g/cm3 and T =5 − 40 and convection in shaping the observable neutrino signal. MeV, respectively [1–6]. While the supernova explosion We model convection as a diffusive process described mechanism and associated fall back of material are ex- by standard time dependent mixing length theory [2]. In pected to influence the neutrino emission at early time Ref. [15] it was shown that hydrodynamic simulations in (i.e. t . 1 s post bounce) the late time neutrino signal is two-dimensions were well reproduced by a simple mixing shaped by the properties of the PNS, such as the nuclear scheme in one-dimension, suggesting that this is a rea- equation of state (EoS), neutrino opacities in dense mat- sonable approximation. Further, since convection pro- ter, and other microphysical properties that affect the ceeds efficiently throughout our simulations, our results cooling timescale by influencing either neutrino diffusion are not particularly sensitive to the chosen mixing length. or convection [4, 7–9]. The linear growth rates are obtained from the standard Here, we present one-dimensional hydrostatic models Ledoux stability analysis as of PNS evolution out to late times for two EoSs. We 2 − g ∇ ∇ include approximate convective transport along with dif- ω = (γs ln(s)+ γYL ln(YL)) , (1) γnB fusive neutrino transport, both consistent with the un- derlying EoS. This allows us to gauge the importance of where arXiv:1112.0335v1 [astro-ph.HE] 1 Dec 2011 convection and the effect of medium modifications to the ∂ ln P ∂ ln P ∂ ln P neutrino opacities in dense matter to the temporal char- γnB = γs = γYL = ∂ ln n ∂ ln s ∂ ln Y acteristics of the neutrino signal. The basic framework B s,YL nB ,YL L nB ,s for PNS evolution is similar to that described in [4], ex- and g, P , s, nB, and YL are the local acceleration due to cept that a treatment of convection is included. We find gravity, the pressure, entropy per baryon, baryon num- that the behavior of the nuclear symmetry energy (which ber density and the fraction of leptons in dense matter, is a measure of the energy difference between dense neu- respectively. Convective instability sets in for ω2 > 0. tron matter and symmetric nuclear matter) at high den- The form of the growth rate clearly implies that there sity significantly influences convection and thereby affects will be a strong interplay between the nuclear EoS and the observable neutrino signal. the patterns of convection within the PNS. Large scale convective overturn of material will directly In addition, the PNS may be subject to doubly diffu- transport energy and lepton number in the PNS and al- sive instabilities due to the lateral transport of compo- ters the gradients along which neutrinos diffuse, thereby sition and energy by neutrinos [2, 16, 17]. Early one- strongly affecting the neutrino signal accompanying PNS dimensional work [2] suggested that neutron fingering 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 be- tween 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.