-
Muon-induced collisional flavor instability in core-collapse supernova
Authors:
Jiabao Liu,
Hiroki Nagakura,
Ryuichiro Akaho,
Akira Ito,
Masamichi Zaizen,
Shun Furusawa,
Shoichi Yamada
Abstract:
Neutrinos are known to undergo flavor conversion among their three flavors. In the theoretical modeling of core-collapse supernova (CCSN), there has been a great deal of attention to recent discoveries of a new type of neutrino flavor conversions, namely collisional flavor instability (CFI), in which the instability is induced by the flavor-dependent decoherence due to the disparity of neutrino-ma…
▽ More
Neutrinos are known to undergo flavor conversion among their three flavors. In the theoretical modeling of core-collapse supernova (CCSN), there has been a great deal of attention to recent discoveries of a new type of neutrino flavor conversions, namely collisional flavor instability (CFI), in which the instability is induced by the flavor-dependent decoherence due to the disparity of neutrino-matter interactions among flavors. In this paper, we study how the appearance of on-shell muons and associated neutrino-matter interactions can impact CFIs based on linear stability analysis of flavor conversions. Some striking results emerge from the present study. First, we analytically show that breaking beta- and pair equilibrium is a necessary condition to trigger CFIs. This also indicates that CFIs with on-shell muons could appear in $e τ$ and $μτ$ neutrino mixing sectors in very high-density region ($\gtrsim 10^{13} {\rm g/cm^{3}}$), exhibiting a possibility of large impacts of CFIs on CCSN. Second, resonance-like CFIs, having a much higher growth rate than normal CFIs, can be triggered by muons. The resonance point of CFIs is different between $e τ$ and $μτ$ sectors; the former (latter) occurs at $μ_{e (μ)} = μ_{n} - μ_{p}$, where $μ_{i}$ denotes the chemical potential of $i$ constitute ($n$ and $p$ represent neutrons and protons, respectively). Our result suggests that the non-linear evolution of CFI with on-shell muons would induce flavor conversions with the complex interplay among all three different neutrino-mixing sectors.
△ Less
Submitted 15 July, 2024;
originally announced July 2024.
-
Leptonic and semi-leptonic neutrino interactions with muons in the proto-neutron star cooling
Authors:
Ken'ichi Sugiura,
Shun Furusawa,
Kohsuke Sumiyoshi,
Shoichi Yamada
Abstract:
It is known that muons are scarce just after the birth of a proto-neutron star via a supernova explosion but get more abundant as the proto-neutron star cools via neutrino emissions on the Kelvin-Helmholtz timescale.
In this paper we evaluate all the relevant rates of the neutrino interactions with muons at different times in the proto-neutron star cooling.
We are particularly interested in th…
▽ More
It is known that muons are scarce just after the birth of a proto-neutron star via a supernova explosion but get more abundant as the proto-neutron star cools via neutrino emissions on the Kelvin-Helmholtz timescale.
In this paper we evaluate all the relevant rates of the neutrino interactions with muons at different times in the proto-neutron star cooling.
We are particularly interested in the late phase ($ t \gtrsim 10 \operatorname{s}$), which will be accessible in the next Galactic supernova but has not been studied well so far.
We calculate both leptonic and semi-leptonic processes, for the latter of which we pay attention also to the form factors with their dependence on the transferred momentum as well as to the modification of the dispersion relations for nucleons on the mean field level.
We find that the flavor-exchange reactions $ν_e + μ^- \rightarrow ν_μ + e^-$ and $\barν_μ + μ^- \rightarrow \barν_e + e^-$ can be dominant, particularly at low energies, over the capture of $ν_e$ on neutron and the scatterings of $\barν_μ$ on nucleons as the opacity sources for these species and that the inverse muon decay $ \barν_e + ν_μ + e^- \leftrightarrows μ^- $ can overwhelm the scatterings of $\barν_e$ and $ν_μ$ on nucleons again at low energies.
At high energies, on the other hand, the corrections in the semi-leptonic processes mentioned above are more important.
We also show the non-trivial energy- and angular dependences of the flavor-exchange reactions and the inverse muon decay.
In the study of the diffusion coefficients from these reactions, we find that $\barν_μ$ is most affected.
These pieces of information are indispensable for numerical computations and the interpretation of results thereof for the proto-neutron star cooling particularly at the very late phase.
△ Less
Submitted 13 November, 2022;
originally announced November 2022.
-
Nuclei in Core-Collapse Supernovae Engine
Authors:
Shun Furusawa,
Hiroki Nagakura
Abstract:
Herein, we review the nuclear equations of state (EOSs) %for core-collapse supernova simulations and the constituent nuclei of core-collapse supernovae (CCSNe) and their roles in CCSN simulations. Various nuclei such as deuterons, iron, and extremely neutron-rich nuclei compose in the central engines of CCSNe. The center of a collapsing core is dominated by neutron-rich heavy nuclei prior to the o…
▽ More
Herein, we review the nuclear equations of state (EOSs) %for core-collapse supernova simulations and the constituent nuclei of core-collapse supernovae (CCSNe) and their roles in CCSN simulations. Various nuclei such as deuterons, iron, and extremely neutron-rich nuclei compose in the central engines of CCSNe. The center of a collapsing core is dominated by neutron-rich heavy nuclei prior to the occurrence of core bounce. Their weak interactions significantly affect the neutrino emission and the size of the produced proto-neutron star. After a core bounce, heavy nuclei are dissolved to protons, neutrons, and light nuclei between the expanding shock wave and the newly formed neutron star (NS). Some of the key components in determining the shock-wave dynamics and supernova explosion of outer envelopes are neutrino interactions of nucleons and light nuclei such as deuterons. An EOS provides the relations between thermodynamical properties and the nuclear composition, and is needed to simulate this explosion. Further investigations on uniform and non-uniform nuclear matter are needed to improve the understanding of the mechanism of CCSNe and the properties of supernova nuclei. The knowledge of the EOS for uniform nuclear matter is being continually improved by a combination of microscopic calculations, terrestrial experiments, and NS observations. With reference to various nuclear experiments and current theories, the finite temperature effects on heavy nuclei, formation of light nuclei in dilute nuclear matter, and transition to uniform nuclear matter should be improved in the model of the EOS for non-uniform nuclear matter.
△ Less
Submitted 19 December, 2022; v1 submitted 2 November, 2022;
originally announced November 2022.
-
Effects of nuclear matter and composition in core-collapse supernovae and long-term proto-neutron star cooling
Authors:
Kohsuke Sumiyoshi,
Shun Furusawa,
Hiroki Nagakura,
Akira Harada,
Hajime Togashi,
Ken'ichiro Nakazato,
Hideyuki Suzuki
Abstract:
We study the influence of hot and dense matter in core-collapse supernovae by adopting up-to-date nuclear equation of state (EOS) based on the microscopic nuclear many-body frameworks. We explore effects of EOS based on the Dirac Brueckner Hartree-Fock theory through comparisons with those based on the variational method. We also examine effects of the differences in the composition of nuclei and…
▽ More
We study the influence of hot and dense matter in core-collapse supernovae by adopting up-to-date nuclear equation of state (EOS) based on the microscopic nuclear many-body frameworks. We explore effects of EOS based on the Dirac Brueckner Hartree-Fock theory through comparisons with those based on the variational method. We also examine effects of the differences in the composition of nuclei and nucleons by using the same EOS by the variational method but employing two different treatments in computations of nuclear abundances. We perform numerical simulations of core-collapse supernovae adopting the three EOSs. We also perform numerical simulations of the long-term evolution over 70 s of the proto-neutron star cooling. We show that impacts by different modeling of composition are remarkable as in those by different treatments of uniform matter in the gravitational collapse, bounce, and shock propagation. The cooling of proto-neutron star and the resulting neutrino emission are also affected by the compositional difference even if the same treatment in computing uniform matter of EOS.
△ Less
Submitted 9 December, 2022; v1 submitted 6 September, 2022;
originally announced September 2022.
-
Equation of state in neutron stars and supernovae
Authors:
Kohsuke Sumiyoshi,
Toru Kojo,
Shun Furusawa
Abstract:
Neutron stars and supernovae provide cosmic laboratories of highly compressed matter at supra nuclear saturation density which is beyond the reach of terrestrial experiments. The properties of dense matter is extracted by combining the knowledge of nuclear experiments and astrophysical observations via theoretical frameworks. A matter in neutron stars is neutron rich, and may further accommodate n…
▽ More
Neutron stars and supernovae provide cosmic laboratories of highly compressed matter at supra nuclear saturation density which is beyond the reach of terrestrial experiments. The properties of dense matter is extracted by combining the knowledge of nuclear experiments and astrophysical observations via theoretical frameworks. A matter in neutron stars is neutron rich, and may further accommodate non-nucleonic degrees of freedom such as hyperons and quarks. The structure and composition of neutron stars are determined by equations of state of matter, which are the primary subject in this chapter. In case of supernovae, the time evolution includes several dynamical stages whose descriptions require equations of state at finite temperature and various lepton fractions. Equations of state also play essential roles in neutron star mergers which allow us to explore new conditions of matter not achievable in static neutron stars and supernovae. Several types of hadron-to-quark transitions, from first order transitions to crossover, are reviewed, and their characteristics are summarized.
△ Less
Submitted 30 June, 2022;
originally announced July 2022.
-
Protoneutron Star Convection Simulated with a New General Relativistic Boltzmann Neutrino Radiation-Hydrodynamics Code
Authors:
Ryuichiro Akaho,
Akira Harada,
Hiroki Nagakura,
Wakana Iwakami,
Hirotada Okawa,
Shun Furusawa,
Hideo Matsufuru,
Kohsuke Sumiyoshi,
Shoichi Yamada
Abstract:
We investigate the protoneutron star (PNS) convection using our newly developed general relativistic Boltzmann neutrino radiation-hydrodynamics code. This is a pilot study for more comprehensive investigations later. As such, we take a snapshot of a PNS at 2.3 seconds after bounce from a 1D PNS cooling calculation and run our simulation for $\sim160\,\mathrm{ms}$ in 2D under axisymmetry. The origi…
▽ More
We investigate the protoneutron star (PNS) convection using our newly developed general relativistic Boltzmann neutrino radiation-hydrodynamics code. This is a pilot study for more comprehensive investigations later. As such, we take a snapshot of a PNS at 2.3 seconds after bounce from a 1D PNS cooling calculation and run our simulation for $\sim160\,\mathrm{ms}$ in 2D under axisymmetry. The original PNS cooling calculation neglected convection entirely and the initial condition is linearly unstable to convection. We find in our 2D simulation that convection is instigated there indeed and expands inward after being full-fledged. The convection is then settled to a quasi-steady state in $\sim100\,\mathrm{ms}$, being sustained by the negative $Y_e$ gradient, which is in turn maintained by neutrino emissions. It enhances the luminosities and mean energies of all species of neutrinos compared to 1D. Taking advantage of the Boltzmann solver, we analyze the possible occurrence of the neutrino fast flavor conversion (FFC). We found that FFC is likely to occur in the regions, where $Y_e$ is lower, and that the growth rate can be as high as $\sim 10^{-1}\,{\mathrm{cm}^{-1}}$.
△ Less
Submitted 20 December, 2022; v1 submitted 3 June, 2022;
originally announced June 2022.
-
Principal-Axis Analysis of the Eddington Tensor for the Early Post-Bounce Phase of Rotational Core-Collapse Supernovae
Authors:
Wakana Iwakami,
Akira Harada,
Hiroki Nagakura,
Ryuichiro Akaho,
Hirotada Okawa,
Shun Furusawa,
Hideo Matsufuru,
Kohsuke Sumiyoshi,
Shoichi Yamada
Abstract:
Using full Boltzmann neutrino transport, we performed two-dimensional (2D) core-collapse supernova simulations in axisymmetry for two progenitor models with 11.2M and 15.0M both rotational and non-rotational. We employed the results obtained in the early post-bounce phase (t < 20 ms) to assess performance under rapid rotation of some closure relations commonly employed in the truncated moment meth…
▽ More
Using full Boltzmann neutrino transport, we performed two-dimensional (2D) core-collapse supernova simulations in axisymmetry for two progenitor models with 11.2M and 15.0M both rotational and non-rotational. We employed the results obtained in the early post-bounce phase (t < 20 ms) to assess performance under rapid rotation of some closure relations commonly employed in the truncated moment method. We first made a comparison in 1D under spherical symmetry, though, of the Eddington factor p defined in the fluid rest frame (FR). We confirmed that the maximum entropy closure for the Fermionic distribution (MEFD) performs better than others near the proto-neutron star surface, where p < 1/3 occurs, but does not work well even in 1D when the phase space occupancy satisfies e < 0.5 together with p < 1/3, the condition known to be not represented by MEFD. For the 2D models with the rapid rotation, we employed the principal-axis analysis of the Eddington tensor. We paid particular attention to the direction of the longest principal axis. We observed in FR that it is aligned neither with the radial direction nor with the neutrino flux in 2D, particularly so in convective and/or rapidly rotating regions, the fact not accommodated in the moment method. We repeated the same analysis in the laboratory frame (LB) and found again that the direction of the longest principal axis is not well reproduced by MEFD because the interpolation between the optically thick and thin limits is not very accurate in this frame.
△ Less
Submitted 21 May, 2022; v1 submitted 13 September, 2021;
originally announced September 2021.
-
Multidimensional Boltzmann Neutrino Transport Code in Full General Relativity for Core-collapse Simulations
Authors:
Ryuichiro Akaho,
Akira Harada,
Hiroki Nagakura,
Kohsuke Sumiyoshi,
Wakana Iwakami,
Hirotada Okawa,
Shun Furusawa,
Hideo Matsufuru,
Shoichi Yamada
Abstract:
We develop a neutrino transfer code for core-collapse simulations, that directly solves the multidimensional Boltzmann equations in full general relativity. We employ the discrete ordinate method, which discretizes the six dimensional phase space. The code is an extension of our special relativistic code coupled to a Newtonian hydrodynamics code, which is currently employed for core-collapse super…
▽ More
We develop a neutrino transfer code for core-collapse simulations, that directly solves the multidimensional Boltzmann equations in full general relativity. We employ the discrete ordinate method, which discretizes the six dimensional phase space. The code is an extension of our special relativistic code coupled to a Newtonian hydrodynamics code, which is currently employed for core-collapse supernova simulations. In order to demonstrate our code's capability to treat general relativistic effects, we conduct some tests: we first compute the free streaming of neutrinos in the Schwarzschild and Kerr spacetimes and compare the results with the geodesic curves; in the Schwarzschild case we deploy not only a 1-dimensional grid in space under spherical symmetry but also a 2-dimensional spatial mesh under axisymmetry in order to assess the capability of the code to compute the spatial advection of neutrinos; secondly, we calculate the neutrino transport in a fixed matter background, which is taken from a core-collapse supernova simulation with our general relativistic but spherically symmetric Boltzmann-hydrodynamics code, to obtain a steady neutrino distribution; the results are compared with those given by the latter code.
△ Less
Submitted 19 March, 2021; v1 submitted 21 October, 2020;
originally announced October 2020.
-
Simulations of the Early Post-Bounce Phase of Core-Collapse Supernovae in Three-Dimensional Space with Full Boltzmann Neutrino Transport
Authors:
Wakana Iwakami,
Hirotada Okawa,
Hiroki Nagakura,
Akira Harada,
Shun Furusawa,
Kosuke Sumiyoshi,
Hideo Matsufuru,
Shoichi Yamada
Abstract:
We report on the core-collapse supernova simulation we conducted for a 11.2 M progenitor model in three-dimensional space up to 20 ms after bounce, using a radiation hydrodynamics code with full Boltzmann neutrino transport. We solve the six-dimensional Boltzmann equations for three neutrino species and the three-dimensional compressible Euler equations with Furusawa and Togashi's nuclear equation…
▽ More
We report on the core-collapse supernova simulation we conducted for a 11.2 M progenitor model in three-dimensional space up to 20 ms after bounce, using a radiation hydrodynamics code with full Boltzmann neutrino transport. We solve the six-dimensional Boltzmann equations for three neutrino species and the three-dimensional compressible Euler equations with Furusawa and Togashi's nuclear equation of state. We focus on the prompt convection at 10 ms after bounce and investigate how neutrinos are transported in the convective matter. We apply a new analysis based on the eigenvalues and eigenvectors of the Eddington tensor and make a comparison between the Boltzmann transport results and the M1 closure approximation in the transition regime between the optically thick and thin limits. We visualize the eigenvalues and eigenvectors using an ellipsoid, in which each principal axis is parallel to one of the eigenvectors and has a length proportional to the corresponding eigenvalue. This approach enables us to understand the difference between the Eddington tensor derived directly from the Boltzmann simulation and the one given by the M1 prescription from a new perspective. We find that the longest principal axis of the ellipsoid is almost always nearly parallel to the energy flux in the M1 closure approximation whereas in the Boltzmann simulation it becomes perpendicular in some transition regions, where the mean free path is 0.1 times the radius. In three spatial dimensions, the convective motions make it difficult to predict where this happens and to possibly improve the closure relation there.
△ Less
Submitted 31 August, 2020; v1 submitted 5 April, 2020;
originally announced April 2020.
-
The Boltzmann-radiation-hydrodynamics Simulations of Core-collapse Supernovae with Different Equations of State: the Role of Nuclear Composition and the Behavior of Neutrinos
Authors:
Akira Harada,
Hiroki Nagakura,
Wakana Iwakami,
Hirotada Okawa,
Shun Furusawa,
Kohsuke Sumiyoshi,
Hideo Matsufuru,
Shoichi Yamada
Abstract:
Using the Boltzmann-radiation-hydrodynamics code, which solves the Boltzmann equation for neutrino transport, we present the results of the simulations with the nuclear equations of state (EOSs) of Lattimer and Swesty (LS) and Furusawa and Shen (FS). We extend the simulation time of the LS model and conduct thorough investigations, though our previous paper briefly reported some of the results. On…
▽ More
Using the Boltzmann-radiation-hydrodynamics code, which solves the Boltzmann equation for neutrino transport, we present the results of the simulations with the nuclear equations of state (EOSs) of Lattimer and Swesty (LS) and Furusawa and Shen (FS). We extend the simulation time of the LS model and conduct thorough investigations, though our previous paper briefly reported some of the results. Only the LS model shows the shock revival. This seems to originate from the nuclear composition: the different nuclear composition results in the different energy loss by photodissociation and hence the different strength of the prompt convection and the later neutrino-driven convection. The protoneutron star seen in the FS model is more compact than that in the LS model because the existence of multinuclear species softens the EOS. For the behavior of neutrinos, we examined the flux and the Eddington tensor of neutrinos. In the optically thick region, the diffusion of neutrinos and the dragging by the motion of matter determine the flux. In the optically thin region, the free-streaming determines it. The Eddington tensor is compared with that obtained from the M1-closure relation. The M1-closure scheme overestimates the contribution from the velocity-dependent terms in the semitransparent region.
△ Less
Submitted 17 November, 2020; v1 submitted 19 March, 2020;
originally announced March 2020.
-
Quantum statistical effects in warm nuclear matter with light and heavy clusters
Authors:
Shun Furusawa,
Igor Mishustin
Abstract:
We have investigated the compositions of hot and dense nuclear matter with the focus on the quantum-statistical effects for light clusters. Our main observation is that the formation of heavy nuclei in stellar matter leads to the reduction of the number densities of light clusters. As a result, the quantum statistical effects such as Bose-Einstein condensation of deuterons and $α$-particles are su…
▽ More
We have investigated the compositions of hot and dense nuclear matter with the focus on the quantum-statistical effects for light clusters. Our main observation is that the formation of heavy nuclei in stellar matter leads to the reduction of the number densities of light clusters. As a result, the quantum statistical effects such as Bose-Einstein condensation of deuterons and $α$-particles are suppressed. The deviations in number densities of light clusters between Boltzmann and quantum statistics at sub-nuclear densities and temperatures 1-3 MeV are $\sim$0.2$\%$ at most. The condensation of $α$-particles in iso-symmetric nuclear matter is predicted only under the assumptions that it is composed only of nucleons and light clusters and no heavy nuclei are present. We also found that the Coulomb screening effects hardly affect the critical baryon densities for alpha condensation, although mass and chemical potential of $α$-particles are modified.
△ Less
Submitted 19 December, 2019;
originally announced December 2019.
-
Fast collective neutrino oscillations inside the neutrino sphere in core-collapse supernovae
Authors:
Milad Delfan Azari,
Shoichi Yamada,
Taiki Morinaga,
Hiroki Nagakura,
Shun Furusawa,
Akira Harada,
Hirotada Okawa,
Wakana Iwakami,
Kohsuke Sumiyoshi
Abstract:
Neutrinos are believed to have a key role in the explosion mechanism of core-collapse supernovae as they carry most of the energy released by the gravitational collapse of a massive star. If their flavor is converted fast inside the neutrino sphere, the supernova explosion may be influenced. This paper is reporting the results of the extended work of our previous paper. We perform a thorough surve…
▽ More
Neutrinos are believed to have a key role in the explosion mechanism of core-collapse supernovae as they carry most of the energy released by the gravitational collapse of a massive star. If their flavor is converted fast inside the neutrino sphere, the supernova explosion may be influenced. This paper is reporting the results of the extended work of our previous paper. We perform a thorough survey of the ELN crossing in one of our self-consistent, realistic Boltzmann simulations in two spatial dimensions under axisymmetry for the existence of the crossings between $ν_e$ and $\barν_e$ angular distributions, or the electron lepton number (ELN) crossing. We report for the first time the positive detections deep inside the core of the massive star in the vicinity of neutrino sphere at $r$ $\approx$ 16 - 21 km. We find that low values of the electron fraction $Y_e$ produced by convective motions together with the appearance of light elements are critically important to give rise to the ELN crossing by enhancing the chemical potential difference between proton and neutron, and hence by mitigating the Fermi-degeneracy of $ν_e$. Since the region of positive detection are sustained and, in fact, expanding with time, it may have an impact on the explosion of core-collapse supernovae, observational neutrino astronomy and nucleosynthesis of heavy nuclei.
△ Less
Submitted 15 January, 2020; v1 submitted 14 October, 2019;
originally announced October 2019.
-
New Neutron Star Equation of State with Quark-Hadron Crossover
Authors:
Gordon Baym,
Shun Furusawa,
Tetsuo Hatsuda,
Toru Kojo,
Hajime Togashi
Abstract:
We present a much improved equation of state for neutron star matter, QHC19, with a smooth crossover from the hadronic regime at lower densities to the quark regime at higher densities. We now use the Togashi et al.~equation of state (Togashi:2017), a generalization of the Akmal-Pandharipande-Ravenhall equation of state of uniform nuclear matter, in the entire hadronic regime; the Togashi equation…
▽ More
We present a much improved equation of state for neutron star matter, QHC19, with a smooth crossover from the hadronic regime at lower densities to the quark regime at higher densities. We now use the Togashi et al.~equation of state (Togashi:2017), a generalization of the Akmal-Pandharipande-Ravenhall equation of state of uniform nuclear matter, in the entire hadronic regime; the Togashi equation of state consistently describes non-uniform as well as uniform matter, and matter at beta equilibrium without the need for an interpolation between pure neutron and symmetric nuclear matter. We describe the quark matter regime at higher densities with the Nambu--Jona--Lasinio model, now identifying tight constraints on the phenomenological universal vector repulsion between quarks and the pairing interaction between quarks arising from the requirements of thermodynamic stability and causal propagation of sound. The resultant neutron star properties agree very well with the inferences of the LIGO/Virgo collaboration, from GW170817, of the pressure vs. baryon density, neutron star radii, and tidal deformabilities. The maximum neutron star mass allowed by QHC19 is 2.35 $M_\odot$, consistent with all neutron star mass determinations.
△ Less
Submitted 22 November, 2019; v1 submitted 21 March, 2019;
originally announced March 2019.
-
Comparing treatments of weak reactions with nuclei in simulations of core-collapse supernovae
Authors:
Hiroki Nagakura,
Shun Furusawa,
Hajime Togashi,
Sherwood Richers,
Kohsuke Sumiyoshi,
Shoichi Yamada
Abstract:
We perform an extensive study of the influence of nuclear weak interactions on core-collapse supernovae (CCSNe), paying particular attention to consistency between nuclear abundances in the equation of state (EOS) and nuclear weak interactions. We compute properties of uniform matter based on the variational method. For inhomogeneous nuclear matter, we take a full ensemble of nuclei into account w…
▽ More
We perform an extensive study of the influence of nuclear weak interactions on core-collapse supernovae (CCSNe), paying particular attention to consistency between nuclear abundances in the equation of state (EOS) and nuclear weak interactions. We compute properties of uniform matter based on the variational method. For inhomogeneous nuclear matter, we take a full ensemble of nuclei into account with various finite-density and thermal effects and directly use the nuclear abundances to compute nuclear weak interaction rates. To quantify the impact of a consistent treatment of nuclear abundances on CCSN dynamics, we carry out spherically symmetric CCSN simulations with full Boltzmann neutrino transport, systematically changing the treatment of weak interactions, EOSs, and progenitor models. We find that the inconsistent treatment of nuclear abundances between the EOS and weak interaction rates weakens the EOS dependence of both the dynamics and neutrino signals. We also test the validity of two artificial prescriptions for weak interactions of light nuclei and find that both prescriptions affect the dynamics. Furthermore, there are differences in neutrino luminosities by ~10% and in average neutrino energies by 0.25-1 MeV from those of the fiducial model. We also find that the neutronization burst neutrino signal depends on the progenitor more strongly than on the EOS, preventing a detection of this signal from constraining the EOS.
△ Less
Submitted 23 December, 2018;
originally announced December 2018.
-
Sensitivity of nuclear statistical equilibrium to nuclear uncertainties during stellar core collapse
Authors:
Shun Furusawa
Abstract:
I have systematically investigated the equations of state (EOSs) in nuclear statistical equilibrium under thermodynamic conditions relevant for core collapse of massive stars by varying the bulk properties of nuclear matter, the mass data for neutron-rich nuclei, and the finite-temperature modifications of the nuclear model. It is found that the temperature dependence of the nuclear free energies…
▽ More
I have systematically investigated the equations of state (EOSs) in nuclear statistical equilibrium under thermodynamic conditions relevant for core collapse of massive stars by varying the bulk properties of nuclear matter, the mass data for neutron-rich nuclei, and the finite-temperature modifications of the nuclear model. It is found that the temperature dependence of the nuclear free energies has a significant impact on the entropy and nuclear composition, which affect the dynamics of core-collapse supernovae. There is a little influence from the bulk properties and the mass data. For all models, common nuclei that are likely to contribute to core-deleptonization are those near $Z\approx30$ and $N\approx50$. A model with a semi-empirical expression for internal degrees of freedom, however, overestimates the number densities of magic nuclei with $N\approx50$ and $82$, while a model, in which nuclear shell effects are not considered, underestimates the number densities of heavy nuclei, and especially of the magic nuclei. Other models, which include the temperature dependence of shell effects in the internal degrees of freedom and/or in the nuclear internal free energy, indicate that the difference in population between magic nuclei and non-magic nuclei disappears as the temperature increases. The construction of complete statistical EOS will require further theoretical and experimental studies of medium-mass, neutron-rich nuclei with proton numbers 25-45 and neutron numbers 40-85.
△ Less
Submitted 30 November, 2018; v1 submitted 26 November, 2018;
originally announced November 2018.
-
On the Neutrino Distributions in Phase Space for the Rotating Core-collapse Supernova Simulated with a Boltzmann-neutrino-radiation-hydrodynamics Code
Authors:
Akira Harada,
Hiroki Nagakura,
Wakana Iwakami,
Hirotada Okawa,
Shun Furusawa,
Hideo Matsufuru,
Kohsuke Sumiyoshi,
Shoichi Yamada
Abstract:
With the Boltzmann-radiation-hydrodynamics code, which we have developed to solve numerically the Boltzmann equations for neutrino transfer, the Newtonian hydrodynamics equations, and the Newtonian self-gravity simultaneously and consistently, we simulate the collapse of a rotating core of the progenitor with a zero-age-main-sequence mass of $11.2\,M_\odot$ and a shelluler rotation of…
▽ More
With the Boltzmann-radiation-hydrodynamics code, which we have developed to solve numerically the Boltzmann equations for neutrino transfer, the Newtonian hydrodynamics equations, and the Newtonian self-gravity simultaneously and consistently, we simulate the collapse of a rotating core of the progenitor with a zero-age-main-sequence mass of $11.2\,M_\odot$ and a shelluler rotation of $1\,{\rm rad\,s^{-1}}$ at the center. We pay particular attention in this paper to the neutrino distribution in phase space, which is affected by the rotation. By solving the Boltzmann equations directly, we can assess the rotation-induced distortion of the angular distribution in momentum space, which gives rise to the rotational component of the neutrino flux. We compare the Eddington tensors calculated both from the raw data and from the M1-closure approximation. We demonstrate that the Eddington tensor is determined by complicated interplays of the fluid velocity and the neutrino interactions and that the M1-closure, which assumes that the Eddington factor is determined by the flux factor, fails to fully capture this aspect, especially in the vicinity of the shock. We find that the error in the Eddington factor reaches $\sim 20\%$ in our simulation. This is due not to the resolution but to the different dependence of the Eddington and flux factors on the angular profile of the neutrino distribution function, and hence modification to the closure relation is needed.
△ Less
Submitted 6 March, 2019; v1 submitted 29 October, 2018;
originally announced October 2018.
-
A new equation of state for core-collapse supernovae based on realistic nuclear forces and including a full nuclear ensemble
Authors:
S. Furusawa,
H. Togashi,
H. Nagakura,
K. Sumiyoshi,
S. Yamada,
H. Suzuki,
M. Takano
Abstract:
We have constructed a nuclear equation of state (EOS) that includes a full nuclear ensemble for use in core-collapse supernova simulations. It is based on the EOS for uniform nuclear matter that two of the authors derived recently, applying a variational method to realistic two- and there-body nuclear forces. We utilize an extended liquid drop model of heavy nuclei and a quantum-theoretical mass e…
▽ More
We have constructed a nuclear equation of state (EOS) that includes a full nuclear ensemble for use in core-collapse supernova simulations. It is based on the EOS for uniform nuclear matter that two of the authors derived recently, applying a variational method to realistic two- and there-body nuclear forces. We utilize an extended liquid drop model of heavy nuclei and a quantum-theoretical mass evaluation for light nuclei. In addition to realistic nuclear forces, the inclusion of in-medium effects on the full ensemble of nuclei makes the new EOS one of the most realistic EOS's for supernova simulations. We make comparisons with the FYSS EOS, which is based on the same formulation for the nuclear ensemble but adopts the relativistic mean field (RMF) theory with the TM1 parameter set for uniform nuclear matter. The new EOS is softer than the FYSS EOS around and above nuclear saturation densities. We find that neutron-rich nuclei with small mass numbers are more abundant in the new EOS than in the FYSS EOS because of the larger saturation densities and smaller symmetry energy of nuclei in the former. We apply the two EOS's to 1D supernova simulations and find that the new EOS gives lower electron fractions and higher temperatures in the collapse phase owing to the smaller symmetry energy. As a result, the inner core has smaller masses for the new EOS. It is more compact, on the other hand, due to the softness of the new EOS and bounces at higher densities. The ensuing outward propagations of the shock wave in the outer core are very similar in the two simulations, which may be an artifact, though, caused by the use of the same tabulated electron capture rates for heavy nuclei ignoring differences in the nuclear composition between the two EOS's in these computations.
△ Less
Submitted 20 July, 2017;
originally announced July 2017.
-
Equilibrium nuclear ensembles taking into account vaporization of hot nuclei in dense stellar matter
Authors:
Shun Furusawa,
Igor Mishustin
Abstract:
We investigate the high-temperature effect on the nuclear matter that consists of mixture of nucleons and all nuclei in the dense and hot stellar environment. The individual nuclei are described within the compressible liquid-drop model that is based on Skyrme interactions for bulk energies and that takes into account modifications of the surface and Coulomb energies at finite temperatures and den…
▽ More
We investigate the high-temperature effect on the nuclear matter that consists of mixture of nucleons and all nuclei in the dense and hot stellar environment. The individual nuclei are described within the compressible liquid-drop model that is based on Skyrme interactions for bulk energies and that takes into account modifications of the surface and Coulomb energies at finite temperatures and densities. The free-energy density is minimized with respect to the individual equilibrium densities of all heavy nuclei and the nuclear composition. We find that their optimized equilibrium densities become smaller and smaller at high temperatures because of the increase of thermal contributions to bulk free energies and the reduction of surface energies. The neutron-rich nuclei become unstable and disappear one after another at some temperatures. The calculations are performed for two sets of model parameters leading to different values of the slope parameter in the nuclear symmetry energy. It is found that the larger slope parameter reduces the equilibrium densities and the melting temperatures. We also compare the new model with some other approaches and find that the mass fractions of heavy nuclei in the previous calculations that omit vaporization are underestimated at $T \lesssim 10$MeV and overestimated at $T\gtrsim10$~MeV. The further sophistication of calculations of nuclear vaporization and of light clusters would be required to construct the equation of state for explosive astrophysical phenomena.
△ Less
Submitted 23 February, 2018; v1 submitted 1 July, 2017;
originally announced July 2017.
-
Neutrino emissions in all flavors up to the pre-bounce of massive stars and the possibility of their detections
Authors:
Chinami Kato,
Hiroki Nagakura,
Shun Furusawa,
Koh Takahashi,
Hideyuki Umeda,
Takashi Yoshida,
Koji Ishidoshiro,
Shoichi Yamada
Abstract:
This paper is a sequel to our previous one (Kato et al.2015), which calculated the luminosities and spectra of electron-type anti-neutrinos ($\barν_e$'s) from the progenitors of core-collapse supernovae. Expecting that a capability to detect electron-type neutrinos ($ν_e$'s) will increase dramatically with the emergence of liquid-argon detectors such as DUNE, we broaden the scope in this study to…
▽ More
This paper is a sequel to our previous one (Kato et al.2015), which calculated the luminosities and spectra of electron-type anti-neutrinos ($\barν_e$'s) from the progenitors of core-collapse supernovae. Expecting that a capability to detect electron-type neutrinos ($ν_e$'s) will increase dramatically with the emergence of liquid-argon detectors such as DUNE, we broaden the scope in this study to include all-flavors of neutrinos emitted from the pre-bounce phase. We pick up three progenitor models of an electron capture supernova (ECSN) and iron-core collapse supernovae (FeCCSNe). We find that the number luminosities reach $\sim10^{57} \mathrm{s^{-1}}$ and $\sim10^{53} \mathrm{s^{-1}}$ at maximum for $ν_e$ and $\barν_e$, respectively. We also estimate the numbers of detection events at terrestrial neutrino detectors including DUNE, taking flavor oscillations into account and assuming the distance to the progenitors to be 200 pc. It is demonstrated that $\barν_e$'s from the ECSN-progenitor will be undetected at almost all detectors, whereas we will be able to observe $\gtrsim$15900 $ν_e$'s at DUNE for the inverted mass hierarchy. From the FeCCSN-progenitors, the number of $\barν_e$ events will be largest for JUNO, 200-900 $\barν_e$'s, depending on the mass hierarchy whereas the number of $ν_e$ events at DUNE is $\gtrsim$2100 for the inverted mass hierarchy. These results imply that the detection of $\barν_e$'s is useful to distinguish FeCCSN- from ECSN-progenitors, while $ν_e$'s will provide us with detailed information on the collapse phase regardless of the type and mass of progenitor.
△ Less
Submitted 18 September, 2017; v1 submitted 18 April, 2017;
originally announced April 2017.
-
Simulations of core-collapse supernovae in spatial axisymmetry with full Boltzmann neutrino transport
Authors:
Hiroki Nagakura,
Wakana Iwakami,
Shun Furusawa,
Hirotada Okawa,
Akira Harada,
Kohsuke Sumiyoshi,
Shoichi Yamada,
Hideo Matsufuru,
Akira Imakura
Abstract:
We present the first results of our spatially axisymmetric core-collapse supernova simulations with full Boltzmann neutrino transport, which amount to a time-dependent 5-dimensional (2 in space and 3 in momentum space) problem in fact. Special relativistic effects are fully taken into account with a two-energy-grid technique. We performed two simulations for a progenitor of 11.2M, employing differ…
▽ More
We present the first results of our spatially axisymmetric core-collapse supernova simulations with full Boltzmann neutrino transport, which amount to a time-dependent 5-dimensional (2 in space and 3 in momentum space) problem in fact. Special relativistic effects are fully taken into account with a two-energy-grid technique. We performed two simulations for a progenitor of 11.2M, employing different nuclear equations-of-state (EOS's): Lattimer and Swesty's EOS with the incompressibility of K = 220MeV (LS EOS) and Furusawa's EOS based on the relativistic mean field theory with the TM1 parameter set (FS EOS). In the LS EOS the shock wave reaches ~700km at 300ms after bounce and is still expanding whereas in the FS EOS it stalled at ~200km and has started to recede by the same time. This seems to be due to more vigorous turbulent motions in the former during the entire post-bounce phase, which leads to higher neutrino-heating efficiency in the neutrino-driven convection. We also look into the neutrino distributions in momentum space, which is the advantage of the Boltzmann transport over other approximate methods. We find non-axisymmetric angular distributions with respect to the local radial direction, which also generate off-diagonal components of the Eddington tensor. We find that the r θ-component reaches ~10% of the dominant rr-component and, more importantly, it dictates the evolution of lateral neutrino fluxes, dominating over the θθ-component, in the semi-transparent region. These data will be useful to further test and possibly improve the prescriptions used in the approximate methods.
△ Less
Submitted 12 February, 2018; v1 submitted 6 February, 2017;
originally announced February 2017.
-
Dependence of weak interaction rates on the nuclear composition during stellar core collapse
Authors:
Shun Furusawa,
Hiroki Nagakura,
Kohsuke Sumiyoshi,
Chinami Kato,
Shoichi Yamada
Abstract:
We investigate the influences of the nuclear composition on the weak interaction rates of heavy nuclei during the core collapse of massive stars. The nuclear abundances in nuclear statistical equilibrium (NSE) are calculated by some equation of state (EOS) models including in-medium effects on nuclear masses. We systematically examine the sensitivities of electron capture and neutrino-nucleus scat…
▽ More
We investigate the influences of the nuclear composition on the weak interaction rates of heavy nuclei during the core collapse of massive stars. The nuclear abundances in nuclear statistical equilibrium (NSE) are calculated by some equation of state (EOS) models including in-medium effects on nuclear masses. We systematically examine the sensitivities of electron capture and neutrino-nucleus scattering on heavy nuclei to the nuclear shell effects and the single nucleus approximation. We find that the washout of shell effects at high temperatures brings significant change to weak rates by smoothing the nuclear abundance distribution: the electron capture rate decreases by $\sim$20$\%$ in the early phase and increases by $\sim$40$\%$ in the late phase at most, while the cross section for neutrino-nucleus scattering is reduced by $\sim$15$\%$. This is because the open-shell nuclei become abundant instead of those with closed neutron shells as the shell effects disappear. We also find that the single-nucleus description based on the average values leads to underestimations of weak rates. Electron captures and neutrino coherent scattering on heavy nuclei are reduced by $\sim$80$\%$ in the early phase and by $\sim$5$\%$ in the late phase, respectively. These results indicate that NSE like EOS accounting for shell washout is indispensable for the reliable estimation of weak interaction rates in simulations of core-collapse supernovae.
△ Less
Submitted 18 March, 2017; v1 submitted 29 January, 2017;
originally announced January 2017.
-
Self consistent calculation of the nuclear composition in hot and dense stellar matter
Authors:
Shun Furusawa,
Igor Mishustin
Abstract:
We investigate the mass fractions and in-medium properties of heavy nuclei in stellar matter at characteristic densities and temperatures for supernova (SN) explosions. The individual nuclei are described within the compressible liquid-drop model taking into account modifications of bulk, surface and Coulomb energies. The equilibrium properties of nuclei and the full ensemble of heavy nuclei are c…
▽ More
We investigate the mass fractions and in-medium properties of heavy nuclei in stellar matter at characteristic densities and temperatures for supernova (SN) explosions. The individual nuclei are described within the compressible liquid-drop model taking into account modifications of bulk, surface and Coulomb energies. The equilibrium properties of nuclei and the full ensemble of heavy nuclei are calculated self-consistently. It is found that heavy nuclei in the ensemble are either compressed or decompressed depending on the isospin asymmetry of the system. The compression or decompression has a little influence on the binding energies, total mass fractions and average mass numbers of heavy nuclei, although the equilibrium densities of individual nuclei themselves are changed appreciably above one hundredth of normal nuclear density. We find that nuclear structure in single nucleus approximation deviates from the actual one obtained in the multi-nucleus description, since the density of free nucleons is different between these two descriptions. This study indicates that a multi-nucleus description is required to realistically account for in-medium effects on the nuclear structure in supernova matter.
△ Less
Submitted 18 March, 2017; v1 submitted 6 December, 2016;
originally announced December 2016.
-
Supernova equations of state including full nuclear ensemble with in-medium effects
Authors:
Shun Furusawa,
Kohsuke Sumiyoshi,
Shoichi Yamada,
Hideyuki Suzuki
Abstract:
We construct new equations of state for baryons at sub-nuclear densities for the use in core-collapse supernova simulations. The abundance of various nuclei is obtained together with thermodynamic quantities. The formulation is an extension of the previous model, in which we adopted the relativistic mean field theory with the TM1 parameter set for nucleons, the quantum approach for $d$, $t$, $h$ a…
▽ More
We construct new equations of state for baryons at sub-nuclear densities for the use in core-collapse supernova simulations. The abundance of various nuclei is obtained together with thermodynamic quantities. The formulation is an extension of the previous model, in which we adopted the relativistic mean field theory with the TM1 parameter set for nucleons, the quantum approach for $d$, $t$, $h$ and $α$ as well as the liquid drop model for the other nuclei under the nuclear statistical equilibrium. We reformulate the model of the light nuclei other than $d$, $t$, $h$ and $α$ based on the quasi-particle description. Furthermore, we modify the model so that the temperature dependences of surface and shell energies of heavy nuclei could be taken into account. The pasta phases for heavy nuclei and the Pauli- and self-energy shifts for $d$, $t$, $h$ and $α$ are taken into account in the same way as in the previous model. We find that nuclear composition is considerably affected by the modifications in this work, whereas thermodynamical quantities are not changed much. In particular, the washout of shell effect has a great impact on the mass distribution above $T \sim 1$ MeV. This improvement may have an important effect on the rates of electron captures and coherent neutrino scatterings on nuclei in supernova cores.
△ Less
Submitted 20 July, 2017; v1 submitted 6 December, 2016;
originally announced December 2016.
-
Three-dimensional Boltzmann-Hydro code for core-collapse in massive stars II. The Implementation of moving-mesh for neutron star kicks
Authors:
Hiroki Nagakura,
Wakana Iwakami,
Shun Furusawa,
Kohsuke Sumiyoshi,
Shoichi Yamada,
Hideo Matsufuru,
Akira Imakura
Abstract:
We present a newly developed moving-mesh technique for the multi-dimensional Boltzmann-Hydro code for the simulation of core-collapse supernovae (CCSNe). What makes this technique different from others is the fact that it treats not only hydrodynamics but also neutrino transfer in the language of the 3+1 formalism of general relativity (GR), making use of the shift vector to specify the time evolu…
▽ More
We present a newly developed moving-mesh technique for the multi-dimensional Boltzmann-Hydro code for the simulation of core-collapse supernovae (CCSNe). What makes this technique different from others is the fact that it treats not only hydrodynamics but also neutrino transfer in the language of the 3+1 formalism of general relativity (GR), making use of the shift vector to specify the time evolution of the coordinate system. This means that the transport part of our code is essentially general relativistic although in this paper it is applied only to the moving curvilinear coordinates in the flat Minknowski spacetime, since the gravity part is still Newtonian. The numerical aspect of the implementation is also described in detail. Employing the axisymmetric two-dimensional version of the code, we conduct two test computations: oscillations and runaways of proto-neutron star (PNS). We show that our new method works fine, tracking the motions of PNS correctly. We believe that this is a major advancement toward the realistic simulation of CCSNe.
△ Less
Submitted 2 May, 2016;
originally announced May 2016.
-
A Hydrodynamical Study on the Conversion of Hadronic Matter to Quark Matter: II. Diffusion-Induced Conversion
Authors:
Shun Furusawa,
Takahiro Sanada,
Shoichi Yamada
Abstract:
We study transitions of hadronic matter (HM) to 3-flavor quark matter (3QM), regarding the conversion processes as combustion and describing them hydrodynamically. Under the assumption that HM is metastable with their free energies being larger than those of 3QM but smaller than those of 2-flavor quark matter (2QM), we consider in this paper the conversion induced by diffusions of seed 3QM. This i…
▽ More
We study transitions of hadronic matter (HM) to 3-flavor quark matter (3QM), regarding the conversion processes as combustion and describing them hydrodynamically. Under the assumption that HM is metastable with their free energies being larger than those of 3QM but smaller than those of 2-flavor quark matter (2QM), we consider in this paper the conversion induced by diffusions of seed 3QM. This is a sequel to our previous paper, in which the shock-induced conversion was studied in the same frame work. We not only pay attention to the jump condition on both sides of the conversion front but the structures inside the front are also considered by taking into account what happens during the conversion processes on the time scale of weak interactions. We employ for HM the Shen's EOS, which is based on the relativistic mean field theory, and the bag model-based EOS for QM just as in the previous paper. We demonstrated in that paper that in this combination of EOS's the combustion will occur for a wide range of the bag constant and strong coupling constant in the so-called endothermic regime, in which the Hugoniot curve for combustion runs below the initial state. We find that weak deflagration nearly always occurs and that weak detonation is possible only when the diffusion constant is (unrealistically) large and the critical strange fraction is small. The velocities of the conversion front are ~ $10^3-10^7$ cm/s depending on the initial temperature and density as well as the parameters in the QM EOS and become particularly small when the final state is in the mixed phase. Finally we study linear stability of the laminar weak-deflagration front and find that it is unstable in the exothermic regime (Darrius-Landau instability) but stable in the endothermic regime, which is quite contrary to the ordinary combustions.
△ Less
Submitted 4 February, 2016; v1 submitted 24 November, 2015;
originally announced November 2015.
-
A Hydrodynamical Study on the Conversion of Hadronic Matter to Quark Matter: I. Shock-Induced Conversion
Authors:
Shun Furusawa,
Takahiro Sanada,
Shoichi Yamada
Abstract:
We study transitions of hadronic matter (HM) to 3-flavor quark matter (3QM) locally, regarding the conversion processes as combustion and describing them hydrodynamically. Not only the jump condition on both sides of the conversion front but the structures inside the front are also considered by taking into account what happens during the conversion processes on the time scale of weak interactions…
▽ More
We study transitions of hadronic matter (HM) to 3-flavor quark matter (3QM) locally, regarding the conversion processes as combustion and describing them hydrodynamically. Not only the jump condition on both sides of the conversion front but the structures inside the front are also considered by taking into account what happens during the conversion processes on the time scale of weak interactions as well as equations of state (EOS's) in the mixed phase. Under the assumption that HM is metastable with their free energies being larger than those of 3QM but smaller than those of 2-flavor quark matter (2QM), we consider the transition via 2QM triggered by a rapid density rise in a shock wave. Based on the results, we discuss which combustion modes (strong/weak detonation) may be realized. HM is described by an EOS based on the relativistic mean field theory and 2, 3QM's are approximated by the MIT bag model. We demonstrate for a wide range of bag constant and strong coupling constant in this combination of EOS's that the combustion may occur in the so-called endothermic regime, in which the Hugoniot curve for combustion runs below the one for the shock wave in P-V plane, and which has no terrestrial counter part. We find that strong detonation always occurs. Depending on the EOS of quark matter (QM) as well as the density of HM and the Mach number of the detonation front, deconfinement from HM to 2QM is either completed or not completed in the shock wave. In the latter case, which is more likely if the EOS of QM ensures that deconfinement occurs above the nuclear saturation density and that the maximum mass of cold quark stars is larger than two solar mass, the conversion continues further via the mixing state of HM and 3QM on the time scale of weak interactions.
△ Less
Submitted 30 November, 2015; v1 submitted 24 November, 2015;
originally announced November 2015.
-
The influence of inelastic neutrino reactions with light nuclei on the standing accretion shock instability in core-collapse supernovae
Authors:
Shun Furusawa,
Hiroki Nagakura,
Kohsuke Sumiyoshi,
Shoichi Yamada
Abstract:
We perform numerical experiments to investigate the influence of inelastic neutrino reactions with light nuclei on the standing accretion shock instability (SASI). The time evolution of shock waves is calculated with a simple light-bulb approximation for the neutrino transport and a multi-nuclei equation of state. The neutrino absorptions and inelastic interactions with deuterons, tritons, helions…
▽ More
We perform numerical experiments to investigate the influence of inelastic neutrino reactions with light nuclei on the standing accretion shock instability (SASI). The time evolution of shock waves is calculated with a simple light-bulb approximation for the neutrino transport and a multi-nuclei equation of state. The neutrino absorptions and inelastic interactions with deuterons, tritons, helions and alpha particles are taken into account in the hydrodynamical simulations in addition to the ordinary charged-current interactions with nucleons. Axial symmetry is assumed but no equatorial symmetry is imposed. We show that the heating rates of deuterons reach as high as 10% of those of nucleons around the bottom of the gain region. On the other hand, alpha particles are heated near the shock wave, which is important when the shock wave expands and the density and temperature of matter become low. It is also found that the models with heating by light nuclei have different evolutions from those without it in the non-linear phase of SASI. This results is because matter in the gain region has a varying density and temperature and there appear sub-regions that are locally rich in deuterons and alpha particles. Although the light nuclei are never dominant heating sources and they work favorably for shock revival in some cases and unfavorably in other cases, they are non-negligible and warrant further investigation.
△ Less
Submitted 7 August, 2013; v1 submitted 7 May, 2013;
originally announced May 2013.
-
New equations of state based on the liquid drop model of heavy nuclei and quantum approach to light nuclei for core-collapse supernova simulations
Authors:
Shun Furusawa,
Kohsuke Sumiyoshi,
Shoichi Yamada,
Hideyuki Suzuki
Abstract:
We construct new equations of state for baryons at sub-nuclear densities for the use in core-collapse simulations of massive stars. The abundance of various nuclei is obtained together with thermodynamic quantities. A model free energy is constructed, based on the relativistic mean field theory for nucleons and the mass formula for nuclei with the proton number up to ~ 1000. The formulation is an…
▽ More
We construct new equations of state for baryons at sub-nuclear densities for the use in core-collapse simulations of massive stars. The abundance of various nuclei is obtained together with thermodynamic quantities. A model free energy is constructed, based on the relativistic mean field theory for nucleons and the mass formula for nuclei with the proton number up to ~ 1000. The formulation is an extension of the previous model, in which we adopted the liquid drop model to all nuclei under the nuclear statistical equilibrium. We reformulate the new liquid drop model so that the temperature dependences of bulk energies could be taken into account. Furthermore, we extend the region in the nuclear chart, in which shell affects are included, by using theoretical mass data in addition to experimental ones. We also adopt a quantum theoretical mass evaluation of light nuclei, which incorporates the Pauli- and self-energy shifts that are not included in the ordinary liquid drop model. The pasta phases for heavy nuclei are taken into account in the same way as in the previous model. We find that the abundances of heavy nuclei are modified by the shell effects of nuclei and temperature dependence of bulk energies. These changes may have an important effect on the rates of electron captures and coherent neutrino scatterings on nuclei in supernova cores. The abundances of light nuclei are also modified by the new mass evaluation, which may affect the heating and cooling rates of supernova cores and shocked envelopes.
△ Less
Submitted 25 June, 2013; v1 submitted 7 May, 2013;
originally announced May 2013.
-
A comparative study of statistical models for nuclear equation of state of stellar matter
Authors:
N. Buyukcizmeci,
A. S. Botvina,
I. N. Mishustin,
R. Ogul,
M. Hempel,
J. Schaffner-Bielich,
F. -K. Thielemann,
S. Furusawa,
K. Sumiyoshi,
S. Yamada,
H. Suzuki
Abstract:
We compare three different statistical models for the equation of state (EOS) of stellar matter at subnuclear densities and temperatures (0.5-10 MeV) expected to occur during the collapse of massive stars and supernova explosions. The models introduce the distributions of various nuclear species in nuclear statistical equilibrium, but use somewhat different nuclear physics inputs. It is demonstrat…
▽ More
We compare three different statistical models for the equation of state (EOS) of stellar matter at subnuclear densities and temperatures (0.5-10 MeV) expected to occur during the collapse of massive stars and supernova explosions. The models introduce the distributions of various nuclear species in nuclear statistical equilibrium, but use somewhat different nuclear physics inputs. It is demonstrated that the basic thermodynamical quantities of stellar matter under these conditions are similar, except in the region of high densities and low temperatures. We demonstrate that mass and isotopic distributions have considerable differences related to the different assumptions of the models on properties of nuclei at these stellar conditions. Overall, the three models give similar trends, but the details reflect the uncertainties related to the modeling of medium effects, such as the temperature and density dependence of surface and bulk energies of heavy nuclei, and the nuclear shell structure effects. We discuss importance of new physics inputs for astrophysical calculations from experimental data obtained in intermediate energy heavy-ion collisions, in particular, the similarities of the conditions reached during supernova explosions and multifragmentation reactions.
△ Less
Submitted 22 March, 2013; v1 submitted 26 November, 2012;
originally announced November 2012.
-
Hyperon Matter and Black Hole Formation in Failed Supernovae
Authors:
Ken'ichiro Nakazato,
Shun Furusawa,
Kohsuke Sumiyoshi,
Akira Ohnishi,
Shoichi Yamada,
Hideyuki Suzuki
Abstract:
We investigate the emergence of hyperons in black-hole-forming failed supernovae, which are caused by the dynamical collapse of nonrotating massive stars. We perform neutrino-radiation hydrodynamical simulations in general relativity adopting realistic hyperonic equation-of-state (EOS). Attractive and repulsive cases are examined for the potential of Σ hyperons. Since hyperons soften the EOS, they…
▽ More
We investigate the emergence of hyperons in black-hole-forming failed supernovae, which are caused by the dynamical collapse of nonrotating massive stars. We perform neutrino-radiation hydrodynamical simulations in general relativity adopting realistic hyperonic equation-of-state (EOS). Attractive and repulsive cases are examined for the potential of Σ hyperons. Since hyperons soften the EOS, they shorten the time interval from the bounce to black hole formation, which corresponds to the duration of neutrino emission. This effect is larger for the attractive case than the repulsive case because Σ hyperons appear more easily. In addition, we investigate the impacts of pions to find that they also promotes the recollapse towards the black hole formation.
△ Less
Submitted 11 November, 2011;
originally announced November 2011.
-
A new baryonic equation of state at sub-nuclear densities for core-collapse simulations
Authors:
Shun Furusawa,
Shoichi Yamada,
Kohsuke Sumiyoshi,
Hideyuki Suzuki
Abstract:
We calculate a new equation of state for baryons at sub-nuclear densities meant for the use in core-collapse simulations of massive stars. The abundance of various nuclei is obtained together with the thermodynamic quantities. The formulation is the NSE description and the liquid drop approximation of nuclei. The model free energy to minimize is calculated by relativistic mean field theory for nuc…
▽ More
We calculate a new equation of state for baryons at sub-nuclear densities meant for the use in core-collapse simulations of massive stars. The abundance of various nuclei is obtained together with the thermodynamic quantities. The formulation is the NSE description and the liquid drop approximation of nuclei. The model free energy to minimize is calculated by relativistic mean field theory for nucleons and the mass formula for nuclei with the atomic number up to ~ 1000. We have also taken into account the pasta phase, thanks to which the transition to uniform nuclear matter in our EOS occurs in the conventional manner: nuclei are not dissociated to nucleons but survive right up to the transition to uniform nuclear matter. We find that the free energy and other thermodynamical quantities are not very different from those given in the Shen's EOS, one of the standard EOS's that adopt the single nucleus approximation. The average mass is systematically different, on the other hand, which may have an important ramification to the rates of electron captures and coherent neutrino scatterings on nuclei in supernova cores. It is also interesting that the root mean square of the mass number is not very different from the average mass number, since the former is important for the evaluation of coherent scattering rates on nuclei but has been unavailable so far. The EOS table is currently under construction, which will include the weak interaction rates.
△ Less
Submitted 9 August, 2011; v1 submitted 31 March, 2011;
originally announced March 2011.