Search

Monthly notices of the Royal Astronomical Society 505(1), 755 - 770 (2021). doi:10.1093/mnras/stab1319, Supernova remnants (SNRs) can be rich sources of information on the parent SN (supernova) explosion. Thus, investigating the transition from the phase of SN to that of SNR can be crucial to link these two phases of evolution. Here, we aim to study the early development of SNR in more detail, paying the major attention to the transition from the early expansion stage to the Sedov stage and the role played by magnetic field in this transition. To this end, spherical magnetohydrodynamic simulations of SNRs have been performed to study the evolution of magnetic field in young SNRs and explore a sequence of the SNR evolutionary stages in the pre-radiative epoch. Remnants of three supernova types are considered, namely, SNIa, SNIc, and SNIIP, that covers a wide space of parameters relevant for SNRs. Changes in global characteristics and development of spatial distributions are analysed. It is shown that the radial component of magnetic field rapidly drops downstream of the forward shock. Therefore, the radially aligned polarization patterns observed in few young SNRs cannot be reproduced in the 1D MHD simulations. The period SNR takes for the transition from the earliest ejecta-driven phase to the Sedov phase is long enough, with its distinctive physical features, headed by the energy conversion from mostly kinetic one to a fixed ratio between the thermal and kinetic components. This transition worth to be distinguished as a phase in SNR evolutionary scheme. The updated sequence of stages in SNR evolution could be the free expansion (of gas) – energy conversion – Sedov–Taylor – post-adiabatic – radiative., Published by Oxford Univ. Press, Oxford

We model the time-dependent radio emission from a disc accretion event in a T-Tauri star using 3D, ideal magnetohydrodynamic simulations combined with a gyrosynchrotron emission and radiative transfer model. We predict for the first time, the multifrequency (1–1000 GHz) intensity and circular polarization from a flaring T-Tauri star. A flux tube, connecting the star with its circumstellar disc, is populated with a distribution of non-thermal electrons that is allowed to decay exponentially after a heating event in the disc and the system is allowed to evolve. The energy distribution of the electrons, as well as the non-thermal power-law index and loss rate, are varied to see their effect on the overall flux. Spectra are generated from different lines of sight, giving different views of the flux tube and disc. The peak flux typically occurs around 20–30 GHz and the radio luminosity is consistent with that observed from T-Tauri stars. For all simulations, the peak flux is found to decrease and move to lower frequencies with elapsing time. The frequency-dependent circular polarization can reach 10$-30{{\ \rm per\ cent}}$ but has a complex structure that evolves as the flare evolves. Our models show that observations of the evolution of the spectrum and its polarization can provide important constraints on physical properties of the flaring environment and associated accretion event.

Since the day of its explosion, supernova (SN) 1987A has been closely monitored to study its evolution and to detect its central compact relic. In fact, the formation of a neutron star is strongly supported by the detection of neutrinos from the SN. However, besides the detection in the Atacama Large Millimeter/submillimeter Array (ALMA) data of a feature that is compatible with the emission arising from a proto-pulsar wind nebula (PWN), the only hint for the existence of such elusive compact object is provided by the detection of hard emission in NuSTAR data up to ~ 20 keV. We report on the simultaneous analysis of multi-epoch observations of SN 1987A performed with Chandra, XMM-Newton and NuSTAR. We also compare the observations with a state-of-the-art 3D magnetohydrodynamic (MHD) simulation of SN 1987A. A heavily absorbed power-law, consistent with the emission from a PWN embedded in the heart of SN 1987A, is needed to properly describe the high-energy part of the observed spectra. The spectral parameters of the best-fit power-law are in agreement with the previous estimate, and exclude diffusive shock acceleration as a possible mechanism responsible for the observed non-thermal emission. The information extracted from our analysis are used to infer the physical characteristics of the pulsar and the broad-band emission of its nebula, in agreement with the ALMA data. Analysis of the synthetic spectra also show that, in the near future, the main contribution to Fe K emission line will originate in the outermost shocked ejecta of SN 1987A., Accepted for publication in ApJ

This work was granted access to the HPC resources of MesoPSL financed by the Region Ile de France and the project EquipMeso (reference ANR-10-EQPX29-01) of the programme Investissements d’Avenir supervised by the Agence Nationale pour la Recherche.; International audience; Context. The equatorial accretion scenario, caused by the development of the Rayleigh-Taylor (RT) instability at the disk edge, was suggested by accurate three-dimensional magnetohydrodynamic (MHD) modelling, but no observational or experimental confirmation of such phenomena has been evidenced yet.Aims. We studied the propagation of a laterally extended laser-generated plasma stream across a magnetic field and investigated if this kind of structure can be scaled to the case of equatorial ‘tongue’ accretion channels in young stellar objects (YSOs); if so, this would support the possibility of equatorial accretion in young accreting stars.Methods. We conducted a scaled laboratory experiment at the PEARL laser facility. The experiment consists in an optical laser pulse that is focused onto the surface of a Teflon target. The irradiation of the target leads to the expansion of a hot plasma stream into the vacuum, perpendicularly to an externally applied magnetic field. We used a Mach-Zehnder interferometer to diagnose the plasma stream propagation along two axes, to obtain the three-dimensional distribution of the plasma stream.Results. The laboratory experiment shows the propagation of a laterally extended laser-generated plasma stream across a magnetic field. We demonstrate that: (i) such a stream is subject to the development of the RT instability, and (ii) the stream, decomposed into tongues, is able to efficiently propagate perpendicular to the magnetic field. Based on numerical simulations, we show that the origin of the development of the instability in the laboratory is similar to that observed in MHD models of equatorial tongue accretion in YSOs.Conclusions. As we verify that the laboratory plasma scales favourably to accretion inflows of YSOs, our laboratory results support the argument in favour of the possibility of the RT-instability-caused equatorial tongue accretion scenario in the astrophysical case.

Hot Jupiters (HJs) are massive gaseous planets orbiting close to their host stars. Due to their physical characteristics and proximity to the central star, HJs are the natural laboratories to study the process of star-planet interaction (SPI). Phenomena related to SPI may include the inflation and the evaporation of planetary atmospheres, the formation of cometary tails and bow shocks and magnetospheric interaction between the magnetic field of the planet and that of the star. Several works suggest that some systems show enhanced stellar activity in phase with the planetary rotation period. In this work, we use a 3D magneto-hydrodynamic model that describes a system composed of a star and an HJ and that includes the corresponding planetary and stellar winds. The aim is to investigate whether the material evaporating from the planet interacts with the stellar extended corona, and generates observable features. Our simulation shows that, in some conditions, the planetary wind expands and propagates mainly along the planetary orbit. Moreover, part of the planetary wind collides with the stellar wind and a fraction of the planet's outflow is funnelled by the stellar magnetic field and hits the stellar surface. In both events, the material is heated up to temperatures of a few MK by a shock. These phenomena could manifest in the form of enhanced stellar activity at some orbital phases of the planet., Accepted for publication in Astronomische Nachrichten - Astronomical Notes

Throughout an activity cycle, magnetic structures rise to the stellar surface, evolve and decay. Tracing their evolution on a stellar corona allows us to characterize the X-ray cycles. However, directly mapping magnetic structures is feasible only for the Sun, while such structures are spatially unresolved with present-day X-ray instruments on stellar coronae. We present here a new method, implemented by us, that indirectly reproduces the stellar X-ray spectrum and its variability with solar magnetic structures. The technique converts solar corona observations into a format identical to that of stellar X-ray observations and, specifically, XMM-Newton spectra. From matching these synthetic spectra with those observed for a star of interest, a fractional surface coverage with solar magnetic structures can be associated to each X-ray observation. We apply this method to two young solar-like stars: $\epsilon$ Eri ($\sim 400$ Myr), the youngest star to display a coronal cycle ($\sim 3$ yr), and Kepler 63 ($\sim 200$ Myr), for which the X-ray monitoring did not reveal a cyclic variability. We found that even during the cycle minimum a large portion of $\epsilon$ Eri's corona is covered with active structures. Therefore, there is little space for additional magnetic regions during the maximum, explaining the small observed cycle amplitude ($\Delta f \sim 0.12$) in terms of the X-ray luminosity. Kepler 63 displays an even higher coverage with magnetic structure than the corona of $\epsilon$ Eri. This supports the hypothesis that for stars younger than $, Comment: Submitted to AN for the Proceeding of the XMM-Newton 2021 Workshop. 7 pages, 4 figures

Aims. EXor-type objects are protostars that display powerful UV-optical outbursts caused by intermittent and powerful events of magnetospheric accretion. These objects are not yet well investigated and are quite difficult to characterize. Several parameters, such as plasma stream velocities, characteristic densities, and temperatures, can be retrieved from present observations. As of yet, however, there is no information about the magnetic field values and the exact underlying accretion scenario is also under discussion. Methods. We use laboratory plasmas, created by a high power laser impacting a solid target or by a plasma gun injector, and make these plasmas propagate perpendicularly to a strong external magnetic field. The propagating plasmas are found to be well scaled to the presently inferred parameters of EXor-type accretion event, thus allowing us to study the behaviour of such episodic accretion processes in scaled conditions. Results. We propose a scenario of additional matter accretion in the equatorial plane, which claims to explain the increased accretion rates of the EXor objects, supported by the experimental demonstration of effective plasma propagation across the magnetic field. In particular, our laboratory investigation allows us to determine that the field strength in the accretion stream of EXor objects, in a position intermediate between the truncation radius and the stellar surface, should be of the order of 100 G. This, in turn, suggests a field strength of a few kilogausses on the stellar surface, which is similar to values inferred from observations of classical T Tauri stars.

Charged particles can be accelerated to high energies by collisionless shock waves in astrophysical environments, such as supernova remnants. By interacting with the magnetized ambient medium, these shocks can transfer energy to particles. Despite increasing efforts in the characterization of these shocks from satellite measurements at Earth’s bow shock as well as powerful numerical simulations, the underlying acceleration mechanism or a combination thereof is still widely debated. Here we show that astrophysically relevant super-critical quasi-perpendicular magnetized collisionless shocks can be produced and characterized in the laboratory. We observe the characteristics of super-criticality in the shock profile as well as the energization of protons picked up from the ambient gas to hundreds of kiloelectronvolts. Kinetic simulations modelling the laboratory experiment identified shock surfing as the proton acceleration mechanism. Our observations not only provide direct evidence of early-stage ion energization by collisionless shocks but also highlight the role played by this particular mechanism in energizing ambient ions to feed further stages of acceleration. Furthermore, our results open the door to future laboratory experiments investigating the possible transition to other mechanisms, when increasing the magnetic field strength, or the effect that induced shock front ripples could have on acceleration processes. Proton acceleration by a super-critical collisionless shock is observed in laboratory experiments, and numerical simulations suggest shock surfing as the underlying acceleration mechanism.

The shaping of astrophysical outflows into bright, dense, and collimated jets due to magnetic pressure is here investigated using laboratory experiments. Here we look at the impact on jet collimation of a misalignment between the outflow, as it stems from the source, and the magnetic field. For small misalignments, a magnetic nozzle forms and redirects the outflow in a collimated jet. For growing misalignments, this nozzle becomes increasingly asymmetric, disrupting jet formation. Our results thus suggest outflow/magnetic field misalignment to be a plausible key process regulating jet collimation in a variety of objects from our Sun’s outflows to extragalatic jets. Furthermore, they provide a possible interpretation for the observed structuring of astrophysical jets. Jet modulation could be interpreted as the signature of changes over time in the outflow/ambient field angle, and the change in the direction of the jet could be the signature of changes in the direction of the ambient field., Mass outflow is a common process in astrophysical objects. Here the authors investigate in which conditions an astrophysically-scaled laser-produced plasma flow can be collimated and evolves in the presence of a misaligned external magnetic field.

Since the day of its explosion, SN 1987A (SN87A) was closely monitored with the aim to study its evolution and to detect its central compact relic. The detection of neutrinos from the supernova strongly supports the formation of a neutron star (NS). However, the constant and fruitless search for this object has led to different hypotheses on its nature. Up to date, the detection in the ALMA data of a feature somehow compatible with the emission arising from a proto Pulsar Wind Nebula (PWN) is the only hint of the existence of such elusive compact object. Here we tackle this 33-years old issue by analyzing archived observations of SN87A performed Chandra and NuSTAR in different years. We firmly detect nonthermal emission in the $10-20$ kev energy band, due to synchrotron radiation. The possible physical mechanism powering such emission is twofold: diffusive shock acceleration (DSA) or emission arising from an absorbed PWN. By relating a state-of-the-art magneto-hydrodynamic simulation of SN87A to the actual data, we reconstruct the absorption pattern of the PWN embedded in the remnant and surrounded by cold ejecta. We found that, even though the DSA scenario cannot be firmly excluded, the most likely scenario that well explains the data is the PWN emission., Accepted for publication in ApJ Letters

According to the most popular model for the origin of cosmic rays (CRs), supernova remnants (SNRs) are the site where CRs are accelerated. Observations across the electromagnetic spectrum support this picture through the detection of non-thermal emission that is compatible with being synchrotron or inverse Compton radiation from high energy electrons, or pion decay due to proton-proton interactions. These observations of growing quantity and quality promise to unveil many aspects of CRs acceleration and require more and more accurate tools for their interpretation. Here, we show how multi-dimensional MHD models of SNRs, including the effects on shock dynamics due to back-reaction of accelerated CRs and the synthesis of non-thermal emission, turned out to be very useful to investigate the signatures of CRs acceleration and to put constraints on the acceleration mechanism of high energy particles. These models have been used to interpret accurately observations of SNRs in various bands (radio, X-ray and $\gamma$-ray) and to extract from them key information about CRs acceleration., Comment: 16 pages, 5 figures, invited review accepted for publication on New Astronomy

Aims. Investigating the process of matter accretion onto forming stars through scaled experiments in the laboratory is important in order to better understand star and planetary system formation and evolution. Such experiments can indeed complement observations by providing access to the processes with spatial and temporal resolution. A previous investigation revealed the existence of a two-component stream: a hot shell surrounding a cooler inner stream. The shell was formed by matter laterally ejected upon impact and refocused by the local magnetic field. That laboratory investigation was limited to normal incidence impacts. However, in young stellar objects, the complex structure of magnetic fields causes variability of the incident angles of the accretion columns. This led us to undertake an investigation, using laboratory plasmas, of the consequence of having a slanted accretion impacting a young star. Methods. Here, we used high power laser interactions and strong magnetic field generation in the laboratory, complemented by numerical simulations, to study the asymmetry induced upon accretion structures when columns of matter impact the surface of young stars with an oblique angle. Results. Compared to the scenario where matter accretes perpendicularly to the star surface, we observe a strongly asymmetric plasma structure, strong lateral ejecta of matter, poor confinement of the accreted material, and reduced heating compared to the normal incidence case. Thus, slanted accretion is a configuration that seems to be capable of inducing perturbations of the chromosphere and hence possibly influencing the level of activity of the corona.

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., Comment: Accepted for publication in The Astrophysical Journal, 62 pages, 26 figures, 3 tables

(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., Comment: Accepted for publication in Astronomy & Astrophysics, 20 pages, 13 figures, 5 tables

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., 13 pages, 12 figures, accepted for publication in A&A

Abridged. We aim at exploring to which extent the remnant keeps memory of the asymmetries that develop stochastically in the neutrino-heating layer due to hydrodynamic instabilities (e.g., convective overturn and the standing accretion shock instability) during the first second after core bounce. We coupled a 3D HD model of a neutrino-driven SN explosion with 3D MHD/HD simulations of the remnant formation. The simulations cover 2000 years of expansion and include all physical processes relevant to describe the complexities in the SN evolution and the subsequent interaction of the stellar debris with the wind of the progenitor star. The interaction of large-scale asymmetries left from the earliest phases of the explosion with the reverse shock produces, at the age of $\approx 350$~years, an ejecta structure and a remnant morphology which are remarkably similar to those observed in Cas A. Small-scale structures in the large-scale Fe-rich plumes created during the initial stages of the SN, combined with HD instabilities that develop after the passage of the reverse shock, naturally produce a pattern of ring- and crown-like structures of shocked ejecta. The consequence is a spatial inversion of the ejecta layers with Si-rich ejecta being physically interior to Fe-rich ejecta. The full-fledged remnant shows voids and cavities in the innermost unshocked ejecta resulting from the expansion of Fe-rich plumes and their inflation due to the decay of radioactive species. The asymmetric distributions of $^{44}$Ti and $^{56}$Fe and their abundance ratio are both compatible with those inferred from high-energy observations of Chandra and NuSTAR. The main asymmetries observed in the ejecta distribution of Cas A can be explained by the interaction of the reverse shock with the large-scale asymmetries that developed from stochastic processes that originate during the first seconds of the SN blast., Comment: 33 pages, 29 figures, 1 Table, 4 movies, 3 interactive graphics, accepted for publication on A&A, high quality figures will be available on the published paper

The morphology and the distribution of material observed in SNRs reflect the interaction of the SN blast wave with the ambient environment, the physical processes associated with the SN explosion and the internal structure of the progenitor star. IC 443 is a MM SNR located in a quite complex environment: it interacts with a molecular cloud in the NW and SE areas and with an atomic cloud in the NE. In this work we aim at investigating the origin of the complex morphology and multi-thermal X-ray emission observed in SNR IC 443, through the study of the effect of the inhomogeneous ambient medium in shaping its observed structure, and the exploration of the main parameters characterizing the remnant. We developed a 3D HD model for IC 443, which describes the interaction of the SNR with the environment, parametrized in agreement with the results of the multi-wavelength data analysis. We performed an ample exploration of the parameter space describing the initial blast wave and the environment, and the surrounding clouds. From the simulations, we synthesized the X-ray emission maps and spectra and compared them with actual X-ray data collected by XMM-Newton. Our model explains the origin of the complex X-ray morphology of SNR IC 443 in a natural way, being able to reproduce, for the first time, most of the observed features, including the centrally-peaked X-ray morphology (characteristic of MM SNRs) when considering the origin of the explosion at the position where the PWN CXOU J061705.3+222127 was at the time of the explosion. In the model which best reproduces the observations, the mass of the ejecta and the energy of the explosion are $\sim 7 M_\odot$ and $\sim10^{51}$ erg, respectively. From the exploration of the parameter space, we found that the density of the clouds is $n>300$ cm$^{-3}$ and that the age of SNR IC 443 is $\sim8000$ yr., 16 pages, 12 figures. Accepted for publication in Astronomy and Astrophysics

Context. The manifold phases in the evolution of a core-collapse (CC) supernova (SN) play an important role in determining the physical properties and morphology of the resulting supernova remnant (SNR). Thus, the complex morphology of SNRs is expected to reflect possible asymmetries and structures developed during and soon after the SN explosion. Aims. The aim of this work is to bridge the gap between CC SNe and their remnants by investigating how post-explosion anisotropies in the ejecta influence the structure and chemical properties of the remnant at later times. Methods. We performed three-dimensional magneto-hydrodynamical simulations starting soon after the SN event and following the evolution of the system in the circumstellar medium, which includes the wind of the stellar progenitor, for 5000 yr, obtaining the physical scenario of a SNR. Here we focused the analysis on the case of a progenitor red supergiant of 19.8 M⊙. We also investigated how a post-explosion large-scale anisotropy in the SN affects the ejecta distribution and the matter mixing of heavy elements in the remnant during the first 5000 yr of evolution. Results. In the case of a spherically symmetric SN explosion without large-scale anisotropies, the remnant roughly keeps memory of the original onion-like layering of ejecta soon after the SN event. Nevertheless, as the reverse shock hits the ejecta, the element distribution departs from a homologous expansion because of the slowing down of the outermost ejecta layers due to interaction with the reverse shock. In the case of a large-scale anisotropy developed after the SN, we found that the chemical stratification in the ejecta can be strongly modified and the original onion-like layering is not preserved. The anisotropy may cause spatial inversion of ejecta layers, for instance leading to Fe/Si-rich ejecta outside the O shell, and may determine the formation of Fe/Si-rich jet-like features that may protrude the remnant outline. The level of matter mixing and the properties of the jet-like feature are sensitive to the initial physical (density and velocity) and geometrical (size and position) initial characteristics of the anisotropy.

LBVs are massive evolved stars that suffer sporadic and violent mass-loss events. They have been proposed as the progenitors of some core-collapse SNe, but this idea is still debated due to the lack of direct evidence. Since SNRs can carry in their morphology the fingerprints of the progenitor stars as well as of the inhomogeneous CSM sculpted by the progenitors, the study of SNRs from LBVs could help to place core-collapse SNe in context with the evolution of massive stars. We investigate the physical, chemical and morphological properties of the remnants of SNe originating from LBVs, in order to search for signatures, revealing the nature of the progenitors, in the ejecta distribution and morphology of the remnants. As a template of LBVs, we considered the actual LBV candidate Gal 026.47+0.02. We selected a grid of models, which describe the evolution of a massive star with properties consistent with those of Gal 026.47+0.02 and its final fate as core-collapse SN. We developed a 3D HD model that follows the post-explosion evolution of the ejecta from the breakout of the shock wave at the stellar surface to the interaction of the SNR with a CSM characterized by two dense nested toroidal shells, parametrized in agreement with multi-wavelength observations of Gal 026.47+0.02. Our models show a strong interaction of the blast wave with the CSM which determines an important slowdown of the expansion of the ejecta in the equatorial plane where the two shells lay, determining a high degree of asymmetry in the remnant. After 10000 years of evolution the ejecta show an elongated shape forming a broad jet-like structure caused by the interaction with the shells and oriented along the axis of the toroidal shells., 18 pages, 12 figures; accepted for publication in Astronomy & Astrophysics

Despite intense observational work, the stellar accretion process is insufficiently constrained, and consequently the details in the exchange of mass and momentum between the proto-star and its accretion disk remain approximate.

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 (e.g. HH 154 and DG Tau). We aim at investigating the origin of X-ray emission and the effect of perturbations in X-ray emitting stationary shocks in stellar jets. We performed a set of 2.5-dimensional MHD numerical simulations modelling supersonic pulsed jets ramming into a magnetized medium, exploring different parameters for the model. We consider two cases: a jet less dense than the ambient medium (HH 154), and a jet denser than the ambient (DG Tau). In both cases, we found that the jet is collimated by the magnetic field forming a quasi-stationary shock at the base which emits in X-rays even when perturbations are present.

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., Comment: The final version of this manuscript (including Methods) is published in Nature Astronomy and is publicly available at this link: https://rdcu.be/bhO7Z

CTTSs are young stars accreting mass from their circumstellar disks. The material falls into the star at free fall velocity and hits the stellar surface producing shocks, that heat the plasma at few million degrees. In the last twenty years the X-ray and UV observations of these systems have raised several questions. In particular, the observed X-ray luminosity is systematically below the value predicted by theoretical models, and the UV lines show complex profiles which cannot be easily interpreted with current accretion models based only on magnetohydrodynamical effects. To tackle these problems we modeled the structure and the dynamics of the plasma in the impact region using 3D magnetohydrodynamical simulations, but including the effects of radiative transport, for the first time in the Non Local Thermodynamical Equilibrium (NLTE) regime. We found that the radiation arising from the shocked plasma is partially absorbed by the unshocked accretion column, and this might explain the excess of X-ray flux predicted by MHD models which do not include absorption effects. Moreover, due to the absorption of radiation, the pre-shock downfalling accreted material is gradually heated up to a few 105 K due to irradiation of X-rays arising from the shocked plasma at the impact region. We discuss the implication of this pre-shock heating for the UV and X-ray emission arising from the impact region.

Models and observations indicate that the impact of matter accreting onto the surface of young stars produces regions at the base of accretion columns, in which optically thin and thick plasma components coexist. Thus an accurate description of these impacts requires to account for the effects of absorption and emission of radiation. We study the effects of radiation emerging from shock-heated plasma in impact regions on the structure of the pre-shock downfalling material. We investigate if a significant absorption of radiation occurs and if it leads to a pre-shock heating of the accreting gas. We developed a radiation hydrodynamics model describing an accretion column impacting onto the surface of a Classical T Tauri Star. The model takes into account the stellar gravity, the thermal conduction, and the effects of both radiative losses and absorption of radiation by matter in the non local thermodynamic equilibrium regime. After the impact, a hot slab of post-shock plasma develops at the base of the accretion column. Part of radiation emerging from the slab is absorbed by the pre-shock accreting material. As a result, the pre-shock accretion column gradually heats up to temperatures of $10^5$ K, forming a radiative precursor of the shock. The precursor has a thermal structure with the hottest part at $T \approx 10^5$ K, with size comparable to that of the hot slab, above the post-shock region. At larger distances the temperature gradually decreases to $T \approx 10^4$ K., Letter Accepted for publication in A&A

A new experimental set-up designed to investigate the accretion dynamics in newly born stars is presented. It takes advantage of a magnetically collimated stream produced by coupling a laser-generated expanding plasma to a $2\times 10^{5}~{G}\ (20~{T})$ externally applied magnetic field. The stream is used as the accretion column and is launched onto an obstacle target that mimics the stellar surface. This setup has been used to investigate in details the accretion dynamics, as reported in [G. Revet et al., Science Advances 3, e1700982 (2017), arXiv:1708.02528}. Here, the characteristics of the stream are detailed and a link between the experimental plasma expansion and a 1D adiabatic expansion model is presented. Dimensionless numbers are also calculated in order to characterize the experimental flow and its closeness to the ideal MHD regime. We build a bridge between our experimental plasma dynamics and the one taking place in the Classical T Tauri Stars (CTTSs), and we find that our set-up is representative of a high plasma $\beta$ CTTS accretion case.

Accretion of matter is a process that plays a central role in the physics of young stellar objects. The analysis of the structure by which matter settles on the star can unveil key information about the process of star formation by providing details on mass accretion rates, stellar magnetic field configurations, possible effects of accretion on the stellar coronal activity, etc. Here we review some of the achievements obtained by our group by exploiting a multi-disciplinary approach based on the analysis of multi-dimensional magnetohydrodynamic simulations, multi-wavelength observations, and laboratory experiments of accretion impacts occurring onto the surface of classical T Tauri stars (CTTSs).

Classical T Tauri stars (CTTSs) are young low-mass stellar objects accreting mass from their circumstellar disks. They are characterized by high levels of coronal activity as revealed by X-ray observations. This activity may affect the disk stability and the circumstellar environment. Here we investigate if an intense coronal activity due to flares occurring close to the accretion disk may perturb the inner disk stability, disrupt the inner part of the disk and, possibly, trigger accretion phenomena with rates comparable with those observed. We model a magnetized protostar surrounded by an accretion disk through 3D magnetohydrodinamic simulations. We explore cases characterized by a dipole plus an octupole stellar magnetic field configuration and different density of the disk or by different levels of flaring activity. As a result of the simulated intense flaring activity, we observe the formation of several loops that link the star to the disk; all these loops build up a hot extended corona with an X-ray luminosity comparable with typical values observed in CTTSs. The intense flaring activity close to the disk can strongly perturb the disk stability. The flares trigger overpressure waves which travel through the disk and modify its configuration. Accretion funnels may be triggered by the flaring activity, thus contributing to the mass accretion rate of the star. Accretion rates synthesized from the simulations are in a range between 10^-10 and 10^-9M_sun yr^-1 The accretion columns can be perturbed by the flares and they can interact with each other, possibly merging together in larger streams. As a result, the accretion pattern can be rather complex: the streams are highly inhomogeneous, with a complex density structure, and clumped, Comment: 12 pages, 9 figures

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., 15 pages, 9 Figures, accepted for publication on A&A, high quality figures will be available on the published paper

The Vela supernova remnant (SNR) shows several ejecta fragments protruding beyond the forward shock (shrapnel). Recent studies have revealed high Si abundance in two shrapnel (A and G), located in opposite directions with respect to the SNR center. This suggests the possible existence of a Si-rich jet-counterjet structure. We analyzed an XMM-Newton observation of a bright clump, behind shrapnel G, which lies along the direction connecting A and G. The aim is to study the physical and chemical properties of this clump to ascertain whether it is part of this putative jet-like structure. We produced background-corrected and adaptively-smoothed count-rate images and median photon energy maps, and performed a spatially resolved spectral analysis. We identified two structures with different physical properties. The first one is remarkably elongated along the direction connecting A and G. Its X-ray spectrum is much softer than that of the other two shrapnel, to the point of hindering the determination of the Si abundance, however its physical and chemical properties are consistent with those of shrapnel A and G. The second structure, running along the southeast-northwest direction, has a higher temperature and appears like a thin filament. By analyzing the ROSAT data, we have found that this filament is part of a very large and coherent structure that we identified in the western rim of the shell. We obtained a thorough description of the tail of Shrapnel G. In addition we discovered a coherent and very extended feature that we interpret as a signature of an earlier interaction of the remnant with the stellar wind of its progenitor star. The peculiar Ne/O ratio we found in the wind residual may be suggestive of a Wolf-Rayet progenitor for Vela SNR, though further analysis is required to address this point., 10 pages, 13 figures, Accepted to A&A

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., 16 pages, 17 Figures; accepted for publication on MNRAS

We present a three-dimensional kinematic reconstruction of the optically-emitting, oxygen-rich ejecta of supernova remnant N132D in the Large Magellanic Cloud. Data were obtained with the 6.5 m Magellan telescope in combination with the IMACS+GISMO instrument and survey [O III] $\lambda\lambda$4959,5007 line emission in a ${\sim}$3$^{\prime}~\times$ 3$^{\prime}$ region centered on N132D. The spatial and spectral resolution of our data enable detailed examination of the optical ejecta structure. The majority of N132D's optically bright oxygen ejecta are arranged in a torus-like geometry tilted approximately 28$^{\circ}$ with respect to the plane of the sky. The torus has a radius of 4.4 pc ($D_{\rm LMC}$/50 kpc), exhibits a blue-shifted radial velocity asymmetry of $-3000$ to $+2300$ km s$^{-1}$, and has a conspicuous break in its circumference. Assuming homologous expansion from the geometric center of O-rich filaments, the average expansion velocity of 1745 km s$^{-1}$ translates to an age since explosion of 2450 $\pm$ 195 yr. A faint, spatially-separated "runaway knot" (RK) with total space velocity of 3650 km s$^{-1}$ is nearly perpendicular to the torus plane and coincident with X-ray emission that is substantially enhanced in Si relative to the LMC and N132D's bulk ejecta. These kinematic and chemical signatures suggest that the RK may have had its origin deep within the progenitor star. Overall, the main shell morphology and high-velocity, Si-enriched components of N132D have remarkable similarity with that of Cassiopeia A, which was the result of a Type IIb supernova explosion. Our results underscore the need for further observations and simulations that can robustly reconcile whether the observed morphology is dominated by explosion dynamics or shaped by interaction with the environment., Comment: 18 pages, 12 figures, accepted for publication in ApJ; 3D animation of N132D can found at: http://www.physics.purdue.edu/kaboom/n132d-anim.mp4

High energy emission from young stars with disks, with all their components due to accretion and outflow activity, can have a deep impact on the evolution of their disks and on the formation of exo-planetary systems. An inter-disciplinary approach, which combines multi-wavelength observations, magnetohydrodynamical models, and laboratory experiments, allows us to get a more complete description of the accretion/ejection phenomena characterizing young stars. We discuss the case of the HH 154 jet, its X-ray emission localized at the base of the jet and its complex morphology, comparing observations, models, and laser experiments. We present the comparison between magnetohydrodynamical models prediction and high energy observations (UV and X-ray bands) of TW Hya, a promising object to perform also Doppler shift measurements, pushing to the limit the capabilities of currently available instruments (e.g. Chandra Telescope). We discuss how future missions as Athena will improve our understanding the accretion process in young stars.

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., Accepted for publication in A&A

We present a model for the radio evolution of supernova remnants (SNRs) obtained by using three-dimensional (3D) hydrodynamic simulations, coupled with nonlinear kinetic theory of cosmic ray (CR) acceleration in SNRs. We model the radio evolution of SNRs on a global level, by performing simulations for wide range of the relevant physical parameters, such as the ambient density, the supernova (SN) explosion energy, the acceleration efficiency and the magnetic field amplification (MFA) efficiency. We attribute the observed spread of radio surface brightnesses for corresponding SNR diameters to the spread of these parameters. In addition to our simulations of type Ia SNRs, we also considered SNR radio evolution in denser, nonuniform circumstellar environment, modified by the progenitor star wind. These simulations start with the mass of the ejecta substantially higher than in the case of a type Ia SN and presumably lower shock speed. The magnetic field is understandably seen as very important for the radio evolution of SNRs. In terms of MFA, we include both resonant and non-resonant modes in our large scale simulations, by implementing models obtained from first-principles, particle-in-cell (PIC) simulations and non-linear magnetohydrodynamical (MHD) simulations. We test the quality and reliability of our models on a sample consisting of Galactic and extragalactic SNRs. Our simulations give $\Sigma-D$ slopes between -4 and -6 for the full Sedov regime. Recent empirical slopes obtained for the Galactic samples are around -5, while for the extragalactic samples are around -4., Comment: 17 pages, 6 figures, 1 table, corrected typos, added new references and conclusions, 2 new coauthors, accepted for publication in ApJ

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..., Accepted for publication in Astronomy and Astrophysics

Monthly notices of the Royal Astronomical Society 479(3), 4253 - 4270 (2018). doi:10.1093/mnras/sty1750, We present very-high-resolution 1D MHD simulations of the late-stage supernova remnants (SNRs). 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., Published by Oxford Univ. Press, Oxford

According to the magnetospheric accretion scenario, young low-mass stars are surrounded by circumstellar disks which they interact with through accretion of mass. The accretion builds up the star to its final mass and is also believed to power the mass outflows, which may in turn have a significant role in removing the excess angular momentum from the star-disk system. Although the process of mass accretion is a critical aspect of star formation, some of its mechanisms are still to be fully understood. On the other hand, strong flaring activity is a common feature of young stellar objects (YSOs). In the Sun, such events give rise to perturbations of the interplanetary medium. Similar but more energetic phenomena occur in YSOs and may influence the circumstellar environment. In fact, a recent study has shown that an intense flaring activity close to the disk may strongly perturb the stability of circumstellar disks, thus inducing mass accretion episodes (Orlando et al. 2011). Here we review the main results obtained in the field and the future perspectives., 4 pages, 2 Figures; accepted for publication on Acta Polytechnica (Proceedings of the Frascati Workshop 2013)

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., accepted in MNRAS

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., Comment: 10 pages, 5 figures; submitted to: "SN 1987A, 30 years later", Proceedings IAU Symposium No. 331, 2017

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., Comment: A&A, accepted

Supernova remnant shocks are strong candidates for being the source of energetic cosmic rays and hadron acceleration is expected to increase the shock compression ratio, providing higher post-shock densities. We exploited the deep observations of the XMM-Newton Large Program on SN 1006 to verify this prediction. Spatially resolved spectral analysis led us to detect X-ray emission from the shocked ambient medium in SN 1006 and to find that its density significantly increases in regions where particle acceleration is efficient. Our results provide evidence for the effects of acceleration of cosmic ray hadrons on the post-shock plasma in supernova remnants.

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., Accepted for publication in A&A

Context. It is generally accepted that, in classical T Tauri stars, the plasma from the circumstellar disc accretes onto the stellar surface with free-fall velocity and the impact generates a shock. The impact region is expected to contribute to emission in different spectral bands; many studies have confirmed that the X-rays arise from the post-shock plasma but, otherwise, there are no studies in the literature investigating the origin of the observed UV emission which is apparently correlated to accretion. Aims: We investigated the effect of radiative heating of the infalling material by the post-shock plasma at the base of the accretion stream, with the aim to identify in which region a significant part of the UV emission originates. Methods: We developed a one-dimensional hydrodynamic model describing the impact of an accretion stream onto the stellar surface; the model takes into account the gravity, the radiative cooling of an optically thin plasma, the thermal conduction, and the heating due to absorption of X-ray radiation. The latter term represents the heating of the infalling plasma due to the absorption of X-rays emitted from the post-shock region. Results: We found that the radiative heating of the pre-shock plasma plays a non-negligible role in the accretion phenomenon. In particular, the dense and cold plasma of the pre-shock accretion column is gradually heated up to a few 105K due to irradiation of X-rays arising from the shocked plasma at the impact region. This heating mechanism does not affect significantly the dynamics of the post-shock plasma. On the other hand, a region of radiatively heated gas (that we consider a precursor) forms in the unshocked accretion column and contributes significantly to UV emission. Our model naturally reproduces the luminosity of UV emission lines correlated to accretion and shows that most of the UV emission originates from the precursor....

The recurrent symbiotic nova V745 Sco exploded on 2014 February 6 and was observed on February 22 and 23 by the Chandra X-ray Observatory Transmission Grating Spectrometers. By that time the supersoft source phase had already ended and Chandra spectra are consistent with emission from a hot, shock-heated circumstellar medium with temperatures exceeding 10^7K. X-ray line profiles are more sharply peaked than expected for a spherically-symmetric blast wave, with a full width at zero intensity of approximately 2400 km/s, a full width at half maximum of 1200 +/- 30 km/s and an average net blueshift of 165 +/- 10 km/s. The red wings of lines are increasingly absorbed toward longer wavelengths by material within the remnant. We conclude that the blast wave was sculpted by an aspherical circumstellar medium in which an equatorial density enhancement plays a role, as in earlier symbiotic nova explosions. Expansion of the dominant X-ray emitting material is aligned close to the plane of the sky and most consistent with an orbit seen close to face-on. Comparison of an analytical blast wave model with the X-ray spectra, Swift observations and near-infrared line widths indicates the explosion energy was approximately 10^43 erg, and confirms an ejected mass of approximately 10^-7 Msun. The total mass lost is an order of magnitude lower than the accreted mass required to have initiated the explosion, indicating the white dwarf is gaining mass and is a supernova Type 1a progenitor candidate., To appear in the Astrophysical Journal

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., Comment: Accepted for publication in A&A

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., Comment: Accepted for publication in Astronomy and Astrophysics

(Abridged) The symbiotic-like binary Mira and nova V407 Cyg was observed in outburst on March 2010 and monitored in several wavelength bands. Here we report on multi-dimensional hydrodynamic simulations describing the 2010 outburst of V407 Cyg, exploring the first 60 days of evolution. The model takes into account thermal conduction and radiative cooling; the pre-explosion system conditions included the companion star and a circumbinary density enhancement. The simulations showed that the blast and the ejecta distribution are both aspherical due to the inhomogeneous circumstellar medium in which they expand; in particular they are significantly collimated in polar directions (producing a bipolar shock morphology) if the circumstellar envelope is characterized by an equatorial density enhancement. The blast is partially shielded by the Mira companion, producing a wake with dense and hot post-shock plasma on the rear side of the companion star; most of the X-ray emission produced during the evolution of the blast arises from this plasma structure. The observed X-ray lightcurve can be reproduced, assuming values of outburst energy and ejected mass similar to those of RS Oph and U Sco, if a circumbinary gas density enhancement is included in the model. In this case, the 2010 blast propagated through a circumbinary gas density enhancement with radius of the order of 40 AU and gas density \approx 10^6 cm^{-3} and the mass of ejecta in the outburst was M_{ej} \approx 2\times 10^{-7} M_{\odot} with an explosion energy E_{0} \approx 2\times 10^{44} erg. Alternatively, the model can produce a similar X-ray lightcurve without the need of a circumbinary gas density enhancement only if the outburst energy and ejected mass were similar to those at the upper end of ranges for classical novae, namely M_{ej} \approx 5\times 10^{-5} M_{\odot} and E_{0} \approx 5\times 10^{46} erg.

Young low-mass stars are characterized by ejection of collimated outflows and by circumstellar disks which they interact with through accretion of mass. The accretion builds up the star to its final mass and is also believed to power the mass outflows, which may in turn remove the excess angular momentum from the star-disk system. However, although the process of mass accretion is a critical aspect of star formation, some of its mechanisms are still to be fully understood. A point not considered to date and relevant for the accretion process is the evidence of very energetic and frequent flaring events in these stars. Flares may easily perturb the stability of the disks, thus influencing the transport of mass and angular momentum. Here we report on three-dimensional magnetohydrodynamic modeling of the evolution of a flare with an idealized non–equilibrium initial condition occurring near the disk around a rotating magnetized star. The model takes into account the stellar magnetic field, the gravitational force, the viscosity of the disk, the magnetic-field-oriented thermal conduction (including the effects of heat flux saturation), the radiative losses from optically thin plasma, and the coronal heating. We show that, during its first stage of evolution, the flare gives rise to a hot magnetic loop linking the disk to the star. The disk is strongly perturbed by the flare: disk material evaporates under the effect of the thermal conduction and an overpressure wave propagates through the disk. When the overpressure reaches the opposite side of the disk, a funnel flow starts to develop there, accreting substantial disk material onto the young star from the side of the disk opposite to the flare.

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.

Distributions of nonthermal surface brightness of supernova remnants (SNRs) contain important information about the properties of magnetic field and acceleration of charged particles. In the present paper, the synchrotron radio, X-ray, and inverse-Compton (IC) gamma-ray maps of adiabatic SNRs in uniform interstellar medium and interstellar magnetic field are modeled and their morphology is analyzed, with particular emphasis to comparison of azimuthal and radial variations of brightness in radio, X-rays, and gamma-rays. Approximate analytical formulae for the azimuthal and radial profiles of the synchrotron radio and X-ray as well as the IC gamma-ray brightness are derived. They reveal the main factors which influence the pattern of the surface brightness distribution due to leptonic emission processes in shells of SNRs and can account for some non-linear effects of acceleration if necessary. These approximations provide observers and theorists with a set of simple diagnostic tools for quick analysis of the non-thermal maps of SNRs.