Complex organic molecules in protoplanetary disks

Context. Protoplanetary disks are vital objects in star and planet formation, possessing all the material, gas and dust, which may form a planetary system orbiting the new star. Small, simple molecules have traditionally been detected in protoplanetary disks; however, in the ALMA era, we expect the molecular inventory of protoplanetary disks to significantly increase.

Aims. We investigate the synthesis of complex organic molecules (COMs) in protoplanetary disks to put constraints on the achievable chemical complexity and to predict species and transitions which may be observable with ALMA.

Methods. We have coupled a 2D steady-state physical model of a protoplanetary disk around a typical T Tauri star with a large gas-grain chemical network including COMs. We compare the resulting column densities with those derived from observations and perform ray-tracing calculations to predict line spectra. We compare the synthesised line intensities with current observations and determine those COMs which may be observable in nearby objects. We also compare the predicted grain-surface abundances with those derived from cometary comae observations.

Results. We find COMs are efficiently formed in the disk midplane via grain-surface chemical reactions, reaching peak grain-surface fractional abundances similar to 10(-6)-10(-4) that of the H nuclei number density. COMs formed on grain surfaces are returned to the gas phase via non-thermal desorption; however, gas-phase species reach lower fractional abundances than their grain-surface equivalents, similar to 10(-12)-10(-7). Including the irradiation of grain mantle material helps build further complexity in the ice through the replenishment of grain-surface radicals which take part in further grain-surface reactions. There is reasonable agreement with several line transitions of H2CO observed towards T Tauri star-disk systems. There is poor agreement with HC3(N) lines observed towards LkCa 15 and GO Tau and we discuss possible explanations for these discrepancies. The synthesised line intensities for CH3OH are consistent with upper limits determined towards all sources. Our models suggest CH3OH should be readily observable in nearby protoplanetary disks with ALMA; however, detection of more complex species may prove challenging, even with ALMA "Full Science" capabilities. Our grain-surface abundances are consistent with those derived from cometary comae observations providing additional evidence for the hypothesis that comets (and other planetesimals) formed via the coagulation of icy grains in the Sun's natal disk.

Grain growth in protoplanetary disks

The majority of young, low-mass stars are surrounded by optically thick accretion disks. These circumstellar disks provide large reservoirs of gas and dust that will eventually be transformed into planetary systems. Theory and observations suggest that the earliest stage toward planet formation in a protoplanetary disk is the growth of particles, from sub-micron-sized grains to centimeter- sized pebbles. Theory indicates that small interstellar grains are well coupled into the gas and are incorporated to the disk during the proto-stellar collapse. These dust particles settle toward the disk mid-plane and simultaneously grow through collisional coagulation in a very short timescale. Observationally, grain growth can be inferred by measuring the spectral energy distribution at long wavelengths, which traces the continuum dust emission spectrum and hence the dust opacity. Several observational studies have indicated that the dust component in protoplanetary disks has evolved as compared to interstellar medium dust particles, suggesting at least 4 orders of magnitude in particle- size growth. However, the limited angular resolution and poor sensitivity of previous observations has not allowed for further exploration of this astrophysical process.

As part of my thesis, I embarked in an observational program to search for evidence of radial variations in the dust properties across a protoplanetary disk, which may be indicative of grain growth. By making use of high angular resolution observations obtained with CARMA, VLA, and SMA, I searched for radial variations in the dust opacity inside protoplanetary disks. These observations span more than an order of magnitude in wavelength (from sub-millimeter to centimeter wavelengths) and attain spatial resolutions down to 20 AU. I characterized the radial distribution of the circumstellar material and constrained radial variations of the dust opacity spectral index, which may originate from particle growth in these circumstellar disks. Furthermore, I compared these observational constraints with simple physical models of grain evolution that include collisional coagulation, fragmentation, and the interaction of these grains with the gaseous disk (the radial drift problem). For the parameters explored, these observational constraints are in agreement with a population of grains limited in size by radial drift. Finally, I also discuss future endeavors with forthcoming ALMA observations.

YSOVAR: Mid-infrared Variability in the Star-forming Region Lynds 1688

The emission from young stellar objects (YSOs) in the mid-infrared (mid-IR) is dominated by the inner rim of their circumstellar disks. We present IR data from the Young Stellar Object VARiability (YSOVAR) survey of ~800 objects in the direction of the Lynds 1688 (L1688) star-forming region over four visibility windows spanning 1.6 yr using the Spitzer Space Telescope in its warm mission phase. Among all light curves, 57 sources are cluster members identified based on their spectral energy distribution and X-ray emission. Almost all cluster members show significant variability. The amplitude of the variability is larger in more embedded YSOs. Ten out of 57 cluster members have periodic variations in the light curves with periods typically between three and seven days, but even for those sources, significant variability in addition to the periodic signal can be seen. No period is stable over 1.6 yr. Nonperiodic light curves often still show a preferred timescale of variability that is longer for more embedded sources. About half of all sources exhibit redder colors in a fainter state. This is compatible with time-variable absorption toward the YSO. The other half becomes bluer when fainter. These colors can only be explained with significant changes in the structure of the inner disk. No relation between mid-IR variability and stellar effective temperature or X-ray spectrum is found.

Towards the main sequence: detailed analysis of weak line and post-T Tauri stars

A detailed study was performed for a sample of low-mass pre-main-sequence (PMS) stars, previously identified as weak-line T Tauri stars, which are compared to members of the Tucanae and Horologium Associations. Aiming to verify if there is any pattern of abundances when comparing the young stars at different phases, we selected objects in the range from 1 to 100 Myr, which covers most of PMS evolution. High-resolution optical spectra were acquired at European Southern Observatory and Observatorio do Pico dos Dias. The stellar fundamental parameters effective temperature and gravity were calculated by excitation and ionization equilibria of iron absorption lines. Chemical abundances were obtained via equivalent width calculations and spectral synthesis for 44 per cent of the sample, which shows metallicities within 0.5 dex solar. A classification was developed based on equivalent width of Li I 6708 angstrom and Ha lines and spectral types of the studied stars. This classification allowed a separation of the sample into categories that correspond to different evolutive stages in the PMS. The position of these stars in the Hertzsprung-Russell diagram was also inspected in order to estimate their ages and masses. Among the studied objects, it was verified that our sample actually contains seven weak-line T Tauri stars, three are Classical T Tauri, 12 are Fe/Ge PMS stars and 21 are post-T Tauri or young main-sequence stars. An estimation of circumstellar luminosity was obtained using a disc model to reproduce the observed spectral energy distribution. Most of the stars show low levels of circumstellar emission, corresponding to less than 30 per cent of the total emission.

DYNAMICAL EVOLUTION OF VISCOUS DISKS AROUND Be STARS. I. PHOTOMETRY

Be stars possess gaseous circumstellar disks that modify in many ways the spectrum of the central B star. Furthermore, they exhibit variability at several timescales and for a large number of observables. Putting the pieces together of this dynamical behavior is not an easy task and requires a detailed understanding of the physical processes that control the temporal evolution of the observables. There is an increasing body of evidence that suggests that Be disks are well described by standard alpha-disk theory. This paper is the first of a series that aims at studying the possibility of inferring several disk and stellar parameters through the follow-up of various observables. Here we study the temporal evolution of the disk density for different dynamical scenarios, including the disk buildup as a result of a long and steady mass injection from the star, the disk dissipation that occurs after mass injection is turned off, as well as scenarios in which active periods are followed by periods of quiescence. For those scenarios, we investigate the temporal evolution of continuum photometric observables using a three-dimensional non-LTE radiative transfer code. We show that light curves for different wavelengths are specific of a mass loss history, inclination angle, and alpha viscosity parameter. The diagnostic potential of those light curves is also discussed.

Détection et caractérisation de nouveaux disques circumstellaires autour d’étoiles de faibles masses et naines brunes jeunes.

La présence de disques circumstellaires signale la formation actuelle ou passée de systèmes planétaires, pour lesquels les processus de formation sont encore mal compris. Ce mémoire porte sur la détection et la caractérisation de disques circumstellaires autour d’étoiles de faibles masses (types spectraux > K5) et de naines brunes qui sont candidates ou membres d’associations cinématiques jeunes. Nous présentons ici les résultats de cette recherche ainsi que son implication pour la compréhension des processus de formation et d’évolution des systèmes planétaires. De l’échantillon initial composé de ∼ 1600 objets provenant des relevés BANYAN de Malo et al. ainsi que Gagné et al., dont seulement 600 satisfont nos critères de qualité sur les données, quatre nouveaux candidats de disque ont été découverts en détectant leur excès d’émission infrarouge dans les données d’archive de la mission WISE. Les données du relevé 2MASS ainsi que les spectres synthétiques BT-Settl ont été conjointement utilisés pour modéliser l’émission des étoiles. Les nouveaux candidats, dont les types spectraux sont tardifs (M4.5 à L0) et les masses se situent entre ∼ 13 et 120 M_Jup, ont des températures de disque de ∼ 135–520 K et des luminosités fractionnaires de 0,021–0,15. Pour deux des cibles, nous avons obtenu des spectres dans les longueurs d’onde visibles et infrarouges proches. Ces nouveaux spectres montrent respectivement des signes d’émission en Hα et Paβ, indiquant la présence d’accrétion, et ainsi de gaz, et renforçant l’hypothèse que ces objets sont réellement jeunes. Ces deux objets, vraisemblablement âgés de 40 Ma, pourraient représenter la première détection et caractérisation de disques porteurs de gaz plus vieux que 20 Ma.

Détection et caractérisation de nouveaux disques circumstellaires autour d’étoiles de faibles masses et naines brunes jeunes

La présence de disques circumstellaires signale la formation actuelle ou passée de systèmes planétaires, pour lesquels les processus de formation sont encore mal compris. Ce mémoire porte sur la détection et la caractérisation de disques circumstellaires autour d’étoiles de faibles masses (types spectraux > K5) et de naines brunes qui sont candidates ou membres d’associations cinématiques jeunes. Nous présentons ici les résultats de cette recherche ainsi que son implication pour la compréhension des processus de formation et d’évolution des systèmes planétaires. De l’échantillon initial composé de ∼ 1600 objets provenant des relevés BANYAN de Malo et al. ainsi que Gagné et al., dont seulement 600 satisfont nos critères de qualité sur les données, quatre nouveaux candidats de disque ont été découverts en détectant leur excès d’émission infrarouge dans les données d’archive de la mission WISE. Les données du relevé 2MASS ainsi que les spectres synthétiques BT-Settl ont été conjointement utilisés pour modéliser l’émission des étoiles. Les nouveaux candidats, dont les types spectraux sont tardifs (M4.5 à L0) et les masses se situent entre ∼ 13 et 120 M_Jup, ont des températures de disque de ∼ 135–520 K et des luminosités fractionnaires de 0,021–0,15. Pour deux des cibles, nous avons obtenu des spectres dans les longueurs d’onde visibles et infrarouges proches. Ces nouveaux spectres montrent respectivement des signes d’émission en Hα et Paβ, indiquant la présence d’accrétion, et ainsi de gaz, et renforçant l’hypothèse que ces objets sont réellement jeunes. Ces deux objets, vraisemblablement âgés de 40 Ma, pourraient représenter la première détection et caractérisation de disques porteurs de gaz plus vieux que 20 Ma.

CSI 2264: Simultaneous Optical and Infrared Light Curves of Young Disk-bearing Stars in NGC 2264 with CoRoT and Spitzer—Evidence for Multiple Origins of Variability

We present the Coordinated Synoptic Investigation of NGC 2264, a continuous 30 day multi-wavelength photometric monitoring campaign on more than 1000 young cluster members using 16 telescopes. The unprecedented combination of multi-wavelength, high-precision, high-cadence, and long-duration data opens a new window into the time domain behavior of young stellar objects. Here we provide an overview of the observations, focusing on results from Spitzer and CoRoT. The highlight of this work is detailed analysis of 162 classical T Tauri stars for which we can probe optical and mid-infrared flux variations to 1% amplitudes and sub-hour timescales. We present a morphological variability census and then use metrics of periodicity, stochasticity, and symmetry to statistically separate the light curves into seven distinct classes, which we suggest represent different physical processes and geometric effects. We provide distributions of the characteristic timescales and amplitudes and assess the fractional representation within each class. The largest category (>20%) are optical "dippers" with discrete fading events lasting ~1-5 days. The degree of correlation between the optical and infrared light curves is positive but weak; notably, the independently assigned optical and infrared morphology classes tend to be different for the same object. Assessment of flux variation behavior with respect to (circum)stellar properties reveals correlations of variability parameters with Hα emission and with effective temperature. Overall, our results point to multiple origins of young star variability, including circumstellar obscuration events, hot spots on the star and/or disk, accretion bursts, and rapid structural changes in the inner disk. Based on data from the Spitzer and CoRoT missions. The CoRoT space mission was developed and is operated by the French space agency CNES, with participation of ESA's RSSD and Science Programmes, Austria, Belgium, Brazil, Germany, and Spain.

YSOVAR: Mid-infrared Variability of Young Stellar Objects and Their Disks in the Cluster IRAS 20050+2720

We present a time-variability study of young stellar objects (YSOs) in the cluster IRAS 20050+2720, performed at 3.6 and 4.5 μm with the Spitzer Space Telescope; this study is part of the Young Stellar Object VARiability (YSOVAR) project. We have collected light curves for 181 cluster members over 60 days. We find a high variability fraction among embedded cluster members of ca. 70%, whereas young stars without a detectable disk display variability less often (in ca. 50% of the cases) and with lower amplitudes. We detect periodic variability for 33 sources with periods primarily in the range of 2–6 days. Practically all embedded periodic sources display additional variability on top of their periodicity. Furthermore, we analyze the slopes of the tracks that our sources span in the color–magnitude diagram (CMD). We find that sources with long variability time scales tend to display CMD slopes that are at least partially influenced by accretion processes, while sources with short variability timescales tend to display extinction-dominated slopes. We find a tentative trend of X-ray detected cluster members to vary on longer timescales than the X-ray undetected members.

New Insights from Aperiodic Variability of Young Stars

Nearly all young stars are variable, with the variability traditionally divided into two classes: periodic variables and aperiodic or "irregular" variables. Periodic variables have been studied extensively, typically using periodograms, while aperiodic variables have received much less attention due to a lack of standard statistical tools. However, aperiodic variability can serve as a powerful probe of young star accretion physics and inner circumstellar disk structure. For my dissertation, I analyzed data from a large-scale, long-term survey of the nearby North America Nebula complex, using Palomar Transient Factory photometric time series collected on a nightly or every few night cadence over several years. This survey is the most thorough exploration of variability in a sample of thousands of young stars over time baselines of days to years, revealing a rich array of lightcurve shapes, amplitudes, and timescales.

I have constrained the timescale distribution of all young variables, periodic and aperiodic, on timescales from less than a day to ~100 days. I have shown that the distribution of timescales for aperiodic variables peaks at a few days, with relatively few (~15%) sources dominated by variability on tens of days or longer. My constraints on aperiodic timescale distributions are based on two new tools, magnitude- vs. time-difference (Δm-Δt) plots and peak-finding plots, for describing aperiodic lightcurves; this thesis provides simulations of their performance and presents recommendations on how to apply them to aperiodic signals in other time series data sets. In addition, I have measured the error introduced into colors or SEDs from combining photometry of variable sources taken at different epochs. These are the first quantitative results to be presented on the distributions in amplitude and time scale for young aperiodic variables, particularly those varying on timescales of weeks to months.

Volatiles in Protoplanetary Disks

Planets are assembled from the gas, dust, and ice in the accretion disks that encircle young stars. Ices of chemical compounds with low condensation temperatures (<200 K), the so-called volatiles, dominate the solid mass reservoir from which planetesimals are formed and are thus available to build the protoplanetary cores of gas/ice giant planets. It has long been thought that the regions near the condensation fronts of volatiles are preferential birth sites of planets. Moreover, the main volatiles in disks are also the main C-and O-containing species in (exo)planetary atmospheres. Understanding the distribution of volatiles in disks and their role in planet-formation processes is therefore of great interest.

This thesis addresses two fundamental questions concerning the nature of volatiles in planet-forming disks: (1) how are volatiles distributed throughout a disk, and (2) how can we use volatiles to probe planet-forming processes in disks? We tackle the first question in two complementary ways. We have developed a novel super-resolution method to constrain the radial distribution of volatiles throughout a disk by combining multi-wavelength spectra. Thanks to the ordered velocity and temperature profiles in disks, we find that detailed constraints can be derived even with spatially and spectrally unresolved data -- provided a wide range of energy levels are sampled. We also employ high-spatial resolution interferometric images at (sub)mm frequencies using the Atacama Large Millimeter Array (ALMA) to directly measure the radial distribution of volatiles.

For the second question, we combine volatile gas emission measurements with those of the dust continuum emission or extinction to understand dust growth mechanisms in disks and disk instabilities at planet-forming distances from the central star. Our observations and models support the idea that the water vapor can be concentrated in regions near its condensation front at certain evolutionary stages in the lifetime of protoplanetary disks, and that fast pebble growth is likely to occur near the condensation fronts of various volatile species.