Source author record

Ilaria Pascucci

Ilaria Pascucci appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

38works
4topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

38 published item(s)

preprint2026arXiv

Protoplanetary disk cavities with JWST-MIRI: a dichotomy in molecular emission

The evolution of planet-forming regions in protoplanetary disks is of fundamental importance to understanding planet formation. Disks with a central deficit in dust emission, a "cavity", have long attracted interest as potential evidence for advanced disk clearing by protoplanets and/or winds. Before JWST, infrared spectra showed that these disks typically lack the strong molecular emission observed in full disks. In this work, we combine a sample of 12 disks with millimeter cavities of a range of sizes ($\sim2$-70 au) and different levels of millimeter and infrared continuum deficits. We analyze their molecular spectra as observed with MIRI on JWST, homogeneously reduced with the new JDISCS pipeline. This analysis demonstrates a stark dichotomy in molecular emission where "molecule-rich" (MR) cavities follow global trends between water, CO, and OH luminosity and accretion luminosity as in full disks, while "molecule-poor" (MP) cavities are significantly sub-luminous in all molecules except sometimes OH. Disk cavities generally show sub-luminous organic emission, higher OH/H$_2$O ratios, and suggest a lower water column density. The sub-thermal excitation of CO and water vibrational lines suggests a decreased gas density in the emitting layer in all cavities, supporting model expectations for C$_2$H$_2$ photodissociation. We discover a bifurcation in infrared index (lower in MR cavities) suggesting that the molecular dichotomy is linked to residual $μ$m-size dust within millimeter disk cavities. Put together, these results suggest a feedback process between dust depletion, gas density decrease, and molecule dissociation. Disk cavities may have a common evolutionary sequence where MR switch into MP over time.

preprint2022arXiv

Determining dispersal mechanisms of protoplanetary disks using accretion and wind mass loss rates

Understanding the origin of accretion and dispersal of protoplanetary disks is fundamental for investigating planet formation. Recent numerical simulations show that launching winds are unavoidable when disks undergo magnetically driven accretion and/or are exposed to external UV radiation. Observations also hint that disk winds are common. We explore how the resulting wind mass loss rate can be used as a probe of both disk accretion and dispersal. As a proof-of-concept study, we focus on magnetocentrifugal winds, MRI (magnetorotational instability) turbulence, and external photoevapotaion. By developing a simple, yet physically motivated disk model and coupling it with simulation results available in the literature, we compute the mass loss rate as a function of external UV flux for each mechanism. We find that different mechanisms lead to different levels of mass loss rate, indicating that the origin of disk accretion and dispersal can be determined, by observing the wind mass loss rate resulting from each mechanism. This determination provides important implications for planet formation. This work thus shows that the ongoing and future observations of the wind mass loss rate for protoplanetary disks are paramount to reliably constrain how protoplanetary disks evolve with time and how planet formation takes place in the disks.

preprint2022arXiv

Pterodactyls: A Tool to Uniformly Search and Vet for Young Transiting Planets In TESS Primary Mission Photometry

Kepler's short-period exoplanet population has revealed evolutionary features such as the Radius Valley and the Hot Neptune desert that are likely sculpted by atmospheric loss over time. These findings suggest that the primordial planet population is different from the Gyr-old Kepler population, and motivates exoplanet searches around young stars. Here, we present pterodactyls , a data reduction pipeline specifically built to address the challenges in discovering exoplanets around young stars and to work with TESS Primary Mission 30-min cadence photometry, since most young stars were not pre-selected TESS 2-min cadence targets. pterodactyls builds on publicly available and tested tools in order to extract, detrend, search, and vet transiting young planet candidates. We search five clusters with known transiting planets: Tucana-Horologium Association, IC 2602, Upper Centaurus Lupus, Ursa Major and Pisces Eridani. We show that pterodactyls recovers seven out of the eight confirmed planets and one out of the two planet candidates, most of which were initially detected in 2-min cadence data. For these clusters, we conduct injection-recovery tests to characterize our detection efficiency, and compute an intrinsic planet occurrence rate of 49+-20% for sub-Neptunes and Neptunes (1.8-6 Re) within 12.5 days, which is higher than Kepler's Gyr-old occurrence rates of 6.8+-0.3%. This potentially implies that these planets have shrunk with time due to atmospheric mass loss. However, a proper assessment of the occurrence of transiting young planets will require a larger sample unbiased to planets already detected. As such, pterodactyls will be used in future work to search and vet for planet candidates in nearby clusters and moving groups.

preprint2022arXiv

The Role of Disk Winds in the Evolution and Dispersal of Protoplanetary Disks

The assembly and architecture of planetary systems strongly depend on the physical processes governing the evolution and dispersal of protoplanetary disks. Since Protostars and Planets VI, new observations and theoretical insights favor disk winds as being one of those key processes. This chapter provides a comprehensive review of recent observations probing outflowing gas launched over a range of disk radii for a wide range of evolutionary stages, enabling an empirical understanding of how winds evolve. In parallel, we review theoretical advancements in both magnetohydrodynamic and photoevaporative disk wind models and identify predictions that can be confronted with observations. By linking theory and observations we critically assess the role of disk winds in the evolution and dispersal of protoplanetary disks. Finally, we explore the impact of disk winds on planet formation and evolution and highlight theoretical work, observations, and critical tests for future progress.

preprint2020arXiv

An Improved HR Diagram for the Orion Trapezium Cluster

In this paper, we present a study of the Trapezium cluster in Orion. We analyze flux-calibrated VLT/MUSE spectra of 361 stars to simultaneously measure the spectral types, reddening, and the optical veiling due to accretion. We find that the extinction law from Cardelli et al. (1989) with a total-to-selective extinction value of $R_{\rm V}=$5.5 is more suitable for this cluster. For 68% of the sample the new spectral types are consistent with literature spectral types within 2 subclasses, but as expected, we derive systematically later types than the literature by one to two subclasses for the sources with significant accretion levels. Here we present an improved Hertzsprung-Russell (H-R) diagram of the Trapezium cluster, in which the contamination by optical veiling on spectral types and stellar luminosities has been properly removed. A comparison of the locations of the stars in the H-R diagram with the non-magnetic and magnetic pre-main sequence evolutionary tracks indicates an age of 1--2~Myr. The magnetic pre-main sequence evolutionary tracks can better explain the luminosities of the low-mass stars. In the H-R diagram, the cluster exhibits a large luminosity spread ($σ$(Log~$L_{\star}/L_{\odot})\sim$0.3). By collecting a sample of 14 clusters/groups with different ages, we find that the luminosity spread tends to be constant ($σ$(Log~$L_{\star}/L_{\odot})\sim$0.2--0.25) after 2~Myr, which suggests that age spread is not the main cause of the spread. There are $\sim$0.1~dex larger luminosity spreads for the younger clusters, e.g., the Trapezium cluster, than the older clusters, which can be explained by the starspots, accretion history and circumstellar disk orientations.

preprint2020arXiv

Dual-Wavelength ALMA Observations of Dust Rings in Protoplanetary Disks

We present new Atacama Large Millimeter/submillimeter Array (ALMA) observations for three protoplanetary disks in Taurus at 2.9\,mm and comparisons with previous 1.3\,mm data both at an angular resolution of $\sim0.''1$ (15\,au for the distance of Taurus). In the single-ring disk DS Tau, double-ring disk GO Tau, and multiple-ring disk DL Tau, the same rings are detected at both wavelengths, with radial locations spanning from 50 to 120\,au. To quantify the dust emission morphology, the observed visibilities are modeled with a parametric prescription for the radial intensity profile. The disk outer radii, taken as 95\% of the total flux encircled in the model intensity profiles, are consistent at both wavelengths for the three disks. Dust evolution models show that dust trapping in local pressure maxima in the outer disk could explain the observed patterns. Dust rings are mostly unresolved. The marginally resolved ring in DS Tau shows a tentatively narrower ring at the longer wavelength, an observational feature expected from efficient dust trapping. The spectral index ($α_{\rm mm}$) increases outward and exhibits local minima that correspond to the peaks of dust rings, indicative of the changes in grain properties across the disks. The low optical depths ($τ\sim$0.1--0.2 at 2.9\,mm and 0.2--0.4 at 1.3\,mm) in the dust rings suggest that grains in the rings may have grown to millimeter sizes. The ubiquitous dust rings in protoplanetary disks modify the overall dynamics and evolution of dust grains, likely paving the way towards the new generation of planet formation.

preprint2020arXiv

Earths in Other Solar Systems N-body simulations: the Role of Orbital Damping in Reproducing the Kepler Planetary Systems

The population of exoplanetary systems detected by Kepler provides opportunities to refine our understanding of planet formation. Unraveling the conditions needed to produce the observed exoplanets will sallow us to make informed predictions as to where habitable worlds exist within the galaxy. In this paper, we examine using N-body simulations how the properties of planetary systems are determined during the final stages of assembly. While accretion is a chaotic process, trends in the ensemble properties of planetary systems provide a memory of the initial distribution of solid mass around a star prior to accretion. We also use EPOS, the Exoplanet Population Observation Simulator, to account for detection biases and show that different accretion scenarios can be distinguished from observations of the Kepler systems. We show that the period of the innermost planet, the ratio of orbital periods of adjacent planets, and masses of the planets are determined by the total mass and radial distribution of embryos and planetesimals at the beginning of accretion. In general, some amount of orbital damping, either via planetesimals or gas, during accretion is needed to match the whole population of exoplanets. Surprisingly, all simulated planetary systems have planets that are similar in size, showing that the "peas in a pod" pattern can be consistent with both a giant impact scenario and a planet migration scenario. The inclusion of material at distances larger than what Kepler observes has a profound impact on the observed planetary architectures, and thus on the formation and delivery of volatiles to possible habitable worlds.

preprint2020arXiv

Hints on the origins of particle traps in protoplanetary disks given by the $M_{\rm{dust}}-M_{\star}$ relation

Demographic surveys of protoplanetary disks, carried out mainly with ALMA, have provided access to a large range of disk dust masses ($M_{\rm{dust}}$) around stars with different stellar types and in different star-forming regions. These surveys found a power-law relation between $M_{\rm{dust}}$ and $M_{\star}$ that steepens in time, but which is also flatter for transition disks (TDs). We performed dust evolution models, which included perturbations to the gas surface density with different amplitudes to investigate the effect of particle trapping on the $M_{\rm{dust}}-M_{\star}$ relation. These perturbations were aimed at mimicking pressure bumps that originated from planets. We focused on the effect caused by different stellar and disk masses based on exoplanet statistics that demonstrate a dependence of planet mass on stellar mass and metallicity. Models of dust evolution can reproduce the observed $M_{\rm{dust}}-M_{\star}$ relation in different star-forming regions when strong pressure bumps are included and when the disk mass scales with stellar mass (case of $M_{\rm{disk}}=0.05\,M_\star$ in our models). This result arises from dust trapping and dust growth beyond centimeter-sized grains inside pressure bumps. However, the flatter relation of $M_{\rm{dust}}-M_{\star}$ for TDs and disks with substructures cannot be reproduced by the models unless the formation of boulders is inhibited inside pressure bumps. In the context of pressure bumps originating from planets, our results agree with current exoplanet statistics on giant planet occurrence increasing with stellar mass, but we cannot draw a conclusion about the type of planets needed in the case of low-mass stars. This is attributed to the fact that for $M_\star<1\,M_\odot$, the observed $M_{\rm{dust}}$ obtained from models is very low due to the efficient growth of dust particles beyond centimeter-sizes inside pressure bumps.

preprint2020arXiv

Pebble-driven Planet Formation around Very Low-mass Stars and Brown Dwarfs

We conduct a pebble-driven planet population synthesis study to investigate the formation of planets around very low-mass stars and brown dwarfs, in the (sub)stellar mass range between $0.01 \ M_{\odot}$ and $0.1 \ M_{\odot}$. Based on the extrapolation of numerical simulations of planetesimal formation by the streaming instability, we obtain the characteristic mass of the planetesimals and the initial masses of the protoplanets (largest bodies from the planetesimal size distributions), in either the early self-gravitating phase or the later non-self-gravitating phase of the protoplanetary disk evolution. We find that the initial protoplanets form with masses that increase with host mass, orbital distance and decrease with disk age. Around late M-dwarfs of $0.1 \ M_{\odot}$, these protoplanets can grow up to Earth-mass planets by pebble accretion. However, around brown dwarfs of $0.01 \ M_{\odot}$, planets do not grow larger than Mars mass when the initial protoplanets are born early in self-gravitating disks, and their growth stalls at around $0.01$ Earth-mass when they are born late in non-self-gravitating disks. Around these low mass stars and brown dwarfs, we find no channel for gas giant planet formation because the solid cores remain too small. When the initial protoplanets form only at the water-ice line, the final planets typically have ${\gtrsim} 15\%$ water mass fraction. Alternatively, when the initial protoplanets form log-uniformly distributed over the entire protoplanetary disk, the final planets are either very water-rich (water mass fraction ${\gtrsim}15\%$) or entirely rocky (water mass fraction ${\lesssim}5\%$).

preprint2020arXiv

Size and structures of disks around very low mass stars in the Taurus star-forming region

We aim to estimate if structures, such as cavities, rings, and gaps, are common in disks around VLMS and to test models of structure formation in these disks. We also aim to compare the radial extent of the gas and dust emission in disks around VLMS, which can give us insight about radial drift. We studied six disks around VLMS in the Taurus star-forming region using ALMA Band 7 ($\sim 340\,$GHz) at a resolution of $\sim0.1''$. The targets were selected because of their high disk dust content in their stellar mass regime. Our observations resolve the disk dust continuum in all disks. In addition, we detect the $^{12}$CO ($J=3-2$) emission line in all targets and $^{13}$CO ($J=3-2$) in five of the six sources. The angular resolution allows the detection of dust substructures in three out of the six disks, which we studied by using UV-modeling. Central cavities are observed in the disks around stars MHO\,6 (M5.0) and CIDA\,1 (M4.5), while we have a tentative detection of a multi-ringed disk around J0433. Single planets of masses $0.1\sim0.4\,M_{\rm{Jup}}$ would be required. The other three disks with no observed structures are the most compact and faintest in our sample. The emission of $^{12}$CO and $^{13}$CO is more extended than the dust continuum emission in all disks of our sample. When using the $^{12}$CO emission to determine the gas disk extension $R_{\rm{gas}}$, the ratio of $R_{\rm{gas}}/R_{\rm{dust}}$ in our sample varies from 2.3 to 6.0, which is consistent with models of radial drift being very efficient around VLMS in the absence of substructures. Our observations do not exclude giant planet formation on the substructures observed. A comparison of the size and luminosity of VLMS disks with their counterparts around higher mass stars shows that they follow a similar relation.

preprint2020arXiv

The evolution of dust-disk sizes from a homogeneous analysis of 1-10 Myr-old stars

We utilize ALMA archival data to estimate the dust disk size of 152 protoplanetary disks in Lupus (1-3 Myr), Chamaeleon I (2-3 Myr), and Upper-Sco (5-11 Myr). We combine our sample with 47 disks from Tau/Aur and Oph whose dust disk radii were estimated, as here, through fitting radial profile models to visibility data. We use these 199 homogeneously derived disk sizes to identify empirical disk-disk and disk-host property relations as well as to search for evolutionary trends. In agreement with previous studies, we find that dust disk sizes and millimeter luminosities are correlated, but show for the first time that the relationship is not universal between regions. We find that disks in the 2-3 Myr-old Cha I are not smaller than disks in other regions of similar age, and confirm the Barenfeld et al. (2017) finding that the 5-10 Myr USco disks are smaller than disks belonging to younger regions. Finally, we find that the outer edge of the Solar System, as defined by the Kuiper Belt, is consistent with a population of dust disk sizes which have not experienced significant truncation.

preprint2019arXiv

The Exoplanet Population Observation Simulator. II -- Population Synthesis in the Era of Kepler

The collection of planetary system properties derived from large surveys such as Kepler provides critical constraints on planet formation and evolution. These constraints can only be applied to planet formation models, however, if the observational biases and selection effects are properly accounted for. Here we show how epos, the Exoplanet Population Observation Simulator, can be used to constrain planet formation models by comparing the Bern planet population synthesis models to the Kepler exoplanetary systems. We compile a series of diagnostics, based on occurrence rates of different classes of planets and the architectures of multi-planet systems, that can be used as benchmarks for future and current modeling efforts. Overall, we find that a model with 100 seed planetary cores per protoplanetary disk provides a reasonable match to most diagnostics. Based on these diagnostics we identify physical properties and processes that would result in the Bern model more closely matching the known planetary systems. These are: moving the planet trap at the inner disk edge outward; increasing the formation efficiency of mini-Neptunes; and reducing the fraction of stars that form observable planets. We conclude with an outlook on the composition of planets in the habitable zone, and highlight that the majority of simulated planets smaller than 1.7 Earth radii have substantial hydrogen atmospheres. The software used in this paper is available online for public scrutiny at https://github.com/GijsMulders/epos

preprint2016arXiv

A candidate planetary-mass object with a photoevaporating disk in Orion

In this work, we report the discovery of a candidate planetary-mass object with a photoevaporating protoplanetary disk, Proplyd 133-353, which is near the massive star $θ^{1}$ Ori C at the center of the Orion Nebula Cluster (ONC). The object was known to have extended emission pointing away from $θ^{1}$ Ori C, indicating ongoing external photoevaporation. Our near-infrared spectroscopic data suggests that the central source of Proplyd 133-353 is substellar ($\sim$M9.5), might have a mass probably less than 13 Jupiter mass and an age younger than 0.5 Myr. Proplyd 133-353 shows a similar ratio of X-ray luminosity to stellar luminosity to other young stars in the ONC with a similar stellar luminosity, and has a similar proper motion to the mean one of confirmed ONC members. We propose that Proplyd 133-353 was formed in a very low-mass dusty cloud near $θ^{1}$ Ori C as a second-generation of star formation, which can explain both its young age and the presence of its disk.

preprint2016arXiv

M Stars in the TW Hya Association: Stellar X-rays and Disk Dissipation

To investigate the potential connection between the intense X-ray emission from young, low-mass stars and the lifetimes of their circumstellar, planet-forming disks, we have compiled the X-ray luminosities ($L_X$) of M stars in the $\sim$8 Myr-old TW Hya Association (TWA) for which X-ray data are presently available. Our investigation includes analysis of archival Chandra data for the TWA binary systems TWA 8, 9, and 13. Although our study suffers from poor statistics for stars later than M3, we find a trend of decreasing $L_X/L_{bol}$ with decreasing $T_{eff}$ for TWA M stars wherein the earliest-type (M0--M2) stars cluster near $\log{(L_X/L_{bol})} \approx -3.0$ and then $\log{(L_X/L_{bol})}$ decreases, and its distribution broadens, for types M4 and later. The fraction of TWA stars that display evidence for residual primordial disk material also sharply increases in this same (mid-M) spectral type regime. This apparent anticorrelation between the relative X-ray luminosities of low-mass TWA stars and the longevities of their circumstellar disks suggests that primordial disks orbiting early-type M stars in the TWA have dispersed rapidly as a consequence of their persistent large X-ray fluxes. Conversely, the disks orbiting the very lowest-mass pre-MS stars and pre-MS brown dwarfs in the Association may have survived because their X-ray luminosities and, hence, disk photoevaporation rates are very low to begin with, and then further decline relatively early in their pre-MS evolution.

preprint2015arXiv

A stellar-mass-dependent drop in planet occurrence rates

The Kepler Spacecraft has discovered a large number of planets up to one-year periods and down to terrestrial sizes. While the majority of the target stars are main-sequence dwarfs of spectral type F, G, and K, Kepler covers stars with effective temperature as low as 2500 K, which corresponds to M stars. These cooler stars allow characterization of small planets near the habitable zone, yet it is not clear if this population is representative of that around FGK stars. In this paper, we calculate the occurrence of planets around stars of different spectral types as a function of planet radius and distance from the star, and show that they are significantly different from each other. We further identify two trends: First, the occurrence of Earth to Neptune-sized planets is successively higher toward later spectral types at all orbital periods probed by Kepler; Planets around M stars occur twice as frequently as around G stars, and thrice as frequently as around F stars. Second, a drop in planet occurrence is evident at all spectral types inward of a 10 day orbital period, with a plateau further out. By assigning to each spectral type a median stellar mass, we show that the distance from the star where this drop occurs is stellar mass dependent, and scales with semi-major axis as the cube root of stellar mass. By comparing different mechanisms of planet formation, trapping and destruction, we find that this scaling best matches the location of the pre-main-sequence co-rotation radius, indicating efficient trapping of migrating planets or planetary building blocks. These results demonstrate the stellar-mass dependence of the planet population, both in terms of occurrence rate and of orbital distribution. The prominent stellar-mass dependence of the inner boundary of the planet population shows that the formation or migration of planets is sensitive to the stellar parameters.

preprint2015arXiv

An Increase in the Mass of Planetary Systems around Lower-Mass Stars

Trends in the planet population with host star mass provide an avenue to constrain planet formation theories. We derive the planet radius distribution function for Kepler stars of different spectral types, sampling a range in host star masses. We find that M dwarf stars have 3.5 times more small planets (1.0-2.8 R_Earth) than main-sequence FGK stars, but two times fewer Neptune-sized and larger planets (>2.8 R_Earth). We find no systematic trend in the planet size distribution between spectral types F, G, and K to explain the increasing occurrence rates. Taking into account the mass-radius relationship and heavy-element mass of observed exoplanets, and assuming those are independent of spectral type, we derive the inventory of the heavy-element mass locked up in exoplanets at short orbits. The overall higher planet occurrence rates around M stars are not consistent with the redistribution of the same mass into more, smaller planets. At the orbital periods and planet radii where Kepler observations are complete for all spectral types, the average heavy-element mass locked up in exoplanets increases roughly inversely with stellar mass from 4 M_Earth in F stars to 5 M_Earth in G and K stars to 7 M_Earth in M stars. This trend stands in stark contrast with observed protoplanetary disk masses that decrease towards lower mass stars, and provides a challenge for current planet formation models. Neither models of in situ formation nor migration of fully-formed planets are consistent with these results. Instead, these results are indicative of large-scale inward migration of planetary building blocks --- either through type-I migration or radial drift of dust grains --- that is more efficient for lower mass stars, but does not result in significantly larger or smaller planets.

preprint2015arXiv

The snow line in viscous disks around low-mass stars: implications for water delivery to terrestrial planets in the habitable zone

The water ice or snow line is one of the key properties of protoplanetary disks that determines the water content of terrestrial planets in the habitable zone. Its location is determined by the properties of the star, the mass accretion rate through the disk, and the size distribution of dust suspended in the disk. We calculate the snow line location from recent observations of mass accretion rates and as a function of stellar mass. By taking the observed dispersion in mass accretion rates as a measure of the dispersion in initial disk mass, we find that stars of a given mass will exhibit a range of snow line locations. At a given age and stellar mass, the observed dispersion in mass accretion rates of 0.4 dex naturally leads to a dispersion in snow line locations of 0.2 dex. For ISM-like dust sizes, the one-sigma snow line location among solar mass stars of the same age ranges from 2 to 5 au. For more realistic dust opacities that include larger grains, the snow line is located up to two times closer to the star. We use these locations and the outcome of N-body simulations to predict the amount of water delivered to terrestrial planets that formed in situ in the habitable zone. We find that the dispersion in snow line locations leads to a large range in water content. For ISM-like dust sizes, a significant fraction of habitable-zone terrestrial planets around sun-like stars remain dry, and no water is delivered to the habitable zones of low-mass M stars (less than half a solar mass) as in previous works. The closer-in snow line in disks with larger grains enables water delivery to the habitable zone for a significant fraction of M stars and all FGK stars. Considering their larger numbers and higher planet occurrence, M stars may host most of the water-rich terrestrial planets in the galaxy if these planets are able to hold on to their water in their subsequent evolution.

preprint2015arXiv

Volatile Delivery to Planets from Water-rich Planetesimals around Low Mass Stars

Most models of volatile delivery to accreting terrestrial planets assume that the carriers for water are similar in water content to the carbonaceous chondrites in our Solar System. Here we suggest that the water content of primitive bodies in many planetary systems may actually be much higher, as carbonaceous chondrites have lost some of their original water due to heating from short-lived radioisotopes that drove parent body alteration. Using N-body simulations, we explore how planetary accretion would be different if bodies beyond the water line contained a water mass fraction consistent with chemical equilibrium calculations, and more similar to comets, as opposed to the more traditional water-depleted values. We apply this model to consider planet formation around stars of different masses and identify trends in the properties of Habitable Zone planets and planetary system architecture which could be tested by ongoing exoplanet census data collection. Comparison of such data with the model predicted trends will serve to evaluate how well the N-body simulations and the initial conditions used in studies of planetary accretion can be used to understand this stage of planet formation.

preprint2014arXiv

Constraints on photoevaporation models from (lack of) radio emission in the Corona Australis protoplanetary disks

Photoevaporation due to high-energy stellar photons is thought to be one of the main drivers of protoplanetary disk dispersal. The fully or partially ionized disk surface is expected to produce free-free continuum emission at centimeter (cm) wavelengths that can be routinely detected with interferometers such as the upgraded Very Large Array (VLA). We use deep (rms noise down to 8 $μ$Jy beam$^{-1}$ in the field of view center) 3.5 cm maps of the nearby (130 pc) Corona Australis (CrA) star formation (SF) region to constrain disk photoevaporation models. We find that the radio emission from disk sources in CrA is surprisingly faint. Only 3 out of 10 sources within the field of view are detected, with flux densities of order $10^2$ $μ$Jy. However, a significant fraction of their emission is non-thermal. Typical upper limits for non-detections are $3σ\sim 60~μ$Jy beam$^{-1}$. Assuming analytic expressions for the free-free emission from extreme-UV (EUV) irradiation, we derive stringent upper limits to the ionizing photon luminosity impinging on the disk surface $Φ_\mathrm{EUV}<1-4\times10^{41}$ s$^{-1}$. These limits constrain $Φ_\mathrm{EUV}$ to the low end of the values needed by EUV-driven photoevaporation models to clear protoplanetary disks in the observed few Myr timescale. Therefore, at least in CrA, EUV-driven photoevaporation is not likely to be the main agent of disk dispersal. We also compare the observed X-ray luminosities $L_X$ of disk sources with models in which photoevaporation is driven by such photons. Although predictions are less specific than for the EUV case, most of the observed fluxes (upper limits) are roughly consistent with the (scaled) predictions. Deeper observations, as well as predictions spanning a wider parameter space, are needed to properly test X-ray driven photoevaporation.

preprint2014arXiv

LEECH: A 100 Night Exoplanet Imaging Survey at the LBT

In February 2013, the LEECH (LBTI Exozodi Exoplanet Common Hunt) survey began its 100-night campaign from the Large Binocular Telescope atop Mount Graham in Arizona. LEECH nearly complements other high-contrast planet imaging efforts by observing stars in L' band (3.8 microns) as opposed to the shorter wavelength near-infrared bands (1-2.3 microns). This part of the spectrum offers deeper mass sensitivity for intermediate age (several hundred Myr-old) systems, since their Jovian-mass planets radiate predominantly in the mid-infrared. In this proceedings, we present the science goals for LEECH and a preliminary contrast curve from some early data.

preprint2013arXiv

Protoplanetary Disk Masses from Stars to Brown Dwarfs

We present SCUBA-2 850um observations for 7 very low mass stars (VLMS) and brown dwarfs (BDs): 3 in Taurus, 4 in the TWA, and all classical T Tauri (cTT) analogs. We detect 2 of the 3 Taurus disks, but none of the TWA ones. Our 3sigma limits correspond to a dust mass of 1.2 MEarth in Taurus and a mere 0.2 MEarth in the TWA (3--10x deeper than previous work). We combine our data with other sub-mm/mm surveys of Taurus, rho Oph and the TWA to investigate trends in disk mass and grain growth during the cTT phase. We find : (1) The minimum disk outer radius required to explain the upper envelope of sub-mm/mm fluxes is 100 AU for intermediate-mass stars, solar-types and VLMS, and 20 AU for BDs. (2) While the upper envelope of disk masses increases with Mstar from BDs to VLMS to solar-types, no increase is seen from solar-type to intermediate-mass stars. We propose this is due to enhanced photoevaporation around intermediate masses. (3) Many disks around Taurus and rho Oph intermediate-mass and solar-type stars evince an opacity index beta of 0--1, indicating large grains. Of the only four VLMS/BDs in these regions with multi-wavelength data, three are consistent with large grains, though optically thick disks are not ruled out. (4) For the TWA VLMS (TWA 30A, B), combining our fluxes with Mdot and ages suggests substantial grain growth by 10 Myr. The degree of grain growth in the TWA BDs (2M1207A, SSPM1102) remains largely unknown. (5) A Bayesian analysis shows that mean(log[Mdisk/Mstar]) = -2.4, roughly constant all the way from intermediate-mass stars to VLMS/BDs, and (6) the disk mass in close solar-type Taurus binaries is significantly lower than in singles (by a factor of 10), while that in wide solar-type Taurus binaries is closer to that in singles (lower by a factor of 3). (7) We discuss the implications for planet formation, and for the dependence of Mdot on Mstar.

preprint2013arXiv

The Atomic and Molecular Content of Disks Around Very Low-mass Stars and Brown Dwarfs

There is growing observational evidence that disk evolution is stellar-mass dependent. Here, we show that these dependencies extend to the atomic and molecular content of disk atmospheres. We analyze a unique dataset of high-resolution Spitzer/IRS spectra from 8 very low-mass star and brown dwarf disks. We report the first detections of Ne+, H2, CO2, and tentative detections of H2O toward these faint and low-mass disks. Two of our [NeII] 12.81 micron emission lines likely trace the hot (>5,000 K) disk surface irradiated by X-ray photons from the central stellar/sub-stellar object. The H2 S(2) and S(1) fluxes are consistent with arising below the fully or partially ionized surface traced by the [NeII] emission, in gas at about 600 K. We confirm the higher C2H2/HCN flux and column density ratio in brown dwarf disks previously noted from low-resolution IRS spectra. Our high-resolution spectra also show that the HCN/H2O fluxes of brown dwarf disks are on average higher than those of T Tauri disks. Our LTE modeling hints that this difference extends to column density ratios if H2O lines trace warm > 600 K disk gas. These trends suggest that the inner regions of brown dwarf disks have a lower O/C ratio than those of T Tauri disks which may result from a more efficient formation of non-migrating icy planetesimals. A O/C=1, as inferred from our analysis, would have profound implications on the bulk composition of rocky planets that can form around very low-mass stars and brown dwarfs.

preprint2013arXiv

The Dispersal of Protoplanetary Disks

Protoplanetary disks are the sites of planet formation, and the evolution and eventual dispersal of these disks strongly influences the formation of planetary systems. Disk evolution during the planet-forming epoch is driven by accretion and mass-loss due to winds, and in typical environments photoevaporation by high-energy radiation from the central star is likely to dominate final gas disk dispersal. We present a critical review of current theoretical models, and discuss the observations that are used to test these models and inform our understanding of the underlying physics. We also discuss the role disk dispersal plays in shaping planetary systems, considering its influence on both the process(es) of planet formation and the architectures of planetary systems. We conclude by presenting a schematic picture of protoplanetary disk evolution and dispersal, and discussing prospects for future work.

preprint2013arXiv

The VLTI/MIDI survey of massive young stellar objects - Sounding the inner regions around intermediate- and high-mass young stars using mid-infrared interferometry

We aim to characterize the distribution and composition of circumstellar material around young massive stars, and to investigate exactly which physical structures in these objects are probed by long-baseline mid-infrared interferometric observations. We used the two-telescope interferometric instrument MIDI of the Very Large Telescope Interferometer of the European Southern Observatory to observe a sample of 24 intermediate- and high-mass young stellar objects in the N band (8-13 micron). We had successful fringe detections for 20 objects, and present spectrally-resolved correlated fluxes and visibility levels for projected baselines of up to 128 m. We fit the visibilities with geometric models to derive the sizes of the emitting regions, as well as the orientation and elongation of the circumstellar material. Fourteen objects in the sample show the 10 micron silicate feature in absorption in the total and correlated flux spectra. For 13 of these objects, we were able to fit the correlated flux spectra with a simple absorption model, allowing us to constrain the composition and absorptive properties of the circumstellar material. Nearly all of the massive young stellar objects observed show significant deviations from spherical symmetry at mid-infrared wavelengths. In general, the mid-infrared emission can trace both disks and outflows, and in many cases it may be difficult to disentangle these components on the basis of interferometric data alone, because of the sparse spatial frequency coverage normally provided by current long-baseline interferometers. For the majority of the objects in this sample, the absorption occurs on spatial scales larger than those probed by MIDI. Finally, the physical extent of the mid-infrared emission around these sources is correlated with the total luminosity, albeit with significant scatter.

preprint2013arXiv

Understanding the origin of the [OI] low-velocity component from T Tauri stars

The formation time, masses, and location of planets are strongly impacted by the physical mechanisms that disperse protoplanetary disks and the timescale over which protoplanetary material is cleared out. Accretion of matter onto the central star, protostellar winds/jets, magnetic disk winds, and photoevaporative winds operate concurrently. Hence, disentangling their relative contribution to disk dispersal requires identifying diagnostics that trace different star-disk environments. Here, we analyze the low velocity component (LVC) of the Oxygen optical forbidden lines, which is found to be blueshifted by a few km/s with respect to the stellar velocity. We find that the [OI] LVC profiles are different from those of [NeII] at 12.81μm and CO at 4.7μm lines pointing to different origins for these gas lines. We report a correlation between the luminosity of the [OI] LVC and the accretion luminosity Lacc. We do not find any correlation with the X-ray luminosity, while we find that the higher is the stellar FUV luminosity, the higher is the luminosity of the [OI] LVC. In addition, we show that the [OI]6300Å/5577Å ratio is low (ranging between 1 and 8). These findings favor an origin of the [OI] LVC in a region where OH is photodissociated by stellar FUV photons and argue against thermal emission from an X-ray-heated layer. Detailed modeling of two spectra with the highest S/N and resolution shows that there are two components within the LVC: a broad, centrally peaked component that can be attributed to gas arising in a warm disk surface in Keplerian rotation (with FWHM between ~40 and ~60 km/s), and a narrow component (with FWHM ~10 km/s and small blueshifts of ~2 km/s) that may arise in a cool (<1,000 K) molecular wind.

preprint2012arXiv

A Herschel Survey of Cold Dust in Disks Around Brown Dwarfs and Low-Mass Stars

We report the complete photometric results from our Herschel study which is the first comprehensive program to search for far-infrared emission from cold dust around young brown dwarfs. We surveyed 50 fields containing 51 known or suspected brown dwarfs and very low mass stars that have evidence of circumstellar disks based on Spitzer photometry and/or spectroscopy. The objects with known spectral types range from M3 to M9.5. Four of the candidates were subsequently identified as extragalactic objects. Of the remaining 47 we have successfully detected 36 at 70micron and 14 at 160micron with S/N greater than 3, as well as several additional possible detections with low S/N. The objects exhibit a range of [24]--[70] micron colors suggesting a range in mass and/or structure of the outer disk. We present modeling of the spectral energy distributions of the sample and discuss trends visible in the data. Using two Monte Carlo radiative transfer codes we investigate disk masses and geometry. We find a very wide range in modeled total disk masses from less than 1e-6 solar masses up to 1e-3 solar masses with a median disk mass of order 3e-5 solar masses, suggesting that the median ratio of disk mass to central object mass may be lower than for T Tauri stars. The disk scale heights and flaring angles, however, cover a range consistent with those seen around T Tauri stars. The host clouds in which the young brown dwarfs and low-mass stars are located span a range in estimated age from ~1-3 Myr to ~10 Myr and represent a variety of star-forming environments. No obvious dependence on cloud location or age is seen in the disk properties, though the statistical significance of this conclusion is not strong.

preprint2012arXiv

Observational Constraints on the Stellar Radiation Field Impinging on Transitional Disk Atmospheres

Mid-infrared atomic and ionic line ratios measured in spectra of pre-main sequence stars are sensitive indicators of the hardness of the radiation field impinging on the disk surface. We present a low-resolution Spitzer IRS search for [Ar II] at 6.98 $μ$m, [Ne II] at 12.81 $μ$m, and [Ne III] 15.55 $μ$m lines in 56 transitional disks. These objects, characterized by reduced near-infrared but strong far-infrared excess emission, are ideal targets to set constraints on the stellar radiation field onto the disk because their spectra are not contaminated by shock emission from jets/outflows or by molecular emission lines. After demonstrating that we can detect [Ne II] lines and recover their fluxes from the low-resolution spectra, here we report the first detections of [Ar II] lines towards protoplanetary disks. We did not detect [Ne III] emission in any of our sources. Our [Ne II]/[Ne III] line flux ratios combined with literature data suggest that a soft-EUV or X-ray spectrum produces these gas lines. Furthermore, the [Ar II]/[Ne II] line flux ratios point to a soft X-ray and/or soft-EUV stellar spectrum as the ionization source of the [Ar II] and [Ne II] emitting layer of the disk. If the soft X-ray component dominates over the EUV than we would expect larger photoevaporation rates hence a reduction of the time available to form planets.

preprint2011arXiv

A Herschel Search For Cold Dust in Brown Dwarf Disks: First Results

We report initial results from a {\it Herschel} program to search for far-infrared emission from cold dust around a statistically significant sample of young brown dwarfs. The first three objects in our survey are all detected at 70\micron, and we report the first detection of a brown dwarf at 160\micron. The flux densities are consistent with the presence of substantial amounts of cold dust in the outer disks around these objects. We modeled the SED's with two different radiative transfer codes. We find that a broad range of model parameters provides a reasonable fit to the SED's, but that the addition of our 70\micron, and especially the 160\micron\ detection enables strong lower limits to be placed on the disk masses since most of the mass is in the outer disk. We find likely disk masses in the range of a few $\times 10^{-6}$ to $10^{-4}$ \msun. Our models provide a good fit to the SED's and do not require dust settling.

preprint2011arXiv

Disk Imaging Survey of Chemistry with SMA: II. Southern Sky Protoplanetary Disk Data and Full Sample Statistics

This is the second in a series of papers based on data from DISCS, a Submillimeter Array observing program aimed at spatially and spectrally resolving the chemical composition of 12 protoplanetary disks. We present data on six Southern sky sources - IM Lup, SAO 206462 (HD 135344b), HD 142527, AS 209, AS 205 and V4046 Sgr - which complement the six sources in the Taurus star forming region reported previously. CO 2-1 and HCO+ 3-2 emission are detected and resolved in all disks and show velocity patterns consistent with Keplerian rotation. Where detected, the emission from DCO+ 3-2, N2H+ 3-2, H2CO 3-2 and 4-3,HCN 3-2 and CN 2-1 are also generally spatially resolved. The detection rates are highest toward the M and K stars, while the F star SAO 206462 has only weak CN and HCN emission, and H2CO alone is detected toward HD 142527. These findings together with the statistics from the previous Taurus disks, support the hypothesis that high detection rates of many small molecules depend on the presence of a cold and protected disk midplane, which is less common around F and A stars compared to M and K stars. Disk-averaged variations in the proposed radiation tracer CN/HCN are found to be small, despite two orders of magnitude range of spectral types and accretion rates. In contrast, the resolved images suggest that the CN/HCN emission ratio varies with disk radius in at least two of the systems. There are no clear observational differences in the disk chemistry between the classical/full T Tauri disks and transitional disks. Furthermore, the observed line emission does not depend on measured accretion luminosities or the number of infrared lines detected, which suggests that the chemistry outside of 100 AU is not coupled to the physical processes that drive the chemistry in the innermost few AU.

preprint2011arXiv

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 sub-millimeter 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 $\sim10^{-4}-10^{-5}$\ms\ for $0.06{\rm AU} <r<3.5$AU and $\sim 0.06$\ms\ for $ 3.5 {\rm AU} <r<200$AU. We find that the inner dust hole ($r<3.5$AU) in the disk must be depleted of gas by $\sim 1-2$ orders of magnitude compared to the extrapolated surface density distribution of the outer disk. Grain growth alone is therefore not a viable explanation for the dust hole. CO vibrational emission arises within $r\sim 0.5$AU from thermal excitation of gas. [OI] 6300Å and 5577Å forbidden lines and OH mid-infrared emission are mainly due to prompt emission following UV photodissociation of OH and water at $r\lesssim0.1$AU and at $r\sim 4$AU. [NeII] emission is consistent with an origin in X-ray heated neutral gas at $r\lesssim 10$AU, and may not require the presence of a significant EUV ($hν>13.6$eV) flux from TW Hya. H$_2$ pure rotational line emission comes primarily from $r\sim 1-30$AU. [OI]63$μ$m, HCO$^+$ and CO pure rotational lines all arise from the outer disk at $r\sim30-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 $\sim 4-7$M$_J$ situated at $\sim 3$AU. Using our photoevaporation models and the best surface density profile match to observations, we estimate a current photoevaporative mass loss rate of $4\times10^{-9}$\ms\ yr$^{-1}$ and a remaining disk lifetime of $\sim 5$ million years.

preprint2011arXiv

Evidence Against an Edge-On Disk Around the Extrasolar Planet 2MASS 1207 b and a New Thick Cloud Explanation for its Under-Luminosity

(Abridged) Since the discovery of the first directly-imaged, planetary-mass object, 2MASS 1207 b, several works have sought to explain a disparity between its observed temperature and luminosity. Given its known age, distance, and spectral type, 2MASS 1207 b is under-luminous by a factor of ~10 (~2.5 mags) when compared to standard models of brown-dwarf/giant-planet evolution. In this paper, we study three possible sources of 2MASS 1207 b's under-luminosity. First, we investigate Mohanty et al. (2007)'s hypothesis that a near edge-on disk might be responsible for 2MASS 1207 b's under-luminosity. We conclude that the hypothesis is unlikely due to the lack of variability seen in multi-epoch photometry and unnecessary due to the increasing sample of under-luminous brown-dwarfs/giant-exoplanets that cannot be explained by an edge-on disk. Next, we test the analogous possibility that a spherical shell of dust, could explain 2MASS 1207 b's under-luminosity. Models containing enough dust to create ~2.5 mags of extinction, placed at reasonable radii, are ruled out by our new Gemini/T-ReCS 8.7 micron photometric upper-limit for 2MASS 1207 b. Finally, we investigate the possibility that 2MASS 1207 b is intrinsically cooler than the commonly used AMES-DUSTY fits to its spectrum, and thus it is not, in fact, under-luminous. New, thick cloud model grids by Madhusudhan et al. (2011) fit 2MASS 1207 b's 1-10 micron SED well, but they do not quite fit its near-infrared spectrum. However, we suggest that with some "tuning", they might be capable of simultaneously reproducing 2MASS 1207 b's spectral shape and luminosity. In this case, the whole class of young, under-luminous brown-dwarfs/giant-exoplanets might be explained by atmospheres that are able to suspend thick, dusty clouds in their photospheres at lower temperatures than field brown-dwarfs.

preprint2011arXiv

Measuring Organic Molecular Emission in Disks with Low Resolution Spitzer Spectroscopy

We explore the extent to which Spitzer IRS spectra taken at low spectral resolution can be used in quantitative studies of organic molecular emission from disks surrounding low mass young stars. We use Spitzer IRS spectra taken in both the high and low resolution modules for the same sources to investigate whether it is possible to define line indices that can measure trends in the strength of the molecular features in low resolution data. We find that trends in HCN emission strength seen in the high resolution data can be recovered in low resolution data. In examining the factors that influence the HCN emission strength, we find that the low resolution HCN flux is modestly correlated with stellar accretion rate and X-ray luminosity. Correlations of this kind are perhaps expected based on recent observational and theoretical studies of inner disk atmospheres. Our results demonstrate the potential of using the large number of low resolution disk spectra that reside in the Spitzer archive to study the factors that influence the strength of molecular emission from disks. Such studies would complement results for the much smaller number of circumstellar disks that have been observed at high resolution with IRS.

preprint2010arXiv

Spitzer Spectroscopy of the Transition Object TW Hya

We report sensitive Spitzer IRS spectroscopy in the 10-20 micron region of TW Hya, a nearby T Tauri star. The unusual spectral energy distribution of the source, that of a transition object, indicates that the circumstellar disk in the system has experienced significant evolution, possibly as a result of planet formation. The spectrum we measure is strikingly different from that of other classical T Tauri stars reported in the literature, displaying no strong emission features of water, acetylene, or HCN. The difference indicates that the inner planet formation region (within 5 AU) of the gaseous disk has evolved physically and/or chemically away from the classical T Tauri norm. Nevertheless, TW Hya does show a rich spectrum of emission features of atoms (HI, [NeII], and [NeIII]) and molecules (H2, OH, CO2, HCO+, and possibly CH3), some of which are also detected in classical T Tauri spectra. The properties of the neon emission are consistent with an origin for the emission in a disk irradiated by X-rays (with a possible role for additional irradiation by stellar EUV). The OH emission we detect, which also likely originates in the disk, is hot, arising from energy levels up to 23,000 K above ground, and may be produced by the UV photodissociation of water. The HI emission is surprisingly strong, with relative strengths that are consistent with case B recombination. While the absence of strong molecular emission in the 10-20 micron region may indicate that the inner region of the gaseous disk has been partly cleared by an orbiting giant planet, chemical and/or excitation effects may be responsible instead. We discuss these issues and how our results bear on our understanding of the evolutionary state of the TW Hya disk.

preprint2010arXiv

Stellar-Mass-Dependent Disk Structure in Coeval Planet-Forming Disks

Previous studies suggest that the planet-forming disks around very-low-mass stars/brown dwarfs may be flatter than those around more massive stars, in contrast to model predictions of larger scale heights for gas-disks around lower-mass stars. We conducted a statistically robust study to determine whether there is evidence for stellar-mass-dependent disk structure in planet-forming disks. We find a statistically significant difference in the Spitzer/IRAC color distributions of disks around very-low-mass and low-mass stars all belonging to the same star-forming region, the Chamaeleon I star-forming region. We show that self consistently calculated disk models cannot fit the median spectral energy distributions (SEDs) of the two groups. These SEDs can be only explained by flatter disk models, consistent with the effect of dust settling in disks. We find that relative to the disk structure predicted for flared disks the required reduction in disk scale height is anti-correlated with the stellar mass, i.e. disks around lower-mass stars are flatter. Our results show that the initial and boundary conditions of planet formation are stellar-mass-dependent, an important finding that must be considered in planet formation models.

preprint2009arXiv

Anachronistic Grain Growth and Global Structure of the Protoplanetary Disk Associated with the Mature Classical T Tauri Star, PDS 66

We present ATCA interferometric observations of the old (13 Myr), nearby (86pc) classical T Tauri star, PDS 66. Unresolved 3 and 12 mm continuum emission is detected towards PDS 66, and upper limits are derived for the 3 and 6 cm flux densities. The mm-wave data show a spectral slope flatter than that expected for ISM-sized dust particles, which is evidence of grain growth. We also present HST/NICMOS 1.1 micron PSF-subtracted coronagraphic imaging of PDS 66. The HST observations reveal a bilaterally symmetric circumstellar region of dust scattering about 0.32% of the central starlight, declining radially in surface brightness. The light-scattering disk of material is inclined 32 degrees from face-on, and extends to a radius of 170 AU. These data are combined with published optical and longer wavelength observations to make qualitative comparisons between the median Taurus and PDS 66 spectral energy distributions (SEDs). By comparing the near-infrared emission to a simple model, we determine that the location of the inner disk radius is consistent with the dust sublimation radius (1400 K at 0.1 AU). We place constraints on the total disk mass using a flat-disk model and find that it is probably too low to form gas giant planets according to current models. Despite the fact that PDS 66 is much older than a typical classical T Tauri star (< 5 Myr), its physical properties are not much different.

preprint2009arXiv

Do we really know the dust? Systematics and uncertainties of the mid-infrared spectral analysis methods

The spectral region around 10 micrometer, showing prominent dust emission bands, is commonly used to derive the chemical composition of protoplanetary dust. Different analysis methods have been proposed for this purpose, but so far, no comparative test has been performed to test their validity. We calculated model spectra of disk models with different geometries and central sources, using a 2D radiative transfer code. These spectra were then fitted in a blind test using four different spectral decomposition methods. We studied the effect of disk structure (flared vs. flat), inclination angle, size of the inner disk hole and stellar luminosity on the fitted chemical composition. Our results show that the derived dust compositions by all methods deviate systematically from the real chemical composition. Out of the four tested spectral decomposition methods, our new two-layer temperature distribution method, differs the least from the input dust composition and the results show the weakest systematic effects. The reason for the deviations of the results given by other methods lies in their simplifying assumptions (e.g. single average grain temperature or one component continuum). We also tested the influence of different noise levels on the results of the spectral decomposition methods. We find that, for ground-based observations (8-13 micrometer) the expected uncertainty in the value of the crystallinity is about 11% for a signal-to-noise ratio of 100, while for space-based observations (7-17 micrometer) the uncertainty is about 5%. On the basis of our results, we propose a recipe for the analysis and interpretation of dust spectroscopy data in the mid-infrared which should be especially valuable for analysing Spitzer data and ground-based infrared spectroscopy data in the 10 micrometer window.

preprint2009arXiv

Evidence for Disk Photoevaporation Driven by the Central Star

The lifetime of isolated protoplanetary disks is thought to be set by the combination of viscous accretion and photoevaporation driven by stellar high-energy photons. Observational evidence for magnetospheric accretion in young sun-like stars is robust. Here we report the first observational evidence for disk photoevaporation driven by the central star. We acquired high-resolution (R~30,000) spectra of the [NeII] 12.81 micron line from 7 circumstellar disks using VISIR on Melipal/VLT. We show that the 3 transition disks in the sample all have [NeII] line profiles consistent with those predicted by a photoevaporative flow driven by stellar extreme UV photons. The ~6 km/s blue-shift of the line from the almost face-on disk of TW Hya is clearly inconsistent with emission from a static disk atmosphere and convincingly points to the presence of a photoevaporative wind. We do not detect any [NeII] line close to the stellar velocity from the sample of classical optically thick (non-transition) disks. We conclude that most of the spectrally unresolved [NeII] emission in these less evolved systems arises from jets/outflows rather than from the disk. The pattern of the [NeII] detections and non-detections suggests that extreme UV-driven photoevaporation starts only at a later stage in the disk evolution.

preprint2009arXiv

Far-Ultraviolet H2 Emission from Circumstellar Disks

We analyze the far-ultraviolet (FUV) spectra of 33 classical T Tauri stars (CTTS), including 20 new spectra obtained with the Advanced Camera for Surveys Solar Blind Channel (ACS/SBC) on the Hubble Space Telescope. Of the sources, 28 are in the ~1 Myr old Taurus-Auriga complex or Orion Molecular Cloud, 4 in the 8-10 Myr old Orion OB1a complex and one, TW Hya, in the 10 Myr old TW Hydrae Association. We also obtained FUV ACS/SBC spectra of 10 non-accreting sources surrounded by debris disks with ages between 10 and 125 Myr. We use a feature in the FUV spectra due mostly to electron impact excitation of \h2 to study the evolution of the gas in the inner disk. We find that the \h2 feature is absent in non-accreting sources, but is detected in the spectra of CTTS and correlates with accretion luminosity. Since all young stars have active chromospheres which produce strong X-ray and UV emission capable of exciting \h2 in the disk, the fact that the non-accreting sources show no \h2 emission implies that the \h2 gas in the inner disk has dissipated in the non-accreting sources, although dust (and possibly gas) remains at larger radii. Using the flux at 1600 Å, we estimate that the column density of \h2 left in the inner regions of the debris disks in our sample is less than ~ 3x10^-6 g cm^-2, nine orders of magnitude below the surface density of the minimum mass solar nebula at 1 AU.