Supernova Blast Wave Research Articles

The Role of Inhomogeneities in Supernova Shock Breakout Emission

Abstract The breakout of a supernova blast wave from its progenitor star provides strong constraints on the star and its immediate surroundings. These surroundings are shaped by mass loss from the star and can include a wide variety of inhomogeneities. Here we present results of multidimensional radiation-hydrodynamics calculations of the interactions of the supernova blast wave with inhomogeneities in the immediate surroundings of a massive Wolf–Rayet star, calculating the effect these interactions have on the shock breakout signal from supernovae.

Supernovae within Pre-existing Wind-blown Bubbles: Dust Injection versus Ambient Dust Destruction

Abstract By means of 3D hydrodynamical simulations, we evaluate here the impact that supernova (SN) explosions occurring within wind-driven bubbles have on the survival or destruction of dust grains. We consider both the dust generated within the ejecta and the dust initially present in the ambient gas and later locked up in the surrounding wind-driven shell (WDS). The collision of the SN blast wave with the WDS leads to a transmitted shock that moves into the shell and a reflected shock that moves into the ejecta. The transmitted shock is capable of destroying large amounts of the dust locked in the shell, but only if the mass of the WDS is small, less than a few tens the ejected mass. Conversely, massive WDSs, with several times the ejected mass, lead upon the interaction to strong radiative cooling, which inhibits the Sedov–Taylor phase and weakens the transmitted shock, making it unable to traverse the WDS. In such a case, the destruction/disruption of the ambient dust is largely inhibited. On the other hand, the SN remnants grow rapidly in the very tenuous region excavated by the stellar winds, and thus a large fraction of the dust generated within the ejecta is not efficiently destroyed by the SN reverse shock, nor by the reflected shock. Our calculations favor a scenario in which core-collapse SNe within sufficiently massive WDSs supply more dust to the interstellar medium than they are able to destroy.

How do velocity structure functions trace gas dynamics in simulated molecular clouds?

Context. Supersonic disordered flows accompany the formation and evolution of molecular clouds (MCs). It has been argued that this is turbulence that can support against gravitational collapse and form hierarchical sub-structures. Aims. We examine the time evolution of simulated MCs to investigate: What physical process dominates the driving of turbulent flows? How can these flows be characterised? Are they consistent with uniform turbulence or gravitational collapse? Do the simulated flows agree with observations? Methods. We analysed three MCs that have formed self-consistently within kiloparsec-scale numerical simulations of the interstellar medium (ISM). The simulated ISM evolves under the influence of physical processes including self-gravity, stratification, magnetic fields, supernova-driven turbulence, and radiative heating and cooling. We characterise the flows using velocity structure functions (VSFs) with and without density weighting or a density cutoff, and computed in one or three dimensions. However, we do not include optical depth effects that can hide motions in the densest gas, limiting comparison of our results with observations. Results. In regions with sufficient resolution, the density-weighted VSFs initially appear to follow the expectations for uniform turbulence, with a first-order power-law exponent consistent with Larson’s size-velocity relationship. Supernova blast wave impacts on MCs produce short-lived coherent motions at large scales, increasing the scaling exponents for a crossing time. Gravitational contraction drives small-scale motions, producing scaling coefficients that drop or even turn negative as small scales become dominant. Removing the density weighting eliminates this effect as it emphasises the diffuse ISM. Conclusions. We conclude that two different effects coincidentally reproduce Larson’s size velocity relationship. Initially, uniform turbulence dominates, so the energy cascade produces VSFs that are consistent with Larson’s relationship. Later, contraction dominates and the density-weighted VSFs become much shallower or even inverted, but the relationship of the global average velocity dispersion of the MCs to their radius follows Larson’s relationship, reflecting virial equilibrium or free-fall collapse. The injection of energy by shocks is visible in the VSFs, but decays within a crossing time.

High-energy cosmic rays from compact galactic star clusters: Particle fluxes and anisotropy

It has been shown that supernova blast waves interacting with winds from massive stars in compact star clusters may be capable of producing cosmic-ray (CR) protons to above 1017 eV. We give a brief description of the colliding-shock-flow mechanism and look at generalizations of the diffusion of ∼100 PeV CRs in the turbulent galactic magnetic field present in the galactic disk. We calculate the temporal evolution of the CR anisotropy from a possible distribution of young compact massive star clusters assuming the sources are intermittent on time scales of a few million years, i.e., comparable to their residence time in the Milky Way. Within the confines of our model, we determine the galactic/extra-galactic fraction of high-energy CRs resulting in anisotropies consistent with observed values. We find that galactic star clusters may contribute a substantial fraction of ∼100 PeV CRs without producing anisotropies above observed limits.

Multiple populations in globular clusters: New insights from chemical evolution and horizontal-branch models

AbstractIn order to investigate the origin of multiple populations in globular clusters (GCs), we have constructed new chemical evolution models for proto-GCs where the supernova blast waves undergo blowout without expelling the ambient gas. Chemical enrichments in our models are then dictated by the winds of massive stars together with the asymptotic-giant-branch stars ejecta. We find that the observed Na-O anti-correlation can be reproduced when multiple episodes of starburst and enrichment are allowed to continue in proto-GCs. The “mass budget problem” is mostly resolved by our models without ad-hoc assumptions on star formation efficiency, initial mass function, and significant loss of first-generation stars. Interestingly, ages and chemical abundances predicted by this chemical evolution model are in good agreements with those independently obtained from our stellar evolution model for the horizontal-branch. We also discuss observational evidence for the GC-like multiple populations in the Milky Way bulge.

A Radio Source Coincident with the Superluminous Supernova PTF10hgi: Evidence for a Central Engine and an Analog of the Repeating FRB 121102?

Abstract We present the detection of an unresolved radio source coincident with the position of the Type I superluminous supernova (SLSN) PTF10hgi (z = 0.098) about 7.5 yr post-explosion, with a flux density of F ν (6 GHz) ≈ 47.3 μJy and a luminosity of L ν (6 GHz) ≈ 1.1 × 1028 erg s−1 Hz−1. This represents the first detection of radio emission coincident with an SLSN on any timescale. We investigate various scenarios for the origin of the radio emission: star formation activity, an active galactic nucleus, and a non-relativistic supernova blastwave. While any of these would be quite novel if confirmed, none appear likely when considered within the context of the other properties of the host galaxy, previous radio observations of SLSNe, and the general population of hydrogen-poor supernovae (SNe). Instead, the radio emission is reminiscent of the quiescent radio source associated with the repeating FRB 121102, which has been argued to be powered by a magnetar born in a SLSN or long gamma-ray burst explosion several decades ago. We show that the properties of the radio source are consistent with a magnetar wind nebula or an off-axis jet, indicating the presence of a central engine. Our directed search for fast radio bursts from the location of PTF10hgi using 40 minutes of Very Large Array phased-array data reveals no detections to a limit of 22 mJy (10σ; 10 ms duration). We outline several follow-up observations that can conclusively establish the origin of the radio emission.

Bell-instability-mediated Spectral Modulation of Hadronic Gamma-Rays from a Supernova Remnant Interacting with a Molecular Cloud

Abstract Supernova remnants (SNRs) are believed to be the site of galactic cosmic-ray acceleration. However, the details of the cosmic-ray acceleration are still not well understood. Gamma-ray observation is a promising method to study cosmic-ray acceleration in the SNRs, because a hadronic gamma-ray can trace high-energy cosmic-rays above ∼GeV energy. Conventional theory predicts that the hadronic gamma-ray shows a flat νF ν spectrum from the pion-creation threshold energy to the maximum energy of diffusive shock acceleration. In this paper, by employing numerical simulations that solve a hybrid system of the magnetohydrodynamics of a molecular cloud and diffusive propagation of cosmic-rays, we demonstrate that the hadronic gamma-ray spectrum can be harder than the conventional one and that the modulated spectrum becomes consistent with observations. The modification mechanism is explained as follows: The cosmic-rays accelerated at the supernova blast wave shock propagate into a clump of a molecular cloud. The cosmic-ray streaming in the cloud induces the so-called Bell instability that induces Alfvén waves in the cloud. The induced magnetic field fluctuations prevent further cosmic-ray incursion by diminishing the diffusion coefficient for the cosmic-rays below ∼1 TeV energy. This mechanism reinforces recent claims of a similar spectral modification by magnetic field amplification around a molecular cloud by Inoue et al. and Celli et al.

Explaining the Multiple Populations in Globular Clusters by Multiple Episodes of Star Formation and Enrichment without Gas Expulsion from Massive Star Feedback

Abstract In order to investigate the origin of multiple stellar populations found in globular clusters (GCs) in the halo and bulge of the Milky Way, we have constructed chemical evolution models for their putative low-mass progenitors. In light of recent theoretical developments, we assume that supernova blast waves undergo blowout without expelling the pre-enriched ambient gas, while relatively slow winds of massive stars (WMSs), together with the winds and ejecta from low- to high-mass asymptotic giant branch stars, are all locally retained in these less massive systems. Interestingly, we find that the observed Na–O anti-correlations in metal-poor GCs can be reproduced when multiple episodes of starburst and enrichment are allowed to continue in these subsystems. A specific form of star formation history with decreasing time intervals between the successive stellar generations, however, is required to obtain this result, which is in good agreement with the parameters obtained from synthetic horizontal branch models. The “mass budget problem” is also much alleviated by our models without ad hoc assumptions on star formation efficiency, initial mass function, and the preferential loss of first-generation stars. We also apply these models to investigate the origin of super-He-rich red clump stars in the metal-rich bulge suggested by Lee et al. We find that chemical enrichment by the WMSs can naturally reproduce the required strong He enhancement in metal-rich subsystems. Our results further underscore that gas expulsion or retention is a key factor in understanding the multiple populations in GCs.

SOFIA mid-infrared observations of Supernova 1987A in 2016 – forward shocks and possible dust re-formation in the post-shocked region

The equatorial ring of Supernova (SN) 1987A has been exposed to forward shocks from the SN blast wave, and it has been suggested that these forward shocks have been causing on-going destruction of dust in the ring. We obtained SOFIA FORCAST 11.1, 19.7 and 31.5 micron photometry of SN\,1987A in 2016. Compared with Spitzer measurements 10 years earlier, the 31.5 micron flux has significantly increased. The excess at 31.5 micron appears to be related to the Herschel 70 micron excess, which was detected 5 years earlier. The dust mass needed to account for the the 31.5--70 micron excess is 3--7x10^-4 Msun, more than ten times larger than the ring dust mass (1x10^-5 Msun) estimate from the data 10-years earlier. We argue that dust grains are re-formed or grown in the post-shock regions in the ring after forward shocks have destroyed pre-existing dust grains in the ring and released refractory elements into gas. In the post-shock region, atoms can stick to surviving dust grains, and the dust mass may have increased (grain growth), or dust grains might have condensed directly from the gas. An alternative possibility is that the outer part of the expanding ejecta dust might have been heated by X-ray emission from the circumstellar ring. The future development of this excess could reveal whether grains are reformed in the post-shocked region of the ring or eject dust is heated by X-ray.

The consequences of a nearby supernova on the early solar system

If the Sun was born in a relatively compact open cluster (half-mass radius ≲3 pc) with ≳103stars, it is quite likely that a massive (≳10M⊙) star was nearby when it exploded in a supernova. The repercussions of a supernova can be rather profound, and the current solar system may still bear the scars of this traumatic event. The truncation of the Kuiper belt and the tilt of the ecliptic plane with respect to the Sun’s rotation axis could represent such signatures. We simulated the effect of a nearby supernova on the young solar system using the Astrophysical Multipurpose Software Environment (AMUSE). Our calculations are realized in two subsequent steps in which we study the effect of the supernova irradiation on the circumstellar disk and the effect of the impact of the nuclear blast-wave which arrives a few decades later. We find that the blast wave of our adopted supernova exploding at a distance of 0.15–0.40 pc and at an angle of 35°–65° with respect to the angular-momentum axis of the circumsolar disk would induce a misalignment between the Sun’s equator and its disk to 5.°6 ± 1.°2, consistent with the current value. The blast of a supernova truncates the disk at a radius between 42 and 55 au, which is consistent with the current edge of the Kuiper belt. For the most favored parameters, the irradiation by the supernova as well as the blast wave heat the majority of the disk to ≳1200 K, which is sufficiently hot to melt chondrules in the circumstellar disk. The majority of the planetary system may have been affected by a nearby supernova. Some of its repercussions, such as truncation and tilting of the disk, may still be visible in the current planetary system topology. The amount of material from the supernova blast wave that is accreted by the circumstellar disk is too small by several orders of magnitude to explain the current abundance of the short-lived radionuclide26Al.

Probing maximum energy of cosmic rays in SNR through gamma rays and neutrinos from the molecular clouds around SNR W28

The galactic cosmic rays are generally believed to be originated in supernova remnants (SNRs), produced in diffusive shock acceleration (DSA) process in supernova blast waves driven by expanding SNRs. One of the key unsettled issue in SNR origin of cosmic ray model is the maximum attainable energy by a cosmic ray particle in the supernova shock. Recently it has been suggested that an amplification of effective magnetic field strength at the shock may take place in young SNRs due to growth of magnetic waves induced by accelerated cosmic rays and as a result the maximum energy achieved by cosmic rays in SNR may reach the knee energy instead of ∼ 200 TeV as predicted earlier under normal magnetic field situation. In the present work we investigate the implication of such maximum energy scenarios on TeV gamma rays and neutrino fluxes from the molecular clouds interacting with the SNR W28. The authors compute the gamma-ray and neutrino flux assuming two different values for the maximum energy reached by cosmic rays in the SNR, from CR interaction in nearby molecular clouds. Both protons and nuclei are considered as accelerated particles and as target material. Our findings suggest that the issue of the maximum energy of cosmic rays in SNRs can be observationally settled by the upcoming gamma-ray experiment the Large High Altitude Air Shower Observatory (LHAASO). The estimated neutrino fluxes from the molecular clouds are , however, out of reach of the present/near future generation of neutrino telescopes.

The effect of cosmic ray acceleration on supernova blast wave dynamics

Non-relativistic shocks accelerate ions to highly relativistic energies provided that the orientation of the magnetic field is closely aligned with the shock normal (quasi-parallel shock configuration). In contrast, quasi-perpendicular shocks do not efficiently accelerate ions. We model this obliquity-dependent acceleration process in a spherically expanding blast wave setup with the moving-mesh code {\sc arepo} for different magnetic field morphologies, ranging from homogeneous to turbulent configurations. A Sedov-Taylor explosion in a homogeneous magnetic field generates an oblate ellipsoidal shock surface due to the slower propagating blast wave in the direction of the magnetic field. This is because of the efficient cosmic ray (CR) production in the quasi-parallel polar cap regions, which softens the equation of state and increases the compressibility of the post-shock gas. We find that the solution remains self-similar because the ellipticity of the propagating blast wave stays constant in time. This enables us to derive an effective ratio of specific heats for a composite of thermal gas and CRs as a function of the maximum acceleration efficiency. We finally discuss the behavior of supernova remnants expanding into a turbulent magnetic field with varying coherence lengths. For a maximum CR acceleration efficiency of about 15 per cent at quasi-parallel shocks (as suggested by kinetic plasma simulations), we find an average efficiency of about 5 per cent, independent of the assumed magnetic coherence length.

Contribution of supernovae-driven galactic fountain on fuel rate

The interstellar medium in our galaxy is a complex system which involves multiphase gas in a wide range of temperatures and densities. Accompanied by physical processes such as cooling, gravity and magnetic fields it is known that the supernovae blast waves can reshuffle the ISM by generating a galactic fountain. In our study, we mainly focus on estimating the fuel rate of the fountain which is defined as the supply rate of material from the gas falling back into the galactic disk. We take into account the effect of cooling by combining the use of non-equilibrium ionization (NEI) and collisional ionization equilibrium (CIE). Dealing with cooling in this manner, we aim to reduce computation time by using a prebuilt CIE cooling table at applicable physical conditions instead of constantly considering time dependent ion fractions for NEI cooling.

Dust in Supernovae and Supernova Remnants II: Processing and Survival

Observations have recently shown that supernovae are efficient dust factories, as predicted for a long time by theoretical models. The rapid evolution of their stellar progenitors combined with their efficiency in precipitating refractory elements from the gas phase into dust grains make supernovae the major potential suppliers of dust in the early Universe, where more conventional sources like Asymptotic Giant Branch (AGB) stars did not have time to evolve. However, dust yields inferred from observations of young supernovae or derived from models do not reflect the net amount of supernova-condensed dust able to be expelled from the remnants and reach the interstellar medium. The cavity where the dust is formed and initially resides is crossed by the high velocity reverse shock which is generated by the pressure of the circumstellar material shocked by the expanding supernova blast wave. Depending on grain composition and initial size, processing by the reverse shock may lead to substantial dust erosion and even complete destruction. The goal of this review is to present the state of the art about processing and survival of dust inside supernova remnants, in terms of theoretical modelling and comparison to observations.

The Acceleration of Charged Particles at a Spherical Shock Moving through an Irregular Magnetic Field

Abstract We investigate the physics of charged-particle acceleration at spherical shocks moving into a uniform plasma containing a turbulent magnetic field with a uniform mean. This has applications to particle acceleration at astrophysical shocks, most notably, to supernovae blast waves. We numerically integrate the equations of motion of a large number of test protons moving under the influence of electric and magnetic fields determined from a kinematically defined plasma flow associated with a radially propagating blast wave. Distribution functions are determined from the positions and velocities of the protons. The unshocked plasma contains a magnetic field with a uniform mean and an irregular component having a Kolmogorov-like power spectrum. The field inside the blast wave is determined from Maxwell’s equations. The angle between the average magnetic field and unit normal to the shock varies with position along its surface. It is quasi-perpendicular to the unit normal near the sphere’s equator, and quasi-parallel to it near the poles. We find that the highest intensities of particles, accelerated by the shock, are at the poles of the blast wave. The particles “collect” at the poles as they approximately adhere to magnetic field lines that move poleward from their initial encounter with the shock at the equator, as the shock expands. The field lines at the poles have been connected to the shock the longest. We also find that the highest-energy protons are initially accelerated near the equator or near the quasi-perpendicular portion of the shock, where the acceleration is more rapid.

Numerical studies on the link between radioisotopic signatures on Earth and the formation of the Local Bubble

The discovery of radionuclides like 60Fe with half-lives of million years in deep-sea crusts and sediments offers the unique possibility to date and locate nearby supernovae. We want to quantitatively establish that the 60Fe enhancement is the result of several supernovae which are also responsible for the formation of the Local Bubble, our Galactic habitat. We performed three-dimensional hydrodynamic adaptive mesh refinement simulations (with resolutions down to subparsec scale) of the Local Bubble and the neighbouring Loop I superbubble in different homogeneous, self-gravitating environments. For setting up the Local and Loop I superbubble, we took into account the time sequence and locations of the generating core-collapse supernova explosions, which were derived from the mass spectrum of the perished members of certain stellar moving groups. The release of 60Fe and its subsequent turbulent mixing process inside the superbubble cavities was followed via passive scalars, where the yields of the decaying radioisotope were adjusted according to recent stellar evolution calculations. The models are able to reproduce both the timing and the intensity of the 60Fe excess observed with rather high precision, provided that the external density does not exceed 0.3 cm-3 on average. Thus the two best-fit models presented here were obtained with background media mimicking the classical warm ionised and warm neutral medium. We also found that 60Fe (which is condensed onto dust grains) can be delivered to Earth via two physical mechanisms: either through individual fast-paced supernova blast waves, which cross the Earth's orbit sometimes even twice as a result of reflection from the Local Bubble's outer shell, or, alternatively, through the supershell of the Local Bubble itself, injecting the 60Fe content of all previous supernovae at once, but over a longer time range.

Variable interstellar radiation fields in simulated dwarf galaxies: supernovae versus photoelectric heating

We present high-resolution hydrodynamical simulations of isolated dwarf galaxies including self-gravity, non-equilibrium cooling and chemistry, interstellar radiation fields (IRSF) and shielding, star formation, and stellar feedback. This includes spatially and temporally varying photoelectric (PE) heating, photoionization, resolved supernova (SN) blast waves and metal enrichment. A new flexible method to sample the stellar initial mass function allows us to follow the contribution to the ISRF, the metal output and the SN delay times of individual massive stars. We find that SNe play the dominant role in regulating the global star formation rate, shaping the multi-phase interstellar medium (ISM) and driving galactic outflows. Outflow rates (with mass-loading factors of a few) and hot gas fractions of the ISM increase with the number of SNe exploding in low-density environments where radiative energy losses are low. While PE heating alone can suppress star formation slightly more (a factor of a few) than SNe alone can do, it is unable to drive outflows and reproduce the multi-phase ISM that emerges naturally when SNe are included. These results are in conflict with recent results of Forbes et al. who concluded that PE heating is the dominant process suppressing star formation in dwarfs, about an order of magnitude more efficient than SNe. Potential origins for this discrepancy are discussed. In the absence of SNe and photoionization (mechanisms to disperse dense clouds), the impact of PE heating is highly overestimated owing to the (unrealistic) proximity of dense gas to the radiation sources. This leads to a substantial boost of the infrared continuum emission from the UV-irradiated dust and a far infrared line-to-continuum ratio too low compared to observations. Though sub-dominant in regulating star formation, the ISRF controls the abundance of molecular hydrogen via photodissociation.

New insights on the origin of multiple stellar populations in globular clusters

AbstractIn order to investigate the origin of multiple stellar populations in the halo and bulge of the Milky Way, we have constructed chemical evolution models for the low-mass proto-Galactic subsystems such as globular clusters (GCs). Unlike previous studies, we assume that supernova blast waves undergo blowout without expelling the pre-enriched gas, while relatively slow winds of massive stars (WMS), together with the winds and ejecta from low and intermediate mass asymptotic-giant-branch stars (AGBs), are all locally retained in these less massive systems. We find that the observed Na-O anti-correlations in metal-poor GCs can be reproduced, when multiple episodes of starbursts are allowed to continue in these subsystems. A specific form of star formation history (SFH) with decreasing time intervals between the stellar generations, however, is required to obtain this result, which is in good agreement with the parameters obtained from our stellar evolution models for the horizontal-branch. The “mass budget problem” is also much alleviated by our models without ad-hoc assumptions on star formation efficiency (SFE) and initial mass function (IMF). We also applied these models to investigate the origin of super-helium-rich red clump stars in the metal-rich bulge as recently suggested by Lee et al. (2015). We find that chemical enrichments by the WMS can naturally reproduce the required helium enhancement (ΔY/ΔZ = 6) for the second generation stars. Disruption of proto-GCs in a hierarchical merging paradigm would have provided helium enhanced stars to the bulge field.

Implications of supernova remnant origin model of galactic cosmic rays on gamma rays from young supernova remnants

It is widely believe that galactic cosmic rays are originated in supernova remnants (SNRs) where they are accelerated by diffusive shock acceleration process at supernova blast waves driven by expanding SNRs. In recent theoretical developments of the diffusive shock acceleration theory in SNRs, protons are expected to accelerate in SNRs at least up to the knee energy. If SNRs are true generator of cosmic rays, they should accelerate not only protons but also heavier nuclei with right proportion and the maximum energy of heavier nuclei should be atomic mass (Z) times that of protons. In this work we investigate the implications of acceleration of heavier nuclei in SNRs on energetic gamma rays those are produced in hadronic interaction of cosmic rays with ambient matter. Our findings suggest that the energy conversion efficiency has to be nearly double for the mixed cosmic ray composition instead of pure protons to explain the observation and secondly the gamma ray flux above few tens of TeV would be significantly higher if cosmic rays particles can attain energies Z times of the knee energy in lieu of 200 TeV, as suggested earlier for non-amplified magnetic fields. The two stated maximum energy paradigm will be discriminated in future by the upcoming gamma ray experiments like Cherenkov Telescope array (CTA).

A “Cosmic Comb” Model of Fast Radio Bursts

Abstract Recent observations of fast radio bursts (FRBs) indicate a perplexing, inconsistent picture. We propose a unified scenario to interpret diverse FRBs observed. A regular pulsar, otherwise unnoticeable at a cosmological distance, may produce a bright FRB if its magnetosphere is suddenly “combed” by a nearby, strong plasma stream toward the anti-stream direction. If the Earth is to the night side of the stream, the combed magnetic sheath would sweep across the direction of Earth and make a detectable FRB. The stream could be an AGN flare, a GRB or supernova blastwave, a tidal disruption event, or even a stellar flare. Since it is the energy flux received by the pulsar rather than the luminosity of the stream origin that defines the properties of the FRB, this model predicts a variety of counterparts of FRBs, including a possible connection between FRB 150418 and an AGN flare, a possible connection between FRB 131104 and a weak GRB, a steady radio nebula associated with the repeating FRB 121102, and probably no bright counterparts for some FRBs.