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. II. On the importance of photochemistry and X-ray ionization

We investigate the impact of photochemistry and X-ray ionization on the molecular composition of, and ionization fraction in, a protoplanetary disk surrounding a typical T Tauri star. We use a sophisticated physical model, which includes a robust treatment of the radiative transfer of UV and X-ray radiation, and calculate the time-dependent chemical structure using a comprehensive chemical network. In previous work, we approximated the photochemistry and X-ray ionization; here, we recalculate the photoreaction rates using the explicit UV wavelength spectrum and wavelength-dependent reaction cross sections. We recalculate the X-ray ionization rate using our explicit elemental composition and X-ray energy spectrum. We find that photochemistry has a larger influence on the molecular composition than X-ray ionization. Observable molecules sensitive to the photorates include OH, HCO+, N2H+, H2O, CO2, and CH3OH. The only molecule significantly affected by the X-ray ionization is N2H+, indicating that it is safe to adopt existing approximations of the X-ray ionization rate in typical T Tauri star-disk systems. The recalculation of the photorates increases the abundances of neutral molecules in the outer disk, highlighting the importance of taking into account the shape of the UV spectrum in protoplanetary disks. A recalculation of the photoreaction rates also affects the gas-phase chemistry due to the adjustment of the H/H2 and C+/C ratios. The disk ionization fraction is not significantly affected by the methods adopted to calculate the photochemistry and X-ray ionization. We determine that there is a probable "dead zone" where accretion is suppressed, present in a layer, Z/R lsim 0.1-0.2, in the disk midplane, within R ˜ 200 AU.

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.

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.

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).

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.

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.

A SuperWASP search for additional transiting planets in 24 known systems

We present results from a search for additional transiting planets in 24 systems already known to contain a transiting planet. We model the transits due to the known planet in each system and subtract these models from light curves obtained with the SuperWASP (Wide Angle Search for Planets) survey instruments. These residual light curves are then searched for evidence of additional periodic transit events. Although we do not find any evidence for additional planets in any of the planetary systems studied, we are able to characterize our ability to find such planets by means of Monte Carlo simulations. Artificially generated transit signals corresponding to planets with a range of sizes and orbital periods were injected into the SuperWASP photometry and the resulting light curves searched for planets. As a result, the detection efficiency as a function of both the radius and orbital period of any second planet is calculated. We determine that there is a good (>50 per cent) chance of detecting additional, Saturn-sized planets in P ~ 10 d orbits around planet-hosting stars that have several seasons of SuperWASP photometry. Additionally, we confirm previous evidence of the rotational stellar variability of WASP-10, and refine the period of rotation. We find that the period of the rotation is 11.91 +/- 0.05 d, and the false alarm probability for this period is extremely low (~10-13).

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).

The Kepler-10 Planetary System Revisited by HARPS-N: A Hot Rocky World and a Solid Neptune-Mass Planet

Kepler-10b was the first rocky planet detected by the Kepler satellite and confirmed with radial velocity follow-up observations from Keck-HIRES. The mass of the planet was measured with a precision of around 30%, which was
insufficient to constrain models of its internal structure and composition in detail. In addition to Kepler-10b, a second planet transiting the same star with a period of 45 days was statistically validated, but the radial velocities were only
good enough to set an upper limit of 20 M⊕ for the mass of Kepler-10c. To improve the precision on the mass for planet b, the HARPS-N Collaboration decided to observe Kepler-10 intensively with the HARPS-N spectrograph
on the Telescopio Nazionale Galileo on La Palma. In total, 148 high-quality radial-velocity measurements were obtained over two observing seasons. These new data allow us to improve the precision of the mass determination for Kepler-10b to 15%. With a mass of 3.33 ± 0.49 M⊕ and an updated radius of 1.47+0.03 −0.02 R⊕, Kepler-10b has a density of 5.8 ± 0.8 g cm−3, very close to the value predicted by models with the same internal structure and composition as the Earth. We were also able to determine a mass for the 45-day period planet Kepler-10c, with an even better precision of 11%. With a mass of 17.2 ± 1.9 M⊕ and radius of 2.35+0.09 −0.04 R⊕, Kepler-10c has a density of 7.1 ± 1.0 g cm−3. Kepler-10c appears to be the first strong evidence of a class of more massive solid planets with longer orbital periods

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.