Unshocked Regions Research Articles

A JWST Survey of the Supernova Remnant Cassiopeia A

We present initial results from a James Webb Space Telescope (JWST) survey of the youngest Galactic core-collapse supernova remnant, Cassiopeia A (Cas A), made up of NIRCam and MIRI imaging mosaics that map emission from the main shell, interior, and surrounding circumstellar/interstellar material (CSM/ISM). We also present four exploratory positions of MIRI Medium Resolution Spectrograph integral field unit spectroscopy that sample ejecta, CSM, and associated dust from representative shocked and unshocked regions. Surprising discoveries include (1) a weblike network of unshocked ejecta filaments resolved to ∼0.01 pc scales exhibiting an overall morphology consistent with turbulent mixing of cool, low-entropy matter from the progenitor’s oxygen layer with hot, high-entropy matter heated by neutrino interactions and radioactivity; (2) a thick sheet of dust-dominated emission from shocked CSM seen in projection toward the remnant’s interior pockmarked with small (∼1″) round holes formed by ≲0.″1 knots of high-velocity ejecta that have pierced through the CSM and driven expanding tangential shocks; and (3) dozens of light echoes with angular sizes between ∼0.″1 and 1′ reflecting previously unseen fine-scale structure in the ISM. NIRCam observations place new upper limits on infrared emission (≲20 nJy at 3 μm) from the neutron star in Cas A’s center and tightly constrain scenarios involving a possible fallback disk. These JWST survey data and initial findings help address unresolved questions about massive star explosions that have broad implications for the formation and evolution of stellar populations, the metal and dust enrichment of galaxies, and the origin of compact remnant objects.

XMM-Newton view of the shock heating in an early merging cluster, CIZA J1358.9−4750

AbstractCIZA J1358.9−4750 is a nearby galaxy cluster in the early phase of a major merger. The two-dimensional temperature map using XMM-Newton EPIC-PN observation confirms the existence of a high-temperature region, which we call the “hot region,” in the “bridge region” connecting the two clusters. The ∼500 kpc wide region between the south-east and north-west boundaries also has higher pseudo-pressure compared to the unshocked regions, suggesting the existence of two shocks. The southern shock front is clearly visible in the X-ray surface brightness image and has already been reported by Kato et al. (2015, PASJ, 67, 71). The northern one, on the other hand, is newly discovered. To evaluate their Mach number, we constructed a three-dimensional toy merger model with overlapping shocked and unshocked components in the line of sight. The unshocked and pre-shock intracluster medium (ICM) conditions are estimated based on those outside the interacting bridge region, assuming point symmetry. The hot-region spectra are modeled with two-temperature thermal components, assuming that the shocked condition follows the Rankin–Hugoniot relation with the pre-shock condition. As a result, the shocked region is estimated to have a line-of-sight depth of ∼1 Mpc with a Mach number of ∼1.3 in the south-east shock and ∼1.7 in the north-west shock. The age of the shock waves is estimated to be ∼260 Myr. This three-dimensional merger model is consistent with the Sunyaev–Zel’dovich signal obtained using the Planck observation within the cosmic microwave background fluctuations. The total flow of the kinetic energy of the ICM through the south-east shock was estimated to be ∼2.2 × 1042 erg s−1. Assuming that $10\%$ of this energy is converted into ICM turbulence, the line–of–sight velocity dispersion is calculated to be ∼200 km s−1, which is basically resolvable via upcoming high spectral resolution observations.

Dust grains from the heart of supernovae

Dust grains are classically thought to form in the winds of asymptotic giant branch (AGB) stars. However, there is increasing evidence today for dust formation in supernovae (SNe). To establish the relative importance of these two classes of stellar sources of dust, it is important to know the fraction of freshly formed dust in SN ejecta that is able to survive the passage of the reverse shock and be injected in the interstellar medium. With this aim, we have developed a new code, GRASH_Rev, that allows following the dynamics of dust grains in the shocked SN ejecta and computing the time evolution of the mass, composition, and size distribution of the grains. We considered four well-studied SNe in the Milky Way and Large Magellanic Cloud: SN 1987A, CasA, the Crab nebula, and N49. These sources have been observed with both Spitzer and Herschel, and the multiwavelength data allow a better assessment the mass of warm and cold dust associated with the ejecta. For each SN, we first identified the best explosion model, using the mass and metallicity of the progenitor star, the mass of 56Ni, the explosion energy, and the circumstellar medium density inferred from the data. We then ran a recently developed dust formation model to compute the properties of freshly formed dust. Starting from these input models, GRASH_Rev self-consistently follows the dynamics of the grains, considering the effects of the forward and reverse shock, and allows predicting the time evolution of the dust mass, composition, and size distribution in the shocked and unshocked regions of the ejecta. All the simulated models aagree well with observations. Our study suggests that SN 1987A is too young for the reverse shock to have affected the dust mass. Hence the observed dust mass of 0.7−0.9 M⊙ in this source can be safely considered as indicative of the mass of freshly formed dust in SN ejecta. Conversely, in the other three SNe, the reverse shock has already destroyed between 10−40% of the initial dust mass. However, the largest dust mass destruction is predicted to occur between 103 and 105 yr after the explosions. Since the oldest SN in the sample has an estimated age of 4800 yr, current observations can only provide an upper limit to the mass of SN dust that will enrich the interstellar medium, the so-called effective dust yields. We find that only between 1−8% of the currently observed mass will survive, resulting in an average SN effective dust yield of (1.55 ± 1.48) × 10-2M⊙. This agrees well with the values adopted in chemical evolution models that consider the effect of the SN reverse shock. We discuss the astrophysical implications of our results for dust enrichment in local galaxies and at high redshift.

Low-mass planets in nearly inviscid disks: numerical treatment

Embedded planets disturb the density structure of the ambient disk and gravitational back-reaction will induce possibly a change in the planet's orbital elements. The accurate determination of the forces acting on the planet requires careful numerical analysis. Recently, the validity of the often used fast orbital advection algorithm (FARGO) has been put into question, and special numerical resolution and stability requirements have been suggested. In this paper we study the process of planet-disk interaction for small mass planets of a few Earth masses, and reanalyze the numerical requirements to obtain converged and stable results. One focus lies on the applicability of the FARGO-algorithm. Additionally, we study the difference of two and three-dimensional simulations, compare global with local setups, as well as isothermal and adiabatic conditions. We study the influence of the planet on the disk through two- and three-dimensional hydrodynamical simulations. To strengthen our conclusions we perform a detailed numerical comparison where several upwind and Riemann-solver based codes are used with and without the FARGO-algorithm. With respect to the wake structure and the torque density acting on the planet we demonstrate that the FARGO-algorithm yields correct results, and that at a fraction of the regular cpu-time. We find that the resolution requirements for achieving convergent results in unshocked regions are rather modest and depend on the pressure scale height of the disk. By comparing the torque densities of 2D and 3D simulations we show that a suitable vertical averaging procedure for the force gives an excellent agreement between the two. We show that isothermal and adiabatic runs can differ considerably, even for adiabatic indices very close to unity.

Shock Geometry and Inverse Compton Emission from the Wind of a Binary Pulsar

AbstractPSR B1259 – 63 is a 47 ms radio pulsar with a high spin-down luminosity which is in a close, highly eccentric 3·5 yr orbit about a bright stellar companion. The binary system may be a detectable source of hard ã γ-rays produced by inverse Compton scattering of photons from the B2e star SS2883 by electrons and positrons in the pulsar wind. The star provides an enormous density of optical photons in the vicinity of the pulsar, particularly at epochs near periastron. We calculate the emission from the unshocked region of the pulsar wind, assuming that it terminates at a shock where it attains pressure balance with the companion’s wind. The spectra and light curves for the inverse Compton emission from the shock-terminated wind are compared with those for an unterminated wind. If the pulsar’s wind is weaker than that from the companion star, the termination of the wind decreases the inverse Compton flux, particularly near periastron. The termination shock geometry has the effect of decreasing the asymmetry of the γ-ray light curve around periastron, which arises because of the asymmetrical variation of the scattering angle.

Nonequilibrium Thermodynamics and Cosmological Pancake Formation

We investigate the influence of non equilibrium thermodynamics on cosmological structure formation. In this paper, we consider the collapse of planar perturbations usually called Zel'dovich pancakes. We have developed for that purpose a new two fluids (gas and dark matter) hydrodynamical code, with three different thermodynamical species: electrons, ions and neutral particles (T_e\ne T_i \ne T_n). We describe in details the complex structure of accretion shock waves. We include several relevant processes for a low density, high temperature, collisional plasma such as non-equilibrium chemical reactions, cooling, shock heating, thermal energy equipartition between electrons, ions and neutral particles and electronic conduction. We find two different regions in the pancake structure: a thermal precursor ahead of the compression front and an equipartition wave after the compression front where electrons and ions temperatures differ significantly. This complex structure may have two interesting consequences: pre-heating of unshocked regions in the vicinity of massive X-ray clusters and ions and electrons temperatures differences in the outer regions of X-rays clusters.

A self‐consistent model for the particles and fields upstream of an outgassing comet: 2. A time‐dependent description

In a previous paper [Flammer et al., 1991], we determined the global variation of the magnetic field and solar wind flow parameters in the unshocked region upstream of an outgassing comet using a kinetic treatment for the cometary ions. Two different assumptions were made concerning the cometary ion distribution function: the pickup cometary ions formed either a velocity space gyrotropic ring distribution or a velocity space isotropic shell distribution in the solar wind frame of reference. In the present paper we consider the general case wherein the newly picked up ions are elastically pitch angle scattered from the initial ring distribution to a shell with some characteristic time scale τ. Using theoretically determined parameters which reflect the conditions expected at comet Kopff (a typical short‐period comet and a possible target for the CRAF/Cassini mission) for various heliocentric distances, we determine how the pitch angle scattering rate affects the global morphology.