Structure Of Ejecta Research Articles

Planetary Nebula Morphologies Indicate a Jet-Driven Explosion of SN 1987A and Other Core-Collapse Supernovae

I demonstrate the usage of planetary nebulae (PNe) to infer that a pair of jets shaped the ejecta of the core-collapse supernova (CCSN) SN 1987A. The main structure of the SN 1987A inner ejecta, the ‘keyhole’, comprises two low-intensity zones. The northern one has a bright rim on its front, while the southern one has an elongated nozzle. An earlier comparison of the SN 1987A ‘keyhole’ with bubbles in the galaxy group NGC 5813 led to its identification as a jet-shaped rim–nozzle structure. Here, I present rim–nozzle asymmetry in planetary nebulae (PNe), thought to be shaped by jets, which solidifies the claim that jets powered the ejecta of SN 1987A and other CCSNe. This finding for the iconic SN 1987A with its unique properties strengthens the jittering-jets explosion mechanism (JJEM) of CCSNe. In a few hundred years, the CCSN 1987A will have a complicated structure with two main symmetry axes, one along the axis of the three circumstellar rings that was shaped by two opposite 20,000-year pre-explosion jets, and the other along the long axis of the ‘keyhole’ that was shaped by the main (but not the only) jet pair of the exploding jets of SN 1987A in the frame of the JJEM.

DART Impact Ejecta Plume Evolution: Implications for Dimorphos

The NASA Double Asteroid Redirection Test (DART) spacecraft impacted the moon Dimorphos of the [65803] Didymos binary system and changed the binary orbit period, demonstrating asteroid deflection by a kinetic impact and indicating that more momentum was transferred to Dimorphos by escaping impact ejecta than was incident with DART. Images of the DART impact ejecta plume were obtained by the Light Italian cubesat for Imaging of Asteroids (LICIACube) in the first few minutes after the DART impact. The ejecta plume imaged by LICIACube 158 s after the DART impact prior to closest approach shows no evidence for plume clearing at low altitude. The ejecta plume imaged 175 s after the DART impact is optically thick up to projected altitudes of 200 m above the surface of Dimorphos. These observations are compared with models of the impact ejecta plume optical depth, structure, and evolution, which are developed from point-source scaling models fitted to numerical simulations of the DART impact into a rubble pile Dimorphos with different material strengths. The observations of the impact plume optical depth and the high momentum transfer from the DART impact are not consistent with impact and ejecta plume models assuming the Dimorphos cohesive strength to be as high as 5000 Pa. Models with 5 and 50 Pa Dimorphos cohesive strength provide the overall best consistency with plume opacity observations and high momentum transfer.

Type Ia supernova explosion models are inherently multidimensional

Theoretical and observational approaches to settling the important questions surrounding the progenitor systems and the explosion mechanism of normal Type Ia supernovae have thus far failed. With its unique capability to obtain continuous spectra through the near- and mid-infrared, JWST now offers completely new insights into Type Ia supernovae. In particular, observing them in the nebular phase allows us to directly see the central ejecta and thereby constrain the explosion mechanism. We aim to understand and quantify differences in the structure and composition of the central ejecta of various Type Ia supernova explosion models. We examined the currently most popular explosion scenarios using self-consistent multidimensional explosion simulations of delayed-detonation and pulsationally assisted, gravitationally confined delayed detonation Chandrasekhar-mass models and double-detonation sub-Chandrasekhar-mass and violent merger models. We find that the distribution of radioactive and stable nickel in the final ejecta, both observable in nebular spectra, are significantly different between different explosion scenarios. Therefore, comparing synthetic nebular spectra with JWST observations should allow us to distinguish between explosion models. We show that the explosion ejecta are inherently multidimensional for all models, and the Chandrasekhar-mass explosions simulated in spherical symmetry in particular lead to a fundamentally unphysical ejecta structure. Moreover, we show that radioactive and stable nickel cover a significant range of densities at a fixed velocity of the homologously expanding ejecta. Any radiation transfer postprocessing has to take these variations into account to obtain faithful synthetic observables; this will likely require multidimensional radiation transport simulations.

1991T-Like Type Ia Supernovae as an Extension of the Normal Population

Type Ia supernovae (SNe) remain poorly understood despite decades of investigation. Massive computationally intensive hydrodynamic simulations have been developed and run to model an ever-growing number of proposed progenitor channels. Further complicating the matter, a large number of subtypes of Type Ia SNe have been identified in recent decades. Due to the massive computational load required, inference of the internal structure of Type Ia SNe ejecta directly from observations using simulations has previously been computationally intractable. However, deep-learning emulators for radiation transport simulations have alleviated such barriers. We perform abundance tomography on 40 Type Ia SNe from optical spectra using the radiative transfer code TARDIS accelerated by the probabilistic DALEK deep-learning emulator. We apply a parametric model of potential outer ejecta structures to comparatively investigate abundance distributions and internal ionization fractions of intermediate-mass elements (IMEs) between normal and 1991T-like Type Ia SNe in the early phases. Our inference shows that the outer ejecta of 1991T-like Type Ia SNe is underabundant in the typical intermediate mass elements that heavily contribute to the spectral line formation seen in normal Type Ia SNe at early times. Additionally, we find that the IMEs present in 1991T-like Type Ia SNe are highly ionized compared to those in the normal Type Ia population. Finally, we conclude that the transition between normal and 1991T-like Type Ia SNe appears to be continuous observationally and that the observed differences come out of a combination of both abundance and ionization fractions in these SNe populations.

A sequence of Type Ib, IIb, II-L, and II-P supernovae from binary-star progenitors with varying initial separations

Over the last decade, evidence has accumulated that massive stars do not typically evolve in isolation but instead follow a tumultuous journey with a companion star on their way to core collapse. While Roche-lobe overflow appears instrumental for the production of a large fraction of Type Ib and Ic supernovae (SNe), variations in the initial orbital period, $P_ init $, of massive interacting binaries may also produce a wide diversity of case B, BC, or C systems, with pre-SN stars endowed from minute to massive H-rich envelopes. Focusing here on the explosion of the primary donor star, originally 12.6\ we used radiation hydrodynamics and nonlocal thermodynamic equilibrium time-dependent radiative transfer to document the gas and radiation properties of such SNe, covering Types Ib, IIb, II-L, and II-P. Variations in init $ are the root cause of the wide diversity of our SN light curves, which present single-peak, double-peak, fast-declining, or plateau-like morphologies in the $V$ band. The different ejecta structures, expansion rates, and relative abundances (e.g., H, He, and can lead to a great deal of diversity in terms of spectral line shapes (absorption versus emission strength and width) and evolution. We emphasize that Halpha is a key tracer of these modulations, and that He is an enduring optical diagnostic for the presence of He. Our grid of simulations fares well against representative Type Ib, IIb, and II-P SNe, but interaction with circumstellar material, which is ignored in this work, is likely at the origin of the tension between our Type II-L SN models and observations (e.g., of SN\,2006Y). Remaining discrepancies in the rise time to bolometric maximum of our models call for a proper account of both small-scale and large-scale structures in core-collapse SN ejecta. Discrepant Type II-P SN models, with a high plateau brightness but small spectral line widths, can be fixed by adopting more compact red-supergiant star progenitors.

X-Ray Diagnostics of Cassiopeia A’s “Green Monster”: Evidence for Dense Shocked Circumstellar Plasma

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 × 1011 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.

Xkn: a semi-analytic framework for the modelling of kilonovae

Abstract After GW170817, kilonovae have become of great interest for the astronomical, astrophysics and nuclear physics communities, due to their potential in revealing key information on the compact binary merger from which they emerge, such as the fate of the central remnant or the composition of the expelled material. Therefore, the landscape of models employed for their analysis is rapidly evolving, with multiple approaches being used for different purposes. In this paper, we present xkn, a semi-analytic framework which predicts and interprets the bolometric luminosity and the broadband light curves of such transients. xkn models the merger ejecta structure accounting for different ejecta components and non-spherical geometries. In addition to light curve models from the literature based on time scale and random-walk arguments, it implements a new model, xkn-diff, which is grounded on a solution of the radiative transfer equation for homologously expanding material. In order to characterize the variety of the ejecta conditions, it employs time and composition dependent heating rates, thermalization efficiencies and opacities. We compare xkn light curves with reference radiative transfer calculations, and we find that xkn-diff significantly improves over previous semi-analytic prescriptions. We view xkn as an ideal tool for extensive parameter estimation data analysis applications.

Modeling a Coronal Mass Ejection from an Extended Filament Channel. II. Interplanetary Propagation to 1 au

We present observations and modeling results of the propagation and impact at Earth of a high-latitude, extended filament channel eruption that commenced on 2015 July 9. The coronal mass ejection (CME) that resulted from the filament eruption was associated with a moderate disturbance at Earth. This event could be classified as a so-called “problem storm” because it lacked the usual solar signatures that are characteristic of large, energetic, Earth-directed CMEs that often result in significant geoeffective impacts. We use solar observations to constrain the initial parameters and therefore to model the propagation of the 2015 July 9 eruption from the solar corona up to Earth using 3D magnetohydrodynamic heliospheric simulations with three different configurations of the modeled CME. We find the best match between observed and modeled arrival at Earth for the simulation run that features a toroidal flux rope structure of the CME ejecta, but caution that different approaches may be more or less useful depending on the CME–observer geometry when evaluating the space weather impact of eruptions that are extreme in terms of their large size and high degree of asymmetry. We discuss our results in the context of both advancing our understanding of the physics of CME evolution and future improvements to space weather forecasting.

GRMHD Simulations of Neutron-star Mergers with Weak Interactions: r-process Nucleosynthesis and Electromagnetic Signatures of Dynamical Ejecta

Fast neutron-rich material ejected dynamically over ≲10 ms during the merger of a binary neutron star (BNS) can give rise to distinctive electromagnetic counterparts to the system’s gravitational-wave emission that serve as a “smoking gun” to distinguish between a BNS and an NS–black hole merger. We present novel ab initio modeling of the kilonova precursor and kilonova afterglow based on 3D general-relativistic magnetohydrodynamic simulations of BNS mergers with nuclear, tabulated, finite-temperature equations of state (EOSs), weak interactions, and approximate neutrino transport. We analyze dynamical mass ejection from 1.35–1.35 M ⊙ binaries, consistent with properties of the first observed BNS merger GW170817, using three nuclear EOSs that span the range of allowed compactness of 1.35 M ⊙-neutron stars. Nuclear reaction network calculations yield a robust second-to-third-peak r-process. We find few ×10−6 M ⊙ of fast (v > 0.6c) ejecta that give rise to broadband synchrotron emission on ∼years timescales, consistent with tentative evidence for excess X-ray/radio emission following GW170817. We find ≈2 × 10−5 M ⊙ of free neutrons that power a kilonova precursor on ≲ hours timescale. A boost in early UV/optical brightness by a factor of a few due to previously neglected relativistic effects, with enhancements up to ≲10 hr post-merger, is promising for future detection with UV/optical telescopes like Swift or ULTRASAT. We find that a recently predicted opacity boost due to highly ionized lanthanides at ≳70,000 K is unlikely to affect the early kilonova based on the obtained ejecta structures. Azimuthal inhomogeneities in dynamical ejecta composition for soft EOSs found here (“lanthanide/actinide pockets”) may have observable consequences for both early kilonova and late-time nebular emission.

Polytropic Behavior in the Structures of Interplanetary Coronal Mass Ejections

The polytropic process characterizes the thermodynamics of space plasma particle populations. The polytropic index, γ, is particularly important as it describes the thermodynamic behavior of the system by quantifying the changes in temperature as the system is compressed or expanded. Using Wind spacecraft plasma and magnetic field data during 1995 February–2015 December, we investigate the thermodynamic evolution in 336 interplanetary coronal mass ejection (ICME) events. For each event, we derive the index γ in the sheath and magnetic ejecta structures, along with the pre- and post-event regions. We then examine the distributions of all γ indices in these four regions and derive the entropic gradient of each, which is indicative of the ambient heating. We find that in the ICME sheath region, where wave turbulence is expected to be highest, the thermodynamics takes longest to recover into the original quasi-adiabatic process, while it recovers faster in the quieter ejecta region. This pattern creates a thermodynamic cycle, featuring a near adiabatic value γ ∼ γ a (=5/3) upstream of the ICMEs, γ a − γ ∼ 0.26 in the sheaths, γ a − γ ∼ 0.13 in the ICME ejecta, and recovers again to γ ∼ γ a after the passage of the ICME. These results expose the turbulent heating rates in the ICME plasma: the lower the polytropic index from its adiabatic value and closer to its isothermal value, the larger the entropic gradient, and thus, the rate of turbulent heating that heats the ICME plasma.

Magnetar Wind-Driven Shock Breakout Emission after Double Neutron Star Mergers: The Effect of the Anisotropy of the Merger Ejecta

A rapidly rotating and highly magnetized remnant neutron star (NS; magnetar) could survive from a merger of double NSs and drive a powerful relativistic wind. The early interaction of this wind with the previous merger ejecta can lead to shock breakout (SBO) emission mainly in ultraviolet and soft X-ray bands, which provides an observational signature for the existence of the remnant magnetar. Here, we investigate the effect of an anisotropic structure of the merger ejecta on the SBO emission. It is found that the bolometric light curve of the SBO emission can be broadened, since the SBO can occur at different times for different directions. In more detail, the profile of the SBO light curve can be highly dependent on the ejecta structure and, thus, we can in principle use the SBO light curves to probe the structure of the merger ejecta in future.

Chemical Stratification in a Long Gamma-Ray Burst Cocoon and Early-time Spectral Signatures of Supernovae Associated with Gamma-Ray Bursts

Abstract We present the results of 3D hydrodynamic simulations of a gamma-ray burst (GRB) jet emanating from a massive star with a particular focus on the formation of high-velocity quasi-spherical ejecta and the jet-induced chemical mixing. Recent early-time optical observations of supernovae associated with GRBs (e.g., GRB 171205A/SN 2017iuk) indicate a considerable amount of heavy metals in the high-velocity outer layers of the ejecta. Using our jet simulations, we show that the density and chemical structure of the outer ejecta implied by observations can be naturally reproduced by a powerful jet penetrating the progenitor star. We consider three representative jet models with a stripped massive star, a standard jet, a weak jet, and a jet choked by an extended circumstellar medium, to clarify the differences in the dynamical evolution and the chemical properties of the ejected materials. The standard jet successfully penetrates the progenitor star and creates a quasi-spherical ejecta component (cocoon). The jet-induced mixing significantly contaminates the cocoon with heavy elements that have been otherwise embedded in the inner layer of the ejecta. The weak and choked jet models fail to produce an ultrarelativistic jet but produce a quasi-spherical cocoon with different chemical properties. We discuss the impact of the different jet−star interactions on the expected early-time electromagnetic signatures of long GRBs and how to probe the jet dynamics from observations.

An analysis of the spectroscopic signatures of layering in the ejecta of Type Iax supernovae

ABSTRACT Investigations of some Type Iax supernovae (SNe Iax) have led to the suggestion that their ejecta must be layered to some degree. Such an ejecta structure has been argued as inconsistent with the well-mixed composition predicted by pure deflagrations. Based on explosion models, we create toy models in which the ejecta are artificially stratified and progressively mixed until a uniform composition is obtained. We find that models that are heavily mixed, containing burned and unburned material at all velocities, produce reasonably good agreement with SN 2012Z, for which a layered structure has been suggested. We also discuss how existing ejecta compositions determined for SNe Iax do not necessarily contradict pure deflagration models and may be consistent with a steeper density profile. We investigate previous claims that differences in line profile shapes may be due to strong blending, by presenting a series of models with different plasma states. These models indicate that blending could indeed explain differences in the observed profiles. Alternatively, stratification could also explain such differences; however, all of our models indicate that this does not necessarily require stratification in abundance. Sufficient stratification in the ionization state can be achieved even for a well-mixed model. Based on our analysis, we demonstrate that there is insufficient evidence to suggest the ejecta of SNe Iax must be layered and therefore argue the pure deflagration scenario is not ruled out, even for the brightest SNe Iax. Our analysis does not indicate the ejecta cannot be layered to some degree, but observations within days of explosion are necessary to determine the extent to which the outer ejecta could be layered.

ALMA CO Observations of the Mixed-morphology Supernova Remnant W49B: Efficient Production of Recombining Plasma and Hadronic Gamma Rays via Shock–Cloud Interactions

Abstract We carried out new CO(J = 2–1) observations toward the mixed-morphology supernova remnant (SNR) W49B with the Atacama Large Millimeter/submillimeter Array. We found that CO clouds at ∼10 km s−1 show a good spatial correspondence to the synchrotron radio continuum as well as to an X-ray deformed shell. The bulk mass of molecular clouds accounts for the western part of the shell, not the eastern shell, where near-infrared H2 emission is detected. The molecular clouds at ∼10 km s−1 show higher kinetic temperatures of ∼20–60 K, suggesting that modest shock heating occurred. The expanding motion of the clouds with ΔV ∼ 6 km s−1 was formed by strong winds from the progenitor system. We argue that the barrel-like structure of Fe-rich ejecta was possibly formed not only by an asymmetric explosion, but also by interactions with dense molecular clouds. We also found a negative correlation between the CO intensity and the electron temperature of recombining plasma, implying that the origin of the high-temperature recombining plasma in W49B can be understood to be the thermal conduction model. The total energy of accelerated cosmic-ray protons W p is estimated to be ∼2 × 1049 erg by adopting an averaged gas density of ∼650 ± 200 cm−3. The SNR age–W p diagram indicates that W49B shows one of the highest in situ values of W p among gamma-ray-bright SNRs.

Tight multimessenger constraints on the neutron star equation of state from GW170817 and a forward model for kilonova light-curve synthesis

ABSTRACT We present a rapid analytic framework for predicting kilonova light curves following neutron star (NS) mergers, where the main input parameters are binary-based properties measurable by gravitational wave detectors (chirp mass and mass ratio, orbital inclination) and properties dependent on the nuclear equation of state (tidal deformability, maximum NS mass). This enables synthesis of a kilonova sample for any NS source population, or determination of the observing depth needed to detect a live kilonova given gravitational wave source parameters in low latency. We validate this code, implemented in the public mosfit package, by fitting it to GW170817. A Bayes factor analysis overwhelmingly (B > 1010) favours the inclusion of an additional luminosity source in addition to lanthanide-poor dynamical ejecta during the first day. This is well fit by a shock-heated cocoon model, though differences in the ejecta structure, opacity or nuclear heating rate cannot be ruled out as alternatives. The emission thereafter is dominated by a lanthanide-rich viscous wind. We find the mass ratio of the binary is q = 0.92 ± 0.07 (90 per cent credible interval). We place tight constraints on the maximum stable NS mass, MTOV $=2.17^{+0.08}_{-0.11}$ M⊙. For a uniform prior in tidal deformability, the radius of a 1.4-M⊙ NS is R1.4 ∼ 10.7 km. Re-weighting with a prior based on equations of state that support our credible range in MTOV, we derive a final measurement R1.4 $=11.06^{+1.01}_{-0.98}$ km. Applying our code to the second gravitationally detected NS merger, GW190425, we estimate that an associated kilonova would have been fainter (by ∼0.7 mag at 1 d post-merger) and declined faster than GW170817, underlining the importance of tuning follow-up strategies individually for each GW-detected NS merger.

Double detonations of sub-MCh CO white dwarfs: variations in Type Ia supernovae due to different core and He shell masses

Sub-Chandrasekhar mass carbon-oxygen white dwarfs with a surface helium shell have been proposed as progenitors of Type Ia supernovae (SNe Ia). If true, the resulting thermonuclear explosions should be able to account for at least some of the range of SNe Ia observables. To study this, we conducted a parameter study based on three-dimensional simulations of double detonations in carbon-oxygen white dwarfs with a helium shell, assuming different core and shell masses. An admixture of carbon to the shell and solar metallicity are included in the models. The hydrodynamic simulations were carried out using the AREPO code. This allowed us to follow the helium shell detonation with high numerical resolution, and this improves the reliability of predicted nucleosynthetic shell detonation yields. The addition of carbon to the shell leads to a lower production of 56Ni, while including solar metallicity increases the production of intermediate mass elements. The production of higher mass elements is further shifted to stable isotopes at solar metallicity. Moreover, we find different core detonation ignition mechanisms depending on the core and shell mass configuration. This has an influence on the ejecta structure. We present the bolometric light curves predicted from our explosion simulations using the Monte Carlo radiative transfer code ARTIS and make comparisons with bolometric SNe Ia data. The bolometric light curves of our models show a range of brightnesses, which is able to account for subluminous to normal brightness SNe Ia. We show the model bolometric width-luminosity relation compared to data for a range of model viewing angles. We find that, on average, our brighter models lie within the observed data. The ejecta asymmetries produce a wide distribution of observables, which might account for outliers in the data. However, the models overestimate the extent of this compared to data. We also find that the bolometric decline rate over 40 days, Δm40(bol), appears systematically faster than data.

On the Origin of the Asymmetry of the Ejecta Structure and Explosion of G350.1–0.3

Abstract We present X-ray analysis of the ejecta of supernova remnant (SNR) G350.1–0.3 observed with Chandra and Suzaku, clarify the ejecta’s kinematics over a decade, and obtain a new observational clue to understanding the origin of the asymmetric explosion. Two images from Chandra X-ray Observatory taken in 2009 and 2018 are analyzed with several methods and enable us to measure the velocities in the plane of the sky. A maximum velocity is 4640 ± 290 km s−1 (0.218 ± 0.014 arcsec yr−1) in the eastern region in the remnant. These findings trigger us to scrutinize the Doppler effects in the spectra of the thermal emission, and the velocities in the line-of-sight direction are estimated to be 1000 km s−1. The results are confirmed by analyzing the spectra of Suzaku. Combining the proper motions and line-of-sight velocities, the ejecta’s 3D velocities are ∼3000–5000 km s−1. The center of the explosion is more stringently constrained by finding the optimal time to reproduce the observed spatial expansion. Our findings that the age of the SNR is estimated at most to be 655 yr and the CCO is observed as a point source object against the SNR strengthen the “hydrodynamical kick” hypothesis on the origin of the remnant.

The fully developed remnant of a neutrino-driven supernova

Context. The remnants of core-collapse supernovae (SNe) are probes of the physical processes associated with their parent SNe. Aims. Here we aim to explore to which extent the remnant keeps memory of the asymmetries that develop stochastically in the neutrino-heating layer due to hydrodynamic instabilities (e.g., convective overturn and the standing accretion shock instability; SASI) during the first second after core bounce. Methods. We coupled a three-dimensional (3D) hydrodynamic model of a neutrino-driven SN explosion, which has the potential to reproduce the observed morphology of the Cassiopeia A (Cas A) remnant, with 3D (magneto)-hydrodynamic simulations of the remnant formation. The simulations cover ≈2000 yr of expansion and include all physical processes relevant to describe the complexities in the SN evolution and the subsequent interaction of the stellar debris with the wind of the progenitor star. Results. The interaction of large-scale asymmetries left from the earliest phases of the explosion with the reverse shock produces, at the age of ≈350 yr, an ejecta structure and a remnant morphology which are remarkably similar to those observed in Cas A. Small-scale structures in the large-scale Fe-rich plumes that were created during the initial stages of the SN, combined with hydrodynamic instabilities that develop after the passage of the reverse shock, naturally produce a pattern of ring- and crown-like structures of shocked ejecta. The consequence is a spatial inversion of the ejecta layers with Si-rich ejecta being physically interior to Fe-rich ejecta. The full-fledged remnant shows voids and cavities in the innermost unshocked ejecta, which are physically connected with ring-like features of shocked ejecta in the main shell in most cases, resulting from the expansion of Fe-rich plumes and their inflation due to the decay of radioactive species. The asymmetric distributions of 44Ti and 56Fe, which are mostly concentrated in the northern hemisphere, and pointing opposite to the kick velocity of the neutron star, as well as their abundance ratio are both compatible with those inferred from high-energy observations of Chandra and NuSTAR. Finally, the simulations show that the fingerprints of the SN can still be visible ≈2000 yr after the explosion. Conclusions. The main asymmetries and features observed in the ejecta distribution of Cas A can be explained by the interaction of the reverse shock with the initial large-scale asymmetries that developed from stochastic processes (e.g., convective overturn and SASI activity) that originate during the first seconds of the SN blast.

Deriving CME Density From Remote Sensing Data and Comparison to In‐Situ Measurements

AbstractWe determine the three‐dimensional geometry and deprojected mass of 29 well‐observed coronal mass ejections (CMEs) and their interplanetary counterparts (ICMEs) using combined Solar Terrestrial Relations Observatory ‐ Solar and Heliospheric Observatory white‐light data. From the geometry parameters, we calculate the volume of the CME for the magnetic ejecta (flux‐rope type geometry) and sheath structure (shell‐like geometry resembling the (I)CME frontal rim). Working under the assumption that the CME mass is roughly equally distributed within a specific volume, we expand the CME self‐similarly and calculate the CME density for distances close to the Sun (15–30 Rs) and at 1 AU. Specific trends are derived comparing calculated and in‐situ measured proton densities at 1 AU, though large uncertainties are revealed due to the unknown mass and geometry evolution: (1) a moderate correlation for the magnetic structure having a mass that stays rather constant (cc ≈ 0.56 − 0.59), and (2) a weak correlation for the sheath density (cc ≈ 0.26) by assuming the sheath region is an extra mass—as expected for a mass pile‐up process—that is in its amount comparable to the initial CME deprojected mass. High correlations are derived between in‐situ measured sheath density and the solar wind density (cc ≈ −0.73) and solar wind speed (cc ≈ 0.56) as measured 24 h ahead of the arrival of the disturbance. This gives additional confirmation that the sheath‐plasma indeed stems from piled‐up solar wind material. While the CME interplanetary propagation speed is not related to the sheath density, the size of the CME may play some role in how much material could be piled up.

Radiative-transfer modeling of supernovae in the nebular-phase

Supernova (SN) explosions play a pivotal role in the chemical evolution of the Universe and the origin of life through the metals they release. Nebular phase spectroscopy constrains such metal yields, for example through forbidden line emission associated with O I, Ca II, Fe II, or Fe III. Fluid instabilities during the explosion produce a complex 3D ejecta structure, with considerable macroscopic, but no microscopic, mixing of elements. This structure sets a formidable challenge for detailed nonlocal thermodynamic equilibrium radiative transfer modeling, which is generally limited to 1D in grid-based codes. Here, we present a novel and simple method that allows for macroscopic mixing without any microscopic mixing, thereby capturing the essence of mixing in SN explosions. With this new technique, the macroscopically mixed ejecta are built by shuffling the shells from the unmixed coasting ejecta in mass space, or equivalently in velocity space. The method requires no change to the radiative transfer, but it necessitates high spatial resolution to resolve the rapid variation in composition with depth inherent to this shuffled-shell structure. We show the results for a few radiative-transfer simulations for a Type II SN explosion from a 15 M⊙ progenitor star. Our simulations capture the strong variations in temperature or ionization between the various shells that are rich in H, He, O, or Si. Because of nonlocal energy deposition, γ rays permeate through an extended region of the ejecta, making the details of the shell arrangement unimportant. The greater physical consistency of the method delivers spectral properties at nebular times that are more reliable, in particular in terms of individual emission line strengths, which may serve to constrain the SN yields as well as the progenitor mass for core collapse SNe. The method works for all SN types.