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

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., Accepted for publication in ApJ. 29 pages, 8 tables, 10 figures, 2 figuresets

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

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

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

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

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

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

We present our results from a mid-infrared interferometric survey targeted at the planet-forming region in the circumstellar disks around low- and intermediate-mass young stars. Our sample consist of 82 objects, including T Tauri stars, Herbig Ae stars, and young eruptive stars. Our main results are: 1) Disks around T Tauri stars are similar to those around Herbig Ae stars, but are relatively more extended once we account for stellar luminosity. 2) From the distribution of the sizes of the mid-infrared emitting region we find that inner dusty disk holes may be present in roughly half of the sample. 3) Our analysis of the silicate spectral feature reveals that the dust in the inner $\sim$1 au region of disks is generally more processed than that in the outer regions. 4) The dust in the disks of T Tauri stars typically show weaker silicate emission in the N band spectrum, compared to Herbig Ae stars, which may indicate a general difference in the disk structure. Our data products are available at VizieR, and at the following web page: http://konkoly.hu/MIDI_atlas, 4 pages, 2 figures, published in "Origins: From the Protosun to the First Steps of Life. Proceedings of the International Astronomical Union, Volume 345"

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

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.

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

Extreme outbursts in young stars may be a common stage of pre-main-sequence stellar evolution. These outbursts, caused by enhanced accretion and accompanied by increased luminosity, can also strongly impact the evolution of the circumstellar environment. We present ALMA and VLA observations of EX Lupi, a prototypical outburst system, at 100 GHz, 45 GHz, and 15 GHz. We use these data, along with archival ALMA 232 GHz data, to fit radiative transfer models to EX Lupi's circumstellar disk in its quiescent state following the extreme outburst in 2008. The best fit models show a compact disk with a characteristic dust radius of 45 au and a total mass of 0.01 M$_{\odot}$. Our modeling suggests grain growth to sizes of at least 3 mm in the disk, possibly spurred by the recent outburst, and an ice line that has migrated inward to $0.2-0.3$ au post-outburst. At 15 GHz, we detected significant emission over the expected thermal disk emission which we attribute primarily to stellar (gyro)synchrotron and free-free disk emission. Altogether, these results highlight what may be a common impact of outbursts on the circumstellar dust., Accepted to ApJ, 15 pages, 8 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., 7 pages, 3 figures. Published in ApJ Letters

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., 21 pages, 11 figures, accepted by ApJ

V582 Aur is a pre-main sequence FU Orionis type eruptive star, which entered a brightness minimum in 2016 March due to changes in the line-of-sight extinction. Here, we present and analyze new optical $B$, $V$, $R_C$ and $I_C$ band multiepoch observations and new near-infrared $J$, $H$ and $K_S$ band photometric measurements from 2018 January$-$2019 February, as well as publicly available mid-infrared WISE data. We found that the source shows a significant optical$-$near-infrared variability, and the current brightness minimum has not completely finished yet. If the present dimming originates from the same orbiting dust clump that caused a similar brightness variation in 2012, than our results suggest a viscous spreading of the dust particles along the orbit. Another scenario is that the current minimum is caused by a dust structure, that is entering and leaving the inner part of the system. The WISE measurements could be consistent with this scenario. Our long-term data, as well as an accretion disk modeling hint at a general fading of V582 Aur, suggesting that the source will reach the quiescent level in $\sim$80 years., Comment: 8 pages, 4 figures, accepted for publication in ApJ

We studied the disc of the unclassified B[e] star HD 50138, in order to explore its structure, and to find indications for the evolutionary status of this system, whether it is a young Herbig Be or a post-main-sequence star. Using high spatial resolution interferometric measurements from MIDI instrument (N-band) on the Very Large Telescope Interferometer, we analysed the disc size, the time-variability of the disc's thermal emission, and the spectral shape of the $10~\mu$m silicate feature. By fitting simple disc models, we determined the inclination and the mid-infrared size of the disc, confirming earlier results based on a lower number of observations. We searched for mid-infrared temporal variability of different regions of the disc, and concluded that its morphology is not experiencing significant changes over the observed epochs. We characterized the mid-infrared silicate feature by determining the feature amplitude and the $11.3/9.8~\mu$m flux ratio. The latter parameter is a good indicator of the grain size. The shape of the feature suggests the presence of crystalline silicate grains in the disc. The interferometric data revealed a strong radial trend in the mineralogy: while the disc's innermost region seems to be dominated by forsterite grains, at intermediate radii both forsterite and enstatite may be present. The outer disc may predominantly contain amorphous silicate particles. A comparison of the observed spectral shape with that of a sample of intermediate-mass stars (supergiants, Herbig Ae/Be stars, unclassified B[e] stars) implied that the evolutionary state of HD 50138 cannot be unambiguously decided from mid-IR spectroscopy., Comment: 13 pages, 7 figures, to be published in MNRAS

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

We present near-IR imaging polarimetry of five classical FU Ori-type objects (FU Ori, V1057 Cyg, V1515 Cyg, V1735 Cyg, Z CMa) with a $\sim$0\farcs1 resolution observed using HiCIAO+AO188 at Subaru Telescope. We observed scattered light associated with circumstellar dust around four of them (i.e., all but V1515 Cyg). Their polarized intensity distribution shows a variety of morphologies with arms, tails or streams, spikes and fragmented distributions, many of which were reported in our previous paper. The morphologies of these reflection nebulae significantly differ from many other normal young stellar objects (Class I-II objects). These structures are attributed to gravitationally unstable disks, trails of clump ejections, dust blown by a wind or a jet, and a stellar companion. We can consistently explain our results with the scenario that their accretion outbursts (FUor outbursts) are triggered by gravitationally fragmenting disks, and with the hypothesis that many low-mass young stellar objects experience such outbursts., 21 pages, 10 figures, accepted by ApJ; resubmitted to astro-ph on August 22, 2018

An essential step to understanding protoplanetary evolution is the study of disks that contain gaps or inner holes. The pretransitional disk around the Herbig star HD 169142 exhibits multi-gap disk structure, differentiated gas and dust distribution, planet candidates, and near-infrared fading in the past decades, which make it a valuable target for a case study of disk evolution. Using near-infrared interferometric observations with VLTI/PIONIER, we aim to study the dust properties in the inner sub-au region of the disk in the years 2011-2013, when the object is already in its near-infrared faint state. We first performed simple geometric modeling to characterize the size and shape of the NIR-emitting region. We then performed Monte-Carlo radiative transfer simulations on grids of models and compared the model predictions with the interferometric and photometric observations. We find that the observations are consistent with optically thin gray dust lying at Rin ~ 0.07 au, passively heated to T ~ 1500 K. Models with sub-micron optically thin dust are excluded because such dust will be heated to much higher temperatures at similar distance. The observations can also be reproduced with a model consisting of optically thick dust at Rin ~ 0.06 au, but this model is plausible only if refractory dust species enduring ~2400 K exist in the inner disk., Comment: 13 pages, 9 figures

Context. Protoplanetary disks show large diversity regarding their morphology and dust composition. With mid-infrared interferometry the thermal emission of disks can be spatially resolved, and the distribution and properties of the dust within can be studied. Aims. Our aim is to perform a statistical analysis on a large sample of 82 disks around low- and intermediate-mass young stars, based on mid-infrared interferometric observations. We intend to study the distribution of disk sizes, variability, and the silicate dust mineralogy. Methods. Archival mid-infrared interferometric data from the MIDI instrument on the VLTI are homogeneously reduced and calibrated. Geometric disk models are used to fit the observations to get spatial information about the disks. An automatic spectral decomposition pipeline is applied to analyze the shape of the silicate feature. Results. We present the resulting data products in the form of an atlas, containing N band correlated and total spectra, visibilities, and differential phases. The majority of our data can be well fitted with a continuous disk model, except for a few objects, where a gapped model gives a better match. From the mid-infrared size--luminosity relation we find that disks around T Tauri stars are generally colder and more extended with respect to the stellar luminosity than disks around Herbig Ae stars. We find that in the innermost part of the disks ($r \lesssim 1$~au) the silicate feature is generally weaker than in the outer parts, suggesting that in the inner parts the dust is substantially more processed. We analyze stellar multiplicity and find that in two systems (AB Aur and HD 72106) data suggest a new companion or asymmetric inner disk structure. We make predictions for the observability of our objects with the upcoming MATISSE instrument, supporting the practical preparations of future MATISSE observations of T Tauri stars., Comment: 54 pages, 14 figures, 6 tables, accepted for publication in A&A

The FU Ori-type young stellar objects are characterized by a sudden increase in luminosity by 1$-$2 orders of magnitude, followed by slow return to the pre-outburst state on timescales of $\sim$10$-$100 yr. The outburst strongly affects the entire disk, changing its thermal structure and radiation field. In this paper, using a detailed physical-chemical model we study the impact of the FU Ori outburst on the disk chemical inventory. Our main goal is to identify gas-phase molecular tracers of the outburst activity that could be observed after the outburst with modern telescopes such as ALMA and NOEMA. We find that the majority of molecules experience a considerable increase in the total disk gas-phase abundances due to the outburst, mainly due to the sublimation of their ices. Their return to the pre-outburst chemical state takes different amounts of time, from nearly instantaneous to very long. Among the former ones we identify CO, NH$_3$, C$_2$H$_6$, C$_3$H$_4$, etc. Their abundance evolution tightly follows the luminosity curve. For CO the abundance increase does not exceed an order of magnitude, while for other tracers the abundances increase by 2$-$5 orders of magnitude. Other molecules like H$_2$CO and NH$_2$OH have longer retention timescales, remaining in the gas phase for $\sim10-10^3$ yr after the end of the outburst. Thus H$_2$CO could be used as an indicator of the previous outbursts in the post-outburst FU Ori systems. We investigate the corresponding time-dependent chemistry in detail and present the most favorable transitions and ALMA configurations for future observations., Comment: accepted for publication in ApJ

Having disk-to-star accretion rates on the order of 10$^{-4}\,M_{\odot}$/yr, FU Orionis-type stars (FUors) are thought to be the visible examples for episodic accretion. FUors are often surrounded by massive envelopes, which replenish the disk material and enable the disk to produce accretion outbursts. We observed the FUor-type star V346 Nor with ALMA at 1.3~mm continuum and in different CO rotational lines. We mapped the density and velocity structure of its envelope and analyzed the results using channel maps, position-velocity diagrams, and spectro-astrometric methods. We discovered a pseudo-disk and a Keplerian disk around a 0.1 $M_{\odot}$ central star. We determined an infall rate from the envelope onto the disk of 6$\times$10$^{-6}\,M_{\odot}$/yr, a factor of few higher than the quiescent accretion rate from the disk onto the star. This hints for a mismatch between the infall and accretion rates as the cause of the eruption., Comment: 2 pages, 1 figure, published in "Origins: From the Protosun to the First Steps of Life. Proceedings of the International Astronomical Union, Volume 345"

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

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 the last decade many kinematic groups of young stars (, Comment: Accepted for publication in MNRAS, 14 pages, 3 figures, 3 tables

Context: Quantifying the gas content inside the dust gaps of transition disks is important to establish their origin. Aims: We seek to constrain the surface density of warm gas in the disk of HD 139614, a Herbig Ae star with a transition disk exhibiting a dust gap from 2.3 to 6 AU. Methods: We have obtained ESO/VLT CRIRES high-resolution spectra of CO ro-vibrational emission. We derived constraints on the disk's structure by modeling the line-profiles, the spectroastrometric signal, and the rotational diagrams using flat Keplerian disk models. Results: We detected v=1-0 12CO, 2-1 12CO, 1-0 13CO, 1-0 C18O, and 1-0 C17O ro-vibrational lines. 12CO v=1-0 lines have an average width of 14 km/s, Tgas of 450 K and an emitting region from 1 to 15 AU. 13CO and C18O lines are on average 70 and 100 K colder, 1 and 4 km/s narrower, and are dominated by emission at R>6 AU. The 12CO v=1-0 line-profile indicates that if there is a gap in the gas it must be narrower than 2 AU. We find that a drop in the gas surface density (delta_gas) at R2 AU) gas gap suggest the presence of an embedded, Comment: accepted for publication in A&A; 30 pages (including Appendix), 24 figure; v2-v4: text corrections

Z CMa is a complex pre-main sequence binary with a current separation of 100 mas, known to consist of an FU Orionis star (SE component) and an embedded Herbig Be star (NW component). Immediately when the late-2008 outburst of Z CMa was announced to the community, we initiated a high angular resolution imaging campaign with VLT/NaCo, Keck/NIRC2, VLT/SINFONI, and Keck/OSIRIS which aimed at characterizing the outburst of both components of the system in the near-infrared. We confirm that the NW star dominates the system flux in the 1.1-3.8 microns range and is responsible for the photometric outburst. We extract the first medium-resolution (R=2000-4000) near-infrared (1.1-2.4 microns) spectra of the individual components during and after the outburst. The SE component has a spectrum typical of FU Orionis objects. The NW component spectrum is characteristic of embedded outbursting protostars and EX Or objects. It displays numerous emission lines during the outburst whose intensity correlates with the system activity. In particular, we find a correlation between the Brackett gamma equivalent width and the system brightness. The bluing of the continuum of the NW component along with the absolute flux and color-variation of the system during the outburst suggests that the outburst was caused by a complex interplay between a variation of the extinction in the line of sight of the NW component on one hand, and the emission of shocked regions close to the NW component on the other. We confirm the recently reported wiggling of the SE component jet from [Fe II] line emission. We find a point-like structure associated with a peak emission at 2.098 microns coincidental with the clump or arm seen in broadband polarization differential imaging as well as additional diffuse emission along a PA=214 degrees. The origin of these two structures is unclear and deserves further investigation., Comment: 18 pages, 18 figures, 7 tables. Accepted for publication in Astronomy and Astrophysics

FU Orionis-type objects are young, low-mass stars with large outbursts in visible light that last for several years or decades. They are thought to represent an evolutionary phase during the life of every young star when accretion from the circumstellar disk is enhanced during recurring time periods. These outbursts are able to rapidly build up the star while affecting the circumstellar disk and thus the ongoing or future planet formation. In many models infall from a circumstellar envelope seems to be necessary to trigger the outbursts. We observed the J=1$-$0 rotational transition of $^{13}$CO and C$^{18}$O towards eight northern FU Orionis-type stars (V1057 Cyg, V1515 Cyg, V2492 Cyg, V2493 Cyg, V1735 Cyg, V733 Cep, RNO 1B and RNO 1C) and derive temperatures and envelope masses and discuss the morphology and kinematics of the circumstellar material. We detected extended CO emission associated with all our targets. Smaller scale CO clumps were found to be associated with five objects with radii of 2000$-$5000 AU and masses of 0.02$-$0.5 $M_{\odot}$; these are clearly heated by the central stars. Three of these envelopes are also strongly detected in the 2.7 mm continuum. No central CO clumps were detected around V733 Cep and V710 Cas but there are many other clumps in their environments. Traces of outflow activity were observed towards V1735 Cyg, V733 Cep and V710 Cas. The diversity of the observed envelopes enables us to set up an evolutionary sequence between the objects. We find their evolutionary state to range from early, embedded Class I stage to late, Class II-type objects with very low-mass circumstellar material. The results reinforce the idea of FU Orionis-type stars as representatives of a transitory stage between embedded Class I young stellar objects and classical T-Tauri stars., 17 pages, 11 figures; accepted for publication in A&A

Context. DG Tau is a low-mass pre-main sequence star, whose strongly accreting protoplanetary disk exhibits a so-far enigmatic behavior: its mid-infrared thermal emission is strongly time-variable, even turning the 10 $��$m silicate feature from emission to absorption temporarily. Aims. We look for the reason for the spectral variability at high spatial resolution and at multiple epochs. Methods. We study the temporal variability of the mid-infrared interferometric signal, observed with the VLTI/MIDI instrument at six epochs between 2011 and 2014. We fit a geometric disk model to the observed interferometric signal to obtain spatial information about the disk. We also model the mid-infrared spectra by template fitting to characterize the profile and time dependence of the silicate emission. We use physically motivated radiative transfer modeling to interpret the mid-infrared interferometric spectra. Results. The inner disk (r1-3 au) spectra show a crystalline silicate feature in emission, similar to the spectra of comet Hale-Bopp. The striking difference between the inner and outer disk spectral feature is highly unusual among T Tauri stars. The mid-infrared variability is dominated by the outer disk. The strength of the silicate feature changed by more than a factor of two. Between 2011 and 2014 the half-light radius of the mid-infrared-emitting region decreased from 1.15 to 0.7 au. Conclusions. For the origin of the absorption we discuss four possible explanations: a cold obscuring envelope, an accretion heated inner disk, a temperature inversion on the disk surface and a misaligned inner geometry. The silicate emission in the outer disk can be explained by dusty material high above the disk plane, whose mass can change with time, possibly due to turbulence in the disk., 16 pages, 13 figures

Young stellar objects (YSOs) may undergo periods of active accretion (outbursts), during which the protostellar accretion rate is temporarily enhanced by a few orders of magnitude. Whether or not these accretion outburst YSOs possess similar dust/gas reservoirs to each other, and whether or not their dust/gas reservoirs are similar as quiescent YSOs, are issues not yet clarified. The aim of this work is to characterize the millimeter thermal dust emission properties of a statistically significant sample of long and short duration accretion outburst YSOs (i.e., FUors and EXors) and the spectroscopically identified candidates of accretion outbursting YSOs (i.e., FUor-like objects). We have carried out extensive Submillimeter Array (SMA) observations mostly at $\sim$225 GHz (1.33 mm) and $\sim$272 GHz (1.10 mm), from 2008 to 2017. We covered accretion outburst YSOs located at $$3-$\sigma$ significance. Detected sources except for the two cases of V883 Ori and NGC 2071 MM3 were observed with $\sim$1$"$ angular resolution. Overall our observed targets show a systematically higher millimeter luminosity distribution than those of the $M_{*}>$0.3 $M_{\odot}$ Class II YSOs in the nearby ($\lesssim$400 pc) low-mass star-forming molecular clouds (e.g., Taurus, Lupus, Upp Scorpio, and Chameleon I). In addition, at 1 mm our observed confirmed binaries or triple-system sources are systematically fainter than the rest of the sources even though their 1 mm fluxes are broadly distributed. We may have detected $\sim$30-60\% millimeter flux variability from V2494 Cyg and V2495 Cyg, from the observations separated by $\sim$1 year., Comment: 23 pages, 17 figures, accepted to A&A