Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userExplosion Energies for Core-collapse Supernovae I: Analytic, Spherically Symmetric Solutions
171
0
0.0
(
0
)
Added by
Mariam Gogilashvili
Publication date
2020
fields
Physics
and research's language is
English
Ask ChatGPT about the research
No Arabic abstract
Recent multi-dimensional simulations of core-collapse supernovae are producing successful explosions and explosion-energy predictions. In general, the explosion-energy evolution is monotonic and relatively smooth, suggesting a possible analytic solution. We derive analytic solutions for the expansion of the gain region under the following assumptions: spherical symmetry, one-zone shell, and powered by neutrinos and $alpha$ particle recombination. We consider two hypotheses: I) explosion energy is powered by neutrinos and $alpha$ recombination, II) explosion energy is powered by neutrinos alone. Under these assumptions, we derive the fundamental dimensionless parameters and analytic scalings. For the neutrino-only hypothesis (II), the asymptotic explosion energy scales as $E_{infty} approx 1.5 M_g v_0^2 eta^{2/3}$, where $M_g$ is the gain mass, $v_0$ is the free-fall velocity at the shock, and $eta$ is a ratio of the heating and dynamical time scales. Including both neutrinos and recombination (hypothesis I), the asymptotic explosion energy is $E_{infty} approx M_g v_0^2 (1.5eta^{2/3} + beta f(rho_0))$, where $beta$ is the dimensionless recombination parameter. We use Bayesian inference to fit these analytic models to simulations. Both hypotheses fit the simulations of the lowest progenitor masses that tend to explode spherically. The fits do not prefer hypothesis I or II; however, prior investigations suggest that $alpha$ recombination is important. As expected, neither hypothesis fits the higher-mass simulations that exhibit aspherical explosions. In summary, this explosion-energy theory is consistent with the spherical explosions of low progenitor masses; the inconsistency with higher progenitor-mass simulations suggests that a theory for them must include aspherical dynamics.
rate research
Read More
Neutrinos are a guaranteed signal from supernova explosions in the Milky Way, and a most valuable messenger that can provide us with information about the deepest parts of supernovae. In particular, neutrinos will provide us with physical quantities, such as the radius and mass of protoneutron stars (PNS), which are the central engine of supernovae. This requires a theoretical model that connects observables such as neutrino luminosity and average energy with physical quantities. Here, we show analytic solutions for the neutrino-light curve derived from the neutrino radiation transport equation by employing the diffusion approximation and the analytic density solution of the hydrostatic equation for a PNS. The neutrino luminosity and the average energy as functions of time are explicitly presented, with dependence on PNS mass, radius, the total energy of neutrinos, surface density, and opacity. The analytic solutions provide good representations of the numerical models from a few seconds after the explosion and allow a rough estimate of these physical quantities to be made from observational data.
Shock revival in core-collapse supernovae (CCSNe) may be due to the neutrino mechanism. While it is known that in a neutrino-powered CCSN, explosion begins when the neutrino luminosity of the proto-neutron star exceeds a critical value, the physics of this condition in time-dependent, multidimensional simulations are not fully understood. citet{Pejcha2012} found that an `antesonic condition exists for time-steady spherically symmetric models, potentially giving a physical explanation for the critical curve observed in simulations. In this paper, we extend that analysis to time-dependent, spherically symmetric polytropic models. We verify the critical antesonic condition in our simulations, showing that models exceeding it drive transonic winds whereas models below it exhibit steady accretion. In addition, we find that (1) high spatial resolution is needed for accurate determination of the antesonic ratio and shock radius at the critical curve, and that low resolution simulations systematically underpredict these quantities, making explosion more difficult at lower resolution; (2) there is an important physical connection between the critical mass accretion rate at explosion and the mass loss rate of the post-explosion wind: the two are directly proportional at criticality, implying that, at criticality, the wind kinetic power is tied directly to the accretion power; (3) the value of the post-shock adiabatic index $Gamma$ has a large effect on the length and time scales of the post-bounce evolution of the explosion larger values of $Gamma$ result in a longer transition from the accretion to wind phases.
We investigate the post-explosion phase in core-collapse supernovae with 2D hydrodynamical simulations and a simple neutrino treatment. The latter allows us to perform 46 simulations and follow the evolution of the 32 successful explosions during several seconds. We present a broad study based on three progenitors (11.2 $M_odot$, 15 $M_odot$, and 27 $M_odot$), different neutrino-heating efficiencies, and various rotation rates. We show that the first seconds after shock revival determine the final explosion energy, remnant mass, and properties of ejected matter. Our results suggest that a continued mass accretion increases the explosion energy even at late times. We link the late-time mass accretion to initial conditions such as rotation strength and shock deformation at explosion time. Only some of our simulations develop a neutrino-driven wind that survives for several seconds. This indicates that neutrino-driven winds are not a standard feature expected after every successful explosion. Even if our neutrino treatment is simple, we estimate the nucleosynthesis of the exploding models for the 15 $M_odot$ progenitor after correcting the neutrino energies and luminosities to get a more realistic electron fraction.
There are now $sim$20 multi-dimensional core-collapse supernova (CCSN) simulations that explode. However, these simulations have explosion energies that are a few times $10^{50}$ erg, not $10^{51}$ erg. In this manuscript, we compare the inferred explosion energies of these simulations and observations of 38 SN~IIP. Assuming a log-normal distribution, the mean explosion energy for the observations is $mu_{rm obs} = -0.13pm 0.05$ ($log_{10}(E/10^{51}, {rm erg})$) and the width is $sigma_{rm obs} = 0.21^{+0.05}_{-0.04}$. Only three CCSN codes have sufficient simulations to compare with observations: CHIMERA, CoCoNuT-FMT, and FORNAX. Currently, FORNAX has the largest sample of simulations. The two-dimensional FORNAX simulations show a correlation between explosion energy and progenitor mass, ranging from linear to quadratic, $E_{rm sim} propto M^{1-2}$; this correlation is consistent with inferences from observations. In addition, we infer the ratio of the observed-to-simulated explosion energies, $Delta=log_{10}(E_{rm obs}/E_{rm sim})$. For the CHIMERA set, $Delta=0.33pm0.06$; for CoCoNuT-FMT, $Delta=0.62pm0.05$; for FORNAX2D, $Delta=0.73pm0.05$, and for FORNAX3D, $Delta=0.95pm0.06$. On average, the simulations are less energetic than inferred energies from observations ($Delta approx 0.7$), but we also note that the variation among the simulations (max($Delta$)-min($Delta$) $approx 0.6$) is as large as this average offset. This suggests that further improvements to the simulations could resolve the discrepancy. Furthermore, both the simulations and the observations are heavily biased. In this preliminary comparison, we model these biases, but to more reliably compare the explosion energies, we recommend strategies to un-bias both the simulations and observations.
We present modelling of line polarization to study multi-dimensional geometry of stripped-envelope core-collapse supernovae (SNe). We demonstrate that a purely axisymmetric, two-dimensional geometry cannot reproduce a loop in the Stokes Q-U diagram, i.e., a variation of the polarization angles along the velocities associated with the absorption lines. On the contrary, three-dimensional (3D) clumpy structures naturally reproduce the loop. The fact that the loop is commonly observed in stripped-envelope SNe suggests that SN ejecta generally have a 3D structure. We study the degree of line polarization as a function of the absorption depth for various 3D clumpy models with different clump sizes and covering factors. Comparison between the calculated and observed degree of line polarization indicates that a typical size of the clump is relatively large, >~ 25 % of the photospheric radius. Such large-scale clumps are similar to those observed in the SN remnant Cassiopeia A. Given the small size of the observed sample, the covering factor of the clumps is only weakly constrained (~ 5-80 %). The presence of large-scale clumpy structure suggests that the large-scale convection or standing accretion shock instability takes place at the onset of the explosion.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
