Showing 1–43 of 43 results for author: Ciesla, F

Exoplanet Volatile Carbon Content as a Natural Pathway for Haze Formation

Authors: E. A. Bergin, E. Kempton, M. Hirschmann, S. T. Bastelberger, D. J. Teal, G. A. Blake, F. Ciesla, J. Li

Abstract: We explore terrestrial planet formation with a focus on the supply of solid-state organics as the main source of volatile carbon. For the water-poor Earth, the water ice line, or ice sublimation front, within the planet-forming disk has long been a key focal point. We posit that the soot line, the location where solid-state organics are irreversibly destroyed, is also a key location within the dis… ▽ More We explore terrestrial planet formation with a focus on the supply of solid-state organics as the main source of volatile carbon. For the water-poor Earth, the water ice line, or ice sublimation front, within the planet-forming disk has long been a key focal point. We posit that the soot line, the location where solid-state organics are irreversibly destroyed, is also a key location within the disk. The soot line is closer to the host star than the water snowline and overlaps with the location of the majority of detected exoplanets. In this work, we explore the ultimate atmospheric composition of a body that receives a major portion of its materials from the zone between the soot line and water ice line. We model a silicate-rich world with 0.1% and 1% carbon by mass with variable water content. We show that as a result of geochemical equilibrium, the mantle of these planets would be rich in reduced carbon but have relatively low water (hydrogen) content. Outgassing would naturally yield the ingredients for haze production when exposed to stellar UV photons in the upper atmosphere. Obscuring atmospheric hazes appear common in the exoplanetary inventory based on the presence of often featureless transmission spectra (Kreidberg et al. 2014, Knutson et al. 2014, Libby-Roberts et al. 2020). Such hazes may be powered by the high volatile content of the underlying silicate-dominated mantle. Although this type of planet has no solar system counterpart, it should be common in the galaxy with potential impact on habitability. △ Less

Submitted 8 May, 2023; originally announced May 2023.

Comments: 15 pages, 6 figures, accepted by Astrophysical Journal Letters. ApJL version has 5 figures. We include an extra figure (Figure 6) in this submission which is an artist rendering of a young disk including the soot and water ice lines. Image Credit for Fig. 6: Ari Gea/SayoStudio

The Volatile Carbon to Oxygen Ratio as a Tracer for the Formation Locations of Interstellar Comets

Authors: Darryl Z. Seligman, Leslie A. Rogers, Samuel H. C. Cabot, John W. Noonan, Theodore Kareta, Kathleen E. Mandt, Fred Ciesla, Adam McKay, Adina D. Feinstein, W. Garrett Levine, Jacob L. Bean, Thomas Nordlander, Mark R. Krumholz, Megan Mansfield, Devin J. Hoover, Eric Van Clepper

Abstract: Based on the occurrence rates implied by the discoveries of 1I/`Oumuamua and 2I/Borisov, the forthcoming Rubin Observatory Legacy Survey of Space and Time (LSST) should detect $\ge1$ interstellar objects every year (Hoover et al. 2021). We advocate for future measurements of the production rates of H$_2$O, CO$_2$ and CO in these objects to estimate their carbon to oxygen ratios, which traces forma… ▽ More Based on the occurrence rates implied by the discoveries of 1I/`Oumuamua and 2I/Borisov, the forthcoming Rubin Observatory Legacy Survey of Space and Time (LSST) should detect $\ge1$ interstellar objects every year (Hoover et al. 2021). We advocate for future measurements of the production rates of H$_2$O, CO$_2$ and CO in these objects to estimate their carbon to oxygen ratios, which traces formation locations within their original protoplanetary disks. We review similar measurements for Solar System comets, which indicate formation interior to the CO snowline. By quantifying the relative processing in the interstellar medium and Solar System, we estimate that production rates will not be representative of primordial compositions for the majority of interstellar comets. Preferential desorption of CO and CO$_2$ relative to H$_2$O in the interstellar medium implies that measured C/O ratios represent lower limits on the primordial ratios. Specifically, production rate ratios of ${\rm Q}({\rm CO})/{\rm Q}({\rm H_2O})<.2$ and ${\rm Q}({\rm CO})/{\rm Q}({\rm H_2O})>1$ likely indicate formation interior and exterior to the CO snowline, respectively. The high C/O ratio of 2I/Borisov implies that it formed exterior to the CO snowline. We provide an overview of the currently operational facilities capable of obtaining these measurements that will constrain the fraction of ejected comets that formed exterior to the CO snowline. This fraction will provide key insights into the efficiency of and mechanisms for cometary ejection in exoplanetary systems. △ Less

Submitted 1 June, 2022; v1 submitted 27 April, 2022; originally announced April 2022.

Comments: Accepted for publication at PSJ, 33 pages, 14 figures, preprint for reference at Exoplanets IV Program Number 405.03

Chemical Feedbacks of Pebble Growth: Impacts on CO depletion and C/O ratios

Authors: Eric Van Clepper, Jennifer B. Bergner, Arthur D. Bosman, Edwin Bergin, Fred J. Ciesla

Abstract: Observations of protoplanetary disks have revealed them to be complex and dynamic, with vertical and radial transport of gas and dust occurring simultaneously with chemistry and planet formation. Previous models of protoplanetary disks focused primarily on chemical evolution of gas and dust in a static disk, or dynamical evolution of solids in a chemically passive disk. In this paper, we present a… ▽ More Observations of protoplanetary disks have revealed them to be complex and dynamic, with vertical and radial transport of gas and dust occurring simultaneously with chemistry and planet formation. Previous models of protoplanetary disks focused primarily on chemical evolution of gas and dust in a static disk, or dynamical evolution of solids in a chemically passive disk. In this paper, we present a new 1D method for modelling pebble growth and chemistry simultaneously. Gas and small dust particles are allowed to diffuse vertically, connecting chemistry at all elevations of the disk. Pebbles are assumed to form from the dust present around the midplane, inheriting the composition of ices at this location. We present the results of this model after 1 Myr of disk evolution around a 1$M_\odot$ star at various locations both inside and outside of the CO snowline. We find that for a turbulent disk ($α= 10^{-3}$), CO is depleted from the surface layers of the disk by roughly 1-2 orders of magnitude, consistent with observations of protoplanetary disks. This is achieved by a combination of ice sequestration and decreasing UV opacity, both driven by pebble growth. Further, we find the selective removal of ice species via pebble growth and sequestration can increase gas phase C/O ratios to values of approximately unity. However, our model is unable to produce C/O values of $\sim$1.5-2.0 inferred from protoplanetary disk observations, implying selective sequestration of ice is not sufficient to explain C/O ratios $>1$. △ Less

Submitted 1 February, 2022; originally announced February 2022.

Comments: 20 pages, 16 figures, 3 tables, accepted to ApJ

Thermal Processing of Solids Encountering a Young Jovian Core

Authors: Megan N. Barnett, Fred J. Ciesla

Abstract: Jupiter's enhancement in nitrogen relative to hydrogen when compared to the Sun has been interpreted as evidence that its early formation occurred beyond the N$_{2}$ snowline ($\sim$ 20-40 AU). However, the rapid growth necessary to form Jupiter before the dissipation of the solar nebula would lead to the forming planet's core reaching very high temperatures ($>$1000 K), which would lead to it war… ▽ More Jupiter's enhancement in nitrogen relative to hydrogen when compared to the Sun has been interpreted as evidence that its early formation occurred beyond the N$_{2}$ snowline ($\sim$ 20-40 AU). However, the rapid growth necessary to form Jupiter before the dissipation of the solar nebula would lead to the forming planet's core reaching very high temperatures ($>$1000 K), which would lead to it warming its surroundings. Here, we explore the effects of a luminous planetary core on the solids that it ultimately accretes. We find that a critical transition occurs where very hot (rapidly accreting) cores drive off volatiles prior to accretion, while cool cores (slowly accreting) are able to inherit volatile rich solids. Given Jupiter's nitrogen enrichment, if it formed beyond the N$_{2}$ snowline, its core could not have accreted solids at a rate above 10$^{-10}$ M$_{\odot}$ yr$^{-1}$. Our results suggest that either Jupiter formed in more distal regions of the solar nebula, or nitrogen loss was suppressed, either by its incorporation in more refractory carriers or because it was trapped within ices which devolatilized at higher temperatures. △ Less

Submitted 10 January, 2022; v1 submitted 3 January, 2022; originally announced January 2022.

Comments: Accepted to the Astrophysical Journal, 11 pages, 5 figures

Journal ref: ApJ 925, 141 (2022)

Why do M dwarfs have more transiting planets?

Authors: Gijs D. Mulders, Joanna Drążkowska, Nienke van der Marel, Fred J. Ciesla, Ilaria Pascucci

Abstract: We propose a planet formation scenario to explain the elevated occurrence rates of transiting planets around M dwarfs compared to sun-like stars discovered by Kepler. We use a pebble drift and accretion model to simulate the growth of planet cores inside and outside of the snow line. A smaller pebble size interior to the snow line delays the growth of super-Earths, allowing giant planet cores in t… ▽ More We propose a planet formation scenario to explain the elevated occurrence rates of transiting planets around M dwarfs compared to sun-like stars discovered by Kepler. We use a pebble drift and accretion model to simulate the growth of planet cores inside and outside of the snow line. A smaller pebble size interior to the snow line delays the growth of super-Earths, allowing giant planet cores in the outer disk to form first. When those giant planets reach pebble isolation mass they cut off the flow of pebbles to the inner disk and prevent the formation of close-in super-Earths. We apply this model to stars with masses between 0.1 and 2 solar mass and for a range of initial disk masses. We find that the masses of hot super-Earths and of cold giant planets are anti-correlated. The fraction of our simulations that form hot super-Earths is higher around lower-mass stars and matches the exoplanet occurrence rates from Kepler. The fraction of simulations forming cold giant planets is consistent with the stellar mass dependence from radial velocity surveys. A key testable prediction of the pebble accretion hypothesis is that the occurrence rates of super-Earths should decrease again for M dwarfs near the sub-stellar boundary like Trappist-1. △ Less

Submitted 6 October, 2021; originally announced October 2021.

Comments: Published in ApJ Letters

Journal ref: The Astrophysical Journal Letters, Volume 920, Number 1 (2021)

The Mass Budgets and Spatial Scales of Exoplanet Systems and Protoplanetary Disks

Authors: Gijs D. Mulders, Ilaria Pascucci, Fred J. Ciesla, Rachel B. Fernandes

Abstract: Planets are born from disks of gas and dust, and observations of protoplanetary disks are used to constrain the initial conditions of planet formation. However, dust mass measurements of Class II disks with ALMA have called into question whether they contain enough solids to build the exoplanets that have been detected to date. In this paper, we calculate the mass and spatial scale of solid materi… ▽ More Planets are born from disks of gas and dust, and observations of protoplanetary disks are used to constrain the initial conditions of planet formation. However, dust mass measurements of Class II disks with ALMA have called into question whether they contain enough solids to build the exoplanets that have been detected to date. In this paper, we calculate the mass and spatial scale of solid material around Sun-like stars probed by transit and radial velocity exoplanet surveys, and compare those to the observed dust masses and sizes of Class II disks in the same stellar mass regime. We show that the apparent mass discrepancy disappears when accounting for observational selection and detection biases. We find a discrepancy only when the planet formation efficiency is below 100%, or if there is a population of undetected exoplanets that significantly contributes to the mass in solids. We identify a positive correlation between the masses of planetary systems and their respective orbital periods, which is consistent with the trend between the masses and the outer radii of Class II dust disks. This implies that, despite a factor 100 difference in spatial scale, the properties of protoplanetary disks seem to be imprinted on the exoplanet population. △ Less

Submitted 26 July, 2021; originally announced July 2021.

Comments: Accepted for Publication in ApJ

Ice inheritance in dynamical disk models

Authors: Jennifer Bergner, Fred Ciesla

Abstract: The compositions of planet-forming disks are set by a combination of material inherited from the interstellar medium and material reprocessed during disk formation and evolution. Indeed, comets and primitive meteorites exhibit interstellar-like isotopic ratios and/or volatile compositions, supporting that some pristine material was incorporated intact into icy planetesimals in the Solar Nebula. To… ▽ More The compositions of planet-forming disks are set by a combination of material inherited from the interstellar medium and material reprocessed during disk formation and evolution. Indeed, comets and primitive meteorites exhibit interstellar-like isotopic ratios and/or volatile compositions, supporting that some pristine material was incorporated intact into icy planetesimals in the Solar Nebula. To date, the survival of volatile interstellar material in the disk stage has not been modeled using realistic disk physics. Here, we present a modeling framework to track the destruction of interstellar ices on dust grains undergoing transport processes within a disk, with a particular focus on explaining the incorporation of pristine material into icy planetesimals. We find it is difficult to explain inheritance through the local assembly of comets, as ice destruction is rapid for small (<10um) grains in the inner few tens of au. Instead, a plausible pathway to inheritance is to form pebbles at larger disk radii, which then drift inwards to the comet-forming zone with their ices mostly preserved. Small grains beyond ~100 au can experience ice photodissociation at the tens of percent level, however little of the ice is actually lost from the grain, likely making this a robust site for in situ ice chemistry. Our models also indicate that many complex organic species should survive passage through the disk intact. This raises the possibility that organics synthesized in the interstellar medium can be delivered to terrestrial planets by icy body impact and thus potentially participate in origins of life chemistry. △ Less

Submitted 7 July, 2021; originally announced July 2021.

Comments: Accepted to ApJ

Early volatile depletion on planetesimals inferred from C-S systematics of iron meteorite parent bodies

Authors: Marc M. Hirschmann, Edwin A. Bergin, Geoff A. Blake, Fred J. Ciesla, Jie Li

Abstract: During the formation of terrestrial planets, volatile loss may occur through nebular processing, planetesimal differentiation, and planetary accretion. We investigate iron meteorites as an archive of volatile loss during planetesimal processing. The carbon contents of the parent bodies of magmatic iron meteorites are reconstructed by thermodynamic modelling. Calculated solid/molten alloy partition… ▽ More During the formation of terrestrial planets, volatile loss may occur through nebular processing, planetesimal differentiation, and planetary accretion. We investigate iron meteorites as an archive of volatile loss during planetesimal processing. The carbon contents of the parent bodies of magmatic iron meteorites are reconstructed by thermodynamic modelling. Calculated solid/molten alloy partitioning of C increases greatly with liquid S concentration and inferred parent body C concentrations range from 0.0004 to 0.11 wt.%. Parent bodies fall into 2 compositional clusters characterized by cores with medium, and low C/S. Both of these require significant planetesimal degassing, as metamorphic devolatilization on chondrite-like precursors is insufficient to account for their C depletions. Planetesimal core formation models, ranging from closed system extraction to degassing of a wholly molten body, show that significant open system silicate melting and volatile loss is required to match medium and low C/S parent body core compositions. Greater depletion in C relative to S is the hallmark of silicate degassing, indicating that parent body core compositions record processes that affect composite silicate/iron planetesimals. Degassing of bare cores stripped of their silicate mantles would deplete S with negligible C loss, and could not account for inferred parent body core compositions. Devolatilization during small-body differentiation is thus a key process in shaping the volatile inventory of terrestrial planets derived from planetesimals and planetary embryos. △ Less

Submitted 6 April, 2021; originally announced April 2021.

Comments: 50 pages including supplementary materials

Journal ref: Publications of the National Academy of Sciences, March 30, 2021 118 (13)

Earth's carbon deficit caused by early loss through irreversible sublimation

Authors: Jie Li, Edwin A. Bergin, Geoffrey A. Blake, Fred J. Ciesla, Marc M. Hirschmann

Abstract: Carbon is an essential element for life but its behavior during Earth's accretion is not well understood. Carbonaceous grains in meteoritic and cometary materials suggest that irreversible sublimation, and not condensation, governs carbon acquisition by terrestrial worlds. Through astronomical observations and modeling we show that the sublimation front of carbon carriers in the solar nebula, or t… ▽ More Carbon is an essential element for life but its behavior during Earth's accretion is not well understood. Carbonaceous grains in meteoritic and cometary materials suggest that irreversible sublimation, and not condensation, governs carbon acquisition by terrestrial worlds. Through astronomical observations and modeling we show that the sublimation front of carbon carriers in the solar nebula, or the soot line, moved inward quickly so that carbon-rich ingredients would be available for accretion at 1 au after the first million years. On the other hand, geological constraints firmly establish a severe carbon deficit in Earth, requiring the destruction of inherited carbonaceous organics in the majority of its building blocks. The carbon-poor nature of the Earth thus implies carbon loss in its precursor material through sublimation within the first million years. △ Less

Submitted 6 April, 2021; originally announced April 2021.

Comments: 21 pages including main article and supplementary materials

Journal ref: Science Advances published 02 Apr 2021: Vol. 7, no. 14, DOI: 10.1126/sciadv.abd3632

CO Depletion in Protoplanetary Disks: A Unified Picture Combining Physical Sequestration and Chemical Processing

Authors: Sebastiaan Krijt, Arthur D. Bosman, Ke Zhang, Kamber R. Schwarz, Fred J. Ciesla, Edwin A. Bergin

Abstract: The gas-phase CO abundance (relative to hydrogen) in protoplanetary disks decreases by up to 2 orders of magnitude from its ISM value ${\sim}10^{-4}$, even after accounting for freeze-out and photo-dissociation. Previous studies have shown that while local chemical processing of CO and the sequestration of CO ice on solids in the midplane can both contribute, neither of these processes appears cap… ▽ More The gas-phase CO abundance (relative to hydrogen) in protoplanetary disks decreases by up to 2 orders of magnitude from its ISM value ${\sim}10^{-4}$, even after accounting for freeze-out and photo-dissociation. Previous studies have shown that while local chemical processing of CO and the sequestration of CO ice on solids in the midplane can both contribute, neither of these processes appears capable of consistently reaching the observed depletion factors on the relevant timescale of $1{-}3\mathrm{~Myr}$. In this study, we model these processes simultaneously by including a compact chemical network (centered on carbon and oxygen) to 2D ($r+z$) simulations of the outer ($r>20\mathrm{~au}$) disk regions that include turbulent diffusion, pebble formation, and pebble dynamics. In general, we find that the CO/H$_2$ abundance is a complex function of time and location. Focusing on CO in the warm molecular layer, we find that only the most complete model (with chemistry and pebble evolution included) can reach depletion factors consistent with observations. In the absence of pressure traps, highly-efficient planetesimal formation, or high cosmic ray ionization rates, this model also predicts a resurgence of CO vapor interior to the CO snowline. We show the impact of physical and chemical processes on the elemental (C/O) and (C/H) ratios (in the gas and ice phases), discuss the use of CO as a disk mass tracer, and, finally, connect our predicted pebble ice compositions to those of pristine planetesimals as found in the Cold Classical Kuiper Belt and debris disks. △ Less

Submitted 18 July, 2020; originally announced July 2020.

Comments: Accepted for publication in The Astrophysical Journal

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

Authors: Gijs D. Mulders, David P. O'Brien, Fred J. Ciesla, Daniel Apai, Ilaria Pascucci

Abstract: 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 determine… ▽ More 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. △ Less

Submitted 8 May, 2020; originally announced May 2020.

Comments: Resubmitted to ApJ. Planet formation models available online at http://eos-nexus.org/genesis-database/

Tracking Dust Grains During Transport and Growth in Protoplanetary Disks

Authors: William Misener, Sebastiaan Krijt, Fred J. Ciesla

Abstract: Protoplanetary disks are dynamic objects, within which dust grains and gas are expected to be redistributed over large distances. Evidence for this redistribution is seen both in other protoplanetary disks and in our own Solar System, with high-temperature materials thought to originate close to the central star found in the cold, outer regions of the disks. While models have shown this redistribu… ▽ More Protoplanetary disks are dynamic objects, within which dust grains and gas are expected to be redistributed over large distances. Evidence for this redistribution is seen both in other protoplanetary disks and in our own Solar System, with high-temperature materials thought to originate close to the central star found in the cold, outer regions of the disks. While models have shown this redistribution is possible through a variety of mechanisms, these models have generally ignored the possible growth of solids via grain-grain collisions that would occur during transit. Here we investigate the interplay of coagulation and radial and vertical transport of solids in protoplanetary disks, considering cases where growth is limited by bouncing or by fragmentation. We find that in all cases, growth effectively limits the ability for materials to be carried outward or preserved at large distances from the star. This is due to solids being incorporated into large aggregates which drift inwards rapidly under the effects of gas drag. We discuss the implications for mixing in protoplanetary disks, and how the preservation of high temperature materials in outer disks may require structures or outward flow patterns to avoid them being lost via radial drift. △ Less

Submitted 1 October, 2019; originally announced October 2019.

Comments: 16 pages, 8 figures, 1 table, accepted for publication in The Astrophysical Journal

Journal ref: Astrophysical Journal (2019) 885:118

Systematic Variations of CO Gas Abundance with Radius in Gas-rich Protoplanetary Disks

Authors: Ke Zhang, Edwin A. Bergin, Kamber R. Schwarz, Sebastiaan Krijt, Fred Ciesla

Abstract: CO is the most widely used gas tracer of protoplanetary disks. Its abundance is usually assumed to be an interstellar ratio throughout the warm molecular layer of the disk. But recent observations of low CO gas abundance in many protoplanetary disks challenge our understanding of physical and chemical evolutions in disks. Here we investigate the CO abundance structures in four well-studied disks a… ▽ More CO is the most widely used gas tracer of protoplanetary disks. Its abundance is usually assumed to be an interstellar ratio throughout the warm molecular layer of the disk. But recent observations of low CO gas abundance in many protoplanetary disks challenge our understanding of physical and chemical evolutions in disks. Here we investigate the CO abundance structures in four well-studied disks and compare their structures with predictions of chemical processing of CO and transport of CO ice-coated dust grains in disks. We use spatially resolved CO isotopologue line observations and detailed thermo-chemical models to derive CO abundance structures. We find that the CO abundance varies with radius by an order of magnitude in these disks. We show that although chemical processes can efficiently reduce the total column of CO gas within 1 Myr under an ISM level of cosmic-ray ionization rate, the depletion mostly occurs at the deep region of a disk. Without sufficient vertical mixing, the surface layer is not depleted enough to reproduce weak CO emissions observed. The radial profiles of CO depletion in three disks are qualitatively consistent with predictions of pebble formation, settling, and drifting in disks. But the dust evolution alone cannot fully explain the high depletion observed in some disks. These results suggest that dust evolution may play a significant role in transporting volatile materials and a coupled chemical-dynamical study is necessary to understand what raw materials are available for planet formation at different distances from the central star. △ Less

Submitted 8 August, 2019; originally announced August 2019.

Comments: 17 pages, 8 figures, accepted for publication in the ApJ

Dust Transport and Processing in Centrifugally Driven Protoplanetary Disk Winds

Authors: Steven Giacalone, Seth Teitler, Arieh Königl, Sebastiaan Krijt, Fred J. Ciesla

Abstract: There is evidence that protoplanetary disks--including the protosolar one--contain crystalline dust grains on spatial scales where the dust temperature is lower than the threshold value for their formation through thermal annealing of amorphous interstellar silicates. We interpret these observations in terms of an extended, magnetocentrifugally driven disk wind that transports grains from the inne… ▽ More There is evidence that protoplanetary disks--including the protosolar one--contain crystalline dust grains on spatial scales where the dust temperature is lower than the threshold value for their formation through thermal annealing of amorphous interstellar silicates. We interpret these observations in terms of an extended, magnetocentrifugally driven disk wind that transports grains from the inner disk--where they are thermally processed by the stellar radiation after being uplifted from the disk surfaces--to the outer disk regions. For any disk radius $r$ there is a maximum grain size $a_\mathrm{max}(r)$ that can be uplifted from that location: grains of size $a$$\ll$$a_\mathrm{max}$ are carried away by the wind, whereas those with $a$$\lesssim$$a_\mathrm{max}$ reenter the disk at larger radii. A significant portion of the reentering grains converge to--and subsequently accumulate in--a narrow region just beyond $r_\mathrm{max}(a)$, the maximum radius from which grains of size $a$ can be uplifted. We show that this model can account for the inferred crystallinity fractions in classical T Tauri and Herbig Ae disks and for their indicated near constancy after being established early in the disk evolution. It is also consistent with the reported radial gradients in the mean grain size, crystallinity, and crystal composition. In addition, this model yields the properties of the grains that remain embedded in the outflows from protoplanetary disks and naturally explains the inferred persistence of small grains in the surface layers of these disks. △ Less

Submitted 10 July, 2019; originally announced July 2019.

Comments: ApJ, in press

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

Authors: Gijs D. Mulders, Christoph Mordasini, Ilaria Pascucci, Fred J. Ciesla, Alexandre Emsenhuber, Daniel Apai

Abstract: 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 constrai… ▽ More 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 △ Less

Submitted 13 December, 2019; v1 submitted 21 May, 2019; originally announced May 2019.

Comments: Accepted in ApJ

The Efficiency of Noble Gas Trapping in Astrophysical Environments

Authors: Fred Ciesla, Sebastiaan Krijt, Reika Yokochi, Scott Sandford

Abstract: Amorphous ice has long been invoked as a means for trapping extreme volatiles into solids, explaining the abundances of these species in comets and planetary atmospheres. Experiments have shown that such trapping is possible and have been used to estimate the abundances of each species in primitive ices after they formed. However, these experiments have been carried out at deposition rates which e… ▽ More Amorphous ice has long been invoked as a means for trapping extreme volatiles into solids, explaining the abundances of these species in comets and planetary atmospheres. Experiments have shown that such trapping is possible and have been used to estimate the abundances of each species in primitive ices after they formed. However, these experiments have been carried out at deposition rates which exceed those expected in a molecular cloud or solar nebula by many orders of magnitude. Here we develop a numerical model which reproduces the experimental results and apply it to those conditions expected in molecular clouds and protoplanetary disks. We find that two regimes of ice trapping exist: `burial trapping' where the ratio of trapped species to water in the ice reflects that same ratio in the gas and `equilibrium trapping' where the ratio in the ice depends only on the partial pressure of the trapped species in the gas. The boundary between these two regimes is set by both the temperature and rate of ice deposition. Such effects must be accounted for when determining the source of trapped volatiles during planet formation. △ Less

Submitted 14 November, 2018; originally announced November 2018.

Comments: 13 pages, 6 figures

Journal ref: Astrophysical Journal (2018) v867 #2

Transport of CO in Protoplanetary Disks: Consequences of Pebble Formation, Settling, and Radial Drift

Authors: Sebastiaan Krijt, Kamber R. Schwarz, Edwin A. Bergin, Fred J. Ciesla

Abstract: Current models of (exo)planet formation often rely on a large influx of so-called `pebbles' from the outer disk into the planet formation region. In this paper, we investigate how the formation of pebbles in the cold outer regions of protoplanetary disks and their subsequent migration to the inner disk can alter the gas-phase CO distribution both interior and exterior to the midplane CO snowline.… ▽ More Current models of (exo)planet formation often rely on a large influx of so-called `pebbles' from the outer disk into the planet formation region. In this paper, we investigate how the formation of pebbles in the cold outer regions of protoplanetary disks and their subsequent migration to the inner disk can alter the gas-phase CO distribution both interior and exterior to the midplane CO snowline. By simulating the resulting CO abundances in the midplane as well as the warm surface layer, we identify observable signatures of large-scale pebble formation and migration that can be used as `smoking guns' for these important processes. Specifically, we find that after $1\mathrm{~Myr}$, the formation and settling of icy pebbles results in the removal of up to $80\%$ of the CO vapor in the warm ($T>22\mathrm{~K}$) disk layers outside the CO snowline, while the radial migration of pebbles results in the generation of a plume of CO vapor interior the snowline, increasing the CO abundance by a factor ${\sim}2{-}6$ depending on the strength of the turbulence and the sizes of the individual pebbles. The absence of this plume of CO vapor in young nearby disks could indicate efficient conversion of CO into a more refractory species, or a reduction in the radial mass flux of pebbles by, for example, disk inhomogeneities or early planetesimal formation. △ Less

Submitted 6 August, 2018; originally announced August 2018.

Comments: Accepted for publication in The Astrophysical Journal

The Exoplanet Population Observation Simulator. I - The Inner Edges of Planetary Systems

Authors: Gijs D. Mulders, Ilaria Pascucci, Daniel Apai, Fred J. Ciesla

Abstract: The Kepler survey provides a statistical census of planetary systems out to the habitable zone. Because most planets are non-transiting, orbital architectures are best estimated using simulated observations of ensemble populations. Here, we introduce EPOS, the Exoplanet Population Observation Simulator, to estimate the prevalence and orbital architectures of multi-planet systems based on the lates… ▽ More The Kepler survey provides a statistical census of planetary systems out to the habitable zone. Because most planets are non-transiting, orbital architectures are best estimated using simulated observations of ensemble populations. Here, we introduce EPOS, the Exoplanet Population Observation Simulator, to estimate the prevalence and orbital architectures of multi-planet systems based on the latest Kepler data release, DR25. We estimate that at least 42% of sun-like stars have nearly coplanar planetary systems with 7 or more exoplanets. The fraction of stars with at least one planet within 1 au could be as high as 100% depending on assumptions about the distribution of single transiting planets. We estimate an occurrence rate of planets in the habitable zone around sun-like stars of eta_earth=36+-14%. The innermost planets in multi-planet systems are clustered around an orbital period of 10 days (0.1 au), reminiscent of the protoplanetary disk inner edge or could be explained by a planet trap at that location. Only a small fraction of planetary systems have the innermost planet at long orbital periods, with fewer than ~8% and ~3% having no planet interior to the orbit of Mercury and Venus, respectively. These results reinforce the view that the solar system is not a typical planetary system, but an outlier among the distribution of known exoplanetary systems. We predict that at least half of the habitable zone exoplanets are accompanied by (non-transiting) planets at shorter orbital periods, hence knowledge of a close-in exoplanet could be used as a way to optimize the search for Earth-size planets in the Habitable Zone with future direct imaging missions. △ Less

Submitted 21 May, 2018; originally announced May 2018.

Comments: Accepted in AAS journals, code available on github

A comprehensive understanding of planet formation is required for assessing planetary habitability and for the search for life

Authors: Dániel Apai, Fred Ciesla, Gijs D. Mulders, Ilaria Pascucci, Richard Barry, Klaus Pontoppidan, Edwin Bergin, Alex Bixel, Sean Brittain, Shawn D. Domagal-Goldman, Yasuhiro Hasegawa, Hannah Jang-Condell, Renu Malhotra, Michael R. Meyer, Andrew Youdin, Johanna Teske, Neal Turner

Abstract: Dozens of habitable zone, approximately earth-sized exoplanets are known today. An emerging frontier of exoplanet studies is identifying which of these habitable zone, small planets are actually habitable (have all necessary conditions for life) and, of those, which are earth-like. Many parameters and processes influence habitability, ranging from the orbit through detailed composition including v… ▽ More Dozens of habitable zone, approximately earth-sized exoplanets are known today. An emerging frontier of exoplanet studies is identifying which of these habitable zone, small planets are actually habitable (have all necessary conditions for life) and, of those, which are earth-like. Many parameters and processes influence habitability, ranging from the orbit through detailed composition including volatiles and organics, to the presence of geological activity and plate tectonics. While some properties will soon be directly observable, others cannot be probed by remote sensing for the foreseeable future. Thus, statistical understanding of planetary systems' formation and evolution is a key supplement to the direct measurements of planet properties. Probabilistically assessing parameters we cannot directly measure is essential to reliably assess habitability, to prioritizing habitable-zone planets for follow-up, and for interpreting possible biosignatures. △ Less

Submitted 23 March, 2018; originally announced March 2018.

Comments: White paper submitted to the NAS Committee on Exoplanet Science Strategy

Disk evolution and the fate of water

Authors: Lee Hartmann, Fred Ciesla, Oliver Gressel, Richard Alexander

Abstract: We review the general theoretical concepts and observational constraints on the distribution and evolution of water vapor and ice in protoplanetary disks, with a focus on the Solar System. Water is expected to freeze out at distances greater than 1-3 AU from solar-type central stars; more precise estimates are difficult to obtain due to uncertainties in the complex processes involved in disk evolu… ▽ More We review the general theoretical concepts and observational constraints on the distribution and evolution of water vapor and ice in protoplanetary disks, with a focus on the Solar System. Water is expected to freeze out at distances greater than 1-3 AU from solar-type central stars; more precise estimates are difficult to obtain due to uncertainties in the complex processes involved in disk evolution, including dust growth, settling, and radial drift, and the level of turbulence and viscous dissipation within disks. Interferometric observations are now providing constraints on the positions of CO snow lines, but extrapolation to the unresolved regions where water ice sublimates will require much better theoretical understanding of mass and angular momentum transport in disks as well as more refined comparison of observations with sophisticated disk models. △ Less

Submitted 10 August, 2017; originally announced August 2017.

Comments: 26 pages, 8 figures: in the proceedings of the ISSI workshop on "The delivery of water to proto-planets, planets and satellites", to appear in Space Science Reviews

Destruction of Refractory Carbon in Protoplanetary Disks

Authors: Dana E. Anderson, Edwin A. Bergin, Geoffrey A. Blake, Fred J. Ciesla, Ruud Visser, Jeong-Eun Lee

Abstract: The Earth and other rocky bodies in the inner solar system contain significantly less carbon than the primordial materials that seeded their formation. These carbon-poor objects include the parent bodies of primitive meteorites, suggesting that at least one process responsible for solid-phase carbon depletion was active prior to the early stages of planet formation. Potential mechanisms include th… ▽ More The Earth and other rocky bodies in the inner solar system contain significantly less carbon than the primordial materials that seeded their formation. These carbon-poor objects include the parent bodies of primitive meteorites, suggesting that at least one process responsible for solid-phase carbon depletion was active prior to the early stages of planet formation. Potential mechanisms include the erosion of carbonaceous materials by photons or atomic oxygen in the surface layers of the protoplanetary disk. Under photochemically generated favorable conditions, these reactions can deplete the near-surface abundance of carbon grains and polycyclic aromatic hydrocarbons by several orders of magnitude on short timescales relative to the lifetime of the disk out to radii of ~20-100+ au from the central star depending on the form of refractory carbon present. Due to the reliance of destruction mechanisms on a high influx of photons, the extent of refractory carbon depletion is quite sensitive to the disk's internal radiation field. Dust transport within the disk is required to affect the composition of the midplane. In our current model of a passive, constant-alpha disk, where alpha = 0.01, carbon grains can be turbulently lofted into the destructive surface layers and depleted out to radii of ~3-10 au for 0.1-1 um grains. Smaller grains can be cleared out of the planet-forming region completely. Destruction may be more effective in an actively accreting disk or when considering individual grain trajectories in non-idealized disks. △ Less

Submitted 27 July, 2017; originally announced July 2017.

Comments: 15 pages, 9 figures, Accepted for publication in ApJ

Fast litho-panspermia in the habitable zone of the TRAPPIST-1 system

Authors: Sebastiaan Krijt, Timothy J. Bowling, Richard J. Lyons, Fred J. Ciesla

Abstract: With several short-period, Earth-mass planets in the habitable zone, the TRAPPIST-1 system potentially allows litho-panspermia to take place on very short timescales. We investigate the efficiency and speed of inter-planetary material transfer resulting from impacts onto the habitable zone planets. By simulating trajectories of impact ejecta from their moment of ejection until (re-)accretion, we f… ▽ More With several short-period, Earth-mass planets in the habitable zone, the TRAPPIST-1 system potentially allows litho-panspermia to take place on very short timescales. We investigate the efficiency and speed of inter-planetary material transfer resulting from impacts onto the habitable zone planets. By simulating trajectories of impact ejecta from their moment of ejection until (re-)accretion, we find that transport between the habitable zone planets is fastest for ejection velocities around and just above planetary escape velocity. At these ejection velocities, ${\sim}10\%$ of the ejected material reaches another habitable zone planet within $10^2\mathrm{~yr}$, indicating litho-panspermia can be 4 to 5 orders of magnitude faster in TRAPPIST-1 than in the Solar System. △ Less

Submitted 5 April, 2017; originally announced April 2017.

Comments: Accepted for publication in ApJL

Tracing water vapor and ice during dust growth

Authors: Sebastiaan Krijt, Fred J. Ciesla, Edwin A. Bergin

Abstract: The processes that govern the evolution of dust and water (in the form of vapor or ice) in protoplanetary disks are intimately connected. We have developed a model that simulates dust coagulation, dust dynamics (settling, turbulent mixing), vapor diffusion, and condensation/sublimation of volatiles onto grains in a vertical column of a protoplanetary disk. We employ the model to study how dust gro… ▽ More The processes that govern the evolution of dust and water (in the form of vapor or ice) in protoplanetary disks are intimately connected. We have developed a model that simulates dust coagulation, dust dynamics (settling, turbulent mixing), vapor diffusion, and condensation/sublimation of volatiles onto grains in a vertical column of a protoplanetary disk. We employ the model to study how dust growth and dynamics influence the vertical distribution of water vapor and water ice in the region just outside the radial snowline. Our main finding is that coagulation (boosted by the enhanced stickiness of icy grains) and the ensuing vertical settling of solids results in water vapor being depleted, but not totally removed, from the region above the snowline on a timescale commensurate with the vertical turbulent mixing timescale. Depending on the strength of the turbulence and the temperature, the depletion can reach factors of up to ${\sim}50$ in the disk atmosphere. In our isothermal column, this vapor depletion results in the vertical snowline moving closer to the midplane (by up to 2 gas scale heights) and the gas-phase C/O ratio above the vertical snowline increasing. Our findings illustrate the importance of dynamical effects and the need for understanding coevolutionary dynamics of gas and solids in planet-forming environments. △ Less

Submitted 20 October, 2016; originally announced October 2016.

Comments: Accepted for publication in The Astrophysical Journal

A Steeper than Linear Disk Mass-Stellar Mass Scaling Relation

Authors: I. Pascucci, L. Testi, G. J. Herczeg, F. Long, C. F. Manara, N. Hendler, G. D. Mulders, S. Krijt, F. Ciesla, Th. Henning, S. Mohanty, E. Drabek-Maunder, D. Apai, L. Szucs, G. Sacco, J. Olofsson

Abstract: The disk mass is among the most important input parameter for every planet formation model to determine the number and masses of the planets that can form. We present an ALMA 887micron survey of the disk population around objects from 2 to 0.03Msun in the nearby 2Myr-old Chamaeleon I star-forming region. We detect thermal dust emission from 66 out of 93 disks, spatially resolve 34 of them, and ide… ▽ More The disk mass is among the most important input parameter for every planet formation model to determine the number and masses of the planets that can form. We present an ALMA 887micron survey of the disk population around objects from 2 to 0.03Msun in the nearby 2Myr-old Chamaeleon I star-forming region. We detect thermal dust emission from 66 out of 93 disks, spatially resolve 34 of them, and identify two disks with large dust cavities of about 45AU in radius. Assuming isothermal and optically thin emission, we convert the 887micron flux densities into dust disk masses, hereafter Mdust. We find that the Mdust-Mstar relation is steeper than linear with power law indices 1.3-1.9, where the range reflects two extremes of the possible relation between the average dust temperature and stellar luminosity. By re-analyzing all millimeter data available for nearby regions in a self-consistent way, we show that the 1-3 Myr-old regions of Taurus, Lupus, and Chamaeleon I share the same Mdust-Mstar relation, while the 10Myr-old Upper Sco association has a steeper relation. Theoretical models of grain growth, drift, and fragmentation reproduce this trend and suggest that disks are in the fragmentation-limited regime. In this regime millimeter grains will be located closer in around lower-mass stars, a prediction that can be tested with deeper and higher spatial resolution ALMA observations. △ Less

Submitted 11 August, 2016; originally announced August 2016.

Comments: ApJ, in press

Dust diffusion and settling in the presence of collisions: Trapping (sub)micron grains in the midplane

Authors: Sebastiaan Krijt, Fred J. Ciesla

Abstract: In protoplanetary disks, the distribution and abundance of small (sub)micron grains are important for a range of physical and chemical processes. For example, they dominate the optical depth at short wavelengths and their surfaces are the sites of many important chemical reactions such as the formation of water. Based on their aerodynamical properties (i.e., their strong dynamical coupling with th… ▽ More In protoplanetary disks, the distribution and abundance of small (sub)micron grains are important for a range of physical and chemical processes. For example, they dominate the optical depth at short wavelengths and their surfaces are the sites of many important chemical reactions such as the formation of water. Based on their aerodynamical properties (i.e., their strong dynamical coupling with the surrounding gas) it is often assumed that these small grains are well-mixed with the gas. Our goal is to study the vertical (re)distribution of grains taking into account settling, turbulent diffusion, as well as collisions with other dust grains. Assuming a fragmentation-limited background dust population, we developed a Monte Carlo approach that follows single monomers as they move through a vertical column of gas and become incorporated in different aggregates as they undergo sticking and fragmenting collisions. We find that (sub)micron grains are not necessarily well-mixed vertically, but can become trapped in a thin layer with a scale-height that is significantly smaller than that of the gas. This collisional trapping occurs when the timescale for diffusion is comparable to or longer than the collision timescale in the midplane and its effect is strongest when the most massive particles in the size-distribution show significant settling. Based on simulations and analytical considerations we conclude that for typical dust-to-gas ratios and turbulence levels, the collisional trapping of small grains should be a relatively common phenomenon. The absence of trapping could then indicate a low dust-to-gas ratio, possibly because a large portion of the dust mass has been removed through radial drift or is locked up in planetesimals. △ Less

Submitted 8 March, 2016; originally announced March 2016.

Comments: 19 pages, 9 figures; accepted for publication in ApJ

Direct imaging of the water snow line at the time of planet formation using two ALMA continuum bands

Authors: Andrea Banzatti, Paola Pinilla, Luca Ricci, Klaus M. Pontoppidan, Til Birnstiel, Fred Ciesla

Abstract: Molecular snow lines in protoplanetary disks have been studied theoretically for decades because of their importance in shaping planetary architectures and compositions. The water snow line lies in the planet formation region at < 10 AU, and so far its location has been estimated only indirectly from spatially-unresolved spectroscopy. This work presents a proof-of-concept method to directly image… ▽ More Molecular snow lines in protoplanetary disks have been studied theoretically for decades because of their importance in shaping planetary architectures and compositions. The water snow line lies in the planet formation region at < 10 AU, and so far its location has been estimated only indirectly from spatially-unresolved spectroscopy. This work presents a proof-of-concept method to directly image the water snow line in protoplanetary disks through its physical and chemical imprint in the local dust properties. We adopt a physical disk model that includes dust coagulation, fragmentation, drift, and a change in fragmentation velocities of a factor 10 between dry silicates and icy grains as found by laboratory work. We find that the presence of a water snow line leads to a sharp discontinuity in the radial profile of the dust emission spectral index α_mm, due to replenishment of small grains through fragmentation. We use the ALMA simulator to demonstrate that this effect can be observed in protoplanetary disks using spatially-resolved ALMA images in two continuum bands. We explore the model dependence on the disk viscosity and find that the spectral index reveals the water snow line for a wide range of conditions, with opposite trends when the emission is optically thin rather than thick. If the disk viscosity is low (α_visc < 10^-3) the snow line produces a ring-like structure with a minimum at α_mm ~ 2 in the optically thick regime, possibly similar to what has been measured with ALMA in the innermost region of the HL Tau disk. △ Less

Submitted 20 November, 2015; originally announced November 2015.

Comments: 6 pages, 4 figures, accepted for publication in ApJ Letters

Planet Hunters X. KIC 8462852 - Where's the Flux?

Authors: T. S. Boyajian, D. M. LaCourse, S. A. Rappaport, D. Fabrycky, D. A. Fischer, D. Gandolfi, G. M. Kennedy, H. Korhonen, M. C. Liu, A. Moor, K. Olah, K. Vida, M. C. Wyatt, W. M. J. Best, J. Brewer, F. Ciesla, B. Csak, H. J. Deeg, T. J. Dupuy, G. Handler, K. Heng, S. B. Howell, S. T. Ishikawa, J. Kovacs, T. Kozakis , et al. (24 additional authors not shown)

Abstract: Over the duration of the Kepler mission, KIC8462852 was observed to undergo irregularly shaped, aperiodic dips in flux of up to $\sim 20$\%. The dipping activity can last for between 5 and 80 days. We characterize the object with high-resolution spectroscopy, spectral energy distribution fitting, radial velocity measurements, high-resolution imaging, and Fourier analyses of the Kepler light curve.… ▽ More Over the duration of the Kepler mission, KIC8462852 was observed to undergo irregularly shaped, aperiodic dips in flux of up to $\sim 20$\%. The dipping activity can last for between 5 and 80 days. We characterize the object with high-resolution spectroscopy, spectral energy distribution fitting, radial velocity measurements, high-resolution imaging, and Fourier analyses of the Kepler light curve. We determine that KIC8462852 is a typical main-sequence F3 V star that exhibits no significant IR excess, and has no very close interacting companions. In this paper, we describe various scenarios to explain the dipping events observed in the Kepler light curve. We confirm that the dipping signals in the data are not caused by any instrumental or data processing artifact, and thus are astrophysical in origin. We construct scenario-independent constraints on the size and location of a body in the system that is needed to reproduce the observations. We deliberate over several assorted stellar and circumstellar astrophysical scenarios, most of which have problems explaining the data in hand. By considering the observational constraints on dust clumps in orbit around a normal main-sequence star, we conclude that the scenario most consistent with the data in hand is the passage of a family of exocomet or planetesimal fragments, all of which are associated with a single previous break-up event, possibly caused by tidal disruption or thermal processing. The minimum total mass associated with these fragments likely exceeds $10^{-6}$~\mearth, corresponding to an original rocky body of $>100$~km in diameter. We discuss the necessity of future observations to help interpret the system. △ Less

Submitted 25 January, 2016; v1 submitted 11 September, 2015; originally announced September 2015.

Comments: Accepted for publication in MNRAS. 17 pages, 13 figures

Tracing the Ingredients for a Habitable Earth from Interstellar Space through Planet Formation

Authors: Edwin A. Bergin, Geoffrey A. Blake, Fred Ciesla, Marc M. Hirschmann, Jie Li

Abstract: We use the C/N ratio as a monitor of the delivery of key ingredients of life to nascent terrestrial worlds. Total elemental C and N contents, and their ratio, are examined for the interstellar medium, comets, chondritic meteorites and terrestrial planets; we include an updated estimate for the Bulk Silicate Earth (C/N = 49.0 +/- 9.3). Using a kinetic model of disk chemistry, and the sublimation/co… ▽ More We use the C/N ratio as a monitor of the delivery of key ingredients of life to nascent terrestrial worlds. Total elemental C and N contents, and their ratio, are examined for the interstellar medium, comets, chondritic meteorites and terrestrial planets; we include an updated estimate for the Bulk Silicate Earth (C/N = 49.0 +/- 9.3). Using a kinetic model of disk chemistry, and the sublimation/condensation temperatures of primitive molecules, we suggest that organic ices and macro-molecular (refractory or carbonaceous dust) organic material are the likely initial C and N carriers. Chemical reactions in the disk can produce nebular C/N ratios of ~1-12, comparable to those of comets and the low end estimated for planetesimals. An increase of the C/N ratio is traced between volatile-rich pristine bodies and larger volatile-depleted objects subjected to thermal/accretional metamorphism. The C/N ratios of the dominant materials accreted to terrestrial planets should therefore be higher than those seen in carbonaceous chondrites or comets. During planetary formation, we explore scenarios leading to further volatile loss and associated C/N variations owing to core formation and atmospheric escape. Key processes include relative enrichment of nitrogen in the atmosphere and preferential sequestration of carbon by the core. The high C/N BSE ratio therefore is best satisfied by accretion of thermally processed objects followed by large-scale atmospheric loss. These two effects must be more profound if volatile sequestration in the core is effective. The stochastic nature of these processes hints that the surface/atmospheric abundances of biosphere-essential materials will likely be variable. △ Less

Submitted 16 July, 2015; originally announced July 2015.

Comments: Accepted by PNAS per http://www.pnas.org/content/early/2015/07/01/1500954112.abstract?sid=9fd8abea-9d33-46d8-b755-217d10b1c24b

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

Authors: Gijs D. Mulders, Fred J. Ciesla, Michiel Min, Ilaria Pascucci

Abstract: 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 rate… ▽ More 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. △ Less

Submitted 13 May, 2015; originally announced May 2015.

Comments: 8 pages, 5 figures, accepted by ApJ

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

Authors: Fred J. Ciesla, Gijs D. Mulders, Ilaria Pascucci, Daniel Apai

Abstract: 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 radioi… ▽ More 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. △ Less

Submitted 25 February, 2015; originally announced February 2015.

Comments: Accepted to ApJ

Sulfurization of Iron in the Dynamic Solar Nebula and Implications for Planetary Compositions

Authors: Fred J. Ciesla

Abstract: One explanation for the enhanced ratio of volatiles to hydrogen in Jupiter's atmosphere compared to a a gas of solar composition is that the planet accreted volatile-bearing clathrates during its formation. Models, however, suggest that S would be over abundant if clathrates were the primary carrier of Jupiter's volatiles. This led to the suggestion that S was depleted in the outer nebula due to t… ▽ More One explanation for the enhanced ratio of volatiles to hydrogen in Jupiter's atmosphere compared to a a gas of solar composition is that the planet accreted volatile-bearing clathrates during its formation. Models, however, suggest that S would be over abundant if clathrates were the primary carrier of Jupiter's volatiles. This led to the suggestion that S was depleted in the outer nebula due to the formation troilite (FeS). Here, this depletion is quantitatively explored by modeling the coupled dynamical and chemical evolution of Fe grains in the solar nebula. It is found that disks that undergo rapid radial expansion from an initially compact state may allow sufficient production of FeS and carry H$_{2}$S-depleted gas outward where ices would form, providing the conditions needed for S-depleted clathrates to form. However, this expansion would also carry FeS grains to this region, which could also be incorporated into planetesimals. Thus for clathrates to be a viable source of volatiles, models must account for the presence of both H$_{2}$S in FeS in the outer solar nebula. △ Less

Submitted 8 January, 2015; originally announced January 2015.

Comments: 5 pages, 1 figure; Accepted to ApJ Letters

The Phases of Water Ice in the Solar Nebula

Authors: Fred J. Ciesla

Abstract: Understanding the phases of water ice that were present in the solar nebula has implications for understanding cometary and planetary compositions as well as internal evolution of these bodies. Here we show that amorphous ice formed more readily than previously recognized, with formation at temperatures <70 K being possible under protoplanetary disk conditions. We further argue that photodesorptio… ▽ More Understanding the phases of water ice that were present in the solar nebula has implications for understanding cometary and planetary compositions as well as internal evolution of these bodies. Here we show that amorphous ice formed more readily than previously recognized, with formation at temperatures <70 K being possible under protoplanetary disk conditions. We further argue that photodesorption and freeze-out of water molecules near the surface layers of the solar nebula would have provided the conditions needed for amorphous ice to form. This processing would be a natural consequence of ice dynamics, and would allow for the trapping of noble gases and other volatiles in water ice in the outer solar nebula. △ Less

Submitted 21 February, 2014; originally announced February 2014.

Comments: Accepted for Publication in Astrophysical Journal Letters on 2/12/14

The early impact histories of meteorite parent bodies

Authors: Thomas M Davison, David P O'Brien, Fred J Ciesla, Gareth S Collins

Abstract: We have developed a statistical framework that uses collisional evolution models, shock physics modeling and scaling laws to determine the range of plausible collisional histories for individual meteorite parent bodies. It is likely that those parent bodies that were not catastrophically disrupted sustained hundreds of impacts on their surfaces - compacting, heating, and mixing the outer layers; i… ▽ More We have developed a statistical framework that uses collisional evolution models, shock physics modeling and scaling laws to determine the range of plausible collisional histories for individual meteorite parent bodies. It is likely that those parent bodies that were not catastrophically disrupted sustained hundreds of impacts on their surfaces - compacting, heating, and mixing the outer layers; it is highly unlikely that many parent bodies escaped without any impacts processing the outer few kilometers. The first 10 - 20 Myr were the most important time for impacts, both in terms of the number of impacts and the increase of specific internal energy due to impacts. The model has been applied to evaluate the proposed impact histories of several meteorite parent bodies: up to 10 parent bodies that were not disrupted in the first 100 Myr experienced a vaporizing collision of the type necessary to produce the metal inclusions and chondrules on the CB chondrite parent; around 1 - 5% of bodies that were catastrophically disrupted after 12 Myr sustained impacts at times that match the heating events recorded on the IAB/winonaite parent body; more than 75% of 100 km radius parent bodies which survived past 100 Myr without being disrupted sustained an impact that excavates to the depth required for mixing in the outer layers of the H chondrite parent body; and to protect the magnetic field on the CV chondrite parent body, the crust would have had to have been thick (~ 20 km) in order to prevent it being punctured by impacts. △ Less

Submitted 10 September, 2013; v1 submitted 29 July, 2013; originally announced July 2013.

Comments: 30 pages, 11 figures, 3 tables. Accepted for publication in Meteoritics & Planetary Science

Indication of insensitivity of planetary weathering behavior and habitable zone to surface land fraction

Authors: Dorian S. Abbot, Nicolas B. Cowan, Fred J. Ciesla

Abstract: It is likely that unambiguous habitable zone terrestrial planets of unknown water content will soon be discovered. Water content helps determine surface land fraction, which influences planetary weathering behavior. This is important because the silicate weathering feedback determines the width of the habitable zone in space and time. Here a low-order model of weathering and climate, useful for ga… ▽ More It is likely that unambiguous habitable zone terrestrial planets of unknown water content will soon be discovered. Water content helps determine surface land fraction, which influences planetary weathering behavior. This is important because the silicate weathering feedback determines the width of the habitable zone in space and time. Here a low-order model of weathering and climate, useful for gaining qualitative understanding, is developed to examine climate evolution for planets of various land-ocean fractions. It is pointed out that, if seafloor weathering does not depend directly on surface temperature, there can be no weathering-climate feedback on a waterworld. This would dramatically narrow the habitable zone of a waterworld. Results from our model indicate that weathering behavior does not depend strongly on land fraction for partially ocean-covered planets. This is powerful because it suggests that previous habitable zone theory is robust to changes in land fraction, as long as there is some land. Finally, a mechanism is proposed for a waterworld to prevent complete water loss during a moist greenhouse through rapid weathering of exposed continents. This process is named a "waterworld self-arrest," and it implies that waterworlds can go through a moist greenhouse stage and end up as planets like Earth with partial ocean coverage. This work stresses the importance of surface and geologic effects, in addition to the usual incident stellar flux, for habitability. △ Less

Submitted 8 August, 2012; originally announced August 2012.

Comments: 15 pages, 6 figures, accepted at ApJ

Post-Impact Thermal Evolution of Porous Planetesimals

Authors: T. M. Davison, F. J. Ciesla, G. S. Collins

Abstract: Impacts between planetesimals have largely been ruled out as a heat source in the early Solar System, by calculations that show them to be an inefficient heat source and unlikely to cause global heating. However, the long-term, localized thermal effects of impacts on planetesimals have never been fully quantified. Here, we simulate a range of impact scenarios between planetesimals to determine the… ▽ More Impacts between planetesimals have largely been ruled out as a heat source in the early Solar System, by calculations that show them to be an inefficient heat source and unlikely to cause global heating. However, the long-term, localized thermal effects of impacts on planetesimals have never been fully quantified. Here, we simulate a range of impact scenarios between planetesimals to determine the post-impact thermal histories of the parent bodies, and hence the importance of impact heating in the thermal evolution of planetesimals. We find on a local scale that heating material to petrologic type 6 is achievable for a range of impact velocities and initial porosities, and impact melting is possible in porous material at a velocity of > 4 km/s. Burial of heated impactor material beneath the impact crater is common, insulating that material and allowing the parent body to retain the heat for extended periods (~ millions of years). Cooling rates at 773 K are typically 1 - 1000 K/Ma, matching a wide range of measurements of metallographic cooling rates from chondritic materials. While the heating presented here is localized to the impact site, multiple impacts over the lifetime of a parent body are likely to have occurred. Moreover, as most meteorite samples are on the centimeter to meter scale, the localized effects of impact heating cannot be ignored. △ Less

Submitted 12 July, 2012; originally announced July 2012.

Comments: 38 pages, 9 figures, Revised for Geochimica et Cosmochimica Acta (Sorry, they do not accept LaTeX)

Residence Times of Particles in Diffusive Protoplanetary Disk Environments II. Radial Motions and Applications to Dust Annealing

Authors: Fred J. Ciesla

Abstract: The origin of crystalline grains in comets and the outer regions of protoplanetary disks remains a mystery. It has been suggested that such grains form via annealing of amorphous precursors in the hot, inner region of a protoplanetary disk, where the temperatures needed for such transformations were found, and were then transported outward by some dynamical means. Here we develop a means of tracki… ▽ More The origin of crystalline grains in comets and the outer regions of protoplanetary disks remains a mystery. It has been suggested that such grains form via annealing of amorphous precursors in the hot, inner region of a protoplanetary disk, where the temperatures needed for such transformations were found, and were then transported outward by some dynamical means. Here we develop a means of tracking the paths that dust grains would have taken through a diffusive protoplanetary disk and examine the types and ranges of environments that particles would have seen over a 10$^{6}$ year time period in the dynamic disk. We then combine this model with three annealing laws to examine how the dynamic evolution of amorphous grains would have led to their physical restructuring and their delivery to various regions of the disk. It is found that "sibling particles"-- those particles that reside at the same location at a given period of time--take a wide range of unique and independent paths through the disk to arrive there. While high temperatures can persist in the disk for very long time periods, we find that those grains which are delivered to the cold outer regions of the disk are largely annealed in the first few $\times10^{5}$ yrs of disk history. This suggests that the crystallinity of grains in the outer disk would be determined early and remain unchanged for much of disk history, in agreement with recent astronomical observations. △ Less

Submitted 23 August, 2011; originally announced August 2011.

Comments: 33 pages, 13 figures, To be published in ApJ

Heterogeneous distribution of Al-26 at the birth of the Solar System

Authors: Kentaro Makide, Kazuhide Nagashima, Alexander N. Krot, Gary R. Huss, Fred J. Ciesla, Eric Hellebrand, Eric Gaidos, Le Yang

Abstract: It is believed that Al-26, a short-lived (t1/2 = 0.73 Ma) and now extinct radionuclide, was uniformly distributed in the nascent Solar System with the initial Al-26/Al-27 ratio of ~5.2\times10-5, suggesting its external stellar origin. However, the stellar source of Al-26 and the manner in which it was injected into the solar system remain controversial: the Al-26 could have been produced by an as… ▽ More It is believed that Al-26, a short-lived (t1/2 = 0.73 Ma) and now extinct radionuclide, was uniformly distributed in the nascent Solar System with the initial Al-26/Al-27 ratio of ~5.2\times10-5, suggesting its external stellar origin. However, the stellar source of Al-26 and the manner in which it was injected into the solar system remain controversial: the Al-26 could have been produced by an asymptotic giant branch star, a supernova, or a Wolf-Rayet star and injected either into the protosolar molecular cloud or protoplanetary disk. Corundum (Al2O3) is thermodynamically predicted to be the first condensate from a cooling gas of solar composition. Here we show that micron-sized corundum condensates from O-16-rich gas (Big Delta O-17 ~ -25%) of solar composition recorded heterogeneous distribution of Al-26 at the birth of the solar system: the inferred initial Al-26/Al-27 ratio ranges from ~6.5x10-5 to <2x10-6; ~50% of the corundum grains measured are Al-26-poor. Other Al-26-poor, O-16-rich refractory objects include grossite (CaAl4O7)- and hibonite(CaAl12O19)-rich calcium-aluminum-rich inclusions (CAIs) in CH chondrites, platy hibonite crystals in CM chondrites, and FUN (fractionation and unidentified nuclear isotopic anomalies) CAIs in CV, CO, and CR chondrites. Considering the apparently early and short duration (<0.3 Ma) of condensation of refractory O-16-rich solids in the solar system, we infer that Al-26 was injected into the collapsing protosolar molecular cloud and later homogenized in the protoplanetary disk. The apparent lack of correlation between Al-26 abundance and O-isotope compositions of corundum grains put important constraints on the stellar source of Al-26 in the solar system. △ Less

Submitted 13 May, 2011; originally announced May 2011.

Comments: Astrophysical Journal Letters 733, L31

Residence Times of Particles in Diffusive Protoplanetary Disk Environments I. Vertical Motions

Authors: Fred J. Ciesla

Abstract: The chemical and physical evolution of primitive materials in protoplanetary disks are determined by the types of environments they are exposed to and their residence times within each environment. Here a method for calculating representative paths of materials in diffusive protoplanetary disks is developed and applied to understanding how the vertical trajectories that particles take impact their… ▽ More The chemical and physical evolution of primitive materials in protoplanetary disks are determined by the types of environments they are exposed to and their residence times within each environment. Here a method for calculating representative paths of materials in diffusive protoplanetary disks is developed and applied to understanding how the vertical trajectories that particles take impact their overall evolution. The methods are general enough to be applied to disks with uniform diffusivity, the so-called "constant-$α$" cases, and disks with a spatially varying diffusivity, such as expected in "layered-disks." The average long-term dynamical evolution of small particles and gaseous molecules is independent of the specific form of the diffusivity in that they spend comparable fractions of their lifetimes at different heights in the disk. However, the paths that individual particles and molecules take depend strongly on the form of the diffusivity leading to a different range of behavior of particles in terms of deviations from the mean. As temperatures, gas densities, chemical abundances, and photon fluxes will vary with height in protoplanetary disks, the different paths taken by primitive materials will lead to differences in their chemical and physical evolution. Examples of differences in gas phase chemistry and photochemistry are explored here. The methods outlined here provide a means of integrating particle dynamics with chemical models to understand the formation and evolution of primitive materials in protoplanetary disks over timescales of $10^5$-$10^6$ years. △ Less

Submitted 7 October, 2010; originally announced October 2010.

Comments: Accepted to ApJ

Two-Dimensional Transport of Solids in Viscous Protoplanetary Disks

Authors: Fred Ciesla

Abstract: Large-scale radial transport of solids appears to be a fundamental consequence of protoplanetary disk evolution based on the presence of high temperature minerals in comets and the outer regions of protoplanetary disks around other stars. Further, inward transport of solids from the outer regions of the solar nebula has been postulated to be the manner in which short-lived radionuclides were int… ▽ More Large-scale radial transport of solids appears to be a fundamental consequence of protoplanetary disk evolution based on the presence of high temperature minerals in comets and the outer regions of protoplanetary disks around other stars. Further, inward transport of solids from the outer regions of the solar nebula has been postulated to be the manner in which short-lived radionuclides were introduced to the terrestrial planet region and the cause of the variations in oxygen isotope ratios seen in primitive materials. Here, both outward and inward transport of solids are investigated in the context of a two-dimensional, viscously evolving protoplanetary disk. The dynamics of solids are investigated to determine how they depend on particle size and the particular stage of protoplanetary disk evolution, corresponding to different rates of mass transport. It is found that the outward flows that arise around the disk midplane of a protoplanetary disk aid in the outward transport of solids up to the size of CAIs and can increase the crystallinity fraction of silicate dust at 10 AU around a solar mass star to as much as $\sim$40% in the case of rapidly evolving disks, decreasing as the accretion rate onto the star slows. High velocity, inward flows along the disk surface aid in the rapid transport of solids from the outer disk to the inner disk, particularly for small dust. Despite the diffusion that occurs throughout the disk, the large-scale, meridonal flows associated with mass transport prevent complete homogenization of the disk, allowing compositional gradients to develop that vary in intensity for a timescale of one million years. △ Less

Submitted 5 January, 2009; v1 submitted 19 December, 2008; originally announced December 2008.

Comments: 49 pages, 18 figures. Accepted to Icarus. Corrected Reference

Cooling of Dense Gas by H2O Line Emission and an Assessment of its Effects in Chondrule-Forming Shocks

Authors: M. A. Morris, S. J. Desch, F. J. Ciesla

Abstract: We consider gas at densities appropriate to protoplanetary disks and calculate its ability to cool due to line radiation emitted by H2O molecules within the gas. Our work follows that of Neufeld & Kaufman (1993; ApJ, 418, 263), expanding on their work in several key aspects, including use of a much expanded line database, an improved escape probability formulism, and the inclusion of dust grains… ▽ More We consider gas at densities appropriate to protoplanetary disks and calculate its ability to cool due to line radiation emitted by H2O molecules within the gas. Our work follows that of Neufeld & Kaufman (1993; ApJ, 418, 263), expanding on their work in several key aspects, including use of a much expanded line database, an improved escape probability formulism, and the inclusion of dust grains, which can absorb line photons. Although the escape probabilities formally depend on a complicated combination of optical depth in the lines and in the dust grains, we show that the cooling rate including dust is well approximated by the dust-free cooling rate multiplied by a simple function of the dust optical depth. We apply the resultant cooling rate of a dust-gas mixture to the case of a solar nebula shock pertinent to the formation of chondrules, millimeter-sized melt droplets found in meteorites. Our aim is to assess whether line cooling can be neglected in chondrule-forming shocks or if it must be included. We find that for typical parameters, H2O line cooling shuts off a few minutes past the shock front; line photons that might otherwise escape the shocked region and cool the gas will be absorbed by dust grains. During the first minute or so past the shock, however, line photons will cool the gas at rates ~ 10,000 K/hr, dropping the temperature of the gas (and most likely the chondrules within the gas) by several hundred K. Inclusion of H2O line cooling therefore must be included in models of chondrule formation by nebular shocks. △ Less

Submitted 29 September, 2008; originally announced September 2008.

Comments: Accepted for publication in The Astrophysical Journal

Journal ref: Astrophysical Journal 691 (2009) 320-331

Radial Transport in the Solar Nebula: Implications for Moderately Volatile Element Depletions in Chondritic Meteorites

Authors: Fred J. Ciesla

Abstract: In this paper, the possibility that the moderately volatile element depletions observed in chondritic meteorites are the results of planetesimals accreting in a solar nebula that cooled from an initially hot state (temperatures > 1350 K out to ~2-4 AU) is explored. A model is developed to track the chemical inventory of planetesimals that accrete in a viscously evolving protoplanetary disk, acco… ▽ More In this paper, the possibility that the moderately volatile element depletions observed in chondritic meteorites are the results of planetesimals accreting in a solar nebula that cooled from an initially hot state (temperatures > 1350 K out to ~2-4 AU) is explored. A model is developed to track the chemical inventory of planetesimals that accrete in a viscously evolving protoplanetary disk, accounting for the redistribution of solids and vapor by advection, diffusion, and gas drag. It is found that depletion trends similar to those observed in the chondritic meteorites can be reproduced for a small range of model parameters. However, the necessary range of parameters is inconsistent with observations of disks around young stars and other constraints on meteorite parent body formation. Thus, counter to previous work, it is concluded that the global scale evolution of the solar nebula is not the cause for the observed depletion trends. △ Less

Submitted 1 February, 2007; originally announced February 2007.

Comments: 31 pages, 9 Figures, 1 Table. Revised for Meteoritics & Planetary Science (sorry, they do not accept LaTeX)

Dust Coagulation and Settling in Layered Protoplanetary Disks

Authors: Fred J. Ciesla

Abstract: Previous models of dust growth in protoplanetary disks considered either uniformly laminar or turbulent disks. This Letter explores how dust growth occurs in a layered protoplanetary disk in which the magnetorotational instability generates turbulence only in the surface layers of a disk. Two cases are considered: a completely laminar dead zone and a dead zone in which turbulence is ``stirred up… ▽ More Previous models of dust growth in protoplanetary disks considered either uniformly laminar or turbulent disks. This Letter explores how dust growth occurs in a layered protoplanetary disk in which the magnetorotational instability generates turbulence only in the surface layers of a disk. Two cases are considered: a completely laminar dead zone and a dead zone in which turbulence is ``stirred up'' from the MRI acting above. It is found that dust is depleted from high altitudes in layered disks faster than in those cases of a uniformly laminar or turblent disks. This is a result of the accelerated growth of particles in the turbulent regions and their storage in the lower levels where they escape energetic collisions which would result in disruption. Thus the regions of a protoplanetary disk above a dead zone would become rapidly depleted in small dust grains, whereas the outer regions, where the MRI is active throughout, will maintain a small dust poplulation at all heights due to the disruptive collisions and vertical mixing from turbulence. This structure is similar to that which has been inferred for disks around TW Hydra, GM Auriga, and CoKu Tau/4, which are depleted in dust close to the star, but are optically thick at larger heliocentric distances. △ Less

Submitted 27 November, 2006; originally announced November 2006.

Comments: 4 pages, 3 figures, accepted to ApJ Letters

The Evolution of the Water Distribution in a Viscous Protoplanetary Disk

Authors: Fred J. Ciesla, Jeffrey N. Cuzzi

Abstract: (Abridged) Astronomical observations have shown that protoplanetary disks are dynamic objects through which mass is transported and accreted by the central star. Age dating of meteorite constituents shows that their creation, evolution, and accumulation occupied several Myr, and over this time disk properties would evolve significantly. Moreover, on this timescale, solid particles decouple from… ▽ More (Abridged) Astronomical observations have shown that protoplanetary disks are dynamic objects through which mass is transported and accreted by the central star. Age dating of meteorite constituents shows that their creation, evolution, and accumulation occupied several Myr, and over this time disk properties would evolve significantly. Moreover, on this timescale, solid particles decouple from the gas in the disk and their evolution follows a different path. Here we present a model which tracks how the distribution of water changes in an evolving disk as the water-bearing species experience condensation, accretion, transport, collisional destruction, and vaporization. Because solids are transported in a disk at different rates depending on their sizes, the motions will lead to water being concentrated in some regions of a disk and depleted in others. These enhancements and depletions are consistent with the conditions needed to explain some aspects of the chemistry of chondritic meteorites and formation of giant planets. The levels of concentration and depletion, as well as their locations, depend strongly on the combined effects of the gaseous disk evolution, the formation of rapidly migrating rubble, and the growth of immobile planetesimals. We present examples of evolution under a range of plausible assumptions and demonstrate how the chemical evolution of the inner region of a protoplanetary disk is intimately connected to the physical processes which occur in the outer regions. △ Less

Submitted 11 November, 2005; originally announced November 2005.

Comments: 45 pages, 7 figures, revised for publication in Icarus