Search

The maximum energy of electrons accelerated by supernova remnants (SNR) is typically limited by radiative losses. In this scenario, the synchrotron cooling time scale is equal to the acceleration time scale. On the other hand, the low propagation speed of a shock in a dense medium is expected to result in an extended acceleration time scale, thus inducing a decrease in the maximum electron energy for a given SNR age and in the X-ray nonthermal flux. The young Kepler's SNR shows an enhanced efficiency of the acceleration process, which is close to the Bohm limit in the north of its shell, where the shock is slowed down by a dense circumstellar medium. Conversely, in the south, where no interaction with a dense medium is evident and the shock speed is high, the acceleration proceeds with a higher Bohm factor. To investigate this scenario, we studied the temporal evolution of the non-thermal emission, taking advantage of two Chandra X-ray observations of Kepler's SNR (performed in 2006 and 2014). We analyzed the spectra of different filaments both in the north and south of the shell, and measured their proper motion. We found a region with low shock velocity where we measured a significant decrease in flux from 2006 to 2014. This could be the first evidence of fading synchrotron emission in Kepler's SNR. This result suggests that under a certain threshold of shock speed the acceleration process could exit the loss-limited regime., Comment: Presented at the 38th International Cosmic Ray Conference (ICRC 2023), 8 pages, 2 figures, 2 tables

The origin of cosmic rays is a pivotal open issue of high-energy astrophysics. Supernova remnants are strong candidates to be the Galactic factory of cosmic rays, their blast waves being powerful particle accelerators. However, supernova remnants can power the observed flux of cosmic rays only if they transfer a significant fraction of their kinetic energy to the accelerated particles, but conclusive evidence for such efficient acceleration is still lacking. In this scenario, the shock energy channeled to cosmic rays should induce a higher post-shock density than that predicted by standard shock conditions. Here we show this effect, and probe its dependence on the orientation of the ambient magnetic field, by analyzing deep X-ray observations of the Galactic remnant of SN 1006. By comparing our results with state-of-the-art models, we conclude that SN 1006 is an efficient source of cosmic rays and obtain an observational support for the quasi-parallel acceleration mechanism., Comment: Published in Nature Communications. The published version of the paper, including supplementary material, is freely available online here https://doi.org/10.1038/s41467-022-32781-4 or as a PDF here https://www.nature.com/articles/s41467-022-32781-4.pdf

Synchrotron X-ray emission in young supernova remnants (SNRs) is a powerful diagnostic tool to study the population of high energy electrons accelerated at the shock front and the acceleration process. We performed a spatially resolved spectral analysis of NuSTAR and XMM-Newton observations of the young Kepler's SNR, aiming to study in detail its non-thermal emission in hard X-rays. We selected a set of regions all around the rim of the shell and extracted the corresponding spectra. The spectra were analyzed by adopting a model of synchrotron radiation in the loss-limited regime, to constrain the dependence of the cutoff energy of the synchrotron radiation on the shock velocity. We identify two different regimes of particle acceleration, characterized by different Bohm factors. In the north, where the shock interacts with a dense circumstellar medium (CSM), we found a more efficient acceleration than in the south, where the shock velocity is higher and there are no signs of shock interaction with dense CSM. Our results suggest an enhanced efficiency of the acceleration process in regions where the shock-CSM interaction generates an amplified and turbulent magnetic field. By combining hard X-ray spectra with radio and $\gamma-$ray observations of Kepler's SNR, we modelled the spectral energy distribution. In the light of our results we propose that the observed $\gamma-$ray emission is mainly hadronic, and originates in the northern part of the shell., Comment: 16 pages, 7 figures, 5 tables, Accepted for publication in ApJ

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

Ariel, the Atmospheric Remote-sensing Infrared Exoplanet Large-survey, was adopted as the fourth medium-class mission in ESA's Cosmic Vision programme to be launched in 2029. During its 4-year mission, Ariel will study what exoplanets are made of, how they formed and how they evolve, by surveying a diverse sample of about 1000 extrasolar planets, simultaneously in visible and infrared wavelengths. It is the first mission dedicated to measuring the chemical composition and thermal structures of hundreds of transiting exoplanets, enabling planetary science far beyond the boundaries of the Solar System. The payload consists of an off-axis Cassegrain telescope (primary mirror 1100 mm x 730 mm ellipse) and two separate instruments (FGS and AIRS) covering simultaneously 0.5-7.8 micron spectral range. The satellite is best placed into an L2 orbit to maximise the thermal stability and the field of regard. The payload module is passively cooled via a series of V-Groove radiators; the detectors for the AIRS are the only items that require active cooling via an active Ne JT cooler. The Ariel payload is developed by a consortium of more than 50 institutes from 16 ESA countries, which include the UK, France, Italy, Belgium, Poland, Spain, Austria, Denmark, Ireland, Portugal, Czech Republic, Hungary, the Netherlands, Sweden, Norway, Estonia, and a NASA contribution., Comment: Ariel Definition Study Report, 147 pages. Reviewed by ESA Science Advisory Structure in November 2020. Original document available at: https://www.cosmos.esa.int/documents/1783156/3267291/Ariel_RedBook_Nov2020.pdf/

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

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

The goal of this study is to assess the impact of the stellar spots on the extraction of the planetary transmission spectra observed by ARIEL. We develop a method to model the stellar spectrum of a star in the presence of spots by using the out-of-transit observations. It is based on a chi squared minimization procedure of the out-of-transit spectrum on a grid of stellar spectra with different sizes and temperatures of the spots. The approach allows us also to study the temporal evolution of the spots when comparing stellar spectra observed at different epochs. We also present a method to correct the transit depth variations due to non-occulted stellar spots and estimate the error we introduce if we apply the same correction to crossings over the stellar spots. The method is tested on three types of stellar targets that ARIEL will observe in its 4-years mission lifetime. In all the explored cases, the approach allows to reliably recover the spot parameters (size and temperature) from out-of-transit observations and, for non-occulted spots, to confidently recover the planetary atmosphere transmission spectrum within the noise level(with average uncertainty of at most 3.3% of the planetary signal). Conversely, we find systematic biases in the inferred planetary spectra due to the occulted spots, with measurable effects for the brightest targets especially for more contrasted spots., Comment: 16 pages, 13 figures

The aim of this work is to bridge the gap between CC SNe and their remnants by investigating how post-explosion anisotropies in the ejecta influence the structure and chemical properties of the remnant at later times. We performed three-dimensional magneto-hydrodynamical simulations starting soon after the SN event and following the evolution of the system in the circumstellar medium (consisting of the wind of the stellar progenitor), for 5000 years, obtaining the physical scenario of a SNR. Here we focused the analysis on the case of a progenitor red supergiant of 19.8 M_sun. We also investigated how a post-explosion large-scale anisotropy in the SN affects the ejecta distribution and the matter mixing of heavy elements in the remnant, during the first 5000 years of evolution. In the case of a spherically symmetric SN explosion without large-scale anisotropies, the remnant roughly keeps memory of the original onion-like layering of ejecta soon after the SN event. Nevertheless, as the reverse shock hits the ejecta, the element distribution departs from a homologous expansion, because of the slowing down of the outermost ejecta layers due to interaction with the reverse shock. In the case of a large-scale anisotropy developed after the SN, we found that the chemical stratification in the ejecta can be strongly modified and the original onion-like layering is not preserved. The anisotropy may cause spatial inversion of ejecta layers, for instance leading to Fe/Si-rich ejecta outside the O shell, and may determine the formation of Fe/Si-rich jet-like features that may protrude the remnant outline. The level of matter mixing and the properties of the jet-like feature are sensitive to the initial physical (density and velocity) and geometrical (size and position) initial characteristics of the anisotropy., Comment: 16 pages, 14 figures, accepted to A&A

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

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

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

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

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

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

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

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

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

We explore UV and optical variability signatures for several hundred members of NGC 2264 (3 Myr). We performed simultaneous u- and r-band monitoring over two full weeks with CFHT/MegaCam. About 750 young stars are probed; 40% of them are accreting. Statistically distinct variability properties are observed for accreting and non-accreting cluster members. The accretors exhibit a significantly higher level of variability than the non-accretors, especially in the UV. The amount of u-band variability correlates statistically with UV excess in disk-bearing objects, which suggests that accretion and star-disk interaction are the main sources of variability. Cool magnetic spots, several hundred degrees colder than the photosphere and covering from 5 to 30% of the stellar surface, appear to be the leading factor of variability for the non-accreting stars. In contrast, accretion spots, a few thousand degrees hotter than the photosphere and covering a few percent of the stellar surface, best reproduce the variability of accreting objects. The color behavior is also found to be different between accreting and non-accreting stars. Typical variability amplitudes for accreting members rapidly increase from r to u, which indicates a much stronger contrast at short wavelengths; a lower color dependence in the amplitudes is instead measured for diskless stars. We find that u-band variability on hour timescales is typically about 10% of the peak-to-peak variability on day timescales, while longer term (years) variability is consistent with amplitudes measured over weeks. We conclude that for both accreting and non-accreting stars, the mid-term rotational modulation by spots is the leading timescale for a variability of up to several years. In turn, this suggests that the accretion process is essentially stable over years, although it exhibits low-level shorter term variations in single accretion events., Comment: 14 pages, 13 figures, accepted for publication in Astronomy & Astrophysics

We aim at characterizing the accretion properties of several hundred members of the star-forming cluster NGC 2264 (3 Myr). We performed a deep u,g,r,i mapping and a simultaneous u+r monitoring of the region with CFHT/MegaCam in order to directly probe the accretion process from UV excess measurements. Photometric properties and stellar parameters are determined homogeneously for about 750 monitored young objects, spanning the mass range 0.1-2 Mo. About 40% are classical (accreting) T Tauri stars, based on various diagnostics (H_alpha, UV and IR excesses). The remaining non-accreting members define the (photospheric+chromospheric) reference UV emission level over which flux excess is detected and measured. We revise the membership status of cluster members based on UV accretion signatures and report a new population of 50 CTTS candidates. A large range of UV excess is measured for the CTTS population, varying from a few 0.1 to 3 mag. We convert these values to accretion luminosities and obtain mass accretion rates ranging from 1e-10 to 1e-7 Mo/yr. Taking into account a mass-dependent detection threshold for weakly accreting objects, we find a >6sigma correlation between mass accretion rate and stellar mass. A power-law fit, properly accounting for upper limits, yields M_acc $\propto$ M^{1.4+/-0.3}. At any given stellar mass, we find a large spread of accretion rates, extending over about 2 orders of magnitude. The monitoring of the UV excess on a timescale of a couple of weeks indicates that its variability typically amounts to 0.5 dex, much smaller than the observed spread. We suggest that a non-negligible age spread across the cluster may effectively contribute to the observed spread in accretion rates at a given mass. In addition, different accretion mechanisms (like, e.g., short-lived accretion bursts vs. more stable funnel-flow accretion) may be associated to different M_acc regimes., Comment: 24 pages, 21 figures, accepted for publication in Astronomy & Astrophysics

Blind source separation techniques are used to reanalyse two exoplanetary transit lightcurves of the exoplanet HD189733b recorded with the IR camera IRAC on board the Spitzer Space Telescope at 3.6$\mu$m during the "cold" era. These observations, together with observations at other IR wavelengths, are crucial to characterise the atmosphere of the planet HD189733b. Previous analyses of the same datasets reported discrepant results, hence the necessity of the reanalyses. The method we used here is based on the Independent Component Analysis (ICA) statistical technique, which ensures a high degree of objectivity. The use of ICA to detrend single photometric observations in a self-consistent way is novel in the literature. The advantage of our reanalyses over previous work is that we do not have to make any assumptions on the structure of the unknown instrumental systematics. Such "admission of ignorance" may result in larger error bars than reported in the literature, up to a factor $1.6$. This is a worthwhile trade-off for much higher objectivity, necessary for trustworthy claims. Our main results are (1) improved and robust values of orbital and stellar parameters, (2) new measurements of the transit depths at 3.6$\mu$m, (3) consistency between the parameters estimated from the two observations, (4) repeatability of the measurement within the photometric level of $\sim 2 \times 10^{-4}$ in the IR, (5) no evidence of stellar variability at the same photometric level within 1 year., Comment: 43 pages, 18 figures

X-ray polarimetry, sometimes alone, and sometimes coupled to spectral and temporal variability measurements and to imaging, allows a wealth of physical phenomena in astrophysics to be studied. X-ray polarimetry investigates the acceleration process, for example, including those typical of magnetic reconnection in solar flares, but also emission in the strong magnetic fields of neutron stars and white dwarfs. It detects scattering in asymmetric structures such as accretion disks and columns, and in the so-called molecular torus and ionization cones. In addition, it allows fundamental physics in regimes of gravity and of magnetic field intensity not accessible to experiments on the Earth to be probed. Finally, models that describe fundamental interactions (e.g. quantum gravity and the extension of the Standard Model) can be tested. We describe in this paper the X-ray Imaging Polarimetry Explorer (XIPE), proposed in June 2012 to the first ESA call for a small mission with a launch in 2017 but not selected. XIPE is composed of two out of the three existing JET-X telescopes with two Gas Pixel Detectors (GPD) filled with a He-DME mixture at their focus and two additional GPDs filled with pressurized Ar-DME facing the sun. The Minimum Detectable Polarization is 14 % at 1 mCrab in 10E5 s (2-10 keV) and 0.6 % for an X10 class flare. The Half Energy Width, measured at PANTER X-ray test facility (MPE, Germany) with JET-X optics is 24 arcsec. XIPE takes advantage of a low-earth equatorial orbit with Malindi as down-link station and of a Mission Operation Center (MOC) at INPE (Brazil)., Comment: 49 pages, 14 figures, 6 tables. Paper published in Experimental Astronomy http://link.springer.com/journal/10686

One scenario proposed to explain the million degrees solar corona is a finely-stranded corona where each strand is heated by a rapid pulse. However, such fine structure has neither been resolved through direct imaging observations nor conclusively shown through indirect observations of extended superhot plasma. Recently it has been shown that the observed difference in appearance of cool and warm coronal loops (~1 MK, ~2-3 MK, respectively) -- warm loops appearing "fuzzier" than cool loops -- can be explained by models of loops composed of subarcsecond strands, which are impulsively heated up to ~10 MK. That work predicts that images of hot coronal loops (>~6 MK) should again show fine structure. Here we show that the predicted effect is indeed widely observed in an active region with the Solar Dynamics Observatory, thus supporting a scenario where impulsive heating of fine loop strands plays an important role in powering the active corona., Comment: 11 pages, 4 figures, accepted for publication

Coronal active regions are observed to get fuzzier and fuzzier (i.e. more and more confused and uniform) in harder and harder energy bands or lines. We explain this evidence as due to the fine multi-temperature structure of coronal loops. To this end, we model bundles of loops made of thin strands, each heated by short and intense heat pulses. For simplicity, we assume that the heat pulses are all equal and triggered only once in each strand at a random time. The pulse intensity and cadence are selected so as to have steady active region loops ($\sim 3$ MK), on the average. We compute the evolution of the confined heated plasma with a hydrodynamic loop model. We then compute the emission along each strand in several spectral lines, from cool ($\leq 1$ MK), to warm ($2-3$ MK) lines, detectable with Hinode/EIS, to hot X-ray lines. The strands are then put side-by-side to construct an active region loop bundle. We find that in the warm lines ($2-3$ MK) the loop emission fills all the available image surface. Therefore the emission appears quite uniform and it is difficult to resolve the single loops, while in the cool lines the loops are considerably more contrasted and the region is less fuzzy. The main reasons for this effect are that, during their evolution, i.e. pulse heating and slow cooling, each strand spends a relatively long time at temperatures around $2-3$ MK, and that it has a high emission measure during that phase, so the whole region appears more uniform or smudged. We make the prediction that the fuzziness should be reduced in the hot UV and X-ray lines., Comment: 27 pages, 14 figures

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

We have investigated the optical thickness of the coronal plasma through the analysis of high-resolution X-ray spectra of a large sample of active stars observed with the High Energy Transmission Grating Spectrometer on Chandra. In particular, we probed for the presence of significant resonant scattering in the strong Lyman series lines arising from hydrogen-like oxygen and neon ions. The active RS CVn-type binaries II Peg and IM Peg and the single M dwarf EV Lac show significant optical depth. For these active coronae, the Lya/Lyb ratios are significantly depleted as compared with theoretical predictions and with the same ratios observed in similar active stars. Interpreting these decrements in terms of resonance scattering of line photons out of the line-of-sight, we are able to derive an estimate for the typical size of coronal structures, and from these we also derive estimates of coronal filling factors. For all three sources we find that the both the photon path length as a fraction of the stellar radius, and the implied surface filling factors are very small and amount to a few percent at most. The measured Lya/Lyb ratios are in good agreement with APED theoretical predictions, thus indicating negligible optical depth, for the other sources in our sample. We discuss the implications for coronal structuring and heating flux requirements. For the stellar sample as a whole, the data suggest increasing quenching of Lya relative to Lyb as function of both L_x/L_bol and the density-sensitive MgXI forbidden to intercombination line ratio, as might generally be expected., Comment: Accepted for publication on the Astrophysical Journal

This work is prompted by the evidence of sharply peaked emission measure distributions in active stars, and by the claims of isothermal loops in solar coronal observations, at variance with the predictions of hydrostatic loop models with constant cross-section and uniform heating. We address the problem with loops heated at the foot-points. Since steady heating does not allow static loop models solutions, we explore whether pulse-heated loops can exist and appear as steady loops, on a time average. We simulate pulse-heated loops, using the Palermo-Harvard 1-D hydrodynamic code, for different initial conditions corresponding to typical coronal temperatures of stars ranging from intermediate to active ($T \sim 3$--$10 \times 10^6$ K). We find long-lived quasi-steady solutions even for heating concentrated at the foot-points over a spatial region of the order of $\sim 1/5$ of the loop half length and broader.These solutions yield an emission measure distribution with a peak at high temperature, and the cool side of the peak is as steep as $\sim T^{5}$, in contrast to the usual $\sim T^{3/2}$ of hydrostatic models with constant cross-section and uniform heating. Such peaks are similar to those found in the emission measure distribution of active stars around $10^7$ K., Comment: 9 pages, 6 figures, accepted for publication in The Astrophysical Journal

An analysis of Lyman series lines arising from hydrogen-like oxygen and neon ions in the coronae of the active RS CVn-type binaries II Peg and IM Peg, observed using the {\it Chandra} High Resolution Transmission Grating Spectrograph, shows significant decrements in the Ly$\alpha$/Ly$\beta$ ratios as compared with theoretical predictions and with the same ratios observed in similar active binaries. We interpret these decrements in terms of resonance scattering of line photons out of the line-of-sight; these observations present the first strong evidence for this effect in active stellar coronae. The net line photon loss implies a non-uniform and asymmetric surface distribution of emitting structures on these stars. Escape probability arguments, together with the observed line ratios and estimates of the emitting plasma density, imply typical line-of-sight sizes of the coronal structures that dominate the X-ray emission of $10^{10}$ cm at temperatures of $3\times 10^6$ K and of $10^8$ cm at $10^7$ K. These sizes are an order of magnitude larger than predicted by simple quasi-static coronal loops models, but are still very small compared to the several $10^{11}$ cm radii of the underlying stars., Comment: Accepted for publication on ApJ Letters. 14 pages, 2 figures

We have analyzed high-resolution X-ray spectra of a sample of 22 active stars observed with the High Energy Transmission Grating Spectrometer on {\em Chandra} in order to investigate their coronal plasma density. Densities where investigated using the lines of the He-like ions O VII, Mg XI, and Si XIII. While Si XIII lines in all stars of the sample are compatible with the low-density limit, Mg XI lines betray the presence of high plasma densities ($> 10^{12}$ cm$^{-3}$) for most of the sources with higher X-ray luminosity ($> 10^{30}$ erg/s); stars with higher $L_X$ and $L_X/L_{bol}$ tend to have higher densities at high temperatures. Ratios of O VII lines yield much lower densities of a few $10^{10}$ cm$^{-3}$, indicating that the ``hot'' and ``cool'' plasma resides in physically different structures. Our findings imply remarkably compact coronal structures, especially for the hotter plasma emitting the Mg XI lines characterized by coronal surface filling factor, $f_{MgXI}$, ranging from $10^{-4}$ to $10^{-1}$, while we find $f_{OVII}$ values from a few $10^{-3}$ up to $\sim 1$ for the cooler plasma emitting the O VII lines. We find that $f_{OVII}$ approaches unity at the same stellar surface X-ray flux level as solar active regions, suggesting that these stars become completely covered by active regions. At the same surface flux level, $f_{MgXI}$ is seen to increase more sharply with increasing surface flux. These results appear to support earlier suggestions that hot $10^7$ K plasma in active coronae arises from flaring activity, and that this flaring activity increases markedly once the stellar surface becomes covered with active regions., Comment: 53 pages, 19 figures, accepted for publication in Astrophysical Journal. A version of the paper with higher quality figures is available from http://www.astropa.unipa.it/Library/preprint.html