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.

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

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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 {\alpha}_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 ({\alpha}_visc < 10^-3) the snow line produces a ring-like structure with a minimum at {\alpha}_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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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-$\alpha$" 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.

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.

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.

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.

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.

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