Source author record

A. Mezzacappa

A. Mezzacappa appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

12works
9topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

12 published item(s)

preprint2016arXiv

The Multi-Dimensional Character of Core-Collapse Supernovae

Core-collapse supernovae, the culmination of massive stellar evolution, are spectacular astronomical events and the principle actors in the story of our elemental origins. Our understanding of these events, while still incomplete, centers around a neutrino-driven central engine that is highly hydrodynamically unstable. Increasingly sophisticated simulations reveal a shock that stalls for hundreds of milliseconds before reviving. Though brought back to life by neutrino heating, the development of the supernova explosion is inextricably linked to multi-dimensional fluid flows. In this paper, the outcomes of three-dimensional simulations that include sophisticated nuclear physics and spectral neutrino transport are juxtaposed to learn about the nature of the three dimensional fluid flow that shapes the explosion. Comparison is also made between the results of simulations in spherical symmetry from several groups, to give ourselves confidence in the understanding derived from this juxtaposition.

preprint2014arXiv

Planning the Future of U.S. Particle Physics (Snowmass 2013): Chapter 9: Computing

These reports present the results of the 2013 Community Summer Study of the APS Division of Particles and Fields ("Snowmass 2013") on the future program of particle physics in the U.S. Chapter 9, on Computing, discusses the computing challenges for future experiments in the Energy, Intensity, and Cosmic Frontiers, for accelerator science, and for particle theory, as well as structural issues in supporting the intense uses of computing required in all areas of particle physics.

preprint2012arXiv

Turbulence and magnetic field amplification from spiral SASI modes in core-collapse supernovae

The stationary accretion shock instability (SASI) plays a central role in modern simulations of the explosion phase of core-collapse supernovae (CCSNe). It may be key to realizing neutrino powered explosions, and possibly links birth properties of pulsars (e.g., kick, spin, and magnetic field) to supernova dynamics. Using high-resolution magnetohydrodynamic simulations, we study the development of turbulence, and subsequent amplification of magnetic fields in a simplified model of the post-bounce core-collapse supernova environment. Turbulence develops from secondary instabilities induced by the SASI. Our simulations suggest that the development of turbulence plays an important role for the subsequent evolution of the SASI. The turbulence also acts to amplify weak magnetic fields via a small-scale dynamo.

preprint2010arXiv

2D and 3D Core-Collapse Supernovae Simulation Results Obtained with the CHIMERA Code

Much progress in realistic modeling of core-collapse supernovae has occurred recently through the availability of multi-teraflop machines and the increasing sophistication of supernova codes. These improvements are enabling simulations with enough realism that the explosion mechanism, long a mystery, may soon be delineated. We briefly describe the CHIMERA code, a supernova code we have developed to simulate core-collapse supernovae in 1, 2, and 3 spatial dimensions. We then describe the results of an ongoing suite of 2D simulations initiated from a 12, 15, 20, and 25 solar mass progenitor. These have all exhibited explosions and are currently in the expanding phase with the shock at between 5,000 and 20,000 km. We also briefly describe an ongoing simulation in 3 spatial dimensions initiated from the 15 solar mass progenitor.

preprint2010arXiv

Exploring the QCD phase transition in core collapse supernova simulations in spherical symmetry

For finite chemical potential effective models of QCD predict a first order phase transition. In favour for the search of such a phase transition in nature, we construct an equation of state for strange quark matter based on the MIT bag model. We apply this equation of state to highly asymmetric core collapse supernova matter with finite temperatures and large baryon densities. The phase transition is constructed using the general Gibbs conditions, which results in an extended coexistence region between the pure hadronic and pure quark phases in the phase diagram, i.e. the mixed phase. The supernovae are simulated via general relativistic radiation hydrodynamics based on three flavor Boltzmann neutrino transport in spherical symmetry. During the dynamical evolution temperatures above 10 MeV, baryon densities above nuclear saturation density and a proton-to-baryon ratio below 0.2 are obtained. At these conditions the phase transition is triggered which leads to a significant softening of the EoS for matter in the mixed phase. As a direct consequence of the stiffening of the EoS again for matter in the pure quark phase, a shock wave forms at the boundary between the mixed and the pure hadronic phases. This shock is accelerated and propagates outward which releases a burst of neutrinos dominated by electron anti-neutrinos due to the lifted degeneracy of the shock-heated hadronic material. We discuss the radiation-hydrodynamic evolution of the phase transition at the example of several low and intermediate mass Fe-core progenitor stars and illustrate the expected neutrino signal from the phase transition.

preprint2010arXiv

Mechanisms of Core-Collapse Supernovae & Simulation Results from the CHIMERA Code

Unraveling the mechanism for core-collapse supernova explosions is an outstanding computational challenge and the problem remains essentially unsolved despite more than four decades of effort. However, much progress in realistic modeling has occurred recently through the availability of multi-teraflop machines and the increasing sophistication of supernova codes. These improvements have led to some key insights which may clarify the picture in the not too distant future. Here we briefly review the current status of the three explosion mechanisms (acoustic, MHD, and neutrino heating) that are currently under active investigation, concentrating on the neutrino heating mechanism as the one most likely responsible for producing explosions from progenitors in the mass range ~10 to ~25 solar masses. We then briefly describe the CHIMERA code, a supernova code we have developed to simulate core-collapse supernovae in 1, 2, and 3 spatial dimensions. We finally describe the results of an ongoing suite of 2D simulations initiated from a 12, 15, 20, and 25 solar mass progenitor. These have all exhibited explosions and are currently in the expanding phase with the shock at between 5,000 and 10,000 km. We finally very briefly describe an ongoing simulation in 3 spatial dimensions initiated from the 15 solar mass progenitor.

preprint2010arXiv

Protoneutron star evolution and the neutrino driven wind in general relativistic neutrino radiation hydrodynamics simulations

Massive stars end their life in an explosion event with kinetic energies of the order 1 Bethe. Immediately after the explosion has been launched, a region of low density and high entropy forms behind the ejecta which is continuously subject to neutrino heating. The neutrinos emitted from the remnant at the center, the protoneutron star (PNS), heat the material above the PNS surface. This heat is partly converted into kinetic energy and the material accelerates to an outflow that is known as the neutrino driven wind. For the first time, we simulate the collapse, bounce, explosion and the neutrino driven wind phases consistently over more than 20 seconds. Our numerical model is based on spherically symmetric general relativistic radiation hydrodynamics using spectral three flavor Boltzmann neutrino transport. In simulations where no explosions are obtained naturally, we model neutrino driven explosions for low and intermediate mass Fe-core progenitor stars by enhancing the charged current reaction rates. In the case of a special progenitor star, the O-Ne-Mg-core, the explosion in spherical symmetry was obtained without enhanced opacities. The post explosion evolution is in qualitative agreement with static steady-state and parametrized dynamic models of the neutrino driven wind. On the other hand, we find generally smaller neutrino luminosities and mean neutrino energies as well as a different evolutionary behavior of the neutrino luminosities and mean neutrino energies. The neutrino driven wind is proton-rich for more than 10 seconds and the contraction of the PNS differs from the assumptions made for the conditions at the inner boundary in previous neutrino driven wind studies. Despite the moderately large entropies per baryon of about 100 and the fast expansion timescale, the conditions found in our model are unlikely to favor...

preprint2009arXiv

The neutrino signal from protoneutron star accretion and black hole formation

We discuss the formation of stellar mass black holes via protoneutron star (PNS) collapse. In the absence of an earlier explosion, the PNS collapses to a black hole due to the continued mass accretion onto the PNS. We present an analysis of the emitted neutrino spectra of all three flavors during the PNS contraction. Special attention is given to the physical conditions which depend on the input physics, e.g. the equation of state (EoS) and the progenitor model. The PNSs are modeled as the central object in core collapse simulations using general relativistic three-flavor Boltzmann neutrino transport in spherical symmetry. The simulations are launched from several massive progenitors of 40 and 50 solar mass. We analyze the electron-neutrino luminosity dependencies and construct a simple approximation for the electron-neutrino luminosity, which depends only on the physical conditions at the electron-neutrinosphere. In addition, we analyze different mu/tau-neutrino pair-reactions separately and compare the differences during the post-bounce phase of failed core collapse supernova explosions of massive progenitors. We also investigate the connection between the increasing mu/tau-neutrino luminosity and the PNS contraction during the accretion phase before black hole formation. Comparing the different post bounce phase of the progenitor models under investigation, we find large differences in the emitted neutrino spectra. These differences and the analysis of the electron-neutrino luminosity indicate a strong progenitor model dependency of the emitted neutrino signal.

preprint2008arXiv

Magnetic field generation by the stationary accretion shock instability

By adding a weak magnetic field to a spherically symmetric fluid configuration that caricatures a stalled shock in the post-bounce supernova environment, we explore the capacity of the stationary accretion shock instability (SASI) to generate magnetic fields. The SASI develops upon perturbation of the initial condition, and the ensuing flow generates--{\em in the absence of rotation}--dynamically significant magnetic fields ($\sim 10^{15}$ G) on a time scale that is relevant for the explosion mechanism of core-collapse supernovae. We describe our model, present some recent results, and discuss their potential relevance for supernova models.

preprint2005arXiv

Composition of the Innermost Core Collapse Supernova Ejecta

With presently known input physics and computer simulations in 1D, a self-consistent treatment of core collapse supernovae does not yet lead to successful explosions, while 2D models show some promise. Thus, there are strong indications that the delayed neutrino mechanism works combined with a multi-D convection treatment for unstable layers. On the other hand there is a need to provide correct nucleosynthesis abundances for the progressing field of galactic evolution and observations of low metallicity stars. The innermost ejecta is directly affected by the explosion mechanism, i.e. most strongly the yields of Fe-group nuclei for which an induced piston or thermal bomb treatment will not provide the correct yields because the effect of neutrino interactions is not included. We apply parameterized variations to the neutrino scattering cross sections and alternatively, parameterized variations are applied to the neutrino absorption cross sections on nucleons in the ``gain region''. We find that both measures lead to similar results, causing explosions and a Ye>0.5 in the innermost ejected layers, due to the combined effect of a short weak interaction time scale and a negligible electron degeneracy, unveiling the proton-neutron mass difference. We include all weak interactions (electron and positron capture, beta-decay, neutrino and antineutrino capture on nuclei, and neutrino and antineutrino capture on nucleons) and present first nucleosynthesis results for these innermost ejected layers to discuss how they improve predictions for Fe-group nuclei. The proton-rich environment results in enhanced abundances of 45Sc, 49Ti, and 64Zn as requested by chemical evolution studies and observations of low metallicity stars as well as appreciable production of nuclei in the mass range up to A=80.

preprint2004arXiv

Supernova Simulations with Boltzmann Neutrino Transport: A Comparison of Methods

Accurate neutrino transport has been built into spherically symmetric simulations of stellar core collapse and postbounce evolution. The results of such simulations agree that spherically symmetric models with standard microphysical input fail to explode by the delayed, neutrino-driven mechanism. Independent groups implemented fundamentally different numerical methods to tackle the Boltzmann neutrino transport equation. Here we present a direct and detailed comparison of such neutrino radiation-hydrodynamical simulations for two codes, Agile-Boltztran of the Oak Ridge-Basel group and Vertex of the Garching group. The former solves the Boltzmann equation directly by an implicit, general relativistic discrete angle method on the adaptive grid of a conservative implicit hydrodynamics code with second-order TVD advection. In contrast, the latter couples a variable Eddington factor technique with an explicit, moving-grid, conservative high-order Riemann solver with important relativistic effects treated by an effective gravitational potential. The presented study is meant to test both neutrino radiation-hydrodynamics implementations and to provide a data basis for comparisons and verifications of supernova codes to be developed in the future. Results are discussed for simulations of the core collapse and post-bounce evolution of a 13 solar mass star with Newtonian gravity and a 15 solar mass star with relativistic gravity.

preprint2003arXiv

A Finite Difference Representation of Neutrino Radiation Hydrodynamics in Spherically Symmetric General Relativistic Space-Time

We present an implicit finite difference representation for general relativistic radiation hydrodynamics in spherical symmetry. Our code, Agile-Boltztran, solves the Boltzmann transport equation for the angular and spectral neutrino distribution functions in self-consistent simulations of stellar core collapse and postbounce evolution. It implements a dynamically adaptive grid in comoving coordinates. Most macroscopically interesting physical quantities are defined by expectation values of the distribution function. We optimize the finite differencing of the microscopic transport equation for a consistent evolution of important expectation values. We test our code in simulations launched from progenitor stars with 13 solar masses and 40 solar masses. ~0.5 s after core collapse and bounce, the protoneutron star in the latter case reaches its maximum mass and collapses further to form a black hole. When the hydrostatic gravitational contraction sets in, we find a transient increase in electron flavor neutrino luminosities due to a change in the accretion rate. The muon- and tauon-neutrino luminosities and rms energies, however, continue to rise because previously shock-heated material with a non-degenerate electron gas starts to replace the cool degenerate material at their production site. We demonstrate this by supplementing the concept of neutrinospheres with a more detailed statistical description of the origin of escaping neutrinos. We compare the evolution of the 13 solar mass progenitor star to simulations with the MGFLD approximation, based on a recently developed flux limiter. We find similar results in the postbounce phase and validate this MGFLD approach for the spherically symmetric case with standard input physics.