Search

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

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

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., 15 pages, 9 figures, 7 tables Accepted for publication in MNRAS

Stars collect most of their mass during the protostellar stage, yet the accretion luminosity and stellar parameters, which are needed to compute the mass accretion rate, are poorly constrained for the youngest sources. The aim of this work is to fill this gap, computing the stellar properties and the accretion rates for a large sample of Class I protostars located in nearby (−8 and ∼10−4 M ⊙ yr−1 in a stellar mass range between ∼0.1 and 3 M ⊙. We find systematically larger mass accretion rates for our Class I sample than for Class II objects. Although the mass accretion rate we found is high, it still suggests that either stars collect most of their mass before the Class I stage, or eruptive accretion is needed during the overall protostellar phase. Indeed, our results suggest that for a large number of protostars the disk can be unstable, which can result in accretion bursts and disk fragmentation in the past or in the future.

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., 29 pages, 14 figures, 7 tables, Accepted for publication in ApJ

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

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

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., Accepted for publication in ApJ. 18 pages, 10 figures

Extreme debris discs are characterized by unusually strong mid-infrared excess emission, which often proves to be variable. The warm dust in these discs is of transient nature and is likely related to a recent giant collision occurring close to the star in the terrestrial region. Here we present the results of a 877 days long, gap-free photometric monitoring performed by the Spitzer Space Telescope of the recently discovered extreme debris disc around TYC 4209-1322-1. By combining these observations with other time-domain optical and mid-infrared data, we explore the disc variability of the last four decades with particular emphasis on the last 12 yr. During the latter interval the disc showed substantial changes, the most significant was the brightening and subsequent fading between 2014 and 2018 as outlined in WISE data. The Spitzer light curves outline the fading phase and a subsequent new brightening of the disc after 2018, revealing an additional flux modulation with a period of ~39 days on top of the long-term trend. We found that all these variations can be interpreted as the outcome of a giant collision that happened at an orbital radius of ~0.3 au sometime in 2014. Our analysis implies that a collision on a similar scale could have taken place around 2010, too. The fact that the disc was already peculiarly dust rich 40 yr ago, as implied by IRAS data, suggests that these dust production events belong to a chain of large impacts triggered by an earlier even more catastrophic collision., Comment: 18 pages, 12 figures, published in MNRAS

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., 19 pages, 18 figures

Exploring the properties of exoplanets near or inside the radius valley provides insights on the transition from the rocky super-Earths to the larger, hydrogen-rich atmosphere mini-Neptunes. Here, we report the discovery of TOI-1452 b, a transiting super-Earth ($R_{\rm p} = 1.67 \pm 0.07$ R$_{\oplus}$) in an 11.1--day temperate orbit ($T_{\rm eq} = 326 \pm 7$ K) around the primary member ($H = 10.0$, $T_{\rm eff} = 3185 \pm 50$ K) of a nearby visual binary M dwarf. The transits were first detected by TESS, then successfully isolated between the two $3.2^{\prime\prime}$ companions with ground-based photometry from OMM and MuSCAT3. The planetary nature of TOI-1452 b was established through high-precision velocimetry with the near-infrared SPIRou spectropolarimeter as part of the ongoing SPIRou Legacy Survey. The measured planetary mass ($4.8 \pm 1.3$ M$_{\oplus}$) and inferred bulk density ($5.6^{+1.8}_{-1.6}$ g/cm$^3$) is suggestive of a rocky core surrounded by a volatile-rich envelope. More quantitatively, the mass and radius of TOI-1452 b, combined with the stellar abundance of refractory elements (Fe, Mg and Si) measured by SPIRou, is consistent with a core mass fraction of $18\pm6$ % and a water mass fraction of $22^{+21}_{-13}$%. The water world candidate TOI-1452 b is a prime target for future atmospheric characterization with JWST, featuring a Transmission Spectroscopy Metric similar to other well-known temperate small planets such as LHS 1140 b and K2-18 b. The system is located near Webb's northern Continuous Viewing Zone, implying that is can be followed at almost any moment of the year., Published in The Astronomical Journal

Evidence of a relation between the mass accretion rate and the disk mass is established for young, Class II pre-main-sequence stars. This observational result opened an avenue to test theoretical models and constrain the initial conditions of disk formation, fundamental in the understanding of the emergence of planetary systems. However, it is becoming clear that planet formation starts even before the Class II stage, in disks around Class 0 and I protostars. We show for the first time evidence for a correlation between the mass accretion rate and the disk mass for a large sample of Class I young stars located in nearby (

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

We report results of a spectropolarimetric and photometric monitoring of the weak-line T Tauri star V410 Tau based on data collected mostly with SPIRou, the near-infrared (NIR) spectropolarimeter recently installed at the Canada-France-Hawaii Telescope, as part of the SPIRou Legacy Survey large programme, and with TESS between October and December 2019. Using Zeeman-Doppler Imaging (ZDI), we obtained the first maps of photospheric brightness and large-scale magnetic field at the surface of this young star derived from NIR spectropolarimetric data. For the first time, ZDI is also simultaneously applied to high-resolution spectropolarimetric data and very-high-precision photometry. V410 Tau hosts both dark and bright surface features and magnetic regions similar to those previously imaged with ZDI from optical data, except for the absence of a prominent dark polar spot. The brightness distribution is significantly less contrasted than its optical equivalent, as expected from the difference in wavelength. The large-scale magnetic field (~410 G), found to be mainly poloidal, features a dipole of ~390 G, again compatible with previous studies at optical wavelengths. NIR data yield a surface differential rotation slightly weaker than that estimated in the optical at previous epochs. Finally, we measured the radial velocity of the star and filtered out the stellar activity jitter using both ZDI and Gaussian Process Regression down to a precision of ~0.15 and 0.08 $\mathrm{km\,s^{-1}}$ RMS, respectively, confirming the previously published upper limit on the mass of a potential close-in massive planet around V410 Tau., 19 pages, 27 figures, accepted for publication in MNRAS

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 accretion of material from the protoplanetary disk onto the star is essential in the formation of Sun-like stars and in setting the properties of the disk at the time when planet formation happens. Although often described by simple and static models, the accretion process is inherently variable. For this reason, variability of young stars is perhaps the most important source of information on the accretion process and on the behavior of the innermost region of the systems. We designed a program to study six young stars in the Chamaeleon I star-forming region with the aim of characterizing the accretion process in time, to study the accretion diagnostic parameters and to examine the innermost part of the circumstellar disk. We combined the high-precision optical photometric observations of the TESS space telescope with contemporaneous ground-based IJHK-band photometry, and multi-epoch high-resolution optical spectroscopic observations obtained by the VLT/ESPRESSO and the 2.2m /FEROS spectrographs. Besides uncovering the periodic light variations that can be attributed to the stellar rotation, the uninterrupted TESS observations allowed us to examine the shorter timescale fluctuations due to accretion variability and also revealed flare-like events in the light curves. Our high-resolution spectroscopy allowed us to examine the amplitude, timescale and pattern of variability in accretion tracers which carry information on the distribution and the kinematics of the accreting material, the density structure of the inner disk, stellar activity, and the presence of outflows or jets. Where we have contemporaneous photometry and spectroscopy we linked the photometric variability to spectroscopic variations., {"references":["Hussain et al ., 2009, MNRAS, Vol 398, Issue 1, p 189 200","Carpenter et al ., 2001, AJ, Vol 121, Issue 6, p 3160 3190","Costigan et al ., 2014, MNRAS, Vol . 440, Issue 4, p 3444 3461"]}

FU Orionis-type objects (FUors) are low-mass pre-main sequence stars undergoing a temporary, but significant increase of mass accretion rate from the circumstellar disk onto the protostar. It is not yet clear what triggers the accretion bursts and whether the disks of FUors are in any way different from disks of non-bursting young stellar objects. Motivated by this, we conducted a 1.3 mm continuum survey of ten FUors and FUor-like objects with ALMA, using both the 7 m array and the 12 m array in two different configurations to recover emission at the widest possible range of spatial scales. We detected all targeted sources and several nearby objects as well. To constrain the disk structure, we fit the data with models of increasing complexity from 2D Gaussian to radiative transfer, enabling comparison with other samples modeled in a similar way. The radiative transfer modeling gives disk masses that are significantly larger than what is obtained from the measured millimeter fluxes assuming optically thin emission, suggesting that the FUor disks are optically thick at this wavelength. In comparison with samples of regular Class II and Class I objects, the disks of FUors are typically a factor of 2.9-4.4 more massive and a factor of 1.5-4.7 smaller in size. A significant fraction of them (65-70%) may be gravitationally unstable., Comment: 34 pages, 8 tables, 29 figures, accepted for publication in the Astrophysical Journal Supplement Series

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 multiepoch optical and infrared photometry of WX Cha, an M0 binary system, with an almost edge-on disk (i = 87°) 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 multiwavelength light curves exhibit quasi-periodic variability of 0.35–0.53 mag in the near-infrared, and of 1.3 mag in the 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 between L acc ∼ 1.6 L ⊙ − 3.2 L ⊙ and M ̇ acc ∼ 3.31 × 10 − 7 M ⊙ yr − 1 − 7.76 × 10 − 7 M ⊙ yr−1. 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.

V899 Mon is an eruptive young star showing characteristics of both FUors and EXors. It reached a peak brightness in 2010, then briefly faded in 2011, followed by a second outburst. We conducted multi-filter optical photometric monitoring, as well as optical and near-infrared spectroscopic observations of V899 Mon. The light curves and color-magnitude diagrams show that V899 Mon has been gradually fading after its second outburst peak in 2018, but smaller accretion bursts are still happening. Our spectroscopic observations taken with Gemini/IGRINS and VLT/MUSE show a number of emission lines, unlike during the outbursting stage. We used the emission line fluxes to estimate the accretion rate and found that it has significantly decreased compared to the outbursting stage. The mass loss rate is also weakening. Our 2D spectro-astrometric analysis of emission lines recovered jet and disk emission of V899 Mon. We found the emission from permitted metallic lines and the CO bandheads can be modeled well with a disk in Keplerian rotation, which also gives a tight constraint for the dynamical stellar mass of 2 ${M_{\odot}}$. After a discussion of the physical changes that led to the changes in the observed properties of V899 Mon, we suggest this object is finishing its second outburst., 31 pages, 26 figures, accepted for publication in ApJ

The second $Gaia$ Data Release (DR2) contains astrometric and photometric data for more than 1.6 billion objects with mean $Gaia$ $G$ magnitude $, Comment: 19 pages, 12 figures, 3 tables

The evolution of young stars and disks is driven by the interplay of several processes, notably accretion and ejection of material. Critical to correctly describe the conditions of planet formation, these processes are best probed spectroscopically. About five-hundred orbits of the Hubble Space Telescope (HST) are being devoted in 2020-2022 to the ULLYSES public survey of about 70 low-mass (M, Comment: Accepted for publication on Astronomy & Astrophysics. 15 pages + appendix, language edited version

We present results of our study of the close pre-main sequence spectroscopic binary CVSO 104 in Ori OB1, based on data obtained within the PENELLOPE legacy program. We derive, for the first time, the orbital elements of the system and the stellar parameters of the two components. The system is composed of two early M-type stars and has an orbital period of about 5 days and a mass ratio of 0.92, but contrarily to expectations does not appear to have a tertiary companion. Both components have been (quasi-)synchronized, but the orbit is still very eccentric. The spectral energy distribution clearly displays a significant infrared excess compatible with a circumbinary disk. The analysis of HeI and Balmer line profiles, after the removal of the composite photospheric spectrum, reveals that both components are accreting at a similar level. We also observe excess emission in H$\alpha$ and H$\beta$, which appears redshifted or blueshifted by more than 100 km/s with respect to the mass center of the system depending on the orbital phase. This additional emission could be connected with accretion structures, such as funnels of matter from the circumbinary disk. We also analyze the optical companion located at about 2".4 from the spectroscopic binary. This companion, that we named CVSO 104B, turns out to be a background Sun-like star not physically associated with the PMS system and not belonging to Ori OB1., Comment: Accepted for publication in Astronomy and Astrophysics; 16 pages, 18 figures

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

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$��$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$��$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., 39 pages, 11 figures, accepted for publication in ApJ

Full list of authors: Hodgkin, S. T.; Harrison, D. L.; Breedt, E.; Wevers, T.; Rixon, G.; Delgado, A.; Yoldas, A.; Kostrzewa-Rutkowska, Z.; Wyrzykowski, Ł.; van Leeuwen, M.; Blagorodnova, N.; Campbell, H.; Eappachen, D.; Fraser, M.; Ihanec, N.; Koposov, S. E.; Kruszyńska, K.; Marton, G.; Rybicki, K. A.; Brown, A. G. A.; Burgess, P. W.; Busso, G.; Cowell, S.; De Angeli, F.; Diener, C.; Evans, D. W.; Gilmore, G.; Holland, G.; Jonker, P. G.; van Leeuwen, F.; Mignard, F.; Osborne, P. J.; Portell, J.; Prusti, T.; Richards, P. J.; Riello, M.; Seabroke, G. M.; Walton, N. A.; Ábrahám, P.; Altavilla, G.; Baker, S. G.; Bastian, U.; O'Brien, P.; de Bruijne, J.; Butterley, T.; Carrasco, J. M.; Castañeda, J.; Clark, J. S.; Clementini, G.; Copperwheat, C. M.; Cropper, M.; Damljanovic, G.; Davidson, M.; Davis, C. J.; Dennefeld, M.; Dhillon, V. S.; Dolding, C.; Dominik, M.; Esquej, P.; Eyer, L.; Fabricius, C.; Fridman, M.; Froebrich, D.; Garralda, N.; Gomboc, A.; González-Vidal, J. J.; Guerra, R.; Hambly, N. C.; Hardy, L. K.; Holl, B.; Hourihane, A.; Japelj, J.; Kann, D. A.; Kiss, C.; Knigge, C.; Kolb, U.; Komossa, S.; Kóspál, Á.; Kovács, G.; Kun, M.; Leto, G.; Lewis, F.; Littlefair, S. P.; Mahabal, A. A.; Mundell, C. G.; Nagy, Z.; Padeletti, D.; Palaversa, L.; Pigulski, A.; Pretorius, M. L.; van Reeven, W.; Ribeiro, V. A. R. M.; Roelens, M.; Rowell, N.; Schartel, N.; Scholz, A.; Schwope, A.; Sipőcz, B. M.; Smartt, S. J.; Smith, M. D.; Serraller, I.; Steeghs, D.; Sullivan, M.; Szabados, L.; Szegedi-Elek, E.; Tisserand, P.; Tomasella, L.; van Velzen, S.; Whitelock, P. A.; Wilson, R. W.; Young, D. R.-- This is an Open Access article distributed under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited., Context. 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, CE = 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 CI = 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% at G = 13, and 3% at G = 19. The per-Transit astrometry is accurate to 55 mas 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. © S. T. Hodgkin et al. 2021., This work has made use of data from the European Space Agency (ESA) mission Gaia (https://www.cosmos.esa.int/gaia), processed by the Gaia Data Processing and Analysis Consortium (DPAC, https://www.cosmos.esa.int/web/gaia/dpac/consortium). Funding for the DPAC has been provided by national institutions, in particular the institutions participating in the Gaia Multilateral Agreement. Further details of funding authorities and individuals contributing to the success of the mission is shown at https://gea.esac.esa.int/archive/documentation/GEDR3/Miscellaneous/sec_acknowl/. We thank the United Kingdom Particle Physics and Astronomy Research Council (PPARC), the United Kingdom Science and Technology Facilities Council (STFC), and the United Kingdom Space Agency (UKSA) through the following grants to the University of Bristol, the University of Cambridge, the University of Edinburgh, the University of Leicester, the Mullard Space Sciences Laboratory of University College London, and the United Kingdom Rutherford Appleton Laboratory (RAL): PP/D006511/1, PP/D006546/1, PP/D006570/1, ST/I000852/1, ST/J005045/1, ST/K00056X/1, ST/K000209/1, ST/K000756/1, ST/L006561/1, ST/N000595/1, ST/S000623/1, ST/N000641/1, ST/N000978/1, ST/N001117/1, ST/S000089/1, ST/S000976/1, ST/S001123/1, ST/S001948/1, ST/S002103/1, and ST/V000969/1. This paper made use of the Whole Sky Database (WSDB) created by Sergey Koposov and maintained at the Institute of Astronomy, Cambridge with financial support from the Science and Technology Facilities Council (STFC) and the European Research Council (ERC). We thank the William Herschel and Isaac Newton Telescopes on the Roque de los Muchachos Observatory, La Palma, Spain, as well as the Optical Infrared Coordination Network for Astronomy (OPTICON) for their support of this project through telescope time, especially during the commissioning and verification phases. We thank the Copernico 1.82 m telescope (Mt. Ekar, Asiago Italy) operated by INAF Padova for supporting the project through telescope time (under the Large Programme Tomasella-SNe) during the verification phases. We acknowledge observations taken as part of the PESSTO project collected at the European Organisation for Astronomical Research in the Southern Hemisphere under ESO programme 199.D-0143. Authors at the ICCUB were supported by the Spanish Ministry of Science, Innovation and University (MICIU/FEDER, UE) through grant RTI2018-095076-B-C21, and the Institute of Cosmos Sciences University of Barcelona (ICCUB, Unidad de Excelencia ’María de Maeztu’) through grant CEX2019-000918-M. This work is supported by Polish NCN grants: Daina No. 2017/27/L/ST9/03221, Harmonia No. 2018/30/M/ST9/00311, Preludium No. 2017/25/N/ST9/01253 and MNiSW grant DIR/WK/2018/12 as well as the European Commission’s Horizon2020 OPTICON grant No. 730890. The Authors would like to thank the Warsaw University OGLE project for their continuous support in this work. AB acknowledges financial support from the Netherlands Research School for Astronomy (NOVA). AG acknowledges the financial support from the Slovenian Research Agency (grants P1-0031, I0-0033, J1-8136, J1-2460). AH was funded in part by the Leverhulme Trust through grant RPG-2012-541 and by the European Research Council grant 320360. AP acknowledges support from the NCN grant no. 2016/21/B/ST9/01126. CM acknowledges support from Jim and Hiroko Sherwin. DAK acknowledges support from the Spanish research projects AYA 2014-58381-P, AYA2017-89384-P, from Juan de la Cierva Incorporación fellowship IJCI-2015-26153, and from Spanish National Research Project RTI2018-098104-J-I00 (GRBPhot). EB and STH are funded by the Science and Technology Facilities Council grant ST/S000623/1. TW was funded in part by European Research Council grant 320360 and by European Commission grant 730980. GC acknowledges the Agenzia Spaziale Italiana (ASI) for its continuing support through contract 2018-24-HH.0 to the Italian Istituto Nazionale di Astrofisica (INAF). GD acknowledges the observing grant support from the Institute of Astronomy and Rozhen NAO BAS through the bilateral joint research project “Gaia Celestial Reference Frame (CRF) and fast variable astronomical objects” (during 2020-2022, leader is G. Damljanovic), and support by the Ministry of Education, Science and Technological Development of the Republic of Serbia (contract no. 451-03-68/2020-14/200002). G. Marton acknowledges support from the EC Horizon 2020 project OPTICON (730890) and the ESA PRODEX contract no. 4000129910. MF is supported by a Royal Society - Science Foundation Ireland University Research Fellowship NB acknowledges support from the research programme VENI, with project number 016.192.277, which is (partly) financed by the Netherlands Organisation for Scientific Research (NWO). NI is partially supported by Polish NCN DAINA grant no. 2017/27/L/ST9/03221. PAW acknowledges research funding from the South African National Research Foundation. RWW was funded by the Science and Technology Facilities Council grant ST/P000541/1. V.A.R.M.R. acknowledges financial support from Radboud Excellence Initiative, the Fundação para a Ciência e a Tecnologia (FCT) in the form of an exploratory project of reference IF/00498/2015/CP1302/CT0001, FCT and the Ministério da Ciência, Tecnologia e Ensino Superior (MCTES) through national funds and when applicable co-funded EU funds under the project UIDB/EEA/50008/2020, and supported by Enabling Green E-science for the Square Kilometre Array Research Infrastructure (ENGAGE-SKA), POCI-01-0145-FEDER-022217, and PHOBOS, POCI-01-0145-FEDER-029932, funded by Programa Operacional Competitividade e Internacionalização (COMPETE 2020) and FCT, Portugal. ZKR acknowledges funding from the Netherlands Research School for Astronomy (NOVA). ZN acknowledges support from the ESA PRODEX contract nr. 4000129910., With funding from the Spanish government through the Severo Ochoa Centre of Excellence accreditation SEV-2017-0709.

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., Accepted to MNRAS

Among the low-mass pre-main sequence stars, a small group called FU Orionis-type objects (FUors) are notable for undergoing powerful accretion outbursts. V1057 Cyg, a classical example of an FUor, went into outburst around 1969-1970, after which it faded rapidly, making it the fastest fading FUor known. Around 1995, a more rapid increase in fading occurred. Since that time, strong photometric modulations have been present. We present nearly 10 years of source monitoring at Piszk\'estet\H{o} Observatory, complemented with optical/near-infrared photometry and spectroscopy from the Nordic Optical Telescope, Bohyunsan Optical Astronomy Observatory, Transiting Exoplanet Survey Satellite, and the Stratospheric Observatory for Infrared Astronomy. Our light curves show continuation of significant quasi-periodic variability in brightness over the past decade. Our spectroscopic observations show strong wind features, shell features, and forbidden emission lines. All of these spectral lines vary with time. We also report the first detection of [S II], [N II], and [O III] lines in the star., Comment: 37 pages, 20 figures

A protoplanetary disk typically forms a dead zone near its midplane at the distance of a few au from the central protostar. Accretion through such a magnetically layered disk can be intrinsically unstable and has been associated with episodic outbursts in young stellar objects. We present the first investigation into the effects of low metallicity environment on the structure of the dead zone as well as the resulting outbursting behavior of the protoplanetary disk. We conducted global numerical hydrodynamic simulations of protoplanetary disk formation and evolution in the thin-disk limit. The consequences of metallicity were considered via its effects on the gas and dust opacity of the disk, the thickness of the magnetically active surface layer, and the temperature of the prestellar cloud core. We show that the metal poor disks accumulate much more mass in the innermost regions, as compared to the solar metallicity counterparts. The duration of the outbursting phase also varies with metallicity - the low metallicity disks showed more powerful luminosity eruptions with a shorter burst phase, which was confined mostly to the early, embedded stages of the disk evolution. The lowest metallicity disks with the higher cloud core temperature showed the most significant differences. The occurrence of outbursts was relatively rare in the disks around low mass stars and this was especially true at lowest metallicities. We conclude that the metal content of the disk environment can have profound effects on both the disk structure and evolution in terms of episodic accretion., Comment: 24 pages, 11 figures

Characterization of gas component in debris disks is of fundamental importance for understanding its origin. Toward this goal, we have conducted non-LTE (local thermodynamic equilibrium) analyses of the rotational spectral lines of CO including those of rare isotopologues ($^{13}$CO and C$^{18}$O) observed toward the gaseous debris disks of 49 Ceti and HD 21997 with the Atacama Large Millimeter/submillimeter Array (ALMA) and Atacama Compact Array (ACA). The analyses have been carried out for a wide range of the H$_{2}$ density, and the observed line intensities are found to be reproduced, as far as the H$_{2}$ density is higher than 10$^{3}$ cm$^{-3}$. The CO column density and the gas temperature are evaluated to be (1.8-5.9)$\times$10$^{17}$ cm$^{-2}$ and 8-11 K for 49 Ceti and (2.6-15)$\times$10$^{17}$ cm$^{-2}$ and 8-12 K for HD 21997, respectively, where the H$_{2}$ collision is assumed for the rotational excitation of CO. The results do not change significantly even if electron collision is considered. Thus, CO molecules can be excited under environments of no H$_{2}$ or a small number of H$_{2}$ molecules, even where the collision with CO, C, O, and C$^{+}$ would make an important contribution for the CO excitation in addition to H$_{2}$. Meanwhile, our result does not rule out the case of abundant H$_{2}$ molecules. The low gas temperature observed in the debris disks is discussed in terms of inefficient heating by interstellar and stellar UV radiation., 25 pages, 5 figures, 9 tables, accepted for publication in ApJ

Analyzing Spitzer and Herschel archival measurements we identified a debris disk around the young K7/M0 star CP-72 2713. The system belongs to the 24Myr old $\beta$ Pic moving group. Our new 1.33mm continuum observation, obtained with the ALMA 7-m array, revealed an extended dust disk with a peak radius of 140au, probably tracing the location of the planetesimal belt in the system. The disk is outstandingly large compared to known spatially resolved debris disks and also to protoplanetary disks around stars of comparable masses. The dynamical excitation of the belt at this radius is found to be reconcilable with planetary stirring, while self-stirring by large planetesimals embedded in the belt can work only if these bodies form very rapidly, e.g. via pebble concentration. By analyzing the spectral energy distribution we derived a characteristic dust temperature of 43K and a fractional luminosity of 1.1$\times$10$^{-3}$. The latter value is prominently high, we know only four other similarly dust-rich Kuiper-belt analogs within 40pc of the Sun., Comment: Published in AJ, 20 pages, 8 figures; v2:typos corrected, references revised, minor changes in the abstract and the introduction

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., accepted for publication in the ApJ

We conducted global hydrodynamic simulations of protoplanetary disk evolution with an adaptive Shakura-Sunyaev α 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., 2 pages, 2 figures

Full author list: E. Szegedi-Elek, P. Ábrahám, Ł. Wyrzykowski, M. Kun, Á. Kóspál, L. Chen, G. Marton, A. Moór, C. Kiss, A. Pál, L. Szabados, J. Varga, E. Varga-Verebélyi, C. Andreas, E. Bachelet, R. Bischoff, A. Bódi, E. Breedt, U. Burgaz, T. Butterley, J. M. Carrasco, V. Čepas, G. Damljanovic, I. Gezer, V. Godunova, M. Gromadzki, A. Gurgul, L. Hardy, F. Hildebrandt, S. Hoffmann, M. Hundertmark, N. Ihanec, R. Janulis, Cs. Kalup, Z. Kaczmarek, R. Könyves-Tóth, M. Krezinger, K. Kruszyńska, S. Littlefair, M. Maskoliūnas, L. Mészáros, P. Mikołajczyk, M. Mugrauer, H. Netzel, A. Ordasi, E. Pakštienė, K. A. Rybicki, K. Sárneczky, B. Seli, A. Simon, K. Šiškauskaitė, Á. Sódor, K. V. Sokolovsky2, W. Stenglein, R. Street, R. Szakáts, L. Tomasella, Y. Tsapras, K. Vida, J. Zdanavičius, M. Zieliński, P. Zieliński, and O. Ziółkowska, We present optical-infrared photometric and spectroscopic observations of Gaia 18dvy, located in the Cygnus OB3 association at a distance of 1.88 kpc. Gaia 18dvy was noted by the Gaia alerts system when its light curve exhibited a ⪆4 mag rise in 2018-2019. The brightening was also observable at mid-infared wavelengths. The infrared colors of Gaia 18dvy became bluer as the outburst progressed. Its optical and near-infrared spectroscopic characteristics in the outburst phase are consistent with those of bona fide FU Orionis-type young eruptive stars. The progenitor of the outburst is probably a low-mass K-type star with an optical extinction of ∼3 mag. A radiative transfer modeling of the circumstellar structure, based on the quiescent spectral energy distribution, indicates a disk with a mass of 4 10-3 M o˙. Our simple accretion disk modeling implies that the accretion rate had been exponentially increasing for more than 3 yr until mid-2019, when it reached a peak value of 6.9 10-6 M o˙ yr-1. In many respects, Gaia 18dvy is similar to the FU Ori-type object HBC 722., We thank the anonymous referee, whose questions and comments significantly improved the paper. This project has received funding from the European Research Council (ERC) under the European Union’s Horizon 2020 research and innovation program under grant agreement No. 716155 (SACCRED), Lendület LP2018-7/2019 and KEP-7/2018 of the Hungarian Academy of Sciences, GINOP 2.3.2-15-2016-00003 and PD-128360 of the Hungarian National Research, Development and Innovation Office, Polish NCN DAINA grant 2017/27/ L/ST9/03221, European Commission’s Horizon 2020 OPTICON grant 730890, Polish MNiSW grant DIR/WK/2018/12, grant No. S-LL-19-2 of the Research Council of Lithuania, Project No. 176011 “Dynamics and kinematics of celestial bodies and systems” of the Ministry of Education, Science and Technological Development of the Republic of Serbia, DFG priority program SPP 1992 “Exploring the Diversity of Extrasolar Planets” (WA 1047/11-1), the MINECO (Spanish Ministry of Economy) through grant RTI2018-095076-B-C21 (MINECO/FEDER, UE). The Joan Oró Telescope (TJO) of the Montsec Astronomical Observatory (OAdM) is owned by the Catalan Government and is operated by the Institute for Space Studies of Catalonia (IEEC). MG is supported by the Polish NCN MAESTRO grant 2014/14/ A/ST9/00121. We acknowledge ESA Gaia, DPAC, and the Photometric Science Alerts Team. We thank Christina Conner, Megan Davis, Alessandro Dellarovere, Hannah Gallamore, Mira Ghazali, Aaron Kruskie, Dylan Mankel, Jesse Leahy–McGregor, Brandon McIntyre, Barrett Ross, Courtney Wicklund, and Evan Zobel for observing Gaia 18dvy at the Michigan State University Observatory. Based on observations made with the Nordic Optical Telescope, operated by the Nordic Optical Telescope Scientific Association at the Observatorio del Roque de los Muchachos, La Palma, Spain, of the Instituto de Astrofísica de Canarias. Based on observations obtained with telescopes of the University Observatory Jena, which is operated by the Astrophysical Institute of the Friedrich-Schiller-University.

FU Orionis-type stars are young stellar objects showing large outbursts due to highly enhanced accretion from the circumstellar disk onto the protostar. FUor-type outbursts happen in a wide variety of sources from the very embedded ones to those with almost no sign of extended emission beyond the disk. The subsequent eruptions might gradually clear up the obscuring envelope material and drive the protostar on its way to become a disk-only T Tauri star. We used VLT/VISIR to obtain the first spectra that cover the 8-13 $\mu$m mid-infrared wavelength range in low-resolution of five recently discovered FUors. Four objects from our sample show the 10 $\mu$m silicate feature in emission. We study the shape and strength of the silicate feature in these objects and find that they mostly contain large amorphous grains, suggesting that large grains are typically not settled to the midplane in FUor disks. This is a general characteristic of FUors, as opposed to regular T Tauri-type stars whose disks display anything from pristine small grains to significant grain growth. We classify our targets by determining whether the silicate feature is in emission or in absorption, and confront them with the evolutionary scenarios on the dispersal of the envelopes around young stars. In our sample, all Class II objects exhibit silicate emission, while for Class I objects, the appearance of the feature in emission or absorption depends on the viewing angle with respect to the outflow cavity. This highlights the importance of geometric effects when interpreting the silicate feature., Comment: 7 pages, 1 table, 3 figures, accepted for publication in the Astrophysical Journal Letters