Search

The evolution of protoplanetary disks in regions with massive OB stars is influenced by externally driven winds that deplete the outer parts of disks. These winds have previously been studied via forbidden oxygen emission lines, which also arise in isolated disks in low-mass star forming-regions (SFRs) with weak external UV fields in photoevaporative or magnetic (internal) disk winds. It is crucial to determine how to disentangle external winds from internal ones. Here, we report a proxy for unambiguously identifying externally driven winds with a forbidden line of neutral atomic carbon, [C i] 8727 A. We compare for the first time the spatial location of the emission in the [O i] 5577 A, [O i] 6300 A, and [C i] 8727 A lines traced by VLT/MUSE-NFM, with the ALMA Band 7 continuum disk emission in a sample of 12 proplyds in the Orion Nebula Cluster (ONC). We confirm that the [O i] 5577 A emission is co-spatial with the disk emission, whereas the [O i] 6300 A is emitted both on the disk surface and on the ionization front of the proplyds. We show for the first time that the [C i] 8727 A line is also co-spatial with the disk surface in proplyds, as seen in the MUSE and ALMA data comparison. To verify whether the [C i] 8727 A line is detected in regions where external photoevaporation is not expected, we examine VLT/X-Shooter spectra for young stars in low-mass SFRs. Although the [O i] lines are well detected in all these targets, there is <<10% detection rate in the case of the [C i] 8727 A line. This number increases substantially to a ~40% detection rate in sigma-Orionis, a region with intermediate UV radiation. The spatial location of the [C i] 8727 A line emission and the lack of its detection in isolated disks in low-mass SFRs strongly suggest that this line is a tell-trace tracer of externally driven photoevaporative winds, which agrees with recent excitation models., Comment: Accepted to A&A on October 25, 2024

We present a study of the abundance of calcium in the innermost disk of 70 T Tauri stars in the star-forming regions of Chamaeleon I, Lupus and Orion OB1b. We use calcium as a proxy for the refractory material that reaches the inner disk. We used magnetospheric accretion models to analyze the Ca II emission lines and estimate abundances in the accretion flows of the stars, which feed from the inner disks. We find Ca depletion in disks of all three star-forming regions, with 57% of the sample having [Ca/H] < -0.30 relative to the solar abundance. All disks with cavities and/or substructures show depletion, consistent with trapping of refractories in pressure bumps. Significant Ca depletion ([Ca/H] < -0.30) is also measured in 60% of full disks, although some of those disks may have hidden substructures or cavities. We find no correlation between Ca abundance and stellar or disk parameters except for the mass accretion rate onto the star. This could suggest that the inner and outer disks are decoupled, and that the mass accretion rate is related to a mass reservoir in the inner disk, while refractory depletion reflects phenomena in the outer disk related to the presence of structure and forming planets. Our results of refractory depletion and timescales for depletion are qualitatively consistent with expectations of dust growth and radial drift including partitioning of elements and constitute direct evidence that radial drift of solids locked in pebbles takes place in disks., Comment: Accepted for publication in The Astrophysical Journal

Growing evidence suggests that protoplanetary discs may be influenced by late stage infall from the interstellar medium (ISM). It remains unclear the degree to which infall shapes disc populations at ages $\gtrsim 1$~Myr. We explore possible spatial correlations between stellar accretion rates in the Lupus star forming region, which would support the hypothesis that infall can regulate stellar accretion. We consider both the `clustered' stars towards the center of Lupus 3, and the `distributed' stars that are more sparsely distributed across the Lupus complex. We take the observed accretion rates in the literature and explore spatial correlations. In particular, we test whether the clustered stars exhibit a radial gradient in normalised accretion rates, and whether the distributed stars have spatially correlated accretion rates. We find statistically significant correlations for both the clustered and distributed samples. The clustered sample exhibits higher accretion rates in the central region, consistent with the expected Bondi-Hoyle-Lyttleton accretion rate. Stars that are spatially closer among the distributed population also exhibit more similar accretion rates. These results cannot be explained by the stellar mass distribution for either sample. Age gradients are disfavoured, though not discounted, because normalised disc dust masses are not spatially correlated across the region. Spatially correlated stellar accretion rates within the Lupus star forming region argue in favour of an environmental influence on stellar accretion, possibly combined with internal processes in the inner disc. Refined age measurements and searches for evidence of infalling material are potential ways to further test this finding., Comment: 10 pages, 6 figures, accepted for publication in A&A

A successful theory of star formation should predict the number of objects as a function of their mass produced through star-forming events. Previous studies in star-forming regions and the solar neighborhood identify a mass function increasing from the hydrogen-burning limit down to about 10 M$_{J}$. Theory predicts a limit to the fragmentation process, providing a natural turnover in the mass function down to the opacity limit of turbulent fragmentation thought to be 2-10 M$_{J}$. Programs to date have not been sensitive enough to probe the hypothesized opacity limit of fragmentation. Here we present the first identification of a turnover in the initial mass function below 12 M$_{J}$ within NGC 2024, a young star-forming region. With JWST/NIRCam deep exposures across 0.7-5 {\mu}m, we identified several free floating objects down to ~ 3 M$_{J}$ with sensitivity to 0.5 M$_{J}$. We present evidence for a double power law model increasing from about 60 M$_{J}$ to roughly 12 M$_{J}$, consistent with previous studies, followed by a decrease down to 0.5 M$_{J}$. Our results support the predictions of star and brown dwarf formation theory, identifying the theoretical turnover in the mass function and suggest the fundamental limit of turbulent fragmentation near 3 M$_{J}$., Comment: 16 pages, 4 figures, submitted 6 September 2024

The aim of this work is to characterize the accretion process of the classical T Tauri Star RU Lup. We studied optical high-resolution spectroscopic observations from CHIRON and ESPRESSO, obtained simultaneously with photometric data from AAVSO and TESS. We detected a periodic modulation in the narrow component of the He I 5876 line with a period that is compatible with the stellar rotation period, indicating the presence of a compact region on the stellar surface that we identified as the footprint of the accretion shock. We show that this region is responsible for the veiling spectrum, which is made up of a continuum component plus narrow line emission. An analysis of the high-cadence TESS light curve reveals quasi-periodic oscillations on timescales shorter than the stellar rotation period, suggesting that the accretion disk in RU~Lup extends inward of the corotation radius, with a truncation radius at $\sim 2 ~ R_{\star}$. This is compatible with predictions from three-dimensional magnetohydrodynamic models of accretion through a magnetic boundary layer (MBL). In this scenario, the photometric variability of RU Lup is produced by a nonstationary hot spot on the stellar surface that rotates with the Keplerian period at the truncation radius. The analysis of the broad components of selected emission lines reveals the existence of a non-axisymmetric, temperature-stratified flow around the star, in which the gas leaves the accretion disk at the truncation radius and accretes onto the star channeled by the magnetic field lines. The unusually rich metallic emission line spectrum of RU Lup might be characteristic of the MBL regime of accretion. In conclusion, the behavior of RU Lup reveals many similarities to predictions from the MBL accretion scenario. Alternative explanations would require the existence of a hot spot with a complex shape, or a warped structure in the inner disk., Comment: 20 pages, 21 figures. Accepted for publication in A&A

This paper analyses spectropolarimetric observations of the classical T Tauri star (CTTS) GM Aurigae collected with SPIRou, the near-infrared spectropolarimeter at the Canada-France-Hawaii Telescope, as part of the SLS and SPICE Large Programs. We report for the first time results on the large-scale magnetic field at the surface of GM Aur using Zeeman Doppler imaging. Its large-scale magnetic field energy is almost entirely stored in an axisymmetric poloidal field, which places GM Aur close to other CTTSs with similar internal structures. A dipole of about 730 G dominates the large-scale field topology, while higher-order harmonics account for less than 30 per-cent of the total magnetic energy. Overall, we find that the main difference between our three reconstructed maps (corresponding to sequential epochs) comes from the evolving tilt of the magnetic dipole, likely generated by non-stationary dynamo processes operating in this largely convective star rotating with a period of about 6 d. Finally, we report a 5.5$\sigma$ detection of a signal in the activity-filtered radial velocity data of semi-amplitude 110 $\pm$ 20 m/s at a period of 8.745 $\pm$ 0.009 d. If attributed to a close-in planet in the inner accretion disc of GM Aur, it would imply that this planet candidate has a minimum mass of 1.10 $\pm$ 0.30 Mjup and orbits at a distance of 0.082 $\pm$ 0.002 au., Comment: Published in MNRAS. 24 pages, 15 figures

Classical T Tauri Stars (CTTSs) are highly variable stars that possess gas- and dust-rich disks from which planets form. Much of their variability is driven by mass accretion from the surrounding disk, a process that is still not entirely understood. A multi-epoch optical spectral monitoring campaign of four CTTSs (TW Hya, RU Lup, BP Tau, and GM Aur) was conducted along with contemporaneous HST UV spectra and ground-based photometry in an effort to determine accretion characteristics and gauge variability in this sample. Using an accretion flow model, we find that the magnetospheric truncation radius varies between 2.5-5 R* across all of our observations. There is also significant variability in all emission lines studied, particularly Halpha, Hbeta, and Hgamma. Using previously established relationships between line luminosity and accretion, we find that, on average, most lines reproduce accretion rates consistent with accretion shock modeling of HST spectra to within 0.5 dex. Looking at individual contemporaneous observations, however, these relationships are less accurate, suggesting that variability trends differ from the trends of the population and that these empirical relationships should be used with caution in studies of variability., Comment: 39 pages, 10 figures, 12 tables

The dynamical interactions between young binaries can perturb the material distribution of their circumstellar disks, and modify the planet formation process. In order to constrain the impact and nature of the binary interaction in the RW Aur system (bound or unbound), we analyzed the circumstellar material at 1.3 mm wavelengths, as observed at multiple epochs by ALMA. We analyzed the disk properties through parametric visibility modeling, and we used this information to constrain the dust morphology and the binary orbital period. We imaged the dust continuum emission of RW Aur with a resolution of 3 au, and we find that the radius enclosing 90% of the flux (R90%) is 19 au and 14 au for RW Aur A and B, respectively. By modeling the relative distance of the disks at each epoch, we find a consistent trend of movement for the disk of RW Aur B moving away from the disk of RW Aur A at an approximate rate of 3 mas/yr (about 0.5 au/yr in sky-projected distance). By combining ALMA astrometry, historical astrometry, and the dynamical masses of each star, we constrain the RW Aur binary stars to be most likely in a high-eccentricity elliptical orbit with a clockwise prograde orientation relative to RW Aur A, although low-eccentricity hyperbolic orbits are not ruled out by the astrometry. Our analysis does not exclude the possibility of a disk collision during the last interaction, which occurred $295_{-74}^{+21}$ yr ago relative to beginning of 2024. Evidence for the close interaction is found in a tentative warp of 6 deg in the inner 3 au of the disk of RW Aur A, in the brightness temperature of both disks, and in the morphology of the gas emission. A narrow ring that peaks at 6 au around RW Aur B is suggestive of captured material from the disk around RW Aur A., Comment: Accepted in A&A. All self-calibrated data and products are publicly available in https://doi.org/10.5281/zenodo.12825068 . Short video summary in https://youtu.be/IOjEW8Vrj9Q

This study focuses on the validation of high-resolution regional reanalyses to understand their effectiveness in reproducing precipitation patterns over Italy, a climate change hotspot characterized by coastal sea-land interaction and complex orography. Nine reanalysis products were evaluated, with the ECMWF global reanalysis ERA5 serving as a benchmark. These included both European (COSMO-REA6, CERRA) and Italy-specific (BOLAM, MERIDA, MERIDA-HRES, MOLOCH, SPHERA, VHR-REA\_IT) datasets, using different models and parametrizations. The inter-comparison involved determining the effective resolution of daily precipitation fields using wavelet techniques and assessing intense precipitation statistics through frequency distributions. In-situ observations and observational gridded datasets were used to independently validate reanalysis precipitation fields. The capability of reanalyses to depict daily precipitation patterns was assessed, highlighting a maximum radius of precipitation misplacement of about 15 km, with notably lower skills during summer. An overall overestimation of precipitation was identified in the reanalysis climatological fields over the Po Valley and the Alps, whereas multiple products showed an underestimation of precipitations across the North-West coast, the Apennines, and Southern Italy. Finally, a comparison with a time-consistent observational dataset (UniMi/ISAC-CNR) revealed a non-stable deviation from observations in the annual precipitation cumulate of the reanalysis products analyzed. This should be taken into account when interpreting precipitation trends over Italy.

Studies of the stellar and accretion properties of classical T Tauri stars (CTTS) require comparison with photospheric spectral templates. Here we aim at expanding the currently available grid of wide-wavelength coverage observed spectra of non-accreting stars with additional new spectra and an interpolation method that allows us to obtain a continuous grid of low resolution spectra ranging from spectral type G8 to M9.5, while also mitigating observational uncertainties. This interpolated grid is then implemented in the self-consistent method to derive stellar and accretion properties of CTTS. With the new templates, we aim to estimate a lower limit on the accretion luminosities that can be obtained through a study of the UV excess emission using observed templates. We analyse the molecular photospheric features present in the VLT/X-Shooter spectra of the targets to perform a spectral classification, including estimates of their extinction. We apply a non-parametric fitting method to the full grid of observed templates to obtain an interpolated grid of templates. We use the uncertainties on our interpolated grid to estimate a lower limit on the accretion luminosity that we can measure with this method. We find that the measurable accretion luminosities ranges from $\sim 2.7$ dex lower than the stellar luminosity in M5.5 stars to $\sim 1.3$ dex lower for G8 stars. For young stars with masses of $\sim 1M_{\odot}$ and ages of 3-6 Myr this limit translates into an observational limit of mass accretion rate on the order of $10^{-10} \rm M_{\odot}/yr$. The implementation of an interpolated grid of observed templates allows us to better disentangle degenerate solutions, leading to a more reliable estimate of accretion rates in young accreting stars., Comment: Accepted to A&A, Version before changes by the language editor

We present the overview and first results from the North-PHASE Legacy Survey, which follows six young clusters for five years, using the 2 deg$^2$ FoV of the JAST80 telescope from the Javalambre Observatory (Spain). North-PHASE investigates stellar variability on timescales from days to years for thousands of young stars distributed over entire clusters. This allows us to find new YSO, characterise accretion and study inner disk evolution within the cluster context. Each region (Tr37, CepOB3, IC5070, IC348, NGC2264, and NGC1333) is observed in six filters (SDSS griz, u band, and J0660, which covers H$\alpha$), detecting cluster members as well as field variable stars. Tr37 is used to prove feasibility and optimise the variability analysis techniques. In Tr37, variability reveals 50 new YSO, most of them proper motion outliers. North-PHASE independently confirms the youth of astrometric members, efficiently distinguishes accreting and non-accreting stars, reveals the extent of the cluster populations along Tr37/IC1396 bright rims, and detects variability resulting from rotation, dips, and irregular bursts. The proper motion outliers unveil a more complex star formation history than inferred from Gaia alone, and variability highlights previously hidden proper motion deviations in the surrounding clouds. We also find that non-YSO variables identified by North-PHASE cover a different variability parameter space and include long-period variables, eclipsing binaries, RR Lyr, and $\delta$ Scuti stars. These early results also emphasize the power of variability to complete the picture of star formation where it is missed by astrometry., Comment: Accepted by MNRAS

The main goal of this study is to screen the PENELLOPE/UVES targets for outflow activity and find microjets via spectro-astrometry in, e.g., the [OI]$\lambda$6300 line. In total, 34 T\,Tauri stars of the PENELLOPE survey have been observed with the high resolution slit spectrograph UVES in three different slit positions rotated by $120^\text{o}$. Our spectro-astrometric analysis in the [OI]$\lambda$6300 wind line reveals two newly discovered microjets associated with Sz\,103 and XX\,Cha. Both microjets have an extent of about $0.04$ arcseconds, that is, $<10\,\text{au}$, and we confined their orientation by the three slit observations. Furthermore, we confirm the binary nature of VW\,Cha and CVSO\,109. We present (further) evidence that DK\,Tau\,B and CVSO\,104\,A are spectroscopic binaries. Sz\,115 is tentatively a spectroscopic binary. The origin of the LVC, that is, MHD winds versus photoevaporative winds, of the Sz\,103 and XX\,Cha microjets remains unclear., Comment: 12 pages, 8 figures, accepted for publication in A&A

A numerical investigation of the flow evolution over a pitching NACA 0012 airfoil incurring in deep dynamic stall phenomena is presented. The experimental data at Reynolds number Re = 135 000 and reduced frequency k = 0.1, provided by Lee and Gerontakos, are compared to numerical simulations using different turbulence models. After a preliminary space and time convergence study, two- and three-dimensional URANS with different turbulence models are explored, highlighting the advantages and the drawbacks. Then, the turbulence-resolving capabilities of hybrid RANS/LES strategies are exploited to recover and better represent the dynamic stall vortex. In detail, Scale-Adaptive Simulations (SAS) and Stress-Blended Eddy Simulations (SBES) are adopted. Furthermore, the LES resolved portion allows a spectral analysis of the force and moment coefficients to investigate the contribution of frequency lower than the pitching one. Finally, a comparison of the proposed approaches with other numerical simulations is given.

We present the first results of the JWST program PROJECT-J (PROtostellar JEts Cradle Tested with JWST ), designed to study the Class I source HH46 IRS and its outflow through NIRSpec and MIRI spectroscopy (1.66 to 28 micron). The data provide line-images (~ 6.6" in length with NIRSpec, and up to 20" with MIRI) revealing unprecedented details within the jet, the molecular outflow and the cavity. We detect, for the first time, the red-shifted jet within ~ 90 au from the source. Dozens of shock-excited forbidden lines are observed, including highly ionized species such as [Ne III] 15.5 micron, suggesting that the gas is excited by high velocity (> 80 km/s) shocks in a relatively high density medium. Images of H2 lines at different excitations outline a complex molecular flow, where a bright cavity, molecular shells, and a jet-driven bow-shock interact with and are shaped by the ambient conditions. Additional NIRCam 2 micron images resolve the HH46 IRS ~ 110 au binary system and suggest that the large asymmetries observed between the jet and the H2 wide angle emission could be due to two separate outflows being driven by the two sources. The spectra of the unresolved binary show deep ice bands and plenty of gaseous lines in absorption, likely originating in a cold envelope or disk. In conclusion, JWST has unraveled for the first time the origin of the HH46 IRS complex outflow demonstrating its capability to investigate embedded regions around young stars, which remain elusive even at near-IR wavelengths., Comment: 28 pages, 15 figures, Accepted for publication on The Astrophysical Journal (9 April 2024)

In the Orion Nebula Cluster (ONC), protoplanetary disks exhibit ionized gas clouds in the form of a striking teardrop shape as massive stars irradiate the disk material. We present the first spatially and spectrally resolved observations of 12 proplyds, using Integral Field Spectroscopy observations performed with the MUSE instrument in Narrow Field Mode (NFM) on the VLT. We present the morphology of the proplyds in seven emission lines and measure the radius of the ionization front (I-front) of the targets in four tracers, covering transitions of different ionization states for the same element. We also derive stellar masses for the targets. The measurements follow a consistent trend of increasing I-front radius for a decreasing strength of the far-UV radiation as expected from photoevaporation models. By analyzing the ratios of the I-front radii as measured in the emission lines of Ha, [OI] 6300, [OII] 7330, and [OIII] 5007, we observe the ionization stratification, that is, the most ionized part of the flow being the furthest from the disk (and closest to the UV source). The ratios of I-front radii scale in the same way for all proplyds in our sample regardless of the incident radiation. We show that the stratification can help constrain the densities near the I-front by using a 1D photoionization model. We derive the upper limits of photoevaporative mass-loss rates by assuming ionization equilibrium, and estimate values decreasing towards lower impinging radiation. We do not find a correlation between Mloss and stellar mass. The highest mass-loss rate is for the proplyd 244-440. These values of Mloss, combined with estimates of the disk mass with ALMA, confirm previous estimates of the short lifetime of these proplyds. This work demonstrates the potential of this dataset and offers a new set of observables to be used to test current and future models of external photoevaporation., Comment: Submitted to A&A on December 18, 2023. Accepted on April 2, 2024

Eruptive stars are a class of young stellar objects that show an abrupt increase in their luminosity. These burst-like episodes are thought to dominate the stellar accretion process during the class 0/class I stage. We present an overview of a survey of seven episodically accreting protostars aimed at studying their potentially complex circumstellar surroundings. The observations were performed with the instrument SPHERE, mounted at the VLT. We observed the eruptive stars in $H$-band with the near-infrared imager IRDIS and used the polarimeter to extract the polarized light scattered from the stars' surroundings. We produced polarized light images for three FUor objects, Z CMa, V960 Mon, and FU Ori, and four EXor objects, XZ Tau, UZ Tau, NY Ori, and EX Lup. We calculated the intrinsic polarization fraction for all the observed stars. In all systems we registered scattered light from around the primary star. FU Ori and V960 Mon are surrounded by complex structures including spiral-like features. In Z CMa, we detected a point source 0.7 arcsec to the northeast of the primary. Based on the astrometric measurements from archival Keck/NIRC2 data, we find this source to be a third member of the system. Further, Z CMa displays an outflow extending thousands of au. Unlike the other EXor objects in our sample, XZ Tau shows bright, extended scattered light structures, also associated with an outflow on a scale of hundreds of au. The other EXors show relatively faint disk-like structures in the immediate vicinity of the coronagraph. Asymmetric arms were only found around FUor objects, while faint disks seem to predominantly occur around EXors. Importantly, for Z CMa the detection of the faint extended structure questions previous interpretations of the system's dynamic state. The streamer which was associated with a fly-by object turned out to be part of a huge outflow extending 6000 au., Comment: Accepted for publication in A&A. 14 pages, 7 figures, 2 tables

We have performed a comprehensive study of the Orion Nebula Cluster (ONC) combining the photometric data obtained by the two \textit{HST} Treasury programs that targeted this region. To consistently analyze the rich dataset obtained in a wide variety of filters, we adopted a Bayesian approach to fit the Spectral Energy Distribution of the sources, deriving mass, age, extinction, distance, and accretion for each source in the region. The three dimensional study of mass distribution for bona-fide cluster members shows that mass segregation in the ONC extends to sub-solar masses, while the age distribution strongly supports the idea that star formation in the ONC is best described by a major episode of star formation that happened $\sim 1$ Myr ago. For masses $\gtrsim 0.1$ \Msun, our derived empirical initial mass function (IMF) is in good agreement with a Chabrier system IMF. Both the accretion luminosity (\Lacc) and mass accretion rates (\dMacc) are best described by broken power-law relations. This suggests that for the majority of young circumstellar disks in this cluster the excess emission may be dominated by X-ray-driven photoevaporation by the central star rather than external photoevaporation. If this is the case, the slopes of the power-law relations may be largely determined by the initial conditions set at the onset of the star formation process, which may be quite similar between regions that eventually form clusters of different sizes.

The sample of planet-forming disks observed by high-contrast imaging campaigns over the last decade is mature enough to enable the demographical analysis of individual star-forming regions. We present the full census of Taurus sources with VLT/SPHERE polarimetric images available. The whole sample sums up to 43 targets (of which 31 have not been previously published) corresponding to one-fifth of the Class II population in Taurus and about half of such objects that are observable. A large fraction of the sample is apparently made up of isolated faint disks (equally divided between small and large self-shadowed disks). Ambient signal is visible in about one-third of the sample. This probes the interaction with the environment and with companions or the outflow activity of the system. The central portion of the Taurus region almost exclusively hosts faint disks, while the periphery also hosts bright disks interacting with their surroundings. The few bright disks are found around apparently older stars. The overall picture is that the Taurus region is in an early evolutionary stage of planet formation. Yet, some objects are discussed individually, as in an intermediate or exceptional stage of the disk evolution. This census provides a first benchmark for the comparison of the disk populations in different star forming regions., Comment: Accepted by A&A

We present SPHERE/IRDIS H-band data for a sample of 23 stars in the Orion Star forming region observed within the DESTINYS (Disk Evolution Study Through Imaging of Nearby Young Stars) program. We use polarization differential imaging in order to detect scattered light from circumstellar dust. From the scattered light observations we characterize the disk orientation, radius and contrast. We analyse the disks in context of the stellar parameters and the environment of the Orion star-forming region. We use ancillary X-shooter spectroscopic observations to characterize the central stars in the systems. We furthermore use a combination of new and archival ALMA mm-continuum observations to characterize the dust masses present in the circumstellar disks. Within our sample we detect extended circumstellar disks in 10 of 23 systems. Of these, three are exceptionally extended (V351 Ori, V599 Ori and V1012 Ori) and show scattered light asymmetries which may indicate perturbations by embedded planets or (in the case of V599 Ori) by an outer stellar companion. Our high resolution imaging observations are also sensitive to close (sub)stellar companions and we detect 9 such objects in our sample of which 5 were previously unknown. We find in particular a possible sub-stellar companion (either a very low mass star or a high mass brown dwarf) 137 au from the star RY Ori. We find a strong anti-correlation between disk detection and multiplicity, with only 2 of our 10 disk detections located in stellar multiple systems. We also find a correlation between scattered light contrast and the millimetre flux suggesting that disks that have a high dust content are typically bright in near-infrared scattered light. Conversely we do not find significant correlations between scattered light contrast of the disks and the stellar mass or age., Comment: Accepted by A&A

We used VLT/SPHERE to observe 20 systems in the Cha I cloud in polarized scattered light in the near-infrared. We combined the scattered light observations with existing literature data on stellar properties and with archival ALMA continuum data to study trends with system age and dust mass. We also connected resolved near-infrared observations with the spectral energy distributions of the systems. In 13 of the 20 systems included in this study we detected resolved scattered light signals from circumstellar dust. For the CR Cha, CT Cha, CV Cha, SY Cha, SZ Cha, and VZ Cha systems we present the first detailed descriptions of the disks in scattered light. The observations found typically smooth or faint disks, often with little substructure, with the notable exceptions of SZ Cha, which shows an extended multiple-ringed disk, and WW Cha, which shows interaction with the cloud environment. New high S/N K- band observations of the HD 97048 system in our survey reveal a significant brightness asymmetry that may point to disk misalignment and subsequent shadowing of outer disk regions, possibly related to the suggested planet candidate in the disk. We resolve for the first time the stellar binary in the CS Cha system. Multiple wavelength observations of the disk around CS Cha have revealed that the system contains small, compact dust grains that may be strongly settled, consistent with numerical studies of circumbinary disks. We find in our sample that there is a strong anti-correlation between the presence of a (close) stellar companion and the detection of circumstellar material with five of our seven nondetections located in binary systems., Comment: Accepted by A&A

High-angular resolution imaging by ALMA has revealed the near-universality and diversity of substructures in protoplanetary disks. However, disks around M-type pre-main-sequence stars are still poorly sampled, despite the prevalence of M-dwarfs in the galaxy. Here we present high-resolution (~50 mas, 8 au) ALMA Band 6 observations of six disks around mid-M stars in Taurus. We detect dust continuum emission in all six disks, 12CO in five disks, and 13CO line in two disks. The size ratios between gas and dust disks range from 1.6 to 5.1. The ratio of about 5 for 2M0436 and 2M0450 indicates efficient dust radial drift. Four disks show rings and cavities and two disks are smooth. The cavity sizes occupy a wide range: 60 au for 2M0412, and ~10 au for 2M0434, 2M0436 and 2M0508. Detailed visibility modeling indicates that small cavities of 1.7 and 5.7 au may hide in the two smooth disks 2M0450 and CIDA 12. We perform radiative transfer fitting of the infrared SEDs to constrain the cavity sizes, finding that micron-sized dust grains may have smaller cavities than millimeter grains. Planet-disk interactions are the preferred explanation to produce the large 60 au cavity, while other physics could be responsible for the three ~10 au cavities under current observations and theories. Currently, disks around mid-to-late M stars in Taurus show a higher detection frequency of cavities than earlier type stars, although a more complete sample is needed to evaluate any dependence of substructure on stellar mass., Comment: 25 pages, 6 figures. Accepted for publication in ApJ

The time evolution of the dependence of the mass accretion rate with the stellar mass and the disk mass represents a fundamental way to understand the evolution of protoplanetary disks and the formation of planets. In this work, we present observations with X-Shooter of 26 Class II very low-mass stars and brown dwarfs in the Ophiuchus, Cha-I, and Upper Scorpius star-forming regions (SFRs). These new observations extend down to SpT M9 ($\sim$0.02 $M_\odot$) the measurement of the mass accretion rate in Ophiuchus and Cha-I and add 11 very-low-mass stars to the sample of objects studied with broadband spectroscopy in Upper Scorpius. We obtained their SpT, extinction and physical parameters, and we used the intensity of various emission lines to derive their accretion luminosity and mass accretion rates. Combining these new observations with data from the literature, we compare relations between accretion and stellar and disk properties of four different SFRs with different ages: Ophiuchus (1 Myr), Lupus (2 Myr), Cha-I (3 Myr), and Upper Scorpius (5-12 Myr). We find the slopes of the $L_*-L\mathrm{_{acc}}$ and $M_*-\dot{M}\mathrm{_{acc}}$ relationships to steepen between Ophiuchus, Lupus, and Cha-I and that both relationships may be better described with a single power law. We also find the relationship between the disk mass and the mass accretion rate of the stellar population to steepen with time down to the age of Upper Scorpius. Overall, we observe hints of a faster evolution into low accretion rates of low-mass stars and brown dwarfs. We also find that brown dwarfs present higher $M\mathrm{_{disk}}/\dot{M}\mathrm{_{acc}}$ ratios (i.e., longer accretion depletion timescales) than stars in Ophiuchus, Lupus, and Cha-I. This apparently contradictory result may imply that the evolution of protoplanetary disks around brown dwarfs is different from what is seen in the stellar regime., Comment: Accepted for publication at A&A, X-Shooter spectra will be made public on Vizier upon publication

This paper aims at revisit the physical and dynamical properties of the warm atomic gas in the inner disk region of classical T Tauri stars (CTTs) and relate them to the properties of the outer dusty disk. We used the high resolution (R=115,000) spectra of 36 CTTs observed as part of the GHOsT project and analysed the profile and luminosity of the brightest optical forbidden lines, namely [OI]630 and 557nm, [SII]406 and 673nm, and [NII]658nm. We find that in about 40% of sources the so-called narrow low-velocity component (NLVC) display a peak velocity compatible with the stellar velocity. In these sources, that typically show lower mass accretion rates and the absence of a jet, the [OI]630nm profiles are well fitted by a simple Keplerian disk model, indicating that the emission from the disk is dominant with respect to the wind contribution. For transitional disks (TD), no correlation is found between $R_{kep}$, derived from the line HWHM, and the size of the dust cavity. We also see an anti-correlation between the [OI] 557/630 nm ratio and $R_{kep}$, which suggests that the [OI] emitting region expands as the gas cools and becomes less dense. We confirmed previous findings on the density and temperature ranges implied by the line ratios, and additionally constrained the ionisation fraction in the NLVC to be < 0.1. We however discuss the limits of applying this diagnostic to winds that are not spatially resolved. For the outflow component, we estimated the mass-loss for both the disk winds and jets and compared the results with X-ray photoevaporative models. We conclude that without better knowledge of the wind geometry, and given the limitation of the diagnostics, the mass-loss in the wind traced by the LVC cannot be constrained better than a factor of 100, with a mass-loss/mass-accretion ratio spanning between ~ 0.01 and more than 1., Comment: Accepted for publication in Astronomy and Astrophysics

Observations of interstellar material infalling onto star- and planet-forming systems have become increasingly common thanks to recent advancements in radio interferometry. These structures replenish disks with fresh material, have the potential to significantly alter their dynamics, trigger the formation of substructures, induce shocks, and modify their physical and chemical properties. In this study, we combine new ALMA band 3 and archival band 6 observations to characterize the dust content and mass infall rate of a 4,000 au arc-like structure infalling onto M512, a class I young stellar object located in the Lynds 1641 region of the Orion A molecular cloud. We measure for the first time spectral index maps and derive a dust opacity index profile along a streamer, constraining grain properties and its dust mass. We measure a spectral index $\alpha \sim$ 3.2 across the entire structure, and a dust opacity index $\beta \sim$ 1.6. Given grain properties consistent with the measured $\beta$, the structure can host up to 245 M$_{\oplus}$ of dust, being comparable or even exceeding the mass of the inner, unresolved 600 au, which contains the protoplanetary disk of M512. Such a massive streamer can strongly affect the evolution of the star- and planet-forming inner system. Assuming typical ISM dust-to-gas ratio of 1%, free-fall timescales (50 kyr) imply total mass infall rates up to 1.5 $\cdot$ 10$^{-6}$ M$_{\odot}$/yr. M512 has been classified as an outbursting source with multi-epoch photometry, thus representing an interesting case study to explore the possible connection between infalling streamers and accretion outbursts., Comment: Accepted for publication in A&A, 13 pages, 5 figures, 3 tables

Protostellar jets and outflows are signposts of active star formation. In H II regions, molecular tracers like CO only reveal embedded portions of the outflow. Outside the natal cloud, outflows are dissociated, ionized, and eventually completely ablated, leaving behind only the high-density jet core. Before this process is complete, there should be a phase where the outflow is partially molecular and partially ionized. In this paper, we capture the HH 900 outflow while this process is in action. New observations from the ERIS/SPIFFIER near-IR integral field unit (IFU) spectrograph using the K-middle filter ($\lambda$=2.06-2.34 $\mu$m) reveal H$_2$ emission from the dissociating outflow and Br-$\gamma$ tracing its ionized skin. Both lines trace the wide-angle outflow morphology but H$_2$ only extends $\sim$5000 au into the H II region while Br-$\gamma$ extends the full length of the outflow ($\sim$12,650 au), indicating rapid dissociation of the molecules. H$_2$ has higher velocities further from the driving source, consistent with a jet-driven outflow. Diagnostic line ratios indicate that photoexcitation, not just shocks, contributes to the excitation in the outflow. We argue that HH 900 is the first clear example of an evaporating molecular outflow and predict that a large column of neutral material that may be detectable with ALMA accompanies the dissociating molecules. Results from this study will help guide the interpretation of near-IR images of externally irradiated jets and outflows such as those obtained with the James Webb Space Telescope (JWST) in high-mass star-forming regions where these conditions may be common., Comment: MNRAS, accepted

Individuals and organizations cope with an always-growing amount of data, which is heterogeneous in its contents and formats. An adequate data management process yielding data quality and control over its lifecycle is a prerequisite to getting value out of this data and minimizing inherent risks related to multiple usages. Common data governance frameworks rely on people, policies, and processes that fall short of the overwhelming complexity of data. Yet, harnessing this complexity is necessary to achieve high-quality standards. The latter will condition any downstream data usage outcome, including generative artificial intelligence trained on this data. In this paper, we report our concrete experience establishing a simple, cost-efficient framework that enables metadata-driven, agile and (semi-)automated data governance (i.e. Data Governance 4.0). We explain how we implement and use this framework to integrate 25 years of clinical study data at an enterprise scale in a fully productive environment. The framework encompasses both methodologies and technologies leveraging semantic web principles. We built a knowledge graph describing avatars of data assets in their business context, including governance principles. Multiple ontologies articulated by an enterprise upper ontology enable key governance actions such as FAIRification, lifecycle management, definition of roles and responsibilities, lineage across transformations and provenance from source systems. This metadata model is the keystone to data governance 4.0: a semi-automatised data management process that considers the business context in an agile manner to adapt governance constraints to each use case and dynamically tune it based on business changes.

While the most exciting explanation of the observed dust asymmetries in protoplanetary disks is the presence of protoplanets, other mechanisms can also form the dust features. This paper presents dual-wavelength Atacama Large Millimeter/submillimeter Array (ALMA) observations of a large asymmetric dusty ring around the M-type star CIDA 9A. We detect a dust asymmetry in both 1.3 mm and 3.1 mm data. To characterize the asymmetric structure, a parametric model is used to fit the observed visibilities. We report a tentative azimuthal shift of the dust emission peaks between the observations at the two wavelengths. This shift is consistent with a dust trap caused by a vortex, which may be formed by an embedded protoplanet or other hydrodynamical instabilities, such as a dead zone. Deep high-spatial observations of dust and molecular gas are needed to constrain the mechanisms that formed the observed millimeter cavity and dust asymmetry in the protoplanetary disk around CIDA 9A., Comment: Accepted for publication in ApJ, 12 pages, 7 figures with appendix

Quantum Sensing is a rapidly expanding research field that finds one of its applications in Fundamental Physics, as the search for Dark Matter. Devices based on superconducting qubits have already been successfully applied in detecting few-GHz single photons via Quantum Non-Demolition measurement (QND). This technique allows us to perform repeatable measurements, bringing remarkable sensitivity improvements and dark count rate suppression in experiments based on high-precision microwave photon detection, such as for Axions and Dark Photons search. In this context, the INFN Qub-IT project goal is to realize an itinerant single-photon counter based on superconducting qubits that will exploit QND for enhancing Axion search experiments. In this study, we present Qub-IT's status towards the realization of its first superconducting qubit device, illustrating design and simulation procedures and the characterization of fabricated Coplanar Waveguide Resonators (CPWs) for readout. We match target qubit parameters and assess a few-percent level agreement between lumped and distributed element simulation models. We reach a maximum internal quality factor of 9.2x10^5 for -92 dBm on-chip readout power.

Massive star-forming regions are thought to be the most common birth environments in the Galaxy and the only birth places of very massive stars. Their presence in the stellar cluster alters the conditions within the cluster impacting at the same time the evolution of other cluster members. In principle, copious amounts of ultraviolet radiation produced by massive stars can remove material from outer parts of the protoplanetary disks around low- and intermediate-mass stars in the process of external photoevaporation, effectively reducing the planet-formation capabilities of those disks. Here, we present deep VLT/MUSE observations of low-mass stars in Trumpler 14, one of the most massive, young, and compact clusters in the Carina Nebula Complex. We provide spectral and stellar properties of 717 sources and based on the distribution of stellar ages derive the cluster age of $\sim$1~Myr. The majority of the stars in our sample have masses $\leqslant$1~$M_\odot$, what makes our spectroscopic catalogue the most deep to date in term of masses, and proves that detailed investigations of low-mass stars are possible in the massive but distant regions. Spectroscopic studies of low-mass members of the whole Carina Nebula Complex are missing. Our work provides an important step forward towards filling this gap and set the stage for follow-up investigation of accretion properties in Trumpler 14., Comment: Accepted for publication in A&A, 27 pages, 28 figures

HM Lup is a young M-type star that accretes material from a circumstellar disk through a magnetosphere. Our aim is to study the inner disk structure of HM Lup and to characterize its variability. We used spectroscopic data from HST/STIS, X-Shooter, and ESPRESSO taken in the framework of the ULLYSES and PENELLOPE programs, together with photometric data from TESS and AAVSO. The 2021 TESS light curve shows variability typical for young stellar objects of the "accretion burster" type. The spectra cover the temporal evolution of the main burst in the 2021 TESS light curve. We compared the strength and morphology of emission lines from different species and ionization stages. We determined the mass accretion rate from selected emission lines and from the UV continuum excess emission at different epochs, and we examined its relation to the photometric light curves. The emission lines in the optical spectrum of HM Lup delineate a temperature stratification along the accretion flow. While the wings of the H I and He I lines originate near the star, the lines of species such as Na I, Mg I, Ca I, Ca II, Fe I, and Fe II are formed in an outer and colder region. The shape and periodicity of the 2019 and 2021 TESS light curves, when qualitatively compared to predictions from magnetohydrodynamic models, suggest that HM Lup was in a regime of unstable ordered accretion during the 2021 TESS observation due to an increase in the accretion rate. Although HM Lup is not an extreme accretor, it shows enhanced emission in the metallic species during this high accretion state that is produced by a density enhancement in the outer part of the accretion flow., Comment: 15 pages, 14 figures. Accepted for publication in A&A

The evolution of protoplanetary disks is regulated by an interplay of several processes, either internal to the system or related to the environment. As most of the stars and planets have formed in massive stellar clusters, studying the effects of UV radiation on disk evolution is of paramount importance. Here we test the impact of external photoevaporation on the evolution of disks in the $\sigma$ Orionis cluster by conducting the first combined large-scale UV to IR spectroscopic and mm-continuum survey of this region. We study a sample of 50 targets located at increasing distances from the central, OB system $\sigma$ Ori. We combine new VLT/X-Shooter spectra with new and previously published ALMA measurements of disk dust and gas fluxes and masses. We confirm the previously found decrease of $M_{\rm dust}$ in the inner $\sim$0.5 pc of the cluster. This is particularly evident when considering the disks around the more massive stars ($\ge$ 0.4 $M_{\odot}$), where those located in the inner part ($<$ 0.5 pc) have $M_{\rm dust}$ about an order of magnitude lower than the more distant ones. About half of the sample is located in the region of the $\dot{M}_{\rm acc}$ vs $M_{\rm disk}$ expected by models of external photoevaporation, namely showing shorter disk lifetimes. These are observed for all targets with projected separation from $\sigma$ Ori $<$ 0.5 pc, proving that the presence of a massive stellar system affects disk evolution. External photoevaporation is a viable mechanism to explain the observed shorter disk lifetimes and lower $M_{\rm dust}$ in the inner $\sim$0.5 pc of the cluster. Follow-up observations of the low stellar mass targets are crucial to confirm the dependence of the external photoevaporation process with stellar host mass. This work confirms that the effects of external photoevaporation are significant down to impinging radiation as low as $\sim 10^{4}$ G$_0$., Comment: Accepted for publication on Astronomy & Astrophysics. 13 pages, 7 figures + appendix. Abstract abridged to meet arXiv requirements

As the classic viscous paradigm for protoplanetary disk accretion is challenged by the observational evidence of low turbulence, the alternative scenario of MHD disk winds is being explored as potentially able to reproduce the same observed features traditionally explained with viscosity. Although the two models lead to different disk properties, none of them has been ruled out by observations - mainly due to instrumental limitations. In this work, we present a viable method to distinguish between the viscous and MHD framework based on the different evolution of the distribution in the disk mass ($M_{\mathrm{d}}$) - accretion rate ($\dot M$) plane of a disk population. With a synergy of analytical calculations and 1D numerical simulations, performed with the population synthesis code \texttt{Diskpop}, we find that both mechanisms predict the spread of the observed ratio $M_{\mathrm{d}}/\dot M$ in a disk population to decrease over time; however, this effect is much less pronounced in MHD-dominated populations as compared to purely viscous populations. Furthermore, we demonstrate that this difference is detectable with the current observational facilities: we show that convolving the intrinsic spread with the observational uncertainties does not affect our result, as the observed spread in the MHD case remains significantly larger than in the viscous scenario. While the most recent data available show a better agreement with the wind model, ongoing and future efforts to obtain direct gas mass measurements with ALMA and ngVLA will cause a reassessment of this comparison in the near future., Comment: 11 pages, 6 figures, 1 table. Published in ApJL

The Cryogenic IR echelle Spectrometer (CRIRES) instrument at the Very Large Telescope (VLT) was in operation from 2006 to 2014. Great strides in characterizing the inner regions of protoplanetary disks were made using CRIRES observations in the L- and M-band at this time. The upgraded instrument, CRIRES+, became available in 2021 and covers a larger wavelength range simultaneously. Here we present new CRIRES+ Science Verification data of the binary system S Coronae Australis (S CrA). We aim to characterize the upgraded CRIRES+ instrument for disk studies and provide new insight into the gas in the inner disk of the S CrA N and S systems. We analyze the CRIRES+ data taken in all available L- and M-band settings, providing spectral coverage from 2.9 to 5.5 $\mu$m. We detect emission from $^{12}$CO (v=1-0, v=2-1, and v=3-2), $^{13}$CO (v=1-0), hydrogen recombination lines, OH, and H$_2$O in the S CrA N disk. In the fainter S CrA S system, only the $^{12}$CO v=1-0 and the hydrogen recombination lines are detected. The $^{12}$CO v=1-0 emission in S CrA N and S shows two velocity components, a broad component coming from $\sim$0.1 au in S CrA N and $\sim$0.03 au in S CrA S and a narrow component coming from $\sim$3 au in S CrA N and $\sim$5 au in S CrA S. We fit local thermodynamic equilibrium slab models to the rotation diagrams of the two S CrA N velocity components and find that they have similar column densities ($\sim$1-7$\times$10$^{17}$ cm$^{-2}$), but that the broad component is coming from a hotter and narrower region. Two filter settings, M4211 and M4368, provide sufficient wavelength coverage for characterizing CO and H$_2$O at $\sim$5 $\mu$m, in particular covering low- and high-$J$ lines. CRIRES+ provides spectral coverage and resolution that are crucial complements to low-resolution observations, such as those with JWST, where multiple velocity components cannot be distinguished., Comment: 15 pages, 13 figures, accepted to A&A

Accretion plays a central role in the physics that governs the evolution and dispersal of protoplanetary disks. The primary goal of this paper is to analyze the stability over time of the mass accretion rate onto TW Hya, the nearest accreting solar-mass young star. We measure veiling across the optical spectrum in 1169 archival high-resolution spectra of TW Hya, obtained from 1998--2022. The veiling is then converted to accretion rate using 26 flux-calibrated spectra that cover the Balmer jump. The accretion rate measured from the excess continuum has an average of $2.51\times10^{-9}$~M$_\odot$~yr$^{-1}$ and a Gaussian distribution with a FWHM of 0.22 dex. This accretion rate may be underestimated by a factor of up to 1.5 because of uncertainty in the bolometric correction and another factor of 1.7 because of excluding the fraction of accretion energy that escapes in lines, especially Ly$\alpha$. The accretion luminosities are well correlated with He line luminosities but poorly correlated with H$\alpha$ and H$\beta$ luminosity. The accretion rate is always flickering over hours but on longer timescales has been stable over 25 years. This level of variability is consistent with previous measurements for most, but not all, accreting young stars., Comment: Accepted by ApJ. 31 pages

We present VLT/MUSE Narrow Field Mode (NFM) observations of a pair of disc-bearing young stellar objects towards the Orion Bar: 203-504 and 203-506. Both of these discs are subject to external photoevaporation, where winds are launched from their outer regions due to environmental irradiation. Intriguingly, despite having projected separation from one another of only 1.65{\arcsec} (660au at 400pc), 203-504 has a classic teardrop shaped ``proplyd'' morphology pointing towards $\theta^2$Ori A (indicating irradiation by the EUV of that star, rather than $\theta^1$ Ori C) but 203-506 has no ionisation front, indicating it is not irradiated by stellar EUV at all. However, 203-506 does show [CI] 8727{\AA} and [OI] 6300{\AA} in emission, indicating irradiation by stellar FUV. This explicitly demonstrates the importance of FUV irradiation in driving mass loss from discs. We conclude that shielding of 203-506 from EUV is most likely due to its position on the observers side of an ionized layer lying in the foreground of the Huygens Region. We demonstrate that the outflow HH 519, previously thought to be emanating from 203-504 is actually an irradiated cloud edge and identify a new compact outflow from that object approximately along our line of sight with a velocity $\sim130$\,km\,s$^{-1}$., Comment: 14 pages, 8 figures. Accepted for publication in MNRAS

PDS 70 is so far the only young disc where multiple planets have been detected by direct imaging. The disc has a large cavity when seen at sub-mm and NIR wavelengths, which hosts two massive planets. This makes PDS 70 the ideal target to study the physical conditions in a strongly depleted inner disc shaped by two giant planets, and in particular to test whether disc winds can play a significant role in its evolution. Using X-Shooter and HARPS spectra, we detected for the first time the wind-tracing [O I] 6300AA line, and confirm the low-moderate value of mass-accretion rate in the literature. The [O I] line luminosity is high with respect to the accretion luminosity when compared to a large sample of discs with cavities in nearby star-forming regions. The FWHM and blue-shifted peak of the [O I] line suggest an emission in a region very close to the star, favouring a magnetically driven wind as the origin. We also detect wind emission and high variability in the He I 10830AA line, which is unusual for low-accretors. We discuss that, although the cavity of PDS 70 was clearly carved out by the giant planets, the substantial inner disc wind could also have had a significant contribution to clearing the inner-disc., Comment: 13 pages, 7 figures, accepted for publication in ApJ

EX Lupi is the prototype by which EXor-type outbursts were defined. It has experienced multiple accretion-related bursts and outbursts throughout the last decades, whose study have greatly extended our knowledge about the effects of these types of events. This star experienced a new burst in 2022. We used multi-band photometry to create color-color and color-magnitude diagrams to exclude the possibility that the brightening could be explained by a decrease in extinction. We obtained VLT/X-shooter spectra to determine the Lacc and Macc during the peak of the burst and after its return to quiescence using 2 methods: empirical relationships between line luminosity and Lacc, and a slab model of the whole spectrum. We examined the 130 year light curve of EX Lupi to provide statistics on the number of outbursts experienced during this period of time. Our analysis of the data taken during the 2022 burst confirmed that a change in extinction is not responsible for the brightening. Our two approaches in calculating the Macc were in agreement, and resulted in values that are 2 orders of magnitude above what had previously been estimated, thus suggesting that EX Lupi is a strong accretor even when in quiescence. We determined that in 2022 March the Macc increased by a factor of 7 with respect to the quiescent level. We also found hints that even though the Macc had returned to almost its pre-outburst levels, certain physical properties of the gas had not returned to the quiescent values. We found that the mass accreted during this three month event was 0.8 lunar masses, which is approximately half of what is accreted during a year of quiescence. We calculated that if EX Lupi remains as active as it has been for the past 130 years, during which it has experienced at least 3 outbursts and 10 bursts, then it will deplete the mass of its circumstellar material in less than 160000 yr., Comment: 18 pages, 12 figures, 4 tables. Accepted for publication in A&A