Search

The evolution of protoplanetary disks in regions with massive OB stars is influenced by externally driven winds that deplete the outer parts of disks. These winds have previously been studied via forbidden oxygen emission lines, which also arise in isolated disks in low-mass star forming-regions (SFRs) with weak external UV fields in photoevaporative or magnetic (internal) disk winds. It is crucial to determine how to disentangle external winds from internal ones. Here, we report a proxy for unambiguously identifying externally driven winds with a forbidden line of neutral atomic carbon, [C i] 8727 A. We compare for the first time the spatial location of the emission in the [O i] 5577 A, [O i] 6300 A, and [C i] 8727 A lines traced by VLT/MUSE-NFM, with the ALMA Band 7 continuum disk emission in a sample of 12 proplyds in the Orion Nebula Cluster (ONC). We confirm that the [O i] 5577 A emission is co-spatial with the disk emission, whereas the [O i] 6300 A is emitted both on the disk surface and on the ionization front of the proplyds. We show for the first time that the [C i] 8727 A line is also co-spatial with the disk surface in proplyds, as seen in the MUSE and ALMA data comparison. To verify whether the [C i] 8727 A line is detected in regions where external photoevaporation is not expected, we examine VLT/X-Shooter spectra for young stars in low-mass SFRs. Although the [O i] lines are well detected in all these targets, there is <<10% detection rate in the case of the [C i] 8727 A line. This number increases substantially to a ~40% detection rate in sigma-Orionis, a region with intermediate UV radiation. The spatial location of the [C i] 8727 A line emission and the lack of its detection in isolated disks in low-mass SFRs strongly suggest that this line is a tell-tale tracer of externally driven photoevaporative winds, which agrees with recent excitation models., Comment: Accepted to A&A on October 25, 2024

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 multi-epoch optical spectral monitoring campaign of four CTTSs (TW Hya, RU Lup, BP Tau, and GM Aur) was conducted along with contemporaneous 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-5 R* across all of our observations. There is also significant variability in all emission lines studied, particularly Halpha, Hbeta, and Hgamma. 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., Comment: 39 pages, 10 figures, 12 tables

In the Orion Nebula Cluster (ONC), protoplanetary disks exhibit ionized gas clouds in the form of a striking teardrop shape as massive stars irradiate the disk material. We present the first spatially and spectrally resolved observations of 12 proplyds, using Integral Field Spectroscopy observations performed with the MUSE instrument in Narrow Field Mode (NFM) on the VLT. We present the morphology of the proplyds in seven emission lines and measure the radius of the ionization front (I-front) of the targets in four tracers, covering transitions of different ionization states for the same element. We also derive stellar masses for the targets. The measurements follow a consistent trend of increasing I-front radius for a decreasing strength of the far-UV radiation as expected from photoevaporation models. By analyzing the ratios of the I-front radii as measured in the emission lines of Ha, [OI] 6300, [OII] 7330, and [OIII] 5007, we observe the ionization stratification, that is, the most ionized part of the flow being the furthest from the disk (and closest to the UV source). The ratios of I-front radii scale in the same way for all proplyds in our sample regardless of the incident radiation. We show that the stratification can help constrain the densities near the I-front by using a 1D photoionization model. We derive the upper limits of photoevaporative mass-loss rates by assuming ionization equilibrium, and estimate values decreasing towards lower impinging radiation. We do not find a correlation between Mloss and stellar mass. The highest mass-loss rate is for the proplyd 244-440. These values of Mloss, combined with estimates of the disk mass with ALMA, confirm previous estimates of the short lifetime of these proplyds. This work demonstrates the potential of this dataset and offers a new set of observables to be used to test current and future models of external photoevaporation., Comment: Submitted to A&A on December 18, 2023. Accepted on April 2, 2024

We developed a grid of stellar rotation models for low-mass and solar-type Classical T Tauri stars (CTTS) ($0.3M_{\odot}

Observations of debris disks have significantly improved over the past decades, both in terms of sensitivity and spatial resolution. At near-infrared wavelengths, new observing strategies and post-processing algorithms allow us to drastically improve the final images, revealing faint structures in the disks. These structures inform us about the properties and spatial distribution of the small dust particles. We present new $H$-band observations of the disk around HD 129590, which display an intriguing arc-like structure in total intensity but not in polarimetry, and propose an explanation for the origin of this arc. Assuming geometric parameters for the birth ring of planetesimals, our model provides the positions of millions of particles of different sizes to compute scattered light images. We demonstrate that if the grain size distribution is truncated or strongly peaks at a size larger than the radiation pressure blow-out size we are able to produce an arc quite similar to the observed one. If the birth ring is radially narrow, given that particles of a given size have similar eccentricities, they will have their apocenters at the same distance from the star. Since this is where the particles will spend most of their time, this results in a "apocenter pile-up" that can look like a ring. Due to more efficient forward scattering this arc only appears in total intensity observations and remains undetected in polarimetric data. This scenario requires sharp variations either in the grain size distribution or for the scattering efficiencies $Q_\mathrm{sca}$. Alternative possibilities such as a wavy size distribution and a size-dependent phase function are interesting candidates to strengthen the apocenter pile-up. We also discuss why such arcs are not commonly detected in other systems, which can mainly be explained by the fact that most parent belts are usually broad., Comment: Accepted for publication in A&A, abstract shortened

Absorption features in stellar atmospheres are often used to calibrate photocentric velocities for kinematic analysis of further spectral lines. The Li feature at $\sim$ 6708 {\AA} is commonly used, especially in the case of young stellar objects for which it is one of the strongest absorption lines. However, this is a complex line comprising two isotope fine-structure doublets. We empirically measure the wavelength of this Li feature in a sample of young stars from the PENELLOPE/VLT programme (using X-Shooter, UVES and ESPRESSO data) as well as HARPS data. For 51 targets, we fit 314 individual spectra using the STAR-MELT package, resulting in 241 accurately fitted Li features, given the automated goodness-of-fit threshold. We find the mean air wavelength to be 6707.856 {\AA}, with a standard error of 0.002 {\AA} (0.09 km/s) and a weighted standard deviation of 0.026 {\AA} (1.16 km/s). The observed spread in measured positions spans 0.145 {\AA}, or 6.5 km/s, which is up to a factor of six higher than typically reported velocity errors for high-resolution studies. We also find a correlation between the effective temperature of the star and the wavelength of the central absorption. We discuss how exclusively using this Li feature as a reference for photocentric velocity in young stars could potentially be introducing a systematic positive offset in wavelength to measurements of further spectral lines. If outflow tracing forbidden lines, such as [O i] 6300 {\AA}, are actually more blueshifted than previously thought, this then favours a disk wind as the origin for such emission in young stars., Comment: Accepted for publication in A&A

The vertical structure of debris disks provides clues about their dynamical evolution and the collision rate of the unseen planetesimals. Thanks to the ever-increasing angular resolution of contemporary instruments and facilities, we are beginning to constrain the scale height of a handful of debris disks, either at near-infrared or millimeter wavelengths. Nonetheless, this is often done for individual targets only. We present here the geometric modeling of eight disks close to edge-on, all observed with the same instrument (SPHERE) and using the same mode (dual-beam polarimetric imaging). Motivated by the presence of CO gas in two out of the eight disks, we then investigate the impact that gas can have on the scale height by performing N-body simulations including gas drag and collisions. We show that gas can quickly alter the dynamics of particles (both in the radial and vertical directions), otherwise governed by gravity and radiation pressure. We find that, in the presence of gas, particles smaller than a few tens of microns can efficiently settle toward the midplane at the same time as they migrate outward beyond the birth ring. For second generation gas ($M_\mathrm{gas} \leq 0.1$ $M_\oplus$), the vertical settling should be best observed in scattered light images compared to observations at millimeter wavelengths. But if the gas has a primordial origin ($M_\mathrm{gas} \geq 1$ $M_\oplus$), the disk will appear very flat both at near-infrared and sub-mm wavelengths. Finally, far beyond the birth ring, our results suggest that the surface brightness profile can be as shallow as $\sim -2.25$., Comment: Accepted in MNRAS, 18 pages, 11 Figures

Observations have shown that the disk frequency and the fraction of accreting pre-main-sequence stars decrease with the age of the population and that some stars appear to have disks while their accretion has stopped. Still, it is unclear how disk-bearing stars stop their accretion. To provide insight into the last stages of accretion in low-mass young stars, we conducted a survey of disk-bearing stars that are thought to be non-accretors to identify stars still accreting at very low rates. Here we present the first catalog of the survey of 170 disk-bearing non-accreting stars in Chamaeleon I, Orion OB1, Upper Scorpius, $\gamma$ Velorum, and Upper Centaurus Lupus, using He I $\lambda$10830 as a sensitive probe of accretion. We classify the line profiles into six types and argue that those showing redshifted and/or blueshifted absorption are still accreting. Using these classifications, we found that, among disk-bearing stars previously classified as non-accretors, at least 20-30% are still accreting, with a larger fraction of those at younger population ages. While the difference between the outer disk signature and accretion status is unclear, we find a difference between the inner disk excess and accretion status. There is no preference in the mass of the newly identified accretors, suggesting that the processes inhibiting accretion do not directly depend on mass in the typical mass range of T Tauri stars. Lastly, we found that at a low accretion level, the H$\alpha$ width at the 10% height criteria mischaracterizes a larger fraction of accretors than the line's equivalent width., Comment: Published in AJ. Data tables are available on the Journal's website

One of the most important questions in the field of planet formation is how mm-cm sized dust particles overcome the radial drift and fragmentation barriers to form kilometer-sized planetesimals. ALMA observations of protoplanetary disks, in particular transition disks or disks with clear signs of substructures, can provide new constraints on theories of grain growth and planetesimal formation and therefore represent one possibility to progress on this issue. We here present ALMA band 4 (2.1 mm) observations of the transition disk system Sz 91 and combine them with previously obtained band 6 (1.3 mm) and 7 (0.9 mm) observations. Sz 91 with its well defined mm-ring, more extended gas disk, and evidence of smaller dust particles close to the star, is a clear case of dust filtering and the accumulation of mm sized particles in a gas pressure bump. We computed the spectral index (nearly constant at $\sim$3.34), optical depth (marginally optically thick), and maximum grain size ($\sim\,0.61$ mm) in the dust ring from the multi-wavelength ALMA observations and compared the results with recently published simulations of grain growth in disk substructures. Our observational results are in very good agreement with the predictions of models for grain growth in dust rings that include fragmentation and planetesimal formation through the streaming instability., Comment: Accepted for publication in the ApJ 19 pages, 13 figures

We present a study of the evolution of the inner few astronomical units of protoplanetary disks around low-mass stars. We consider nearby stellar groups with ages spanning from 1 to 11 Myr, distributed into four age bins. Combining PANSTARSS photometry with spectral types, we derive the reddening consistently for each star, which we use (1) to measure the excess emission above the photosphere with a new indicator of IR excess and (2) to estimate the mass accretion rate ($\dot{M}$) from the equivalent width of the H$\alpha$ line. Using the observed decay of $\dot{M}$ as a constrain to fix the initial conditions and the viscosity parameter of viscous evolutionary models, we use approximate Bayesian modeling to infer the dust properties that produce the observed decrease of the IR excess with age, in the range between 4.5 and $24\,\mu$m. We calculate an extensive grid of irradiated disk models with a two-layered wall to emulate a curved dust inner edge and obtain the vertical structure consistent with the surface density predicted by viscous evolution. We find that the median dust depletion in the disk upper layers is $\epsilon \sim 3 \times 10^{-3}$ at 1.5 Myr, consistent with previous studies, and it decreases to $\epsilon \sim 3 \times 10^{-4}$ by 7.5 Myr. We include photoevaporation in a simple model of the disk evolution and find that a photoevaporative wind mass-loss rate of $\sim 1 -3 \times 10 ^{-9} \, M_{\odot}yr^{-1}$ agrees with the decrease of the disk fraction with age reasonably well. The models show the inward evolution of the H$_2$O and CO snowlines., Comment: Accepted for publication in the Astrophysical Journal

We present polarized light observations of the transitional disk around Sz 91 acquired with VLT/NaCo at $H$ (1.7$\mu$m) and $K_s$ (2.2$\mu$m) bands. We resolve the disk and detect polarized emission up to $\sim$0."5 ($\sim$80 au) along with a central cavity at both bands. We computed a radiative transfer model that accounts for the main characteristics of the polarized observations. We found that the emission is best explained by small, porous grains distributed in a disk with a $\sim$45 au cavity. Previous ALMA observations have revealed a large sub-mm cavity ($\sim$83 au) and extended gas emission from the innermost (<16 au) regions up to almost 400 au from the star. Dynamical clearing by multiple low-mass planets arises as the most probable mechanism for the origin of Sz 91's peculiar structure. Using new $L'$ band ADI observations we can rule out companions more massive than $M_p$ $\geq$ 8 $M_\mathrm{Jup}$ beyond 45 au assuming hot-start models. The disk is clearly asymmetric in polarized light along the minor axis, with the north side brighter than the south side. Differences in position angle between the disk observed at sub-mm wavelengths with ALMA and our NaCo observations were found. This suggests that the disk around Sz 91 could be highly structured. Higher signal-to-noise near-IR and sub-mm observations are needed to confirm the existence of such structures and to improve the current understanding in the origin of transitional disks., Comment: 14 pages, 8 figures, accepted for publication in MNRAS

Although the magnetospheric accretion model has been extensively applied to T Tauri Stars with typical mass accretion rates, the very low accretion regime is still not fully explored. Here we report multi-epoch observations and modeling of CVSO 1335, a 5 Myr old solar mass star which is accreting mass from the disk, as evidenced by redshifted absorption in the H$\alpha$ profile, but with very uncertain estimates of mass accretion rate using traditional calibrators. We use the accretion shock model to constrain the mass accretion rate from the Balmer jump excess measured with respect to a non-accreting template, and we model the H$\alpha$ profile, observed simultaneously, using magnetospheric accretion models. Using data taken on consecutive nights, we found that the accretion rate of the star is low, $4-9 \times 10^{-10} \,$ M$_{\odot}\,$ yr$^{-1}$, suggesting a variability on a timescale of days. The observed H$\alpha$ profiles point to two geometrically isolated accretion flows, suggesting a complex infall geometry. The systems of redshifted absorptions observed are consistent with the star being a dipper, although multi-band photometric monitoring is needed to confirm this hypothesis., Comment: Accepted for publication in the Astrophysical Journal. 15 pages, 8 figures, 4 tables

We present a detailed study of the disk around the intermediate mass star SO 411, aiming to explain the spectral energy distribution of this star. We show that this is a transitional disk truncated at $\sim$11 au, with $\sim$0.03 lunar masses of optically thin dust inside the cavity. Gas also flows through the cavity, since we find that the disk is still accreting mass onto the star, at a rate of $\sim 5x10^{-9}$ Msun/yr. Until now, SO 411 has been thought to belong to the $\sim$3 Myr old {$\sigma$} Orionis cluster. However, we analyzed the second Gaia Data Release in combination with kinematic data previously reported, and found that SO 411 can be associated with an sparse stellar population located in front of the {$\sigma$} Orionis cluster. If this is the case, then SO 411 is older and even more peculiar, since primordial disks in this stellar mass range are scarce for ages $>$5 Myr. Analysis of the silicate 10$\mu$m feature of SO 411 indicates that the observed feature arises at the edge of the outer disk, and displays a very high crystallinity ratio of $\sim$0.5, with forsterite the most abundant silicate crystal. The high forsterite abundance points to crystal formation in non-equilibrium conditions. The PAH spectrum of SO 411 is consistent with this intermediate state between the hot and luminous Herbig Ae and the less massive and cooler T Tauri stars. Analysis of the 7.7$\mu$m PAH feature indicates that small PAHs still remain in the SO 411 disk., Comment: Accepted in the Astrophysical Journal (17 pages, 9 Figures)

We report FUV, optical, and NIR observations of three T Tauri stars in the Orion OB1b subassociation with H$\alpha$ equivalent widths consistent with low or absent accretion and various degrees of excess flux in the mid-infrared. We aim to search for evidence of gas in the inner disk in HST ACS/SBC spectra, and to probe the accretion flows onto the star using H$\alpha$ and He I $\lambda$10830 in spectra obtained at the Magellan and SOAR telescopes. At the critical age of 5 Myr, the targets are at different stages of disk evolution. One of our targets is clearly accreting, as shown by redshifted absorption at free-fall velocities in the He I line and wide wings in H$\alpha$; however, a marginal detection of FUV H$_2$ suggests that little gas is present in the inner disk, although the spectral energy distribution indicates that small dust still remains close to the star. Another target is surrounded by a transitional disk, with an inner cavity in which little sub-micron dust remains. Still, the inner disk shows substantial amounts of gas, accreting onto the star at a probably low, but uncertain rate. The third target lacks both a He I line or FUV emission, consistent with no accretion or inner gas disk; its very weak IR excess is consistent with a debris disk. Different processes occurring in targets with ages close to the disk dispersal time suggest that the end of accretion phase is reached in diverse ways., Comment: 13 pages, 3 figures. Accepted for publication in The Astrophysical Journal

We present \emph{Herschel} PACS observations of 8 Classical T Tauri Stars in the $\sim 7-10$ Myr old OB1a and the $\sim 4-5$ Myr old OB1b Orion sub-asscociations. Detailed modeling of the broadband spectral energy distributions, particularly the strong silicate emission at 10 $\mu$m, shows that these objects are (pre)transitional disks with some amount of small optically thin dust inside their cavities, ranging from $\sim 4$ AU to $\sim 90$ AU in size. We analyzed \emph{Spitzer} IRS spectra for two objects in the sample: CVSO-107 and CVSO-109. The IRS spectrum of CVSO-107 indicates the presence of crystalline material inside its gap while the silicate feature of CVSO-109 is characterized by a pristine profile produced by amorphous silicates; the mechanisms creating the optically thin dust seem to depend on disk local conditions. Using millimeter photometry we estimated dust disk masses for CVSO-107 and CVSO-109 lower than the minimum mass of solids needed to form the planets in our Solar System, which suggests that giant planet formation should be over in these disks. We speculate that the presence and maintenance of optically thick material in the inner regions of these pre-transitional disks might point to low-mass planet formation., Comment: 20 pages, 10 figures, 8 tables

We present new Herschel PACS observations of 32 T Tauri stars in the young ($\sim$3 Myr) $\sigma$ Ori cluster. Most of our objects are K & M stars with large excesses at 24 $\mu$m. We used irradiated accretion disk models of D'Alessio et al. (2006) to compare their spectral energy distributions with our observational data. We arrive at the following six conclusions. (i) The observed disks are consistent with irradiated accretion disks systems. (ii) Most of our objects (60%) can be explained by significant dust depletion from the upper disk layers. (iii) Similarly, 61% of our objects can be modeled with large disk sizes ($\rm R_{\rm d} \geq$ 100 AU). (iv) The masses of our disks range between 0.03 to 39 $\rm M_{Jup}$, where 35% of our objects have disk masses lower than 1 Jupiter. Although these are lower limits, high mass ($>$ 0.05 M$_{\odot}$) disks, which are present e.g, in Taurus, are missing. (v) By assuming a uniform distribution of objects around the brightest stars at the center of the cluster, we found that 80% of our disks are exposed to external FUV radiation of $300 \leq G_{0} \leq 1000$, which can be strong enough to photoevaporate the outer edges of the closer disks. (vi) Within 0.6 pc from $\sigma$ Ori we found forbidden emission lines of [NII] in the spectrum of one of our large disk (SO662), but no emission in any of our small ones. This suggests that this object may be an example of a photoevaporating disk., Comment: Accepted by the Astrophysical Journal. 37 pages, 20 figures. Section 6 presents estimations of disks masses and radii and evidence of external photoevaporation

We present a study of 15 new brown dwarfs belonging to the $\sim7$ Myr old 25 Orionis group and Orion OB1a sub-association with spectral types between M6 and M9 and estimated masses between $\sim0.07$M$_\odot$ and $\sim0.01$ M$_\odot$. By comparing them through a Bayesian method with low mass stars ($0.8\lesssim$ M/M$_\odot\lesssim0.1$) from previous works in the 25 Orionis group, we found statistically significant differences in the number fraction of classical T Tauri stars, weak T Tauri stars, class II, evolved discs and purely photospheric emitters at both sides of the sub-stellar mass limit. Particularly we found a fraction of $3.9^{+2.4}_{-1.6}~\%$ low mass stars classified as CTTS and class II or evolved discs, against a fraction of $33.3^{+10.8}_{-9.8}~\%$ in the sub-stellar mass domain. Our results support the suggested scenario in which the dissipation of discs is less efficient for decreasing mass of the central object., Comment: Accepted for publication in Monthly Notices of the Royal Astronomical Society. 15 pages, 9 figures and 5 tables