T Tauri Stars Research Articles

Measuring the Spot Variability of T Tauri Stars Using Near-infrared Atomic Fe and Molecular OH Lines

As part of the Young Exoplanets Spectroscopic Survey, this study explores the spot variability of 13 T Tauri Stars (TTSs) in the near-infrared H band, using spectra from the Immersion GRating INfrared Spectrometer. By analyzing effective temperature (T eff) sensitive lines of atomic Fe i at ∼1.56259 μm and ∼1.56362 μm, and molecular OH at ∼1.56310 and ∼1.56317 μm, we develop an empirical equivalent width ratio (EWR) relationship for T eff in the range of 3400–5000 K. This relationship allows for precise relative T eff estimates to within tens of Kelvin and demonstrates compatibility with solar metallicity target models. However, discrepancies between observational data and model predictions limit the extension of the T eff–EWR relationship to a broader parameter space. Our study reveals that both classical and weak-line TTSs can exhibit T eff variations exceeding 150 K over a span of 2 yr. The detection of a quarter-phase delay between the EWR and radial velocity phase curves in TTSs indicates spot-driven signals. A phase delay of 0.06 ± 0.13 for CI Tau, however, suggests additional dynamics, potentially caused by planetary interaction, inferred from a posited 1:1 commensurability between the rotation period and orbital period. Moreover, a positive correlation between T eff variation amplitude and stellar inclination angle supports the existence of high-latitude spots on TTSs, further enriching our understanding of stellar surface activity in young stars.

A Multiwavelength, Multiepoch Monitoring Campaign of Accretion Variability in T Tauri Stars from the ODYSSEUS Survey. III. Optical Spectra* * Based on observations collected at the European Southern Observatory under ESO program 106.20Z8.

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

Simultaneous Far-ultraviolet and Near-ultraviolet Observations of T Tauri Stars with UVIT/AstroSat: Probing the Accretion Process in Young Stars

We present results from simultaneous far-ultraviolet (FUV) and near-ultraviolet (NUV) observations of T Tauri stars (TTSs) in the Taurus molecular cloud with UVIT/AstroSat. This is the very first UVIT study of TTSs. From the spectral energy distribution of TTSs from FUV to IR, we show that classical TTSs (CTTSs) emit significantly higher UV excess compared to weak-line TTSs (WTTSs). The equivalent blackbody temperatures corresponding to the UV excess in CTTSs (>104 K) are also found to be relatively higher than those in WTTSs (<9250 K). From the UV excess, we have reclassified two WTTSs (BS Tau and V836 Tau) as CTTSs, which has been supported by the follow-up optical spectroscopic study using the Himalayan Chandra Telescope, showing strong Hα line emission. We find that CTTSs show strong excess emission in both the FUV (>107) and NUV (>103) bands, while WTTSs show strong excess only in the FUV (≲105), suggesting that excess emission in the NUV can be used as a tool to classify the TTSs. We also find a linear correlation between UV luminosity (a primary indicator of mass accretion) and Hα luminosity (a secondary indicator of mass accretion) with a slope of 1.20 ± 0.22 and intercept of 2.16 ± 0.70.

Magnetic field, magnetospheric accretion and candidate planet of the young star GM Aurigae observed with SPIRou

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

A Multiwavelength, Multiepoch Monitoring Campaign of Accretion Variability in T Tauri Stars from the ODYSSEUS Survey. II. Photometric Light Curves

Classical T Tauri stars (CTTSs) are young, low-mass stars that accrete material from their surrounding protoplanetary disk. To better understand accretion variability, we conducted a multiepoch, multiwavelength photometric monitoring campaign of four CTTSs, TW Hya, RU Lup, BP Tau, and GM Aur, in 2021 and 2022, contemporaneous with Hubble Space Telescope UV and optical spectra. We find that all four targets display significant variability in their light curves, generally on days-long timescales (but in some cases year-to-year), often due to periodicity associated with stellar rotation and to stochastic accretion variability. There is a strong connection between mass accretion and photometric variability in all bands, but the relationship varies per target and epoch. Thus, photometry should be used with caution as a direct measure of accretion in CTTSs.

A Multiwavelength, Multiepoch Monitoring Campaign of Accretion Variability in T Tauri Stars from the ODYSSEUS Survey. I. HST Far-UV and Near-UV Spectra

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

Long-period modulation of the classical T Tauri star CI Tau

Context. Detecting planets within protoplanetary disks around young stars is essential for understanding planet formation and evolution. However, planet detection using the radial velocity method faces challenges due to the strong stellar activity in these early stages. Aims. We detect long-term periodicities in photometric and spectroscopic time series of the classical T Tauri star (CTTS) CI Tau, and retrieve evidence for inner embedded planets in its disk. Methods. The study conducted photometric and spectroscopic analyses using K2 and Las Cumbres Observatory Global Network light curves, and high-resolution spectra from ESPaDOnS and SPIRou. We focused our radial velocity analysis on a wavelength domain less affected by spot activity. To account for spot effects, a quasi-periodic Gaussian process model was applied to K2 light curve, ESPaDOnS, and SPIRou radial velocity data. Additionally, a detailed bisector analysis on cross-correlation functions was carried out to understand the cause of long-term periodicity. Results. We detect coherent periods at ~6.6d, 9d, ~11.5d, ~14.2d, and ~25.2d, the last of which is seen consistently across all datasets. Bisector analysis of the cross-correlation functions provides strong hints for combined activity-induced and Doppler reflex signals in the radial velocities at a period of 25.2 d. Our analysis suggests that this periodicity is best explained by the presence of a 3.6 ± 0.3 MJup eccentric (e ~ 0.58) planet at a semi-major axis of 0.17 au. Conclusions. We report the detection of a massive inner planet in CI Tau. Our study outlines the difficulty of searching for disk-embedded planets in the inner 0.1 au of young and active systems. When searching for planets in actively accreting stars such as CI Tau, we demonstrate that the primary limitation is stellar activity rather than the precision of RV measurements provided by the instrument.

Source of the Fe Kalpha emission in T Tauri stars. Radiation induced by relativistic electrons during flares. An application to RY Tau

T Tauri stars (TTSs) are magnetically active stars that accrete matter from the inner border of the surrounding accretion disk; plasma becomes trapped into the large-scale magnetic structures and falls onto the star, heating the surface through the so-called accretion shocks. The X-ray spectra of the TTSs show prominent Fe Kalpha fluorescence emission at 6.4 keV (hereafter, Fe Kalpha emission) that cannot be explained in a pure accretion scenario because its excitation requires significantly more energy than the maximum available through the well-constrained free-fall velocity. Neither can it be produced by the hot coronal plasma. TTSs display all signs of magnetic activity, and magnetic reconnection events are expected to occur frequently. In these events, electrons may become accelerated to relativistic speeds, and their interaction with the environmental matter may result in Fe Kalpha emission. It is the aim of this work to evaluate the expected Fe Kalpha emission in the context of the TTS research and compare it with the actual Fe Kalpha measurements obtained during the flare detected while monitoring RY Tau with the XMM-Newton satellite. The propagation of high-energy electrons in dense gas generates a cascade of secondary particles that results in an electron shower of random nature, whose evolution and radiative throughput was simulated in this work using the Monte Carlo code PENELOPE. A set of conditions representing the environment of the TTSs where these showers may impinge was taken into account to generate a grid of models that can aid the interpretation of the data. The simulations show that the electron beams produce a hot spot at the point of impact; strong Fe Kalpha emission and X-ray continuum radiation are produced by the spot. This emission is compatible with RY Tau observations. The Fe Kalpha emission observed in TTSs could be produced by beams of relativistic electrons accelerated in magnetic reconnection events during flares.

The Effect of Starspots on Spectroscopic Age and Mass Estimates of Nonaccreting T Tauri Stars in the Taurus–Auriga Star-forming Region

Accurate age and mass determinations for young pre-main-sequence stars are made challenging by the presence of large-scale starspots. We present results from a near-infrared spectroscopic survey of 10 T-Tauri Stars in Taurus-Auriga that characterize spot-filling factors and temperatures, the resulting effects on temperature and luminosity determinations, and the consequences for inferred stellar masses and ages. We constructed composite models of spotted stars by combining BTSettl-CIFIST synthetic spectra of atmospheres to represent the spots and the photosphere along with continuum emission from a warm inner disk. Using a Markov Chain Monte Carlo algorithm, we find the best-fit spot and photospheric temperatures, spot-filling factors, as well as disk-filling factors. This methodology allowed us to reproduce the 0.75–2.40 μm stellar spectra and molecular feature strengths for all of our targets, disentangling the complicated multicomponent emission. For a subset of stars with multiepoch observations spanning an entire stellar rotation, we correlate the spectral variability and changes in the filling factors with rotational periods observed in K2 and AAVSO photometry. Combining spot-corrected effective temperatures and Gaia distances, we calculate luminosities and use the Stellar Parameters of Tracks with Starspots models to infer spot-corrected masses and ages for our sample of stars. Our method of accounting for spots results in an average increase of 60% in mass and a doubling in age with respect to traditional methods using optical spectra that do not account for the effect of spots.

Chandra X-ray analysis of herbig Ae/Be stars

Abstract Herbig Ae/Be (HAeBe) stars are intermediate-mass pre-main sequence stars, characterized by infrared excess and emission lines. They are observed to emit X-rays, whose origin is a matter of discussion and not settled yet. X-ray emission is not expected in HAeBe stars, as they lack the sub-surface convective zone. In this study, we retrieved observations from the Chandra archive for 62 HAeBe stars, among which 44 sources (detection fraction ∼71 per cent) were detected in X-rays, with 7 being new detections. We use this sample as a test bed to conduct a comparative analysis of the X-ray properties of HAeBe stars and their low-mass counterparts, T Tauri Stars (TTSs). Further, we compare the X-ray properties of HAeBe stars and TTSs with optical and IR properties to constrain the X-ray emission mechanism in HAeBe stars. We found no correlation between X-ray emission and disk properties of HAeBe stars, confirming that X-rays are not related to accretion shocks. About 56 % of HAeBe stars without any known sub-arcsec companions have lower plasma temperatures (kT ≤ 2 keV). We observe flaring/variability in HAeBe stars with confirmed low-mass companions. These stars show plasma temperatures &amp;gt; 2 keV, similar to TTSs. Guided by this information we discuss the role of a T Tauri companion for X-ray emission seen in our sample of HAeBe stars. From the results obtained in this paper, we suggest that X-ray emission from HAeBe stars may not be related to accretion shocks or hidden TTS, but rather can be due to magnetically driven coronal emission.

GIARPS High-resolution Observations of T Tauri stars (GHOsT)

Aims. This paper aims to revisit the kinematical and physical properties of the warm (T ~ 5000–10 000 K) atomic gas in the inner disk (&lt;5 au) region of classical T Tauri stars (CTTs) and relate them to the properties of the outer dusty disk resolved with ALMA. We also want to define constraints for the mass-loss in the inner atomic winds and jets to assess their role in the evolution and dispersal of planet-forming disks. Methods. We used the high resolution (R = 115 000, ~2.6 km s−1) spectra of 36 CTTs observed as part of the GIARPS High-resolution Observations of T Tauri stars (GHOsT) project and analysed the profile and luminosity of the brightest optical forbidden lines, namely [O I] 630 nm, [O I] 557 nm, [S II] 406 nm, [S II] 673 nm, and [N II] 658 nm. Results. We decomposed the line profiles into different velocity components, and concentrated our analysis mostly on the so-called narrow low-velocity component (NLVC). We find that about 40% of sources display a NLVC peak velocity (Vp) compatible with the stellar velocity. These include the transitional disks (TD) and typically show a single low velocity component (LVC), lower mass accretion rates, and the absence of a jet. They therefore might represent later evolutionary stages where the emission from the disk is dominant with respect to the wind contribution. No difference in kinematical properties was instead found between sources with full disks and disks with substructures as resolved by ALMA. The [O I] 630 nm profiles peaking at the stellar velocity are well fitted by a simple Keplerian disk model, where the emission line region extends from ~0.01 au up to several tens of au in some cases. The [O I] emission is detected inside the sub-millimetre dust cavities of all the TDs. No correlation is found between Rkep, derived from the line half width at half maximum (HWHM), and the size of the dust cavity. We see an anti-correlation between the [O I] 557/630 nm ratio and Rkep, which suggests that the [O I] emitting region expands as the gas dominating the emission cools and becomes less dense. We confirmed previous findings that the line ratios observed in the LVC, if compared with a thermal single temperature and density model, imply ne ~ 106–108 cm−3 and Te ~ 5000–10 000 K, and additionally constrained the ionisation fraction in the NLVC to be xe &lt; 0.1. We however discuss the limits of applying this diagnostic to winds that are not spatially resolved. Conclusions. The emission from the disk should be considered as an important contribution to the forbidden line emission in CTTs. Also, the clearing of warm atomic gas from the upper disk layers does not seem to follow the dispersal of the bulk of molecular gas and dust during late disk evolution. For the outflow component, we estimated the mass-loss for both the disk winds and jets. We conclude that without better knowledge of the wind geometry and spatial extent, and given the limitation of the diagnostics, the mass-loss rates in the wind traced by the blue shifted LVC cannot be constrained better than a factor of 100, with a Ṁwind/Ṁacc spanning between ~0.01 and more than 1. When compared with synthetic [O I] 630 nm images of X-ray photoevaporation models, the estimated Ṁwind represents a lower limit to the total mass-loss rate of the model, indicating that [O I] 630 nm is likely not the best tracer to probe mass-loss in low-velocity winds.

First detection and modeling of spatially resolved Lyαin TW Hya

Abstract Lyman-α (Lyα) is the strongest emission line in the accretion-generated UV spectra from T-Tauri stars and, as such, plays a critical role in regulating chemistry within the surrounding protoplanetary disks. Due to its resonant nature, the scattering of Lyαphotons along the line-of-sight encodes information about the physical properties of the intervening H i medium. In this work, we present the first spatially resolved spectral images of Lyαemission across a protoplanetary disk in the iconic face-on T-Tauri star TW Hya, observed with HST-STIS at spatial offsets 0″, ±0.2″, and ±0.4″. To comprehensively interpret these Lyαspectra, we utilize a 3D Monte-Carlo Lyαradiative transfer simulation considering the H iwind and protoplanetary disk. From the simulation, we constrain the wind’s properties: the H icolumn density ∼1020 cm−2 and the outflow velocity ∼200 km s−1. Our findings indicate that successfully interpreting the observed spectra necessitates scattering contributions in the H ilayer within the disk. Furthermore, to explore the effect of Lyαradiative transfer on protoplanetary disk chemistry, we compute the radiation field within the scattering medium and reveal that the wind reflection causes more Lyαphotons to penetrate the disk. Our results show the necessity of spatially resolved Lyαobservations of a broad range of targets, which will decode the complex interactions between the winds, protoplanetary disks, and surrounding environments.

Discovery of Two Different Full Disk Evolutionary Patterns of M-type T Tauri Stars with LAMOST DR8

The full disk, full of gas and dust, determines the upper limit of planet masses, and its lifetime is critical for planet formation, especially for giant planets. In this work, we studied the evolutionary timescales of the full disks of T Tauri stars (TTSs) and their relations to accretion. Combined with Gaia EDR3, Two Micron All Sky Survey, and Wide-field Infrared Survey Explorer data, 1077 disk-bearing TTS candidates were found in LAMOST DR8, and stellar parameters were obtained. Among them, 783 are newly classified by spectra as classical T Tauri stars (CTTSs; 169) or weak-lined T Tauri stars (WTTSs). Based on EW and FWHM of Hα, 157 TTSs in accretion were identified, with ∼82% also having full disks. For TTSs with M < 0.35M ☉, about 80% seem to already lose their full disks at ∼0.1 Myr, which may explain their lower mass, while the remaining 20% with full disks evolve at similar rates of non-full disks within 5 Myr, allowing enough time and material to form giant planets. The fraction of accreting TTSs to disk-bearing TTSs is stable at ∼10% and can last ∼5–10 Myr, suggesting that full disks and accretion evolve with similar rates as non-full disks. For TTSs with M > 0.35 M ☉, almost all full disks can survive more than 0.1 Myr, most for 1 Myr and some even for 20 Myr. For TTSs with M > 0.35 M ☉, almost all full disks can survive more than 0.1 Myr, most for 1 Myr, and some even for 20 Myr, which implies planets are more likely to be formed in their disks than those of M < 0.35 M ☉, and thus M dwarfs with M > 0.35 M ☉ can have more planets. The fraction of full-disk TTSs to disk-bearing TTSs decreases with age following the relation f ∝ t −0.35, and similar relations existed in the fraction of accreting TTSs and the fraction of full-disk CTTSs, suggesting faster full disks and accretion evolution than non-full disks. For full-disk stars, the ratio of accretion of lower-mass stars is systematically lower than that of higher-mass stars, confirming the dependence of accretion on stellar mass, which may be reflective of an observational bias in the detection of accretion levels, with the lower-mass stars crossing below the detection threshold earlier than higher-mass stars.

Identifying the population of T-Tauri stars in Taurus: UV–optical synergy

Identifying the population of T-Tauri stars in Taurus: UV–optical synergy

Radio-continuum decrements associated to shadowing from the central warp in transition disc DoAr 44

ABSTRACT Warps have often been used to explain disc properties, but well-characterized examples are important due to their role in disc evolution. Scattered light images of discs with central gaps have revealed sharp warps, such that the outer rings are shadowed by tilted inner discs. The near-infrared intensity drops along the ring around T-Tauri star DoAr 44 have been interpreted in terms of a central warp. We report new ALMA observations of DoAr 44 in the continuum at 230 and 350 GHz (at ∼10 au), along with a new epoch of Spectro Polarimetric High-contrast Exoplanet REsearch (SPHERE)/Infrared Dual-band Imager and Spectrograph differential polarized imaging taken during excellent weather conditions. The ALMA observations resolve the ring and confirm the decrements proposed from deconvolution of coarse 336 GHz data. The scattered light image constrains the dips, which correspond to a misaligned inner disc with a relative inclination ξ = 21.4 $^{+6.7}_{-8.3}$ deg. The SPHERE intensity profile shows a morphological change compared to a previous epoch that may be interpreted as a variable orientation of the inner disc, from ξ ∼ 30 deg to ξ ∼ 20 deg. The intensity dips probably correspond to temperature decrements, as their mm-spectral index, $\alpha ^{230 \textrm {GHz}}_{350 \textrm {GHz}} \sim$ 2.0 ± 0.1, is indicative of optically thick emission. The azimuth of the two temperature decrements are leading clockwise relative to the infrared-dips, by η = 14.95 deg and η = 7.92 deg. For a retrograde disc, such shifts are expected from a thermal lag and imply gas surface densities of Σg = 117 ± 10 g cm−2 and Σg = 48 ± 10 g cm−2. A lopsided disc, with contrast ratio fr = 2.4 ± 0.5, is also consistent with the large continuum crescent.

Mass Accretion, Spectral, and Photometric Properties of T Tauri Stars in Taurus Based on TESS and LAMOST

We present the analysis of 16 classical T Tauri stars (CTTSs) using LAMOST and TESS data, investigating spectral properties, photometric variations, and mass accretion rates. All 16 stars exhibit emissions in Hα lines, from which the average mass accretion rate of 1.76 × 10−9 M ⊙ yr−1 is derived. Two of the stars, DL Tau and Haro 6-13, show mass accretion bursts simultaneously in TESS, ASAS-SN, and/or the ZTF survey. Based on these observations, we find that the mass accretion rates of DL Tau and Haro 6-13 reach their maxima of 2.5 × 10−8 M ⊙ yr−1 and 2 × 10−10 M ⊙ yr−1, respectively, during the TESS observation. We detect 13 flares among these stars. The flare frequency distribution shows that the CTTSs’ flare activity is not only dominated by strong flares with high energy but also much more active than those of solar-type and young low-mass stars. By comparing the variability classes reported in the literature, we find that the transition timescale between different classes of variability in CTTSs, such as from stochastic (S) to bursting (B) or from quasi-periodic symmetric to quasi-periodic dipping, may range from 1.6 to 4 yr. We observe no significant correlation between inclination and mass accretion rates derived from the emission indicators. This suggests that inner disk properties may be more important than those of outer disks. Finally, we find a relatively significant positive correlation between the asymmetric metric M and the cold disk inclination compared to the literature. A weak negative correlation between the periodicity metric Q value and inclination has also been found.

T Tauri stars in the SuperWASP and NSVS surveys II. Spectral modelling

ABSTRACT We present results from long-term spectroscopic monitoring of 21 T-Tauri stars located in the Taurus–Auriga star-forming region (SFR). We combine medium and high-dispersion Echelle spectroscopy obtained at the Stará Lesná, Skalnaté Pleso (both in Slovakia), and Tautenburg (Germany) observatories with low-resolution flux-calibrated spectra from Asiago (Italy) observatory all taken between 2015 and 2018. We extend the coverage by additional medium-resolution spectra from Stará Lesná obtained in 2022. In the previous paper, we measured photometric periods of these targets in a range of 0.7–3.1 d, which could be due to the rotation of a spotted surface or binarity. Here, we use the broadening-function technique to determine the radial and projected rotational velocities to reveal any close binary companion. Our analysis concludes that no such companion is present with an orbital period equal to the photometric period. We focus our analysis primarily on determining atmospheric parameters such as surface gravity log g, effective temperature Teff, and metallicity [Fe/H]. Additionally, we measure the equivalent width of H α, Li i, and interstellar Na i lines. We also investigate the effect of possible reddening on individual targets and construct the HR diagram of our sample. Using pre-main-sequence evolutionary models, we determine the age of our targets. This analysis hints at ages younger than 50 Myr with mean age 5 ± 3 Myr, masses between 0.75 and 2.10 M⊙, and minimum radii in the range 0.60–3.17 R⊙. Altogether, the results are consistent with expected young stars with larger radii than those of main-sequence stars.

The impact of dust evolution on the dead zone outer edge in magnetized protoplanetary disks

Context. The dead zone outer edge corresponds to the transition from the magnetically dead to the magnetorotational instability (MRI) active regions in the outer protoplanetary disk midplane. It has been previously hypothesized to be a prime location for dust particle trapping. A more consistent approach to access such an idea has yet to be developed, since the interplay between dust evolution and MRI-driven accretion over millions of years has been poorly understood. Aims. We provide an important step toward a better understanding of the MRI–dust coevolution in protoplanetary disks. In this pilot study, we present a proof of concept that dust evolution ultimately plays a crucial role in the MRI activity. Methods. First, we study how a fixed power-law dust size distribution with varying parameters impacts the MRI activity, especially the steady-state MRI-driven accretion, by employing and improving our previous 1+1D MRI-driven turbulence model. Second, we relax the steady-state accretion assumption in this disk accretion model, and partially couple it to a dust evolution model in order to investigate how the evolution of dust (dynamics and grain growth processes combined) and MRI-driven accretion are intertwined on million-year timescales, from a more sophisticated modeling of the gas ionization degree. Results. Dust coagulation and settling lead to a higher gas ionization degree in the protoplanetary disk, resulting in stronger MRI-driven turbulence as well as a more compact dead zone. On the other hand, fragmentation has an opposite effect because it replenishes the disk in small dust particles which are very efficient at sweeping up free electrons and ions from the gas phase. Since the dust content of the disk decreases over millions of years of evolution due to radial drift, the MRI-driven turbulence overall becomes stronger and the dead zone more compact until the disk dust-gas mixture eventually behaves as a grain-free plasma. Furthermore, our results show that dust evolution alone does not lead to a complete reactivation of the dead zone. For typical T-Tauri stars, we find that the dead zone outer edge is expected to be located roughly between 10 au and 50 au during the disk lifetime for our choice of the magnetic field strength and configuration. Finally, the MRI activity evolution is expected to be crucially sensitive to the choice made for the minimum grain size of the dust distribution. Conclusions. The MRI activity evolution (hence the temporal evolution of the MRI-induced α parameter) is controlled by dust evolution and occurs on a timescale of local dust growth, as long as there are enough dust particles in the disk to dominate the recombination process for the ionization chemistry. Once that is no longer the case, the MRI activity evolution is expected to be controlled by gas evolution and occurs on a viscous evolution timescale.

Updated Planetary Mass Constraints of the Young V1298 Tau System Using MAROON-X

The early K-type T-Tauri star, V1298 Tau (V = 10 mag, age ≈ 20–30 Myr) hosts four transiting planets with radii ranging from 4.9 to 9.6 R ⊕. The three inner planets have orbital periods of ≈8–24 days while the outer planet’s period is poorly constrained by single transits observed with K2 and the Transiting Exoplanet Survey Satellite (TESS). Planets b, c, and d are proto–sub-Neptunes that may be undergoing significant mass loss. Depending on the stellar activity and planet masses, they are expected to evolve into super-Earths/sub-Neptunes that bound the radius valley. Here we present results of a joint transit and radial velocity (RV) modeling analysis, which includes recently obtained TESS photometry and MAROON-X RV measurements. Assuming circular orbits, we obtain a low-significance (≈2σ) RV detection of planet c, implying a mass of and a conservative 2σ upper limit of <39 M ⊕. For planets b and d, we derive 2σ upper limits of M b < 159 M ⊕ and M d < 41 M ⊕, respectively. For planet e, plausible discrete periods of P e > 55.4 days are ruled out at the 3σ level while seven solutions with 43.3 < P e/d < 55.4 are consistent with the most probable 46.768131 ± 000076 days solution within 3σ. Adopting the most probable solution yields a 2.6σ RV detection with a mass of 0.66 ± 0.26 M Jup. Comparing the updated mass and radius constraints with planetary evolution and interior structure models shows that planets b, d, and e are consistent with predictions for young gas-rich planets and that planet c is consistent with having a water-rich core with a substantial (∼5% by mass) H2 envelope.

Multiple Rings and Asymmetric Structures in the Disk of SR 21

Crescent-like asymmetric dust structures discovered in protoplanetary disks indicate dust aggregations. Thus, researching these structures helps us understand the planet formation process. Here we analyze ALMA data of the protoplanetary disk around the T-Tauri star SR 21, which has asymmetric structures that were detected in previous submillimeter observations. Imaged at ALMA Band 6 (1.3 mm) with a spatial resolution of about 0.″04, the disk is found to consist of two rings and three asymmetric structures, with two of the asymmetric structures being in the same ring. Compared to the Band 6 image, the Band 3 (2.7 mm) image also shows the three asymmetric structures, but with some clumps. The elongated asymmetric structures in the outer ring could be due to the interactions of a growing planet. Based on fitting the Band 3 and Band 6 dust continuum data, two branches of solutions for the maximum dust size in the disk are suggested: one is larger than 1 mm, and the other is smaller than 300 μm. High-resolution continuum observations at longer wavelengths as well as polarization observations can help break the degeneracy. We also suggest that the prominent spiral previously identified in VLT/SPHERE observations south of the star at 0.″25 may be the scattered-light counterpart of the inner arc, and the structure is a dust-trapping vortex in nature. The discovered features in SR 21 make it a good target for studying the evolution of asymmetric structures and planet formation.