Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userModeling core collapse supernovae in 2 and 3 dimensions with spectral neutrino transport
101
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
The overwhelming evidence that the core collapse supernova mechanism is inherently multidimensional, the complexity of the physical processes involved, and the increasing evidence from simulations that the explosion is marginal presents great computational challenges for the realistic modeling of this event, particularly in 3 spatial dimensions. We have developed a code which is scalable to computations in 3 dimensions which couples PPM Lagrangian with remap hydrodynamics [1], multigroup, flux-limited diffusion neutrino transport [2], with many improvements), and a nuclear network [3]. The neutrino transport is performed in a ray-by-ray plus approximation wherein all the lateral effects of neutrinos are included (e.g., pressure, velocity corrections, advection) except the transport. A moving radial grid option permits the evolution to be carried out from initial core collapse with only modest demands on the number of radial zones. The inner part of the core is evolved after collapse along with the rest of the core and mantle by subcycling the lateral evolution near the center as demanded by the small Courant times. We present results of 2-D simulations of a symmetric and an asymmetric collapse of both a 15 and an 11 M progenitor. In each of these simulations we have discovered that once the oxygen rich material reaches the shock there is a synergistic interplay between the reduced ram pressure, the energy released by the burning of the shock heated oxygen rich material, and the neutrino energy deposition which leads to a revival of the shock and an explosion.
rate research
Read More
The proposal that core collapse supernovae are neutrino driven is still the subject of active investigation more than fifty years after the seminal paper by Colgate and White. The modern version of this paradigm, which we owe to Wilson, proposes that the supernova shock wave is powered by neutrino heating, mediated by the absorption of electron-flavor neutrinos and antineutrinos emanating from the proto-neutron star surface, or neutrinosphere. Neutrino weak interactions with the stellar core fluid, the theory of which is still evolving, are flavor and energy dependent. The associated neutrino mean free paths extend over many orders of magnitude and are never always small relative to the stellar core radius. Thus, neutrinos are never always fluid like. Instead, a kinetic description of them in terms of distribution functions that determine the number density of neutrinos in the six-dimensional phase space of position, direction, and energy, for both neutrinos and antineutrinos of each flavor, or in terms of angular moments of these neutrino distributions that instead provide neutrino number densities in the four-dimensional phase-space subspace of position and energy, is needed. In turn, the computational challenge is twofold: (i) to map the kinetic equations governing the evolution of these distributions or moments onto discrete representations that are stable, accurate, and, perhaps most important, respect physical laws such as conservation of lepton number and energy and the Fermi--Dirac nature of neutrinos and (ii) to develop efficient, supercomputer-architecture-aware solution methods for the resultant nonlinear algebraic equations. In this review, we present the current state of the art in attempts to meet this challenge.
We have conducted a series of numerical experiments with the spherically symmetric, general relativistic, neutrino radiation hydrodynamics code Agile-BOLTZTRAN to examine the effects of several approximations used in multidimensional core-collapse supernova simulations. Our code permits us to examine the effects of these approximations quantitatively by removing, or substituting for, the pieces of supernova physics of interest. These approximations include: (1) using Newtonian versus general relativistic gravity, hydrodynamics, and transport; (2) using a reduced set of weak interactions, including the omission of non-isoenergetic neutrino scattering, versus the current state-of-the-art; and (3) omitting the velocity-dependent terms, or observer corrections, from the neutrino Boltzmann kinetic equation. We demonstrate that each of these changes has noticeable effects on the outcomes of our simulations. Of these, we find that the omission of observer corrections is particularly detrimental to the potential for neutrino-driven explosions and exhibits a failure to conserve lepton number. Finally, we discuss the impact of these results on our understanding of current, and the requirements for future, multidimensional models.
We develop a numerical code to calculate the neutrino transfer with multi-energy and multi-angle in three dimensions (3D) for the study of core-collapse supernovae. The numerical code solves the Boltzmann equations for neutrino distributions by the discrete-ordinate (S_n) method with a fully implicit differencing for time advance. The Boltzmann equations are formulated in the inertial frame with collision terms being evaluated to the zeroth order of v/c. A basic set of neutrino reactions for three neutrino species is implemented together with a realistic equation of state of dense matter. The pair process is included approximately in order to keep the system linear. We present numerical results for a set of test problems to demonstrate the ability of the code. The numerical treatments of advection and collision terms are validated first in the diffusion and free streaming limits. Then we compute steady neutrino distributions for a background extracted from a spherically symmetric, general relativistic simulation of 15Msun star and compare them with the results in the latter computation. We also demonstrate multi-D capabilities of the 3D code solving neutrino transfers for artificially deformed supernova cores in 2D and 3D. Formal solutions along neutrino paths are utilized as exact solutions. We plan to apply this code to the 3D neutrino-radiation hydrodynamics simulations of supernovae. This is the first article in a series of reports on the development.
We have made core-collapse supernova simulations that allow oscillations between electron neutrinos (or their anti particles) with right-handed sterile neutrinos. We have considered a range of mixing angles and sterile neutrino masses including those consistent with sterile neutrinos as a dark matter candidate. We examine whether such oscillations can impact the core bounce and shock reheating in supernovae. We identify the optimum ranges of mixing angles and masses that can dramatically enhance the supernova explosion by efficiently transporting electron anti-neutrinos from the core to behind the shock where they provide additional heating leading to much larger explosion kinetic energies. We show that this effect can cause stars to explode that otherwise would have collapsed. We find that an interesting periodicity in the neutrino luminosity develops due to a cycle of depletion of the neutrino density by conversion to sterile neutrinos that shuts off the conversion, followed by a replenished neutrino density as neutrinos transport through the core.
We have explored the impact of sterile neutrino dark matter on core-collapse supernova explosions. We have included oscillations between electron neutrinos or mixed $mu,tau$ neutrinos and right-handed sterile neutrinos into a supernova model. We have chosen sterile neutrino masses and mixing angles that are consistent with sterile neutrino dark matter candidates as indicated by recent x-ray flux measurements. Using these simulations, we have explored the impact of sterile neutrinos on the core bounce and shock reheating. We find that, for ranges of sterile neutrino mass and mixing angle consistent with most dark matter constraints, the shock energy can be significantly enhanced and even a model that does not explode can be made to explode. In addition, we have found that the presence of a sterile neutrino may lead to detectable changes in the observed neutrino luminosities.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
