Molecular line emission from a protoplanetary disk irradiated externally by a nearby massive star

Star formation often occurs within or nearby stellar clusters. Irradiation by nearby massive stars can photoevaporate protoplanetary disks around young stars (so-called proplyds) which raises questions regarding the ability of planet formation to take place in these environments. We investigate the two-dimensional physical and chemical structure of a protoplanetary disk surrounding a low-mass (T Tauri) star which is irradiated by a nearby massive O-type star to determine the survivability and observability of molecules in proplyds. Compared with an isolated star-disk system, the gas temperature ranges from a factor of a few (in the disk midplane) to around two orders of magnitude (in the disk surface) higher in the irradiated disk. Although the UV flux in the outer disk, in particular, is several orders of magnitude higher, the surface density of the disk is sufficient for effective shielding of the disk midplane so that the disk remains predominantly molecular in nature. We also find that non-volatile molecules, such as HCN and H2O, are able to freeze out onto dust grains in the disk midplane so that the formation of icy planetesimals, e.g., comets, may also be possible in proplyds. We have calculated the molecular line emission from the disk assuming LTE and determined that multiple transitions of atomic carbon, CO (and isotopologues, 13CO and C18O), HCO+, CN, and HCN may be observable with ALMA, allowing characterization of the gas column density, temperature, and optical depth in proplyds at the distance of Orion (˜400 pc).

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.

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.

Chemical evolution of protoplanetary disks - the effects of viscous accretion, turbulent mixing and disk winds

We calculate the chemical evolution of protoplanetary disks considering radial viscous accretion, vertical turbulent mixing, and vertical disk winds. We study the effects on the disk chemical structure when different models for the formation of molecular hydrogen on dust grains are adopted. Our gas-phase chemistry is extracted from the UMIST Database for Astrochemistry (Rate06) to which we have added detailed gas-grain interactions. We use our chemical model results to generate synthetic near- and mid-infrared local thermodynamic equilibrium line emission spectra and compare these with recent Spitzer observations. Our results show that if H2 formation on warm grains is taken into consideration, the H2O and OH abundances in the disk surface increase significantly. We find that the radial accretion flow strongly influences the molecular abundances, with those in the cold midplane layers particularly affected. On the other hand, we show that diffusive turbulent mixing affects the disk chemistry in the warm molecular layers, influencing the line emission from the disk and subsequently improving agreement with observations. We find that NH3, CH3OH, C2H2, and sulfur-containing species are greatly enhanced by the inclusion of turbulent mixing. We demonstrate that disk winds potentially affect the disk chemistry and the resulting molecular line emission in a manner similar to that found when mixing is included.

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.

Molecular Hydrogen Emission from Protoplanetary Disks

We have modeled self-consistently the density and temperature profiles of gas and dust in protoplanetary disks, taking into account irradiation from a central star. Making use of this physical structure, we have calculated the level populations of molecular hydrogen and the line emission from the disks. As a result, we can reproduce the observed strong line spectra of molecular hydrogen from protoplanetary disks, both in the ultraviolet (UV) and the near-infrared, but only if the central star has a strong UV excess radiation.

Molecular Hydrogen Emission from Protoplanetary Disks: Effects of X-ray Irradiation and Dust Evolution

Detailed models for the density and temperature profiles of gas and dust in protoplanetary disks are constructed by taking into account X-ray and UV irradiation from a central T Tauri star, as well as dust size growth and settling toward the disk midplane. The spatial and size distributions of dust grains are numerically computed by solving the coagulation equation for settling dust particles, with the result that the mass and total surface area of dust grains per unit volume of the gas in the disks are very small, except at the midplane. The H2 level populations and line emission are calculated using the derived physical structure of the disks. X-ray irradiation is the dominant heating source of the gas in the inner disk and in the surface layer, while the UV heating dominates otherwise. If the central star has strong X-ray and weak UV radiation, the H2 level populations are controlled by X-ray pumping, and the X-rayinduced transition lines could be observable. If the UV irradiation is strong, the level populations are controlled by thermal collisions or UV pumping, depending on the dust properties. As the dust particles evolve in the disks, the gas temperature at the disk surface drops because the grain photoelectric heating becomes less efficient. This makes the level populations change from LTE to non-LTE distributions, which results in changes to the line ratios. Our results suggest that dust evolution in protoplanetary disks could be observable through the H2 line ratios. The emission lines are strong from disks irradiated by strong UV and X-rays and possessing small dust grains; such disks will be good targets in which to observe H2 emission.

Effects of Accretion Flow on the Chemical Structure in the Inner Regions of Protoplanetary Disks

Aims. We study the dependence of the profiles of molecular abundances and line emission on the accretion flow in the hot (100 K) inner region of protoplanetary disks.
Methods. The gas-phase reactions initiated by evaporation of the ice mantle on dust grains are calculated along the accretion flow. We focus on methanol, a molecule that is formed predominantly by the evaporation of warm ice mantles, to demonstrate how its abundance profile and line emission depend on the accretion flow.
Results. Our results indicate that some evaporated molecules retain high abundances only when the accretion velocity is sufficiently high, and that methanol could be useful as a diagnostic of the accretion flow by means of ALMA observations at the disk radius of 10 AU.

Water in protoplanetary disks

Infrared water line emission from protoplanetary disks, recently observed by the Spitzer and Herschel space telescopes, is thought to trace the surface layer of the inner to outer regions of the disks. We have modelled the water abundance profile and line emission, especially focusing on the effects of dust size growth and turbulent mixing. Comparison between model calculations and observations suggests a small grain model with turbulent mixing is preferred. Copyright © International Astronomical Union 2014.

Chemical Processes in Protoplanetary Disks

We have developed a high-resolution combined physical and chemical model of a protoplanetary disk surrounding a typical T Tauri star. Our aims were to use our model to calculate the chemical structure of disks on small scales (submilliarcsecond in the inner disk for objects at the distance of Taurus, ~140 pc) to investigate the various chemical processes thought to be important in disks and to determine potential molecular tracers of each process. Our gas-phase network was extracted from the UMIST Database for Astrochemistry to which we added gas–grain interactions including freezeout and thermal and non-thermal desorption (cosmic-ray-induced desorption, photodesorption, and X-ray desorption), and a grain-surface network. We find that cosmic-ray-induced desorption has the least effect on our disk chemical structure while photodesorption has a significant effect, enhancing the abundances of most gas-phase molecules throughout the disk and affecting the abundances and distribution of HCN, CN, and CS, in particular. In the outer disk, we also see enhancements in the abundances of H2O and CO2. X-ray desorption is a potentially powerful mechanism in disks, acting to homogenize the fractional abundances of gas-phase species across the depth and increasing the column densities of most molecules, although there remain significant uncertainties in the rates adopted for this process. The addition of grain-surface chemistry enhances the fractional abundances of several small complex organic molecules including CH3OH, HCOOCH3, and CH3OCH3 to potentially observable values (i.e., a fractional abundance of greater than 10-11).

Complex organic molecules along the accretion flow in isolated and externally irradiated protoplanetary disks

(Abridged) The birth environment of the Sun will have influenced the conditions in the pre-solar nebula, including the attainable chemical complexity, important for prebiotic chemistry. The formation and distribution of complex organic molecules (COMs) in a disk around a T Tauri star is investigated for two scenarios: (i) an isolated disk, and (ii) a disk irradiated externally by a nearby massive star. The chemistry is calculated along the accretion flow from the outer disk inwards using a comprehensive network. Two simulations are performed, one beginning with complex ices and one with simple ices only. For the isolated disk, COMs are transported without major alteration into the inner disk where they thermally desorb into the gas reaching an abundance representative of the initial assumed ice abundance. For simple ices, COMs efficiently form on grain surfaces under the conditions in the outer disk. Gas-phase COMs are released into the molecular layer via photodesorption. For the irradiated disk, complex ices are also transported inwards; however, they undergo thermal processing caused by the warmer conditions in the irradiated disk which tends to reduce their abundance along the accretion flow. For simple ices, grain-surface chemistry cannot synthesise COMs in the outer disk because the necessary grain-surface radicals, which tend to be particularly volatile, are not sufficiently abundant on the grain surfaces. Gas-phase COMs are formed in the inner region of the irradiated disk via gas-phase chemistry induced by the desorption of strongly bound molecules such as methanol; hence, the abundances are not representative of the initial molecular abundances injected into the outer disk. These results suggest that the composition of comets formed in isolated disks may differ from those formed in externally irradiated disks with the latter composed of more simple ices.

The VLT-FLAMES tarantula survey:XII. Rotational velocities of the single O-type stars

Context. The 30 Doradus (30 Dor) region of the Large Magellanic Cloud, also known as the Tarantula nebula, is the nearest starburst region. It contains the richest population of massive stars in the Local Group, and it is thus the best possible laboratory to investigate open questions on the formation and evolution of massive stars. Aims. Using ground-based multi-object optical spectroscopy obtained in the framework of the VLT-FLAMES Tarantula Survey (VFTS), we aim to establish the (projected) rotational velocity distribution for a sample of 216 presumably single O-type stars in 30 Dor. The sample is large enough to obtain statistically significant information and to search for variations among subpopulations - in terms of spectral type, luminosity class, and spatial location - in the field of view. Methods. We measured projected rotational velocities, 3e sin i, by means of a Fourier transform method and a profile fitting method applied to a set of isolated spectral lines. We also used an iterative deconvolution procedure to infer the probability density, P(3e), of the equatorial rotational velocity, 3e. Results. The distribution of 3e sin i shows a two-component structure: a peak around 80 km s1 and a high-velocity tail extending up to 600 km s^-1 This structure is also present in the inferred distribution P(3e) with around 80% of the sample having 0 <3e ≤ 300 km s^-1 and the other 20% distributed in the high-velocity region. The presence of the low-velocity peak is consistent with what has been found in other studies for late O- and early B-type stars. Conclusions. Most of the stars in our sample rotate with a rate less than 20% of their break-up velocity. For the bulk of the sample, mass loss in a stellar wind and/or envelope expansion is not efficient enough to significantly spin down these stars within the first few Myr of evolution. If massive-star formation results in stars rotating at birth with a large portion of their break-up velocities, an alternative braking mechanism, possibly magnetic fields, is thus required to explain the present-day rotational properties of the O-type stars in 30 Dor. The presence of a sizeable population of fast rotators is compatible with recent population synthesis computations that investigate the influence of binary evolution on the rotation rate of massive stars. Even though we have excluded stars that show significant radial velocity variations, our sample may have remained contaminated by post-interaction binary products. That the highvelocity tail may be populated primarily (and perhaps exclusively) by post-binary interaction products has important implications for the evolutionary origin of systems that produce gamma-ray bursts. © 2013 Author(s).

The Role of Edges and Line-Ends in Illusory Contour Formation

Illusory contours can be induced along directions approximately collinear to edges or approximately perpendicular to the ends of lines. Using a rating scale procedure we explored the relation between the two types of inducers by systematically varying the thickness of inducing elements to result; in varying amounts of "edge-like" or "line-like" induction. Inducers for om illusory figures consisted of concentric rings with arcs missing. Observers judged the clarity and brightness of illusory figures as the number of arcs, their thicknesses, and spacings were parametrically varied. Degree of clarity and amount of induced brightness were both found to be inverted-U functions of the number of arcs. These results mandate that any valid model of illusory contour formation must account for interference effects between parallel lines or between those neural units responsible for completion of boundary signals in directions perpendicular to the ends of thin lines. Line width was found to have an effect on both clarity and brightness, a finding inconsistent with those models which employ only completion perpendicular to inducer orientation.

The Role of Edges and Line-Ends in Illusory Contour Formation

Illusory contours can be induced along direction approximately collinear to edges or approximately perpendicular to the ends of lines. Using a rating scale procedure we explored the relation between the two types of inducers by systematically varying the thickness of inducing elements to result in varying amounts of "edge-like" or "line-like" induction. Inducers for our illusory figures consisted of concentric rings with arcs missing. Observers judged the clarity and brightness of illusory figures as the number of arcs, their thicknesses, and spacings were parametrically varied. Degree of clarity and amount of induced brightness were both found to be inverted-U functions of the number of arcs. These results mandate that any valid model of illusory contour formation must account for interference effects between parallel lines or between those neural units responsible for completion of boundary signals in directions perpendicular to the ends of thin lines. Line width was found to have an efFect on both clarity and brightness, a finding inconsistent with those models which employ only completion perpendicular to inducer orientation.