Spectral analysis of X-ray emission from ejecta in supernova remnants (SNRs) is hampered by the low spectral resolution of CCD cameras, which creates a degeneracy between the best-fit values of abundances and emission measure. The combined contribution of shocked ambient medium and ejecta to the X-ray emission complicates the determination of the ejecta mass and chemical composition, leading to big uncertainties in mass estimates and it can introduce a bias in the comparison between the observed ejecta composition and the yields predicted by explosive nucleosynthesis. We explore the capabilities of present and future spectral instruments with the aim of identifying a spectral feature which may allow us to discriminate between metal-rich and pure-metal plasmas in X-ray spectra of SNRs. We studied the behavior of the most common X-ray emission processes of an optically thin plasma in the high-abundance regime. We investigated spectral features of bremsstrahlung, radiative recombination continua (RRC) and line emission, by exploring a wide range of chemical abundances, temperatures and ionization parameters. We synthesized X-ray spectra from a 3D hydrodynamic (HD) simulation of Cas A, by using the response matrix from the Chandra/ACIS-S CCD detector and that of the XRISM/Resolve X-ray calorimeter. We found that a bright RRC shows up when the plasma is made of pure-metal ejecta, and a high spectral resolution is needed to identify this ejecta signature. We verified the applicability of our novel diagnostic tool and we propose a promising target for the future detection of such spectral feature: the southeastern Fe-rich clump of Cas A. While there is no way to unambiguously reveal pure-metal ejecta emission with CCD detectors, X-ray calorimeters will be able to pinpoint the presence of pure-metal RRC and to recover correctly absolute mass and the chemical composition of the ejecta.

Recent X-ray studies revealed over-ionized recombining plasmas (RPs) in a dozen mixed-morphology (MM) supernova remnants (SNRs). However, the physical process of the over-ionization has not been fully understood yet. Here we report on spatially resolved spectroscopy of X-ray emission from W44, one of the over-ionized MM-SNRs, using XMM-Newton data from deep observations, aiming to clarify the physical origin of the over-ionization. We find that combination of low electron temperature and low recombination timescale is achieved in the region interacting with dense molecular clouds. Moreover, a clear anti-correlation between the electron temperature and the recombining timescale is obtained from each of the regions with and without the molecular clouds. The results are well explained if the plasma was over-ionized by rapid cooling through thermal conduction with the dense clouds hit by the blast wave of W44. Given that a few other over-ionized SNRs show evidence for adiabatic expansion as the major driver of the rapid cooling, our new result indicates that both processes can contribute to over-ionization in SNRs, with the dominant channel depending on the evolutionary stage.

(Abridged) We aim at linking the dynamical and radiative properties of the remnant of SN 1987A to the geometrical and physical characteristics of the parent aspherical SN explosion and to the internal structure of its progenitor star. We performed 3D hydrodynamic simulations which describe the long-term evolution of SN 1987A from the onset of the SN to the full-fledged remnant at the age of 50 years, accounting for the pre-SN structure of the progenitor star. The simulations include all physical processes relevant for the complex phases of SN evolution and for the interaction of the SNR with the highly inhomogeneous ambient environment around SN 1987A. From the simulations, we synthesize observables to be compared with observations. By comparing the model results with observations, we constrained the initial SN anisotropy causing Doppler shifts observed in emission lines of heavy elements from ejecta, and leading to the remnant evolution observed in the X-ray band in the last 30 years. In particular, we found that the high mixing of ejecta unveiled by high redshifts and broadenings of [FeII] and $^{44}$Ti lines require a highly asymmetric SN explosion channeling a significant fraction of energy along an axis almost lying in the plane of the central equatorial ring around SN 1987A, roughly along the line-of-sight but with an offset of 40 deg, with the lobe propagating away from the observer slightly more energetic than the other. We found unambiguously that the observed distribution of ejecta and the dynamical and radiative properties of the SNR can be best reproduced if the structure of the progenitor star was that of a blue supergiant resulted from the merging of two massive stars.

We perform three-dimensional hydrodynamic simulations of aspherical core-collapse supernovae focusing on the matter mixing in SN 1987A. The impacts of four progenitor (pre-supernova) models and parameterized aspherical explosions are investigated. The four pre-supernova models include a blue supergiant (BSG) model based on a slow merger scenario developed recently for the progenitor of SN 1987A (Urushibata et al. 2018). The others are a BSG model based on a single star evolution and two red supergiant (RSG) models. Among the investigated explosion (simulation) models, a model with the binary merger progenitor model and with an asymmetric bipolar-like explosion, which invokes a jetlike explosion, best reproduces constraints on the mass of high velocity $^{56}$Ni, as inferred from the observed [Fe II] line profiles. The advantage of the binary merger progenitor model for the matter mixing is the flat and less extended $\rho \,r^3$ profile of the C+O core and the helium layer, which may be characterized by the small helium core mass. From the best explosion model, the direction of the bipolar explosion axis (the strongest explosion direction), the neutron star (NS) kick velocity, and its direction are predicted. Other related implications and future prospects are also given.

Magnetars are regarded as the most magnetized neutron stars in the Universe. Aiming to unveil what kinds of stars and supernovae can create magnetars, we have performed a state-of-the-art spatially resolved spectroscopic X-ray study of the supernova remnants (SNRs) Kes 73, RCW 103, and N49, which host magnetars 1E 1841-045, 1E 161348-5055, and SGR 0526-66, respectively. The three SNRs are O- and Ne-enhanced and are evolving in the interstellar medium with densities of >1--2 cm$^{-3}$. The metal composition and dense environment indicate that the progenitor stars are not very massive. The progenitor masses of the three magnetars are constrained to be < 20 Msun (11--15 Msun for Kes 73, < 13 Msun for RCW 103, and ~13 --17 Msun for N49). Our study suggests that magnetars are not necessarily made from very massive stars, but originate from stars that span a large mass range. The explosion energies of the three SNRs range from $10^{50}$ erg to ~2$\times 10^{51}$ erg, further refuting that the SNRs are energized by rapidly rotating (millisecond) pulsars. We report that RCW 103 is produced by a weak supernova explosion with significant fallback, as such an explosion explains the low explosion energy (~$10^{50}$ erg), small observed metal masses ($M_{\rm O}\sim 4\times 10^{-2}$ Msun and $M_{\rm Ne}\sim 6\times 10^{-3}$ Msun), and sub-solar abundances of heavier elements such as Si and S. Our study supports the fossil field origin as an important channel to produce magnetars, given the normal mass range ($M_{\rm ZAMS} < 20$ Msun) of the progenitor stars, the low-to-normal explosion energy of the SNRs, and the fact that the fraction of SNRs hosting magnetars is consistent with the magnetic OB stars with high fields.

Coronal mass ejections (CMEs), often associated with flares, are the most powerful magnetic phenomena occurring on the Sun. Stars show magnetic activity levels up to 10^4 times higher, and CME effects on stellar physics and circumstellar environments are predicted to be significant. However, stellar CMEs remain observationally unexplored. Using time-resolved high-resolution X-ray spectroscopy of a stellar flare on the active star HR 9024 observed with Chandra/HETGS, we distinctly detected Doppler shifts in S XVI, Si XIV, and Mg XII lines that indicate upward and downward motions of hot plasmas (~10-25 MK) within the flaring loop, with velocity v~100-400 km/s, in agreement with a model of flaring magnetic tube. Most notably, we also detected a later blueshift in the O VIII line which reveals an upward motion, with v=90+/-30 km/s, of cool plasma (~4 MK), that we ascribe to a CME coupled to the flare. From this evidence we were able to derive a CME mass of 1x10^21 g and a CME kinetic energy of 5x10^34 erg. These values provide clues in the extrapolation of the solar case to higher activity levels, suggesting that CMEs could indeed be a major cause of mass and angular momentum loss.

Astrophysical shocks at all scales, from those in the heliosphere up to the cosmological shock waves, are typically "collisionless", because the thickness of their jump region is much shorter than the collisional mean free path. Across these jumps, electrons, protons, and ions are expected to be heated at different temperatures. Supernova remnants (SNRs) are ideal targets to study collisionless processes because of their bright post-shock emission and fast shocks. Although optical observations of Balmer-dominated shocks in young SNRs showed that the post-shock proton temperature is higher than the electron temperature, the actual dependence of the post-shock temperature on the particle mass is still widely debated. We tackle this longstanding issue through the analysis of deep multi-epoch and high-resolution observations of the youngest nearby supernova remnant, SN 1987A, made with the Chandra X-ray telescope. We introduce a novel data analysis method by studying the observed spectra in close comparison with a dedicated full 3-D hydrodynamic simulation. The simulation is able to reproduce self-consistently the whole broadening of the spectral lines of many ions altogether. We can therefore measure the post shock temperature of protons and selected ions through comparison of the model with observations. We have obtained information about the heating processes in collisional shocks by finding that the ion to proton temperature ratio is always significantly higher than one and increases linearly with the ion mass for a wide range of masses and shock parameters.

Young stars exhibit short-term photometric variability caused by mass accretion events from circumstellar disks, the presence of dusty warps within the inner disks, starspots that rotate across the stellar surfaces, and flares. Long-term variability also occurs owing to starspot longevity and cycles, and from changes in stellar angular momenta and activity as the stars age. We propose to observe the Carina star-forming region in different bands with a cadence of 30 minutes every night for one week per year to clarify the nature of both the short-term and long-term variability of the thousands of young stars in this region. By obtaining well-sampled multicolor lightcurves of this dense young cluster, LSST would acquire the first statistically significant data on how these objects vary on both short and long timescales. This information will allow us to relate the observed variability to stellar properties such as mass, age, binarity, and to environmental properties such as location within or exterior to the H II region, and to the presence or absence of a circumstellar disk.

Aims. We investigate the role played by a pre-supernova (SN) ambient magnetic field on the dynamics of the expanding remnant of SN 1987A and the origin and evolution of the radio emission from the remnant, in particular, during the interaction of the blast wave with the nebula surrounding the SN. Methods. We model the evolution of SN 1987A from the breakout of the shock wave at the stellar surface to the expansion of its remnant through the surrounding nebula by 3D MHD simulations. The model considers the radiative cooling, the deviations from equilibrium of ionization, the deviation from temperature-equilibration between electrons and ions, and a plausible configuration of the pre-SN ambient magnetic field. From the simulations, we synthesize both thermal X-ray and non-thermal radio emission and compare the model results with observations. Results. The presence of an ambient magnetic field reduces the erosion and fragmentation of the dense equatorial ring after the impact of the SN blast wave. As a result, the ring survives the passage of the blast, at least, during the time covered by the simulations (40 years). Our model is able to reproduce the morphology and lightcurves of SN 1987A in both X-ray and radio bands. The model reproduces the observed radio emission if the flux originating from the reverse shock is heavily suppressed. In this case, the radio emission originates mostly from the forward shock traveling through the H II region. Possible mechanisms for the suppression of emission from the reverse shock are investigated. We find that synchrotron self-absorption and free-free absorption have negligible effects on the emission during the interaction with the nebula. We suggest that the emission from the reverse shock at radio frequencies might be limited by highly magnetized ejecta.

Vela X is the prototypical example of a pulsar wind nebula whose morphology and detailed structure have been affected by the interaction with the reverse shock of its host supernova remnant. The resulting complex of filamentary structure and mixed-in ejecta embedded in a nebula that is offset from the pulsar provides the best example we have of this middle-age state that characterizes a significant fraction of composite SNRs, and perhaps all of the large-diameter PWNe seen as TeV sources. Here we report on an XMM-Newton Large Project study of Vela X, supplemented by additional Chandra observations. Through broad spectral modeling as well as detailed spectral investigations of discrete emission regions, we confirm previous studies that report evidence for ejecta material within Vela X, and show that equivalent width variations of O VII and O VIII are consistent with temperature maps within the PWN that show low-temperature regions where the projected SNR emission appears to dominate emission from the ejecta. We identify spectral variations in the nonthermal emission, with hard emission being concentrated near the pulsar. We carry out investigations of the Vela X "cocoon" structure and, with hydrodynamical simulations, show that its overall properties are consistent with structures formed in the late-phase evolution of a composite SNR expanding into a surrounding medium with a density gradient, with ejecta material being swept beyond the pulsar and compressed into an elongated structure in the direction opposite the high external density.

We present ~400ks NuSTAR observations of the northeast (NE) and southwest (SW) non-thermal limbs of the Galactic SNR SN1006. We discovered three sources with X-ray emission detected at >50keV. Two of them are identified as background AGN. We extract the NuSTAR spectra from a few regions along the non-thermal limbs and jointly analyze them with the XMM-Newton spectra and the radio data. The broad-band radio/X-ray spectra can be well described with a synchrotron emission model from a single population of CR electrons with a power law energy distribution and an exponential cutoff. The power law index of the electron particle distribution function (PDF) is ~1.88-1.95 for both the NE and SW limbs, and we do not find significant evidence for a variation of this index at different energy (curvature). There are significant spatial variations of the synchrotron emission parameters. The highest energy electrons are accelerated in regions with the lowest expansion velocity, which is opposite to what has been found in the Tycho's SNR. In addition to a gradual steepening of synchrotron emission from the center of the non-thermal limbs to larger azimuthal angles, we also find that both the emission spectrum and the PDF are significantly flatter in three regions in the SW limb where the shock encounters higher density ambient medium. The NE limb also shows significantly higher cutoff energy in the PDF than the SW limb. By comparing with the roughly symmetric TeV emission and largely asymmetric GeV emission from the two non-thermal limbs, we conclude that the asymmetry in the ambient medium and magnetic fields may have largely modified the acceleration and emission of CR leptons.

We present very-high-resolution 1D MHD simulations of the late-stage supernova remnants (SNR). In the post-adiabatic stage the magnetic field has an important and significant dynamical effect on the shock dynamics, the flow structure, and hence the acceleration and emission of cosmic rays. We find that the tangential component of the magnetic field provides pressure support that to a fair degree prevents the collapse of the radiative shell and thus limits the total compression ratio of the partially or fully radiative forward shock. A consequence is that the spectra of cosmic rays would not be as hard as in hydrodynamic simulations. We also investigated the effect on the flow profiles of the magnetic-field inclination and a large-scale gradient in the gas density and/or the magnetic field. A positive density gradient shortens the evolutionary stages whereas a shock obliquity lowers the shock compression. The compression of the tangential component of the magnetic field leads to its dominance in the downstream region of post-adiabatic shocks for a wide range of orientation of the upstream field, which may explain why one preferentially observes tangential radio polarization in old SNRs. As most cosmic rays are produced at late stages of SNR evolution, the post-adiabatic phase and the influence of the magnetic field during it are most important for modeling the cosmic-ray acceleration at old SNRs and the gamma-ray emission from late-stage SNRs interacting with clouds.

IC443 is a supernova remnant located in a quite complex environment since it interacts with nearby clouds. Indications for the presence of overionized plasma have been found though the possible physical causes of overionization are still debated. Moreover, because of its peculiar position and proper motion, it is not clear if the pulsar wind nebula (PWN) within the remnant is the relic of the IC443 progenitor star or just a rambling one seen in projection on the remnant. Here we address the study of IC443 plasma in order to clarify the relationship PWN-remnant, the presence of overionization and the origin of the latter. We analyzed two \emph{XMM-Newton} observations producing background-subtracted, vignetting-corrected and mosaicked images in two different energy bands and we performed a spatially resolved spectral analysis of the X-ray emission. We identified an elongated (jet-like) structure with Mg-rich plasma in overionization. The head of the jet is interacting with a molecular cloud and the jet is aligned with the position of the PWN at the instant of the supernova explosion. Interestingly, the direction of the jet of ejecta is somehow consistent with the direction of the PWN jet. Our discovery of a jet of ejecta in IC443 enlarge the sample of core-collapse SNRs with collimated ejecta structures. IC443's jet is the first one which shows overionized plasma, possibly associated with the adiabatic expansion of ejecta. The match between the jet's direction and the original position of the PWN strongly supports the association between the neutron star and IC443.

Several observations of stellar jets show evidence of X-ray emitting shocks close to the launching site. In some cases, the shocked features appear to be stationary, also for YSOs at different stages of evolution. We study the case of HH 154, the jet originating from the embedded binary Class 0/I protostar IRS 5, and the case of the jet associated to DG Tau, a more evolved Class II disk-bearing source or Classical T Tauri star (CTTS), both located in the Taurus star-forming region. We aim at investigating the effect of perturbations in X-ray emitting stationary shocks in stellar jets; the stability and detectability in X-rays of these shocks; and explore the differences in jets from Class 0 to Class II sources. We performed a set of 2.5-dimensional magnetohydrodynamic numerical simulations that modelled supersonic jets ramming into a magnetized medium. The jet is formed by two components: a continously driven component that forms a quasi-stationary shock at the base of the jet; and a pulsed component constituted by blobs perturbing the shock. We explored different parameters for both components. We studied two cases: a jet less dense than the ambient medium (light jet), representing the case of HH 154; and a jet denser than the ambient (heavy jet), associated with DG Tau. We synthesized the count rate from the simulations and compared with available Chandra observations. Our model explains the formation of X-ray emitting quasi-stationary shocks observed at the base of jets in a natural way, being able to reproduce the observed jet properties at different evolutionary phases (in particular, for HH 154 and DG Tau). The jet is collimated by the magnetic field forming a quasi-stationary shock at the base which emits in X-rays even when perturbations formed by a train of blobs are present. We found similar collimation mechanisms dominating in both heavy and light jets...

Molecular clouds are predicted to emit nonthermal X-rays when they are close to particle-accelerating supernova remnants (SNRs), and the hard X-ray source AX J1714.1-3912, near the SNR RX J1713.7-3946, has long been considered a candidate for diffuse nonthermal emission associated with cosmic rays diffusing from the remnant to a closeby molecular cloud. We aim at ascertaining the nature of this source by analyzing two dedicated X-ray observations performed with Suzaku and Chandra. We extracted images from the data in various energy bands, spectra, and light curves and studied the long-term evolution of the X-ray emission on the basis of the ~4.5 yr time separation between the two observations. We found that there is no diffuse emission associated with AX J1714.1-3912, which is instead the point-like source CXOU J171343.9-391205. We discovered rapid time variability (timescale ~1 ks), together with a high intrinsic absorption and a hard nonthermal spectrum (power law with photon index Gamma~1.4). We also found that the X-ray flux of the source drops down by 1-2 orders of magnitude on a timescale of a few years. Our results suggest a possible association between AX J1714.1-3912 and a previously unknown supergiant fast X-ray transient, although further follow-up observations are necessary to prove this association definitively.

High resolution spectroscopy, providing constraints on plasma motions and temperatures, is a powerful means to investigate the structure of accretion streams in CTTS. In particular, the accretion shock region, where the accreting material is heated to temperatures of a few MK as it continues its inward bulk motion, can be probed by X-ray spectroscopy. To attempt to detect for the first time the motion of this X-ray-emitting post-shock material, we searched for a Doppler shift in the deep Chandra/HETGS observation of the CTTS TW Hya. This test should unveil the nature of this X-ray emitting plasma component in CTTS, and constrain the accretion stream geometry. We searched for a Doppler shift in the X-ray emission from TW Hya with two different methods, by measuring the position of a selected sample of emission lines, and by fitting the whole TW Hya X-ray spectrum, allowing the line-of-sight velocity to vary. We found that the plasma at T~2-4 MK has a line-of-sight velocity of 38.3+/-5.1 km/s with respect to the stellar photosphere. This result definitively confirms that this X-ray-emitting material originates in the post-shock region, at the base of the accretion stream, and not in coronal structures. The comparison of the observed velocity along the line of sight, 38.3+/-5.1 km/s, with the inferred intrinsic velocity of the post shock of TW Hya, v_post~110-120 km/s, indicates that the footpoints of the accretion streams on TW Hya are located at low latitudes on the stellar surface. Our results indicate that complex magnetic field geometries, such as that of TW Hya, permit low-latitude accretion spots. Moreover, since on TW Hya the redshift of the soft X-ray emission is very similar to that of the narrow component of the CIV resonance doublet at 1550 Ang, as found by Ardila et al. (2013), then the plasma at 2-4 MK and that at 0.1 MK likely originate in the same post-shock regions.

The Vela supernova remnant displays several ejecta, which are fragment-like features protruding beyond the front of its primary blast shock wave. They appear to be "shrapnel", bowshock-shaped relics of the supernova explosion. One of these pieces of shrapnel (A), located in the northeastern edge of the remnant, is peculiar because its X-ray spectrum exhibits a high Si abundance, in contrast to the other observed ejecta fragments, which show enhanced O, Ne, and Mg abundances. We present the analysis of another fragment located opposite to shrapnel A with respect to the center of the shell, in the southwestern boundary of the remnant, named shrapnel G. We thoroughly analyzed a dedicated XMM-Newton observation of shrapnel G by producing background-subtracted and exposure-corrected maps in different energy ranges, which we complemented with a spatially resolved spectral analysis of the X-ray emission. The fragment presents a bowshock-like shape with its anti-apex pointing to the center of the remnant. Its X-ray spectrum is best fit by a thermal plasma out of equilibrium of ionization with low O and Fe, roughly solar Ne and Mg, and a significantly high Si abundance, which is required to fit a very clear Si line at ~1.85 keV. Its chemical composition and spectral properties are compatible with those of shrapnel A, which is located on the opposite side of the remnant. As a consequence of the nucleosynthesis, pieces of Si-rich shrapnel are expected to originate in deeper layers of the progenitor star compared to ejecta with lower-Z elements. A high velocity and density contrast with respect to the surrounding ejecta are necessary to make shrapnel A and G overtake the forward shock. The line connecting shrapnel A and G crosses almost exactly the expansion center of the remnant, strongly suggesting a Si-rich jet-counterjet structure, reminiscent of that observed in the young remnant Cas A.

The acceleration times of the highest-energy particles which emit gamma-rays in young and middle-age SNRs are comparable with SNR age. If the number of particles starting acceleration was varying during early times after the supernova explosion then this variation should be reflected in the shape of the gamma-ray spectrum. We use the solution of the non-stationary equation for particle acceleration in order to analyze this effect. As a test case, we apply our method to describe gamma-rays from IC443. As a proxy of the IC443 parent supernova we consider SN1987A. First, we infer the time dependence of injection efficiency from evolution of the radio spectral index in SN1987A. Then, we use the inferred injection behavior to fit the gamma-ray spectrum of IC443. We show that the break in the proton spectrum needed to explain the gamma-ray emission is a natural consequence of the early variation of the cosmic ray injection, and that the very-high energy gamma-rays originate from particles which began acceleration during the first months after the supernova explosion. We conclude that the shape of the gamma-ray spectrum observed today in SNRs critically depends on the time variation of the cosmic ray injection process in the immediate post explosion phases. With the same model, we estimate also the possibility in the future to detect gamma-rays from SN 1987A.

Polarized radio emission has been mapped with great detail in several Galactic supernova remnants (SNRs), but has not yet been exploited to the extent it deserves. We have developed a method to model maps of the Stokes parameters for shell-like SNRs during their Sedov evolution phase. At first, 3-dimensional structure of a SNR has been computed, by modeling the distribution of the magnetohydrodynamic parameters and of the accelerated particles. The generation and dissipation of the turbulent component of magnetic field everywhere in SNR are also considered taking into account its interaction with accelerated particles. Then, in order to model the emission, we have used a generalization of the classical synchrotron theory, valid for the case in which the magnetic field has ordered and disordered components. Finally, 2-dimensional projected maps have been derived, for different orientations of SNR and of interstellar magnetic field with respect to the observer. An important effect to consider is the Faraday rotation of the polarization planes inside the SNR interior. In this paper we present details of the model, and describe general properties of the images.

Supernova remnants (SNRs) are diffuse extended sources characterized by a complex morphology and a non-uniform distribution of ejecta. Such a morphology reflects pristine structures and features of the progenitor supernova (SN) and the early interaction of the SN blast wave with the inhomogeneous circumstellar medium (CSM). Deciphering the observations of SNRs might open the possibility to investigate the physical properties of both the interacting ejecta and the shocked CSM. This requires accurate numerical models which describe the evolution from the SN explosion to the remnant development and which connect the emission properties of the remnants to the progenitor SNe. Here we show how multi-dimensional SN-SNR hydrodynamic models have been very effective in deciphering observations of SNR Cassiopeia A and SN 1987A, thus unveiling the structure of ejecta in the immediate aftermath of the SN explosion and constraining the 3D pre-supernova structure and geometry of the environment surrounding the progenitor SN.

The Galactic supernova remnant Kes 78 is surrounded by dense molecular clouds, whose projected position overlaps with the extended HESS gamma-ray source HESS J1852-000. The X-ray emission from the remnant has been recently revealed by Suzaku observations, which have shown indications for a hard X-ray component in the spectra, possibly associated with synchrotron radiation. We aim at describing the spatial distribution of the physical properties of the X-ray emitting plasma and at revealing the effects of the interaction of the remnant with the inhomogeneous ambient medium. We also aim at investigating the origin of the gamma-ray emission, which may be Inverse Compton radiation associated with X-ray synchrotron emitting electrons or hadronic emission originating from the impact of high energy protons on the nearby clouds. We analyzed an XMM-Newton EPIC observation of Kes 78 by performing image analysis and spatially resolved spectral analysis on a set of three regions. We tested our findings against the observations of the 12CO and 13CO emission in the environment of the remnant. We revealed the complex X-ray morphology of Kes 78 and found variations of the spectral properties of the plasma, with significantly denser and cooler material at the eastern edge of the remnant, which we interpret as a signature of interaction with a molecular cloud. We also exclude the presence of narrow filaments emitting X-ray synchrotron radiation. Assuming that the very high energy gamma-ray emission is associated with Kes 78, the lack of synchrotron emission rules out a leptonic origin. A hadronic origin is further supported by evidence for interaction of the remnant with a dense molecular cloud in its eastern limb.

The symbiotic nova V745 Sco was observed in outburst on 2014 February 6. Its observations by the Chandra X-ray Observatory at days 16 and 17 have revealed a spectrum characterized by asymmetric and blue-shifted emission lines. Here we investigate the origin of these asymmetries through three-dimensional hydrodynamic simulations describing the outburst during the first 20 days of evolution. The model takes into account thermal conduction and radiative cooling and assumes a blast wave propagates through an equatorial density enhancement. From the simulations, we synthesize the X-ray emission and derive the spectra as they would be observed with Chandra. We find that both the blast wave and the ejecta distribution are efficiently collimated in polar directions due to the presence of the equatorial density enhancement. The majority of the X-ray emission originates from the interaction of the blast with the equatorial density enhancement and is concentrated on the equatorial plane as a ring-like structure. Our "best-fit" model requires a mass of ejecta in the outburst $M_{ej} \approx 3\times 10^{-7}\,M_{\odot}$ and an explosion energy $E_b \approx 3 \times 10^{43}$ erg and reproduces the distribution of emission measure vs temperature and the evolution of shock velocity and temperature inferred from the observations. The model predicts asymmetric and blue-shifted line profiles similar to those observed and explains their origin as due to substantial X-ray absorption of red-shifted emission by ejecta material. The comparison of predicted and observed Ne and O spectral line ratios reveals no signs of strong Ne enhancement and suggests the progenitor is a CO white dwarf.

X-ray observations of protostellar jets show evidence of strong shocks heating the plasma up to temperatures of a few million degrees. In some cases, the shocked features appear to be stationary. They are interpreted as shock diamonds. We aim at investigating the physics that guides the formation of X-ray emitting stationary shocks in protostellar jets, the role of the magnetic field in determining the location, stability, and detectability in X-rays of these shocks, and the physical properties of the shocked plasma. We performed a set of 2.5-dimensional magnetohydrodynamic numerical simulations modelling supersonic jets ramming into a magnetized medium and explored different configurations of the magnetic field. The model takes into account the most relevant physical effects, namely thermal conduction and radiative losses. We compared the model results with observations, via the emission measure and the X-ray luminosity synthesized from the simulations. Our model explains the formation of X-ray emitting stationary shocks in a natural way. The magnetic field collimates the plasma at the base of the jet and forms there a magnetic nozzle. After an initial transient, the nozzle leads to the formation of a shock diamond at its exit which is stationary over the time covered by the simulations (~ 40 - 60 yr; comparable with time scales of the observations). The shock generates a point-like X-ray source located close to the base of the jet with luminosity comparable with that inferred from X-ray observations of protostellar jets. For the range of parameters explored, the evolution of the post-shock plasma is dominated by the radiative cooling, whereas the thermal conduction slightly affects the structure of the shock.

We used combined data from XMM-Newton and Chandra observatories to study the X-ray morphology of SNR G306.3-0.9. A spatially-resolved spectral analysis was used to obtain physical and geometrical parameters of different regions of the remnant. Spitzer infrared observations were also used to constrain the progenitor supernova and study the environment in which the SNR evolved. The X-ray morphology of the remnant displays a non-uniform structure of semi-circular appearance, with a bright southwest region and very weak or almost negligible X-ray emission in its northern part. These results indicate that the remnant is propagating in a non-uniform environment as the shock fronts are encountering a high-density medium, where enhanced infrared emission is detected. The X-ray spectral analysis of the selected regions shows distinct emission-line features of several metal elements, confirming the thermal origin of the emission. The X-ray spectra are well represented by a combination of two absorbed thermal plasma models: one in equilibrium ionization with a mean temperature of ~0.19 keV, and another out of equilibrium ionization at a higher temperature of ~1.1 or 1.6-1.9 keV. For regions located in the northeast, central, and southwest part of the SNR, we found elevated abundances of Si, S, Ar, Ca, and Fe, typical of ejecta material. The outer regions located northwest and south show values of the abundances above solar but lower than to those found in the central regions. This suggests that the composition of the emitting outer parts of the SNR is a combination of ejecta and shocked material of the interstellar medium. The comparison between the S/Si, Ar/Si, and Ca/Si abundances ratios (1.75, 1.27, and 2.72 in the central region, respectively), favor a Type Ia progenitor for this SNR, a result that is also supported by an independent morphological analysis using X-ray and IR data.

The supernova remnant SN 1006 is a source of high-energy particles and its southwestern limb is interacting with a dense ambient cloud, thus being a promising region for gamma-ray hadronic emission. We aim at describing the physics and the nonthermal emission associated with the shock-cloud interaction to derive the physical parameters of the cloud (poorly constrained by the data analysis), to ascertain the origin of the observed spatial variations in the spectral properties of the X-ray synchrotron emission, and to predict spectral and morphological features of the resulting gamma-ray emission. We performed 3-D magnetohydrodynamic simulations modeling the evolution of SN 1006 and its interaction with the ambient cloud, and explored different model setups. By applying the REMLIGHT code on the model results, we synthesized the synchrotron X-ray emission, and compared it with actual observations, to constrain the parameters of the model. We also synthesized the leptonic and hadronic gamma-ray emission from the models, deriving constraints on the energy content of the hadrons accelerated at the southwestern limb. We found that the impact of the SN 1006 shock front with a uniform cloud with density 0.5 cm^-3 can explain the observed morphology, the azimuthal variations of the cutoff frequency of the X-ray synchrotron emission, and the shock proper motion in the interaction region. Our results show that the current upper limit for the total hadronic energy in the southwestern limb is 2.5e49 erg.

Based on the XMM-Newton large program on SN1006 and our newly developed spatially resolved spectroscopy tools (Paper~I), we study the thermal emission from ISM and ejecta of SN1006 by analyzing the spectra extracted from 583 tessellated regions dominated by thermal emission. With some key improvements in spectral analysis as compared to Paper~I, we obtain much better spectral fitting results with less residuals. The spatial distributions of the thermal and ionization states of the ISM and ejecta show different features, which are consistent with a scenario that the ISM (ejecta) is heated and ionized by the forward (reverse) shock propagating outward (inward). Different elements have different spatial distributions and origins, with Ne mostly from the ISM, Si and S from the ejecta, and O and Mg from both ISM and ejecta. Fe L-shell lines are only detected in a small shell-like region SE to the center of SN1006, indicating that most of the Fe-rich ejecta has not yet or just recently been reached by the reverse shock. The overall ejecta abundance patterns for most of the heavy elements, except for Fe and sometimes S, are consistent with typical Type~Ia SN products. The NW half of the SNR interior probably represents a region with turbulently mixed ISM and ejecta, so has enhanced emission from O, Mg, Si, S, lower ejecta temperature, and a large diversity of ionization age. In addition to the asymmetric ISM distribution, an asymmetric explosion of the progenitor star is also needed to explain the asymmetric ejecta distribution.

The remnants of core-collapse supernovae (SNe) have complex morphologies that may reflect asymmetries and structures developed during the progenitor SN explosion. Here we investigate how the morphology of the SNR Cassiopeia A (Cas A) reflects the characteristics of the progenitor SN with the aim to derive the energies and masses of the post-explosion anisotropies responsible for the observed spatial distribution of Fe and Si/S. We model the evolution of Cas A from the immediate aftermath of the progenitor SN to the three-dimensional interaction of the remnant with the surrounding medium. The post-explosion structure of the ejecta is described by small-scale clumping of material and larger-scale anisotropies. The hydrodynamic multi-species simulations consider an appropriate post-explosion isotopic composition of the ejecta. The observed average expansion rate and shock velocities can be well reproduced by models with ejecta mass $M_{\rm ej}\approx 4M_{\odot}$ and explosion energy $E_{\rm SN}\approx 2.3\times 10^{51}$ erg. The post-explosion anisotropies (pistons) reproduce the observed distributions of Fe and Si/S if they had a total mass of $\approx 0.25\,M_{\odot}$ and a total kinetic energy of $\approx 1.5\times 10^{50}$ erg. The pistons produce a spatial inversion of ejecta layers at the epoch of Cas A, leading to the Si/S-rich ejecta physically interior to the Fe-rich ejecta. The pistons are also responsible for the development of bright rings of Si/S-rich material which form at the intersection between the reverse shock and the material accumulated around the pistons during their propagation. Our result supports the idea that the bulk of asymmetries observed in Cas A are intrinsic to the explosion.

Runaway O- and early B-type stars passing throughout the interstellar medium at supersonic velocities and characterized by strong stellar winds may produce bow shocks that can serve as particle acceleration sites. Previous theoretical models predict the production of high energy photons by non-thermal radiative processes, but their efficiency is still debated. We aim to test and explain the possibility of emission from the bow shocks formed by runaway stars traveling through the interstellar medium by using previous theoretical models. We apply our model to AE Aurigae, the first reported star with an X-ray detected bow shock, to BD+43 3654, in which the observations failed in detecting high energy emission, and to the transition phase of a supergiant star in the late stages of its life.From our analysis, we confirm that the X-ray emission from the bow shock produced by AE Aurigae can be explained by inverse Compton processes involving the infrared photons of the heated dust. We also predict low high energy flux emission from the bow shock produced by BD+43 3654, and the possibility of high energy emission from the bow shock formed by a supergiant star during the transition phase from blue to red supergiant.Bow shock formed by different type of runaway stars are revealed as a new possible source of high energy photons in our neighbourhood.

We present an X-ray analysis of the central region of supernova remnant (SNR) G332.5-5.6 through an exhaustive analysis of XMM-Netwon observations with complementary infrared observations. We characterize and discuss the origin of the observed X-ray morphology, which presents a peculiar plane edge over the west side of the central region. The morphology and spectral properties of the X-ray supernova remnant were studied using a single full frame XMM-Newton observation in the 0.3 to 10.0 keV energy band. Archival infrared WISE observations at 8, 12 and 24 \mu m were also used to investigate the properties of the source and its surroundings at different wavelengths. The results show that the extended X-ray emission is predominantly soft (0.3-1.2 keV) and peaks around 0.5 keV, which shows that it is an extremely soft SNR. X-ray emission correlates very well with central regions of bright radio emission. On the west side the radio/X-ray emission displays a plane-like feature with a terminal wall where strong infrared emission is detected. Our spatially resolved X-ray spectral analysis confirms that the emission is dominated by weak atomic emission lines of N, O, Ne, and Fe, all of them undetected in previous X-ray studies. These characteristics suggest that the X-ray emission is originated in an optically thin thermal plasma, whose radiation is well fitted by a non-equilibrium ionization collisional plasma (VNEI) X-ray emission model. Our study favors a scenario where G332.5-5.6 is expanding in a medium with an abrupt density change (the wall), likely a dense infrared emitting region of dust on the western side of the source.

Based on our newly developed methods and the XMM-Newton large program of SN1006, we extract and analyze the spectra from 3596 tessellated regions of this SNR each with 0.3-8 keV counts $>10^4$. For the first time, we map out multiple physical parameters, such as the temperature ($kT$), electron density ($n_e$), ionization parameter ($n_et$), ionization age ($t_{ion}$), metal abundances, as well as the radio-to-X-ray slope ($\alpha$) and cutoff frequency ($\nu_{cutoff}$) of the synchrotron emission. We construct probability distribution functions of $kT$ and $n_et$, and model them with several Gaussians, in order to characterize the average thermal and ionization states of such an extended source. We construct equivalent width (EW) maps based on continuum interpolation with the spectral model of each regions. We then compare the EW maps of OVII, OVIII, OVII K$\delta-\zeta$, Ne, Mg, SiXIII, SiXIV, and S lines constructed with this method to those constructed with linear interpolation. We further extract spectra from larger regions to confirm the features revealed by parameter and EW maps, which are often not directly detectable on X-ray intensity images. For example, O abundance is consistent with solar across the SNR, except for a low-abundance hole in the center. This "O Hole" has enhanced OVII K$\delta-\zeta$ and Fe emissions, indicating recently reverse shocked ejecta, but also has the highest $n_et$, indicating forward shocked ISM. Therefore, a multi-temperature model is needed to decompose these components. The asymmetric metal distributions suggest there is either an asymmetric explosion of the SN or an asymmetric distribution of the ISM.

The emission of supernova remnants reflects the properties of both the progenitor supernovae and the surrounding environment. The complex morphology of the remnants, however, hampers the disentanglement of the two contributions. Here we aim at identifying the imprint of SN 1987A on the X-ray emission of its remnant and at constraining the structure of the environment surrounding the supernova. We performed high-resolution hydrodynamic simulations describing SN 1987A soon after the core-collapse and the following three-dimensional expansion of its remnant between days 1 and 15000 after the supernova. We demonstrated that the physical model reproducing the main observables of SN 1987A during the first 250 days of evolution reproduces also the X-ray emission of the subsequent expanding remnant, thus bridging the gap between supernovae and supernova remnants. By comparing model results with observations, we constrained the explosion energy in the range $1.2-1.4\times 10^{51}$~erg and the envelope mass in the range $15-17 M_{\odot}$. We found that the shape of X-ray lightcurves and spectra at early epochs (<15 years) reflects the structure of outer ejecta: our model reproduces the observations if the outermost ejecta have a post-explosion radial profile of density approximated by a power law with index $\alpha = -8$. At later epochs, the shapes of X-ray lightcurves and spectra reflect the density structure of the nebula around SN 1987A. This enabled us to ascertain the origin of the multi-thermal X-ray emission, to disentangle the imprint of the supernova on the remnant emission from the effects of the remnant interaction with the environment, and to constrain the pre-supernova structure of the nebula.

Young supernova remnants show a characteristic ejecta-dominated X-ray emission that allows us to probe the products of the explosive nucleosynthesis processes and to ascertain important information about the physics of the supernova explosions. Hard X-ray observations have recently revealed the radioactive decay lines of 44Ti at ~67.9 keV and ~78.4 keV in the Tycho's SNR. We here analyze the set of XMM-Newton archive observations of the Tycho's SNR. We produce equivalent width maps of the Fe K and Ca XIX emission lines and find indications for a stratification of the abundances of these elements and significant anisotropies. We then perform a spatially resolved spectral analysis by identifying five different regions characterized by high/low values of the Fe K equivalent width. We find that the spatial distribution of the Fe K emission is correlated with that of the Cr XXII. We also detect the Ti K-line complex in the spectra extracted from the two regions with the highest values of the Fe and Cr equivalent widths. The Ti line emissions remains undetected in regions where the Fe and Cr equivalent widths are low. Our results indicate that the post-shock Ti is spatially co-located with other iron-peak nuclei in Tycho's SNR, in agreement with the predictions of multi-D models of Type Ia supernovae.

We investigate the plausibility of detecting X-ray emission from a stellar jet that impacts against a dense molecular cloud. This scenario may be usual for classical T Tauri stars with jets in dense star-forming complexes. We first model the impact of a jet against a dense cloud by 2D axisymmetric hydrodynamic simulations, exploring different configurations of the ambient environment. Then, we compare our results with XMM-Newton observations of the Herbig-Haro object HH 248, where extended X-ray emission aligned with the optical knots is detected at the edge of the nearby IC 434 cloud. Our simulations show that a jet can produce plasma with temperatures up to 10 MK, consistent with production of X-ray emission, after impacting a dense cloud. We find that jets denser than the ambient medium but less dense than the cloud produce detectable X-ray emission only at the impact onto the cloud. From the exploration of the model parameter space, we constrain the physical conditions (jet density and velocity, cloud density) that reproduce well the intrinsic luminosity and emission measure of the X-ray source possibly associated with HH 248. Thus, we suggest that the extended X-ray source close to HH 248 corresponds to the jet impacting on a dense cloud.

We investigate the properties of X-ray emission from accretion shocks in classical T Tauri stars (CTTSs), generated where the infalling material impacts the stellar surface. Both observations and models of the accretion process reveal several aspects that are still unclear: the observed X-ray luminosity in accretion shocks is below the predicted value, and the density versus temperature structure of the shocked plasma, with increasing densities at higher temperature, deduced from the observations, is at odds with that proposed in the current picture of accretion shocks. To address these open issues we investigate whether a correct treatment of the local absorption by the surrounding medium is crucial to explain the observations. To this end, we describe the impact of an accretion stream on a CTTS by considering a magnetohydrodynamic model. From the model results we synthesize the X-ray emission from the accretion shock by producing maps and spectra. We perform density and temperature diagnostics on the synthetic spectra, and we directly compare the results with the observations. Our model shows that the X-ray fluxes inferred from the emerging spectra are lower than expected because of the complex local absorption by the optically thick material of the chromosphere and of the unperturbed stream. Moreover, our model including the effects of local absorption explains in a natural way the apparently puzzling pattern of density versus temperature observed in the X-ray emission from accretion shocks.

The supernova remnant SN 1006 is a powerful source of high-energy particles and evolves in a relatively tenuous and uniform environment, though interacting with an atomic cloud in its northwestern limb. The X-ray image of SN 1006 reveals an indentation in the southwestern part of the shock front and the HI maps show an isolated cloud (southwestern cloud) having the same velocity as the northwestern cloud and whose morphology fits perfectly in the indentation. We performed spatially resolved spectral analysis of a set of small regions in the southwestern nonthermal limb and studied the deep X-ray spectra obtained within the XMM-Newton SN 1006 Large Program. We also analyzed archive HI data, obtained combining single dish and interferometric observations. We found that the best-fit value of the N_H derived from the X-ray spectra significantly increases in regions corresponding to the southwestern cloud, while the cutoff energy of the synchrotron emission decreases. The amount of the N_H variations corresponds perfectly with the HI column density of the southwestern cloud, as measured from the radio data. The decrease in the cutoff energy at the indentation clearly reveals that the back side of the cloud is actually interacting with the remnant. The southwestern limb therefore presents a unique combination of efficient particle acceleration and high ambient density, thus being the most promising region for gamma-ray hadronic emission in SN 1006. We estimate that such emission will be detectable with the Fermi telescope within a few years.

(abridged) AIMS. We investigate the dynamics and stability of post-shock plasma streaming along nonuniform stellar magnetic fields at the impact region of accretion columns. We study how the magnetic field configuration and strength determine the structure, geometry, and location of the shock-heated plasma. METHODS. We model the impact of an accretion stream onto the chromosphere of a CTTS by 2D axisymmetric magnetohydrodynamic simulations. Our model takes into account the gravity, the radiative cooling, and the magnetic-field-oriented thermal conduction. RESULTS. The structure, stability, and location of the shocked plasma strongly depend on the configuration and strength of the magnetic field. For weak magnetic fields, a large component of B may develop perpendicular to the stream at the base of the accretion column, limiting the sinking of the shocked plasma into the chromosphere. An envelope of dense and cold chromospheric material may also develop around the shocked column. For strong magnetic fields, the field configuration determines the position of the shock and its stand-off height. If the field is strongly tapered close to the chromosphere, an oblique shock may form well above the stellar surface. In general, a nonuniform magnetic field makes the distribution of emission measure vs. temperature of the shocked plasma lower than in the case of uniform magnetic field. CONCLUSIONS. The initial strength and configuration of the magnetic field in the impact region of the stream are expected to influence the chromospheric absorption and, therefore, the observability of the shock-heated plasma in the X-ray band. The field strength and configuration influence also the energy balance of the shocked plasma, its emission measure at T > 1 MK being lower than expected for a uniform field. The above effects contribute in underestimating the mass accretion rates derived in the X-ray band.

Protostellar jets appear at all stages of star formation when the accretion process is still at work. Jets travel at velocities of hundreds of km/s, creating strong shocks when interacting with interstellar medium. Several cases of jets have been detected in X-rays, typically showing soft emission. For the first time, we report evidence of hard X-ray emission possibly related to non-thermal processes not explained by previous models of the post-shock emission predicted in the jet/ambient interaction scenario. HH 80 is located at the south head of the jet associated to the massive protostar IRAS 18162-2048. It shows soft and hard X-ray emission in regions that are spatially separated, with the soft X-ray emission region situated behind the region of hard X-ray emission. We propose a scenario for HH 80 where soft X-ray emission is associated to thermal processes from the interaction of the jet with denser ambient matter and the hard X-ray emission is produced by synchrotron radiation at the front shock.

The spectral shape of the synchrotron X-ray emission from SN 1006 reveals the fundamental role played by radiative losses in shaping the high-energy tail of the electron spectrum. We analyze data from the XMM-Newton SN 1006 Large Program and confirm that in both nonthermal limbs the loss-limited model correctly describes the observed spectra. We study the physical origin of the observed variations of the synchrotron cutoff energy across the shell. We investigate the role played by the shock velocity and by the electron gyrofactor. We found that the cutoff energy of the syncrotron X-ray emission reaches its maximum value in regions where the shock has experienced its highest average speed. This result is consistent with the loss-limited framework. We also find that the electron acceleration in both nonthermal limbs of SN 1006 proceeds close to the Bohm diffusion limit, the gyrofactor being in the range 1.5-4. We finally investigate possible explanations for the low values of cutoff energy measured in thermal limbs.

Synchrotron X-ray emission from the rims of young supernova remnants allows us to study the high-energy tail of the electrons accelerated at the shock front. The analysis of X-ray spectra can provide information on the physical mechanisms that limit the energy achieved by the electrons in the acceleration process. We aim at verifying whether the maximum electron energy in SN 1006 is limited by synchrotron losses and at obtaining information on the shape of the cutoff in the X-ray synchrotron emission. We analyzed the deep observations of the XMM-Newton SN 1006 Large Program. We performed spatially resolved spectral analysis of a set of small regions in the nonthermal limbs and studied the X-ray spectra by adopting models that assume different electron spectra. We found out that a loss-limited model provides the best fit to all the spectra and this indicates that the shape of the cutoff in the electron momentum (p) distribution has the form exp[-(p/p_cut)^2]. We also detected residual thermal emission from shocked ambient medium and confirmed the reliability of previous estimates of the post-shock density. Our results indicate that radiative losses play a fundamental role in shaping the electron spectrum in SN 1006.

The study of both supernova remnants and the hot and cold phases of the interstellar medium are essential for understanding the final stages of stellar evolution and their feedback on the evolution of galaxies through injection of energy and heavy elements. These studies are also crucial for understanding the physics of supernovae, their cosmological implication, and the origin of galactic cosmic rays. The unique capabilities of Athena+ will allow us to explore a new parameter space. Spatially-resolved high-resolution spectroscopy using Athena+ X-IFU of young remnants will allow to characterize individual parcels of ejected material in the line of sight in terms of kinematics, ionization and composition, providing access to the three dimensional geometry of the explosion. Athena+ will also allow studying shock physics and particle acceleration in supernova remnants, as well as their interaction with their environment. Athena+ X-IFU will also characterize the ionization mechanisms competing in forming the complex structures of the hot interstellar medium, likely to keep the echo of past star formation activity, both in our Galaxy and nearby ones. For the first time the dust and gas of the densest cold medium, like in the Galactic Centre environment, will be studied. Athena+ X-IFU will observe, along with the Mg K and Si K edges, which are the main tracers of the silicates content of the ISM, the Fe K edge with unprecedented sensitivity and energy-resolution. This will allow us to study for the first time the nature of Fe-bearing dust in such regions.

This White Paper, submitted to the recent ESA call for science themes to define its future large missions, advocates the need for a transformational leap in our understanding of two key questions in astrophysics: 1) How does ordinary matter assemble into the large scale structures that we see today? 2) How do black holes grow and shape the Universe? Hot gas in clusters, groups and the intergalactic medium dominates the baryonic content of the local Universe. To understand the astrophysical processes responsible for the formation and assembly of these large structures, it is necessary to measure their physical properties and evolution. This requires spatially resolved X-ray spectroscopy with a factor 10 increase in both telescope throughput and spatial resolving power compared to currently planned facilities. Feedback from supermassive black holes is an essential ingredient in this process and in most galaxy evolution models, but it is not well understood. X-ray observations can uniquely reveal the mechanisms launching winds close to black holes and determine the coupling of the energy and matter flows on larger scales. Due to the effects of feedback, a complete understanding of galaxy evolution requires knowledge of the obscured growth of supermassive black holes through cosmic time, out to the redshifts where the first galaxies form. X-ray emission is the most reliable way to reveal accreting black holes, but deep survey speed must improve by a factor ~100 over current facilities to perform a full census into the early Universe. The Advanced Telescope for High Energy Astrophysics (Athena+) mission provides the necessary performance (e.g. angular resolution, spectral resolution, survey grasp) to address these questions and revolutionize our understanding of the Hot and Energetic Universe. These capabilities will also provide a powerful observatory to be used in all areas of astrophysics.

Many supernova remnants (SNRs) are characterized by a knotty ejecta structure. The Vela SNR is an excellent example of remnant in which detached clumps of ejecta are visible as X-ray emitting bullets that have been observed and studied in great detail. We aim at modelling the evolution of ejecta shrapnel in the Vela SNR, investigating the role of their initial parameters (position and density) and addressing the effects of thermal conduction and radiative losses. We performed a set of 2-D hydrodynamic simulations describing the evolution of a density inhomogeneity in the ejecta profile. We explored different initial setups. We found that the final position of the shrapnel is very sensitive to its initial position within the ejecta, while the dependence on the initial density contrast is weaker. Our model also shows that moderately overdense knots can reproduce the detached features observed in the Vela SNR. Efficient thermal conduction produces detectable effects by determining an efficient mixing of the ejecta knot with the surrounding medium and shaping a characteristic elongated morphology in the clump.

Aims: We want to probe the physics of fast collision-less shocks in supernova remnants. In particular, we are interested in the non-equilibration of temperatures and particle acceleration. Specifically, we aim to measure the oxygen temperature with regards to the electron temperature. In addition, we search for synchrotron emission in the northwestern thermal rim. Methods: This study is part of a dedicated deep observational project of SN 1006 using XMM-Newton, which provides us with currently the best resolution spectra of the bright northwestern oxygen knot. We aim to use the reflection grating spectrometer to measure the thermal broadening of the O vii line triplet by convolving the emission profile of the remnant with the response matrix. Results: The line broadening was measured to be {\sigma}_e = 2.4 \pm 0.3 eV, corresponding to an oxygen temperature of 275$^{+72}_{-63}$ keV. From the EPIC spectra we obtain an electron temperature of 1.35 \pm 0.10 keV. The difference in temperature between the species provides further evidence of non-equilibration of temperatures in a shock. In addition, we find evidence for a bow shock that emits X-ray synchrotron radiation, which is at odds with the general idea that due to the magnetic field orientation only in the NE and SW region X-ray synchrotron radiation should be emitted. We find an unusual H{\alpha} and X-ray synchrotron geometry, in that the H{\alpha} emission peaks downstream of the synchrotron emission. This may be an indication for a peculiar H{\alpha} shock, in which the density is lower and neutral fraction are higher than in other supernova remnants, resulting in a peak in H{\alpha} emission further downstream of the shock.

This contribution summarizes the splinter session "Non-thermal processes in coronae and beyond" held at the Cool Stars 17 workshop in Barcelona in 2012. It covers new developments in high energy non-thermal effects in the Earth's exosphere, solar and stellar flares, the diffuse emission in star forming regions and reviews the state and the challenges of the underlying atomic databases.

Runaway stars produce shocks when passing through interstellar medium at supersonic velocities. Bow shocks have been detected in the mid-infrared for several high-mass runaway stars and in radio waves for one star. Theoretical models predict the production of high-energy photons by non-thermal radiative processes in a number sufficiently large to be detected in X-rays. To date, no stellar bow shock has been detected at such energies. We present the first detection of X-ray emission from a bow shock produced by a runaway star. The star is AE Aur, which was likely expelled from its birthplace by the encounter of two massive binary systems and now is passing through the dense nebula IC 405. The X-ray emission from the bow shock is detected at 30" to the northeast of the star, coinciding with an enhancement in the density of the nebula. From the analysis of the observed X-ray spectrum of the source and our theoretical emission model, we confirm that the X-ray emission is produced mainly by inverse Compton upscattering of infrared photons from dust in the shock front.

Shock fronts in young supernova remnants are the best candidates for being sites of cosmic ray acceleration up to a few PeV, though conclusive experimental evidence is still lacking. Hadron acceleration is expected to increase the shock compression ratio, providing higher postshock densities, but X-ray emission from shocked ambient medium has not firmly been detected yet in remnants where particle acceleration is at work. We exploited the deep observations of the XMM-Newton Large Program on SN 1006 to verify this prediction. We performed spatially resolved spectral analysis of a set of regions covering the southeastern rim of SN 1006. We studied the spatial distribution of the thermodynamic properties of the ambient medium and carefully verified the robustness of the result with respect to the analysis method. We detected the contribution of the shocked ambient medium. We also found that the postshock density of the interstellar medium significantly increases in regions where particle acceleration is efficient. Under the assumption of uniform preshock density, we found that the shock compression ratio reaches a value of ~6 in regions near the nonthermal limbs. Our results support the predictions of shock modification theory and indicate that effects of acceleration of cosmic ray hadrons on the postshock plasma can be observed in supernova remnants.

The detection of very hot plasma in the quiescent corona is important for diagnosing heating mechanisms. The presence and the amount of such hot plasma is currently debated. The SphinX instrument on-board CORONAS-PHOTON mission is sensitive to X-ray emission well above 1 keV and provides the opportunity to detect the hot plasma component. We analyzed the X-ray spectra of the solar corona collected by the SphinX spectrometer in May 2009 (when two active regions were present). We modelled the spectrum extracted from the whole Sun over a time window of 17 days in the 1.34-7 keV energy band by adopting the latest release of the APED database. The SphinX broadband spectrum cannot be modelled by a single isothermal component of optically thin plasma and two components are necessary. In particular, the high statistics and the accurate calibration of the spectrometer allowed us to detect a very hot component at ~7 million K with an emission measure of ~2.7 x 10^44 cm^-3. The X-ray emission from the hot plasma dominates the solar X-ray spectrum above 4 keV. We checked that this hot component is invariably present both at high and low emission regimes, i.e. even excluding resolvable microflares. We also present and discuss a possible non-thermal origin (compatible with a weak contribution from thick-target bremsstrahlung) for this hard emission component. Our results support the nanoflare scenario and might confirm that a minor flaring activity is ever-present in the quiescent corona, as also inferred for the coronae of other stars.

Recent X-ray observations of the supernova remnant IC443 interacting with molecular clouds have shown the presence of a new population of hard X-ray sources related to the remnant itself, which has been interpreted in terms of fast ejecta fragment propagating inside the dense environment. Prompted by these studies, we have obtained a deep {\sl XMM-Newton} observation of the supernova remnant (SNR) Kes 69, which also shows signs of shock-cloud interaction. We report on the detection of 18 hard X-ray sources in the field of Kes 69, a significant excess of the expected galactic source population in the field, spatially correlated with CO emission from the cloud in the remnant environment. The spectra of 3 of the 18 sources can be described as hard power laws with photon index <2 plus line emission associated to K-shell transitions. We discuss the two most promising scenarios for the interpretation of the sources, namely fast ejecta fragments (as in IC443) and cataclysmic variables. While most of the observational evidences are consistent with the former interpretation, we cannot rule out the latter.

We investigate the role played by initial clumping of ejecta and by efficient acceleration of cosmic rays (CRs) in determining the density structure of the post-shock region of a Type Ia supernova remnant (SNR) through detailed 3D MHD modeling. Our model describes the expansion of a SNR through a magnetized interstellar medium (ISM), including the initial clumping of ejecta and the effects on shock dynamics due to back-reaction of accelerated CRs. The model predictions are compared to the observations of SN 1006. We found that the back-reaction of accelerated CRs alone cannot reproduce the observed separation between the forward shock (FS) and the contact discontinuity (CD) unless the energy losses through CR acceleration and escape are very large and independent of the obliquity angle. On the contrary, the clumping of ejecta can naturally reproduce the observed small separation and the occurrence of protrusions observed in SN 1006, even without the need of accelerated CRs. We conclude that FS-CD separation is a probe of the ejecta structure at the time of explosion rather than a probe of the efficiency of CR acceleration in young SNRs.

The heating of the solar corona is one of the big questions in astrophysics. Rapid pulses called nanoflares are among the best candidate mechanisms. The analysis of the time variability of coronal X-ray emission is potentially a very useful tool to detect impulsive events. We analyze the small-scale variability of a solar active region in a high cadence Hinode/XRT observation. The dataset allows us to detect very small deviations of emission fluctuations from the distribution expected for a constant rate. We discuss the deviations in the light of the pulsed-heating scenario.

The thermal structure of the post-shock region of a young supernova remnant (SNR) is heavily affected by two main physical effects, the back-reaction of accelerated cosmic rays (CRs) and the Rayleigh-Taylor (RT) instabilities developing at the contact discontinuity between the ejecta and the shocked interstellar medium (ISM). Here, we investigate the role played by both physical mechanisms in the evolution of SNRs through detailed 3D MHD modeling. Our model describes the expansion of the remnant through a magnetized ISM, including consistently the initial ejecta clumping and the effects on shock dynamics due to back-reaction of accelerated CRs. We discuss the role of the initial ejecta clumpiness in developing strong instabilities at the contact discontinuity which may extend upstream to the main shock and beyond.

Mixed-morphology SNRs are characterized by a shell-like radio emission, a centrally peaked X-ray morphology, and by interaction with molecular clouds. Many models have been proposed to explain these peculiar remnants, but their physical origin is still unclear. The recent discovery of over-ionized (i. e. recombining) ejecta in 3 mixed-morphology SNRs has dramatically challenged all the previous models and opened up new, unexpected scenarios. I review the main properties of these remnants and their peculiar X-ray spectral properties. I also discuss the hydrodynamic model developed to explain the presence of over-ionized ejecta in W49B and present a list of open issues that still need to be clarified.

X-ray emission from about ten protostellar jets has been discovered and it appears as a feature common to the most energetic jets. Although X-ray emission seems to originate from shocks internal to jets, the mechanism forming these shocks remains controversial. One of the best studied X-ray jet is HH 154 that has been observed by Chandra over a time base of about 10 years. We analyze the Chandra observations of HH 154 by investigating the evolution of its X-ray source. We show that the X-ray emission consists of a bright stationary component and a faint elongated component. We interpret the observations by developing a hydrodynamic model describing a protostellar jet originating from a nozzle and compare the X-ray emission synthesized from the model with the X-ray observations. The model takes into account the thermal conduction and radiative losses and shows that the jet/nozzle leads to the formation of a diamond shock at the nozzle exit. The shock is stationary over the period covered by our simulations and generates an X-ray source with luminosity and spectral characteristics in excellent agreement with the observations. We conclude that the X-ray emission from HH 154 is consistent with a diamond shock originating from a nozzle through which the jet is launched into the ambient medium. We suggest that the physical origin of the nozzle could be related to the dense gas in which the HH 154 driving source is embedded and/or to the magnetic field at the jet launching/collimation region.

The synchrotron radio morphology of bilateral supernova remnants depends on the mechanisms of particle acceleration and on the viewing geometry. However, unlike X-ray and $\gamma$-ray morphologies, the radio emission does not depend on the cut-off region of the parent electron population, making it a simpler and more straightforward tool to investigate the physics of cosmic ray production in SNRs. Our aim is to derive from the radio morphology tight constraints on the direction of the local magnetic field and its gradient, and on the obliquity dependence of the electron injection efficiency. We perform a set of 3D MHD simulations describing the expansion of a spherical SNR through a magnetized medium with a non-uniform ISMF. From the simulations, we derive non-thermal radio maps and compare them with observations of the SN1006 remnant. We find that the radio morphology of SN1006 at 1 GHz is best-fitted by a model with quasi-parallel injection efficiency, a magnetic field aspect angle of $38^\circ\pm 4^\circ$ with the line of sight, and a gradient of the field strength toward the galactic plane, higher then the expected variations of the large scale field of the Galaxy. We conclude that the radio limbs of SN1006 are polar caps, not lying in the plane of sky}. The study of the synchrotron radio emission of SNRs is of crucial importance to derive information on the galactic magnetic field in the vicinity of the remnants, and more hints on the actual injection efficiency scenario.

It is generally agreed that small impulsive energy bursts called nanoflares are responsible for at least some of the Sun's hot corona, but whether they are the explanation for most of the multi-million degree plasma has been a matter of ongoing debate. We here present evidence that nanoflares are widespread in an active region observed by the X-Ray Telescope on-board the Hinode mission. The distributions of intensity fluctuations have small but important asymmetries, whether taken from individual pixels, multi-pixel subregions, or the entire active region. Negative fluctuations (corresponding to reduced intensity) are greater in number but weaker in amplitude, so that the median fluctuation is negative compared to a mean of zero. Using MonteCarlo simulations, we show that only part of this asymmetry can be explained by Poisson photon statistics. The remainder is explainable with a tendency for exponentially decreasing intensity, such as would be expected from a cooling plasma produced from a nanoflare. We suggest that nanoflares are a universal heating process within active regions.

Vela X is a large, 3x2 degrees, radio-emitting pulsar wind nebula (PWN) powered by the Vela pulsar in the Vela supernova remnant. Using four Suzaku/XIS observations pointed just outside Vela X, we find hard X-ray emission extending throughout the fields of view. The hard X-ray spectra are well represented by a power-law. The photon index is measured to be constant at Gamma~2.4, similar to that of the southern outer part of Vela X. The power-law flux decreases with increasing distance from the pulsar. These properties lead us to propose that the hard X-ray emission is associated with the Vela PWN. The larger X-ray extension found in this work strongly suggests that distinct populations relativistic electrons form the X-ray PWN and Vela X, as was recently inferred from multiwavelength spectral modeling of Vela X.

W49B is a mixed-morphology supernova remnant with thermal X-ray emission dominated by the ejecta. In this remnant, the presence of overionized plasma has been directly established, with information about its spatial structure. However, the physical origin of the overionized plasma in W49B has not yet been understood. We investigate this intriguing issue through a 2D hydrodynamic model that takes into account, for the first time, the mixing of ejecta with the inhomogeneous circumstellar and interstellar medium, the thermal conduction, the radiative losses from optically thin plasma, and the deviations from equilibrium of ionization induced by plasma dynamics. The model was set up on the basis of the observational results. We found that the thermal conduction plays an important role in the evolution of W49B, in- ducing the evaporation of the circumstellar ring-like cloud (whose presence has been deduced from previous observations) that mingles with the surrounding hot medium, cooling down the shocked plasma, and pushes the ejecta backwards to the center of the remnant, forming there a jet-like structure. During the evolution, a large region of overionized plasma forms within the remnant. The overionized plasma originates from the rapid cooling of the hot plasma originally heated by the shock reflected from the dense ring-like cloud. In particular, we found two different ways for the rapid cooling of plasma to appear: i) the mixing of relatively cold and dense material evaporated from the ring with the hot shocked plasma and ii) the rapid adiabatic expansion of the ejecta. The spatial distribution of the radiative recombination continuum predicted by the numerical model is in good agreement with that observed.

Experimental spectra and images of the supernova remnant SN1006 have been reported for radio, X-ray and TeV gamma-ray bands. Several comparisons between models and observations have been discussed in the literature, showing that the broad-band spectrum from the whole remnant as well as a sharpest radial profile of the X-ray brightness can be both fitted by adopting a model of SN1006 which strongly depends on the non-linear effects of the accelerated cosmic rays; these models predict post-shock magnetic field (MF) strengths of the order of 150 micro G. Here we present a new way to compare models and observations, in order to put constraints on the physical parameters and mechanisms governing the remnant. In particular, we show that a simple model based on the classic MHD and cosmic rays acceleration theories allows us to investigate the spatially distributed characteristics of SN1006 and to put observational constraints on the kinetics and MF. Our method includes modelling and comparison of the azimuthal and radial profiles of the surface brightness in radio, hard X-rays and TeV gamma-rays as well as the azimuthal variations of the electron maximum energy. In addition, this simple model also provides good fits to the radio-to-gamma-ray spectrum of SN1006. We find that our best-fit model predicts an effective MF strength inside SN1006 of 32 micro G, in good agreement with the `leptonic' model suggested by the HESS Collaboration (2010). Finally, some difficulties in both the `classic' and the non-linear models are discussed. A number of evidences about non-uniformity of MF around SN1006 are noted.

Observations of SNRs in X-ray and gamma-ray bands promise to contribute with important information in our understanding on the nature of galactic cosmic rays. The analysis of SNRs images collected in different energy bands requires the support of theoretical modeling of synchrotron and inverse Compton (IC) emission. We develop a numerical code (REMLIGHT) to synthesize, from MHD simulations, the synchrotron radio, X-ray and IC gamma-ray emission from SNRs expanding in non-uniform interstellar medium (ISM) and/or non-uniform interstellar magnetic field (ISMF). As a first application, the code is used to investigate the effects of non-uniform ISMF on the SNR morphology in the non-thermal X-ray and gamma-ray bands. We perform 3D MHD simulations of a spherical SNR shock expanding through a magnetized ISM with a gradient of ambient magnetic field strength. The model includes an approximate treatment of upstream magnetic field amplification and the effect of shock modification due to back reaction of accelerated cosmic rays. From the simulations, we synthesize the synchrotron radio, X-ray and IC gamma-ray emission with REMLIGHT, making different assumptions about the details of acceleration and injection of relativistic electrons. A gradient of the ambient magnetic field strength induces asymmetric morphologies in radio, X-ray and gamma-ray bands independently from the model of electron injection if the gradient has a component perpendicular to the line-of-sight. The degree of asymmetry of the remnant morphology depends on the details of the electron injection and acceleration and is different in the radio, X-ray, and gamma-ray bands. The non-thermal X-ray morphology is the most sensitive to the gradient, showing the highest degree of asymmetry. The IC gamma-ray emission is weakly sensitive to the non-uniform ISMF, the degree of asymmetry of the SNR morphology being the lowest in this band.

Protostellar jets are known to emit in a wide range of bands, from radio to IR to optical bands, and to date also about ten X-ray emitting jets have been detected, with a rate of discovery of about one per year. We aim at investigating the mechanism leading to the X-ray emission detected in protostellar jets and at constraining the physical parameters that describe the jet/ambient interaction by comparing our model predictions with observations. We perform 2D axisymmetric hydrodynamic simulations of the interaction between a supersonic jet and the ambient. The jet is described as a train of plasma blobs randomly ejected by the stellar source along the jet axis. We explore the parameter space by varying the ejection rate, the initial jet Mach number, and the initial density contrast between the ambient and the jet. We synthesized from the model the X-ray emission as it would be observed with the current X-ray telescopes. The mutual interactions among the ejected blobs and of the blobs with the ambient medium lead to complex X-ray emitting structures within the jet: irregular chains of knots; isolated knots with measurable proper motion; apparently stationary knots; reverse shocks. The predicted X-ray luminosity strongly depends on the ejection rate and on the initial density contrast between the ambient and the jet, with a weaker dependence on the jet Mach number. Our model represents the first attempt to describe the X-ray properties of all the X-ray emitting protostellar jets. The comparison between our model predictions and the observations can provide a useful diagnostic tool necessary for a proper interpretation of the observations. In particular, we suggest that the observable quantities derived from the spectral analysis of X-ray observations can be used to constrain the ejection rate, a parameter explored in our model that is not measurable by current observations.

Recent Suzaku X-ray observations of the ejecta-dominated supernova remnant W49B have shown that in the global spectrum there is a clear indication for the presence of overionized plasma whose physical origin is still under debate. In order to ascertain the physical origin of such a rapidly cooling plasma, we focus on the study of its spatial localization within the X-ray emitting ejecta. We confirm the presence of a saw-edged excess (interpreted as a strong radiative recombination continuum) in the global spectrum above 8 keV, emerging above the ionization-equilibrium model. We produce a hardness ratio map to determine where the plasma is overionized and we perform a spectral analysis of the regions with and without strong overionization. We find that the overionized plasma is localized in the center of the remnant and in its western jet, while it is not detected in the bright eastern jet, where the expansion of the ejecta is hampered by their interaction with a dense interstellar cloud. The location of overionized plasma suggests that the inner ejecta are rapidly cooling by expansion, unlike the outer ejecta, for which expansion is hampered by interstellar clouds seen in H2

The Cygnus Loop is interacting with a protrusion of the cavity wall in its eastern edge (the XA region), where the X-ray emission is very bright. The complexity of the environment and the non-linear physical processes of the shock-cloud interaction make the origin of the X-ray emission still not well understood. Our purpose is to understand the physical origin of the X-ray emission in the XA region, addressing, in particular, the role of thermal conduction in the interaction process. We analyzed two XMM-Newton data sets, performing image analysis and spatially resolved spectral analysis on a set of homogeneous regions. We applied a recently developed diagnostic tool to compare spectral analysis results with predictions of theoretical models, and to estimate the efficiency of thermal conduction on the X-ray emitting shocked plasma. We found that the inhomogeneous cavity wall contains both large clumps (the protrusion) and small isolated clumps with different densities. A large indentation bent over to the south is detected. The abundance of the surrounding ISM is ~0.2 times solar value. We confirmed the important role of thermal conduction in the evolution of X-ray emitting plasma during shock-cloud interaction.

X-ray emitting features originating from the interaction of supernova shock waves with small interstellar gas clouds are revealed in many X-ray observations of evolved supernova remnants (e.g. Cygnus Loop and Vela), but their interpretation is not straightforward. We develop a self-consistent method for the analysis and interpretation of shock-cloud interactions in middle-aged supernova remnants, which can provide the key parameters of the system and the role of relevant physical effects like the thermal conduction, without the need to run ad-hoc numerical simulations and to bother of morphology details. We explore all the possible values of the shock speed and cloud density contrast relevant to middle-aged SNRs with a set of hydrodynamic simulations of shock-cloud interaction, including the effects of thermal conduction and radiative cooling. From the simulations, we synthesize spatially and spectrally resolved focal-plane data as they would be collected with XMM-Newton/EPIC, an X-ray instrument commonly used in these studies. We devise and tune up two diagnostic tools, the first based on the mean-photon energy vs. count rate scatter plot and the second on the spectral analysis of the interaction region, that can be used to highlight the effects of thermal conduction and to derive the shock speed in case of efficient conduction at work. These tools can be used to ascertain information from X-ray observations, without the need to develop detailed and ad-hoc numerical models for the interpretation of the data.

HH objects are characterized by a complex knotty morphology detected mainly along the axis of protostellar jets in a wide range of bands. Evidence of interactions between knots formed in different epochs have been found, suggesting that jets may result from the ejection of plasma blobs from the source. We aim at investigating the physical mechanism leading to the irregular knotty structure observed in jets in different bands and the complex interactions occurring among blobs of plasma ejected from the stellar source. We perform 2D axisymmetric HD simulations of a randomly ejected pulsed jet. The jet consists of a train of blobs which ram with supersonic speed into the ambient medium. The initial random velocity of each blob follows an exponential distribution. We explore the ejection rate parameter to derive constraints on the physical properties of protostellar jets by comparison of model results with observations. Our model takes into account radiative losses and thermal conduction. We find that the mutual interactions of blobs ejected at different epochs and with different speed lead to a variety of plasma components not described by current models. The main features characterizing the random pulsed jet scenario are: single high speed knots, showing a measurable proper motion in nice agreement with observations; irregular chains of knots aligned along the jet axis and possibly interacting with each other; reverse shocks interacting with outgoing knots; oblique shocks produced by the reflection of shocks at the jet cocoon. All these structures concur to determine the morphology of the jet in different bands. We also find that the thermal conduction plays a crucial role in damping out HD instabilities that would develop within the cocoon and that contribute to the jet breaking.

We propose a method to synthesize the inverse Compton (IC) gamma-ray image of a supernova remnant starting from the radio (or hard X-ray) map and using results of the spatially resolved X-ray spectral analysis. The method is successfully applied to SN 1006. We found that synthesized IC gamma-ray images of SN 1006 show morphology in nice agreement with that reported by the H.E.S.S. collaboration. The good correlation found between the observed very-high energy gamma-ray and X-ray/radio appearance can be considered as an evidence that the gamma-ray emission of SN 1006 observed by H.E.S.S. is leptonic in origin, though the hadronic origin may not be excluded.

Efficient particle acceleration can modify the structure of supernova remnants. In this context we present the results of the combined analysis of the XMM-Newton EPIC archive observations of SN 1006. We aim at describing the spatial distribution of the physical and chemical properties of the X-ray emitting plasma at the shock front. We investigate the contribution of thermal and non-thermal emission to the X-ray spectrum at the rim of the remnant to study how the acceleration processes affect the X-ray emitting plasma. We perform a spatially resolved spectral analysis on a set of regions covering the whole rim of the shell and we exploit the results of the spectral analysis to produce a count-rate image of the "pure" thermal emission of SN 1006 in the 0.5-0.8 keV energy band (subtracting the non-thermal contribution). This image significantly differs from the total image in the same band, specially near the bright limbs. We find that thermal X-ray emission can be associated with the ejecta and we study the azimuthal variation of their physical and chemical properties finding anisotropies in temperature and chemical composition. Thanks to our thermal image we trace the position of the contact discontinuity over the whole shell and we compare it with that expected from 3-D MHD models of SNRs with unmodified shock. We conclude that the shock is modified everywhere in the rim and that the aspect angle between the interstellar magnetic field and the line of sight is significantly lower than 90 degrees

Recent developments on the study of mixed morphology supernova remnants (MMSNRs) have revealed the presence of metal rich X-ray emitting plasma inside a fraction of these remnant, a feature not properly addressed by traditional models for these objects. Radial profiles of thermodynamical and chemical parameters are needed for a fruitful comparison of data and model of MMSNRs, but these are available only in a few cases. We analyze XMM-Newton data of two MMSNRs, namely IC443 and G166.0+4.3, previously known to have solar metal abundances, and we perform spatially resolved spectral analysis of the X-ray emission. We detected enhanced abundances of Ne, Mg and Si in the hard X-ray bright peak in the north of IC443, and of S in the outer regions of G166.0+4.3. The metal abundances are not distributed uniformly in both remnants. The evaporating clouds model and the radiative SNR model fail to reproduce consistently all the observational results. We suggest that further deep X-ray observations of MMSNRs may reveal more metal rich objects. More detailed models which include ISM-ejecta mixing are needed to explain the nature of this growing subclass of MMSNRs.

A number of important processes taking place around strong shocks in supernova remnants (SNRs) depend on the shock obliquity. The measured synchrotron flux is a function of the aspect angle between interstellar magnetic field (ISMF) and the line of sight. Thus a model of non-thermal emission from SNRs should account for the orientation of the ambient magnetic field. We develop a new method for the estimation of the aspect angle, based on the comparison between observed and synthesized radio maps of SNRs, making different assumptions about the dependence of electron injection efficiency on the shock obliquity. The method uses the azimuthal profile of radio surface brightness as a probe for orientation of ambient magnetic field because it is almost insensitive to the downstream distribution of magnetic field and emitting electrons. We apply our method to a new radio image of SN 1006 produced on the basis of archival VLA and Parkes data. The image recovers emission from all spatial structures with angular scales from few arcseconds to 15 arcmin. We explore different models of injection efficiency and find the following best-fitting values for the aspect angle of SN 1006: phi=70 +/- 4.2 deg if the injection is isotropic, phi=64 +/- 2.8 deg for quasi-perpendicular injection (SNR has an equatorial belt in both cases) and phi=11 +/- 0.8 deg for quasi-parallel injection (polar-cap model of SNR). In the last case, SN 1006 is expected to have a centrally-peaked morphology contrary to what is observed. Therefore, our analysis provides some indication against the quasi-parallel injection model.

We investigate the spatial distribution of the physical and chemical properties of the hot X-ray emitting plasma of the supernova remnant IC 443, in order to get important constraints on its ionization stage, on the progenitor supernova explosion, on the age of the remnant, and its physical association with a close pulsar wind nebula. The hard X-ray thermal emission (1.4-5.0 keV) of IC 443 displays a centrally-peaked morphology, its brightness peaks being associated with hot (kT>1 keV) X-ray emitting plasma. A ring-shaped structure, characterized by high values of equivalent widths and median photon energy, encloses the PWN. Its hard X-ray emission is spectrally characterized by a collisional ionization equilibrium model, and strong emission lines of Mg, Si, and S, requiring oversolar metal abundances. Dynamically, the location of the ejecta ring suggests an SNR age of ~4,000 yr. The presence of overionized plasma in the inner regions of IC 443, addressed in previous works, is much less evident in our observations.

We present the results of the combined analysis of three XMM-Newton EPIC observations of the northern rim of the Vela SNR. The three pointings cover an area of ~10 pc^2 (at 250 pc) behind the main shock front and we aim at studying with high resolution the spatial distribution of the physical and chemical properties of the X-ray emitting plasma on this large scale. We produce count-rate images and equivalent width maps of the Ne IX and Mg XI emission blends. We also perform a spatially resolved spectral analysis of a set of physically homogeneous regions. We reveal physical and chemical inhomogeneities in the X-ray emitting plasma. In particular, we find large variations of the O, Ne, Mg, and Fe abundances. In some bright knots we also find unexpectedly enhanced Ne and Mg abundances, with values significantly larger than solar. Our results support a possible association of a few X-ray emitting knots with previously undetected residuals of stellar fragments (i. e. shrapnels) observed, in projection, inside the Vela shell.

In the framework of the study of supernova remnants and their complex interaction with the interstellar medium and the circumstellar material, we focus on the galactic supernova remnant W49B. Its morphology exhibits an X-ray bright elongated nebula, terminated on its eastern end by a sharp perpendicular structure aligned with the radio shell. The X-ray spectrum of W49B is characterized by strong K emission lines from Si, S, Ar, Ca, and Fe. There is a variation of the temperature in the remnant with the highest temperature found in the eastern side and the lowest one in the western side. The analysis of the recent observations of W49B indicates that the remnant may be the result of an asymmetric bipolar explosion where the ejecta are collimated along a jet-like structure and the eastern jet is hotter and more Fe-rich than the western one. Another possible scenario associates the X-ray emission with a spherical explosion where parts of the ejecta are interacting with a dense belt of ambient material. To overcome this ambiguity we present new results of the analysis of an XMM-Newton observation and we perform estimates of the mass and energy of the remnant. We conclude that the scenario of an anisotropic jet-like explosion explains quite naturally our observation results, but the association of W49B with a hypernova and a gamma-ray burst, although still possible, is not directly supported by any evidence.

In the framework of the study of the X-ray and optical emission in supernova remnants we focus on an isolated X-ray knot in the northern rim of the Vela SNR (Vela FilD), whose X-ray emission has been studied and discussed in Paper I. We aim at understanding the physical origin of the X-ray and optical emission in FilD, at understanding the role of the different physical processes at work, and at obtaining a key for the interpretation of future X-ray observations of SNRs. To this end we have pursued an accurate ``forward'' modeling of the interaction of the Vela SNR shock with an ISM cloud. We perform hydrodynamic simulations and we directly compare the observables synthesized from the simulations with the data. We explore four different model setups, choosing the values of the physical parameters on the basis of our preliminary analysis of the X-ray data. We synthesize X-ray emission maps and spectra filtered through the XMM-Newton EPIC-MOS instrumental response. The impact of a shock front at 6 million Kelvin on an elliptical cloud 30 times denser than the ambient medium describes well the shock-cloud interaction processes in the Vela FilD region in terms of spectral properties and morphology of the X-ray and optical emission. The bulk of the X-ray emission in the FilD knot originates in the cloud material heated by the transmitted shock front, but significant X-ray emission is also associated to the cloud material, which evaporates, as an effect of thermal conduction, in the intercloud medium. The physical origin of the FilD optical emission is associated to thermal instabilities. In the FilD knot the X-ray emission associated to the reflected shock front is negligible.

In the framework of the study of supernova remnants and their complex interaction with the interstellar medium, we report on an XMM-Newton EPIC observation of the Galactic supernova remnant W49B. We investigate the spatial distribution of the chemical and physical properties of the plasma, so as to get important constraints on the physical scenario, on the dynamics of the supernova explosion, and on the interaction of the supernova remnant with the ambient interstellar clouds. We present line images, equivalent width maps, and a spatially resolved spectral analysis of a set of homogeneous regions. The X-ray spectrum of W49B is characterized by strong K emission lines from Si, S, Ar, Ca and Fe. In all the regions studied, the X-ray spectrum is dominated by the emission from the ejecta and there is no indication of radial stratification of the elements. A high overabundance of Ni (Ni/Ni_sun=10^{+2}_{-1}) is required in the bright central region and the previous detection of Cr and Mn line emission is confirmed. Spectra are well described by two thermal components in collisional ionization equilibrium. We observe spatial variations in the temperature, with the highest temperature found in the east and the lowest in the west. Our results support a scenario in which the X-ray emission comes from ejecta interacting with a dense belt of ambient material, but another possibility is that the remnant is the result of an asymmetric bipolar explosion with the eastern jet being hotter and more Fe-rich than the western jet. The eastern jet is confined by interaction with ambient molecular clouds. Comparison of the observed abundances with yields for hypernova and supernova nucleosynthesis does not directly support the association of W49B with a gamma-ray burst, although it remains possible.

We analyzed an XMM-Newton EPIC observation of a bright knot, named FilD, in the northern rim of the Vela SNR, where the shock has encountered a cloud. The good combination of sensitivity, spectral, and spatial resolution allowed us to describe the internal structure of the observed ISM clouds and to obtain estimates of their temperature, density, O, Ne, and Fe abundances, and of their extension along the line of sight. We also examined the interaction of the shock with the FilD knot and estimated that the time elapsed from the shock impact is about one cloud crushing time. Our analysis allowed us to conclude that the observed X-ray emission is best explained by the propagation of transmitted shocks through ISM inhomogeneities with an inward increasing density profile. We found that the interaction with the shock determines the heating of the inner part of the clouds and the evaporation of the outer part.