Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userPrincipal-Axis Analysis of the Eddington Tensor for the Early Post-Bounce Phase of Rotational Core-Collapse Supernovae
89
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
One- (1D) and two-dimensional (2D) core-collapse supernova simulations using full Boltzmann neutrino transport for 11.2M and 15.0M progenitor models have been performed to verify the closure relation for the moment method used in the approximate radiation transfer. This study finds areas where the results of the closure relation are inconsistent with those of Boltzmann transport, even for rotational models. In 1D simulations, the Eddington factors p defined in the fluid rest frame (FR) are compared to evaluate the maximum entropy closure for the Fermi-Dirac distribution (MEFD), confirming that MEFD closure performs better than other closures if p < 1/3 and phase space occupancy e > 0.5. In 2D simulations for non-rotating progenitor models, similar results are obtained from the principal-axis analysis of the Eddington tensor kij measured in FR. However, for rotating progenitor models, the principal axes of kij for Boltzmann transport tilt toward oblique directions where matter and neutrinos move relatively fast in azimuthal directions, while the principal axes of kij for MEFD closure are always parallel or perpendicular to the neutrino flux. Thus, the assumption of axisymmetric angular distribution to the flux direction in the closure relation does not hold in the strongly rotating supernova core in the early post-bounce phase. It is also shown that the deviation of the principal axes of kij from the flux direction increases when evaluated in a laboratory frame (LB). The optically thin and thick terms of the pressure tensor in LB negatively impact results in optically thicker and thinner regions, respectively.
rate research
Read 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 Togashis 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.
We explore explosions of massive stars, which are triggered via the quark-hadron phase transition during the early post bounce phase of core-collapse supernovae. We construct a quark equation of state, based on the bag model for strange quark matter. The transition between the hadronic and the quark phases is constructed applying Gibbs conditions. The resulting quark-hadron hybrid equations of state are used in core-collapse supernova simulations, based on general relativistic radiation hydrodynamics and three flavor Boltzmann neutrino transport in spherical symmetry. The formation of a mixed phase reduces the adiabatic index, which induces the gravitational collapse of the central protoneutron star. The collapse halts in the pure quark phase, where the adiabatic index increases. A strong accretion shock forms, which propagates towards the protoneutron star surface. Due to the density decrease of several orders of magnitude, the accretion shock turns into a dynamic shock with matter outflow. This moment defines the onset of the explosion in supernova models that allow for a quark-hadron phase transition, where otherwise no explosions could be obtained. The shock propagation across the neutrinospheres releases a burst of neutrinos. This serves as a strong observable identification for the structural reconfiguration of the stellar core. The ejected matter expands on a short timescale and remains neutron-rich. These conditions might be suitable for the production of heavy elements via the r-process. The neutron-rich material is followed by proton-rich neutrino-driven ejecta in the later cooling phase of the protoneutron star where the vp-process might occur.
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.
We trained deep neural networks (DNNs) as a function of the neutrino energy density, flux, and the fluid velocity to reproduce the Eddington tensor for neutrinos obtained in our first-principles core-collapse supernova (CCSN) simulations. Although the moment method, which is one of the most popular approximations for neutrino transport, requires a closure relation, none of the analytical closure relations commonly employed in the literature captures all aspects of the neutrino angular distribution in momentum space. In this paper, we developed a closure relation by using the DNN that takes the neutrino energy density, flux, and the fluid velocity as the input and the Eddington tensor as the output. We consider two kinds of DNNs: a conventional DNN named a component-wise neural network (CWNN) and a tensor-basis neural network (TBNN). We found that the diagonal component of the Eddington tensor is reproduced better by the DNNs than the M1-closure relation especially for low to intermediate energies. For the off-diagonal component, the DNNs agree better with the Boltzmann solver than the M1 closure at large radii. In the comparison between the two DNNs, the TBNN has slightly better performance than the CWNN. With the new closure relations at hand based on the DNNs that well reproduce the Eddington tensor with much smaller costs, we opened up a new possibility for the moment method.
We explore the abundance of light clusters in core-collapse supernovae at post-bounce stage in a quantum statistical approach. Adopting the profile of a supernova core from detailed numerical simulations, we study the distribution of light bound clusters up to alpha particles (A=2-4) as well as heavy nuclei (A > 4) in dense matter at finite temperature. Within the frame of a cluster-mean field approach, the abundances of light clusters are evaluated accounting for self-energy, Pauli blocking and effects of continuum correlations. We find that deuterons and tritons, in addition to 3He and 4He, appear abundantly in a wide region from the surface of the proto-neutron star to the position of the shock wave. The appearance of light clusters may modify the neutrino emission in the cooling region and the neutrino absorption in the heating region, and, thereby, influence the supernova mechanism.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
