Low EUV Luminosities Impinging on Protoplanetary Disks

The amount of high-energy stellar radiation reaching the surface of protoplanetary disks is essential to determine their chemistry and physical evolution. Here, we use millimetric and centimetric radio data to constrain the EUV luminosity impinging on 14 disks around young (~2-10Myr) sun-like stars. For each object we identify the long-wavelength emission in excess to the dust thermal emission, attribute that to free-free disk emission, and thereby compute an upper limit to the EUV reaching the disk. We find upper limits lower than 10$^{42}$ photons/s for all sources without jets and lower than $5 \times 10^{40}$ photons/s for the three older sources in our sample. These latter values are low for EUV-driven photoevaporation alone to clear out protoplanetary material in the timescale inferred by observations. In addition, our EUV upper limits are too low to reproduce the [NeII] 12.81 micron luminosities from three disks with slow [NeII]-detected winds. This indicates that the [NeII] line in these sources primarily traces a mostly neutral wind where Ne is ionized by 1 keV X-ray photons, implying higher photoevaporative mass loss rates than those predicted by EUV-driven models alone. In summary, our results suggest that high-energy stellar photons other than EUV may dominate the dispersal of protoplanetary disks around sun-like stars.Comment: Accepted for publication to The Astrophysical Journa

Tracing Slow Winds from T Tauri Stars via Low Velocity Forbidden Line Emission

Using Keck/HIRES spectra {\Delta}v ~ 7 km/s, we analyze forbidden lines of [O I] 6300 {\AA}, [O I] 5577 {\AA} and [S II] 6731 {\AA} from 33 T Tauri stars covering a range of disk evolutionary stages. After removing a high velocity component (HVC) associated with microjets, we study the properties of the low velocity component (LVC). The LVC can be attributed to slow disk winds that could be magnetically (MHD) or thermally (photoevaporative) driven. Both of these winds play an important role in the evolution and dispersal of protoplanetary material. LVC emission is seen in all 30 stars with detected [O I] but only in 2 out of eight with detected [S II] , so our analysis is largely based on the properties of the [O I] LVC. The LVC itself is resolved into broad (BC) and narrow (NC) kinematic components. Both components are found over a wide range of accretion rates and their luminosity is correlated with the accretion luminosity, but the NC is proportionately stronger than the BC in transition disks. The FWHM of both the BC and NC correlates with disk inclination, consistent with Keplerian broadening from radii of 0.05 to 0.5 AU and 0.5 to 5 AU, respectively. The velocity centroids of the BC suggest formation in an MHD disk wind, with the largest blueshifts found in sources with closer to face-on orientations. The velocity centroids of the NC however, show no dependence on disk inclination. The origin of this component is less clear and the evidence for photoevaporation is not conclusive

Emission Lines from the Gas Disk Around TW Hydra and the Origin of the Inner Hole

We compare line emission calculated from theoretical disk models with optical to submillimeter wavelength observational data of the gas disk surrounding TW Hya and infer the spatial distribution of mass in the gas disk. The model disk that best matches observations has a gas mass ranging from approx.10(exp 4) to 10(exp 5) M for 0.06AU 13.6 eV) flux from TW Hya. H2 pure rotational line emission comes primarily from r approx. 1 to 30 AU. [Oi] 63microns, HCO+, and CO pure rotational lines all arise from the outer disk at r approx. 30-120 AU. We discuss planet formation and photoevaporation as causes for the decrease in surface density of gas and dust inside 4 AU. If a planet is present, our results suggest a planet mass approx. 4-7MJ situated at 3 AU. Using our photoevaporation models and the best surface density profile match to observations, we estimate a current photoevaporative mass loss rate of 4x10(exp 9M)/yr and a remaining disk lifetime of approx.5 million years

The Photoevaporative Wind from the Disk of TW Hya

Photoevaporation driven by the central star is expected to be a ubiquitous and important mechanism to disperse the circumstellar dust and gas from which planets form. Here, we present a detailed study of the circumstellar disk surrounding the nearby star TW Hya and provide observational constraints to its photoevaporative wind. Our new high-resolution (R ~ 30,000) mid-infrared spectroscopy in the [Ne II] 12.81 {\mu}m line confirms that this gas diagnostic traces the unbound wind component within 10AU from the star. From the blueshift and asymmetry in the line profile, we estimate that most (>80%) of the [Ne II] emission arises from disk radii where the midplane is optically thick to the redshifted outflowing gas, meaning beyond the 1 or 4AU dust rim inferred from other observations. We re-analyze high-resolution (R ~ 48, 000) archival optical spectra searching for additional transitions that may trace the photoevaporative flow. Unlike the [Ne II] line, optical forbidden lines from OI, SII, and MgI are centered at the stellar velocity and have symmetric profiles. The only way these lines could trace the photoevaporative flow is if they arise from a disk region physically distinct from that traced by the [Ne II] line, specifically from within the optically thin dust gap. However, the small (~10 km/s) FWHM of these lines suggest that most of the emitting gas traced at optical wavelengths is bound to the system rather than unbound. We discuss the implications of our results for a planet-induced versus a photoevaporation-induced gap.Comment: Accepted for publication in Ap

Narrow Na and K Absorption Lines Toward T Tauri Stars: Tracing the Atomic Envelope of Molecular Clouds

We present a detailed analysis of narrow Na i and K i absorption resonance lines toward nearly 40 T Tauri stars in Taurus with the goal of clarifying their origin. The Na i λ5889.95 line is detected toward all but one source, while the weaker K i λ7698.96 line is detected in about two-thirds of the sample. The similarity in their peak centroids and the significant positive correlation between their equivalent widths demonstrate that these transitions trace the same atomic gas. The absorption lines are present toward both disk and diskless young stellar objects, which excludes cold gas within the circumstellar disk as the absorbing material. A comparison of Na i and CO detections and peak centroids demonstrates that the atomic gas and molecular gas are not co-located, the atomic gas being more extended than the molecular gas. The width of the atomic lines corroborates this finding and points to atomic gas about an order of magnitude warmer than the molecular gas. The distribution of Na i radial velocities shows a clear spatial gradient along the length of the Taurus molecular cloud filaments. This suggests that absorption is associated with the Taurus molecular cloud. Assuming that the gradient is due to cloud rotation, the rotation of the atomic gas is consistent with differential galactic rotation, whereas the rotation of the molecular gas, although with the same rotation axis, is retrograde. Our analysis shows that narrow Na i and K i absorption resonance lines are useful tracers of the atomic envelope of molecular clouds. In line with recent findings from giant molecular clouds, our results demonstrate that the velocity fields of the atomic and molecular gas are misaligned. The angular momentum of a molecular cloud is not simply inherited from the rotating Galactic disk from which it formed but may be redistributed by cloud–cloud interactions

Probing stellar accretion with mid-infrared hydrogen lines

In this paper we investigate the origin of the mid-infrared (IR) hydrogen recombination lines for a sample of 114 disks in different evolutionary stages (full, transitional and debris disks) collected from the {\it Spitzer} archive. We focus on the two brighter {H~{\sc i}} lines observed in the {\it Spitzer} spectra, the {H~{\sc i}}(7-6) at 12.37$\mu$m and the {H~{\sc i}}(9-7) at 11.32$\mu$m. We detect the {H~{\sc i}}(7-6) line in 46 objects, and the {H~{\sc i}}(9-7) in 11. We compare these lines with the other most common gas line detected in {\it Spitzer} spectra, the {[Ne~{\sc iii}]} at 12.81$\mu$m. We argue that it is unlikely that the {H~{\sc i}} emission originates from the photoevaporating upper surface layers of the disk, as has been found for the {[Ne~{\sc iii}]} lines toward low-accreting stars. Using the {H~{\sc i}}(9-7)/{H~{\sc i}}(7-6) line ratios we find these gas lines are likely probing gas with hydrogen column densities of 10$^{10}$-10$^{11}$~cm$^{-3}$. The subsample of objects surrounded by full and transitional disks show a positive correlation between the accretion luminosity and the {H~{\sc i}} line luminosity. These two results suggest that the observed mid-IR {H~{\sc i}} lines trace gas accreting onto the star in the same way as other hydrogen recombination lines at shorter wavelengths. A pure chromospheric origin of these lines can be excluded for the vast majority of full and transitional disks.We report for the first time the detection of the {H~{\sc i}}(7-6) line in eight young (< 20~Myr) debris disks. A pure chromospheric origin cannot be ruled out in these objects. If the {H~{\sc i}}(7-6) line traces accretion in these older systems, as in the case of full and transitional disks, the strength of the emission implies accretion rates lower than 10$^{-10}$M$_{\odot}$/yr. We discuss some advantages of extending accretion indicators to longer wavelengths

Diagnostic Line Emission from EUV and X-ray Illuminated Disks and Shocks around Low Mass stars

Extreme ultraviolet (EUV, 13.6 eV < h\nu \lta 100 eV) and X-rays in the 0.1-2 keV band can heat the surfaces of disks around young, low mass stars to thousands of degrees and ionize species with ionization potentials greater than 13.6 eV. Shocks generated by protostellar winds can also heat and ionize the same species close to the star/disk system. These processes produce diagnostic lines (e.g., [NeII] 12.8 $\mu$m and [OI] 6300 \AA) that we model as functions of key parameters such as EUV luminosity and spectral shape, X-ray luminosity and spectral shape, and wind mass loss rate and shock speed. Comparing our models with observations, we conclude that either internal shocks in the winds or X-rays incident on the disk surfaces often produce the observed [NeII] line, although there are cases where EUV may dominate. Shocks created by the oblique interaction of winds with disks are unlikely [NeII] sources because these shocks are too weak to ionize Ne. Even if [NeII] is mainly produced by X-rays or internal wind shocks, the neon observations typically place upper limits of \lta 10^{42} s$^{-1}$ on the EUV photon luminosity of these young low mass stars. The observed [OI] 6300 \AA line has both a low velocity component (LVC) and a high velocity component. The latter likely arises in internal wind shocks. For the former we find that X-rays likely produce more [OI] luminosity than either the EUV layer, the transition layer between the EUV and X-ray layer, or the shear layer where the protostellar wind shocks and entrains disk material in a radial flow across the surface of the disk. Our soft X-ray models produce [OI] LVCs with luminosities up to $10^{-4}$ L$_\odot$, but may not be able to explain the most luminous LVCs.Comment: 51 pages, 10 figures, accepted to Ap

Tracing Slow Winds From T Tauri Stars via Low-velocity Forbidden Line Emission

Using Keck/HIRES spectra (Δ v ∼ 7 km s−1) we analyze forbidden lines of [O I] 6300 Å, [O I] 5577 Å and [S II] 6731 Å from 33 T Tauri stars covering a range of disk evolutionary stages. After removing a high-velocity component (HVC) associated with microjets, we study the properties of the low-velocity component (LVC). The LVC can be attributed to slow disk winds that could be magnetically (magnetohydrodynamic) or thermally (photoevaporative) driven. Both of these winds play an important role in the evolution and dispersal of protoplanetary material. LVC emission is seen in all 30 stars with detected [O I] but only in two out of eight with detected [S II], so our analysis is largely based on the properties of the [O I] LVC. The LVC itself is resolved into broad (BC) and narrow (NC) kinematic components. Both components are found over a wide range of accretion rates and their luminosity is correlated with the accretion luminosity, but the NC is proportionately stronger than the BC in transition disks. The full width at half maximum of both the BC and NC correlates with disk inclination, consistent with Keplerian broadening from radii of 0.05 to 0.5 au and 0.5 to 5 au, respectively. The velocity centroids of the BC suggest formation in an MHD disk wind, with the largest blueshifts found in sources with closer to face-on orientations. The velocity centroids of the NC, however, show no dependence on disk inclination. The origin of this component is less clear and the evidence for photoevaporation is not conclusive

Parallel Mapper

The construction of Mapper has emerged in the last decade as a powerful and effective topological data analysis tool that approximates and generalizes other topological summaries, such as the Reeb graph, the contour tree, split, and joint trees. In this paper, we study the parallel analysis of the construction of Mapper. We give a provably correct parallel algorithm to execute Mapper on multiple processors and discuss the performance results that compare our approach to a reference sequential Mapper implementation. We report the performance experiments that demonstrate the efficiency of our method