Source author record

M. Miceli

M. Miceli appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

39works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

39 published item(s)

preprint2025arXiv

Predicting the X-ray signatures of the imminent T Coronae Borealis outburst through 3D hydrodynamic modeling

T Coronae Borealis (T CrB) is a symbiotic recurrent nova with eruptions in 1866 and 1946. Mounting evidence suggests an imminent outburst, offering a rare opportunity to observe a nearby nova in detail. We constrain the circumbinary medium (CBM) by modeling inter-eruption radio observations and simulate the hydrodynamic evolution of the upcoming outburst to predict its X-ray signatures, focusing on the roles of the red giant companion, accretion disk, and equatorial density enhancement (EDE). We model thermal radio emission from a CBM composed of a spherical wind and a torus-like EDE to estimate its density. We then perform 3D hydrodynamic simulations of the nova outburst, varying explosion energy, ejecta mass, and CBM configuration. From these, we synthesize X-ray light curves and spectra as they would appear to XMM-Newton and XRISM. The CBM in T CrB is significantly less dense than in other symbiotic novae, with a mass-loss rate of $\dot{M} \approx 4 \times 10^{-9}$ M$_{\odot}$ yr$^{-1}$ for a 10 km s$^{-1}$ wind. Despite the low-density EDE, the blast is collimated along the poles by the accretion disk and EDE, producing a bipolar shock. The red giant partially shields the ejecta, forming a bow shock and hot wake. X-ray evolution proceeds through three phases: an early phase (first few hours) dominated by shocked disk material; an intermediate phase ($\sim 1$ week-1 month) driven by reverse-shocked ejecta; and a late phase dominated by shocked EDE. Soft X-rays trace shocked ejecta, hard X-rays arise from shocked ambient gas, and synthetic spectra show asymmetric, blueshifted lines due to absorption by expanding ejecta. The X-ray evolution resembles that of RS Oph and V745 Sco, with a peak luminosity of $L_\mathrm{X} \approx 10^{36}$ erg s$^{-1}$, but features a more prolonged soft X-ray phase, reflecting the lower CBM density and distinct ejecta-environment interaction.

preprint2022arXiv

High energy gamma-ray detection of supernova remnants in the Large Magellanic Cloud

We present the results of a cluster search in the gamma-ray sky images of the Large Magellanic Cloud (LMC) region by means of the Minimum Spanning Tree (MST) and DBSCAN algorithms, at energies higher than 6 and 10 GeV, using 12 years of Fermi-LAT data. Several significant clusters were found, the majority of which associated with previously known gamma-ray sources. We confirm our previous detection of the Supernova Remnants N 49B and N 63A and found new significant clusters associated with the SNRs N 49, N 186D and N 44. These sources are among the brightest X-ray remnants in the LMC and corresponds to core-collapse supernovae interacting with dense HII regions, indicating that an hadronic origin of high energy photons is the most likely process.

preprint2022arXiv

Negative and Positive Feedback from a Supernova Remnant with SHREC: A detailed Study of the Shocked Gas in IC443

Supernova remnants (SNRs) contribute to regulate the star formation efficiency and evolution of galaxies. As they expand into the interstellar medium (ISM), they transfer vast amounts of energy and momentum that displace, compress and heat the surrounding material. Despite the extensive work in galaxy evolution models, it remains to be observationally validated to what extent the molecular ISM is affected by the interaction with SNRs. We use the first results of the ESO-ARO Public Spectroscopic Survey SHREC, to investigate the shock interaction between the SNR IC443 and the nearby molecular clump G. We use high sensitivity SiO(2-1) and H$^{13}$CO$^+$(1-0) maps obtained by SHREC together with SiO(1-0) observations obtained with the 40m telescope at the Yebes Observatory. We find that the bulk of the SiO emission is arising from the ongoing shock interaction between IC443 and clump G. The shocked gas shows a well ordered kinematic structure, with velocities blue-shifted with respect to the central velocity of the SNR, similar to what observed toward other SNR-cloud interaction sites. The shock compression enhances the molecular gas density, n(H$_2$), up to $>$10$^5$ cm$^{-3}$, a factor of >10 higher than the ambient gas density and similar to values required to ignite star formation. Finally, we estimate that up to 50\% of the momentum injected by IC443 is transferred to the interacting molecular material. Therefore the molecular ISM may represent an important momentum carrier in sites of SNR-cloud interactions.

preprint2021arXiv

Detailed characterization of laboratory magnetized super-critical collisionless shock and of the associated proton energization

Collisionless shocks are ubiquitous in the Universe and are held responsible for the production of non-thermal particles and high-energy radiation. In the absence of particle collisions in the system, theoretical works show that the interaction of an expanding plasma with a pre-existing electromagnetic structure (as in our case) is able to induce energy dissipation and allow for shock formation. Shock formation can alternatively take place when two plasmas interact, through microscopic instabilities inducing electromagnetic fields which are able in turn to mediate energy dissipation and shock formation. Using our platform where we couple a fast-expanding plasma induced by high-power lasers (JLF/Titan at LLNL and LULI2000) with high-strength magnetic fields, we have investigated the generation of magnetized collisionless shock and the associated particle energization. We have characterized the shock to be collisionless and super-critical. We report here on measurements of the plasma density, temperature, the electromagnetic field structures, and particle energization in the experiments, under various conditions of ambient plasma and B-field. We have also modeled the formation of the shocks using macroscopic hydrodynamic simulations and the associated particle acceleration using kinetic particle-in-cell simulations. As a companion paper of \citet{yao2020laboratory}, here we show additional results of the experiments and simulations, providing more information to reproduce them and demonstrating the robustness of our interpreted proton energization mechanism to be shock surfing acceleration.

preprint2021arXiv

Sensitivity of the Cherenkov Telescope Array to a dark matter signal from the Galactic centre

We provide an updated assessment of the power of the Cherenkov Telescope Array (CTA) to search for thermally produced dark matter at the TeV scale, via the associated gamma-ray signal from pair-annihilating dark matter particles in the region around the Galactic centre. We find that CTA will open a new window of discovery potential, significantly extending the range of robustly testable models given a standard cuspy profile of the dark matter density distribution. Importantly, even for a cored profile, the projected sensitivity of CTA will be sufficient to probe various well-motivated models of thermally produced dark matter at the TeV scale. This is due to CTA's unprecedented sensitivity, angular and energy resolutions, and the planned observational strategy. The survey of the inner Galaxy will cover a much larger region than corresponding previous observational campaigns with imaging atmospheric Cherenkov telescopes. CTA will map with unprecedented precision the large-scale diffuse emission in high-energy gamma rays, constituting a background for dark matter searches for which we adopt state-of-the-art models based on current data. Throughout our analysis, we use up-to-date event reconstruction Monte Carlo tools developed by the CTA consortium, and pay special attention to quantifying the level of instrumental systematic uncertainties, as well as background template systematic errors, required to probe thermally produced dark matter at these energies. "Full likelihood tables complementing our analysis are provided here [ https://doi.org/10.5281/zenodo.4057987 ]"

preprint2020arXiv

Modeling particle acceleration and non-thermal emission in supernova remnants

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

preprint2020arXiv

The fully developed remnant of a neutrino-driven supernova: Evolution of ejecta structure and asymmetries in SNR Cassiopeia A

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

preprint2019arXiv

Hydrodynamic simulations unravel the progenitor-supernova-remnant connection in SN 1987A

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

preprint2016arXiv

Detailed study of SNR G306.3-0.9 using XMM-Newton and Chandra observations

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

preprint2016arXiv

Formation of X-ray emitting stationary shocks in magnetized protostellar jets

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.

preprint2016arXiv

Modeling non-thermal emission from stellar bow shocks

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

preprint2016arXiv

Modeling SNR Cassiopeia A from the Supernova Explosion to its Current Age: The role of post-explosion anisotropies of ejecta

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.

preprint2015arXiv

Spatial distribution of X-ray emitting ejecta in Tycho's SNR: indications of shocked Titanium

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.

preprint2015arXiv

Supernova 1987A: a Template to Link Supernovae to their Remnants

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 $α= -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.

preprint2014arXiv

Magnetohydrodynamic modeling of the accretion shocks in classical T Tauri stars: the role of local absorption on the X-ray emission

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.

preprint2014arXiv

Shock-cloud interaction and particle acceleration in SN 1006

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.

preprint2013arXiv

Evidence of non-thermal X-ray emission from HH 80

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

preprint2013arXiv

Hydrodynamic modelling of ejecta shrapnel in the Vela supernova remnant

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.

preprint2013arXiv

Radiative accretion shocks along nonuniform stellar magnetic fields in classical T Tauri stars

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

preprint2013arXiv

The Hot and Energetic Universe: The astrophysics of supernova remnants and the interstellar medium

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

preprint2013arXiv

The loss-limited electron energy in SN 1006: effects of the shock velocity and of the diffusion process

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.

preprint2013arXiv

The shape of the cutoff in the synchrotron emission of SN 1006 observed with XMM-Newton

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.

preprint2012arXiv

A population of isolated hard X-ray sources near the supernova remnant Kes 69

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

preprint2012arXiv

AE Aurigae: first detection of non-thermal X-ray emission from a bow shock produced by a runaway star

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

preprint2012arXiv

Nanoflare Evidence from Analysis of the X-Ray Variability of an Active Region Observed with Hinode/XRT

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

preprint2012arXiv

Role of ejecta clumping and back-reaction of accelerated cosmic rays in the evolution of Type Ia 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.

preprint2012arXiv

X-raying hot plasma in solar active regions with the SphinX spectrometer

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

preprint2012arXiv

XMM-Newton evidence of shocked ISM in SN 1006: indications of hadronic acceleration

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.

preprint2011arXiv

Constraints on local interstellar magnetic field from non-thermal emission of SN1006

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

preprint2011arXiv

Role of ejecta clumping and back-reaction of accelerated cosmic rays in the evolution of supernova remnants

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.

preprint2011arXiv

X-ray emission from protostellar jet HH 154: the first evidence of a diamond shock?

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.

preprint2010arXiv

Effects of non-uniform interstellar magnetic field on synchrotron X-ray and inverse-Compton gamma-ray morphology of SNRs

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.

preprint2010arXiv

Generation of radiative knots in a randomly pulsed protostellar jet II. X-ray emission

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.

preprint2010arXiv

Observability and diagnostics in the X-ray band of shock-cloud interactions in supernova remnants

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.

preprint2010arXiv

Observational constraints on the modeling of SN1006

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.

preprint2010arXiv

Spatial identification of the overionized plasma in W49B

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

preprint2009arXiv

Predicted gamma-ray image of SN 1006 due to inverse Compton emission

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.

preprint2009arXiv

Thermal emission, shock modification, and X-ray emitting ejecta in SN 1006

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