BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

S. Orlando

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

ResearcherUnclaimed source record
astro-ph.HEastro-ph.SRastro-ph.GAphysics.plasm-phastro-phhep-ph

Catalog footprint

What is connected

49works
6topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

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

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

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

The X-ray activity of the young solar-like star Kepler-63 and the structure of its corona

The X-ray satellite XMM-Newton has so far revealed coronal cycles in seven solar-like stars. In this sample, the youngest stars $ε$ Eridani (400 Myr) and $ι$ Horologii (600 Myr) display the shortest X-ray cycles and the smallest cycle amplitudes. The corona of $ε$ Eridani was modelled in terms of solar magnetic structures (active regions, cores of active regions and flares) at varying filling factors. The study revealed that 65-95% of its corona is covered with magnetic structures, and this was held responsible for the low X-ray cycle amplitude. It was also hypothesized that the basal surface coverage with magnetic structures may be higher on the corona of younger solar-like stars. To investigate this hypothesis, we study the solar-like star Kepler-63 in the X-rays. With an age of 210 Myr and a photospheric cycle of 1.27 yr, it is so far the youngest star observed in X-rays with the aim of revealing a coronal cycle. In the long-term X-ray lightcurve we do not reveal a periodic variation of the X-ray luminosity, but a factor two change would be possible. As for the case of $ε$ Eridani, we modelled the corona of Kepler-63 with magnetic structures observed on the Sun. The study suggests that 100% of the corona is composed of cores and flares of Class M, justifying the absence of an X-ray cycle and confirming the analogous results derived for $ε$ Eridani. Finally, we establish an empirical relation between the cycle amplitude and the X-ray surface flux . From the absence of a coronal cycle in Kepler-63 we infer that stars with higher X-ray flux than Kepler-63 must host a significant fraction of higher-energetic flares than those necessary to model the corona of Kepler-63. Our study opens new ground for studies of the solar-stellar analogy and the joint exploration of resolved and unresolved variability in stellar X-ray lightcurves.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Characterizing X-ray activity cycles of young solar-like stars with solar observations

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

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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 ]"

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

An X-ray activity cycle on the young solar-like star $ε\ \rm Eridani$

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Laboratory disruption of scaled astrophysical outflows by a misaligned magnetic field

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 from our Sun's outflows to extragalatic jets. Furthermore, they provide a possible interpretation for the observed structuring of astrophysical jets. Jet modulation could be interpreted as the signature of changes over time in the outflow/ambient field angle, and the change in the direction of the jet could be the signature of changes in the direction of the ambient field.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Predicting the time variation of radio emission from MHD simulations of a flaring T-Tauri star

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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

3D MHD modeling of twisted coronal loops

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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Bright hot impacts by erupted fragments falling back on the Sun: magnetic channelling

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

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Bright hot impacts by erupted fragments falling back on the Sun: UV redshifts in stellar accretion

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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Mass Accretion Processes in Young Stellar Objects: Role of Intense Flaring Activity

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

MHD modeling of coronal loops: injection of high-speed chromospheric flows

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

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

MHD modeling of coronal loops: the transition region throat

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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

YSO accretion shocks: magnetic, chromospheric or stochastic flow effects can suppress fluctuations of X-ray emission

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

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Mass accretion to young stars triggered by flaring activity in circumstellar disks

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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Modeling the 2010 blast wave of the symbiotic-like nova V407 Cygni

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

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Post-flare UV light curves explained with thermal instability of loop plasma

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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Multi-wavelength diagnostics of accretion in an X-ray selected sample of CTTSs

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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Radio, X-ray and gamma-ray surface brightness profiles as powerful diagnostic tools for non-thermal SNR shells

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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

The X-ray emission from Z CMa during an FUor-like outburst and the detection of its X-ray jet

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.

0 citations0 reviewsNo code linkedOpen access
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

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

X-ray emitting MHD accretion shocks in classical T Tauri stars. Case for moderate to high plasma-beta values

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.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

The X-ray cycle in the solar-type star HD 81809

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

0 citations0 reviewsNo code linkedOpen access