No Arabic abstract
We investigate a method to construct parametrized progenitor models for core-collapse supernova simulations. Different from all modern core-collapse supernova studies, which rely on progenitor models from stellar evolution calculations, we follow the methodology of Baron & Cooperstein (1990) to construct initial models. Choosing parametrized spatial distributions of entropy and electron fraction as a function of mass coordinate and solving the equation of hydrostatic equilibrium, we obtain the initial density structures of our progenitor models. First, we calculate structures with parameters fitting broadly the evolutionary model s11.2 of Woosley et al. (2002). We then demonstrate the reliability of our method by performing general relativistic hydrodynamic simulations in spherical symmetry with the isotropic diffusion source approximation to solve the neutrino transport. Our comprehensive parameter study shows that initial models with a small central entropy ($lesssim 0.4,k_B$ nucleon$^{-1}$) can explode even in spherically symmetric simulations. Models with a large entropy ($gtrsim 6,k_B$ nucleon$^{-1}$) in the Si/O layer have a rather large explosion energy ($sim 4times 10^{50}$ erg) at the end of the simulations, which is still rapidly increasing.
We construct the equation of state (EOS) of dense matter covering a wide range of temperature, proton fraction, and density for the use of core-collapse supernova simulations. The study is based on the relativistic mean-field (RMF) theory, which can provide an excellent description of nuclear matter and finite nuclei. The Thomas--Fermi approximation in combination with assumed nucleon distribution functions and a free energy minimization is adopted to describe the non-uniform matter, which is composed of a lattice of heavy nuclei. We treat the uniform matter and non-uniform matter consistently using the same RMF theory. We present two sets of EOS tables, namely EOS2 and EOS3. EOS2 is an update of our earlier work published in 1998 (EOS1), where only the nucleon degree of freedom is taken into account. EOS3 includes additional contributions from $Lambda$ hyperons. The effect of $Lambda$ hyperons on the EOS is negligible in the low-temperature and low-density region, whereas it tends to soften the EOS at high density. In comparison with EOS1, EOS2 and EOS3 have an improved design of ranges and grids, which covers the temperature range $T=0.1$--$10^{2.6}$ MeV with the logarithmic grid spacing $Delta log_{10}(T/rm{[MeV]})=0.04$ (92 points including T=0), the proton fraction range $Y_p=0$--0.65 with the linear grid spacing $Delta Y_p = 0.01$ (66 points), and the density range $rho_B=10^{5.1}$--$10^{16},rm{g,cm^{-3}}$ with the logarithmic grid spacing $Delta log_{10}(rho_B/rm{[g,cm^{-3}]}) = 0.1$ (110 points).
Core-collapse supernovae are among Natures grandest explosions. They are powered by the energy released in gravitational collapse and include a rich set of physical phenomena involving all fundamental forces and many branches of physics and astrophysics. We summarize the current state of core-collapse supernova theory and discuss the current set of candidate explosion mechanisms under scrutiny as core-collapse supernova modeling is moving towards self-consistent three-dimensional simulations. Recent work in nuclear theory and neutron star mass and radius measurements are providing new constraints for the nuclear equation of state. We discuss these new developments and their impact on core-collapse supernova modeling. Neutrino-neutrino forward scattering in the central regions of core-collapse supernovae can lead to collective neutrino flavor oscillations that result in swaps of electron and heavy-lepton neutrino spectra. We review the rapid progress that is being made in understanding these collective oscillations and their potential impact on the core-collapse supernova explosion mechanism.
We present multi-dimensional core-collapse supernova simulations using the Isotropic Diffusion Source Approximation (IDSA) for the neutrino transport and a modified potential for general relativity in two different supernova codes: FLASH and ELEPHANT. Due to the complexity of the core-collapse supernova explosion mechanism, simulations require not only high-performance computers and the exploitation of GPUs, but also sophisticated approximations to capture the essential microphysics. We demonstrate that the IDSA is an elegant and efficient neutrino radiation transfer scheme, which is portable to multiple hydrodynamics codes and fast enough to investigate long-term evolutions in two and three dimensions. Simulations with a 40 solar mass progenitor are presented in both FLASH (1D and 2D) and ELEPHANT (3D) as an extreme test condition. It is found that the black hole formation time is delayed in multiple dimensions and we argue that the strong standing accretion shock instability before black hole formation will lead to strong gravitational waves.
We write down and apply the linearized fluid and gravitational equations consistent with pseudo-Newtonian simulations, whereby Newtonian hydrodynamics is used with a pseudo-Newtonian monopole and standard Newtonian gravity for higher multipoles. We thereby eliminate the need to use mode function matching to identify the active non-radial modes in pseudo-Newtonian core-collapse supernova simulations, in favor of the less complex and less costly mode frequency matching method. In doing so, we are able to measure appropriate boundary conditions for a mode calculation.
We present a comparison between several simulation codes designed to study the core-collapse supernova mechanism. We pay close attention to controlling the initial conditions and input physics in order to ensure a meaningful and informative comparison. Our goal is three-fold. First, we aim to demonstrate the current level of agreement between various groups studying the core-collapse supernova central engine. Second, we desire to form a strong basis for future simulation codes and methods to compare to. Lastly, we want this work to be a stepping stone for future work exploring more complex simulations of core-collapse supernovae, i.e., simulations in multiple dimensions and simulations with modern neutrino and nuclear physics. We compare the early (first ~500ms after core bounce) spherically-symmetric evolution of a 20 solar mass progenitor star from six different core-collapse supernovae codes: 3DnSNe-IDSA, AGILE-BOLTZTRAN, FLASH, F{sc{ornax}}, GR1D, and PROMETHEUS-VERTEX. Given the diversity of neutrino transport and hydrodynamic methods employed, we find excellent agreement in many critical quantities, including the shock radius evolution and the amount of neutrino heating. Our results provide an excellent starting point from which to extend this comparison to higher dimensions and compare the development of hydrodynamic instabilities that are crucial to the supernova explosion mechanism, such as turbulence and convection.