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