Source author record

John M. Blondin

John M. Blondin 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

15works
5topics
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

15 published item(s)

preprint2021arXiv

Type Ia Supernova Models: Asymmetric Remnants and Supernova Remnant G1.9+0.3

The youngest Galactic supernova remnant G1.9+0.3, probably the result of a Type Ia supernova, shows surprising anomalies in the distribution of its ejecta in space and velocity. In particular, high-velocity shocked iron is seen in several locations far from the remnant center, in some cases beyond prominent silicon and sulfur emission. These asymmetries strongly suggest a highly asymmetric explosion. We present high-resolution hydrodynamic simulations in two and three dimensions of the evolution from ages of 100 seconds to hundreds of years of two asymmetric Type Ia models, expanding into a uniform medium. At the age of G1.9+0.3 (about 100 years), our 2D model shows almost no iron shocked to become visible in X-rays. Only in a much higher-density environment could significant iron be shocked, at which time the model's expansion speed is completely inconsistent with the observations of G1.9+0.3. Our 3D model, evolving the most asymmetric of a suite of Type Ia SN models from Seitenzahl et al.~(2013), shows some features resembling G1.9+0.3. We characterize its evolution with images of composition in three classes: C and O, intermediate-mass elements (IMEs), and iron-group elements (IGEs). From ages of 13 to 1800 years, we follow the evolution of the highly asymmetric initial remnant as the explosion asymmetries decrease in relative strength to be replaced by asymmetries due to evolutionary hydrodynamic instabilities. At an age of about 100 years, our 3D model has comparable shocked masses of C+O, IMEs, and IGEs, with about 0.03 $M_\odot$ each. Evolutionary changes appear to be rapid enough that continued monitoring with the Chandra X-ray Observatory may show significant variations.

preprint2020arXiv

Chimera: A massively parallel code for core-collapse supernova simulation

We provide a detailed description of the Chimera code, a code developed to model core collapse supernovae in multiple spatial dimensions. The core collapse supernova explosion mechanism remains the subject of intense research. Progress to date demonstrates that it involves a complex interplay of neutrino production, transport, and interaction in the stellar core, three-dimensional stellar core fluid dynamics and its associated instabilities, nuclear burning, and the foundational physics of the neutrino-stellar core weak interactions and the equations of state of all stellar core constituents -particularly, the nuclear equation of state associated with nucleons, both free and bound in nuclei. Chimera, by incorporating detailed neutrino transport, realistic neutrino-matter interactions, three-dimensional hydrodynamics, realistic nuclear, leptonic, and photonic equations of state, and a nuclear reaction network, along with other refinements, can be used to study the role of neutrino radiation, hydrodynamic instabilities, and a variety of input physics in the explosion mechanism itself. It can also be used to compute observables such as neutrino signatures, gravitational radiation, and the products of nucleosynthesis associated with core collapse supernovae. The code contains modules for neutrino transport, multidimensional compressible hydrodynamics, nuclear reactions, a variety of neutrino interactions, equations of state, and modules to provide data for post-processing observables such as the products of nucleosynthesis, and gravitational radiation. Chimera is an evolving code, being updated periodically with improved input physics and numerical refinements. We detail here the current version of the code, from which future improvements will stem, which can in turn be described as needed in future publications.

preprint2020arXiv

Gravitational-wave signal of a core-collapse supernova explosion of a 15 Solar mass star

We report on the gravitational wave signal computed in the context of a three-dimensional simulation of a core collapse supernova explosion of a 15 Solar mass star. The simulation was performed with our neutrino hydrodynamics code Chimera. We detail the gravitational wave strains as a function of time, for both polarizations, and discuss their physical origins. We also present the corresponding spectral signatures. Gravitational wave emission in our model has two key features: low-frequency emission (< 200 Hz) emanates from the gain layer as a result of neutrino-driven convection and the SASI and high-frequency emission (> 600 Hz) emanates from the proto-neutron star due to Ledoux convection within it. The high-frequency emission dominates the gravitational wave emission in our model and emanates largely from the convective layer itself, not from the convectively stable layer above it, due to convective overshoot. Moreover, the low-frequency emission emanates from the gain layer itself, not from the proto-neutron star, due to accretion onto it. We provide evidence of the SASI in our model and demonstrate that the peak of our low-frequency gravitational wave emission spectrum corresponds to it. Given its origin in the gain layer, we classify the SASI emission in our model as p-mode emission and assign a purely acoustic origin, not a vortical-acoustic origin, to it. Our dominant proto-neutron star gravitational wave emission is not well characterized by emission from surface g-modes, complicating the relationship between peak frequencies observed and the mass and radius of the proto-neutron star expressed by analytic estimates under the assumption of surface g-mode emission. We present our frequency normalized characteristic strain along with the sensitivity curves of current- and next-generation gravitational wave detectors.

preprint2016arXiv

An X-ray and Radio Study of the Varying Expansion Velocities in Tycho's Supernova Remnant

We present newly obtained X-ray and radio observations of Tycho's supernova remnant using {\it Chandra} and the Karl G. Jansky Very Large Array in 2015 and 2013/14, respectively. When combined with earlier epoch observations by these instruments, we now have time baselines for expansion measurements of the remnant of 12-15 year in the X-rays and 30 year in the radio. The remnant's large angular size allows for proper motion measurements at many locations around the periphery of the blast wave. We find, consistent with earlier measurements, a clear gradient in the expansion velocity of the remnant, despite its round shape. The proper motions on the western and southwestern sides of the remnant are about a factor of two higher than those in the east and northeast. We showed in an earlier work that this is related to an offset of the explosion site from the geometric center of the remnant due to a density gradient in the ISM, and using our refined measurements reported here, we find that this offset is $\sim 23"$ towards the northeast. An explosion center offset in such a circular remnant has implications for searches for progenitor companions in other remnants.

preprint2015arXiv

A Neutrino-Driven Core Collapse Supernova Explosion of a 15 M Star

We present results from an ab initio three-dimensional, multi-physics core collapse supernova simulation for the case of a 15 M progenitor. Our simulation includes multi-frequency neutrino transport with state-of-the-art neutrino interactions in the "ray-by-ray" approximation, and approximate general relativity. Our model exhibits a neutrino-driven explosion. The shock radius begins an outward trajectory at approximately 275 ms after bounce, giving the first indication of a developing explosion in the model. The onset of this shock expansion is delayed relative to our two-dimensional counterpart model, which begins at approximately 200 ms after core bounce. At a time of 441 ms after bounce, the angle-averaged shock radius in our three-dimensional model has reached 751 km. Further quantitative analysis of the outcomes in this model must await further development of the post-bounce dynamics and a simulation that will extend well beyond 1 s after stellar core bounce, based on the results for the same progenitor in the context of our two-dimensional, counterpart model. This more complete analysis will determine whether or not the explosion is robust and whether or not observables such as the explosion energy, 56Ni mass, etc. are in agreement with observations. Nonetheless, the onset of explosion in our ab initio three-dimensional multi-physics model with multi-frequency neutrino transport and general relativity is encouraging.

preprint2015arXiv

Gravitational Wave Signatures of Ab Initio Two-Dimensional Core Collapse Supernova Explosion Models for 12-25 Solar Masses Stars

We present the gravitational waveforms computed in ab initio two-dimensional core collapse supernova models evolved with the Chimera code for progenitor masses between 12 and 25 solar masses. All models employ multi-frequency neutrino transport in the ray-by-ray approximation, state-of-the-art weak interaction physics, relativistic transport corrections such as the gravitational redshift of neutrinos, two-dimensional hydrodynamics with the commensurate relativistic corrections, Newtonian self-gravity with a general relativistic monopole correction, and the Lattimer-Swesty equation of state with 220 MeV compressibility, and begin with the most recent Woosley-Heger nonrotating progenitors in this mass range. All of our models exhibit robust explosions. Therefore, our waveforms capture all stages of supernova development: 1) a relatively short and weak prompt signal, 2) a quiescent stage, 3) a strong signal due to convection and SASI activity, 4) termination of active accretion onto the proto-neutron star, and 5) a slowly increasing tail that reaches a saturation value. Fourier decomposition shows that the gravitational wave signals we predict should be observable by AdvLIGO for Galactic events across the range of progenitors considered here. The fundamental limitation of these models is in their imposition of axisymmetry. Further progress will require counterpart three-dimensional models, which are underway.

preprint2015arXiv

Recent Progress on Ascertaining the Core Collapse Supernova Explosion Mechanism

We have been working within the fundamental paradigm that core collapse supernovae (CCSNe) may be neutrino driven, since the first suggestion of this by Colgate and White nearly five decades ago. Computational models have become increasingly sophisticated, first in one spatial dimension assuming spherical symmetry, then in two spatial dimensions assuming axisymmetry, and now in three spatial dimensions with no imposed symmetries. The increase in the number of spatial dimensions has been accompanied by an increase in the physics included in the models, and an increase in the sophistication with which this physics has been modeled. Computation has played an essential role in the development of CCSN theory, not simply for the obvious reason that such multidimensional, multi-physics, nonlinear events cannot possibly be fully captured analytically, but for its role in discovery. In particular, the discovery of the standing accretion shock instability (SASI) through computation about a decade ago has impacted all simulations performed since then. Today, we appear to be at a threshold, where neutrinos, neutrino-driven convection, and the SASI, working together over time scales significantly longer than had been anticipated in the past, are able to generate explosions, and in some cases, robust explosions, in a number of axisymmetric models. But how will this play out in three dimensions? Early results from the first three-dimensional (3D), multi-physics simulation of the "Oak Ridge" group are promising. I will discuss the essential components of today's models and the requirements of realistic CCSN modeling, present results from our one-, two-, and three-dimensional models, place our models in context with respect to other efforts around the world, and discuss short- and long-term next steps.

preprint2015arXiv

The Development of Explosions in Axisymmetric Ab Initio Core-Collapse Supernova Simulations of 12-25 $M_\odot$ Stars

We present four ab initio axisymmetric core-collapse supernova simulations for 12, 15, 20, and 25 $M_\odot$ progenitors. All of the simulations yield explosions and have been evolved for at least 1.2 seconds after core bounce and 1 second after material first becomes unbound. Simulations were computed with our Chimera code employing spectral neutrino transport, special and general relativistic transport effects, and state-of-the-art neutrino interactions. Continuing the evolution beyond 1 second allows explosions to develop more fully and the processes powering the explosions to become more clearly evident. We compute explosion energy estimates, including the binding energy of the stellar envelope outside the shock, of 0.34, 0.88, 0.38, and 0.70 B ($10^{51}$ ergs) and increasing at 0.03, 0.15, 0.19, and 0.52 B s$^{-1}$, respectively, for the 12, 15, 20, and 25 $M_\odot$ models. Three models developed pronounced prolate shock morphologies, while the 20 $M_\odot$ model, though exhibiting lobes and accretion streams like the other models, develops an approximately spherical, off-center shock as the explosion begins and then becomes moderately prolate $\sim$600 ms after bounce. This reduces the explosion energy relative to the other models by reducing mass accretion during the critical explosion power-up phase. We examine the growth of the explosion energy in our models through detailed analyses of the energy sources and flows. We find that the 12 and 20 $M_\odot$ models have explosion energies comparable to that of the lower range of observed explosion energies while the 15 and 25 $M_\odot$ models are within the range of observed explosion energies, particularly considering the rate at which their explosion energies are increasing. The ejected $^{56}$Ni masses given by our models are all within observational limits as are the proto-neutron star masses and kick velocities. (Truncated)

preprint2015arXiv

Three-dimensional core-collapse supernova simulated using a 15 $M_\odot$ progenitor

We have performed ab initio neutrino radiation hydrodynamics simulations in three and two spatial dimensions (3D and 2D) of core-collapse supernovae from the same 15 $M_\odot$ progenitor through 440 ms after core bounce. Both 3D and 2D models achieve explosions, however, the onset of explosion (shock revival) is delayed by $\sim$100 ms in 3D relative to the 2D counterpart and the growth of the diagnostic explosion energy is slower. This is consistent with previously reported 3D simulations utilizing iron-core progenitors with dense mantles. In the $\sim$100 ms before the onset of explosion, diagnostics of neutrino heating and turbulent kinetic energy favor earlier explosion in 2D. During the delay, the angular scale of convective plumes reaching the shock surface grows and explosion in 3D is ultimately lead by a single, large-angle plume, giving the expanding shock a directional orientation not dissimilar from those imposed by axial symmetry in 2D simulations. We posit that shock revival and explosion in the 3D simulation may be delayed until sufficiently large plumes form, whereas such plumes form more rapidly in 2D, permitting earlier explosions.

preprint2014arXiv

Two- and Three-Dimensional Multi-Physics Simulations of Core Collapse Supernovae: A Brief Status Report and Summary of Results from the "Oak Ridge" Group

We summarize the results of core collapse supernova theory from one-, two-, and three-dimensional models and provide a snapshot of the field at this time. We also present results from the "Oak Ridge" group in this context. Studies in both one and two spatial dimensions define the necessary} physics that must be included in core collapse supernova models: a general relativistic treatment of gravity (at least an approximate one), spectral neutrino transport, including relativistic effects such as gravitational redshift, and a complete set of neutrino weak interactions that includes state-of-the-art electron capture on nuclei and energy-exchanging scattering on electrons and nucleons. Whether or not the necessarily approximate treatment of this physics in current models that include it is sufficient remains to be determined in the context of future models that remove the approximations. We summarize the results of the Oak Ridge group's two-dimensional supernova models. In particular, we demonstrate that robust neutrino-driven explosions can be obtained. We also demonstrate that our predictions of the explosion energies and remnant neutron star masses are in agreement with observations, although a much larger number of models must be developed before more confident conclusions can be made. We provide preliminary results from our ongoing three dimensional model with the same physics. Finally, we speculate on future outcomes and directions.

preprint2013arXiv

Azimuthal Density Variations Around the Rim of Tycho's Supernova Remnant

{\it Spitzer} images of Tycho's supernova remnant in the mid-infrared reveal limb-brightened emission from the entire periphery of the shell and faint filamentary structures in the interior. As with other young remnants, this emission is produced by dust grains, warmed to $\sim 100$ K in the post-shock environment by collisions with energetic electrons and ions. The ratio of the 70 to 24 $μ$m fluxes is a diagnostic of the dust temperature, which in turn is a sensitive function of the plasma density. We find significant variations in the 70/24 flux ratio around the periphery of Tycho's forward shock, implying order-of-magnitude variations in density. While some of these are likely localized interactions with dense clumps of the interstellar medium, we find an overall gradient in the ambient density surrounding Tycho, with densities 3-10 times higher in the NE than in the SW. This large density gradient is qualitatively consistent with the variations in the proper motion of the shock observed in radio and X-ray studies. Overall, the mean ISM density around Tycho is quite low ($\sim 0.1-0.2$ cm$^{-3}$), consistent with the lack of thermal X-ray emission observed at the forward shock. We perform two-dimensional hydrodynamic simulations of a Type Ia SN expanding into a density gradient in the ISM, and find that the overall round shape of the remnant is still easily acheivable, even for explosions into significant gradients. However, this leads to an offset of the center of the explosion from the geometric center of the remnant of up to 20%, although lower values of 10% are preferred. The best match with hydrodynamical simulations is achieved if Tycho is located at a large (3-4 kpc) distance in a medium with a mean preshock density of $\sim 0.2$ cm$^{-3}$. Such preshock densities are obtained for highly ($\ga 50$%) porous ISM grains.

preprint2012arXiv

Hoyle-Lyttleton Accretion in Three Dimensions

We investigate the stability of gravitational accretion of an ideal gas onto a compact object moving through a uniform medium at Mach 3. Previous three-dimensional simulations have shown that such accretion is not stable, and that strong rotational 'disk-like' flows are generated and accreted on short time scales. We re-address this problem using overset spherical grids that provide a factor of seven improvement in spatial resolution over previous simulations. With our higher spatial resolution we found these 3D accretion flows remained remarkably axisymmetric. We examined two cases of accretion with different sized accretors. The larger accretor produced very steady flow, with the mass accretion rate varying by less than 0.02% over 30 flow times. The smaller accretor exhibited an axisymmetric breathing mode that modulated the mass accretion rate by a constant 20%. Nonetheless, the flow remained highly axisymmetric with only negligible accretion of angular momentum in both cases.

preprint2012arXiv

Turbulent Magnetic Field Amplification from Spiral SASI Modes: Implications for Core-Collapse Supernovae and Proto-Neutron Star Magnetization

We extend our investigation of magnetic field evolution in three-dimensional flows driven by the stationary accretion shock instability (SASI) with a suite of higher-resolution idealized models of the post-bounce core-collapse supernova environment. Our magnetohydrodynamic simulations vary in initial magnetic field strength, rotation rate, and grid resolution. Vigorous SASI-driven turbulence inside the shock amplifies magnetic fields exponentially; but while the amplified fields reduce the kinetic energy of small-scale flows, they do not seem to affect the global shock dynamics. The growth rate and final magnitude of the magnetic energy are very sensitive to grid resolution, and both are underestimated by the simulations. Nevertheless our simulations suggest that neutron star magnetic fields exceeding $10^{14}$ G can result from dynamics driven by the SASI, \emph{even for non-rotating progenitors}.

preprint2012arXiv

X-ray Emission from Strongly Asymmetric Circumstellar Material in the Remnant of Kepler's Supernova

Kepler's supernova remnant resulted from a thermonuclear explosion, but is interacting with circumstellar material (CSM) lost from the progenitor system. We describe a statistical technique for isolating X-ray emission due to CSM from that due to shocked ejecta. Shocked CSM coincides well in position with 24 $μ$m emission seen by {\sl Spitzer}. We find most CSM to be distributed along the bright north rim, but substantial concentrations are also found projected against the center of the remnant, roughly along a diameter with position angle $\sim 100^\circ$. We interpret this as evidence for a disk distribution of CSM before the SN, with the line of sight to the observer roughly in the disk plane. We present 2-D hydrodynamic simulations of this scenario, in qualitative agreement with the observed CSM morphology. Our observations require Kepler to have originated in a close binary system with an AGB star companion.

preprint2011arXiv

RCW 86: A Type Ia Supernova in a Wind-Blown Bubble

We report results from a multi-wavelength analysis of the Galactic SNR RCW 86, the proposed remnant of the supernova of 185 A.D. We report new infrared observations from {\it Spitzer} and {\it WISE}, where the entire shell is detected at 24 and 22 $μ$m. We fit the infrared flux ratios with models of collisionally heated ambient dust, finding post-shock gas densities in the non-radiative shocks of 2.4 and 2.0 cm$^{-3}$ in the SW and NW portions of the remnant, respectively. The Balmer-dominated shocks around the periphery of the shell, large amount of iron in the X-ray emitting ejecta, and lack of a compact remnant support a Type Ia origin for this remnant. From hydrodynamic simulations, the observed characteristics of RCW 86 are successfully reproduced by an off-center explosion in a low-density cavity carved by the progenitor system. This would make RCW 86 the first known case of a Type Ia supernova in a wind-blown bubble. The fast shocks ($> 3000$ km s$^{-1}$) observed in the NE are propagating in the low-density bubble, where the shock is just beginning to encounter the shell, while the slower shocks elsewhere have already encountered the bubble wall. The diffuse nature of the synchrotron emission in the SW and NW is due to electrons that were accelerated early in the lifetime of the remnant, when the shock was still in the bubble. Electrons in a bubble could produce gamma-rays by inverse-Compton scattering. The wind-blown bubble scenario requires a single-degenerate progenitor, which should leave behind a companion star.