We model the time dependent radio emission from a disk 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 multi-frequency (1$-$1000 GHz) intensity and circular polarisation from a flaring T-Tauri star. A flux tube, connecting the star with its circumstellar disk, is populated with a distribution of non-thermal electrons which is allowed to decay exponentially after a heating event in the disk 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 disk. 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 polarisation can reach 10$-$30$\%$ but has a complex structure which evolves as the flare evolves. Our models show that observations of the evolution of the spectrum and its polarisation can provide important constraints on physical properties of the flaring environment and associated accretion event.

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.

Investigating in the laboratory the process of matter accretion onto forming stars through scaled experiments is important in order to better understand star and planetary systems formation and evolution. Such experiments can indeed complement observations by providing access to the processes with spatial and temporal resolution. A first step has been made in [G. Revet et al., Science Advances 3, e1700982 (2017), arXiv:1708.02528]. in allowing such investigations. It revealed the existence of a two components 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 tonormal incidence impacts. However, in young stellar objects, complex structure of magnetic fields causes variability of the incidentangles 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. Here we use 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. Compared to the scenario where matter accretes normal to the star surface, we observe 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 influence the level of activity of the corona.

The shaping of astrophysical outflows into bright, dense and collimated jets due to magnetic pressure is here investigated using laboratory experiments. We notably 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. 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.

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.

In 2015 we started the XMM-Newton monitoring of the young solar-like star Epsilon Eridani (440 Myr), one of the youngest solar-like stars with a known chromospheric CaII cycle. By analyzing the most recent Mount Wilson S-index CaII data of this star, we found that the chromospheric cycle lasts 2.92 +/- 0.02 yr, in agreement with past results. From the long-term X-ray lightcurve, we find clear and systematic X-ray variability of our target, consistent with the chromospheric CaII cycle. The average X-ray luminosity results to be 2 x 10^28 erg/s, with an amplitude that is only a factor 2 throughout the cycle. We apply a new method to describe the evolution of the coronal emission measure distribution of Epsilon Eridani in terms of solar magnetic structures: active regions, cores of active regions and flares covering the stellar surface at varying filling fractions. Combinations of these magnetic structures can describe the observed X-ray emission measure of Epsilon Eridani only if the solar flare emission measure distribution is restricted to events in the decay phase. The interpretation is that flares in the corona of Epsilon Eridani last longer than their solar counterparts. We ascribe this to the lower metallicity of Epsilon Eridani. Our analysis revealed also that the X-ray cycle of Epsilon Eridani is strongly dominated by cores of active regions. The coverage fraction of cores throughout the cycle changes by the same factor as the X-ray luminosity. The maxima of the cycle are characterized by a high percentage of covering fraction of the flares, consistent with the fact that flaring events are seen in the corresponding short-term X-ray lightcurves predominately at the cycle maxima. The high X-ray emission throughout the cycle of Epsilon Eridani is thus explained by the high percentage of magnetic structures on its surface.

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.

In this research note, we present 3DMAP-VR,(3-Dimensional Modeling of Astrophysical Phenomena in Virtual Reality), a project aimed at visualizing 3D MHD models of astrophysical simulations, using virtual reality sets of equipment. The models account for all the relevant physical processes in astrophysical phenomena: gravity, magnetic-field-oriented thermal conduction, energy losses due to radiation, gas viscosity, deviations from proton-electron temperature equilibration, deviations from the ionization equilibrium, cosmic rays acceleration, etc.. We realized an excellent synergy between our 3DMAP-VR project and Sketchfab (one of the largest open access platforms to publish and share 3D virtual reality and augmented reality content) to promote a wide dissemination of results for both scientific and public outreach purposes.

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.

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.

Modeling the dynamics of most astrophysical structures requires an adequate description of the radiation-matter interaction. Several numerical (magneto)hydrodynamics codes were upgraded with a radiation module to fulfill this request. However, those among them that use either the flux-limited diffusion (FLD) or the M1 radiation moment approaches are restricted to the local thermodynamic equilibrium (LTE). This assumption may be not valid in some astrophysical cases. We present an upgraded version of the LTE radiation-hydrodynamics module implemented in the PLUTO code, originally developed by Kolb et al. (2013), which we have extended to handle non-LTE regimes. Starting from the general frequency-integrated comoving-frame equations of radiation hydrodynamics (RHD), we have justified all the assumptions made to obtain the non-LTE equations actually implemented in the module, under the FLD approximation. An operator-split method is employed, with two substeps: the hydrodynamic part is solved with an explicit method by the solvers already available in PLUTO, the non-LTE radiation diffusion and energy exchange part is solved with an implicit method. The module is implemented in the PLUTO environment. It uses databases of radiative quantities that can be provided independently by the user: the radiative power loss, the Planck and Rosseland mean opacities. Our implementation has been validated through different tests, in particular radiative shock tests. The agreement with the semi-analytical solutions (when available) is good, with a maximum error of 7%. Moreover, we have proved that non-LTE approach is of paramount importance to properly model accretion shock structures. Our radiation FLD module represents a step toward the general non-LTE RHD modeling. The module is available, under request, for the community.

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.

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

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.

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

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.

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.

Accretion dynamics in the forming of young stars is still object of debate because of limitations in observations and modelling. Through scaled laboratory experiments of collimated plasma accretion onto a solid in the presence of a magnetic field, we open first window on this phenomenon by tracking, with spatial and temporal resolution, the dynamics of the system and simultaneously measuring multiband emissions. We observe in these experiments that matter, upon impact, is laterally ejected from the solid surface, then refocused by the magnetic field toward the incoming stream. Such ejected matter forms a plasma shell that envelops the shocked core, reducing escaped X-ray emission. This demonstrates one possible structure reconciling current discrepancies between mass accretion rates derived from X-ray and optical observations.

Aims. The data set of the long-term XMM-Newton monitoring program of HD 81809 is analyzed to study its X-ray cycle, to investigate if the latter is related to the chromospheric one, to infer the structure of the corona of HD 81809, and to explore if the coronal activity of HD 81809 can be ascribed to phenomena similar to the solar ones and, therefore, considered an extension of the solar case. Methods. We analyze the observations of HD 81809 performed with XMM-Newton with a regular cadence of 6 months from 2001 to 2016 and representing one of the longest available observational baseline ($\sim 15$~yr) for a solar-like star with a well-studied chromospheric cycle (with a period of $\sim 8$~yr). We investigate the modulation of coronal luminosity and temperature and its relation with the chromospheric cycle. We interpret the data in terms of a mixture of solar-like coronal regions, adopting a methodology originally proposed to study the Sun as an X-ray star. Results. The observations show a well-defined regular cyclic modulation of the X-ray luminosity that reflects the activity level of HD 81809. The data covers approximately two cycles of coronal activity; the modulation has an amplitude of a factor of $\sim 5$ (excluding evident flares, as in the June 2002 observation) and a period of $7.3\pm 1.5$~yr, consistent with that of the chromospheric cycle. We demonstrate that the corona of HD 81809 can be interpreted as an extension of the solar case and it can be modeled with a mixture of solar-like coronal regions along the whole cycle. The activity level is mainly determined by a varying coverage of very bright active regions, similar to cores of active regions observed in the Sun. Evidence of unresolved significant flaring activity is present especially in proximity of cycle maxima.

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.

IPCVs contain a magnetic, rotating white dwarf surrounded by a magnetically truncated accretion disk. To explain their strong flickering X-ray emission, accretion has been successfully taken into account. Nevertheless, observations suggest that accretion phenomena could not be the only process behind it. An intense flaring activity occurring on the surface of the disk may generate a corona, contribute to the thermal X-ray emission and influence the system stability. Our purposes are: investigating the formation of an extended corona above the accretion disk, due to an intense flaring activity occurring on the disk surface; studying its effects on the disk and stellar magnetosphere; assessing its contribution to the observed X-ray flux. We have developed a 3D MHD model of a IPCV. The model takes into account gravity, disk viscosity, thermal conduction, radiative losses and coronal flare heating. To perform a parameter space exploration, several system conditions have been considered, with different magnetic field intensity and disk density values. From the results of the evolution of the model, we have synthesized the thermal X-ray emission. The simulations show the formation of an extended corona, linking disk and star. The flaring activity is capable of strongly influencing the disk configuration and its stability, effectively deforming the magnetic field lines. Hot plasma evaporation phenomena occur in the layer immediately above the disk. The flaring activity gives rise to a thermal X-ray emission in both the [0.1-2.0] keV and the [2.0-10] keV bands. An intense coronal activity occurring on the disk surface of an IPCV can affect the structure of the disk depending noticeably on the density of the disk and the magnetic field of the central object. Moreover, the synthesis of the thermal X-ray fluxes shows that this flaring activity may contribute to the observed thermal X-ray emission.

Channelled fragmented downflows are ubiquitous in magnetized atmospheres, and have been recently addressed from an observation after a solar eruption. We study the possible back-effect of the magnetic field on the propagation of confined flows. We compare two 3D MHD simulations of dense supersonic plasma blobs downfalling along a coronal magnetic flux tube. In one, the blobs move strictly along the field lines; in the other, the initial velocity of the blobs is not perfectly aligned to the magnetic field and the field is weaker. The aligned blobs remain compact while flowing along the tube, with the generated shocks. The misaligned blobs are disrupted and merged by the chaotic shuffling of the field lines, and structured into thinner filaments; Alfven wave fronts are generated together with shocks ahead of the dense moving front. Downflowing plasma fragments can be chaotically and efficiently mixed if their motion is misaligned to field lines, with broad implications, e.g., disk accretion in protostars, coronal eruptions and rain.

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.

Dense plasma fragments were observed to fall back on the solar surface by the Solar Dynamics Observatory after an eruption on 7 June 2011, producing strong EUV brightenings. Previous studies investigated impacts in regions of weak magnetic field. Here we model the $\sim~300$ km/s impact of fragments channelled by the magnetic field close to active regions. In the observations, the magnetic channel brightens before the fragment impact. We use a 3D-MHD model of spherical blobs downfalling in a magnetized atmosphere. The blob parameters are constrained from the observation. We run numerical simulations with different ambient density and magnetic field intensity. We compare the model emission in the 171\AA~ channel of the Atmospheric Imaging Assembly with the observed one. We find that a model of downfall channelled in a $\sim~1$MK coronal loop confined by a magnetic field of $\sim~10-20$G, best explains qualitatively and quantitatively the observed evolution. The blobs are highly deformed, further fragmented, when the ram pressure becomes comparable to the local magnetic pressure and they are deviated to be channelled by the field, because of the differential stress applied by the perturbed magnetic field. Ahead of them, in the relatively dense coronal medium, shock fronts propagate, heat and brighten the channel between the cold falling plasma and the solar surface. This study shows a new mechanism which brightens downflows channelled by the magnetic field, such as in accreting young stars, and also works as a probe of the ambient atmosphere, providing information about the local plasma density and magnetic field.

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 1D 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 few $10^5$\;K 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.

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 perform MHD modeling of a single bright coronal loop to include the interaction with a non-uniform magnetic field. The field is stressed by random footpoint rotation in the central region and its energy is dissipated into heating by growing currents through anomalous magnetic diffusivity that switches on in the corona above a current density threshold. We model an entire single magnetic flux tube, in the solar atmosphere extending from the high-beta chromosphere to the low-beta corona through the steep transition region. The magnetic field expands from the chromosphere to the corona. The maximum resolution is ~30 km. We obtain an overall evolution typical of loop models and realistic loop emission in the EUV and X-ray bands. The plasma confined in the flux tube is heated to active region temperatures (~3 MK) after ~2/3 hr. Upflows from the chromosphere up to ~100 km/s fill the core of the flux tube to densities above 10^9 cm^-3. More heating is released in the low corona than the high corona and is finely structured both in space and time.

Context. The accretion process in Classical T Tauri Stars (CTTSs) can be studied through the analysis of some UV and X-ray emission lines which trace hot gas flows and act as diagnostics of the post-shock downfalling plasma. In the UV band, where higher spectral resolution is available, these lines are characterized by rather complex profiles whose origin is still not clear. Aims. We investigate the origin of UV and X-ray emission at impact regions of density structured (fragmented) accretion streams.We study if and how the stream fragmentation and the resulting structure of the post-shock region determine the observed profiles of UV and X-ray emission lines. Methods. We model the impact of an accretion stream consisting of a series of dense blobs onto the chromosphere of a CTTS through 2D MHD simulations. We explore different levels of stream fragmentation and accretion rates. From the model results, we synthesize C IV (1550 {\AA}) and OVIII (18.97 {\AA}) line profiles. Results. The impacts of accreting blobs onto the stellar chromosphere produce reverse shocks propagating through the blobs and shocked upflows. These upflows, in turn, hit and shock the subsequent downfalling fragments. As a result, several plasma components differing for the downfalling velocity, density, and temperature are present altoghether. The profiles of C IV doublet are characterized by two main components: one narrow and redshifted to speed $\approx$ 50 km s$^{-1}$ and the other broader and consisting of subcomponents with redshift to speed in the range 200 $\approx$ 400 km s$^{-1}$. The profiles of OVIII lines appear more symmetric than C IV and are redshifted to speed $\approx$ 150 km s$^{-1}$. Conclusions. Our model predicts profiles of C IV line remarkably similar to those observed and explains their origin in a natural way as due to stream fragmentation.

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.

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.

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.

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.

A solar eruption after a flare on 7 Jun 2011 produced EUV-bright impacts of fallbacks far from the eruption site, observed with the Solar Dynamics Observatory. These impacts can be taken as a template for the impact of stellar accretion flows. Broad red-shifted UV lines have been commonly observed in young accreting stars. Here we study the emission from the impacts in the Atmospheric Imaging Assembly's UV channels and compare the inferred velocity distribution to stellar observations. We model the impacts with 2D hydrodynamic simulations. We find that the localised UV 1600A emission and its timing with respect to the EUV emission can be explained by the impact of a cloud of fragments. The first impacts produce strong initial upflows. The following fragments are hit and shocked by these upflows. The UV emission comes mostly from the shocked front shell of the fragments while they are still falling, and is therefore redshifted when observed from above. The EUV emission instead continues from the hot surface layer that is fed by the impacts. Fragmented accretion can therefore explain broad redshifted UV lines (e.g. C IV 1550A) to speeds around 400 km/s observed in accreting young stellar objects.

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.

Observations reveal a correspondence between chromospheric type II spicules and bright upwardly moving fronts in the corona observed in the EUV band. However, theoretical considerations suggest that these flows are unlikely to be the main source of heating in coronal magnetic loops. We investigate the propagation of high-speed chromospheric flows into coronal magnetic flux tubes, and the possible production of emission in the EUV band. We simulate the propagation of a dense $10^4$ K chromospheric jet upwards along a coronal loop, by means of a 2-D cylindrical MHD model, including gravity, radiative losses, thermal conduction and magnetic induction. The jet propagates in a complete atmosphere including the chromosphere and a tenuous cool ($\sim 0.8$ MK) corona, linked through a steep transition region. In our reference model, the jet's initial speed is 70 km/s, its initial density is $10^{11}$ cm$^{-3}$, and the ambient uniform magnetic field is 10 G. We explore also other values of jet speed and density in 1-D, and of magnetic field in 2-D, and the jet propagation in a hotter ($\sim 1.5$ MK) background loop. While the initial speed of the jet does not allow it to reach the loop apex, a hot shock front develops ahead of it and travels to the other extreme of the loop. The shock front compresses the coronal plasma and heats it to about $10^6$ K. As a result, a bright moving front becomes visible in the 171 \AA\ channel of the SDO/AIA mission. This result generally applies to all the other explored cases, except for the propagation in the hotter loop. For a cool, low-density initial coronal loop, the post-shock plasma ahead of upward chromospheric flows might explain at least part of the observed correspondence between type II spicules and EUV emission excess.

The expansion of coronal loops in the transition region may considerably influence the diagnostics of the plasma emission measure. The cross sectional area of the loops is expected to depend on the temperature and pressure, and might be sensitive to the heating rate. The approach here is to study the area response to slow changes in the coronal heating rate, and check the current interpretation in terms of steady heating models. We study the area response with a time-dependent 2D MHD loop model, including the description of the expanding magnetic field, coronal heating and losses by thermal conduction and radiation from optically thin plasma. We run a simulation for a loop 50 Mm long and quasi-statically heated to about 4 MK. We find that the area can change substantially with the quasi-steady heating rate, e.g. by ~40% at 0.5 MK as the loop temperature varies between 1 and 4 MK, and, therefore, affects the interpretation of DEM(T) curves.

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.

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.

(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.

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.

Context. Theoretical arguments and numerical simulations of radiative shocks produced by the impact of the accreting gas onto young stars predict quasi-periodic oscillations in the emitted radiation. However, observational data do not show evidence of such periodicity. Aims. We investigate whether physically plausible perturbations in the accretion column or in the chromosphere could disrupt the shock structure influencing the observability of the oscillatory behavior. Methods. We performed local 2D magneto-hydrodynamical simulations of an accretion shock impacting a chromosphere, taking optically thin radiation losses and thermal conduction into account. We investigated the effects of several perturbation types, such as clumps in the accretion stream or chromospheric fluctuations, and also explored a wide range of plasma-\beta values. Results. In the case of a weak magnetic field, the post-shock region shows chaotic motion and mixing, smoothing out the perturbations and retaining a global periodic signature. On the other hand, a strong magnetic field confines the plasma in flux tubes, which leads to the formation of fibrils that oscillate independently. Realistic values for the amplitude, length, and time scales of the perturbation are capable of bringing the fibril oscillations out of phase, suppressing the periodicity of the emission. Conclusions. The strength of a locally uniform magnetic field in YSO accretion shocks determines the structure of the post-shock region, namely, whether it will be somewhat homogeneous or if it will split up to form a collection of fibrils. In the second case, the size and shape of the fibrils is found to depend strongly on the plasma-\beta value but not on the perturbation type. Therefore, the actual value of the protostellar magnetic field is expected to play a critical role in the time dependence of the observable emission.

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.

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.

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.

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.

In the present work we study the C8 flare occurred on September 26, 2000 at 19:49 UT and observed by the SOHO/SUMER spectrometer from the beginning of the impulsive phase to well beyond the disappearance in the X-rays. The emission first decayed progressively through equilibrium states until the plasma reached 2-3 MK. Then, a series of cooler lines, i.e. Ca x, Ca vii, Ne vi, O iv and Si iii (formed in the temperature range log T = 4.3 - 6.3 under equilibrium conditions), are emitted at the same time and all evolve in a similar way. Here we show that the simultaneous emission of lines with such a different formation temperature is due to thermal instability occurring in the flaring plasma as soon as it has cooled below ~ 2 MK. We can qualitatively reproduce the relative start time of the light curves of each line in the correct order with a simple (and standard) model of a single flaring loop. The agreement with the observed light curves is greatly improved, and a slower evolution of the line emission is predicted, if we assume that the model loop consists of an ensemble of subloops or strands heated at slightly different times. Our analysis can be useful for flare observations with SDO/EVE.

(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.

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.

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.

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.

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.

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.

High resolution X-ray spectroscopy has revealed soft X-rays from high density plasma in Classical T-Tauri stars (CTTSs), probably arising from the accretion shock region. However, the mass accretion rates derived from the X-ray observations are consistently lower than those derived from UV/optical/NIR studies. We aim to test the hypothesis that the high density soft X-ray emission is from accretion by analysing optical accretion tracers from an X-ray selected sample of CTTSs in a homogeneous manner. We analyse optical spectra of a sample of CTTSs and calculate the accretion rates based on measuring optical emission lines. These are then compared to the accretion rates derived from the X-ray spectroscopy. We find that, for each CTTS in our sample, the different optical tracers predict mass accretion rates that agree within the errors, albeit with a spread of ~1 order of magnitude. Typically, mass accretion rates derived from Halpha and HeI 5876 Ang are larger than those derived from Hbeta, Hgamma and OI. When comparisons of the optical mass accretion rates are made to the X-ray derived mass accretion rates, we find that: a) the latter are always lower (but by varying amounts); b) the latter range within a factor of ~2 around 2x10^{-10} M_odot yr^{-1}, despite the fact that the former span a range of ~3 orders of magnitude. We suggest that the systematic underestimation of the X-ray derived mass accretion rates could depend on the density distribution inside the accretion streams, where the densest part of the stream is not visible in the X-ray band because of the absorption by the stellar atmosphere. We also suggest that a non-negligible optical depth of X-ray emission lines produced by post-shock accreting plasma may explain the almost constant mass accretion rates derived in X-rays if the effect is larger in stars with larger optical mass accretion rates.

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.

Three-dimensional hydrodynamic simulations exploring the first 18 hours of the 2010 January 28 outburst of the recurrent nova U Scorpii have been performed. Special emphasis was placed on capturing the enormous range in spatial scales in the blast. The pre-explosion system conditions included the secondary star and a flared accretion disk. These conditions can have a profound influence on the evolving blast wave. The blast itself is shadowed by the secondary star, which itself gives rise to a low-temperature bow-shock. The accretion disk is completely destroyed in the explosion. A model with a disk gas density of 10^{15} cm^{-3} produced a blast wave that is collimated and with clear bipolar structures, including a bipolar X-ray emitting shell. The degree of collimation depends on the initial mass of ejecta, energy of explosion, and circumstellar gas density distribution. It is most pronounced for a model with the lowest explosion energy (10^{43} erg) and mass of ejecta (10^{-8} M_{\odot}). The X-ray luminosities of three of six models computed are close to, but consistent with, an upper limit to the early blast X-ray emission obtained by the Swift satellite, the X-ray luminosity being larger for higher circumstellar gas density and higher ejecta mass. The latter consideration, together with estimates of the blast energy from previous outbursts, suggests that the mass of ejecta in the 2010 outburst was not larger than 10^{-7} M_{\odot}.

Context. High resolution X-ray observations of classical T Tauri stars (CTTSs) show a soft X-ray excess due to high density plasma (n_e=10^11-10^13 cm^-3). This emission has been attributed to shock-heated accreting material impacting onto the stellar surface. Aims. We investigate the observability of the shock-heated accreting material in the X-ray band as a function of the accretion stream properties (velocity, density, and metal abundance) in the case of plasma-beta<<1 in the post-shock zone. Methods. We use a 1-D hydrodynamic model describing the impact of an accretion stream onto the chromosphere, including the effects of radiative cooling, gravity and thermal conduction. We explore the space of relevant parameters and synthesize from the model results the X-ray emission in the [0.5-8.0] keV band and in the resonance lines of O VII (21.60 Ang) and Ne IX (13.45 Ang), taking into account the absorption from the chromosphere. Results. The accretion stream properties influence the temperature and the stand-off height of the shocked slab and its sinking in the chromosphere, determining the observability of the shocked plasma. Our model predicts that X-ray observations preferentially detect emission from low density and high velocity shocked accretion streams due to the large absorption of dense post-shock plasma. In all the cases examined, the post-shock zone exhibits quasi-periodic oscillations due to thermal instabilities, but in the case of inhomogeneous streams and beta<<1, the shock oscillations are hardly detectable. Conclusions. We suggest that, if accretion streams are inhomogeneous, the selection effect introduced by the absorption on observable plasma components may explain the discrepancy between the accretion rate measured by optical and X-ray data as well as the different densities measured using different He-like triplets in the X-ray 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.

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.

AIMS. We investigate the stability and dynamics of accretion shocks in CTTSs, considering the case of beta >= 1 in the post-shock region. In these cases the 1D approximation is not valid and a multi-dimensional MHD approach is necessary. METHODS. We model an accretion stream propagating through the atmosphere of a CTTS and impacting onto its chromosphere, by performing 2D axisymmetric MHD simulations. The model takes into account the stellar magnetic field, the gravity, the radiative cooling, and the thermal conduction (including the effects of heat flux saturation). RESULTS. The dynamics and stability of the accretion shock strongly depends on the plasma beta. In the case of shocks with beta > 10, violent outflows of shock-heated material (and possibly MHD waves) are generated at the base of the accretion column and strongly perturb the surrounding stellar atmosphere and the accretion column itself (modifying, therefore, the dynamics of the shock). In shocks with beta ~ 1, the post-shock region is efficiently confined by the magnetic field. The shock oscillations induced by cooling instability are strongly influenced by beta: for beta > 10, the oscillations may be rapidly dumped by the magnetic field, approaching a quasi-stationary state, or may be chaotic with no obvious periodicity due to perturbation of the stream induced by the post-shock plasma itself; for beta ~ 1 the oscillations are quasi-periodic, although their amplitude is smaller and the frequency higher than those predicted by 1D models.

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.

Accretion shocks have been recognized as important X-ray emission mechanism for pre-main sequence stars. Yet the X-ray properties of FUor outbursts, events that are caused by violent accretion, have been given little attention. We have observed the FUor object Z CMa during optical outburst and quiescence with Chandra. No significant changes in X-ray brightness and spectral shape are found, suggesting that the X-ray emission is of coronal nature. Due to the binary nature of Z CMa the origin of the X-ray source is ambiguous. However, the moderate hydrogen column density derived from our data makes it unlikely that the embedded primary star is the X-ray source. The secondary star, which is the FUor object, is thus responsible for both the X-ray emission and the presently ongoing accretion outburst, which seem however to be unrelated phenomena. The secondary is also known to drive a large outflow and jet, that we detect here for the first time in X-rays. The distance of the X-ray emitting outflow source to the central star is higher than in jets of low-mass stars.

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

The synchrotron radio maps of supernova remnants (SNRs) in uniform interstellar medium and interstellar magnetic field (ISMF) are analyzed, allowing different `sensitivity' of injection efficiency to the shock obliquity. The very-high energy gamma-ray maps due to inverse Compton process are also synthesized. The properties of images in these different wavelength bands are compared, with particular emphasis on the location of the bright limbs in bilateral SNRs. Recent H.E.S.S. observations of SN 1006 show that the radio and IC gamma-ray limbs coincide, and we found that this may happen if: i) injection is isotropic but the variation of the maximum energy of electrons is rather quick to compensate for differences in magnetic field; ii) obliquity dependence of injection (either quasi-parallel or quasi-perpendicular) and the electron maximum energy is strong enough to dominate magnetic field variation. In the latter case, the obliquity dependence of the injection and the maximum energy should not be opposite. We argue that the position of the limbs alone and even their coincidence in radio, X-rays and gamma-rays, as it is discovered by H.E.S.S. in SN 1006, cannot be conclusive about the dependence of the electron injection efficiency, the compression/amplification of ISMF and the electron maximum energy on the obliquity angle.

Chandra/HETG observations of the recurrent nova RS Ophiuchi at day 13.9 of its 2006 outburst reveal a spectrum covering a large range in plasma temperature and characterized by asymmetric and blue-shifted emission lines. We investigate the origin of these asymmetries and broadening of emission lines. We perform 3-D hydrodynamic simulations of the blast wave from the 2006 outburst, propagating through the inhomogeneous CSM. The model takes into account the thermal conduction (including the effects of heat flux saturation) and the radiative cooling. From the simulations, we synthesize the X-ray emission and derive the spectra as they would be observed with Chandra/HETG. Our model reproduces the observed X-ray emission in a natural way if the CSM in which the outburst occurred is characterized by an equatorial density enhancement. Our ``best-fit'' model predicts that most of the early X-ray emission originates from a small region propagating in the direction perpendicular to the line-of-sight and localized just behind the interaction front between the blast wave and the equatorial density enhancement. The model predicts asymmetric and blue-shifted line profiles remarkably similar to those observed. These asymmetries are due to substantial X-ray absorption of red-shifted emission by ejecta material. The comparison of high quality data of Chandra/HETG with detailed hydrodynamic modeling has allowed us to unveil, for the first time, the details of the structure emitting in the X-ray band in early phases of the outburst evolution, contributing to a better understanding of the physics of interactions between nova blasts and CSM in recurrent novae. This may have implications for whether or not RS Ophiuchi is a Type Ia SN progenitor system.

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.

High spectral resolution X-ray observations of CTTSs demonstrate the presence of plasma at T~2-3X10^6 K and n_e~10^11-10^13 cm^-3, unobserved in non-accreting stars. Stationary models suggest that this emission is due to shock-heated accreting material, but they do not allow to analyze the stability of such material and its position in the stellar atmosphere. We investigate the dynamics and the stability of shock-heated accreting material in CTTSs and the role of the stellar chromosphere in determining the position and the thickness of the shocked region. We perform 1-D HD simulations of the impact of the accretion flow onto chromosphere of a CTTS, including the effects of gravity, radiative losses from optically thin plasma, thermal conduction and a well tested detailed model of the stellar chromosphere. Here we present the results of a simulation based on the parameters of the CTTS MP Mus. We find that the accretion shock generates an hot slab of material above the chromosphere with a maximum thickness of 1.8X10^9 cm, density n_e~10^11-10^2 cm^-3, temperature T~3X10^6 K and uniform pressure equal to the ram pressure of the accretion flow (~450 dyn cm^-2). The base of the shocked region penetrates the chromosphere and stays where the ram pressure is equal to the thermal pressure. The system evolves with quasi-periodic instabilities of the material in the slab leading to cyclic disappearance and re-formation of the slab. For an accretion rate of ~10^-10 M_sun yr^-1, the shocked region emits a time-averaged X-ray luminosity L_X~7X10^29 erg s^-1, which is comparable to the X-ray luminosity observed in CTTSs of the same mass. Furthermore, the X-ray spectrum synthesized from the simulation matches in detail all the main features of the O VIII and O VII lines of the star MP Mus.

(abridged) Our long-term XMM-Newton program of long-term monitoring of a solar-like star with a well-studied chromospheric cycle, HD 81809 aims to study whether an X-ray cycle is present, along with studying its characteristics and its relation to the chromospheric cycle. Regular observations of HD 81809 were performed with XMM-Newton, spaced by 6 months from 2001 to 2007. We studied the variations in the resulting coronal luminosity and temperature, and compared them with the chromospheric CaII variations. We also modeled the observations in terms of a mixture of active regions, using a methodology originally developed to study the solar corona. Our observations show a well-defined cycle with an amplitude exceeding 1 dex and an average luminosity approximately one order of magnitude higher than in the Sun. The behavior of the corona of HD 81809 can be modeled well in terms of varying coverage of solar-like active regions, with a larger coverage than for the Sun, showing it to be compatible with a simple extension of the solar case.

Impulsive nanoflares are expected to transiently heat the plasma confined in coronal loops to temperatures of the order of 10 MK. Such hot plasma is hardly detected in quiet and active regions, outside flares. During rapid and short heat pulses in rarified loops the plasma can be highly out of equilibrium of ionization. Here we investigate the effects of the non-equilibrium of ionization (NEI) on the detection of hot plasma in coronal loops. Time-dependent loop hydrodynamic simulations are specifically devoted to this task, including saturated thermal conduction, and coupled to the detailed solution of the equations of ionization rate for several abundant elements. In our simulations, initially cool and rarified magnetic flux tubes are heated to 10 MK by nanoflares deposited either at the footpoints or at the loop apex. We test for different pulse durations, and find that, due to NEI effects, the loop plasma may never be detected at temperatures above ~5 MK for heat pulses shorter than about 1 min. We discuss some implications in the framework of multi-stranded nanoflare-heated coronal loops.

The solar corona is a template to understand stellar activity. The Sun is a moderately active star, and its corona differs from active stars: active stellar coronae have a double-peaked EM(T) with the hot peak at 8-20 MK, while the non flaring solar corona has one peak at 1-2 MK. We study the average contribution of flares to the solar EM(T) to investigate indirectly the hypothesis that the hot peak of the EM(T) of active stellar coronae is due to a large number of unresolved solar-like flares, and to infer properties on the flare distribution from nano- to macro-flares. We measure the disk-integrated time-averaged emission measure, EM_F(T), of an unbiased sample of solar flares analyzing uninterrupted GOES/XRS light curves over time intervals of one month. We obtain the EM_Q(T) of quiescent corona for the same time intervals from the Yohkoh/SXT data. To investigate how EM_F(T) and EM_Q(T) vary with the solar cycle, we evaluate them at different phases of the cycle (from Dec. 1991 to Apr. 1998). Irrespective of the solar cycle phase, EM_F(T) appears like a peak of the distribution significantly larger than the values of EM_Q(T) for T~5-10 MK. As a result the time-averaged EM(T) of the whole solar corona is double-peaked, with the hot peak, due to time-averaged flares, located at temperature similar of that of active stars, but less enhanced. The EM_F(T) shape supports the hypothesis that the hot EM(T) peak of active coronae is due to unresolved solar-like flares. If this is the case, quiescent and flare components should follow different scaling laws for increasing stellar activity. In the assumption that the heating of the corona is entirely due to flares, from nano- to macro-flares, then either the flare distribution or the confined plasma response to flares, or both, are bimodal.

We provide a guideline to interpret the UVCS emission lines (in particular O VI and Si XII) during shock wave propagation in the outer solar corona. We use a numerical MHD model performing a set of simulations of shock waves generated in the corona and from the result we compute the plasma emission for the O VI and Si XII including the effects of NEI. We analyze the radiative and spectral properties of our model with the support of a detailed radiation model including Doppler dimming and an analytical model for shocks, and, finally, we synthesize the expected O VI 1032A line profile. We explain several spectral features of the observations like the absence of discontinuities in the O VI emission during the shock passage, the brightening of Si XII emission and the width of the lines. We use our model also to give very simple and general predictions for the strength of the line wings due to the ions shock heating and on the line shape for Limb CMEs or Halo CMEs. The emission coming from post-shock region in the solar corona roughly agrees with the emission from a simple planar and adiabatic shock, but the effect of thermal conduction and the magnetic field may be important depending on the event parameters. Doppler dimming significantly influences the O VI emission while Si XII line brightens mainly because of the shock compression. Significant shock heating is responsible for the wide and faint component of the O VI line usually observed which may be taken as a shock signature in the solar corona.

We explore the importance of magnetic-field-oriented thermal conduction in the interaction of supernova remnant (SNR) shocks with radiative gas clouds and in determining the mass and energy exchange between the clouds and the hot surrounding medium. We perform 2.5D MHD simulations of a shock impacting on an isolated gas cloud, including anisotropic thermal conduction and radiative cooling; we consider the representative case of a Mach 50 shock impacting on a cloud ten-fold denser than the ambient medium. We consider different configurations of the ambient magnetic field and compare MHD models with or without the thermal conduction. The efficiency of the thermal conduction in the presence of magnetic field is, in general, reduced with respect to the unmagnetized case. The reduction factor strongly depends on the initial magnetic field orientation, and it is minimum when the magnetic field is initially aligned with the direction of shock propagation. The thermal conduction contributes to suppress hydrodynamic instabilities, reducing the mass mixing of the cloud and preserving the cloud from complete fragmentation. Depending on the magnetic field orientation, the heat conduction may determine a significant energy exchange between the cloud and the hot surrounding medium which, while remaining always at levels less than those in the unmagnetized case, leads to a progressive heating and evaporation of the cloud. This additional heating may contrast the radiative cooling of some parts of the cloud, preventing the onset of thermal instabilities.

AIMS: We investigate whether the morphology of bilateral supernova remnants (BSNRs) observed in the radio band is determined mainly either by a non-uniform interstellar medium (ISM) or by a non-uniform ambient magnetic field. METHODS: We perform 3-D MHD simulations of a spherical SNR shock propagating through a magnetized ISM. Two cases of shock propagation are considered: 1) through a gradient of ambient density with a uniform ambient magnetic field; 2) through a homogeneous medium with a gradient of ambient magnetic field strength. From the simulations, we synthesize the synchrotron radio emission, making different assumptions about the details of acceleration and injection of relativistic electrons. RESULTS: We find that asymmetric BSNRs are produced if the line-of-sight is not aligned with the gradient of ambient plasma density or with the gradient of ambient magnetic field strength. We derive useful parameters to quantify the degree of asymmetry of the remnants that may provide a powerful diagnostic of the microphysics of strong shock waves through the comparison between models and observations. CONCLUSIONS: BSNRs with two radio limbs of different brightness can be explained if a gradient of ambient density or, most likely, of ambient magnetic field strength is perpendicular to the radio limbs. BSNRs with converging similar radio arcs can be explained if the gradient runs between the two arcs.

Detailed hydrodynamic modeling explained several features of a fragment of the core of a Coronal Mass Ejection observed with SoHO/UVCS at 1.7 Ro on 12 December 1997, but some questions remained unsolved. We investigate the role of the magnetic fields in the thermal insulation and the expansion of an ejected fragment (cloud) traveling upwards in the outer corona. We perform MHD simulations including the effects of thermal conduction and radiative losses of a dense spherical or cylindrical cloud launched upwards in the outer corona, with various assumptions on the strength and topology of the ambient magnetic field; we also consider the case of a cylindrical cloud with an internal magnetic field component along its axis. We find that a weak ambient magnetic field (beta~20) with open topology provides both significant thermal insulation and large expansion. The cylindrical cloud expands more than the spherical one.

Recently X-ray emission from protostellar jets has been detected with both XMM-Newton and Chandra satellites, but the physical mechanism which can give rise to this emission is still unclear. We performed an extensive exploration of a wide space of the main parameters influencing the jet/ambient interaction. Aims include: 1) to constrain the jet/ambient interaction regimes leading to the X-ray emission observed in Herbig-Haro objects in terms of the emission by a shock forming at the interaction front between a continuous supersonic jet and the surrounding medium; 2) to derive detailed predictions to be compared with optical and X-ray observations of protostellar jets; 3) to get insight into the protostellar jet's physical conditions. We performed a set of bidimensional hydrodynamic numerical simulations, in cylindrical coordinates, modeling supersonic jets ramming a uniform ambient medium. The model takes into account the most relevant physical effects, namely the thermal conduction and the radiative losses. Our model explains the observed X-ray emission from protostellar jets in a natural way. In particular we find that the case of a protostellar jet less dense than the ambient medium reproduces well the observations of the nearest Herbig-Haro object, HH154, and allows us to make detailed predictions of a possible X-ray source proper motion (vsh = 500 km/s), detectable with Chandra. Furthermore our results suggest that the simulated protostellar jets which best reproduce the X-rays observations cannot drive molecular outflows.

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.

AIMS. We study and discuss the time-dependent X-ray emission predicted by hydrodynamic modeling of the interaction of a SNR shock wave with an interstellar gas cloud. The scope includes: 1) to study the correspondence between modeled and X-ray emitting structures, 2) to explore two different physical regimes in which either thermal conduction or radiative cooling plays a dominant role, and 3) to investigate the effects of the physical processes at work on the emission of the shocked cloud in the two different regimes. METHODS. We use a detailed hydrodynamic model, including thermal conduction and radiation, and explore two cases characterized by different Mach numbers of the primary shock: M = 30 in which the cloud dynamics is dominated by radiative cooling and M = 50 dominated by thermal conduction. From the simulations, we synthesize the expected X-ray emission, using available spectral codes. RESULTS. The morphology of the X-ray emitting structures is significantly different from that of the flow structures originating from the shock-cloud interaction. The hydrodynamic instabilities are never clearly visible in the X-ray band. Shocked clouds are preferentially visible during the early phases of their evolution. Thermal conduction and radiative cooling lead to two different phases of the shocked cloud: a cold cooling dominated core emitting at low energies and a hot thermally conducting corona emitting in the X-ray band. The thermal conduction makes the X-ray image of the cloud smaller, more diffuse, and shorter-lived than that observed when thermal conduction is neglected.

Magnetic field dependent transient solar observations are suggestive for axion-photon oscillations with light axion(-like) particle involvement. Novel dark-moon measurements with the SMART X-ray detectors can be conclusive for radiatively decaying massive exotica like the generic solar Kaluza-Klein axions. Furthermore, the predicted intrinsic strong solar magnetic fields could be the reason of enhanced low energy axion production. Such an axion component could be the as yet unknown origin of the strong quiet Sun X-ray luminosity at energies below 1 keV. Solar axion telescopes should lower their threshold, aiming to copy processes that might occur near the solar surface, be it due to spontaneous or magnetically induced radiative decay of axion(-like) particles. This is motivated also by the recent claim of an axion-like particle detection by the laser experiment PVLAS.

We have performed, in October 2005, a deep Chandra X-ray observation of HH 154. Comparison with the previous (2001) Chandra observation allows to detect proper motion down to the level predicted by models of X-ray emitting shocks in the jet. The 2005 Chandra observation of HH 154 shows unexpected morphological changes of the X-ray emission in comparison with the 2001 data. Two components are present: a stronger, point-like component with no detectable motion and a weaker component which has expanded in size by approximately 300 AU over the 4 years time base of the two observations. This expansion corresponds to approximately 500 km/s, very close to the velocity of the X-ray emitting shock in the simple theoretical models. The 2005 data show a more complex system than initially thought (and modeled), with multiple components with different properties. The observed morphology is possibly indicating a pulsed jet propagating through a non-homogeneous medium, likely with medium density decreasing with distance from the driving source. Detailed theoretical modeling and deeper X-ray observations will be needed to understand the physics of this fascinating class of sources.

We model the hydrodynamic interaction of a shock wave of an evolved supernova remnant with a small interstellar gas cloud like the ones observed in the Cygnus loop and in the Vela SNR. We investigate the interplay between radiative cooling and thermal conduction during cloud evolution and their effect on the mass and energy exchange between the cloud and the surrounding medium. Through the study of two cases characterized by different Mach numbers of the primary shock (M = 30 and 50, corresponding to a post-shock temperature $T\approx 1.7\times 10^6$ K and $\approx 4.7\times 10^6$ K, respectively), we explore two very different physical regimes: for M = 30, the radiative losses dominate the evolution of the shocked cloud which fragments into cold, dense, and compact filaments surrounded by a hot corona which is ablated by the thermal conduction; instead, for M = 50, the thermal conduction dominates the evolution of the shocked cloud, which evaporates in a few dynamical time-scales. In both cases we find that the thermal conduction is very effective in suppressing the hydrodynamic instabilities that would develop at the cloud boundaries.

We show that stellar coronae can be composed of X-ray emitting structures like those in the solar corona, using a large set of ROSAT/PSPC observations of late-type-stars, and a large set of solar X-ray data collected with Yohkoh/SXT. We have considered data on the solar corona at various phases of the cycle and various kinds of X-ray coronal structures, from flares to the background corona. The surface flux (F_x) vs. spectral hardness ratio (HR) diagram is a fundamental tool for our study. We find that $F_X$ is strongly correlated to HR in stellar coronae, in the solar corona at all phases of the cycle, and in the individual solar coronal structures; all follow the same law. We therefore claim that coronae of late type stars are formed with X-ray structures very similar to the solar ones. In this scenario, the fraction of the stellar surface covered with active regions and with their bright cores increases with activity; the most active stars are brighter and hotter than if they were entirely covered with active regions so they can be explained only with the additional presence of several flares (or flare-like structures) at any time. Using the F_x vs. HR correlation (F_X ~ T^6) we derive new laws relating the temperature, pressure, volumetric heating and characteristic loop length of the coronal plasma, on all the late type stars Individual solar coronal structures and the whole solar corona, follow the same laws. We claim that the strong correlation between surface flux and temperature and the laws mentioned above are just the effect of more fundamental physical mechanisms driving the coronal structures of all the late-type stars from the emergence of new magnetic structures to their dissipation.

The aim of this paper is to validate a methodology for connecting the emission measure of individual solar coronal loops to the integrated emission measure of the entire solar corona, and using this connection to deduce the energetic properties of the corona, and then to show how this methodology can be applied to observations of solar-like stellar coronae. The solar validation is carried out by using spatially resolved X-ray observations of the Sun obtained from the Yohkoh satellite. This work is a further step in our effort to place the "solar-stellar connection" on a quantitative footing. In particular, we show how this analysis procedure can be used in the context of archival Einstein, ROSAT and EUVE data, as well as Chandra and XMM Newton data, as a complementary analysis tool to existing multi-thermal component models.

In previous works we have developed a method to convert solar X-ray data, collected with the Yohkoh/SXT, into templates of stellar coronal observations. Here we apply the method to several solar flares, for comparison with stellar X-ray flares. Eight flares, from weak (GOES class C5.8) to very intense ones (X9) are selected as representative of the flaring Sun. The emission measure distribution vs. temperature, EM(T), of the flaring regions is derived from Yohkoh/SXT observations in the rise, peak and decay of the flares. The EM(T) is rather peaked and centered around $T \approx 10^7$ K for most of the time. Typically, it grows during the rise phase of the flare, and then it decreases and shifts toward lower temperatures during the decay, more slowly if there is sustained heating. The most intense flare we studied shows emission measure even at very high temperature ($T \approx 10^8$ K). Time-resolved X-ray spectra both unfiltered and filtered through the instrumental responses of the non-solar instruments ASCA/SIS and ROSAT/PSPC are then derived. Synthesized ASCA/SIS and ROSAT/PSPC spectra are generally well fitted with single thermal components at temperatures close to that of the EM(T) maximum, albeit two thermal components are needed to fit some flare decays. ROSAT/PSPC spectra show that solar flares are in a two-orders of magnitude flux range ($10^6 - 10^8$ erg cm$^{-2}$ s$^{-1}$) and a narrow PSPC hardness ratio range, however higher than that of typical non-flaring solar-like stars.