Open Access
Issue
A&A
Volume 696, April 2025
Article Number A207
Number of page(s) 10
Section Planets, planetary systems, and small bodies
DOI https://doi.org/10.1051/0004-6361/202453207
Published online 25 April 2025

© The Authors 2025

Licence Creative CommonsOpen Access article, published by EDP Sciences, under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

This article is published in open access under the Subscribe to Open model. Subscribe to A&A to support open access publication.

1 Introduction

The lifetime and structure of a protoplanetary disk (PPD) set the tempo and conditions for planet formation. The gas disk drives the starward drift and concentration of solids and dampens planetesimal motion, accelerating planet formation (Birnstiel 2024). It exerts torques that cause the migration of nascent planets (Nelson 2018). Radial structures in a disk are important in this context (e.g., Lau et al. 2022). Observations of young cluster stars show that disk lifetime decreases with increasing stellar mass (Bayo et al. 2012; Ribas et al. 2015) but with an order-of-magnitude scatter (Pfalzner & Dincer 2024). There are exceptionally long-lived disks (>10 Myr) around some very low-mass stars (VLMSs; ≲0.2 M ; Murphy et al. 2018; Flaherty et al. 2019; Silverberg et al. 2020; Gaidos et al. 2022). The size of dust disks and the occurrence of disk structure (cavities or gaps) positively correlate with stellar mass (van der Marel & Mulders 2021; Kurtovic et al. 2021), although the latter could be in part a selection effect brought about by the former.

While disks can be truncated and their dissipation hastened by close companion stars (Kraus et al. 2012; Akeson et al. 2019; Zurlo et al. 2020, 2021) and far-ultraviolet radiation from neighboring O and B stars (Winter & Haworth 2022), these factors seem to only be important for a small fraction of stars (e.g, Zurlo et al. 2021; Mann et al. 2015). In the canonical scenario for single and wide binary stars, disk lifetime is thought to be set by a competition between viscous or magnetohydrodynamic (MHD) wind-driven accretion versus “internal” photoevaporation (PE) driven by high-energy emission from the central star (Armitage 2011), specifically soft X-rays (Ercolano & Pascucci 2017; Ercolano et al. 2021, 2023; Lin et al. 2025). When the accretion rate falls below the local PE rate, a gap appears; this occurs at a fraction of the gravitational radius (a few au, Liffman 2003). The flow of gas into the inner disk and onto the star then declines or halts (Drake et al. 2009) and a cavity develops. Depletion of dust in the inner disk leads to a concomitant decline in infrared emission that is the signature of “transition” disks (Ercolano & Pascucci 2017).

Incident stellar radiation on the PPD, including X-rays, can be attenuated by intervening structures in the inner disk. Many T Tauri stars (TTSs) exhibit optical variability due to occulting dust in accretion streams (Bouvier et al. 2003), winds (Tambovtseva & Grinin 2008), or instabilities near the corotation radius (Ansdell et al. 2016; Roggero et al. 2021). “See-saw” variability in infrared emission has been explained by shadowing of the disk by interior structures (Muzerolle et al. 2009; Espaillat et al. 2011; Flaherty et al. 2012; Fernandes et al. 2018; Gaidos et al. 2024).

For an interstellar-like composition (e.g., Güver & Özel 2009), the neutral hydrogen column density of the gas corresponding to the observed extinction would be 1021 × 1022 cm−2, but some measurements find much higher values (Güdel et al. 2005, 2007b; Robrade & Schmitt 2007; Schneider et al. 2015; Günther et al. 2018). This intervening gas can absorb X-rays from the star that would otherwise reach the disk. Accreting TTSs have systematically lower observed X-ray emission than their non-accreting coeval counterparts due to circumstellar gas, an effect that is correlated with the accretion rate (Telleschi et al. 2007; Güdel et al. 2007a; Bustamante et al. 2016). Some of the most extreme accretors (e.g., FU Orionis-type objects) host dust- depleted gas that is optically thick to 0.1–1 keV X-rays (Güdel et al. 2008; Skinner et al. 2010; Liebhart et al. 2014). Attenuation along shallower lines of sight to the disk surface (elevations of a few degrees) could be even higher.

Güdel et al. (2010) proposed that X-ray-absorbing gas in the immediate environment of accreting TTSs means that only ≳1 keV X-rays, not more readily absorbed extreme-ultraviolet (EUV) radiation, reach the outer disk. Pascucci et al. (2020) suggest that absorption of soft X-rays by a magnetized inner disk wind (IDW) at 1 au) could suppress any PE wind at ≳ a few astronomical units. The appearance of an IDW can, in turn, depend sensitively on stellar rotation, the large-scale stellar magnetic field, and inner disk accretion (Romanova & Owocki 2015). This motivates the question of whether this link between conditions around the central star and PE-triggered disk dissipation could explain the disparate fates of disks and trends of the disk lifetime and structure. Here, we combine an analytical model of disk evolution with models of pre-main-sequence (pre-MS) stars to investigate this possibility.

2 Model

In the inner disk, gas at >1000 K is partially ionized and influenced by the large-scale stellar magnetic field. Interior to the corotation radius Rcor=(GMP24π2)1/3,R_{\rm cor} = \left(\frac{GM_{\star}P_{\star}^2}{4\pi^2}\right)^{1/3},(1) where G is the gravitational constant, M* is the stellar mass, and P* is the rotation period, the field exerts a negative torque on any disk gas, driving accretion onto the star. Beyond the corotation radius, the torque is positive and disk gas can be accelerated and ejected along magnetic field lines. An ionized disk is expected to be truncated where magnetic pressure exceeds the ram pressure of inflowing gas. Three-dimensional MHD simulations (Takasao et al. 2022; Zhu et al. 2024) support a relation between the truncation radius, Rmag, the strength of the large-scale dipole magnetic field at the stellar surface, B*, and the disk accretion rate, M˙acc$\dot{M}_{\rm acc}$, originally derived for spherical flow onto a dipole field (Ghosh & Lamb 1979): Rmag=(B4R122GMM˙acc2)1/7.R_{\rm mag} = \left(\frac{B_{\star}^4R_{\star}^{12}}{2GM_{\star}\dot{M}_{\rm acc}^2}\right)^{1/7}.(2) Rmag is also where any magnetic field advected inward by disk accretion will accumulate, increasing magnetic pressure.

The structure of the innermost disk and flow in it depends sensitively on the ratio RRmag/Rcor$\rr \equiv R_{\rm mag}/R_{\rm cor}$. If B* is in kilogauss, R* and M* are in solar units, P* is in days, and M˙acc$\dot{M}_{\rm acc}$ is in units of 10−9 M yr−1, R=1.8B4/7R12/7M˙acc2/7M10/21P2/3.\mathcal{R} = 1.8 B_{\star}^{4/7} R_{\star}^{12/7} \dot{M}_{\rm acc}^{-2/7} M_{\star}^{-10/21} P_{\star}^{-2/3}(3)

If R<0.8$\mathcal{R} < 0.8$ then magnetospheric accretion occurs in an unsteady flow that remains close to the disk midplane. At intermediate R$\mathcal{R}$ values of 0.8–1, magnetospheric accretion takes the form of a highly non-axisymmetric “funnel” flow out of the midplane and along field lines onto the star (Bouvier et al. 2007). However, if R1$\mathcal{R} \ge 1$, “propeller accretion” occurs; some disk gas is accelerated in a magneto-centripetal IDW, and the remainder is accreted onto the star (Romanova & Owocki 2015), with the ejection/accretion ratio, f , depending on loading of disk gas onto stellar field lines and the wind’s ratio of Alfven to launch radius. Evidence for such winds includes highly broadened forbidden lines of O I (Fang et al. 2023). The critical M˙acc$\dot{M}_{\rm acc}$ at which Rmag and Rcor intersect (R=1$\mathcal{R} = 1$), in other words the transition between magnetospheric and propeller accretion, is: M˙crit=8.2×109B2R6M5/3P7/3Myr1.\dot{M}_{\rm crit} = 8.2 \times 10^{-9} B_{\star}^2 R_{\star}^6 M_{\star}^{-5/3} P_{\star}^{-7/3} M_{\odot} {\rm yr}^{-1}.(4)

The X-ray emission from the star is calculated using the standard broken power law dependent on the Rossby number, Ro = P*c (Wright et al. 2011), where τc, the convective turnover time in the stellar interior, is taken to be 12.6L1/2$12.6\, L_{\star}^{-1/2}$ days and L* is the stellar luminosity in solar units (Jeffries et al. 2011). LX=7.4×104L1+(Ro/0.13)2.18.L_{\rm X} = \frac{7.4 \times 10^{-4} L_{\star}}{1 + \left(Ro/0.13\right)^{2.18}}.(5)

A geometrically thin wind1 launched from both sides of a disk at Rmag and an angle, θ, with respect to the disk normal will impose an atomic column density, NH, between the star and disk of NHfM˙acc4πRmagcosθμγmpkBT,N_{\rm H} \sim \frac{f \dot{M}_{\rm acc}}{4\pi R_{\rm mag} \cos \theta \sqrt{\mu \gamma m_p k_B T}},(6) where we assume that lines of sight to the disk pass through the transonic flow near the base of the wind, and where μ = 2.38, γ = 1.4, and T are, respectively, the molecular weight, adiabatic constant, and temperature of the flow. W✓e take T to be that of a passively heated inner disk wall: TR/Rtrunc${\approx} T_{\star}\sqrt{R_{\star}/R_{\rm trunc}}$ (Dullemond et al. 2001). The ejection/accretion ratio, f , is taken to be 0.2, near typical values inferred by Watson et al. (2016) and Serna et al. (2024). This yields NH4×1021cm2M˙917/14R43/28M3/28B3/7(T/4000)1/4.N_H \approx 4 \times 10^{21} {\rm cm}^{-2} \dot{M}_{-9}^{17/14} R_{\star}^{-43/28} M_{\star}^{3/28}B_{\star}^{-3/7} (T_{\star}/4000)^{1/4}.(7)

We calculated the photoelectric absorption of 0.1–1 keV X-rays by intervening solar-metallicity gas with NH using the XSPEC package (Arnaud 1996). We modeled coronal X-ray emission to be a two-temperature (0.5 and 2 keV), equal emission measure plasma like that typically used to describe TTSs (e.g., Telleschi et al. 2007), and we adopted coronal abundances found in Taurus stars (Güdel & Telleschi 2007). Unit optical depth is reached at NH 1021 cm−2, weakly dependent on the coronal X-ray spectrum (Appendix B).

We used the relation between PE wind mass loss, M˙PE$\dot{M}_{\rm PE}$, and (attenuated) stellar 0.1–1 keV X-ray luminosity, LX derived from 2D MHD simulations Eq. (9) in Ercolano et al. (2021). Their relation is nonlinear and convex with LX due to self-shielding by the PE wind itself (see also Gorti & Hollenbach 2009). At Lx < 1028 ergs s−1, we used a matched linear relation, M˙PE=2.56×1010(LX/1028)$\dot{M}_{\rm PE} = 2.56 \times 10^{-10} (L_X/10^{28})$ M yr−1.

We adopted Dartmouth stellar evolution models for pre-MS values of R*, T*, and L* (Dotter et al. 2008; Feiden et al. 2015, see our Appendix B.2 for an investigation into the effects of model choice). We used the self-similar solution for time- dependent accretion in a viscous disk with a uniform viscosity (Lynden-Bell & Pringle 1974; Hartmann et al. 1998): M˙acc=M˙0(1+t/t0)5/4.\dot{M}_{\rm acc} = \dot{M}_0 \left(1 + t/t_0\right)^{-5/4}.(8)

The initial accretion rate was taken to be M˙0=2.5×107M1.6$\dot{M}_0 = 2.5\times 10^{-7} M_{\star}^{1.6}$ M yr−1, based on observations of “flat-spectrum” young stellar objects considered to be in the Class I–II transition (Eq. (3) in Gehrig et al. 2023) and a viscous timescale of t0 = 1 × 105 yr.

For given star and prescribed M˙acc$\dot{M}_{\rm acc}$, the remaining parameters are B* and P*. We adopted a nominal value for B* of 0.44 kG as 22% of the total field strength (Lavail et al. 2019) of 2 kG, typical for magnetically saturated stars (Reiners et al. 2022), but considered a range of 0.1–1 kG (see Fig. 9 in Gehrig et al. 2022). Although protostars should spin up as they contract, observations suggest that the P* of Class I and II objects remain relatively constant due to torques from the disk (Königl 1991; Ostriker & Shu 1995, see also the discussion in Mueller et al. 2024). We assume that P* is fixed to a value between 1 and 10 days while a disk is present (Smith et al. 2023; Serna et al. 2024), but consider different scalings of P* with other parameters in Appendix B.

3 Results

3.1 Fiducial and solar-mass scenarios

At the transition to the Class II phase, low-mass protostars are still inflated but already rapidly rotating (P* ∼ 1−10 days, Smith et al. 2023). As a result, M˙acc$\dot{M}_{\rm acc}$ is initially below M˙crit$\dot{M}_{\rm crit}$, R > 1, and for typical P* and B* (see below), accretion will be accompanied by an IDW. Based on Eq. (7), for solar-mass stars the optical depth to 0.1–1 keV X-rays will be 300 and the disk is completely shielded at the beginning of the Class II phase2. If P* is nearly constant due to the star’s exchange of angular with the disk, then the pre-MS star and its magnetic field will contract with time, M˙crit$\dot{M}_{\rm crit}$ will decline faster than M˙acc$\dot{M}_{\rm acc}$, and R will decrease. Eventually, reaches unity, disk accretion enters the magnetospheric regime, the IDW disappears, and stellar X-rays can reach the disk unimpeded. If M˙PE$\dot{M}_{\rm PE}$ also exceeds M˙acc$\dot{M}_{\rm acc}$ then PE opens a gap in the disk (at time tgap), causing the disk to more rapidly dissipate. The dependence of Mcrit on B* and inverse dependence on P* (Eq. (4)) means that disks around more rapidly rotating stars with stronger fields will tend to be longer-lived.

Figure 1 plots M˙acc$\dot{M}_{\rm acc}$, the M˙PE$\dot{M}_{\rm PE}$, and the critical accretion rate at which R = 1 for four scenarios. Time is Dartmouth stellar model time in millions of years3. In Scenario A (M* = 1 M, B* = 0.44 kG and P* = 3 days), M˙crit$\dot{M}_{\rm crit}$ falls below M˙acc$\dot{M}_{\rm acc}$, the IDW ceases, and a gap opens at tgap = 2.3 Myr. In scenario B, with a stronger field (B* = 0.6 kG) and faster rotation (P* = 2 days), M˙acc$\dot{M}_{\rm acc}$ remains below M˙crit$\dot{M}_{\rm crit}$ for the entire simulation, and the IDW persists. A gap could eventually open once M˙acc$\dot{M}_{\rm acc}$ is insufficient for the wind to absorb most soft X-rays (see Sect. 4.1). In Scenario C, the slower rotation (P* = 7 days) causes the IDW to disappear by 0.5 Myr, but X-ray emission is also lower due to slow rotation and more time elapses (5.7 Myr) before M˙acc$\dot{M}_{\rm acc}$ falls below M˙PE$\dot{M}_{\rm PE}$. Scenario D is the same as A except for a lower-mass (0.3 M) star; PE is suppressed by an IDW for longer (13.3 Myr), mostly due to lower M˙acc$\dot{M}_{\rm acc}$. For a given M* and M˙acc$\dot{M}_{\rm acc}$ history, tgap depends only on B* and P*. Figure 2 plots this dependence for M* = 1M over the ranges B* ∈ [0.1, 1] kG and P* ∈ [1, 10] days, showing the locations of Scenarios A, B, and C in this parameter space.

3.2 Disk lifetime and stellar mass

For fixed M* = 1 M but varying B* and P*, the gap opening time, tgap, spans more than an order of magnitude (Fig. 2), consistent with the variation inferred from disk statistics (Pfalzner & Dincer 2024). Superposed on this variation is a trend of longer disk life with lower stellar mass. This pattern is largely due to the correlation between M˙acc$\dot{M}_{\rm acc}$ and M*, which in turn is based on observations of flat-spectrum protostars (Gehrig et al. 2023) but which is also widely observed among TTSs (Manara et al. 2023). A longer elapsed time is required for Mcrit to decline to reach a lower M˙acc$\dot{M}_{\rm acc}$ and the IDW to disappear (compare scenarios A and D in Fig. 1). Figure 3 shows that tgap systematically increases with lower M* for a fixed B* = 0.44 kG, except for the slowest rotators. This is consistent with observations showing longer disk lifetimes around lower-mass stars (Bayo et al. 2012; Ribas et al. 2015; Silverberg et al. 2020). tgap increases with decreasing P*: the tendency of VLMS TTSs to be more rapidly rotating (e.g., Herbst et al. 2002; Kounkel et al. 2022; Smith et al. 2023) means that disks around these objects will be even longer lived. The strong R* and weaker M* dependence of Eq. (7) also indicate that, at a given M˙acc$\dot{M}_{\rm acc}$, a more compact IDW around a lower-mass star will more strongly attenuate X-rays.

For a more realistic treatment of the star-disk angular momentum exchange that affects P*, and thus disk dissipation, we performed disk evolution simulations using the selfconsistent hydrodynamic model described in Gehrig et al. (2022; Cecil et al. 2024). This includes star-disk torques as well as the effects of a PE wind. The results are described in Appendix A and support the dependence of tgap on M*, P*, and B* illustrated by Figs. 2 and 3.

thumbnail Fig. 1

Four scenarios of disk evolution. (A) With M* = 1 M, B* = 0.44 kG, and P* = 3 days, M˙acc$\dot{M}_{\rm acc}$ exceeds M˙crit$\dot{M}_{\rm crit}$ at 2.3 Myr, Rmag retreats inside Rcor, an IDW ceases, and stellar X-rays open a disk gap. (B) In the case of a more rapidly rotating (P* = 2 days) star with a stronger (0.6 kG) field, M˙acc$\dot{M}_{\rm acc}$ never exceeds M˙crit$\dot{M}_{\rm crit}$ and no gap opens. (C) For a more slowly rotating (P* = 7 day) star, there is no IDW but lower LX and PE delays tgap to 5.7 Myr. (D) In the 0.3 M case, lower M˙acc$\dot{M}_{\rm acc}$ remains below M˙crit$\dot{M}_{\rm crit}$ for longer, delaying tgap to 13.3 Myr.

thumbnail Fig. 2

tgap (Myr) as a function of B* and P* for M* = 1 M. Labeled stars correspond to scenarios A, B, and C in Fig. 1 In the shaded region, M˙acc$\dot{M}_{\rm acc}$ remains below M˙crit$\dot{M}_{\rm crit}$ and no gap opens for the entire simulation.

thumbnail Fig. 3

Same as Fig. 2, but with tgap (Myr) as a function of M* and P* for B* = 0.44 kG. Stars correspond to scenarios A, C, and D in Fig. 1.

thumbnail Fig. 4

Threshold accretion rate below which a wind cannot effectively shield the disk from PE-driving X-ray radiation, vs. stellar mass. Curves are calculations for different rotation periods and ages of 1–10 Myr, and points are values of M˙acc$\dot{M}_{\rm acc}$ from Manara et al. (2023).

4 Discussion

4.1 Threshold accretion rate

Decline in M˙acc$\dot{M}_{\rm acc}$ with time means that the gas column in any IDW (see Eq. (7)) will eventually be insufficient (NH ≪ 1021 cm−2) to shield the disk from stellar X-rays. Thus, PE will inevitably trigger accelerated dissipation of the disk and rapid decline in M˙acc$\dot{M}_{\rm acc}$. This effect should appear as a baseline M˙acc$\dot{M}_{\rm acc}$ level, below which few accreting TTSs are found even in sensitive surveys. Thanathibodee et al. (2023) report that the lowest T Tauri accretors have ~ 10−10 M yr1, well above the detection limit of 1 - 5 10−11 M yr1. We calculated the baseline by assuming Rmag = Rcor (most compact IDW) in Eq. (6), adopting f = 0.2 as before, and setting NH to 1021 cm−2, yielding: M˙thresh=1.5×1010Myr1P1/2R1/4M1/4(T/4000)1/4.\dot{M}_{\rm thresh} = 1.5 \times 10^{-10} M_{\odot} {\rm yr}^{-1} P^{1/2} R_{\star}^{1/4} M_{\star}^{1/4} \left(T_{\star}/4000\right)^{1/4}(9)

Figure 4 compares this M*-dependent threshold for two values of P* and a range of ages to M˙acc$\dot{M}_{\rm acc}$ values compiled by Manara et al. (2023). Rather than minimum M˙acc$\dot{M}_{\rm acc}$ values reflecting the rate of PE driven by EUV radiation from the star (Thanathibodee et al. 2023), we interpret this threshold as the minimum M˙acc$\dot{M}_{\rm acc}$ at which an IDW can quench X-ray-driven PE (see also Ercolano et al. 2023).

4.2 Accretion variability from a limit cycle

Feedback between accretion in the inner disk, the attenuation of stellar X-rays, and PE, combined with a delay between the removal of mass by PE and inner disk accretion (i.e., the viscous transport time) raises the possibility of an oscillatory limit cycle (Strogatz 2000) that produces variability on that delay timescale. We demonstrated this phenomenon using a minimalist “toy” model that captures these feedbacks (rather than the complex model described in Sect. 2) and that is represented by time-dependent equations for the inner disk accretion and PE rates: M˙acc(t)=M˙0M˙PE(tΔt),\dot{M}_{\rm acc}(t) = \dot{M}_0 - \dot{M}_{\rm PE}\left(t - \Delta t\right), \label{eqn:toy1}(10a) M˙PE(t)=M˙PE0exp[afM˙acc(t)δ(M˙acc(t)<M˙crit)].\dot{M}_{\rm PE}(t) = \dot{M}_{\rm PE}^0 \exp \left[-a f \dot{M}_{\rm acc}(t) \delta\left(\dot{M}_{\rm acc}(t) &lt; \dot{M}_{\rm crit}\right)\right].(10b)

Equation (10a) describes the decrease in inner disk accretion at time, t, from a baseline rate, M˙0$\dot{M}_0$, due to PE at time t - △t, with △t representing the transport time (e.g., viscous timescale) from the location of the PE wind to the inner disk edge. Equation (10b) describes the exponential attenuation of X-ray-driven PE below an unattenuated rate, M˙PE0,$\dot{M}_{PE}^0$ if an IDW is present (a condition represented by the Boolean integer, δ). The optical depth of the IDW to X-rays is assumed to scale with IDW mass flow times a fixed attenuation coefficient, a, the IDW mass flow is M˙acc$\dot{M}_{\rm acc}$ times a fixed ejection/accretion ratio, f , and PE is assumed to scale linearly with X-ray irradiation.

The equations are iteratively solved, with each iteration representing an elapsed time, △t. If the feedback is sufficiently strong – that is, if M˙PE0>M˙0M˙crit$\dot{M}_{PE}^0 > \dot{M}_0 - \dot{M}_{\rm crit}$ and the product a × f exceeds a threshold value – the solutions of Eqs. (10a) and (10b) become multivalued or “multifurcate” (Fig. 5a). Physically, this could manifest itself as unsteady accretion over the accretion timescale, △t, of the inner disk. Based on the observed dispersion of T Tauri disk masses and accretion rates (e.g, Almendros-Abad et al. 2024), △t could be ~104 - 105 yr. Figure 5b shows the domain of multifurcated solutions if M˙crit$\dot{M}_{\rm crit}$ is 60% of M˙0$\dot{M}_0$. For M˙acc$\dot{M}_{\rm acc}$ [0.4M˙0,M˙0]$\in$ [$0.4\,\dot{M}_0$,$\dot{M}_0$] and high a × f , there is a regime of unsteady accretion (shaded gray region in Fig. 5b with contours representing the ratio of the maximum to minimum values in dex). For low a × f , M˙acc$\dot{M}_{\rm acc}$ is steady but attenuated by an IDW, if M˙PE0<0.4M˙0$\dot{M}_{PE}^0 &lt; 0.4\,\dot{M}_0$, there is no IDW and M˙acc$\dot{M}_{\rm acc}$ = M˙0$\dot{M}_0$, PE completely halts accretion. The unsteady mode dominates as M˙0$\dot{M}_0$ approaches M˙crit$\dot{M}_{\rm crit}$. Variation in M˙acc$\dot{M}_{\rm acc}$ on timescales of ∼104−105 yr, in addition to shorter-term variability (e.g., Claes et al. 2022) could explain some of the observed two- orders-of-magnitude scatter around the relation between M˙acc$\dot{M}_{\rm acc}$ and M* (Manara et al. 2023).

thumbnail Fig. 5

Panel a: above a critical value of the product of the X-ray attenuation coefficient, a, times the ejection/accretion ratio, f (Eq. (10b)), solutions for the inner disk accretion rate multifurcate. Panel b: domains of solutions to Eqs. (10a) and (10b) in af˙$a \dot f$ vs. the maximum PE rate. Contours are the ratio of the maximum to minimum values of M˙acc$\dot{M}_{\rm acc}$ expressed in dex. The dotted line is the parameter range for the case shown in the top panel.

4.3 Disk lifetime and chemistry

Infrared spectroscopy has revealed gas chemistry indicative of a super-solar C/O (~ 1) in the inner disks of some VLMSs (Najita et al. 2013; Tabone et al. 2023; Arabhavi et al. 2024; Kanwar et al. 2024), but not all (Xie et al. 2023). This difference between VLMS and TTS disks has been explained as the result of a more compact snowline, a shorter viscous timescale, and the faster inward transport of gas depleted in oxygen by the condensation of H2O- and CO2-rich solids (Mah et al. 2023). However, it could also be explained by more pervasive suppression of PE that would otherwise form a gap and inhibit later inward transport of O-depleted gas from the outer disk (Gasman et al. 2023). Of four VLMS disks studied with JWST and compared by Kanwar et al. (2024), all have low C/O chemistry except the one (Sz 124) with the highest M˙acc$\dot{M}_{\rm acc}$ and apparent LX, suggesting a lack of an X-ray attenuating IDW. Sz 124 is also the only one with detectable [Ne II] emission indicative of X-ray-driven PE, perhaps sufficient to open a gap such as the one suggested by the highest-resolution analysis (Jennings et al. 2022) and inhibit later transport of O-poor gas to the inner disk.

4.4 Disk structure and planet formation

Diverging disk evolutionary pathways could end in markedly different planetary outcomes (Fig. 6). If PE opens a gap early in the disk history,a typical scenario among solar-mass stars (Fig. 3), a pressure maximum immediately exterior to the gap can trap migrating solids; that is, millimeter-to-centimeter-sized “pebbles” (e.g., Pinilla et al. 2021). This starves the inner disk of condensible solids, especially volatiles (Kalyaan et al. 2021), while accumulation in the pressure maximum could cause the formation of massive cores and giant planets (e.g., Lau et al. 2022). A statistical link between transition disks – a possible stage after gap opening – and giant planets has been previously described (van der Marel & Mulders 2021). Furthermore, the inward migration of giant planets formed further out could be slowed or halted by the evacuation of the inner disk (Alexander & Pascucci 2012; Monsch et al. 2021)4. In longer-lived disks where gap opening is delayed, a likely outcome around lower-mass stars (Fig. 3), solids will drift unimpeded close to the star (Pinilla et al. 2020), explaining why the continuum emission from dust VLMS disks tends to be more compact (Ansdell et al. 2018; Kurtovic et al. 2021). Planet formation and migration in denser compact disks could yield multiple short-period planets in near-resonant chains, which are observed around M dwarfs more frequently than solar-type stars (Gaidos et al. 2016; Hardegree-Ullman et al. 2019; Ment & Charbonneau 2023).

4.5 Caveats

Sensitivity to key parameters is investigated in Appendix B. A disk’s outcome depends on the stellar period (Appendix B.1), the radius (Appendix B.2), and the time-dependence of M˙acc$\dot{M}_{\rm acc}$ (Appendix B.3), as well as B*. We emphasize that while this sensitivity makes exact predictions more difficult, it is the underlying mechanism by which diversity in disk evolution can arise. In the analytic model, we prescribe the values or time-dependence of these parameters, but it is derived self-consistently in our numerical modeling (Appendix A). The degree of shielding of the disk from X-rays for a given M˙acc$\dot{M}_{\rm acc}$ is sensitive to the poorly known geometry and kinematics of the IDW and ejection- accretion ratio (see Appendix B.4). The sonic speed at the base of the wind could vary by a factor of three for temperatures between 103 and 104K. Shielding could also occur from a more static wall or envelope supported by magnetic pressure (Zhu et al. 2024). The PE rates are based entirely on 2D wind models and have only been tested by observations in a limited way (e.g., Ercolano et al. 2021). A model that includes detailed chemistry suggests that PE rates could be lower due to more efficient cooling (Sellek et al. 2024).

5 Conclusions

Photoevaporation driven by soft X-rays from the central star is a pathway by which PPDs develop gaps and inner cavities and eventually dissipate. An IDW or other vertical structure near the inner edge of a disk can shield the outer disk from X-rays, suppressing PE and extending its lifetime. Such a wind is present when the partially ionized disk is truncated by the stellar magnetic field beyond the corotation radius, but will disappear when the contraction of the protostar, its field, and the truncation radius outpaces the effect of a declining accretion rate, allowing X-rays to reach the disk. The timing of these events depends sensitively on the stellar mass, large-scale dipole field strength, and rotation. An analytic model of star and disk evolution shows that this inner-outer disk connection explains:

  • the observed trend of increased disk lifetime and decreased occurrence of disk structure around lower-mass stars;

  • the 1 dex range in disk lifetime for a given mass, with rapidly rotating stars with stronger magnetic fields having longer disk lifetimes;

  • the paucity of disks with accretion rates lower than ~10−10 M yr−1, as below that level the IDW is too diffuse to significantly attenuate X-rays and prevent disk dispersal;

  • some of the observed scatter in disk accretion rates may be a manifestation of kiloyear-timescale variability, specifically a limit cycle due to feedback between PE, inner disk accretion, and the attenuation of X-rays by an accretion-fed wind.

A dichotomous outcome in disk evolution suggests a corresponding divergence in disk chemistry and planet formation outcomes. Disks with early-developing gaps and cavities tend to generate giant planets on 1–10 au orbits. Longer-lived disks where suppression of PE has delayed gap formation are more likely to have higher C/O due to the unimpeded inward transport of O-poor gas, and spawn systems of multiple close-in super-Earths and Neptunes.

Direct evidence that the attenuation of X-rays and suppression of PE prolongs disk life is needed. Laos et al. (2022) found that 0.8–3 keV X-ray emission from six M dwarf stars hosting long-lived disks is not systematically lower than that from their diskless coeval counterparts. However, this energy range is less sensitive to attenuation by gas and less relevant to PE. Moreover, if the wind has a wide opening angle (e.g., ~60 deg, Nemer & Goodman 2024) then attenuation would not be observed along our line of sight among half of the disks in a randomly oriented sample. Observations of more long-lived disks at softer energies (<1 keV) are warranted. Key to more robust tests of this paradigm are near-simultaneous measurements of the large- scale magnetic field strength, accretion rate, and winds as might be accomplished with synergistic combinations of ground- and space-based spectroscopy in the UV and infrared.

thumbnail Fig. 6

Illustration of two scenarios of disk evolution and planet formation. Left: if Rmag < Rcor, unimpeded X-ray irradiation by the central star opens a gap in the disk that eventually depletes the inner disk and stymies planet formation there. A pressure bump outside the cavity can trap solids and enable the formation of one or more giant planets. Right: if Rmag > Rcor, an IDW shields the disk further out from X-rays, slowing or halting PE and delaying or preventing the opening of a gap. Solids drift accumulate closer to the star, leading to more compact dust disks and close-in multi-planet systems.

Acknowledgements

The authors thank Kees Dullemond for a thoughtful and careful review, and Ilaria Pascucci for helpful comments. E.G. was supported as a Gauss Professor at the University of Göttingen by the Niedersächsische Akaemie der Wissenschaften.

Appendix A Numerical modeling

We compared the analytic results represented in Fig. 2 to the output of a numerical model that combines stellar and disk evolution self-consistently. The implicit hydrodynamic TAPIR code (e.g., Ragossnig et al. 2020) includes a detailed calculation of the disk truncation radius (Steiner et al. 2021) and a stellar spin model (Gehrig et al. 2022), and calculates photoevaporative mass loss from the disk (Cecil et al. 2024).

There are several key differences between the numerical and analytic models. First, the stellar rotation period is not fixed but evolves according to changes in moment of inertia and torques from the disk (Gehrig et al. 2022). Second, the code calculates the truncation radius based on the stellar magnetic pressure, the ram pressure of the infalling material, and the static gas pressure in the innermost disk region (e.g., Romanova et al. 2002; Bessolaz et al. 2008; Steiner et al. 2021). Finally, disk accretion is calculated continuously and self-consistently with an α viscosity, a “dead” zone, and PE as a sink term. Feedbacks between PE, inner disk accretion, the location of the magnetic truncation radius, the operation of an IDW, and the attenuation of stellar X- rays are included. We use the model as described in Cecil et al. (2024), but included a “switch” that turns off PE when reaches unity. To avoid numerical difficulties, PE is reduced smoothly starting at R = 0.9 to zero at R = 0.9.

In all runs, the α viscosity parameter is taken to be 5 × 10 (2 × 10−4) outside of (within) the dead zone. In these simulations, disk lifetime refers to the time at which M˙acc$\dot{M}_{\rm acc}$ falls below 10−11 M yr−1. We assume an initial stellar age of 1 Myr, and stellar radii and luminosities are taken from Baraffe et al. (2015). The initial accretion rates are 1.3 10−9 M /yr and 7.7 10−9 M /yr and the X-ray luminosities (in ergs sec−1) are taken to be log LX = 29.5 and 30.3 for 0.3 and 1.0 M stars, respectively.

Figures A.1 and A.2 plot disk lifetime predicted by the model vs. B* and initial (not fixed) P* for 1.0 and 0.3M stars, respectively. These TAPIR results support the overall prediction of the analytic model that PE is suppressed and disk lifetimes are longer around more rapidly rotating stars with strong magnetic fields (compare to Fig. 2), and that this effect is more pervasive around stars with lower masses (compare to Fig. 3). Differences between the numerical and analytic predictions are attributable to the addition of stellar spin evolution in the former. For a given magnetic field strength, stellar rotation tends to evolve toward an equilibrium state that is different from the initial condition (e.g., Ireland et al. 2021; Gehrig et al. 2022; Serna et al. 2024). One notable difference is shorter disk lifetime for stars with the strongest B* and shortest initial P* (upper right-hand corner of Figs. A.1A.2). This is a result of a stronger field exerting torques on both the star and disk, enhancing accretion towards the star and moving the disk truncation radius inward, while spinning down the star and moving the corotation radius outward. This tends to suppress the formation of an IDW which would otherwise protect the disk from PE. However, the condition of both rapid rotation and a strong field would probably be transient in the Class I phase and not a realistic state at the beginning of the Class II phase.

thumbnail Fig. A.1

Disk lifetime vs. B* and initial (not fixed) P* for a 1-M star predicted by the numerical TAPIR star-disk interaction model. Compare to Fig. 2.

thumbnail Fig. A.2

Same as Fig. A.1 but for a 0.3-M star.

Appendix B Sensitivity analyses

B.1 Rotational evolution

Models of star-disk interaction predict that rotation periods do not remain constant as assumed in our analytical model (Sect. 2), but instead that they evolve slightly over 1–10 Myr. The “stellar magnetospheric” (large-scale dipole; Königl 1991) and “trapped magnetic flux” models of star-disk interaction (Ostriker & Shu 1995; Johns-Krull & Gafford 2002) predict rotation periods that scale as PR18/7M˙3/7$P_{\star} \propto R_{\star}^{18/7}\dot{M}^{-3/7}$ and PR4M˙1$P_{\star} \propto R_{\star}^{4}\dot{M}^{-1}$, respectively. (See Mueller et al. 2024 for an observational test of these models.). We re-calculated gap opening times for the solar-mass case using these relations. Figure B.1 plots these results, where the ordinate is now the minimum value of P* that occurs during the simulation. The overall distribution of gap opening times is similar (compare to Fig. 2) except that in the trapped flux scenario, late (rather than no) gap opening is predicted for rapidly rotating stars with high B*, and no gap opens for very slow rotators. We also find that the scaling relations do not predict an interval of quasiconstant rotation for the 0.3 M case, which could reflect either reality or the limitations of our prescription for the accretion rate history of VLMSs.

thumbnail Fig. B.1

tgap vs. B* and P* for 1 M. Same as Fig. 2 except using the stellar magnetosphere (Königl 1991) (top) and trapped flux (Ostriker & Shu 1995) models to calculate evolution of P*. The ordinate is now the minimum rotation period over the disk lifetime.

thumbnail Fig. B.2

Stellar radius evolution for 1- and 0.3-Mstars predicted by the Dartmouth standard model (Dotter 2016), the Dartmouth magnetic model , which includes the inflationary effect of magnetic fields via spots and inhibition of convection (Feiden 2016), and the Baraffe et al. (2015) model. In the Dartmouth models a different set of abundances and opacities results in prolonged deuterium burning compared to the Baraffe et al. (2015) model (G. Feiden, pers. comm.). This delays stellar contraction and potential affects disk evolution via the radius dependence of Eqn. 4.

B.2 Stellar model

For a given stellar mass, rotation period, and dipole magnetic field strength, the critical value of accretion rate at which an IDW appears is highly sensitive to stellar radius (Eqn. 4). Predictions of stellar radius by models depend on the physics that is included (Fig. B.2). Deuterium burning increases the entropy and can delay the early contraction of pre-MS stars (Stahler 1988). The timescale of D-burning in fully convective pre-MS stars is sensitive to the central temperature, which in turn depends on the adopted abundances and opacities (compare Darmouth and Baraffe et al. (2015) models in Fig. B.2). In most cases, deuterium burning occurs too early to matter, with the possible exception of slowly rotating M dwarfs, where it is not completely burned for 3 Myr (Fig. B.2). For these cases it could prolong the disk by delaying contraction and increasing M˙crit$\dot{M}_{\rm crit}$.

Internal magnetic fields can significantly inflate stellar radii over a range of masses and ages via suppression of convection or spots (Feiden 2016). This would increase M˙crit$\dot{M}_{\rm crit}$ and, all else being equal, prolong the lifetime of disks, except for slow rotators where stellar X-ray emission is the limiting factor (e.g., Fig. 1c).

B.3 Accretion rate

Whether PE opens a cavity or is attenuated by an IDW depends on the disk accretion history. We demonstrated this sensitively by two departures from the nominal case for a 1 M star (Eqn. 8), first by tripling the initial value M˙acc$\dot{M}_{\rm acc}$(0), and second by changing the power law index from −5/4 to −2. The index value of -2 corresponds to the case where the disk viscosity depends on radius as r3/2, the maximum dependence consistent with observed surface density profiles and disk demographics (an index greater than −5/4 and negative radius dependence is also excluded by these observations Alexander et al. 2023). These changes slightly decrease and significantly increase the B*P* parameter range where PE is completely suppressed, respectively (top and bottom panels of Fig. B.3)

B.4 Ejection/accretion ratio

For a given accretion rate, the attenuation of X-rays by an IDW depends on the ejection/accretion ratio f , which is poorly constrained but thought to range between ~0.01 and 1 (Watson et al. 2016; Serna et al. 2024). Figure B.4 plots the X-ray attenuation at the end of the IDW phase (~ = 1) vs. M* using Eqn. 6, our nominal two-temperature model of coronal X-ray emission, and different values of f . Attenuation is also sensitive to the flow velocity at the foot of the wind (here assumed to be the sound speed) and the angle of the wind with respect to the disk normal (see Sect. 4.5).

B.5 Stellar coronal X-ray spectrum

X-ray absorption is strongly energy-dependent, soft photons being absorbed more efficiently than hard photons, and thus the effective attenuation over 0.1–1 keV by intervening gas depends on the spectrum of the intrinsic emission from the star. We compared attenuation for different coronal spectra typically used for TTSs, i.e., isothermal plasmas with kT = 0.5, 0.7, 1.0, and 2.0 keV (Telleschi et al. 2007, see our Sect. 2 for details). Figure B.4 plots the calculated attenuation for single-temperature (0.7 and 2 keV) model plasmas, plus a mix of two plasma components at 0.5 keV plus 1 keV with an emission measure ratio of 2:1. The effective attenuation depends only weakly on the temperature because the different spectra have a similar shape over 0.1–1 keV. This is a consequence of the spectra being dominated by a flat bremsstrahlung continuum plus widely distributed emission lines. The denser line forest between 0.7–1 keV is present and similar for all spectra in this temperature range. The strongly temperature-dependent turnover of the bremsstrahlung spectrum from flat to exponentially decreasing occurs around kT and therefore affects spectra at energies above 1 keV, a range not relevant to PE.

thumbnail Fig. B.3

tgap vs. B* and P* for 1 M. Same as Fig. 2 except with a trebled initial accretion rate (top) or an accretion rate power law index of −2 instead of −5/4.

thumbnail Fig. B.4

0.1–1 keV X-ray attenuation at ℛ = 1 vs. stellar mass at the end of the IDW phase for different values of ejection/accretion ratio f.

thumbnail Fig. B.5

Attenuation of 0.1–1 keV X-rays for different coronal emission models.

References

  1. Akeson, R. L., Jensen, E. L. N., Carpenter, J., et al. 2019, ApJ, 872, 158 [Google Scholar]
  2. Alexander, R. D., & Pascucci, I. 2012, MNRAS, 422, L82 [NASA ADS] [CrossRef] [Google Scholar]
  3. Alexander, R., Rosotti, G., Armitage, P. J., et al. 2023, MNRAS, 524, 3948 [NASA ADS] [CrossRef] [Google Scholar]
  4. Almendros-Abad, V., Manara, C. F., Testi, L., et al. 2024, A&A, 685, A118 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  5. Ansdell, M., Gaidos, E., Rappaport, S. A., et al. 2016, ApJ, 816, 69 [Google Scholar]
  6. Ansdell, M., Williams, J. P., Trapman, L., et al. 2018, ApJ, 859, 21 [NASA ADS] [CrossRef] [Google Scholar]
  7. Arabhavi, A. M., Kamp, I., Henning, T., et al. 2024, Science, 384, 1086 [NASA ADS] [CrossRef] [Google Scholar]
  8. Armitage, P. J. 2011, ARA&A, 49, 195 [Google Scholar]
  9. Arnaud, K. A. 1996, in ASP Conf. Ser., 101, 17 [Google Scholar]
  10. Baraffe, I., Homeier, D., Allard, F., & Chabrier, G. 2015, A&A, 577, A42 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  11. Bayo, A., Barrado, D., Huélamo, N., et al. 2012, A&A, 547, A80 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  12. Bessolaz, N., Zanni, C., Ferreira, J., Keppens, R., & Bouvier, J. 2008, A&A, 478, 155 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  13. Birnstiel, T. 2024, ARA&A, 62, 157 [NASA ADS] [CrossRef] [Google Scholar]
  14. Bouvier, J., Grankin, K. N., Alencar, S. H. P., et al. 2003, A&A, 409, 169 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  15. Bouvier, J., Alencar, S. H. P., Harries, T. J., Johns-Krull, C. M., & Romanova, M. M. 2007, in Protostars and Planets V, eds. Reipurth, B., Jewitt, D., & Keil, K. (Tucson: University of Arizona Press), 479 [Google Scholar]
  16. Bustamante, I., Merín, B., Bouy, H., et al. 2016, A&A, 587, A81 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  17. Cecil, M., Gehrig, L., & Steiner, D. 2024, A&A, 687, A136 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  18. Claes, R. A. B., Manara, C. F., Garcia-Lopez, R., et al. 2022, A&A, 664, L7 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  19. Doppmann, G. W., Greene, T. P., Covey, K. R., & Lada, C. J. 2005, AJ, 130, 1145 [Google Scholar]
  20. Dotter, A. 2016, ApJS, 222, 8 [Google Scholar]
  21. Dotter, A., Chaboyer, B., Jevremovic, D., et al. 2008, ApJS, 178, 89 [NASA ADS] [CrossRef] [Google Scholar]
  22. Drake, J. J., Ercolano, B., Flaccomio, E., & Micela, G. 2009, ApJ, 699, L35 [NASA ADS] [CrossRef] [Google Scholar]
  23. Dullemond, C. P., Dominik, C., & Natta, A. 2001, ApJ, 560, 957 [NASA ADS] [CrossRef] [Google Scholar]
  24. Ercolano, B., & Pascucci, I. 2017, Royal Soc. Open Sci., 4, 170114 [NASA ADS] [CrossRef] [Google Scholar]
  25. Ercolano, B., Picogna, G., Monsch, K., Drake, J. J., & Preibisch, T. 2021, MNRAS, 508, 1675 [NASA ADS] [CrossRef] [Google Scholar]
  26. Ercolano, B., Picogna, G., & Monsch, K. 2023, MNRAS, 526, L105 [NASA ADS] [CrossRef] [Google Scholar]
  27. Espaillat, C., Furlan, E., D’Alessio, P., et al. 2011, ApJ, 728, 49 [NASA ADS] [CrossRef] [Google Scholar]
  28. Fang, M., Wang, L., Herczeg, G. J., et al. 2023, Nat. Astron., 7, 905 [Google Scholar]
  29. Feiden, G. A. 2016, A&A, 593, A99 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  30. Feiden, G. A., Jones, J., & Chaboyer, B. 2015, in Cambridge Workshop on Cool Stars, Stellar Systems, and the Sun (Berlin: Springer), 18, 171 [Google Scholar]
  31. Fernandes, R. B., Long, Z. C., Pikhartova, M., et al. 2018, ApJ, 856, 103 [NASA ADS] [CrossRef] [Google Scholar]
  32. Flaherty, K. M., Muzerolle, J., Rieke, G., et al. 2012, ApJ, 748, 71 [NASA ADS] [Google Scholar]
  33. Flaherty, K., Hughes, A. M., Mamajek, E. E., & Murphy, S. J. 2019, ApJ, 872, 92 [Google Scholar]
  34. Gaidos, E., Mann, A. W., Kraus, A. L., & Ireland, M. 2016, MNRAS, 457, 2877 [Google Scholar]
  35. Gaidos, E., Mann, A. W., Rojas-Ayala, B., et al. 2022, MNRAS, 514, 1386 [CrossRef] [Google Scholar]
  36. Gaidos, E., Thanathibodee, T., Hoffman, A., et al. 2024, ApJ, 966, 167 [NASA ADS] [CrossRef] [Google Scholar]
  37. Gasman, D., van Dishoeck, E. F., Grant, S. L., et al. 2023, A&A, 679, A117 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  38. Gehrig, L., Steiner, D., Vorobyov, E. I., & Güdel, M. 2022, A&A, 667, A46 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  39. Gehrig, L., Gaidos, E., & Güdel, M. 2023, A&A, 675, A179 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  40. Ghosh, P., & Lamb, F. K. 1979, ApJ, 232, 259 [Google Scholar]
  41. Gorti, U., & Hollenbach, D. 2009, ApJ, 690, 1539 [Google Scholar]
  42. Güdel, M., & Telleschi, A. 2007, A&A, 474, L25 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  43. Güdel, M., Skinner, S. L., Briggs, K. R., et al. 2005, ApJ, 626, L53 [Google Scholar]
  44. Güdel, M., Skinner, S. L., Mel’Nikov, S. Y., et al. 2007a, A&A, 468, 529 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  45. Güdel, M., Telleschi, A., Audard, M., et al. 2007b, A&A, 468, 515 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  46. Güdel, M., Skinner, S. L., Audard, M., Briggs, K. R., & Cabrit, S. 2008, A&A, 478, 797 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  47. Güdel, M., Lahuis, F., Briggs, K. R., et al. 2010, A&A, 519, A113 [CrossRef] [EDP Sciences] [Google Scholar]
  48. Günther, H. M., Birnstiel, T., Huenemoerder, D. P., et al. 2018, AJ, 156, 56 [Google Scholar]
  49. Güver, T., & Özel, F. 2009, MNRAS, 400, 2050 [Google Scholar]
  50. Hardegree-Ullman, K. K., Cushing, M. C., Muirhead, P. S., & Christiansen, J. L. 2019, AJ, 158, 75 [Google Scholar]
  51. Hartmann, L., Calvet, N., Gullbring, E., & D’Alessio, P. 1998, ApJ, 495, 385 [Google Scholar]
  52. Herbst, W., Bailer-Jones, C. A. L., Mundt, R., Meisenheimer, K., & Wackermann, R. 2002, A&A, 396, 513 [EDP Sciences] [Google Scholar]
  53. Ireland, L. G., Zanni, C., Matt, S. P., & Pantolmos, G. 2021, ApJ, 906, 4 [Google Scholar]
  54. Jeffries, R. D., Jackson, R. J., Briggs, K. R., Evans, P. A., & Pye, J. P. 2011, MNRAS, 411, 2099 [Google Scholar]
  55. Jennings, J., Booth, R. A., Tazzari, M., Clarke, C. J., & Rosotti, G. P. 2022, MNRAS, 509, 2780 [NASA ADS] [Google Scholar]
  56. Johns-Krull, C. M., & Gafford, A. D. 2002, ApJ, 573, 685 [Google Scholar]
  57. Kalyaan, A., Pinilla, P., Krijt, S., Mulders, G. D., & Banzatti, A. 2021, ApJ, 921, 84 [CrossRef] [Google Scholar]
  58. Kanwar, J., Kamp, I., Jang, H., et al. 2024, A&A, 689, A231 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  59. Königl, A. 1991, ApJ, 370, L39 [CrossRef] [Google Scholar]
  60. Kounkel, M., Stassun, K. G., Bouma, L. G., et al. 2022, AJ, 164, 137 [NASA ADS] [CrossRef] [Google Scholar]
  61. Kraus, A. L., Ireland, M. J., Hillenbrand, L. A., & Martinache, F. 2012, ApJ, 745, 19 [Google Scholar]
  62. Kurtovic, N. T., Pinilla, P., Long, F., et al. 2021, A&A, 645, A139 [EDP Sciences] [Google Scholar]
  63. Laos, S., Wisniewski, J. P., Kuchner, M. J., et al. 2022, ApJ, 935, 111 [Google Scholar]
  64. Lau, T. C. H., Drazkowska, J., Stammler, S. M., Birnstiel, T., & Dullemond, C. P. 2022, A&A, 668, A170 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  65. Lavail, A., Kochukhov, O., & Hussain, G. A. J. 2019, A&A, 630, A99 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  66. Liebhart, A., Güdel, M., Skinner, S. L., & Green, J. 2014, A&A, 570, L11 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  67. Liffman, K. 2003, PASA, 20, 337 [NASA ADS] [CrossRef] [Google Scholar]
  68. Lin, Y., Wang, L., Fang, M., Nemer, A., & Goodman, J. 2025, ApJ, 981, 180 [Google Scholar]
  69. Lynden-Bell, D., & Pringle, J. E. 1974, MNRAS, 168, 603 [Google Scholar]
  70. Mah, J., Bitsch, B., Pascucci, I., & Henning, T. 2023, A&A, 677, L7 [CrossRef] [EDP Sciences] [Google Scholar]
  71. Manara, C. F., Ansdell, M., Rosotti, G. P., et al. 2023, in ASP Conf. Ser., 534, 539 [NASA ADS] [Google Scholar]
  72. Mann, R. K., Andrews, S. M., Eisner, J. A., et al. 2015, ApJ, 802, 77 [Google Scholar]
  73. Ment, K., & Charbonneau, D. 2023, AJ, 165, 265 [NASA ADS] [CrossRef] [Google Scholar]
  74. Monsch, K., Picogna, G., Ercolano, B., & Kley, W. 2021, A&A, 646, A169 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  75. Mueller, M. A., Johns-Krull, C. M., Stassun, K. G., & Dixon, D. M. 2024, ApJ, 964, 1 [Google Scholar]
  76. Murphy, S. J., Mamajek, E. E., & Bell, C. P. M. 2018, MNRAS, 476, 3290 [NASA ADS] [CrossRef] [Google Scholar]
  77. Muzerolle, J., Flaherty, K., Balog, Z., et al. 2009, ApJ, 704, L15 [Google Scholar]
  78. Najita, J. R., Carr, J. S., Pontoppidan, K. M., et al. 2013, ApJ, 766, 134 [Google Scholar]
  79. Nelson, R. P. 2018, in Handbook of Exoplanets, eds. Deeg, H. J., & Belmonte, J. A. (Berlin: Springer), 139 [Google Scholar]
  80. Nemer, A., & Goodman, J. 2024, ApJ, 961, 122 [Google Scholar]
  81. Ostriker, E. C., & Shu, F. H. 1995, ApJ, 447, 813 [NASA ADS] [CrossRef] [Google Scholar]
  82. Pascucci, I., Banzatti, A., Gorti, U., et al. 2020, ApJ, 903, 78 [Google Scholar]
  83. Pascucci, I., Beck, T. L., Cabrit, S., et al. 2024, Nat.Astron., 9, 81 [Google Scholar]
  84. Pfalzner, S., & Dincer, F. 2024, ApJ, 963, 122 [Google Scholar]
  85. Pinilla, P., Pascucci, I., & Marino, S. 2020, A&A, 635, A105 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  86. Pinilla, P., Lenz, C. T., & Stammler, S. M. 2021, A&A, 645, A70 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  87. Ragossnig, F., Dorfi, E. A., Ratschiner, B., et al. 2020, Comp. Phys. Commun., 256, 107437 [Google Scholar]
  88. Reiners, A., Shulyak, D., Käpylä, P. J., et al. 2022, A&A, 662, A41 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  89. Ribas, Á., Bouy, H., & Merín, B. 2015, A&A, 576, A52 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  90. Robrade, J., & Schmitt, J. H. M. M. 2007, A&A, 473, 229 [CrossRef] [EDP Sciences] [Google Scholar]
  91. Roggero, N., Bouvier, J., Rebull, L. M., & Cody, A. M. 2021, A&A, 651, A44 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  92. Romanova, M. M., & Owocki, S. P. 2015, Space Sci. Rev., 191, 339 [Google Scholar]
  93. Romanova, M. M., Ustyugova, G. V., Koldoba, A. V., & Lovelace, R. V. E. 2002, ApJ, 578, 420 [Google Scholar]
  94. Schneider, P. C., France, K., Günther, H. M., et al. 2015, A&A, 584, A51 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  95. Sellek, A. D., Grassi, T., Picogna, G., et al. 2024, A&A, 690, A296 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  96. Serna, J., Pinzón, G., Hernández, J., et al. 2024, ApJ, 968, 68 [NASA ADS] [CrossRef] [Google Scholar]
  97. Silverberg, S. M., Wisniewski, J. P., Kuchner, M. J., et al. 2020, ApJ, 890, 106 [NASA ADS] [CrossRef] [Google Scholar]
  98. Skinner, S. L., Güdel, M., Briggs, K. R., & Lamzin, S. A. 2010, ApJ, 722, 1654 [Google Scholar]
  99. Smith, G. D., Gillen, E., Hodgkin, S. T., et al. 2023, MNRAS, 523, 169 [Google Scholar]
  100. Stahler, S. W. 1988, ApJ, 332, 804 [NASA ADS] [CrossRef] [Google Scholar]
  101. Steiner, D., Gehrig, L., Ratschiner, B., et al. 2021, A&A, 655, A110 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  102. Strogatz, S. H. 2000, Nonlinear Dynamics and Chaos: With Applications to Physics, Biology, Chemistry and Engineering (Colorado: Westview Press) [Google Scholar]
  103. Tabone, B., Bettoni, G., van Dishoeck, E. F., et al. 2023, Nat.Astron., 7, 805 [Google Scholar]
  104. Takasao, S., Tomida, K., Iwasaki, K., & Suzuki, T. K. 2022, ApJ, 941, 73 [NASA ADS] [CrossRef] [Google Scholar]
  105. Tambovtseva, L. V., & Grinin, V. P. 2008, Astron. Lett., 34, 231 [Google Scholar]
  106. Telleschi, A., Güdel, M., Briggs, K. R., Audard, M., & Palla, F. 2007, A&A, 468, 425 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  107. Thanathibodee, T., Molina, B., Serna, J., et al. 2023, ApJ, 944, 90 [NASA ADS] [CrossRef] [Google Scholar]
  108. van der Marel, N., & Mulders, G. D. 2021, AJ, 162, 28 [NASA ADS] [CrossRef] [Google Scholar]
  109. Watson, D. M., Calvet, N. P., Fischer, W. J., et al. 2016, ApJ, 828, 52 [NASA ADS] [CrossRef] [Google Scholar]
  110. Weber, M. L., Picogna, G., & Ercolano, B. 2024, A&A, 686, A53 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  111. Winter, A. J., & Haworth, T. J. 2022, Euro. Phys. J. Plus, 137, 1132 [NASA ADS] [CrossRef] [Google Scholar]
  112. Wright, J. T., Fakhouri, O., Marcy, G. W., et al. 2011, PASP, 123, 412 [NASA ADS] [CrossRef] [Google Scholar]
  113. Xie, C., Pascucci, I., Long, F., et al. 2023, ApJ, 959, L25 [NASA ADS] [CrossRef] [Google Scholar]
  114. Zhu, Z., Stone, J. M., & Calvet, N. 2024, MNRAS, 528, 2883 [NASA ADS] [CrossRef] [Google Scholar]
  115. Zurlo, A., Cugno, G., Montesinos, M., et al. 2020, A&A, 633, A119 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  116. Zurlo, A., Cieza, L. A., Ansdell, M., et al. 2021, MNRAS, 501, 2305 [NASA ADS] [CrossRef] [Google Scholar]

1

Winds may actually be more radially extended (e.g., Pascucci et al. 2024) but to first order this does not affect the column density.

2

Infall during the Class I phase will shield the disk from X-rays at earlier times, regardless of the inner disk accretion mode.

3

For a discussion of the complexities of pre-MS model time for young stellar objects, see Sect. 4 of Doppmann et al. (2005).

4

Embedded giant planets could also influence the opening of a gap by a PE wind (Weber et al. 2024).

All Figures

thumbnail Fig. 1

Four scenarios of disk evolution. (A) With M* = 1 M, B* = 0.44 kG, and P* = 3 days, M˙acc$\dot{M}_{\rm acc}$ exceeds M˙crit$\dot{M}_{\rm crit}$ at 2.3 Myr, Rmag retreats inside Rcor, an IDW ceases, and stellar X-rays open a disk gap. (B) In the case of a more rapidly rotating (P* = 2 days) star with a stronger (0.6 kG) field, M˙acc$\dot{M}_{\rm acc}$ never exceeds M˙crit$\dot{M}_{\rm crit}$ and no gap opens. (C) For a more slowly rotating (P* = 7 day) star, there is no IDW but lower LX and PE delays tgap to 5.7 Myr. (D) In the 0.3 M case, lower M˙acc$\dot{M}_{\rm acc}$ remains below M˙crit$\dot{M}_{\rm crit}$ for longer, delaying tgap to 13.3 Myr.

In the text
thumbnail Fig. 2

tgap (Myr) as a function of B* and P* for M* = 1 M. Labeled stars correspond to scenarios A, B, and C in Fig. 1 In the shaded region, M˙acc$\dot{M}_{\rm acc}$ remains below M˙crit$\dot{M}_{\rm crit}$ and no gap opens for the entire simulation.

In the text
thumbnail Fig. 3

Same as Fig. 2, but with tgap (Myr) as a function of M* and P* for B* = 0.44 kG. Stars correspond to scenarios A, C, and D in Fig. 1.

In the text
thumbnail Fig. 4

Threshold accretion rate below which a wind cannot effectively shield the disk from PE-driving X-ray radiation, vs. stellar mass. Curves are calculations for different rotation periods and ages of 1–10 Myr, and points are values of M˙acc$\dot{M}_{\rm acc}$ from Manara et al. (2023).

In the text
thumbnail Fig. 5

Panel a: above a critical value of the product of the X-ray attenuation coefficient, a, times the ejection/accretion ratio, f (Eq. (10b)), solutions for the inner disk accretion rate multifurcate. Panel b: domains of solutions to Eqs. (10a) and (10b) in af˙$a \dot f$ vs. the maximum PE rate. Contours are the ratio of the maximum to minimum values of M˙acc$\dot{M}_{\rm acc}$ expressed in dex. The dotted line is the parameter range for the case shown in the top panel.

In the text
thumbnail Fig. 6

Illustration of two scenarios of disk evolution and planet formation. Left: if Rmag < Rcor, unimpeded X-ray irradiation by the central star opens a gap in the disk that eventually depletes the inner disk and stymies planet formation there. A pressure bump outside the cavity can trap solids and enable the formation of one or more giant planets. Right: if Rmag > Rcor, an IDW shields the disk further out from X-rays, slowing or halting PE and delaying or preventing the opening of a gap. Solids drift accumulate closer to the star, leading to more compact dust disks and close-in multi-planet systems.

In the text
thumbnail Fig. A.1

Disk lifetime vs. B* and initial (not fixed) P* for a 1-M star predicted by the numerical TAPIR star-disk interaction model. Compare to Fig. 2.

In the text
thumbnail Fig. A.2

Same as Fig. A.1 but for a 0.3-M star.

In the text
thumbnail Fig. B.1

tgap vs. B* and P* for 1 M. Same as Fig. 2 except using the stellar magnetosphere (Königl 1991) (top) and trapped flux (Ostriker & Shu 1995) models to calculate evolution of P*. The ordinate is now the minimum rotation period over the disk lifetime.

In the text
thumbnail Fig. B.2

Stellar radius evolution for 1- and 0.3-Mstars predicted by the Dartmouth standard model (Dotter 2016), the Dartmouth magnetic model , which includes the inflationary effect of magnetic fields via spots and inhibition of convection (Feiden 2016), and the Baraffe et al. (2015) model. In the Dartmouth models a different set of abundances and opacities results in prolonged deuterium burning compared to the Baraffe et al. (2015) model (G. Feiden, pers. comm.). This delays stellar contraction and potential affects disk evolution via the radius dependence of Eqn. 4.

In the text
thumbnail Fig. B.3

tgap vs. B* and P* for 1 M. Same as Fig. 2 except with a trebled initial accretion rate (top) or an accretion rate power law index of −2 instead of −5/4.

In the text
thumbnail Fig. B.4

0.1–1 keV X-ray attenuation at ℛ = 1 vs. stellar mass at the end of the IDW phase for different values of ejection/accretion ratio f.

In the text
thumbnail Fig. B.5

Attenuation of 0.1–1 keV X-rays for different coronal emission models.

In the text

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.