Search

Context. Young stars accrete matter from their surrounding protoplanetary disk and drive powerful outflows. These two processes shape the final system architecture, and studying how these processes interact is the goal of the ESO Large programme PENELLOPE. PENELLOPE complements the ULLYSES legacy programme on the Hubble Space Telescope (HST) by providing ground-based – Very Large Telescope (VLT) – optical and near-infrared spectroscopy of more than 80 low-mass young stars. Aims. The main goal of this study is to screen the PENELLOPE/UVES targets for outflow activity and find microjets. A spectro-astrometric analysis in the [OI]λ6300 line in the velocity components of the microjet can give insights into the origin of the line emission, that is, if they originate from a magneto-hydrodynamical (MHD) wind or a photoevaporative wind. Methods. In total, 34 T Tauri stars of the PENELLOPE survey have been observed with the high-resolution optical slit spectrograph UVES (R ~ 65 000, λ = 3300–6800 Å). We formulated four criteria to rank the targets according to their outflow activity. Most of the targets have been observed in three different slit positions rotated by 120° with UVES. Using spectro-astrometric techniques in the [OI]λ6300 and Hα emission lines in each slit position of each target, we searched for outflow signatures, that is, an offset emission with respect to the continuum contribution of the associated T Tauri star. We checked all spectra for the presence of other wind line emission of [SII]λ4068, [SII]λ4076, [OI]λ5577, [OI]λ6300, [OI]λ6363, [SII]λ6716, [SII]λ6731, [NII]λ6548, and [NII]λ6583. Line profiles of Hα were inspected for a P Cygni signature. All [OI]λ6300 line profiles were decomposed into their constituent high-velocity component (HVC) and low-velocity component (LVC). Results. Our spectro-astrometric analysis in the [OI]λ6300 wind line reveals two newly discovered microjets associated with Sz 103 and XX Cha. Both microjets have an extent of about 0″.04, that is, <10 au, and we confined their orientation by the three slit observations. We identified two other interesting targets for which all four outflow criteria are fulfilled: Sz 98 and Sz 99. These targets display peculiar wind lines in their spectra with multiple velocity components, however, with the lack of a spectro-astrometric outflow signature. Furthermore, we confirm the binary nature of VWCha and CVSO 109. We present (further) evidence that DK Tau B and CVSO 104 A are spectroscopic binaries. Sz 115 is tentatively a spectroscopic binary. We find that the P Cygni line profile in the Hα line is not a robust indicator for the presence of outflows. Conclusions. The utilised observing strategy (rotating the UVES slit in three different positions) is very powerful in detecting micro-jets in T Tauri stars. The three slit positions can confine the spatial extend of the forbidden emission line regions. The introduced metric to rank targets according to their outflow activity is useful for follow-up observations. The origin of the LVC, that is, MHD winds versus photoevaporative winds, of the Sz 103 and XX Cha microjets remains unclear. [ABSTRACT FROM AUTHOR]

Context. Most ground-based direct-imaging planet search campaigns use angular differential imaging that distorts the signal from extended sources, such as protoplanetary disks. In the case of the young system PDS 70, for which two planets were detected within the cavity of a protoplanetary disk, obtaining a reliable image of both planets and the disk is essential to understanding planet-disk interactions. Aims. Our goals are to reveal the true intensity of the planets and disk without self-subtraction effects for the first time, search for new giant planets beyond separations of 0.1″, and to study the morphology of the disk shaped by two massive planets. Methods. We present YJHK-band imaging, polarimetry, and spatially resolved spectroscopy of PDS 70 using near-simultaneous reference star differential imaging, also known as star-hopping. We created a radiative transfer model of the system to try to match the near-infrared imaging and polarimetric data within measurement errors. Sub-millimeter imaging data from ALMA were also modeled. Furthermore, we extracted the spectra of the planets and the disk and compared them Results. With strong constraints, we find that the disk is quite flared, with a scale height of ∼15 at the outer edge of the disk at ∼90 au, similar to some disks in the literature. The gap inside ∼50 au is estimated to have ∼1 of the dust density of the outer disk. The northeast outer disk arc seen in previous observations is likely the outer lip of the flared disk. Abundance ratios of sub-micron, micron, and grains estimated by the modeling indicate a shallow grain-size index greater than −2.7, instead of the canonical –3.5. There is both vertical and radial segregation of grains. Planet c is well separated from the disk and has a spectrum similar to planet b, and it is clearly redder than the disk spectra. Planet c is possibly associated with the sudden flaring of the disk starting at ∼50 au. We found no new planets in the system. If we assume DUSTY models and an age of 5 Myr, this indicates no new planets more massive than 5 outside a 12 au separation. [ABSTRACT FROM AUTHOR]

Context. To understand how planets form in protoplanetary disks, it is necessary to characterize their gas and dust distribution and masses. This requires a combination of high-resolution dust continuum and molecular line interferometric observations, coupled with advanced theoretical models of protoplanetary disk physics, chemical composition, and radiative transfer. Aims. We aim to constrain the gas density and temperature distributions as well as gas masses in several T Tauri protoplanetary disks located in Taurus. We use the 12CO, 13CO, and C18O (2–1) isotopologue emission observed at 0.9″ with the IRAM NOrthern Extended Millimeter Array (NOEMA) as part of the MPG-IRAM Observatory Program PRODIGE (PROtostars and DIsks: Global Evolution PIs: P. Caselli & Th. Henning). Our sample consists of Class II disks with no evidence of strong radial substructures. We use these data to constrain the thermal and chemical structure of these disks through theoretical models for gas emission. Methods. To fit the combined optically thick and thin CO line data in Fourier space, we developed the DiskCheF code, which includes the parameterized disk physical structure, machine-learning (ML) accelerated chemistry, and the RADMC-3D line radiative transfer module. A key novelty of DiskCheF is the fast and feasible ML-based chemistry trained on the extended grid of the disk physical-chemical models precomputed with the ANDES2 code. This ML approach allows complex chemical kinetics models to be included in a time-consuming disk fitting without the need to run a chemical code. Results. We present a novel approach to incorporate chemistry into disk modeling without the need to explicitly calculate a chemical network every time. Using this new disk modeling tool, we successfully fit the 12CO, 13CO, and C18O (2-1) data from the CI, CY, DL, DM, DN, and IQ Tau disks. The combination of optically thin and optically thick CO lines allows us to simultaneously constrain the disk temperature and mass distribution, and derive the CO-based gas masses. The best-fit disk gas masses range between 0.005 and 0.04 M⊙. These values are in reasonable agreement with the disk dust masses rescaled by a factor of 100 as well as with other indirect gas measurements via, for example, modeling of the wavelength dependence of the dust continuum emission radii, and HD and CO isotopologue emission. [ABSTRACT FROM AUTHOR]

Context. The physics and chemistry of planet-forming disks are far from being fully understood. To make further progress, both broad line surveys and observations of individual tracers in a statistically significant number of disks are required. Aims. Our aim is to perform a line survey of eight planet-forming Class II disks in Taurus with the IRAM NOrthern Extended Millimeter Array (NOEMA), as a part of the MPG-IRAM Observatory Program PRODIGE (PROtostars and DIsks: Global Evolution; PIs: P. Caselli and Th. Henning). Methods. Compact and extended disks around T Tauri stars CI, CY, DG, DL, DM, DN, IQ Tau, and UZ Tau E are observed in ~80 lines from > 20 C-, O,- N-, and S-bearing species. The observations in four spectral settings at 210–280 GHz with a 1σ rms sensitivity of ~8–12 mJy beam−1 at a 0.9″ and 0.3 km s−1 resolution will be completed in 2024. The uv visibilities are fitted with the DiskFit model to obtain key stellar and disk properties. Results. In this first paper, the combined 12CO, 13CO, and C18O J = 2–1 data are presented. We find that the CO fluxes and disk masses inferred from dust continuum tentatively correlate with the CO emission sizes. We constrained dynamical stellar masses, geometries, temperatures, the CO column densities, and gas masses for each disk. The best-fit temperatures at 100 au are ~ 17–37 K, and decrease radially with the power-law exponent q ~ 0.05–0.76. The inferred CO column densities decrease radially with the power-law exponent p ~ 0.2–3.1. The gas masses estimated from 13CO (2–1) are ~0.001–0.2 M⊙. Conclusions. Using NOEMA, we confirm the presence of temperature gradients in our disk sample. The best-fit CO column densities point to severe CO freeze-out in these disks. The DL Tau disk is an outlier, and has either stronger CO depletion or lower gas mass than the rest of the sample. The CO isotopologue ratios are roughly consistent with the observed values in disks and the low-mass star-forming regions. The high 13CO/C18O ratio of ~23 in DM Tau could be indicative of strong selective photodissociation of C18O in this disk. [ABSTRACT FROM AUTHOR]

Aims. We aim to investigate the origin of the HI Brγ emission in young stars by using GRAVITY to image the innermost region of circumstellar disks, where important physical processes such as accretion and winds occur. With high spectral and angular resolution, we focus on studying the continuum and the HI Brγ-emitting area of the Herbig star HD 58647. Methods. Using VLTI-GRAVITY, we conducted observations of HD 58647 with both high spectral and high angular resolution. Thanks to the extensive uv coverage, we were able to obtain detailed images of the circumstellar environment at a sub-au scale, specifically capturing the continuum and the Brγ-emitting region. Through the analysis of velocity-dispersed images and photocentre shifts, we were able to investigate the kinematics of the HI Brγ-emitting region. Results. The recovered continuum images show extended emission where the disk major axis is oriented along a position angle of 14°. The size of the continuum emission at 5-σ levels is ~1.5 times more extended than the sizes reported from geometrical fitting (3.69 mas ± 0.02 mas). This result supports the existence of dust particles close to the stellar surface, screened from the stellar radiation by an optically thick gaseous disk. Moreover, for the first time with GRAVITY, the hot gas component of HD 58647 traced by the Brγ has been imaged. This allowed us to constrain the size of the Brγ-emitting region and study the kinematics of the hot gas; we find its velocity field to be roughly consistent with gas that obeys Keplerian motion. The velocity-dispersed images show that the size of the hot gas emission is from a more compact region than the continuum (2.3 mas ± 0.2 mas). Finally, the line phases show that the emission is not entirely consistent with Keplerian rotation, hinting at a more complex structure in the hot gaseous disk. [ABSTRACT FROM AUTHOR]

Context. The total disk mass sets the formation potential for exoplanets. Obtaining the disk mass is however not an easy feat, as one needs to consider the optical thickness, temperature, photodissociation, and freeze-out of potential mass tracers. Carbon-monoxide (CO) has been used as a gas mass tracer in T Tauri disks, but was found to be less abundant than expected due to the freeze-out and chemical conversion of CO on the surfaces of cold dust grains. The disks around more massive intermediate mass pre-main sequence stars called Herbig disks are likely to be warmer, allowing for the possibility of using CO as a more effective total gas mass tracer. Aims. This work aims to obtain the gas mass and size of Herbig disks observed with ALMA and compare these to previous works on T Tauri disks and debris disks. Methods. Using ALMA archival data and new NOEMA data of 12CO, 13CO, and C 18O transitions of 35 Herbig disks within 450 pc, the masses were determined using the thermo-chemical code Dust And Lines (DALI). A grid of models was run spanning five orders of magnitude in disk mass, for which the model CO line luminosities could be linked to the observed luminosities. Survival analysis was used to obtain cumulative distributions of the resulting disk masses. These were compared with dust masses from previous work to obtain gas-to-dust ratios for each disk. In addition, radii for all three isotopologues were obtained. Results. The majority of Herbig disks for which 13CO and C18O were detected are optically thick in both. For these disks, the line flux essentially only traces the disk size and only lower limits to the mass can be obtained. Computing the gas mass using a simple optically thin relation between line flux and column density results in an underestimate of the gas mass of at least an order of magnitude compared to the masses obtained with DALI. The inferred gas masses with DALI are consistent with a gas-to-dust ratio of at least 100. These gas-to-dust ratios are two orders of magnitude higher compared to those found for T Tauri disks using similar techniques, even over multiple orders of magnitude in dust mass, illustrating the importance of the chemical conversion of CO in colder T Tauri disks. Similar high gas-to-dust ratios are found for Herbig group I and II disks. Since group II disks have dust masses comparable to T Tauri disks, their higher CO gas masses illustrate the determining role of temperature. Compared to debris disks, Herbig disks have gas masses higher by four orders of magnitude. At least one Herbig disk, HD 163296, has a detected molecular disk wind, but our investigation has not turned up other detections of the CO disk wind in spite of similar sensitivities. Conclusions. Herbig disks are consistent with a gas-to-dust ratio of at least 100 over multiple orders of magnitude in dust mass. This indicates a fundamental difference between CO emission from Herbig disks and T Tauri disks, which is likely linked to the warmer temperature of the Herbig disks. [ABSTRACT FROM AUTHOR]

Context. A characteristic feature of young stellar objects is their variability, which is caused by a variety of different physical processes. High-resolution interferometric observations in the near- and mid-infrared wavelength ranges spanning multiple epochs allow the detailed study of these processes. Aims. We aim at investigating the expected variations of the interferometric observables connected to changes in the measured photometric fluxes of a typical variable accreting central young stellar object with a circumstellar disk. Methods. We calculated visibilities and closure phases as well as the photometric flux of brightness distributions obtained using 3D Monte Carlo radiative transfer simulations for a model of a circumstellar disk with an accreting central star. Results. Changes in the accretion luminosity of the central object, that is, an accreting pre-main-sequence star, can lead to significant variations in the visibility and closure phase of the star-disk system measured with instruments at the Very Large Telescope Interferometer (VLTI) that can be related to changes in the photometric flux. Taking into account additional effects due to baseline variation, interferometric observations can provide valuable contributions to the understanding of the underlying processes. Additionally, we provide the web application VLTI B-VAR that allows the impact of the hour angle on the visibility and closure phase for customized intensity maps to be estimated. [ABSTRACT FROM AUTHOR]

Context. Most studies devoted to Herbig Ae/Be stars (HAeBes) assume solar metallicity. However, the stellar metallicity, [M/H], is a fundamental parameter that can strongly differ depending on the source and may have important implications for planet formation. It has been proposed that the deficit of refractory elements observed in the surfaces of some HAeBes may be linked to the presence of cavities in their disks and is likely caused by Jovian planets that trap the metal-rich content. Aims. This work aims to provide a robust test on the previous proposal by analyzing the largest sample of HAeBes characterized by homogeneously derived [M/H] values and stellar and circumstellar properties. Methods. The spectra of 67 HAeBes, along with their well-known properties drawn from our previous work, have been collected from the ESO Science Archive Facility. Their [M/H] values were derived based on the comparison with Kurucz synthetic models. Statistical analyses were carried out with the aim to test the potential relation between [M/H] and the Meeus group I sources, with spectral energy distributions (SEDs) associated with the presence of cavities potentially carved by giant planets. We critically analyzed the eventual link between [M/H], the SED groups, and the presence of such planets. Results. Our statistical study robustly confirms that group I sources tend to have a lower [M/H] (typically ~ −0.10) than that of group II HAeBes (~ +0.14). A similar analysis involving SED-based transitional disks, with infrared excess only at wavelengths of ≥2.2 µm, does not reveal such a relation with [M/H]. This result indicates that not all processes capable of creating holes in the inner dust disks end up having an effect on the stellar abundances. The spatial distributions of group I and II sources are similar, at least within the available range of distances to the galactic centre and the galactic plane, for which the observed [M/H] differences are not driven by environmental effects. In addition, group I sources tend to have stronger (sub-) mm continuum emission presumably related to the presence of giant planets. Indeed, literature results indicate that disk substructures probably associated with the presence of giant planets are up to ten times more frequent in group I HAeBes than in group II. Finally, along with the metallicities derived for the whole sample, surface gravities and projected rotational velocities are additional outcomes reported in this work. Conclusions. We provide indirect evidence to suggest that giant planets are more frequent around group I/low [M/H] stars than around the rest of the HAeBes. However, a direct test of the previous hypothesis requires multiple detections of forming planets in their disks. Such detections have so far been limited to the candidate around the metal depleted ([M/H] = −0.35 ± −0.25) group I HAeBe star AB Aur, which is consistent with our findings. [ABSTRACT FROM AUTHOR]

Aims. We aim to assess the complementarity between spectroscopic and interferometric observations in the characterisation of the inner star-disc interaction region of young stars. Methods. We used the MCFOST code to solve the non-local thermodynamic equilibrium problem of line formation in non-axisymmetric accreting magnetospheres. We computed the Brγ line profile originating from accretion columns for models with different magnetic obliquities. We also derived monochromatic synthetic images of the Brγ line-emitting region across the line profile. This spectral line is a prime diagnostic of magnetospheric accretion in young stars and is accessible with the long baseline near-infrared interferometer GRAVITY installed at the ESO Very Large Telescope Interferometer. Results. We derive Brγ line profiles as a function of rotational phase and compute interferometric observables, visibilities, and phases, from synthetic images. The line profile shape is modulated along the rotational cycle, exhibiting inverse P Cygni profiles at the time the accretion shock faces the observer. The size of the line's emission region decreases as the magnetic obliquity increases, which is reflected in a lower line flux. We apply interferometric models to the synthetic visibilities in order to derive the size of the line-emitting region. We find the derived interferometric size to be more compact than the actual size of the magnetosphere, ranging from 50 to 90% of the truncation radius. Additionally, we show that the rotation of the non-axisymmetric magnetosphere is recovered from the rotational modulation of the Brγ-to-continuum photo-centre shifts, as measured by the differential phase of interferometric visibilities. Conclusions. Based on the radiative transfer modelling of non-axisymmetric accreting magnetospheres, we show that simultaneous spectroscopic and interferometric measurements provide a unique diagnostic to determine the origin of the Brγ line emitted by young stellar objects and are ideal tools to probe the structure and dynamics of the star-disc interaction region. [ABSTRACT FROM AUTHOR]

This paper is devoted to the study of gas-and-dust disks around young AB stars in the association Orion OB1. The observed spectral flux values of these stars in the wavelength range of 0.36–100 μm are determined. We used the results of multicolor photometric observations in the -bands, as well as data from the WISE and IRAS catalogs. It is shown that many program stars have radiation excesses in the far IR range. It is assumed that such an excess may be the result of the formation of a gas-dust disk around the program stars similar to those for Vega-type stars. [ABSTRACT FROM AUTHOR]

We discuss a technique for studying local gas inhomogeneities in the envelopes of young stars based on multicolor polarimetric observations. The method considered is based on the comparison of the observed time tracks of the variations in the Stokes parameters in the plane with the theoretical tracks computed in terms of a simple inhomogeneity model in the form of a gaseous stream with local magnetic field in the rotating envelope of a young star. The computations take into account the process of single scattering of radiation in the electron gas and the Faraday rotation of the polarization plane during propagation of radiation in a medium with magnetic field. The results of our computations depend strongly on the orbital inclination of the inhomogeneity to the line of sight and on the parameters of its magnetic field, whose influence on the polarization parameters increases with wavelength of radiation. The application of the method to the study of inhomogeneities in the envelopes of Herbig stars AB Aur and HD 163296 is described. [ABSTRACT FROM AUTHOR]

Context. It has been hypothesized that the location of Herbig Ae/Be stars (HAeBes) within the empirical relation between the inner disk radius (rin), inferred from K-band interferometry, and the stellar luminosity (L*), is related to the presence of the innermost gas, the disk-to-star accretion mechanism, the dust disk properties inferred from the spectral energy distributions (SEDs), or a combination of these effects. However, no general observational confirmation has been provided to date. Aims. This work aims to test whether the previously proposed hypotheses do, in fact, serve as a general explanation for the distribution of HAeBes in the size–luminosity diagram. Methods. GRAVITY/VLTI spectro-interferometric observations at ~2.2 μm have been obtained for five HBes representing two extreme cases concerning the presence of innermost gas and accretion modes. V590 Mon, PDS 281, and HD 94509 show no excess in the near-ultraviolet, Balmer region of the spectra (ΔDB), indicative of a negligible amount of inner gas and disk-to-star accretion, whereas DG Cir and HD 141926 show such strong ΔDB values that cannot be reproduced from magnetospheric accretion, but probably come from the alternative boundary layer mechanism. In turn, the sample includes three Group I and two Group II stars based on the Meeus et al. SED classification scheme. Additional data for these and all HAeBes resolved through K-band interferometry have been compiled from the literature and updated using Gaia EDR3 distances, almost doubling previous samples used to analyze the size–luminosity relation. Results. We find no general trend linking the presence of gas inside the dust destruction radius or the accretion mechanism with the location of HAeBes in the size–luminosity diagram. Similarly, our data do not support the more recent hypothesis linking such a location and the SED groups. Underlying trends are present and must be taken into account when interpreting the size–luminosity correlation. In particular, it cannot be statistically ruled out that this correlation is affected by dependencies of both L* and rin on the wide range of distances to the sources. Still, it is argued that the size–luminosity correlation is most likely to be physically relevant in spite of the previous statistical warning concerning dependencies on distance. Conclusions. Different observational approaches have been used to test the main scenarios proposed to explain the scatter of locations of HAeBes in the size–luminosity diagram. However, none of these scenarios have been confirmed as a fitting general explanation and this issue remains an open question. [ABSTRACT FROM AUTHOR]

Context. Young stars interact with their accretion disk through their strong magnetosphere. Aims. We aim to investigate the magnetospheric accretion/ejection process in the young stellar system DoAr 44 (V2062 Oph). Methods. We monitored the system over several rotational cycles, combining high-resolution spectropolarimetry at both optical and near-IR wavelengths with long-baseline near-IR inteferometry and multicolor photometry. Results. We derive a rotational period of 2.96 d from the system's light curve, which is dominated by stellar spots. We fully characterize the central star's properties from the high signal-to-noise, high-resolution optical spectra we obtained during the campaign. DoAr 44 is a young 1.2 M⊙ star, moderately accreting from its disk (Ṁacc = 6.5 10−9M⊙ yr−1), and seen at a low inclination (i ≃ 30°). Several optical and near-IR line profiles probing the accretion funnel flows (Hα, Hβ, HeI 1083 nm, Paβ) and the accretion shock (HeI 587.6 nm) are modulated at the stellar rotation period. The most variable line profile is HeI 1083 nm, which exhibits modulated redshifted wings that are a signature of accretion funnel flows, as well as deep blueshifted absorptions indicative of transient outflows. The Zeeman-Doppler analysis suggests the star hosts a mainly dipolar magnetic field, inclined by about 20° onto the spin axis, with an intensity reaching about 800 G at the photosphere, and up to 2 ± 0.8 kG close to the accretion shock. The magnetic field appears strong enough to disrupt the inner disk close to the corotation radius, at a distance of about 4.6 R⋆ (0.043 au), which is consistent with the 5 R⋆ (0.047 au) upper limit we derived for the size of the magnetosphere in our Paper I from long baseline interferometry. Conclusions. DoAr 44 is a pre-transitional disk system, exhibiting a 25–30 au gap in its circumstellar disk, with the inner and outer disks being misaligned. On a scale of 0.1 au or less, our results indicate that the system is steadily accreting from its inner disk through its tilted dipolar magnetosphere. We conclude that in spite of a highly structured disk on the large scale, perhaps the signature of ongoing planetary formation, the magnetospheric accretion process proceeds unimpeded at the star-disk interaction level. [ABSTRACT FROM AUTHOR]

Aims. We use optical spectroscopy to investigate the disk, wind, and accretion during the 2008 Z CMa NW outburst. Methods. Emission lines were used to constrain the locations, densities, and temperatures of the structures around the star. Results. More than 1000 optical emission lines reveal accretion, a variable, multicomponent wind, and double-peaked lines of disk origin. The variable, non-axisymmetric, accretion-powered wind has slow (~0 km s−1), intermediate (approximately −100 km s−1), and fast (≥−400 km s−1) components. The fast components are of stellar origin and disappear in quiescence, while the slow component is less variable and could be related to a disk wind. The changes in the optical depth of the lines between outburst and quiescence reveal that increased accretion is responsible for the observed outburst. We derive an accretion rate of 10−4M⊙ yr−1 in outburst. The Fe I and weak Fe II lines arise from an irradiated, flared disk at ~0.5–3 × M*/16 M⊙ au with asymmetric upper layers, revealing that the energy from the accretion burst is deposited at scales below 0.5 au. Some line profiles have redshifted asymmetries, but the system is unlikely to be sustained by magnetospheric accretion, especially in outburst. The accretion-related structures extend over several stellar radii and, like the wind, are likely to be non-axisymmetric. The stellar mass may be ~6–8 M⊙, lower than previously thought (~16 M⊙). Conclusions. Emission line analysis is found to be a powerful tool to study the innermost regions and accretion in stars within a very large range of effective temperatures. The density ranges in the disk and accretion structures are higher than in late-type stars, but the overall behavior, including the innermost disk emission and variable wind, is very similar for stars with different spectral types. Our work suggests a common outburst behavior for stars with spectral types ranging from M type to intermediate mass. [ABSTRACT FROM AUTHOR]

Context. AB Aur is a Herbig Ae star that hosts a prototypical transition disk. The disk shows a plethora of features connected with planet formation mechanisms, such as spiral arms, dust cavities, and dust traps. Understanding the physical and chemical characteristics of these features is crucial to advancing our knowledge of the planet formation processes. Aims. We aim to characterize the gaseous disk around the Herbig Ae star AB Aur. A complete spectroscopic study was performed using NOEMA to determine the physical and chemical conditions with high spatial resolution. Methods. We present new NOrthern Extended Millimeter Array (NOEMA) interferometric observations of the continuum and 12CO, 13CO, C18O, H2CO, and SO lines obtained at high resolution. We used the integrated intensity maps and stacked spectra to derive reliable estimates of the disk temperature. By combining our 13CO and C18O observations, we computed the gas-to-dust ratio along the disk. We also derived column density maps for the different species and used them to compute abundance maps. The results of our observations were compared with a set of Nautilus astrochemical models to obtain insight into the disk properties. Results. We detected continuum emission in a ring that extends from 0.6′′ to ~2.0′′, peaking at 0.97′′ and with a strong azimuthal asymmetry. The molecules observed show different spatial distributions, and the peaks of the distributions are not correlated with the binding energy. Using H2CO and SO lines, we derived a mean disk temperature of 39 K. We derived a gas-to-dust ratio that ranges from 10 to 40 along the disk. Abundance with respect to 13CO for SO (~2 × 10−4) is almost one order of magnitude greater than the value derived for H2CO (1.6 × 10−5). The comparison with Nautilus models favors a disk with a low gas-to-dust ratio (40) and prominent sulfur depletion. Conclusions. From a very complete spectroscopic study of the prototypical disk around AB Aur, we derived, for the first time, the gas temperature and the gas-to-dust ratio along the disk, providing information that is essential to constraining hydrodynamical simulations. Moreover, we explored the gas chemistry and, in particular, the sulfur depletion. The derived sulfur depletion is dependent on the assumed C/O ratio. Our data are better explained with C/O ~ 0.7 and S/H = 8 × 10−8. [ABSTRACT FROM AUTHOR]

Context. Classical T Tauri stars are pre-main sequence stars surrounded by an accretion disk. They host a strong magnetic field, and both magnetospheric accretion and ejection processes develop as the young magnetic star interacts with its disk. Studying this interaction is a major goal toward understanding the properties of young stars and their evolution. Aims. The goal of this study is to investigate the accretion process in the young stellar system HQ Tau, an intermediate-mass T Tauri star (1.9 M⊙). Methods. The time variability of the system is investigated both photometrically, using Kepler-K2 and complementary light curves, and from a high-resolution spectropolarimetric time series obtained with ESPaDOnS at CFHT. Results. The quasi-sinusoidal Kepler-K2 light curve exhibits a period of 2.424 d, which we ascribe to the rotational period of the star. The radial velocity of the system shows the same periodicity, as expected from the modulation of the photospheric line profiles by surface spots. A similar period is found in the red wing of several emission lines (e.g., HI, CaII, NaI), due to the appearance of inverse P Cygni components, indicative of accretion funnel flows. Signatures of outflows are also seen in the line profiles, some being periodic, others transient. The polarimetric analysis indicates a complex, moderately strong magnetic field which is possibly sufficient to truncate the inner disk close to the corotation radius, rcor ∼ 3.5 R⋆. Additionally, we report HQ Tau to be a spectroscopic binary candidate whose orbit remains to be determined. Conclusions. The results of this study expand upon those previously reported for low-mass T Tauri stars, as they indicate that the magnetospheric accretion process may still operate in intermediate-mass pre-main sequence stars, such as HQ Tauri. [ABSTRACT FROM AUTHOR]

Context. Studies of the dust distribution, composition, and evolution of protoplanetary disks provide clues for understanding planet formation. However, little is known about the innermost regions of disks where telluric planets are expected to form. Aims. We aim constrain the geometry of the inner disk of the T Tauri star RY Lup by combining spectro-photometric data and interferometric observations in the near-infrared (NIR) collected at the Very Large Telescope Interferometer. We use PIONIER data from the ESO archive and GRAVITY data that were obtained in June 2017 with the four 8m telescopes. Methods. We use a parametric disk model and the 3D radiative transfer code MCFOST to reproduce the spectral energy distribution (SED) and match the interferometric observations. MCFOST produces synthetic SEDs and intensity maps at different wavelengths from which we compute the modeled interferometric visibilities and closure phases through Fourier transform. Results. To match the SED from the blue to the millimetric range, our model requires a stellar luminosity of 2.5 L⊙, higher than any previously determined values. Such a high value is needed to accommodate the circumstellar extinction caused by the highly inclined disk, which has been neglected in previous studies. While using an effective temperature of 4800 K determined through high-resolution spectroscopy, we derive a stellar radius of 2.29 R⊙. These revised fundamental parameters, when combined with the mass estimates available (in the range 1.3–1.5 M⊙), lead to an age of 0.5–2.0 Ma for RY Lup, in better agreement with the age of the Lupus association than previous determinations. Our disk model (that has a transition disk geometry) nicely reproduces the interferometric GRAVITY data and is in good agreement with the PIONIER ones. We derive an inner rim location at 0.12 au from the central star. This model corresponds to an inclination of the inner disk of 50°, which is in mild tension with previous determinations of a more inclined outer disk from SPHERE (70° in NIR) and ALMA (67 ± 5°) images, but consistent with the inclination determination from the ALMA CO spectra (55 ± 5°). Increasing the inclination of the inner disk to 70° leads to a higher line-of-sight extinction and therefore requires a higher stellar luminosity of 4.65 L⊙ to match the observed flux levels. This luminosity would translate to a stellar radius of 3.13 R⊙, leading to an age of 2–3 Ma, and a stellarmass of about 2 M⊙, in disagreement with the observed dynamical mass estimate of 1.3–1.5 M⊙. Critically, this high-inclination inner disk model also fails to reproduce the visibilities observed with GRAVITY. Conclusions. The inner dust disk, as traced by the GRAVITY data, is located at a radius in agreement with the dust sublimation radius. An ambiguity remains regarding the respective orientations of the inner and outer disk, coplanar and mildly misaligned, respectively.As our datasets are not contemporary and the star is strongly variable, a deeper investigation will require a dedicated multi-technique observing campaign. [ABSTRACT FROM AUTHOR]

Aims. We aim to study vertical settling and radial drift of dust in protoplanetary disks from a different perspective: an edge-on view. An estimation of the amplitude of settling and drift is highly relevant to understanding planet formation. Methods. We analyze a sample of 12 HST-selected edge-on protoplanetary disks (i.e., seen with high inclinations) for which the vertical extent of the emission layers can be constrained directly. We present ALMA high angular resolution continuum images (~0.1′′) of these disks at two wavelengths, 0.89 and 2.06 mm (respectively ALMA bands 7 and 4), supplemented with archival band 6 data (1.33 mm) where available. Results. Several sources show constant brightness profiles along their major axis with a steep drop at the outer edges. Two disks have central holes with additional compact continuum emission at the location of the central star. For most sources, the millimeter continuum emission is more compact than the scattered light, both in the vertical and radial directions. Six sources are resolved along their minor axis in at least one millimetric band, providing direct information on the vertical distribution of the millimeter grains. For the second largest disk of the sample, Tau 042021, the significant difference in vertical extent between band 7 and band 4 suggests efficient size-selective vertical settling of large grains. Furthermore, the only Class I object in our sample shows evidence of flaring in the millimeter. Along the major axis, all disks are well resolved. Four of them are larger in band 7 than in band 4 in the radial direction, and three have a similar radial extent in all bands. These three disks are also the ones with the sharpest apparent edges (between 80% and 20% of the peak flux, Δr∕r ~ 0.3), and two of them are binaries. For all disks, we also derive the millimeter brightness temperature and spectral index maps. We find that all edge-on disks in our sample are likely optically thick and that the dust emission reveals low brightness temperatures in most cases (brightness temperatures ≤10 K). The integrated spectral indices are similar to those of disks at lower inclination. Conclusions. The comparison of a generic radiative transfer disk model with our data shows that at least three disks are consistent with a small millimeter dust scale height, of a few au (measured at r = 100 au). This is in contrast with the more classical value of hg ~ 10 au derived from scattered light images and from gas line measurements. These results confirm, by direct observations, that large (millimeter) grains are subject to significant vertical settling in protoplanetary disks. [ABSTRACT FROM AUTHOR]

Context. Optical and infrared spatially unresolved multi-epoch observations have revealed the variability of pre-main sequence stars and/or their environment. Moreover, structures in orbital motion around the central star, resulting from planet-disk interaction, are predicted to cause temporal variations in the brightness distributions of protoplanetary disks. Through repeated observations of pre-main sequence stars with the Very Large Telescope Interferometer (VLTI) over nearly two decades, the ESO Archive has become a treasure chest containing unprecedented high-resolution multi-epoch near- and mid-infrared observations of the potential planet-forming regions in protoplanetary disks. Aims. We aim to investigate whether the existing multi-epoch observations provide evidence for the variability of the brightness distributions of the innermost few astronomical units of protoplanetary disks and to quantify any variations detected. Methods. We present different approaches to search for evidence of temporal variations based on multi-epoch observations obtained with the VLTI instruments PIONIER, AMBER, and MIDI for 68 pre-main sequence stars. Results. For nine objects in our sample, multi-epoch data obtained using equal baselines are available that allow us to directly detect variations in the visibilities due to temporally variable brightness distributions. Significant variations of the near-infrared visibilities obtained in different epochs with PIONIER and/or AMBER for HD 50138, DX Cha, HD 142527, V856 Sco, HD 163296, and R CrA were found. HD 37806, TW Hya, and CPD-36 6759 show no significant variations. By estimating the impact of a small variation of the baseline on the measured squared visibilities, we are able to compare the data of another 12 pre-main sequence stars. Thereby, we find evidence for temporal variations of the brightness distribution of one additional object, AK Sco. Besides the two binaries DX Cha and AK Sco, HD 50138 and V856 Sco also show signs of variability caused by variations of asymmetric structures in the brightness distribution. [ABSTRACT FROM AUTHOR]

At the time of stellar evolution, young stellar objects go through processes of increased activity and instability. Star formation takes place in several stages during which the star accumulates enough mass to initiate thermonuclear reactions in the nucleus. A significant percentage of the mass of Sun-like stars accumulates during periods of increased accretion known as FUor outbursts. Since we know only about two dozen stars of this type, the study of each new object is very important for our knowledge. In this paper, we present data from photometric monitoring on a FUor object V2493 Cyg discovered in 2010. Our data were obtained in the optical region with BVRI Johnson–Cousins set of filters during the period from November 2016 to February 2021. The results of our observations show that during this period no significant changes in the brightness of the star were registered. We only detect variations with a small amplitude around the maximum brightness value. Thus, since 2013 V2493 Cyg remains at its maximum brightness, without a decrease in brightness. Such photometric behavior is not typical of other stars from FUor type. Usually, the light curves of FUors are asymmetrical, with a very rapid rise and gradual decline of the brightness. V2493 Cyg remains unique in this respect with a very rapid rise in brightness and prolonged retention in maximum light. Our period analysis made for the interval February 2013–February 2021 reveals a well-defined period of 914 ± 10 days. Such periodicity can be explained by dust structures remaining from star formation in orbit around the star.

Context. Pre-main-sequence variability characteristics can be used to probe the physical processes leading to the formation and initial evolution of both stars and planets. Aims. The photometric variability of pre-main-sequence stars is studied at optical wavelengths to explore star–disk interactions, accretion, spots, and other physical mechanisms associated with young stellar objects. Methods. We observed a field of 16′ × 16′ in the star-forming region Pelican Nebula (IC 5070) at BVRI wavelengths for 90 nights spread over one year in 2012−2013. More than 250 epochs in the VRI bands are used to identify and classify variables up to V ∼ 21 mag. Their physical association with the cluster IC 5070 is established based on the parallaxes and proper motions from the Gaia second data release (DR2). Multiwavelength photometric data are used to estimate physical parameters based on the isochrone fitting and spectral energy distributions. Results. We present a catalog of optical time-series photometry with periods, mean magnitudes, and classifications for 95 variable stars including 67 pre-main-sequence variables towards star-forming region IC 5070. The pre-main-sequence variables are further classified as candidate classical T Tauri and weak-line T Tauri stars based on their light curve variations and the locations on the color-color and color-magnitude diagrams using optical and infrared data together with Gaia DR2 astrometry. Classical T Tauri stars display variability amplitudes up to three times the maximum fluctuation in disk-free weak-line T Tauri stars, which show strong periodic variations. Short-term variability is missed in our photometry within single nights. Several classical T Tauri stars display long-lasting (≥10 days) single or multiple fading and brightening events of up to two magnitudes at optical wavelengths. The typical mass and age of the pre-main-sequence variables from the isochrone fitting and spectral energy distributions are estimated to be ≤1 M⊙ and ∼2 Myr, respectively. We do not find any correlation between the optical amplitudes or periods with the physical parameters (mass and age) of pre-main-sequence stars. Conclusions. The low-mass pre-main-sequence stars in the Pelican Nebula region display distinct variability and color trends and nearly 30% of the variables exhibit strong periodic signatures attributed to cold spot modulations. In the case of accretion bursts and extinction events, the average amplitudes are larger than one magnitude at optical wavelengths. These optical magnitude fluctuations are stable on a timescale of one year. [ABSTRACT FROM AUTHOR]

Context. T Tauri stars are low-mass young stars whose disks provide the setting for planet formation. Despite this, their structure is poorly understood. We present new infrared interferometric observations of the SU Aurigae circumstellar environment that offer resolution that is three times higher and a better baseline position angle coverage than previous observations. Aims. We aim to investigate the characteristics of the circumstellar material around SU Aur, constrain the disk geometry, composition and inner dust rim structure. Methods. The CHARA array offers unique opportunities for long baseline observations, with baselines up to 331 m. Using the CLIMB three-telescope combiner in the K-band allows us to measure visibilities as well as closure phase. We undertook image reconstruction for model-independent analysis, and fitted geometric models such as Gaussian and ring distributions. Additionally, the fitting of radiative transfer models constrain the physical parameters of the disk. For the first time, a dusty disk wind is introduced to the radiative transfer code TORUS to model protoplanetary disks. Our implementation is motivated by theoretical models of dusty disk winds, where magnetic field lines drive dust above the disk plane close to the sublimation zone. Results. Image reconstruction reveals an inclined disk with slight asymmetry along its minor-axis, likely due to inclination effects obscuring the inner disk rim through absorption of incident star light on the near-side and thermal re-emission and scattering of the far-side. Geometric modelling of a skewed ring finds the inner rim at 0.17 ± 0.02 au with an inclination of 50.9 ± 1.0° and minor axis position angle 60.8 ± 1.2°. Radiative transfer modelling shows a flared disk with an inner radius at 0.18 au which implies a grain size of 0.4 μm assuming astronomical silicates and a scale height of 15.0 at 100 au. Among the tested radiative transfer models, only the dusty disk wind successfully accounts for the K-band excess by introducing dust above the mid-plane. [ABSTRACT FROM AUTHOR]

Context. In recent years, the disk populations in a number of young star-forming regions have been surveyed with the Atacama Large Millimeter/submillimeter Array (ALMA). Understanding the disk properties and their correlation with the properties of the central star is critical to understanding planet formation. In particular, a decrease of the average measured disk dust mass with the age of the region has been observed, consistent with grain growth and disk dissipation. Aims. We aim to compare the general properties of disks and their host stars in the nearby (d = 160 pc) Corona Australis (CrA) star forming region to those of the disks and stars in other regions. Methods. We conducted high-sensitivity continuum ALMA observations of 43 Class II young stellar objects in CrA at 1.3 mm (230 GHz). The typical spatial resolution is ~0.3′′. The continuum fluxes are used to estimate the dust masses of the disks, and a survival analysis is performed to estimate the average dust mass. We also obtained new VLT/X-shooter spectra for 12 of the objects in our sample for which spectral type (SpT) information was missing. Results. Twenty-four disks were detected, and stringent limits have been put on the average dust mass of the nondetections. Taking into account the upper limits, the average disk mass in CrA is 6 ± 3 M⊕. This value is significantly lower than that of disks in other young (1–3 Myr) star forming regions (Lupus, Taurus, Chamaeleon I, and Ophiuchus) and appears to be consistent with the average disk mass of the 5–10 Myr-old Upper Sco. The position of the stars in our sample on the Herzsprung-Russel diagram however seems to confirm that CrA has an age similar to Lupus. Neither external photoevaporation nor a lower-than-usual stellar mass distribution can explain the low disk masses. On the other hand, a low-mass disk population could be explained if the disks were small, which could happen if the parent cloud had a low temperature or intrinsic angular momentum, or if the angular momentum of the cloud were removed by some physical mechanism such as magnetic braking. Even in detected disks, none show clear substructures or cavities. Conclusions. Our results suggest that in order to fully explain and understand the dust mass distribution of protoplanetary disks and their evolution, it may also be necessary to take into consideration the initial conditions of star- and disk-formation process. These conditions at the very beginning may potentially vary from region to region, and could play a crucial role in planet formation and evolution. [ABSTRACT FROM AUTHOR]

Context. The complex system HD 100453 AB with a ring-like circumprimary disk and two spiral arms, one of which is pointing to the secondary, is a good laboratory in which to test spiral formation theories. Aims. We aim to resolve the dust and gas distribution in the disk around HD 100453 A and to quantify the interaction of HD 100453 B with the circumprimary disk. Methods. Using ALMA band 6 dust continuum and CO isotopologue observations we have studied the HD 100453 AB system with a spatial resolution of 0′′.09 × 0′′.17 at 234 GHz. We used smoothed particle hydrodynamics (SPH) simulations and orbital fitting to investigate the tidal influence of the companion on the disk. Results. We resolve the continuum emission around HD 100453 A into a disk between 0′′.22 and 0′′.40 with an inclination of 29.5° and a position angle of 151.0°, an unresolved inner disk, and excess mm emission cospatial with the northern spiral arm which was previously detected using scattered light observations. We also detect CO emission from 7 au (well within the disk cavity) out to 1′′.10, overlapping with HD 100453 B at least in projection. The outer CO disk position angle (PA) and inclination differ by up to 10° from the values found for the inner CO disk and the dust continuum emission, which we interpret as due to gravitational interaction with HD 100453 B. Both the spatial extent of the CO disk and the detection of mm emission at the same location as the northern spiral arm are in disagreement with the previously proposed near co-planar orbit of HD 100453 B. Conclusions. We conclude that HD 100453 B has an orbit that is significantly misaligned with the circumprimary disk. Because it is unclear whether such an orbit can explain the observed system geometry we highlight an alternative scenario that explains all detected disk features where another, (yet) undetected, low mass close companion within the disk cavity, shepherds a misaligned inner disk whose slowly precessing shadows excite the spiral arms. [ABSTRACT FROM AUTHOR]

Context. The disks around some Herbig Be stars have been observed to be more compact than the expected dust sublimation radius for such objects. Highly refractory dust grains and optically thick gas emission have been proposed as possible explanations for this phenomenon. Aims. Previously, the "undersized" Herbig Be star MWC 147 was observed with interferometry, and the results indicated a presence of a compact gaseous disk based on the measured wavelength-dependence of near-infrared or mid-infrared visibilities. Our aim is to search for direct evidence for the presence of hot gas inside of the expected dust sublimation radius of MWC 147. Methods. By combining VLTI/AMBER spectro-interferometry (R = 12 000) with CRIRES spectroscopy (R = 100 000) we can both spectrally and spatially resolve the Brγ line-emitting gas around MWC 147. Additionally, using CHARA/CLIMB enables us to achieve baseline lengths up to 330 m, offering ~2 times higher angular resolution (and a better position angle coverage) than has previously been achieved with interferometry for MWC 147. To model the continuum we fit our AMBER and CLIMB measurements with a geometric model of an inclined Gaussian distribution as well as a ring model. We fit our high-resolution spectra and spectro-interferometric data with a kinematic model of a disk in Keplerian rotation. Results. Our interferometric visibility modelling of MWC 147 indicates the presence of a compact continuum disk with a close to face-on orientation. We model the continuum with an inclined Gaussian and a ring with a radius of 0.60 mas (0.39 au), which is well within the expected dust sublimation radius of 1.52 au. We detect no significant change in the measured visibilities across the Brγ line, indicating that the line-emitting gas is located in the same region as the continuum-emitting disk. Using our differential phase data we construct photocentre displacement vectors across the Brγ line, revealing a velocity profile consistent with a rotating disk. We fit our AMBER spectro-interferometry data with a kinematic model of a disk in Keplerian rotation, where both the line-emitting and continuum-emitting components of the disk originate from the same compact region close to the central star. The presence of line-emitting gas in the same region as the K-band continuum supports the interpretation that the K-band continuum traces an optically thick gas disk. Conclusions. Our spatially and spectrally resolved observations of MWC 147 reveal that the K-band continuum and Brγ emission both originate from a similar region which is 3.9 times more compact than the expected dust sublimation radius for the star; Brγ is emitted from the accretion disk or disk wind region and exhibits a rotational velocity profile. We conclude that we detect the presence of a compact, gaseous accretion disk in Keplerian rotation around MWC 147. [ABSTRACT FROM AUTHOR]

Context. Infrared and (sub-)millimeter observations of disks around T Tauri and Herbig Ae/Be stars point to a chemical differentiation, with a lower detection rate of molecules in disks around hotter stars. Aims. We aim to investigate the underlying causes of the chemical differentiation indicated by observations and perform a comparative study of the chemistry of T Tauri and Herbig Ae/Be disks. This is one of the first studies to compare the chemistry in the outer regions of these two types of disk. Methods. We developed a model to compute the chemical composition of a generic protoplanetary disk, with particular attention to the photochemistry, and applied it to a T Tauri and a Herbig Ae/Be disk. We compiled cross sections and computed photodissociation and photoionization rates at each location in the disk by solving the far-ultraviolet (FUV) radiative transfer in a 1+1D approach using the Meudon PDR code and adopting observed stellar spectra. Results. The warmer disk temperatures and higher ultraviolet flux of Herbig stars compared to T Tauri stars induce some differences in the disk chemistry. In the hot inner regions, H2O, and simple organic molecules like C2H2, HCN, and CH4 are predicted to be very abundant in T Tauri disks and even more in Herbig Ae/Be disks, in contrast with infrared observations that find a much lower detection rate of water and simple organics toward disks around hotter stars. In the outer regions, the model indicates that the molecules typically observed in disks, like HCN, CN, C2H, H2CO, CS, SO, and HCO+, do not have drastic abundance differences between T Tauri and Herbig Ae disks. Some species produced under the action of photochemistry, like C2H and CN, are predicted to have slightly lower abundances around Herbig Ae stars due to a narrowing of the photochemically active layer. Observations indeed suggest that these radicals are somewhat less abundant in Herbig Ae disks, although in any case, the inferred abundance differences are small, of a factor of a few at most. A clear chemical differentiation between both types of disks concerns ices. Owing to the warmer temperatures of Herbig Ae disks, one expects snow lines lying farther away from the star and a lower mass of ices compared to T Tauri disks. Conclusions. The global chemical behavior of T Tauri and Herbig Ae/Be disks is quite similar. The main differences are driven by the warmer temperatures of the latter, which result in a larger reservoir or water and simple organics in the inner regions and a lower mass of ices in the outer disk. [ABSTRACT FROM AUTHOR]

We are totally immersed in the Big Data era and reliable algorithms and methods for data classification are instrumental for astronomical research. Random Forest and Support Vector Machines algorithms have become popular over the last few years and they are widely used for different stellar classification problems. In this article, we explore an alternative supervised classification method scarcely exploited in astronomy, Logistic Regression, that has been applied successfully in other scientific areas, particularly biostatistics. We have applied this method in order to derive membership probabilities for potential T Tauri star candidates from ultraviolet-infrared colour-colour diagrams. [ABSTRACT FROM AUTHOR]

This is the final photometric study of TW Hya based on new MOST satellite observations. During 2014 and 2017, the light curves showed stable 3.75- and  3.69-d quasi-periodic oscillations, respectively. Both values appear to be closely related to the stellar rotation period, as they might be created by changing visibility of a hotspot formed near the magnetic pole directed towards the observer. These major light variations were superimposed on a chaotic, flaring-type activity caused by hotspots resulting from unstable accretion – a situation reminiscent of that in 2011, when TW Hya showed signs of a moderately stable accretion state. In 2015, only drifting quasi-periods were observed, similar to those present in 2008–2009 data and typical for magnetized stars accreting in a strongly unstable regime. A rich set of multicolour data was obtained during 2013–2017 with the primary aim of characterizing the basic spectral properties of the mysterious occultations in TW Hya. Although several possible occultation-like events were identified, they are not as well defined as in the 2011 MOST data. The new ground-based and MOST data show a dozen previously unnoticed flares, as well as small-amplitude 11 min–3 h brightness variations, associated with 'accretion bursts'. We cannot exclude the possibility that the shortest 11–15 min variations could also be caused by thermal instability oscillations in an accretion shock. [ABSTRACT FROM AUTHOR]

Context. Eruptive variables are young stars that show episodic variations of brightness: EXors/FUors variations are commonly associated with enhanced accretion outbursts occurring at intermittent cadence of months/years (EXors) and decades/centuries (FUors). Variations that can be ascribed to a variable extinction along their line of sight are instead classified as UXors. Aims. We aim at investigating the long-term photometric behaviour of three sources classified as eruptive variables. We present data from the archival plates of the Asiago Observatory relative to the fields where the targets are located. For the sake of completeness we have also analysed the Harvard plates of the same regions that cover a much longer historical period, albeit at a lower sensitivity, however we are only able to provide upper limits. Methods. A total of 273 Asiago plates were investigated, providing a total of more than 200 magnitudes for the three stars, which cover a period of about 34 yr between 1958 and 1991. We have compared our data with more recently collected literature data. Results. Our plates analysis of V2492 Cyg provides historical upper limits that seem not to be compatible with the level of the activity monitored during the last decade. Therefore, recently observed accretion phenomena could be associated with the outbursting episodes, more than repetitive obscuration. While a pure extinction does not seem the only mechanism responsible for the ASASSN-15qi fluctuations, it can account quite reasonably for the recent V350 Cep variations. [ABSTRACT FROM AUTHOR]

Disk warps around classical T Tauri stars (CTTSs) can periodically obscure the central star for some viewing geometries. For these so- called AA Tau-like variables, the obscuring material is located in the inner disk and absorption spectroscopy allows one to characterize its dust and gas content. Since the observed emission from CTTSs consists of several components (photospheric, accretion, jet, and disk emission), which can all vary with time, it is generally challenging to disentangling disk features from emission variability. Multi- epoch, flux-calibrated, broadband spectra provide us with the necessary information to cleanly separate absorption from emission variability. We applied this method to three epochs of VLT/X-shooter spectra of the CTTS V 354 Mon (CSI Mon-660) located in NGC 2264 and find that: (a) the accretion emission remains virtually unchanged between the three epochs; (b) the broadband flux evolution is best described by disk material obscuring part of the star, and (c) the Na and K gas absorption lines show only a minor increase in equivalent width during phases of high dust extinction. The limits on the absorbing gas column densities indicate a low gas-to-dust ratio in the inner disk, less than a tenth of the ISM value. We speculate that the evolutionary state of V 354 Mon, rather old with a low accretion rate, is responsible for the dust excess through an evolution toward a dust dominated disk or through the fragmentation of larger bodies that drifted inward from larger radii in a still gas dominated disk. [ABSTRACT FROM AUTHOR]

We have quantified the potential capabilities of detecting local brightness asymmetries in circumstellar disks with the Very Large Telescope Interferometer (VLTI) in the mid-infrared wavelength range. The study is motivated by the need to evaluate theoretical models of planet formation by direct observations of protoplanets at early evolutionary stages, when they are still embedded in their host disk. Up to now, only a few embedded candidate protoplanets have been detected with semi-major axes of 20–50 au. Due to the small angular separation from their central star, only long-baseline interferometry provides the angular resolving power to detect disk asymmetries associated to protoplanets at solar system scales in nearby star-forming regions. In particular, infrared observations are crucial to observe scattered stellar radiation and thermal re-emission in the vicinity of embedded companions directly. For this purpose we performed radiative transfer simulations to calculate the thermal re-emission and scattered stellar flux from a protoplanetary disk hosting an embedded companion. Based on that, visibilities and closure phases are calculated to simulate observations with the future beam combiner MATISSE, operating at the L, M and N bands at the VLTI. We find that the flux ratio of the embedded source to the central star can be as low as 0.5 to 0.6% for a detection at a feasible significance level due to the heated dust in the vicinity of the embedded source. Furthermore, we find that the likelihood for detection is highest for sources at intermediate distances r ≈ 2–5 au and disk masses not higher than ≈10−4 M⊙. [ABSTRACT FROM AUTHOR]

Context. V2492 Cyg is a young pre-main sequence star presenting repetitive brightness variations of significant amplitude (ΔR ⩾ 5 mag) whose physical origin has been ascribed to both extinction (UXor-type) and accretion (EXor-type) variability, although their mutual proportion has not been clarified yet. Recently, V2492 Cyg has reached a level of brightness ever registered in the period of its documented activity. Aims. We aim to derive the variation of the mass accretion rate between low- and high-state and to get new insights on the origin of the variability of V2492 Cyg. Methods. Optical and near-infrared (NIR) photometry and spectroscopy have been obtained in October 2016 and between March and July 2017. The source has remained bright until the end of May 2017, then it started to rapidly fade since the beginning of June at a rate of ∼0.08 mag/day. On mid-July 2017 the source has reached the same low-brightness level as two years before. Extinction and mass accretion rate were derived by means of the luminosity of the brightest lines, in particular Hα and Hβ A couple of optical high-resolution spectra are also presented to derive information on the gas kinematics. Results. Visual extinction variations do not exceed a few magnitudes, while the mass accretion rate is estimated to vary from less than 10-8 up to a few 10-7 M☉ yr-1. This latter is comparable to that estimated on the previous high-state in 2010, likely occurred under more severe extinction conditions. Conclusions. The combined analysis of the optical and NIR observations extends to the present event the original suggestion that the V2492 Cyg variability is a combination of changing extinction and accretion. [ABSTRACT FROM AUTHOR]

We report on the discovery of correlations between dust and CO gas tracers of the 0.1-10 au region in planet-forming disks around young intermediate-mass stars. The abundance of refractory elements on stellar photospheres decreases as the location of hot CO gas emission recedes to larger disk radii, and as the near-infrared excess emission from hot dust in the inner disk decreases. The linked behavior between these observables demonstrates that the recession of infrared CO emission to larger disk radii traces an inner disk region where dust is being depleted. We also find that Herbig disk cavities have either low (~5-10%) or high (~20-35%) near-infrared excess, a dichotomy that has not been captured by the classic definition of "pre-transitional" disks. [ABSTRACT FROM AUTHOR]

Protoplanets can produce structures in protoplanetary disks via gravitational disk–planet interactions. Once detected, such structures serve as signposts of planet formation. Here we investigate the kinematic signatures in disks produced by multi-Jupiter mass ( M _J ) planets using 3D hydrodynamics and radiative transfer simulations. Such a planet opens a deep gap, and drives transonic vertical motions inside. Such motions include both a bulk motion of the entire half-disk column, and turbulence on scales comparable to and smaller than the scale height. They significantly broaden molecular lines from the gap, producing double-peaked line profiles at certain locations, and a kinematic velocity dispersion comparable to thermal after azimuthal averaging. The same planet does not drive fast vertical motions outside the gap, except at the inner spiral arms and the disk surface. Searching for line broadening induced by multi- M _J planets inside gaps requires an angular resolution comparable to the gap width, an assessment of the gap gas temperature to within a factor of 2, and a high sensitivity needed to detect line emission from the gap.

Having compared images of a jet of the young star RWAurA obtained with an interval of 21.3 yr, we have found that the outermost knots of the jet have emerged approximately 350 years ago. We come up with arguments that the jet itself has appeared at the same time, and intensive accretion onto the star has begun due to rearrangement of its protoplanetary disk structure caused by the tidal effect of the companion RWAur B. More precisely suppose that intensification of accretion is a response to changing conditions in the outer-disk regions which has followed after the sound wave, generated by these changes, has passed the disk in the radial direction. In our opinion difference in the parameters of blue and red lobes of the RWAurA jet is a result of the asymmetric distribution of the circumstellar matter above and below the disk due to companion's passage. It was found from the analysis of the RWAur historical light curve that deep and long-term (Δ t > 150 days) light attenuations of RWAurA observed after 2010 had no precedents in the previous 110 years.We also associate the change in the character of photometric variability of the star with the rearrangement of the structure of inner ( r < 1 AU) regions of its protoplanetary disk, and discuss why these changes have begun only 350 years after the beginning of the active accretion phase. [ABSTRACT FROM AUTHOR]

We explore how measurements of protoplanetary disc masses and accretion rates provided by surveys of star-forming regions can be analysed via the dimensionless accretion parameter, which we define as the product of the accretion rate and stellar age divided by the disc mass. By extending and generalizing the study of Jones et al., we demonstrate that this parameter should be less than or of order unity for a wide range of evolutionary scenarios, rising above unity only during the final stages of outside-in clearing by external photoevaporation. We use this result to assess the reliability of disc mass estimates derived from CO isotopologues and sub-mm continuum emission by examining the distribution of accretion efficiencies in regions that are not subject to external photoevaporation.We find that while dust-based mass estimates produce results compatible with theoretical expectations assuming a canonical dust-to-gas ratio, the systematically lower CO-based estimates yield accretion efficiencies significantly above unity in contrast with the theory. This finding provides additional evidence that CO-based disc masses are an underestimate, in line with arguments that have been made on the basis of chemical modelling of relatively small samples. On the other hand, we demonstrate that dust-based mass estimates are sufficiently accurate to reveal distinctly higher accretion efficiencies in the Trapezium cluster, where this result is expected, given the evident importance of external photoevaporation. We therefore propose the dimensionless accretion parameter as a new diagnostic of external photoevaporation in other star-forming regions. [ABSTRACT FROM AUTHOR]

Context: EXor objects are young variables that show episodic variations of brightness commonly associated to enhanced accretion outbursts. Aims: With the aim of investigating the long-term photometric behaviour of a few EXor sources, we present here data from the archival plates of the Asiago observatory, showing the Orion field where the three EXors V1118, V1143, and NY are located. Methods: A total of 484 plates were investigated, providing a total of more than 1000 mag for the three stars, which cover a period of about 35 yr between 1959 to 1993. We then compared our data with literature data. Results: Apart from a newly discovered flare-up of V1118, we identify the same outbursts already known, but we provide two added values: (i) a long-term sampling of the quiescence phase; and (ii) repeated multi-colour observations (BVRI bands). The former allows us to give a reliable characterisation of the quiescence, which represents a unique reference for studies that will analyze future outbursts and the physical changes induced by these events. The latter is useful for confirming whether the intermittent increases of brightness are accretion-driven (as in the case of V1118), or extinction-driven (as in the case of V1143). Accordingly, doubts arise about the V1143 classification as a pure EXor object. Finally, although our plates do not separate NY Ori and the star very close to it, they indicate that this EXor did not undergo any major outbursts during our 40 yr of monitoring. [ABSTRACT FROM AUTHOR]

Aims: The aim of this work is to constrain properties of the disk around the archetype FU Orionis object, FU Ori, with as good as ~25 au resolution. Methods. We resolved FU Ori at 29-37 GHz using the Karl G. Jansky Very Large Array (JVLA) in the A-array configuration, which provided the highest possible angular resolution to date at this frequency band (~007). We also performed complementary JVLA 8-10 GHz observations, Submillimeter Array (SMA) 224 GHz and 272 GHz observations, and compared these with archival Atacama Large Millimeter Array (ALMA) 346 GHz observations to obtain the spectral energy distributions (SEDs). Results: Our 8-10 GHz observations do not find evidence for the presence of thermal radio jets, and constrain the radio jet/wind flux to at least 90 times lower than the expected value from the previously reported bolometric luminosity-radio luminosity correlation. The emission at frequencies higher than 29 GHz may be dominated by the two spatially unresolved sources, which are located immediately around FU Ori and its companion FU Ori S, respectively. Their deconvolved radii at 33 GHz are only a few au, which is two orders of magnitude smaller in linear scale than the gaseous disk revealed by the previous Subaru-HiCIAO 1.6 µm coronagraphic polarization imaging observations. We are struck by the fact that these two spatially compact sources contribute to over 50% of the observed fluxes at 224 GHz, 272 GHz, and 346 GHz. The 8-346 GHz SEDs of FU Ori and FU Ori S cannot be fit by constant spectral indices (over frequency), although we cannot rule out that it is due to the time variability of their (sub)millimeter fluxes. Conclusions: The more sophisticated models for SEDs considering the details of the observed spectral indices in the millimeter bands suggest that the >29 GHz emission is contributed by a combination of free-free emission from ionized gas and thermal emission from optically thick and optically thin dust components. We hypothesize that dust in the innermost parts of the disks (≤0.1 au) has been sublimated, and thus the disks are no longer well shielded against the ionizing photons. The estimated overall gas and dust mass based on SED modeling, can be as high as a fraction of a solar mass, which is adequate for developing disk gravitational instability. Our present explanation for the observational data is that the massive inflow of gas and dust due to disk gravitational instability or interaction with a companion/intruder, was piled up at the few-au scale due to the development of a deadzone with negligible ionization. The piled up material subsequently triggered the thermal instability and the magnetorotational instability when the ionization fraction in the inner sub-au scale region exceeded a threshold value, leading to the high protostellar accretion rate. [ABSTRACT FROM AUTHOR]

Context. Pre-main sequence stars are variable sources. The main mechanisms responsible for their variability are variable extinction, unsteady accretion, and rotational modulation of both hot and dark photospheric spots and X-ray-active regions. In stars with disks, this variability is related to the morphology of the inner circumstellar region (≤0.1AU) and that of the photosphere and corona, all impossible to be spatially resolved with present-day techniques. This has been the main motivation for the Coordinated Synoptic Investigation of NGC 2264, a set of simultaneous observations of NGC 2264 with 15 different telescopes. Aims. In this paper, we focus on the stars with disks. We analyze the X-ray spectral properties extracted during optical bursts and dips in order to unveil the nature of these phenomena. Stars without disks are studied in a companion paper. Methods. We analyze simultaneous CoRoT and Chandra/ACIS-I observations to search for coherent optical and X-ray flux variability in stars with disks. Then, stars are analyzed in two different samples. In stars with variable extinction, we look for a simultaneous increase of optical extinction and X-ray absorption during the optical dips; in stars with accretion bursts, we search for soft X-ray emission and increasing X-ray absorption during the bursts. Results. We find evidence for coherent optical and X-ray flux variability among the stars with variable extinction. In 9 of the 24 stars with optical dips, we observe a simultaneous increase of X-ray absorption and optical extinction. In seven dips, it is possible to calculate the NH/AVV ratio in order to infer the composition of the obscuring material. In 5 of the 20 stars with optical accretion bursts, we observe increasing soft X-ray emission during the bursts that we associate to the emission of accreting gas. It is not surprising that these properties are not observed in all the stars with dips and bursts, since favorable geometric configurations are required. Conclusions. The observed variable absorption during the dips is mainly due to dust-free material in accretion streams. In stars with accretion bursts, we observe, on average, a larger soft X-ray spectral component not observed in non-accreting stars. [ABSTRACT FROM AUTHOR]

We assess the evolutionary status of EXors. These low-mass, pre-main-sequence stars repeatedly undergo sharp luminosity increases, each a year or so in duration. We place into the HR diagram all EXors that have documented quiescent luminosities and effective temperatures, and thus determine their masses and ages. Two alternate sets of pre-main-sequence tracks are used, and yield similar results. Roughly half of EXors are embedded objects, i.e., they appear observationally as Class I or flat-spectrum infrared sources. We find that these are relatively young and are located close to the stellar birthline in the HR diagram. Optically visible EXors, on the other hand, are situated well below the birthline. They have ages of several Myr, typical of classical T Tauri stars. Judging from the limited data at hand, we find no evidence that binarity companions trigger EXor eruptions; this issue merits further investigation. We draw several general conclusions. First, repetitive luminosity outbursts do not occur in all pre-main-sequence stars, and are not in themselves a sign of extreme youth. They persist, along with other signs of activity, in a relatively small subset of these objects. Second, the very existence of embedded EXors demonstrates that at least some Class I infrared sources are not true protostars, but very young pre-main-sequence objects still enshrouded in dusty gas. Finally, we believe that the embedded pre-main-sequence phase is of observational and theoretical significance, and should be included in a more complete account of early stellar evolution. [ABSTRACT FROM AUTHOR]

The protoplanetary system HD 169142 is one of the few cases where a potential candidate protoplanet has recently been detected by direct imaging in the near-infrared. To study the interaction between the protoplanet and the disk itself, observations of the gas and dust surface density structure are needed. This paper reports new ALMA observations of the dust continuum at 1.3 mm, 12CO, 13CO, and C18O J = 2-1 emission from the system HD 169142 (which is observed almost face-on) at an angular resolution of ~ (~35 × 20 au). The dust continuum emission reveals a double-ring structure with an inner ring between (~20-35 au) and an outer ring between (~56-83 au). The size and position of the inner ring is in good agreement with previous polarimetric observations in the near-infrared and is consistent with dust trapping by a massive planet. No dust emission is detected inside the inner dust cavity (R ܐ 20 au) or within the dust gap (~35-56 au) down to the noise level. In contrast, the channel maps of the J = 2-1 line of the three CO isotopologs reveal gas inside the dust cavity and dust gap. The gaseous disk is also much larger than the compact dust emission; it extends to ~1 (~180 au) in radius. This difference and the sharp drop of the continuum emission at large radii point to radial drift of large dust grains (>µm size). Using the thermo-chemical disk code dali, we modeled the continuum and the CO isotopolog emission to quantitatively measure the gas and dust surface densities. The resulting gas surface density is reduced by a factor of ~30-40 inward of the dust gap. The gas and dust distribution indicate that two giant planets shape the disk structure through dynamical clearing (dust cavity and gap) and dust trapping (double-ring dust distribution). [ABSTRACT FROM AUTHOR]

RW Aur is a young binary star that experienced a deep dimming in 2010-2011 in component A and a second even deeper dimming from summer 2014 to summer 2016. We present new unresolved multiband photometry during the 2014-2016 eclipse, new emission line spectroscopy before and during the dimming, archive infrared photometry between 2014 and 2015, as well as an overview of literature data. Spectral observations were carried out with the Fibre-fed RObotic Dual-beam Optical Spectrograph on the Liverpool Telescope. Photometric monitoring was done with the Las Cumbres Observatory Global Telescope Network and James Gregory Telescope. Our photometry shows that RWAur dropped in brightness to R = 12.5 in 2016 March. In addition to the long-term dimming trend, RW Aur is variable on time-scales as short as hours. The short-term variation is most likely due to an unstable accretion flow. This, combined with the presence of accretion-related emission lines in the spectra suggest that accretion flows in the binary system are at least partially visible during the eclipse. The equivalent width of [O I] increases by a factor of 10 in 2014, coinciding with the dimming event, confirming previous reports. The blueshifted part of the Ha profile is suppressed during the eclipse. In combination with the increase in mid-infrared brightness during the eclipse reported in the literature and seen in WISE archival data, and constraints on the geometry of the disc around RW Aur A we arrive at the conclusion that the obscuring screen is part of a wind emanating from the inner disc. [ABSTRACT FROM AUTHOR]

The astronomy literature contains a wealth of photometric measurements of the brightness variations of weak-lined T Tauri stars. Comparing observations through two different filters, the brightness changes of the vast majority of these stars trace straight lines in magnitude- magnitude plots. It is plausible that the slope of the line is determined by a fixed mean starspot temperature, and that changing filling factors cause the linear magnitude-magnitude relations to be described. This is exploited to derive an estimator for the starspot temperature in terms of a linear regression slope, valid for modest starspot filling factors. Simulations are used to demonstrate that such regression slopes first need to be corrected for measurement errors, in order to avoid biased results. The theory is applied to a collection of 45 extensive sets of BVR (and in some cases UBVR) observations of weak-lined T Tauri stars taken from the literature. The results are examined critically, and it is pointed out that there are systematic differences between spot temperatures determined respectively from BV and VR data pairs. The reason for this is not known; possible causes are briefly speculated about. There is generally poor agreement with published spot temperatures, which also vary considerably. The simplest explanation is that the mean temperature of starspots is variable. [ABSTRACT FROM AUTHOR]

We present long-term photometric observations of the young open cluster IC 348 with a baseline time-scale of 2.4 yr. Our study was conducted with several telescopes from the Young Exoplanet Transit Initiative (YETI) network in the Bessel R band to find periodic variability of young stars. We identified 87 stars in IC 348 to be periodically variable; 33 of them were unreported before. Additionally, we detected 61 periodic non-members of which 41 are new discoveries. Our wide field of view was the key to those numerous newly found variable stars. The distribution of rotation periods in IC 348 has always been of special interest. We investigate it further with our newly detected periods but we cannot find a statistically significant bimodality. We also report the detection of a close eclipsing binary in IC 348 composed of a low-mass stellar component (M ≳ 0.09 M⊙) and a K0 pre-main-sequence star (M ≈ 2.7 M⊙). Furthermore, we discovered three detached binaries among the background stars in our field of view and confirmed the period of a fourth one. [ABSTRACT FROM AUTHOR]

Spatially resolved structures in protoplanetary disks hint at unseen planets. Previous imaging observations of the transitional disk around MWC 758 revealed an inner cavity, a ring-like outer disk, emission clumps, and spiral arms, all possibly generated by companions. We present ALMA dust continuum observations of MWC 758 at 0.87 mm wavelength with 43 × 39 mas angular resolution (6.9 × 6.2 au) and 20 μ Jy beam ^−1 rms. The central submillimeter emission cavity is revealed to be eccentric; once deprojected, its outer edge can be well fitted by an ellipse with an eccentricity of 0.1 and one focus on the star. The broad ring-like outer disk is resolved into three narrow rings with two gaps in between. The outer two rings tentatively show the same eccentricity and orientation as the innermost ring bounding the inner cavity. The two previously known dust emission clumps are resolved in both the radial and azimuthal directions, with radial widths equal to ∼4× the local scale height. Only one of the two spiral arms previously imaged in near-infrared (NIR) scattered light is revealed in ALMA dust emission, at a slightly larger stellocentric distance owing to projection effects. We also submit evidence of disk truncation at ∼100 au based on comparing NIR imaging observations with models. The spirals, the north clump, and the truncated disk edge are all broadly consistent with the presence of one companion exterior to the spirals at roughly 100 au.

We report high spatial resolution i' band imaging of the multiple T Tauri system LkHα 262/LkHα 263 obtained during the first commissioning period of the Adaptive Optics Lucky Imager (AOLI) at the 4.2-m William Herschel Telescope, using its Lucky Imaging mode. AOLI images have provided photometry for each of the two components LkHα 263 A and B (0.41 arcsec separation) and marginal evidence for an unresolved binary or a disc in LkHα 262. The AOLI data combined with previously available and newly obtained optical and infrared imaging show that the three components of LkHα 263 are comoving, that there is orbital motion in the AB pair, and, remarkably, that LkHα 262-263 is a common proper motion system with less than 1 mas yr-1 relative motion. We argue that this is a likely five-component gravitationally bounded system. According to BT-settl models the mass of each of the five components is close to 0.4 M⊙ and the age is in the range 1-2 Myr. The presence of discs in some of the components offers an interesting opportunity to investigate the formation and evolution of discs in the early stages of multiple very low mass systems. In particular, we provide tentative evidence that the disc in 263C could be coplanar with the orbit of 263AB. [ABSTRACT FROM AUTHOR]

Given that a majority of stars form in multiple systems, in order to fully understand the star- and planet-formation processes we must seek to understand them in multiple stellar systems. With this in mind, we present an analysis of the enigmatic binary T-Tauri system VV Corona Australis, in which both components host discs, but only one is visible at optical wavelengths. We seek to understand the peculiarities of this system by searching for a model for the binary which explains all the available continuum observations of the system. We present new mid-infrared interferometry and near-infrared (NIR) spectroscopy along with archival millimetre-wave observations, which resolve the binary at 1.3 mm for the first time. We compute a grid of pre-main-sequence radiative transfer models and calculate their posterior probabilities given the observed spectral energy distributions and mid-infrared interferometric visibilities of the binary components, beginning with the assumption that the only differences between the two components are their inclination and position angles. Our best-fitting solution corresponds to a relatively low-luminosity T-tauri binary, with each component's disc having a large scaleheight and viewed at moderate inclination (~50°), with the infrared companion inclined by ~5° more than the primary. Comparing the results of our model to evolutionary models suggests stellar masses ~1.7 M⊙ and an age for the system of 3.5 Myr, towards the upper end of previous estimates. Combining these results with accretion indicators from NIR spectroscopy, we determine an accretion rate of 4.0 × 10-8 M⊙ yr-1 for the primary. We suggest that future observations of VV Corona Australis and similar systems should prioritize high angular resolution sub-mm and NIR imaging of the discs and high-resolution optical/NIR spectroscopy of the central stars. [ABSTRACT FROM AUTHOR]

Results of the time variability monitoring of the two classical T Tauri stars, RU Lup and IM Lup, are presented. Three photometric data sets were utilized: (1) simultaneous (same field) MOST satellite observations over four weeks in each of the years 2012 and 2013, (2) multicolour observations at the South African Astronomical Observatory in April-May of 2013, (3) archival V-filter All Sky Automated Survey (ASAS) data for nine seasons, 2001- 2009. They were augmented by an analysis of high-resolution, public-domain VLT-UT2 Ultraviolet Visual Echelle Spectrograph spectra from the years 2000 to 2012. From the MOST observations, we infer that irregular light variations of RU Lup are caused by stochastic variability of hotspots induced by unstable accretion. In contrast, the MOST light curves of IM Lup are fairly regular and modulated with a period of about 7.19-7.58 d, which is in accord with ASAS observations showing a well-defined 7.247 ± 0.026 d periodicity. We propose that this is the rotational period of IM Lup and is due to the changing visibility of two antipodal hotspots created near the stellar magnetic poles during the stable process of accretion. Reanalysis of RU Lup high-resolution spectra with the broadening function approach reveals signs of a large polar coldspot, which is fairly stable over 13 years. As the star rotates, the spot-induced depression of intensity in the broadening function profiles changes cyclically with period 3.710 58 d, which was previously found by the spectral cross-correlation method. [ABSTRACT FROM AUTHOR]