Search

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

NOTT (formerly Hi-5) is the L'-band (3.5-4.0~microns) nulling interferometer of Asgard, an instrument suite in preparation for the VLTI visitor focus. The primary scientific objectives of NOTT include characterizing (i) young planetary systems near the snow line, a critical region for giant planet formation, and (ii) nearby main-sequence stars close to the habitable zone, with a focus on detecting exozodiacal dust that could obscure Earth-like planets. In 2023-2024, the final warm optics have been procured and assembled in a new laboratory at KU Leuven. First fringes and null measurements were obtained using a Gallium Lanthanum Sulfide (GLS) photonic chip that was also tested at cryogenic temperatures. In this paper, we present an overall update of the NOTT project with a particular focus on the cold mechanical design, the first results in the laboratory with the final NOTT warm optics, and the ongoing Asgard integration activities. We also report on other ongoing activities such as the characterization of the photonic chip (GLS, LiNbO3, SiO), the development of the exoplanet science case, the design of the dispersion control module, and the progress with the self-calibration data reduction software., Comment: 11 pages (incl. 5 figures); Proc. SPIE Astronomical Telescopes + Instrumentation 2024 (Yokohama; Japan), Optical and Infrared Interferometry and Imaging VI

The Classical T Tauri Star (CTTS) stage is a critical phase of the star and planet formation process. In an effort to better understand the mass accretion process, which can dictate further stellar evolution and planet formation, a multi-epoch, multi-wavelength photometric and spectroscopic monitoring campaign of four CTTSs (TW Hya, RU Lup, BP Tau, and GM Aur) was carried out in 2021 and 2022/2023 as part of the Outflows and Disks Around Young Stars: Synergies for the Exploration of ULYSSES Spectra (ODYSSEUS) program. Here we focus on the HST UV spectra obtained by the HST Director's Discretionary Time UV Legacy Library of Young Stars as Essential Standards (ULLYSES) program. Using accretion shock modeling, we find that all targets exhibit accretion variability, varying from short increases in accretion rate by up to a factor of 3 within 48 hours, to longer decreases in accretion rate by a factor of 2.5 over the course of 1 year. This is despite the generally consistent accretion morphology within each target. Additionally, we test empirical relationships between accretion rate and UV luminosity and find stark differences, showing that these relationships should not be used to estimate the accretion rate for individual target. Our work reinforces that future multi-epoch and simultaneous multi-wavelength studies are critical in our understanding of the accretion process in low-mass star formation., Comment: 37 pages, 14 figures

Context. Gaia18cjb is one of the Gaia-alerted eruptive young star candidates which has been experiencing a slow and strong brightening during the last 13 years, similar to some FU Orionis-type objects. Aims. The aim of this work is to derive the young stellar nature of Gaia18cjb, determine its physical and accretion properties to classify its variability. Methods. We conducted monitoring observations using multi-filter optical and near-infrared photometry, as well as near-infrared spectroscopy. We present the analysis of pre-outburst and outburst optical and infrared light curves, color-magnitude diagrams in different bands, the detection of near-IR spectral lines, and estimates of both stellar and accretion parameters during the burst. Results. The optical light curve shows an unusually long (8 years) brightening event of 5 mag in the last 13 years, before reaching a plateau indicating that the burst is still on-going, suggesting a FUor-like nature. The same outburst is less strong in the infrared light curves. The near-infrared spectra, obtained during the outburst, exhibit emission lines typical of highly accreting low-intermediate mass young stars with typical EXor features. The spectral index of Gaia18cjb SED classifies it as a Class I in the pre-burst stage and a Flat Spectrum young stellar object (YSO) during the burst. Conclusions. Gaia18cjb is an eruptive YSO which shows FUor-like photometric features (in terms of brightening amplitude and length of the burst) and EXor-like spectroscopic features and accretion rate, as V350 Cep and V1647 Ori, classified as objects in between FUors and EXors, Comment: 21 pages, 13 figures

We show the first SPHERE/IRDIS and IFS data of the CO-rich debris disc around HD 131488. We use N-body simulations to model both the scattered light images and the SED of the disc in a self-consistent way. We apply the Henyey-Greenstein approximation, Mie theory, and the Discrete Dipole Approximation to model the emission of individual dust grains. Our study shows that only when gas drag is taken into account can we find a model that is consistent with scattered light as well as thermal emission data of the disc. The models suggest a gas surface density of $2\times10^{-5}\,M_\oplus/$au$^2$ which is in agreement with estimates from ALMA observations. Thus, our modelling procedure allows us to roughly constrain the expected amount of gas in a debris disc without actual gas measurements. We also show that the shallow size distribution of the dust leads to a significant contribution of large particles to the overall amount of scattered light. The scattering phase function indicates a dust porosity of $\sim0.2\ldots 0.6$ which is in agreement with a pebble pile scenario for planetesimal growth., Comment: 27 pages, 25 figures, 2 tables, accepted for publication in MNRAS

Gaia21bty, a pre-main sequence star that previously had shown aperiodic dips in its light curve, underwent a considerable $\Delta G\approx2.9$ mag brightening that occurred over a few months between 2020 October - 2021 February. The Gaia lightcurve shows that the star remained near maximum brightness for about $4-6$ months, and then started slowly fading over the next 2 years, with at least three superimposed $\sim$1 mag sudden rebrightening events. Whereas the amplitude and duration of the maximum is typical for EXors, optical and near-infrared spectra obtained at the maximum are dominated by features which are typical for FUors. Modelling of the accretion disc at the maximum indicates that the disc bolometric luminosity is 43 L$_{\odot}$ and the mass accretion rate is $2.5\times10^{-5}$ M$_{\odot}$ yr$^{-1}$, which are typical values for FUors even considering the large uncertainty in the distance ($1.7_{-0.4}^{+0.8}$ kpc). Further monitoring is necessary to understand the cause of the quick brightness decline, the rebrightening, and the other post-outburst light changes, as our multi-colour photometric data suggest that they could be caused by a long and discontinuous obscuration event. We speculate that the outburst might have induced large-scale inhomogeneous dust condensations in the line of sight leading to such phenomena, whilst the FUor outburst continues behind the opaque screen., Comment: Accepted to MNRAS

FU Orionis objects (FUors) are eruptive young stars, which exhibit outbursts that last from decades to a century. Due to the duration of their outbursts, and to the fact that only about two dozens of such sources are known, information on the end of their outbursts is limited. Here we analyse follow-up photometry and spectroscopy of Gaia21elv, a young stellar object, which had a several decades long outburst. It was reported as a Gaia science alert due to its recent fading by more than a magnitude. To study the fading of the source and look for signatures characteristic of FUors, we have obtained follow-up near infrared (NIR) spectra using Gemini South/IGRINS, and both optical and NIR spectra using VLT/X-SHOOTER. The spectra at both epochs show typical FUor signatures, such as a triangular shaped $H$-band continuum, absorption-line dominated spectrum, and P Cygni profiles. In addition to the typical FUor signatures, [OI], [FeII], and [SII] were detected, suggesting the presence of a jet or disk wind. Fitting the spectral energy distributions with an accretion disc model suggests a decrease of the accretion rate between the brightest and faintest states. The rapid fading of the source in 2021 was most likely dominated by an increase of circumstellar extinction. The spectroscopy presented here confirms that Gaia21elv is a classical FUor, the third such object discovered among the Gaia science alerts., Comment: Accepted to MNRAS

Absorption features in stellar atmospheres are often used to calibrate photocentric velocities for kinematic analysis of further spectral lines. The Li feature at $\sim$ 6708 {\AA} is commonly used, especially in the case of young stellar objects for which it is one of the strongest absorption lines. However, this is a complex line comprising two isotope fine-structure doublets. We empirically measure the wavelength of this Li feature in a sample of young stars from the PENELLOPE/VLT programme (using X-Shooter, UVES and ESPRESSO data) as well as HARPS data. For 51 targets, we fit 314 individual spectra using the STAR-MELT package, resulting in 241 accurately fitted Li features, given the automated goodness-of-fit threshold. We find the mean air wavelength to be 6707.856 {\AA}, with a standard error of 0.002 {\AA} (0.09 km/s) and a weighted standard deviation of 0.026 {\AA} (1.16 km/s). The observed spread in measured positions spans 0.145 {\AA}, or 6.5 km/s, which is up to a factor of six higher than typically reported velocity errors for high-resolution studies. We also find a correlation between the effective temperature of the star and the wavelength of the central absorption. We discuss how exclusively using this Li feature as a reference for photocentric velocity in young stars could potentially be introducing a systematic positive offset in wavelength to measurements of further spectral lines. If outflow tracing forbidden lines, such as [O i] 6300 {\AA}, are actually more blueshifted than previously thought, this then favours a disk wind as the origin for such emission in young stars., Comment: Accepted for publication in A&A

EX Lup is a low-mass pre-main sequence star that occasionally shows accretion-related outbursts. Here, we present JWST/MIRI medium resolution spectroscopy obtained for EX Lup fourteen years after its powerful outburst. EX Lup is now in quiescence and displays a Class II spectrum. We detect a forest of emission lines from molecules previously identified in infrared spectra of classical T Tauri disks: H2O, OH, H2, HCN, C2H2, and CO2. The detection of organic molecules demonstrates that they are back after disappearing during the large outburst. Spectral lines from water and OH are for the first time de-blended and will provide a much improved characterization of their distribution and density in the inner disk. The spectrum also shows broad emission bands from warm, sub-micron size amorphous silicate grains at 10 and 18 um. During the outburst, in 2008, crystalline forsterite grains were annealed in the inner disk within 1 au, but their spectral signatures in the 10 um silicate band later disappeared. With JWST we re-discovered these crystals via their 19.0, 20.0, and 23.5 um emission, whose strength implies that the particles are at ~3 au from the star. This suggests that crystalline grains formed in 2008 were transported outwards and now approach the water snowline, where they may be incorporated into planetesimals. Containing several key tracers of planetesimal and planet formation, EX Lup is an ideal laboratory to study the effects of variable luminosity on the planet-forming material and may provide explanation for the observed high crystalline fraction in solar system comets., Comment: 9 pages, 4 figures, accepted for publication in ApJL. JWST/MIRI spectrum is available for download at https://tinyurl.com/spexodisksJWST

FU Orionis-type objects (FUors) are low-mass pre-main-sequence objects which go through a short-lived phase (~100 years) of increased mass accretion rate (from 10^-8 to 10^-4 M_sun yr^-1). These eruptive young stars are in the early stages of stellar evolution and, thus, still deeply embedded in a massive envelope that feeds material to the circumstellar disk that is then accreted onto the star. Some FUors drive molecular outflows, i.e. low-velocity wide-angle magneto-hydrodynamical winds, that inject energy and momentum back to the surrounding envelopes, and help clear the material surrounding the young star. Here we present a 12CO (3--2), 13CO (3--2) and 12CO (4--3) survey of 20 FUor-type eruptive young stars observed with APEX. We use our 13CO (3--2) observations to measure the masses of the envelopes surrounding each FUor and find an agreement with the FUor evolutionary trend found from the 10um silicate feature. We find outflows in 11 FUors, calculate their masses and other kinematic properties, and compare these with those of outflows found around quiescent young stellar objects gathered from the literature. This comparison indicates that outflows in FUors are more massive than outflows in quiescent sources, and that FUor outflows have a higher ratio outflow mass with respect to the envelope than the quiescent sample, indicating that the eruptive young stars have lower star-forming efficiencies. Finally, we found that the outflow forces in FUors are similar to those of quiescent young stellar objects, indicating that their accretion histories are similar or that the FUor outflows have lower velocities., Comment: Accepted for publication in ApJ. 29 pages, 8 tables, 10 figures, 2 figuresets

OB associations, birthplaces of the most luminous stars, are key objects for understanding the formation of high-mass stars and their effects on their environments. The aim of this work is to explore the structure and kinematics of the Cepheus OB2 association and characterize the history of star formation in the region -- in particular, the role of the Cepheus Bubble, surrounding Cepheus OB2. Based on \gaia\ DR3 data we study the spatial and age distribution and kinematics of young stars in the region. We select candidate pre-main-sequence stars in the $M_\mathrm{G}$ vs. \gbp$-$\grp\ colour-magnitude diagram, and using a clustering algorithm, we identify 13 stellar groups belonging to Cep OB2. Four groups, consisting of 10-13 Myr old low- and intermediate-mass stars, are located in the interior of the bubble, and are part of the oldest subsystem of the association Cep OB2a. Younger groups are found on the periphery. The tangential velocities suggest that some groups on the periphery were born in an expanding system of star-forming clouds, whereas others have been formed due to the collision of their parent cloud with the expanding bubble., Comment: Accepted in MNRAS, 17 pages, 17 figures

T Tauri stars produce broad Lyman-alpha emission lines that contribute $\sim$88% of the total UV flux incident on the inner circumstellar disks. Lyman-alpha photons are generated at the accretion shocks and in the protostellar chromospheres and must travel through accretion flows, winds and jets, the protoplanetary disks, and the interstellar medium before reaching the observer. This trajectory produces asymmetric, double-peaked features that carry kinematic and opacity signatures of the disk environments. To understand the link between the evolution of Lyman-alpha emission lines and the disks themselves, we model HST-COS spectra from targets included in Data Release 3 of the Hubble UV Legacy Library of Young Stars as Essential Standards (ULLYSES) program. We find that resonant scattering in a simple spherical expanding shell is able to reproduce the high velocity emission line wings, providing estimates of the average velocities within the bulk intervening H I. The model velocities are significantly correlated with the K band veiling, indicating a turnover from Lyman-alpha profiles absorbed by outflowing winds to emission lines suppressed by accretion flows as the hot inner disk is depleted. Just 30% of targets in our sample have profiles with red-shifted absorption from accretion flows, many of which have resolved dust gaps. At this stage, Lyman-alpha photons may no longer intersect with disk winds along the path to the observer. Our results point to a significant evolution of Lyman-alpha irradiation within the gas disks over time, which may lead to chemical differences that are observable with ALMA and JWST., Comment: accepted to ApJ

Gaia DR3 contains 1.8 billion sources with G-band photometry, 1.5 billion of which with BP and RP photometry, complemented by positions on the sky, parallax, and proper motion. The median number of field-of-view transits in the three photometric bands is between 40 and 44 measurements per source and covers 34 months of data collection. We pursue a classification of Galactic and extra-galactic objects that are detected as variable by Gaia across the whole sky. Supervised machine learning (eXtreme Gradient Boosting and Random Forest) was employed to generate multi-class, binary, and meta-classifiers that classified variable objects with photometric time series in the G, BP, and RP bands. Classification results comprise 12.4 million sources (selected from a much larger set of potential variable objects) and include about 9 million variable stars classified into 22 variability types in the Milky Way and nearby galaxies such as the Magellanic Clouds and Andromeda, plus thousands of supernova explosions in distant galaxies, 1 million active galactic nuclei, and almost 2.5 million galaxies. The identification of galaxies was made possible by the artificial variability of extended objects as detected by Gaia, so they were published in the galaxy_candidates table of the Gaia DR3 archive, separate from the classifications of genuine variability (in the vari_classifier_result table). The latter contains 24 variability classes or class groups of periodic and non-periodic variables (pulsating, eclipsing, rotating, eruptive, cataclysmic, stochastic, and microlensing), with amplitudes from a few milli-magnitudes to several magnitudes., Comment: 106 pages, 513 figure panels, 4 tables. Accepted by Astronomy & Astrophysics (in press)

Classical T Tauri stars (CTTSs) are young stellar objects that accrete materials from their accretion disc influenced by their strong magnetic field. The magnetic pressure truncates the disc at a few stellar radii and forces the material to leave the disc plane and fall onto the stellar surface by following the magnetic field lines. However, this global scheme may be disturbed by the presence of a companion interacting gravitationally with the accreting component. This work is aiming to study the accretion and the magnetic field of the tight eccentric binary DQ Tau, composed of two equal-mass ($\sim$ 0.6 \msun ) CTTSs interacting at different orbital phases. We investigated the variability of the system using a high-resolution spectroscopic and spectropolarimetric monitoring performed with ESPaDOnS at the CFHT. We provide the first ever magnetic field analysis of this system, the Zeeman-Doppler imaging revealed a stronger magnetic field for the secondary than the primary (1.2 kG and 0.5 kG, respectively), but the small-scale fields analysed through Zeeman intensification yielded similar strengths (about 2.5 kG). The magnetic field topology and strengths are compatible with the accretion processes on CTTSs. Both components of this system are accreting, with a change of the main accretor during the orbital motion. In addition, the system displays a strong enhancement of the mass accretion rate at periastron and apastron. We also discovered, for the first time in this system, the apsidal motion of the orbital ellipse., Comment: 18 pages, 20 figures. Accepted for publication in MNRAS

Gaia19fct is one of the Gaia-alerted eruptive young stars that has undergone several brightening events. We conducted monitoring observations using multi-filter optical and near-infrared photometry, as well as near-infrared spectroscopy, to understand the physical properties of Gaia19fct and investigate whether it fits into the historically defined two classes. We present the analyses of light curves, color variations, spectral lines, and CO modeling. The light curves show at least five brightening events since 2015, and the multi-filter color evolutions are mostly gray. The gray evolution indicates that bursts are triggered by mechanisms other than extinction. Our near-infrared spectra exhibit both absorption and emission lines and show time-variability throughout our observations. We found lower rotational velocity and lower temperature from the near-infrared atomic absorption lines than from the optical lines, suggesting that Gaia19fct has a Keplerian rotating disk. The CO overtone features show a superposition of absorption and emission components, which is unlike other young stellar objects. We modeled the CO lines, and the result suggests that the emission and absorption components are formed in different regions. We found that although Gaia19fct exhibits characteristics of both types of eruptive young stars, FU Orionis-type objects (FUors) and EX Lupi-type objects (EXors), it shows more similarity with EXors in general., Comment: 29 pages, 14 figures, 7 tables, Accepted for publication in ApJ

Light curves of young star systems show photometric variability due to different kinematic, and physical processes. One of the main contributors to the photometric variability is the changing mass accretion rate, which regulates the interplay between the forming young star and the protoplanetary disk. We collected high-resolution spectroscopy in eight different epochs, as well as ground-based and space-borne multi-epoch optical and infrared photometry of WX Cha, an M0 binary system, with an almost edge-on disk (i = 87degrees) in the Chamaeleon I star-forming region. Spectroscopic observations cover 72 days, the ground-based optical monitoring covers 42 days while space-borne TESS photometry extends for 56 days. The multi-wavelength light curves exhibit quasi-periodic variability of 0.35 - 0.53 mag in the near-infrared, and of 1.3 mag in g band. We studied the variability of selected emission lines that trace the accretion, computed the accretion luminosity and the mass accretion rate using empirical relations and obtained values of the accretion luminosity between 1.6 and 3.2 Lsun and mass accretion rate between 3.31x10{-7} Msun/yr and 7.76x10^{-7} Msun/yr. Our results show that WX Cha is accreting at a rate larger than what is typical for T Tauri stars in the same star-forming region with the same stellar parameters. We theorize that this is due to the higher disk mass of WX Cha than what is usual for stars with similar stellar mass, and to the binary nature of the system. Daily changes in the accretion luminosity and in the extinction can explain the photometric variability.

DI Cha A is K0-type pre-main sequence star, the brightest component of a quadruple stellar system. Here we report on a detailed study of this star based on archival VLTI/MIDI and VLTI/PIONIER infrared interferometric observations, as well as optical--infrared photometric monitoring from ground-based and space-born instruments. We determined the structure of the circumstellar disk by fitting simultaneously the interferometric visibilities and the spectral energy distribution, using both analytical models and the radiative transfer code RADMC-3D. The modeling revealed that the radial density distribution of the disk appears to have a gap between 0.21 and 3.0 au. The inner ring, whose inner size coincides with the sublimation radius, is devoid of small, submicrometer-sized dust grains. The inner edge of the outer disk features a puffed-up rim, typically seen in intermediate-mass stars. Grain growth, although less progressed, was also detected in the outer disk. The inner ring is variable at mid-infrared wavelengths on both daily and annual timescales, while the star stays remarkably constant in the optical, pointing to geometrical or accretion changes in the disk as possible explanation for the flux variations., Comment: 19 pages, 18 figures

The coevolution of T Tauri stars and their surrounding protoplanetary disks dictates the timescales of planet formation. In this paper, we present magnetospheric accretion and inner disk wall model fits to NUV-NIR spectra of nine classical T Tauri stars in Orion OB1b as part of the Outflows and Disks around Young Stars: Synergies for the Exploration of ULLYSES Spectra (ODYSSEUS) Survey. Using NUV-optical spectra from the Hubble UV Legacy Library of Young Stars as Essential Standards (ULLYSES) Director's Discretionary Program and optical-NIR spectra from the PENELLOPE VLT Large Programme, we find that the accretion rates of these targets are relatively high for the region's intermediate age of 5.0 Myr; rates range from $0.5-17.2 \times 10^{-8}$ M$_{\odot}$/yr, with a median value of $1.2\times 10^{-8}$ M$_{\odot}$/yr. The NIR excesses can be fit with 1200-1800 K inner disk walls located at 0.05-0.10 AU from the host stars. We discuss the significance of the choice in extinction law, as the measured accretion rate depends strongly on the adopted extinction value. This analysis will be extended to the complete sample of T Tauri stars being observed through ULLYSES to characterize accretion and inner disks in star-forming regions of different ages and stellar populations., Comment: 22 pages, 3 figures, 8 tables. Accepted for publication in AJ

Young stars show variability on different time-scales from hours to decades, with a range of amplitudes. We studied two young stars, which triggered the Gaia Science Alerts system due to brightenings on a time-scale of a year. Gaia20bwa brightened by about half a magnitude, whereas Gaia20fgx brightened by about two and half magnitudes. We analyzed the Gaia light curves, additional photometry, and spectra taken with the Telescopio Nazionale Galileo and the Gran Telescopio Canarias. Several emission lines were detected toward Gaia20bwa, including hydrogen lines from H$\alpha$ to H$\delta$, Pa$\beta$, Br$\gamma$, and lines of Ca II, O I, and Na I. The H$\alpha$ and Br$\gamma$ lines were detected toward Gaia20fgx in emission in its bright state, with additional CO lines in absorption, and the Pa$\beta$ line with an inverse P Cygni profile during its fading. Based on the Br$\gamma$ lines the accretion rate was $(2.4-3.1)\times10^{-8}$ $M_\odot$ yr$^{-1}$ for Gaia20bwa and $(4.5-6.6)\times10^{-8}$ $M_\odot$ yr$^{-1}$ for Gaia20fgx during their bright state. The accretion rate of Gaia20fgx dropped by almost a factor of 10 on a time-scale of half a year. The accretion parameters of both stars were found to be similar to those of classical T Tauri stars, lower than those of young eruptive stars. However, the amplitude and time-scale of these brightenings place these stars to a region of the parameter space, which is rarely populated by young stars. This suggests a new class of young stars, which produce outbursts on a time-scale similar to young eruptive stars, but with smaller amplitudes., Comment: Accepted to MNRAS

Context. The Gaia third Data Release (DR3) presents the first catalogue of full-sky variable Young Stellar Object (YSO) candidates observed by the Gaia space telescope during the initial 34 months of science operations. Aims. Numerous types of variable stars were classified using photometric data collected by Gaia. One of the new classes presented in the Gaia DR3 is the class of YSOs showing brightness variability. We analysed 79 375 sources classified as YSO candidates to validate their young nature and investigate the completeness and purity of the sample. Methods. We cross-matched the Gaia DR3 YSO sample with numerous catalogues from the literature, including YSO catalogues based on optical and infrared data, as well as catalogues of extragalactic sources and Galactic variable stars. YSO catalogues were used to quantify the completeness of the Gaia DR3 YSO sample, while others were inspected to calculate the contamination. Results. Among the 79 375 potential Young Stellar Object candidates published in the Gaia DR3 variable star catalogue, the majority of these objects are distributed along the line of sight of well-known Star Forming Regions and the Galactic mid-plane. We found that the contamination level is well below 1%, while the completeness is also at the percent level, taking into account that the Gaia DR3 YSO sample is based on sources that showed significant variability during the data collection period. The number of sources in our sample that were never catalogued before as YSO candidates is in the order of $\sim$40 000 objects., Comment: 18 pages, 34 figures, 6 tables, accepted in A&A, Gaia DR3 paper

Classical T Tauri stars are low-mass objects, which are still accreting material from the surrounding circumstellar disk. The accretion process is essential in the formation of Sun-like stars and in setting the properties of the disk at the time when planet formation occurs. We constructed a complex dataset in order to examine the accretion process of VW Cha, a classical T Tauri multiple system with the aim of studying the physical origin of the photometric and spectroscopic variability of the system. The TESS Space Telescope observed VW Cha between 2019 April 22 and June 19, and we complemented these data with contemporaneous ground-based $I_CJHK$ band photometric measurements. In addition, we obtained high-resolution optical spectra with the VLT/ESPRESSO and the 2.2\,m/FEROS instruments. Analyzing these data, we found that the TESS light curve shows photometric variations on timescales from minutes to weeks with a peak-to-peak amplitude of $\sim$0.8 mag. The near-infrared light curves follow the shape of the optical measurements, however, the peak-to-peak amplitudes are slightly increasing towards the longer wavelengths. We took spectra in both fainter and brighter photometric states of the system, allowing us to examine the origin of a photometric brightening event. Our results show that this brightening event can be explained by increased accretion. In addition, our spectroscopic data also suggest that the primary component of VW Cha is a spectroscopic binary, as it was proposed in earlier works., Comment: Accepted for publication in ApJ. 18 pages, 10 figures

Close encounters between young stellar objects in star forming clusters are expected to dramatically perturb circumstellar disks. Such events are witnessed in numerical simulations of star formation, but few direct observations of ongoing encounters have been made. Here we report sub-0".1 resolution Atacama Large Millimeter Array (ALMA) and Jansky Very Large Array (JVLA) observations towards the million year old binary protostar Z CMa in dust continuum and molecular line emission. A point source ~4700 au from the binary has been discovered at both millimeter and centimeter wavelengths. It is located along the extension of a ~2000 au streamer structure previously found in scattered light imaging, whose counterpart in dust and gas emission is also newly identified. Comparison with simulations shows signposts of a rare flyby event in action. Z CMa is a "double burster", as both binary components undergo accretion outbursts, which may be facilitated by perturbations to the host disk by flybys., Comment: Published in Nature Astronomy. Here is the authors' version with the Supplementary Information integrated into the Methods section

Mass accretion from the circumstellar disk onto the protostar is a fundamental process during star formation. Measuring the mass accretion rate is particularly challenging for stars belonging to binary systems, because it is often difficult to discriminate which component is accreting. DQ~Tau is an almost equal-mass spectroscopic binary system where the components orbit each other every 15.8~days. The system is known to display pulsed accretion, i.e., the periodic modulation of the accretion by the components on eccentric orbit. We present multiepoch ESO/VLT X-Shooter observations of DQ~Tau, with the aim to determine which component of this system is the main accreting source. We use the absorption lines in the spectra to determine the radial velocity of the two components, and measure the continuum veiling as a function of wavelength and time. We fit the observed spectra with non-accreting templates to correct for the photospheric and chromospheric contribution. In the corrected spectra we study in details the profiles of the emission lines and calculate mass accretion rates for the system as a function of orbital phase. In accordance with previous findings, we detect elevated accretion close to periastron. We measure that the accretion rate varies between $10^{-8.5}$ and $10^{-7.3}$ Msun/yr. The emission line profiles suggest that both stars are actively accreting and the dominant accretor is not always the same component, varying in few orbits., Comment: 22 pages, 19 figures

We present follow-up photometric and spectroscopic observations, and subsequent analysis of Gaia20eae. This source triggered photometric alerts during 2020 after showing a $\sim$3 mag increase in its brightness. Its Gaia Alert light curve showed the shape of a typical eruptive young star. We carried out observations to confirm Gaia20eae as an eruptive young star and classify it. Its pre-outburst spectral energy distribution shows that Gaia20eae is a moderately embedded Class II object with $L_\mathrm{bol} = 7.22$ L$_\odot$. The color-color and color-magnitude diagrams indicate that the evolution in the light curve is mostly gray. Multiple epochs of the H$\alpha$ line profile suggest an evolution of the accretion rate and winds. The near-infrared spectra display several emission lines, a feature typical of EXor-type eruptive young stars. We estimated the mass accretion rate during the dimming phase to be $\dot{M} = 3-8 \times 10^{-7}$ M$_\odot$ yr$^{-1}$, higher than typical T Tauri stars of similar mass and comparable to other EXors. We conclude Gaia20eae is a new EXor-type candidate., Comment: Accepted for publication in ApJ. 20 pages, 11 figures

The nature of the gas in CO-rich debris discs remains poorly understood, as it could either be a remnant from the earlier Class II phase or of secondary origin, driven by the destruction of icy planetesimals. The aim of this paper was to elucidate the origin of the gas content in the debris discs via various simple molecules that are often detected in the less-evolved Class II discs. We present millimetre molecular line observations of nine circumstellar discs around A-type stars: four CO-rich debris discs (HD 21997, HD 121617, HD 131488, HD 131835) and five old Herbig Ae protoplanetary discs (HD 139614, HD 141569, HD 142666, HD 145718, HD 100453). The sources were observed with the Atacama Large Millimeter/submillimeter Array (ALMA) in Bands 5 and 6 with 1-2'' resolution. The Herbig Ae discs are detected in the CO isotopologues, CN, HCN, HCO$^+$, C2H, and CS lines. In contrast, only CO isotopologues are detected in the debris discs, showing a similar amount of CO to that found in the Herbig Ae protoplanetary discs. Using chemical and radiative transfer modelling, we show that the abundances of molecules other than CO in debris discs are expected to be very low. We consider multiple sets of initial elemental abundances with various degrees of H2 depletion. We find that the HCO+ lines should be the second brightest after the CO lines and that their intensities strongly depend on the overall CO/H$_2$ ratio of the gas. However, even in the ISM-like scenario, the simulated HCO+ emission remains weak as required by our non-detections., Comment: 15 pages, 9 figures, 7 tables Accepted for publication in MNRAS

Surprisingly strong CO emission has been observed from more than a dozen debris disks around nearby main-sequence stars. The origin of this CO is unclear, in particular whether it is left over from the protoplanetary disk phase or is second-generation material released from collisions between icy bodies like debris dust. The primary unexplored avenue for distinguishing the origin of the material is understanding its molecular composition. Here we present a deep search for five molecules (CN, HCN, HCO+, SiO, and CH3OH) in the debris disk around 49 Ceti. We take advantage of the high sensitivity of the Atacama Large Millimeter/submillimeter Array (ALMA) at Band 7 to integrate for 3.2 hours at modest spatial (1") and spectral (0.8 km/s) resolution. Our search yields stringent upper limits on the flux of all surveyed molecular lines, which imply abundances relative to CO that are orders of magnitude lower than those observed in protoplanetary disks and Solar System comets, and also those predicted in outgassing models of second-generation material. However, if CI shielding is responsible for extending the lifetime of any CO produced in second-generation collisions, as proposed by Kral et al. (2018), then the line ratios do not reflect true ice phase chemical abundances, but rather imply that CO is shielded by its own photodissociation product, CI, but other molecules are rapidly photodissociated by the stellar and interstellar radiation field., Comment: 26 pages, 4 figures, 4 tables, accepted for publication in ApJ

Millimetre continuum observations of debris discs can provide insights into the physical and dynamical properties of the unseen planetesimals that these discs host. The material properties and collisional models of planetesimals leave their signature on the grain size distribution, which can be traced through the millimetre spectral index. We present 8.8 mm observations of the debris discs HD 48370, CPD 72 2713, HD 131488, and HD 32297 using the Australian Telescope Compact Array (ATCA) as part of the PLanetesimals Around TYpicalPre-main seqUence Stars (PLATYPUS) survey. We detect all four targets with a characteristic beam size of 5 arcseconds and derive a grain size distribution parameter that is consistent with collisional cascade models and theoretical predictions for parent planetesimal bodies where binding is dominated by self-gravity. We combine our sample with 19 other millimetre-wavelength detected debris discs from the literature and calculate a weighted mean grain size power law index which is close to analytical predictions for a classical steady state collisional cascade model. We suggest the possibility of two distributions of q in our debris disc sample; a broad distribution (where q is approximately 3.2 to 3.7) for "typical" debris discs (gas-poor/non-detection), and a narrow distribution (where q is less than 3.2) for bright gas-rich discs. Or alternatively, we suggest that there exists an observational bias between the grain size distribution parameter and absolute flux which may be attributed to the detection rates of faint debris discs at cm wavelengths., Comment: 10 pages, 4 figures, accepted

Since July 2014, the Gaia mission has been engaged in a high-spatial-resolution, time-resolved, precise, accurate astrometric, and photometric survey of the entire sky. Aims: We present the Gaia Science Alerts project, which has been in operation since 1 June 2016. We describe the system which has been developed to enable the discovery and publication of transient photometric events as seen by Gaia. Methods: We outline the data handling, timings, and performances, and we describe the transient detection algorithms and filtering procedures needed to manage the high false alarm rate. We identify two classes of events: (1) sources which are new to Gaia and (2) Gaia sources which have undergone a significant brightening or fading. Validation of the Gaia transit astrometry and photometry was performed, followed by testing of the source environment to minimise contamination from Solar System objects, bright stars, and fainter near-neighbours. Results: We show that the Gaia Science Alerts project suffers from very low contamination, that is there are very few false-positives. We find that the external completeness for supernovae, $C_E=0.46$, is dominated by the Gaia scanning law and the requirement of detections from both fields-of-view. Where we have two or more scans the internal completeness is $C_I=0.79$ at 3 arcsec or larger from the centres of galaxies, but it drops closer in, especially within 1 arcsec. Conclusions: The per-transit photometry for Gaia transients is precise to 1 per cent at $G=13$, and 3 per cent at $G=19$. The per-transit astrometry is accurate to 55 milliarcseconds when compared to Gaia DR2. The Gaia Science Alerts project is one of the most homogeneous and productive transient surveys in operation, and it is the only survey which covers the whole sky at high spatial resolution (subarcsecond), including the Galactic plane and bulge., Comment: 26 pages, 26 figures, accepted for publication in Astronomy & Astrophysics

We present a new census of candidate pre-main-sequence stars in the Cepheus flare star-forming region, based on $\textit{Gaia}$ EDR3 parallaxes, proper motions, and colour$-$magnitude diagrams. We identified new candidate members of the previously known young stellar groups associated with NGC$\,$7023, L1177, L1217/L1219, L1228, L1235, and L1251. We studied the 3D structure of the star-forming complex and the distribution of tangential velocities of the young stars. The young stellar groups are located between 330 and 368$\,$pc from the Sun, divide into three kinematic subgroups, and have ages between 1$-$5 million years. The results confirm the scenario of propagating star formation, suggested by previous studies. In addition to the bulk pre-main-sequence star population between 330 and 370$~$pc, there is a scattered and more evolved pre-main-sequence population around 150$-$180$\,$pc. We found new candidate members of the nearby Cepheus Association, and identified a new moving group of 46, 15$-$20 million years old pre-main-sequence stars located at a distance of 178$\,$pc, around the A0-type star HD$\,$190833. A few pre-main-sequence stars are located at 800$-$900$\,$pc, indicative of star-forming regions associated with the Galactic local arm above the Galactic latitude of +10$^{\circ}$., Comment: Accepted in MNRAS, 20 pages, 16 figures

V555 Ori is a T Tauri star, whose 1.5 mag brightening was published as a Gaia science alert in 2017. We carried out optical and near-infrared photometric, and optical spectroscopic observations to understand the light variations. The light curves show that V555 Ori was faint before 2017, entered a high state for about a year, and returned to the faint state by mid-2018. In addition to the long-term flux evolution, quasi-periodic brightness oscillations were also evident, with a period of about 5 days. At optical wavelengths both the long-term and short-term variations exhibited colourless changes, while in the near-infrared they were consistent with changing extinction. We explain the brightness variations as the consequence of changing extinction. The object has a low accretion rate whose variation in itself would not be enough to reproduce the optical flux changes. This behaviour makes V555 Ori similar to the pre-main sequence star AA Tau, where the light changes are interpreted as periodic eclipses of the star by a rotating inner disc warp. The brightness maximum of V555 Ori was a moderately obscured ($A_V$=2.3 mag) state, while the extinction in the low state was $A_V$=6.4 mag. We found that while the Gaia alert hinted at an accretion burst, V555 Ori is a standard dipper, similar to the prototype AA Tau. However, unlike in AA Tau, the periodic behaviour was also detectable in the faint phase, implying that the inner disc warp remained stable in both the high and low states of the system., Comment: Accepted to MNRAS

To investigate the feasibility of ancillary target observations with ESA's Ariel mission, we compiled a list of potentially interesting young stars: FUors, systems harbouring extreme debris discs and a larger sample of young stellar objects showing strong near/mid-infrared excess. These objects can be observed as additional targets in the waiting times between the scheduled exoplanet transit and occultation observations. After analyzing the schedule for Ariel an algorithm was constructed to find the optimal target to be observed in each gap. The selection was mainly based on the slew and stabilization time needed to observe the selected YSO, but it also incorporated the scientific importance of the targets and whether they have already been sufficiently measured. After acquiring an adequately large sample of simulation data, it was concluded that approximately 99.2% of the available -- at least one hour long -- gaps could be used effectively. With an average slewing and stabilization time of about 16.7 minutes between scheduled exoplanet transits and ancillary targets, this corresponds to an additional $2881 \pm 56$ hours of active data gathering. When this additional time is used to observe our selected 200 ancillary targets, a typical signal-to-noise ratio of $\sim$10$^4$ can be achieved along the whole spectral window covered by Ariel., Comment: Accepted for publication in Experimental Astronomy

Extreme debris disks (EDDs) are rare systems with peculiarly large amounts of warm dust that may stem from recent giant impacts between planetary embryos during the final phases of terrestrial planet growth. Here we report on the identification and characterization of six new EDDs. These disks surround F5-G9 type main-sequence stars with ages >100 Myr, have dust temperatures higher than 300K and fractional luminosities between 0.01 and 0.07. Using time-domain photometric data at 3.4 and 4.6$\mu$m from the WISE all sky surveys, we conclude that four of these disks exhibited variable mid-infrared emission between 2010 and 2019. Analyzing the sample of all known EDDs, now expanded to 17 objects, we find that 14 of them showed changes at 3-5$\mu$m over the past decade suggesting that mid-infrared variability is an inherent characteristic of EDDs. We also report that wide-orbit pairs are significantly more common in EDD systems than in the normal stellar population. While current models of rocky planet formation predict that the majority of giant collisions occur in the first 100 Myr, we find that the sample of EDDs is dominated by systems older than this age. This raises the possibility that the era of giant impacts may be longer than we think, or that some other mechanism(s) can also produce EDDs. We examine a scenario where the observed warm dust stems from the disruption and/or collisions of comets delivered from an outer reservoir into the inner regions, and explore what role the wide companions could play in this process., Comment: 39 pages, 11 figures, accepted for publication in ApJ

While accreting through a circumstellar disk, young stellar objects are observed to undergo sudden and powerful accretion events known as FUor or EXor outbursts. Although such episodic accretion is considered to be an integral part of the star formation process, the triggers and mechanisms behind them remain uncertain. We conducted global numerical hydrodynamics simulations of protoplanetary disk formation and evolution in the thin-disk limit, assuming both magnetically layered and fully magnetorotational instability (MRI)-active disk structure. In this paper, we characterize the nature of the outbursts occurring in these simulations. The instability in the dead zone of a typical layered disk results in "MRI outbursts". We explore their progression and their dependence on the layered disk parameters as well as cloud core mass. The simulations of fully MRI-active disks showed an instability analogous to the classical thermal instability. This instability manifested at two temperatures--above approximately 1400 K and 3500 K--due to the steep dependence of Rosseland opacity on the temperature. The origin of these thermally unstable regions is related to the bump in opacity resulting from molecular absorption by water vapor and may be viewed as a novel mechanism behind some of the shorter duration accretion events. Although we demonstrated local thermal instability in the disk, more investigations are needed to confirm that a large-scale global instability will ensue. We conclude that the magnetic structure of a disk, its composition, as well as the stellar mass, can significantly affect the nature of episodic accretion in young stellar objects., Comment: 16 figures

We conducted global hydrodynamic simulations of protoplanetary disk evolution with an adaptive Shakura-Sunyaev {\alpha} prescription to represent the layered disk structure, and starting with the collapse phase of the molecular cloud. With the canonical values of model parameters, self-consistent dead zones formed at the scale of a few au. The instabilities associated with the dead zone and corresponding outbursts, similar to FUor eruptions, were also observed in the simulations.

Luminosity outbursts of the FUOri type stars, which have a magnitude of $\sim100 L_{\odot}$ and last for decades, may affect chemical composition of the surrounding protoplanetary disk. Using astrochemical modeling we analyze the changes induced by the outburst and search for species sensitive to the luminosity rise. Some changes in the disk molecular composition appear not only during the outburst itself but can also retain for decades after the end of the outburst. We analyze main chemical processes responsible for these effects and assess timescales at which chemically inert species return to the pre-outburst abundances., Comment: 2 pages, 2 figures

The Herbig Ae star HD 169142 is known to have a gaseous disk with a large inner hole, and also a photometrically variable inner dust component in the sub-au region. Following up our previous analysis, we further studied the temporal evolution of inner dust around HD 169142, which may provide information on the evolution from late-stage protoplanetary disks to debris disks. We used near-infrared interferometric observations obtained with VLTI/PIONIER to constrain the dust distribution at three epochs spanning six years. We also studied the photometric variability of HD 169142 using our optical-infrared observations and archival data. Our results indicate that a dust ring at ~0.3 au formed at some time between 2013 and 2018, and then faded (but did not completely disappear) by 2019. The short-term variability resembles that observed in extreme debris disks, and is likely related to short-lived dust of secondary origin, though variable shadowing from the inner ring could be an alternative interpretation. If confirmed, this is the first direct detection of secondary dust production inside a protoplanetary disk., Comment: 7 pages, 3 figures. Published in ApJ Letters

EX Lup is the prototype of a class of pre-main sequence eruptive stars defined by their repetitive outbursts lasting several months. In 2008 January-September EX Lup underwent its historically largest outburst, brightening by about 4 magnitudes in visual light. In previous studies we discovered on-going silicate crystal formation in the inner disk during the outburst, but also noticed that the measured crystallinity fraction started decreasing after the source returned to the quiescent phase. Here we present new observations of the 10 $\mu$m silicate feature, obtained with the MIDI and VISIR instruments at Paranal Observatory. The observations demonstrate that within five years practically all crystalline forsterite disappeared from the surface of the inner disk. We reconstruct this process by presenting a series of parametric axisymmetric radiative transfer models of an expanding dust cloud that transports the crystals from the terrestrial zone to outer disk regions where comets are supposed to form. Possibly the early Sun also experienced similar flare-ups, and the forming planetesimals might have incorporated crystalline silicate material produced by such outbursts. Finally, we discuss how far the location of the dust cloud could be constrained by future JWST observations., Comment: 12 pages, 4 figures, accepted for publication in ApJ

We present ALMA observations of 12CO, 13CO, and C18O J=2--1 lines and the 230 GHz continuum for the FU Ori-type object (FUor) V900 Mon (d~1.5 kpc), for which the accretion burst was triggered between 1953 and 2009. We identified CO emission associated with a molecular bipolar outflow extending up to a ~10^4 au scale and a rotating molecular envelope extending over >10^4 au. The interaction with the hot energetic FUor wind, which was observed using optical spectroscopy, appears limited to a region within ~400 au of the star. The envelope mass and the collimation of the extended CO outflow suggest that the progenitor of this FUor is a low-mass Class I young stellar object (YSO). These parameters for V900 Mon, another FUor, and a few FUor-like stars are consistent with the idea that FUor outbursts are associated with normal YSOs. The continuum emission is marginally resolved in our observations with a 0."2x0."15 (~300x225 au) beam, and a Gaussian model provides a deconvolved FWHM of ~90 au. The emission is presumably associated with a dusty circumstellar disk, plus a possible contribution from a wind or a wind cavity close to the star. The warm compact nature of the disk continuum emission could be explained with viscous heating of the disk, while gravitational fragmentation in the outer disk and/or a combination of grain growth and their inward drift may also contribute to its compact nature., Comment: 21 pages, 11 figures, accepted by ApJ

Previous observations revealed the existence of CO gas at nearly protoplanetary level in several dust-rich debris disks around young A-type stars. Here we used the ALMA 7m-array to measure $^{13}$CO and C$^{18}$O emission toward two debris disks, 49 Cet and HD 32297, and detected similarly high CO content ($>$0.01M$_\oplus$). These high CO masses imply a highly efficient shielding of CO molecules against stellar and interstellar ultraviolet photons. Adapting a recent secondary gas disk model that considers both shielding by carbon atoms and self-shielding of CO, we can explain the observed CO level in both systems. Based on the derived gas densities we suggest that, in the HD 32297 disk, dust and gas are coupled and the dynamics of small grains is affected by the gaseous component. For 49 Cet, the question of coupling remains undecided. We found that the main stellar and disk propertiesof 49 Cet and HD 32297 are very similar to those of previously identified debris disks with high CO content. These objects constitute together the first known representatives of shielded debris disks., Comment: submitted to The Astrophysical Journal

L1551 IRS 5 is a FUor-like object located in the Taurus star forming region. We present ALMA 1.3 mm continuum observations using a wide range of baselines. The observations recovered the two circumstellar disks composing the system and, for the first time, resolved the circumbinary ring. We determined the geometry and estimated lower mass limits for the circumstellar disks using simple models. We calculated lower limits for the total mass of both circumstellar disks. After subtracting the two circumstellar disk models from the image, the residuals show a clearly resolved circumbinary ring. Using a radiative transfer model, we show that geometrical effects can explain some of the brightness asymmetries found in the ring. The remaining features are interpreted as enhancements in the dust density., Comment: 8 pages, 3 figure, Accepted for publication in ApJL

Global numerical simulations of protoplanetary disk formation and evolution were conducted in thin-disk limit, where the model included magnetically layered disk structure, a self-consistent treatment for the infall from cloud core as well as the smallest possible inner computational boundary. We compared the evolution of a layered disk with a fully magnetically active disk. We also studied how the evolution depends on the parameters of the layered disk model - the MRI triggering temperature and active layer thickness - as well as the mass of the prestellar cloud core. With the canonical values of parameters a dead zone formed within the inner $\approx$ 15 au region of the magnetically layered disk. The dead zone was not a uniform structure and long-lived, axisymmetric, gaseous rings ubiquitously formed within this region due to the action of viscous torques. The rings showed a remarkable contrast in the disk environment as compared to a fully magnetically active disk and were characterized by high surface density and low effective viscosity. Multiple gaseous rings could form simultaneously in the dead zone region which were highly dynamical and showed complex, time-dependent behavior such as inward migration, vortices, gravitational instability and large-scale spiral waves. An increase in MRI triggering temperature had only marginal effects, while changes in active layer thickness as well as the initial cloud core mass had significant effects on the structure and evolution of the inner disk. Dust with large fragmentation barrier could be trapped in the rings, which may play a key role in planet formation., Comment: 23 pages, 13 figures

In a multi-wavelength study of thermal emission and scattered light images we analyse the dust properties and structure of the debris disc around the A1-type main sequence star 49~Cet. As a basis for this study, we present new scattered light images of the debris disc known to possess both a high amount of dust and gas. The outer region of the disc is revealed in former coronagraphic H-band and our new Y-band images from the Very Large Telescope SPHERE instrument. We use the knowledge of the disc's radial extent inferred from ALMA observations and the grain size distribution found by SED fitting to generate semi-dynamical dust models of the disc. We compare the models to scattered light and thermal emission data and find that a disc with a maximum of the surface density at 110~au and shallow edges can describe both thermal emission and scattered light observations. This suggests that grains close to the blow-out limit and large grains stem from the same planetesimal population and are mainly influenced by radiation pressure. The influence of inwards transport processes could not be analysed in this study., Comment: Accepted by MNRAS

Gas has been detected in a number of debris disks. It is likely secondary, i.e. produced by colliding solids. Here, we report ALMA Band 8 observations of neutral carbon in the CO-rich debris disk around the 15--30 Myr old A-type star HD 32297. We find that C$^0$ is located in a ring at $\sim$110 au with a FWHM of $\sim$80 au, and has a mass of $(3.5\pm0.2)\times10^{-3}$ M$_\oplus$. Naively, such a surprisingly small mass can be accumulated from CO photo-dissociation in a time as short as $\sim$10$^4$ yr. We develop a simple model for gas production and destruction in this system, properly accounting for CO self-shielding and shielding by neutral carbon, and introducing a removal mechanism for carbon gas. We find that the most likely scenario to explain both C$^0$ and CO observations, is one where the carbon gas is rapidly removed on a timescale of order a thousand years and the system maintains a very high CO production rate of $\sim$15 M$_\oplus$ Myr$^{-1}$, much higher than the rate of dust grind-down. We propose a possible scenario to meet these peculiar conditions: the capture of carbon onto dust grains, followed by rapid CO re-formation and re-release. In steady state, CO would continuously be recycled, producing a CO-rich gas ring that shows no appreciable spreading over time. This picture might be extended to explain other gas-rich debris disks., Comment: accepted for publication in the ApJ