Disk Photoevaporation Research Articles

Broadening the Canonical Picture of EUV-driven Photoevaporation of Accretion Disks

Photoevaporation driven by hydrogen-ionizing extreme-ultraviolet (EUV) radiation profoundly shapes the lives of diverse astrophysical objects. We construct an analytical model accounting for the finite timescales of photoheating and photoionization and apply it to the dispersal of protoplanetary disks. The model yields improved estimates for the ionization, temperature, and velocity versus distance from the central source when compared to the classical picture of fully ionized and isothermal winds with temperatures ≈104 K and speeds ≈10 km s−1. In contrast to the classical picture, the photoevaporative winds take on several distinct hydrodynamical and thermochemical states depending on the central star’s EUV emission rate and spectral hardness: T Tauri stars with EUV luminosities around 1030 erg s−1 drive nonisothermal ionized disk winds at lower temperatures than the classical value if the spectrum is soft, with an average deposited energy per photoionization less than about 3.7 eV. If, however, the spectrum is hard, the winds tend to be atomic and isothermal at most disk radii. For lower EUV intensities, even with soft spectra, atomic winds can emerge beyond ∼10 au through advection. We show that these predictions are in general agreement with detailed radiation hydrodynamics calculations. The model furthermore illustrates how the energy efficiency of photoevaporation varies with the intensity and spectral hardness of the EUV illumination, as well as addressing discrepancies in the literature around the effectiveness of X-ray photoevaporation. The findings highlight the importance of the photoheating and photoionization timescales both for modeling and for understanding winds’ observed behavior.

Photoevaporation of protoplanetary discs with PLUTO+PRIZMO I. Lower X-ray--driven mass-loss rates due to enhanced cooling

Photoevaporation is an important process for protoplanetary disc dispersal, but there has so far been a lack of consensus from simulations over the mass-loss rates and the most important part of the high-energy spectrum involved in driving the wind. We aim to isolate the origins of these discrepancies through carefully benchmarked hydrodynamic simulations of X-ray photoevaporation with time-dependent thermochemistry calculated on the fly. We conducted hydrodynamic simulations with pluto where the thermochemistry is calculated using prizmo . We explored the contribution of certain key microphysical processes and the impact of employing different spectra previously used in literature studies. We find that additional cooling results from the excitation of O by neutral H, which leads to dramatically reduced mass-loss across the disc compared to previous X-ray photoevaporation models, with an integrated rate of $ M sun yr^ $. Such rates would allow for longer-lived discs than previously expected from population synthesis. An alternative spectrum with less soft X-ray produces mass-loss rates around a factor of two to three times lower. The chemistry is significantly out of equilibrium, with the survival of H into the wind being aided by advection. This leads to H becoming the dominant coolant at 10s au, thus stabilising a larger radial temperature gradient across the wind as well as providing a possible wind tracer.

The need for spatially resolved observations of polycyclic aromatic hydrocarbons in protoplanetary discs

Context. The signatures of polycyclic aromatic hydrocarbons (PAHs) have been observed in protoplanetary discs, and their emission features obtained from spectral energy distributions (SED) have been used in the literature to characterise their size and determine their abundance. Aims. Two simple disc models (uniform PAH distribution against a PAH gap in the inner disc) are compared to investigate the difference of their SED and obtainable information. Methods. We used the radiative transfer code RADMC-3D to model the SED of two protoplanetary discs orbiting a typical Herbig star, one of which features a depletion of PAHs in the inner disc. We further created artificial images of the discs at face-on view to extract radial profiles of the PAH emission in the infrared. Results. We find that the extracted PAH features from an SED provide limited information about the PAHs in protoplanetary disc environments, except for the ionisation state. The distribution of PAHs in a protoplanetary disc influences the total observed PAH luminosity in a non-linear fashion and alters the relative strength between the 3.3 µm and 11.3 µm features. Furthermore, we produced radial profiles at the 3 µm, 6 µm and, 11 µm PAH emission features and find that they follow a double power-law profile where the slope reflects the radiative environment (single photon regime vs. multi-photon regime) in which the PAHs lie. Conclusions. Using spatially resolved techniques such as IFU or imaging in the era of the James Webb Space Telescope, we find that multi-wavelength radial emission profiles will not only provide information on the spatial distribution of the PAHs, but may also provide information on their size and underlying UV environment, which is crucial for photo-evaporative disc wind models.

The interplay between forming planets and photoevaporating discs

Context. Disc winds and planet–disc interactions are two crucial mechanisms that define the structure, evolution, and dispersal of protoplanetary discs. While winds are capable of removing material from discs, eventually leading to their dispersal, massive planets can shape their disc by creating sub-structures such as gaps and spiral arms. Aims. We studied the interplay between an X-ray photoevaporative disc wind and the sub-structures generated due to planet–disc interactions to determine how their mutual interactions affect the disc’s and the planet’s evolution. Methods. We performed 3D hydrodynamic simulations of viscous discs (α = 6.9 × 10−4) that host a Jupiter-like planet and undergo X-ray photoevaporation. We traced the gas flows within the disc and wind and measured the rate of accretion onto the planet, as well as the gravitational torque that is acting on it. Results. Our results show that the planetary gap removes the wind’s pressure support, allowing wind material to fall back into the gap. This opens new pathways for material from the inner disc (and part of the outer disc) to be redistributed through the wind towards the gap. Consequently, the gap becomes shallower and the flow of mass across the gap in both directions is significantly increased, as is the planet’s mass-accretion rate (by factors of ≈5 and ≈2, respectively). Moreover, the wind-driven redistribution results in a denser inner disc and a less dense outer disc, which, combined with the recycling of a significant portion of the inner wind, leads to longer lifetimes for the inner disc, contrary to the expectation in a planet-induced photoevaporation scenario that has been proposed in the past.

Constraining Free–Free Emission and Photoevaporative Mass-loss Rates for Known Proplyds and New VLA–identified Candidate Proplyds in NGC 1977

We present NSF's Karl G. Jansky Very Large Array observations covering the NGC 1977 region at 3.0, 6.4, and 15.0 GHz. We search for compact radio sources and detect continuum emission from 34 NGC 1977 cluster members and 37 background objects. Of the 34 radio-detected cluster members, 3 are associated with known proplyds in NGC 1977, 22 are associated with additional young stellar objects in NGC 1977, and 9 are newly identified cluster members. We examine the radio spectral energy distributions, circular polarization, and variability of the detected NGC 1977 sources and identify 10 new candidate proplyds whose radio fluxes are dominated by optically thin free–free emission. We use measurements of free–free emission to calculate the mass-loss rates of known proplyds and new candidate proplyds in NGC 1977, and find values ∼10−9 to 10−8 M ⊙ yr−1, which are lower than the mass-loss rates measured toward proplyds in the Orion Nebula Cluster but consistent with the mass-loss rates predicted by external photoevaporation models for spatially extended disks that are irradiated by the typical external ultraviolet (UV) fields encountered in NGC 1977. Finally, we show that photoevaporative disk winds in NGC 1977 may be illuminated by internal or external sources of ionization, depending on their positions within the cluster. This study provides new constraints on disk properties in a clustered star-forming region with a weaker UV environment than the Orion Nebula Cluster but a stronger UV environment than low-mass star-forming regions like Taurus. Such intermediate UV environments represent the typical conditions of Galactic star and planet formation.

Time-dependent, long-term hydrodynamic simulations of the inner protoplanetary disk. III. The influence of photoevaporation

The final stages of a protoplanetary disk are essential for our understanding of the formation and evolution of planets. Photoevaporation is an important mechanism that contributes to the dispersal of an accretion disk and has significant consequences for the disk's lifetime. However, the combined effects of photoevaporation and star-disk interaction have not been investigated in previous studies. A photoevaporative disk evolution model including a detailed formulation of the inner star-disk interaction region will improve the understanding of the final stages of disk evolution. We combined an implicit disk evolution model with a photoevaporative mass-loss profile. By including the innermost disk regions down to 0.01 AU, we could calculate the star-disk interaction, the stellar spin evolution, and the transition from an accreting disk to the propeller regime self-consistently. Starting from an early Class II star-disk system (with an age of 1 Myr), we calculated the long-term evolution of the system until the disk becomes almost completely dissolved. Photoevaporation has a significant effect on disk structure and evolution. The radial extent of the dead zone decreases, and the number of episodic accretion events (outbursts) is reduced by high stellar X-ray luminosities. Reasonable accretion rates ($> 10^ M_ / yr $) in combination with photoevaporative gaps are possible for a dead zone that is still massive enough to develop episodic accretion events. Furthermore, the stellar spin evolution during the Class II evolution is less affected by the star-disk interaction in the case of high X-ray luminosities. Our results suggest that the formation of planets, especially habitable planets, in the dead zone is strongly impaired in the case of strong X-ray luminosities. Additionally, the importance of the star-disk interaction during the Class II phase with respect to the stellar spin evolution is reduced.

JWST Detects Neon Line Variability in a Protoplanetary Disk

We report the first detection of variability in the mid-infrared neon line emission of a protoplanetary disk by comparing a JWST Mid-InfraRed Instrument Medium Resolution Spectrometer spectrum of SZ Cha taken in 2023 with a Spitzer Infrared Spectrograph Short-High spectrum of this object from 2008. We measure the [Ne iii]-to-[Ne ii] line flux ratio, which is a diagnostic of the high-energy radiation field, to distinguish between the dominance of EUV- or X-ray-driven disk photoevaporation. We find that the [Ne iii]-to-[Ne ii] line flux ratio changes significantly from ∼1.4 in 2008 to ∼0.2 in 2023. This points to a switch from EUV-dominated to X-ray-dominated photoevaporation of the disk. We present contemporaneous ground-based optical spectra of the Hα emission line that show the presence of a strong wind in 2023. We propose that this strong wind prevents EUV radiation from reaching the disk surface while the X-rays permeate the wind and irradiate the disk. We speculate that at the time of the Spitzer observations, the wind was suppressed and EUV radiation reached the disk. These observations confirm that the MIR neon emission lines are sensitive to changes in high-energy radiation reaching the disk surface. This highlights the [Ne iii]-to-[Ne ii] line flux ratio as a tool to gauge the efficiency of disk photoevaporation in order to provide constraints on the planet formation timescale. However, multiwavelength observations are crucial to interpret the observations and properly consider the star–disk connection.

The distribution of accretion rates as a diagnostic of protoplanetary disc evolution

ABSTRACT We show that the distribution of observed accretion rates is a powerful diagnostic of protoplanetary disc physics. Accretion due to turbulent (‘viscous’) transport of angular momentum results in a fundamentally different distribution of accretion rates than accretion driven by magnetized disc winds. We find that a homogeneous sample of ≳300 observed accretion rates would be sufficient to distinguish between these two mechanisms of disc accretion at high confidence, even for pessimistic assumptions. Current samples of T Tauri star accretion rates are not this large, and also suffer from significant inhomogeneity, so both viscous and wind-driven models are broadly consistent with the existing observations. If accretion is viscous, the observed accretion rates require low rates of disc photoevaporation (≲10−9 M⊙ yr−1). Uniform, homogeneous surveys of stellar accretion rates can therefore provide a clear answer to the long-standing question of how protoplanetary discs accrete.

Spatial and dynamical structure of the NGC 2264 star-forming region

Context. The formation of stars within molecular clouds and the early stages of stellar evolution (e.g., mass accretion and disk dispersal) are all active research topics. The target of this study, NGC 2264, is a benchmark star-forming region in which these issues can be profitably studied. Aims. We revisit the structure, dynamics, and star-forming history of NGC 2264 in order to advance our understanding of the processes that lead from molecular clouds to protostars, stellar associations, and the evolution of both. Methods. We assembled a new extensive sample of NGC 2264 members. To this end we used new X-ray data obtained with the XMM-Newton telescope, Gaia eDR3 data, and an extensive collection of public and published catalogs. Following a previous suggestion that the star-forming region might extend significantly beyond the better studied areas, our search covers a wide 2.5×2.5 degrees region in the sky. Results. Our catalog comprises more than 2200 candidate members, which is a ∼100% increase over previous determinations. We analyze their spatial distribution and define new substructures. Using Gaia parallaxes we estimate a new average distance to NGC 2264 of 722±2 pc and suggest that the embedded Spokes subregion is ∼20 pc farther away within the molecular cloud. A complex dynamics is unveiled by the available proper motions and radial velocities: we observe signs of global expansion and rotation. At the same time, we observe the collapse and coalescence of two substructures in a region where active star formation is taking place. The fraction of stars with disks and of those undergoing circumstellar accretion varies significantly across the field, suggesting that star formation has been occurring for several million years. A particularly low accretion disk fraction around the O VII star S Mon might be attributed to external disk photoevaporation or to an older age of the stars in the region. Conclusions. NGC 2264 is not dynamically relaxed and its present configuration is the result of multiple dynamical processes. The cloud has been forming stars for several million years and we identify the process that is likely responsible for the ongoing formation activity.

O i] 6300 Å emission as a probe of external photoevaporation of protoplanetary discs

ABSTRACT We study the utility of the [O i] 6300 Å forbidden line for identifying and interpreting externally driven photoevaporative winds in different environments and at a range of distances. Thermally excited [O i] 6300 Å is a well-known tracer of inner disc winds, so any external contribution needs to be distinguishable. In external winds, the line is not thermally excited and instead results from the dissociation of OH, and we study how the line luminosity resulting from that process scales with the disc/environmental parameters. We find that the line luminosity increases dramatically with FUV radiation field strength above around 5000 G0. The predicted luminosities from our models are consistent with measurements of the line luminosity of proplyds in the Orion Nebula Cluster. The high luminosity in strong UV environments alone may act as a diagnostic, but a rise in the [O i]-to-accretion luminosity ratio is predicted to better separate the two contributions. This could provide a means of identifying external photoevaporation in distant clusters where the proplyd morphology of evaporating discs cannot be spatially resolved.

Monitoring accretion rate variability in the Orion Nebula Cluster with the Wendelstein Wide Field Imager

Context. The understanding of the accretion process has a central role in the understanding of star and planet formation. Aims. We aim to test how accretion variability influences previous correlation analyses of the relation between X-ray activity and accretion rates, which is important for understanding the evolution of circumstellar disks and disk photoevaporation. Methods. We monitored accreting stars in the Orion Nebula Cluster from November 24, 2014, until February 17, 2019, for 42 epochs with the Wendelstein Wide Field Imager in the Sloan Digital Sky Survey u′g′r′ filters on the 2 m Fraunhofer Telescope on Mount Wendelstein. Mass accretion rates were determined from the measured ultraviolet excess. The influence of the mass accretion rate variability on the relation between X-ray luminosities and mass accretion rates was analyzed statistically. Results. We find a typical interquartile range of ∼0.3 dex for the mass accretion rate variability on timescales from weeks to ∼2 yr. The variability has likely no significant influence on a correlation analysis of the X-ray luminosity and the mass accretion rate observed at different times when the sample size is large enough. Conclusions. The observed anticorrelation between the X-ray luminosity and the mass accretion rate predicted by models of photoevaporation-starved accretion is likely not due to a bias introduced by different observing times.

The dispersal of protoplanetary discs – III. Influence of stellar mass on disc photoevaporation

ABSTRACT The strong X-ray irradiation from young solar-type stars may play a crucial role in the thermodynamics and chemistry of circumstellar discs, driving their evolution in the last stages of disc dispersal as well as shaping the atmospheres of newborn planets. In this paper, we study the influence of stellar mass on circumstellar disc mass-loss rates due to X-ray irradiation, extending our previous study of the mass-loss rate’s dependence on the X-ray luminosity and spectrum hardness. We focus on stars with masses between 0.1 and 1 M⊙, which are the main target of current and future missions to find potentially habitable planets. We find a linear relationship between the mass-loss rates and the stellar masses when changing the X-ray luminosity accordingly with the stellar mass. This linear increase is observed also when the X-ray luminosity is kept fixed because of the lower disc aspect ratio which allows the X-ray irradiation to reach larger radii. We provide new analytical relations for the mass-loss rates and profiles of photoevaporative winds as a function of the stellar mass that can be used in disc and planet population synthesis models. Our photoevaporative models correctly predict the observed trend of inner-disc lifetime as a function of stellar mass with an increased steepness for stars smaller than 0.3 M⊙, indicating that X-ray photoevaporation is a good candidate to explain the observed disc dispersal process.

The influence of the environment on the spin evolution of low-mass stars – I. External photoevaporation of circumstellar discs

ABSTRACT Massive stars are strong sources of far-ultraviolet radiation that can be hostile to the evolution of protoplanetary discs, driving mass-loss by external photoevaporation and shortening disc-dissipation time-scales. Their effect may also reduce the time-scale of angular momentum exchanges between the disc and host star during the early pre-main-sequence phase. To improve our understanding of the environmental influence on the rotational history of stars, we developed a model that considers the influence of the local far-ultraviolet radiation on the spin evolution of low mass stars. Our model includes an assumption of disc locking, which fixes the rotation rate during the star-disc-interaction phase, with the duration of this phase parametrized as a function of the local far-ultraviolet radiation and stellar mass (in the range of 0.1–1.3 M⊙). In this way, we demonstrate how the feedback from massive stars can significantly influence the spin evolution of stars and explain the mass dependence observed in period-mass distributions of young regions like Upper Sco and NGC 2264. The high far-ultraviolet environments of high-mass stars can skew the period distribution of surrounding stars towards fast-rotation, explaining the excess of fast-rotating stars in the open cluster h Per. The proposed link between rotation and the pre-main-sequence environment opens new avenues for interpreting the rotational distributions of young stars. For example, we suggest that stellar rotation may be used as a tracer for the primordial ultraviolet irradiation for stars up to ∼1 Gyr, which offers a potential method to connect mature planetary systems to their birth environment.

A slim disc approach to external photoevaporation of discs

ABSTRACT The photoevaporation of protoplanetary discs by nearby massive stars present in their birth cluster plays a vital role in their evolution. Previous modelling assumes that the disc behaves like a classical Keplerian accretion disc out to a radius where the photoevaporative outflow is launched. There is then an abrupt change in the angular velocity profile, and the outflow is modelled by forcing the fluid parcels to conserve their specific angular momenta. Instead, we model externally photoevaporating discs using the slim disc formalism. The slim disc approach self-consistently includes the advection of radial and angular momentum as well as angular momentum redistribution by internal viscous torques. Our resulting models produce a smooth transition from a rotationally supported Keplerian disc to a photoevaporative driven outflow, where this transition typically occurs over ∼4–5 scale heights. The penetration of ultraviolet photons predominately sets the radius of the transition and the viscosity’s strength plays a minor role. By studying the entrainment of dust particles in the outflow, we find a rapid change in the dust size and surface density distribution in the transition region due to the steep gas density gradients present. This rapid change in the dust properties leaves a potentially observable signature in the continuum spectral index of the disc at mm wavelengths. Using the slim disc formalism in future evolutionary calculations will reveal how both the gas and dust evolve in their outer regions and the observable imprints of the external photoevaporation process.

X-Ray Superflares from Pre-main-sequence Stars: Flare Energetics and Frequency

Abstract Solar-type stars exhibit their highest levels of magnetic activity during their early convective pre-main-sequence (PMS) phase of evolution. The most powerful PMS flares, superflares and megaflares, have peak X-ray luminosities of log ( L X ) = 30.5 – 34.0 erg s−1 and total energies log ( E X ) = 34 – 38 erg. Among >24,000 X-ray-selected young (t ≲ 5 Myr) members of 40 nearby star-forming regions from our earlier Chandra MYStIX and SFiNCs surveys, we identify and analyze a well-defined sample of 1086 X-ray superflares and megaflares, the largest sample ever studied. Most are considerably more powerful than optical/X-ray superflares detected on main-sequence stars. This study presents energy estimates of these X-ray flares and the properties of their host stars. These events are produced by young stars of all masses over evolutionary stages ranging from protostars to diskless stars, with the occurrence rate positively correlated with stellar mass. Flare properties are indistinguishable for disk-bearing and diskless stars indicating star–disk magnetic fields are not involved. A slope of α ≃ 2 in the flare energy distributions dN / dE X ∝ E X − α is consistent with that of optical/X-ray flaring from older stars and the Sun. Megaflares ( log ( E X ) > 36.2 erg) from solar-mass stars have an occurrence rate of 1.7 − 0.6 + 1.0 flares/star/year and contribute at least 10%–20% to the total PMS X-ray energetics. These explosive events may have important astrophysical effects on protoplanetary disk photoevaporation, ionization of disk gas, production of spallogenic radionuclides in disk solids, and hydrodynamic escape of young planetary atmospheres. Our following paper details plasma and magnetic loop modeling of the >50 brightest X-ray megaflares.

Testing Photoevaporation and MHD Disk Wind Models through Future High-angular Resolution Radio Observations: The Case of TW Hydrae

Abstract We present theoretical predictions for the free–free emission at centimeter wavelengths obtained from photoevaporation and magnetohydrodynamic (MHD) wind disk models adjusted to the case of the TW Hydrae young stellar object. For this system, disk photoevaporation with heating due to the high-energy photons from the star has been proposed as a possible mechanism to open the gap observed in the dust emission with the Atacama Large Millimeter/submillimeter Array. We show that the photoevaporation disk model predicts a radial profile for the free–free emission that is made of two main spatial components, one originated from the bound disk atmosphere at 0.5–1 au from the star, and another more extended component from the photoevaporative wind at larger disk radii. We also show that the stellar X-ray luminosity has a significant impact on both these components. The predicted radio emission from the MHD wind model has a smoother radial distribution which extends to closer distances to the star than the photoevaporation case. We also show that a future radio telescope such as the Next Generation Very Large Array would have enough sensitivity and angular resolution to spatially resolve the main structures predicted by these models.

The imprint of X-ray photoevaporation of planet-forming discs on the orbital distribution of giant planets

Context. Numerical models have shown that disc dispersal via internal photoevaporation driven by the host star can successfully reproduce the observed pile-up of warm Jupiters near 1–2 au. However, since a range of different mechanisms have been proposed to cause the same feature, clear observational diagnostics of disc dispersal leaving an imprint in the observed distribution of giant planets could help in constraining the dominant mechanisms. Aims. We aim to assess the impact of disc dispersal via X-ray-driven photoevaporation (XPE) on giant planet separations in order to provide theoretical constraints on the location and size of any possible features related to this process within the observed semi-major axis distribution of giant planets. Methods. For this purpose, we perform a set of 1D planet population syntheses with varying initial conditions and correlate the gas giants’ final parking locations with the X-ray luminosities of their host stars in order to quantify observables of this process within the semi-major axis versus host star X-ray luminosity plane of these systems. Results. We find that XPE does create an under-density of gas giants near the gravitational radius, with corresponding pile-ups inside and/or outside this location. However, the size and location of these features are strongly dependent on the choice of initial conditions in our model, such as the assumed formation location of the planets. Conclusions. XPE can strongly affect the migration process of giant planets and leave potentially observable signatures within the observed orbital separations of giant planets. However, due to the simplistic approach employed in our model, which lacks a self-consistent treatment of planet formation within an evolving disc, a quantitative analysis of the final planet population orbits is not possible. Our results, however, should strongly motivate future studies to include realistic disc dispersal mechanisms in global planet population synthesis models with self-consistent planet formation modules.

Testing the models of X-ray driven photoevaporation with accreting stars in the Orion Nebula Cluster

Context. Recent works highlight the importance of stellar X-rays on the evolution of the circumstellar disks of young stellar objects, especially for disk photoevaporation. Aims. A signature of this process may be seen in the so far tentatively observed dependence of stellar accretion rates on X-ray luminosities. According to models of X-ray driven photoevaporation, stars with higher X-ray luminosities should show lower accretion rates, on average, in a sample with similar masses and ages. Methods. To this aim, we have analyzed X-ray properties of young stars in the Orion Nebula Cluster determined with Chandra during the COUP observation as well as accretion data obtained from the photometric catalog of the HST Treasury Program. With these data, we have performed a statistical analysis of the relation between X-ray activity and accretion rates using partial linear regression analysis. Results. The initial anticorrelation found with a sample of 332 young stars is considerably weaker compared to previous studies. However, excluding flaring activity or limiting the X-ray luminosity to the soft band (0.5−2.0 keV) leads to a stronger anticorrelation, which is statistically more significant. Furthermore, we have found a weak positive correlation between the higher component of the plasma temperature gained in the X-ray spectral fitting and the accretion rates, indicating that the hardness of the X-ray spectra may influence the accretion process. Conclusions. There is evidence for a weak anticorrelation, as predicted by theoretical models, suggesting that X-ray photoevaporation modulates the accretion rate through the inner disk at late stages of disk evolution, leading to a phase of photoevaporation-starved accretion.

Exploring the consequences of external photoevaporation of circumstellar disks on the rotational evolution of low mass stars

Exploring the consequences of external photoevaporation of circumstellar disks on the rotational evolution of low mass stars

External Photoevaporation of Disks around Low Mass Young Stellar and Sub-Stellar Objects

External Photoevaporation of Disks around Low Mass Young Stellar and Sub-Stellar Objects