Recent far-infrared (IR) observations of supernova remnants (SNRs) have revealed significantly large amounts of newly condensed dust in their ejecta, comparable to the total mass of available ...refractory elements. The dust masses derived from these observations assume that all the grains of a given species radiate at the same temperature, regardless of the dust heating mechanism or grain radius. In this paper, we derive the dust mass in the ejecta of the Crab Nebula, using a physical model for the heating and radiation from the dust. We adopt a power-law distribution of grain sizes and two different dust compositions (silicates and amorphous carbon), and calculate the heating rate of each dust grain by the radiation from the pulsar wind nebula. We find that the grains attain a continuous range of temperatures, depending on their size and composition. The total mass derived from the best-fit models to the observed IR spectrum is 0.019-0.13 Solar Mass, depending on the assumed grain composition. We find that the power-law size distribution of dust grains is characterized by a power-law index of 3.5-4.0 and a maximum grain size larger than 0.1 micron. The grain sizes and composition are consistent with what is expected for dust grains formed in a Type IIP supernova (SN). Our derived dust mass is at least a factor of two less than the mass reported in previous studies of the Crab Nebula that assumed more simplified two-temperature models. These models also require a larger mass of refractory elements to be locked up in dust than was likely available in the ejecta. The results of this study show that a physical model resulting in a realistic distribution of dust temperatures can constrain the dust properties and affect the derived dust masses. Our study may also have important implications for deriving grain properties and mass estimates in other SNRs and for the ultimate question of whether SNe are major sources of dust in the Galactic interstellar medium and in external galaxies.
We analyze and model the infrared spectrum of the Cassiopeia A supernova remnant with the aim of determining the masses of various elements in the unshocked ejecta. In this way, we complement the ...survey of the X-ray-emitting ejecta to provide a complete census of the elemental composition of the Cas A ejecta. We calculate photoionization-recombination equilibria to determine the ionization balance of various elements in the ejecta as a function of density using the X-ray and UV emission from the forward and reverse shocks as the ionizing radiation. With the assumption that all emission lines are principally excited at the ejecta density that maximizes their emission, we can convert observed line intensities into element masses. We find that the majority of the ∼3 M ejecta have already been through the reverse shock and are seen today in X-rays. A minority, ∼0.47 0.05 M , with uncertainties quoted here coming from the data fitting procedure only, are still expanding inside the reverse shock and emitting in the infrared. This component is comprised mainly of O, Si, and S, with no Fe readily detectable. Incorporating uncertainties estimated to come from our modeling, we quote M . We speculate that up to a further 0.07 M of Fe may be present in diffuse gas in the inner ejecta, depending on the Fe charge state.
ABSTRACT
We outline a model of the Crab pulsar wind nebula with two different populations of synchrotron emitting particles, arising from two different acceleration mechanisms: (i) Component-I due to ...Fermi-I acceleration at the equatorial portion of the termination shock, with particle spectral index pI ≈ 2.2 above the injection break corresponding to γwindσwind ∼ 105, peaking in the ultraviolet (UV, γwind ∼ 102 is the bulk Lorentz factor of the wind, σwind ∼ 103 is wind magnetization); and (ii) Component-II due to acceleration at reconnection layers in the bulk of the turbulent Nebula, with particle index pII ≈ 1.6. The model requires relatively slow but highly magnetized wind. For both components, the overall cooling break is in the infrared at ∼0.01 eV, so that the Component-I is in the fast cooling regime (cooling frequency below the peak frequency). In the optical band, Component-I produces emission with the cooling spectral index of αo ≈ 0.5, softening towards the edges due to radiative losses. Above the cooling break, in the optical, UV, and X-rays, Component-I mostly overwhelms Component-II. We hypothesize that acceleration at large-scale current sheets in the turbulent nebula (Component-II) extends to the synchrotron burn-off limit of ϵs ∼ 100 MeV. Thus in our model acceleration in turbulent reconnection (Component-II) can produce both hard radio spectra and occasional gamma-ray flares. This model may be applicable to a broader class of high-energy astrophysical objects, like active galactic nuclei and gamma-ray burst jets, where often radio electrons form a different population from the high-energy electrons.
Supernova explosions are inherently asymmetric and can accelerate new-born neutron stars (NSs) to hundreds of km s−1. Two prevailing theories to explain NS kicks are ejecta asymmetries (e.g., ...conservation of momentum between NS and ejecta) and anisotropic neutrino emission. Observations of supernova remnants (SNRs) can give us insights into the mechanism that generates these NS kicks. In this paper, we investigate the relationship between NS kick velocities and the X-ray morphologies of 18 SNRs observed with the Chandra X-ray Observatory and the Röntgen Satellite (ROSAT). We measure SNR asymmetries using the power-ratio method (a multipole expansion technique), focusing on the dipole, quadrupole, and octupole power ratios. Our results show no correlation between the magnitude of the power ratios and NS kick velocities, but we find that for Cas A and G292.0+1.8, whose emission traces the ejecta distribution, their NSs are preferentially moving opposite to the bulk of the X-ray emission. In addition, we find a similar result for PKS 1209-51, CTB 109, and Puppis A; however, their emission is dominated by circumstellar/interstellar material, so their asymmetries may not reflect their ejecta distributions. Our results are consistent with the theory that NS kicks are a consequence of ejecta asymmetries as opposed to anisotropic neutrino emission. In the future, additional observations to measure NS proper motions within ejecta-dominated SNRs are necessary to robustly constrain the NS kick mechanism.
While theoretical models of dust condensation predict that most refractory elements produced in core-collapse supernovae (SNe) efficiently condense into dust, a large quantity of dust has so far only ...been observed in SN 1987A. We present an analysis of observations from the Spitzer Space Telescope, Herschel Space Observatory, Stratospheric Observatory for Infrared Astronomy, and AKARI of the infrared shell surrounding the pulsar wind nebula in the supernova remnant G54.1+0.3. We attribute a distinctive spectral feature at 21 m to a magnesium silicate grain species that has been invoked in modeling the ejecta-condensed dust in Cas A, which exhibits the same spectral signature. If this species is responsible for producing the observed spectral feature and accounts for a significant fraction of the observed infrared continuum, we find that it would be the dominant constituent of the dust in G54.1+0.3, with possible secondary contributions from other compositions, such as carbon, silicate, or alumina grains. The total mass of SN-formed dust required by this model is at least 0.3 M . We discuss how these results may be affected by varying dust grain properties and self-consistent grain heating models. The spatial distribution of the dust mass and temperature in G54.1+0.3 confirms the scenario in which the SN-formed dust has not yet been processed by the SN reverse shock and is being heated by stars belonging to a cluster in which the SN progenitor exploded. The dust mass and composition suggest a progenitor mass of 16-27 M and imply a high dust condensation efficiency, similar to that found for Cas A and SN 1987A. The study provides another example of significant dust formation in a Type IIP SN explosion and sheds light on the properties of pristine SN-condensed dust.
Abstract
The evolution of a pulsar wind nebula (PWN) depends on properties of the progenitor star, supernova, and surrounding environment. As some of these quantities are difficult to measure, ...reproducing the observed dynamical properties and spectral energy distribution (SED) with an evolutionary model is often the best approach to estimating their values. G21.5−0.9, powered by the pulsar J1833−1034, is a well observed PWN for which previous modeling efforts have struggled to reproduce the observed SED. In this study, we reanalyze archival infrared (IR; Herschel, Spitzer) and X-ray (Chandra, NuSTAR, Hitomi) observations. The similar morphology observed between IR line and continuum images of this source indicates that a significant portion of this emission is generated by surrounding dust and gas, and not synchrotron radiation from the PWN. Furthermore, we find that the broadband X-ray spectrum of this source is best described by a series of power laws fit over distinct energy bands. For all X-ray detectors, we find significant softening and decreasing unabsorbed flux in higher energy bands. Our model for the evolution of a PWN is able to reproduce the properties of this source when the supernova ejecta has a low initial kinetic energy
E
sn
≈ 1.2 × 10
50
erg and the spectrum of particles injected into the PWN at the termination shock is softer at low energies. Lastly, our hydrodynamical modeling of the supernova remnant can reproduce its morphology if there is a significant increase in the density of the ambient medium ∼1.8 pc north of the explosion center.
Abstract The recent survey of the core-collapse supernova remnant Cassiopeia A (Cas A) with the MIRI instrument on board the James Webb Space Telescope (JWST) revealed a large structure in the ...interior region, referred to as the “Green Monster.” Although its location suggests that it is an ejecta structure, the infrared properties of the “Green Monster” hint at a circumstellar medium (CSM) origin. In this companion paper to the JWST Cas A paper, we investigate the filamentary X-ray structures associated with the “Green Monster” using Chandra X-ray Observatory data. We extracted spectra along the “Green Monster” as well as from shocked CSM regions. Both the extracted spectra and a principal component analysis show that the “Green Monster” emission properties are similar to those of the shocked CSM. The spectra are well fit by a model consisting of a combination of a nonequilibrium ionization model and a power-law component, modified by Galactic absorption. All the “Green Monster” spectra show a blueshift corresponding to a radial velocity of around −2300 km s −1 , suggesting that the structure is on the near side of Cas A. The ionization age is around n e t ≈ 1.5 × 10 11 cm −3 s. This translates into a preshock density of ∼12 cm −3 , higher than previous estimates of the unshocked CSM. The relatively high n e t and relatively low radial velocity suggest that this structure has a relatively high density compared to other shocked CSM plasma. This analysis provides yet another piece of evidence that the CSM around Cas A’s progenitor was not that of a smooth steady wind profile.
Abstract
We studied the PeVatron nature of the pulsar wind nebula (PWN) G75.2+0.1 (“Dragonfly”) as part of our NuSTAR observational campaign of energetic PWNe. The Dragonfly is spatially coincident ...with LHAASO J2018+3651, whose maximum photon energy is 0.27 PeV. We detected a compact (radius
1
′
) inner nebula of the Dragonfly without a spectral break in 3–20 keV using NuSTAR. A joint analysis of the inner nebula with archival Chandra and XMM-Newton (XMM) observations yields a power-law spectrum with Γ = 1.49 ± 0.03. Synchrotron burnoff is observed from the shrinkage of the NuSTAR nebula at higher energies, from which we infer the magnetic field in the inner nebula of 24
μ
G at 3.5 kpc. Our analysis of archival XMM data and 13 yr of Fermi-LAT data confirms the detection of an extended (
∼
10
′
) outer nebula in 2–6 keV (Γ = 1.82 ± 0.03) and the nondetection of a GeV nebula, respectively. Using the VLA, XMM, and HAWC data, we modeled a multiwavelength spectral energy distribution of the Dragonfly as a leptonic PeVatron. The maximum injected particle energy of 1.4 PeV from our model suggests that the Dragonfly is likely a PeVatron. Our model prediction of the low magnetic field (2.7
μ
G) in the outer nebula and recent interaction with the host supernova remnant’s reverse shock (4 kyr ago) align with common features of PeVatron PWNe. The origin of its highly asymmetric morphology, pulsar proper motion, PWN–supernova remnant (SNR) interaction, and source distance will require further investigations in the future, including a multiwavelength study using radio, X-ray, and gamma-ray observations.
We present moderately deep (125 ks) XMM-Newton observations of supernova remnant G330.2+1.0. This remnant is one of only a few known that fall into the "synchrotron-dominated" category, with the ...emission almost entirely dominated by a nonthermal continuum. Previous X-ray observations could only characterize the spectra of a few regions. Here, we examine the spectra from 14 regions surrounding the entire rim, finding that the spectral properties of the nonthermal emission do not vary significantly in any systematic way from one part of the forward shock to another, unlike several other remnants of this class. We confirm earlier findings that the power-law index, Γ, ranges from about 2.1-2.5, while the absorbing column density is generally between (2.0-2.6) × 1022 cm−2. Fits with the srcut model find values of the roll-off frequency in the range of 1017.1-1017.5 Hz, implying energies of accelerated electrons of ∼100 TeV. These values imply a high shock velocity of ∼4600 km s−1, favoring a young age of the remnant. Diffuse emission from the interior is nonthermal in origin as well, and fits to these regions yield similar values to those along the rim, also implying a young age. Thermal emission is present in the east, and the spectrum is consistent with a ∼650 km s−1 shock wave encountering interstellar or circumstellar material with a density of ∼1 cm−3.
Abstract We present an imaging survey of the Spitzer I star-forming region in NGC 6822 conducted with the NIRCam and MIRI instruments on board JWST. Located at a distance of 490 kpc, NGC 6822 is the ...nearest non-interacting low-metallicity (∼0.2 Z ⊙ ) dwarf galaxy. It hosts some of the brightest known H ii regions in the local universe, including recently discovered sites of highly embedded active star formation. Of these, Spitzer I is the youngest and most active, and houses 90 color-selected candidate young stellar objects (YSOs) identified from Spitzer Space Telescope observations. We revisit the YSO population of Spitzer I with these new JWST observations. By analyzing color–magnitude diagrams constructed with NIRCam and MIRI data, we establish color selection criteria and construct spectral energy distributions to identify candidate YSOs and characterize the full population of young stars, from the most embedded phase to the more evolved stages. In this way, we have identified 140 YSOs in Spitzer I. Comparing to previous Spitzer studies of the NGC 6822 YSO population, we find that the YSOs we identify are fainter and less massive, indicating that the improved resolution of JWST allows us to resolve previously blended sources into multiple objects.