Source author record

Fred J. Ciesla

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

16works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

16 published item(s)

preprint2022arXiv

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

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

preprint2022arXiv

Thermal Processing of Solids Encountering a Young Jovian Core

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.

preprint2020arXiv

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

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.

preprint2020arXiv

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

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

preprint2019arXiv

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

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

preprint2016arXiv

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

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.

preprint2016arXiv

Tracing water vapor and ice during dust growth

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.

preprint2015arXiv

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

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.

preprint2015arXiv

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

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

preprint2015arXiv

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

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

preprint2014arXiv

The Phases of Water Ice in the Solar Nebula

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.

preprint2012arXiv

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

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.

preprint2011arXiv

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

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.

preprint2011arXiv

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

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.

preprint2010arXiv

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

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.

preprint2007arXiv

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

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.