Source author record

Hilke E. Schlichting

Hilke E. Schlichting appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

astro-ph.EP astro-ph.SR

Catalog footprint

What is connected

21works

2topics

4close collaborators

Actions

Connect this record

Open graph Browse works

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.

preprint2024arXiv

The Kepler Giant Planet Search. I: A Decade of Kepler Planet-host Radial Velocities from W. M. Keck Observatory

Despite the importance of Jupiter and Saturn to Earth's formation and habitability, there has not yet been a comprehensive observational study of how giant exoplanets correlate with the architectural properties of close-in, sub-Neptune-sized exoplanets. This is largely because transit surveys are particularly insensitive to planets at orbital separations >1 au, and so their census of Jupiter-like planets is incomplete, inhibiting our study of the relationship between Jupiter-like planets and the small planets that do transit. To investigate the relationship between close-in, small and distant, giant planets, we conducted the Kepler Giant Planet Survey (KGPS). Using the W. M. Keck Observatory High Resolution Echelle Spectrometer, we spent over a decade collecting 2844 radial velocities (RVs; 2167 of which are presented here for the first time) of 63 Sunlike stars that host 157 transiting planets. We had no prior knowledge of which systems would contain giant planets beyond 1 au, making this survey unbiased with respect to previously detected Jovians. We announce RV-detected companions to 20 stars from our sample. These include 13 Jovians (0.3 MJ < M sin i < 13 MJ, 1 au < a < 10 au), eight nontransiting sub-Saturns, and three stellar-mass companions. We also present updated masses and densities of 84 transiting planets. The KGPS project leverages one of the longest-running and most data-rich collections of RVs of the NASA Kepler systems yet, and it will provide a basis for addressing whether giant planets help or hinder the growth of sub-Neptune-sized and terrestrial planets. Future KGPS papers will examine the relationship between small, transiting planets and their long-period companions.

preprint2022arXiv

Chemical equilibrium between Cores, Mantles, and Atmospheres of Super-Earths and Sub-Neptunes, and Implications for their Compositions, Interiors and Evolution

We investigate equilibrium chemistry between molten metal and silicate, and a hydrogen-rich envelope using 18 independent reactions among 25 phase components for sub-Neptune-like exoplanets. Both reactive metal and unreactive metal sequestered in an isolated core are modeled. The overarching effects of equilibration are oxidation of the envelope and reduction of the mantle and core. Hydrogen and oxygen typically comprise significant fractions of metal cores at chemical equilibrium, leading to density deficits that offer a possible alternative explanation for the low densities of the Trappist-1 planets. Reactions with the magma ocean produce significant amounts of SiO and H$_2$O in the envelopes directly above the magma ocean. Molar concentrations in the envelopes of planets with reactive metal are H$_2>\rm{SiO}>CO\sim\rm{Na}\sim\rm{Mg}>\rm{H}_2\rm{O}>>\rm{CO}_2\sim\rm{CH}_4>>\rm{O}_2$ while for the unreactive metal case H$_2$O becomes the second most abundant species, after H$_2$, providing an arbiter for the two scenarios amenable to observation. The water abundances in the atmospheres exceed those in the mantles by at least an order of magnitude in both scenarios. The water concentrations in the silicate mantles are $\sim 0.01$ wt\% and $\sim 0.1$ wt\% in the reactive and unreactive metal core cases, respectively, limiting the H$_2$O that might be outgassed in a future super-Earth. Less dissolved water in the reactive core case is due to sequestration of H and O in the Fe-rich metal. The total hydrogen budget of most sub-Neptunes can be, to first order, estimated from their atmospheres alone, as the atmospheres typically contain more than 90\% of all H.

preprint2022arXiv

Properties of the radius valley around low mass stars: Predictions from the core-powered mass-loss mechanism

In recent years, analyzing the bimodality in the size distribution of small planets, i.e., the `radius valley', has given us unprecedented insight into the planet formation process. Here we explore the properties of the radius valley for low mass stars, assuming that the core-powered mass-loss is the dominant process shaping the small exoplanet population. We show that the slope of radius valley in the planet size-orbital period space, to first-order, does not vary with stellar mass and has a negative slope of $\text{d log}R_p/\text{d log}P \simeq -0.11$ even for stars as small as 0.1 $M_\odot$, as observed in latest studies. Furthermore, we find that the slope of the radius valley in the planet size-stellar mass space is $\text{d log}R_p/\text{d log}M_\ast \simeq (3 ζ- 2)/36$ where $ζ$ is given by the stellar mass-luminosity relation $L_\ast \propto M_\ast^ζ$. Because $ζ$ is $\gtrsim$ 2 and increases with stellar mass, we predict that the radius valley has a positive slope in the planet size-stellar mass space across FGKM dwarfs. This slope, however, decreases (increases) in magnitude towards lower (higher) mass stars, due to the variation of $ζ$ with stellar mass. While around 1.0 $M_\odot$ stars the slope is $\text{d log}R_p/\text{d log}M_\ast \sim 0.37$, it is as low as $\sim 0.13$ around 0.1 $M_\odot$ stars. In addition, we find that the radius valley is narrower and less empty around lower mass stars. Finally, we show that predictions for the radius valley for core-powered mass-loss and photoevaporation become increasingly distinct for lower mass stars.

preprint2022arXiv

The importance of silicate vapor in determining the structure, radii, and envelope mass fractions of sub-Neptunes

Substantial silicate vapor is expected to be in chemical equilibrium at temperature conditions typical of the silicate-atmosphere interface of sub-Neptune planets, which can exceed 5000 K. Previous models of the atmospheric structure and evolution of these exoplanets, which have been used to constrain their atmospheric mass fractions, have neglected this compositional coupling. In this work, we show that silicate vapor in a hydrogen-dominated atmosphere acts as a condensable species, decreasing in abundance with altitude. The resultant mean molecular weight gradient inhibits convection at temperatures above $\sim 4000$ K, inducing a near-surface radiative layer. This radiative layer decreases the planet's total radius compared to a planet with the same base temperature and a convective, pure H/He atmosphere. Therefore, we expect silicate vapor to have major effects on the inferred envelope mass fraction and thermal evolution of sub-Neptune planets. We demonstrate that differences in radii, and hence in inferred atmospheric masses, are largest for planets which have larger masses, equilibrium temperatures, and atmospheric mass fractions. The effects are largest for younger planets, but differences can persist on gigayear time-scales for some sub-Neptunes. For a $10 M_\oplus$ planet with $T_\mathrm{eq}=1000$ K and an age of $\sim 300$ Myr, an observed radius consistent with an atmospheric mass fraction of 10% when accounting for silicate vapor would be misinterpreted as indicating an atmospheric mass fraction of 2% if a H/He-only atmosphere were assumed. The presence of silicate vapor in the atmosphere is also expected to have important implications for the accretion and loss of primordial hydrogen atmospheres.

preprint2021arXiv

Atmosphere Loss by Aerial Bursts

We present a simple analytic description of atmospheric mass loss by aerial bursts and demonstrate that mass loss from aerial bursts becomes significant when the maximum impactor size that leads to an aerial burst rather than a ground explosion, $r_o$, is larger than the minimum impactor size needed to achieve atmospheric loss, $r_{min}$. For vertical trajectories, which give the most stringent limit, this condition is approximately satisfied when $ρ_o/ρ_i \gtrsim 0.4 v_e/v_\infty$, which implies atmospheric densities need to be comparable to impactor densities for impactor velocities that are a few times the escape velocity of the planet. The range of impactor radii resulting in aerial burst-induced mass loss, $r_o-r_{min}$, increases with the ratio of the atmosphere to the impactor density and with the trajectory angle of the impactor. The range of impactor radii that result in aerial burst-induced mass loss and the atmospheric mass lost is larger in adiabatic atmospheres than isothermal atmospheres of equivalent total mass, scale height, and atmospheric surface density. Our results imply that aerial bursts are not expected to significantly contribute to the atmospheric mass-loss history of Earth, but are expected to play an important role for planets and exoplanets similar to Neptune with significant atmospheres. For Neptune-like atmospheres, the atmospheric mass ejected per impactor mass by aerial bursts is comparable to that lost by ground explosions, which implies that, for impactors following a Dohnanyi size distribution, overall loss by aerial busts is expected to exceed that by ground explosions by a factor of $(r_{ground}/r_{aerial})^{0.5}$.

preprint2020arXiv

Signatures of the Core-Powered Mass-Loss Mechanism in the Exoplanet Population: Dependence on Stellar Properties and Observational Predictions

Recent studies have shown that atmospheric mass-loss powered by the cooling luminosity of a planet's core can explain the observed radius valley separating super-Earths and sub-Neptunes, even without photoevaporation. In this work, we investigate the dependence of this core-powered mass-loss mechanism on stellar mass ($M_\ast$), metallicity ($Z_\ast$) and age ($τ_\ast$). Without making any changes to the underlying planet population, we find that the core-powered mass-loss model yields a shift in the radius valley to larger planet sizes around more massive stars with a slope given by $\text{d log}R_p/\text{d log}M_\ast \simeq 0.35$, in agreement with observations. To first order, this slope is driven by the dependence of core-powered mass-loss on the bolometric luminosity of the host star and is given by $\text{d log}R_p/\text{d log}M_\ast \simeq (3α-2)/36 \simeq 0.33$, where $(L_\ast/L_\odot) = (M_\ast/M_\odot)^α$ is the stellar mass-luminosity relation and $α\simeq 4.6$ for the CKS dataset. We therefore find, in contrast to photoevaporation models, no evidence for a linear correlation between planet and stellar mass, but can't rule it out either. In addition, we show that the location of the radius valley is, to first order, independent of stellar age and metallicity. Since core-powered mass-loss proceeds over Gyr timescales, the abundance of super-Earths relative to sub-Neptunes increases with age but decreases with stellar metallicity. Finally, due the dependence of the envelope's cooling timescale on metallicity, we find that the radii of sub-Neptunes increase with metallicity and decrease with age with slopes given by $\text{d log}R_p/\text{d log}Z_\ast \simeq 0.1$ and $\text{d log}R_p/\text{d log}τ_\ast \simeq -0.1$, respectively. We conclude with a series of observational tests that can differentiate between core-powered mass-loss and photoevaporation models.

preprint2020arXiv

Where are the Extrasolar Mercuries?

We utilize observations of 16 white dwarf stars to calculate and analyze the oxidation states of the parent bodies accreting onto the stars. Oxygen fugacity, a measure of overall oxidation state for rocks, is as important as pressure and temperature in determining the structure of a planet. We find that most of the extrasolar rocky bodies formed under oxidizing conditions, but approximately 1/4 of the polluted white dwarfs have compositions consistent with more reduced parent bodies. The difficulty in constraining the oxidation states of relatively reduced bodies is discussed and a model for the time-dependent evolution of the apparent oxygen fugacity for a hypothetical reduced body engulfed by a WD is investigated. Differences in diffusive fluxes of various elements through the WD envelope yield spurious inferred bulk elemental compositions and oxidation states of the accreting parent bodies under certain conditions. The worst case for biasing against detection of reduced bodies occurs for high effective temperatures. For moderate and low effective temperatures, evidence for relatively reduced parent bodies is preserved under most circumstances for at least several characteristic lifetimes of the debris disk.

preprint2016arXiv

Super-Earth Atmospheres: Self-Consistent Gas Accretion and Retention

Some recently discovered short-period Earth to Neptune sized exoplanets (super Earths) have low observed mean densities which can only be explained by voluminous gaseous atmospheres. Here, we study the conditions allowing the accretion and retention of such atmospheres. We self-consistently couple the nebular gas accretion onto rocky cores and the subsequent evolution of gas envelopes following the dispersal of the protoplanetary disk. Specifically, we address mass-loss due to both photo-evaporation and cooling of the planet. We find that planets shed their outer layers (dozens of percents in mass) following the disk's dispersal (even without photo-evaporation), and their atmospheres shrink in a few Myr to a thickness comparable to the radius of the underlying rocky core. At this stage, atmospheres containing less particles than the core (equivalently, lighter than a few % of the planet's mass) can be blown away by heat coming from the cooling core, while heavier atmospheres cool and contract on a timescale of Gyr at most. By relating the mass-loss timescale to the accretion time, we analytically identify a Goldilocks region in the mass-temperature plane in which low-density super Earths can be found: planets have to be massive and cold enough to accrete and retain their atmospheres, while not too massive or cold, such that they do not enter runaway accretion and become gas giants (Jupiters). We compare our results to the observed super-Earth population and find that low-density planets are indeed concentrated in the theoretically allowed region. Our analytical and intuitive model can be used to investigate possible super-Earth formation scenarios.

preprint2016arXiv

The Formation of Super-Earths and Mini-Neptunes with Giant Impacts

The majority of discovered exoplanetary systems harbour a new class of planets, bodies typically several times more massive than Earth but orbiting their host stars well inside the orbit of Mercury. The origin of these close-in super-Earths and mini-Neptunes is a major unanswered question in planet formation. Unlike Earth, whose atmosphere contains $<10^{-6}$ its total mass, a large fraction of close-in planets have significant gaseous envelopes, containing $1 -10\%$ or more of their total mass. It has been proposed that these close-in planets formed in situ either by delivery of $50-100M_{\oplus}$ of rocky material to the inner disc, or in a disc enhanced relative to the MMSN. In both cases, final assembly of the planets occurs by giant impacts (GIs). Here we test the viability of these scenarios. We show that atmospheres accreted by isolation masses are small ($10^{-3}-10^{-2}$ the core mass) and that atmospheric mass-loss during GIs is significant, with typical post-GI atmospheres that are $8 \times 10^{-4}$ the core mass. Such values are consistent with terrestrial planet atmospheres but more than an order of magnitude below atmospheric masses of $1-10\%$ inferred for many close-in exoplanets. In the most optimistic scenario with no core luminosity, post-GI envelope accretion from a depleted gas disc yields atmospheric masses that are several per cent the core mass. If the gravitational potential energy due to the last mass doubling of the planet by GIs is released over the disc dissipation time-scale as core luminosity, then envelope masses are reduced by about an order of magnitude. Finally we show that radial drift time-scales due to gas drag for many isolation masses are shorter than typical disc lifetimes. Given these challenges, we conclude that most observed close-in planets with envelopes larger than several per cent likely formed at larger separations from their host stars.

preprint2015arXiv

A Search for Ringed Exoplanets using Kepler Photometry

Models are developed to simulate lightcurves of stars dimmed by transiting exoplanets with and without rings. These models are then applied to \textit{Kepler} photometry to search for planetary rings in a sample of 21 exoplanets, mostly hot Jupiters, chosen to offer the best observational opportunity for discovering potential rings. We also examine what kinds of rings might be expected for these planets, in terms of both size and orientation, based on arguments involving the host planet's equilibrium temperature, its likely obliquities, and the formation and stability of possible ring systems. Finding no evidence for rings, for each of the 21 studied planets it is determined on an observational basis which potential rings can be rejected out of a representative set of fiducial rings, varying in both size and orientation. For 12 of the 21 planets, we determined that Saturn-like rings could be ruled out for at least certain orientations. Additionally, the detectability of rings is studied, and it is found that ringed planets with small obliquities (roughly $5^{\circ}-10^{\circ}$) can yield large signals, which is encouraging for future work, since such small obliquities are expected for hot Jupiters.

preprint2015arXiv

Stealing the Gas: Giant Impacts and the Large Diversity in Exoplanet Densities

Although current sensitivity limits are such that true Solar System analogs remain challenging to detect, numerous planetary systems have been discovered that are very different from our own Solar System. The majority of systems harbor a new class of planets, bodies that are typically several times more massive than the Earth but that orbit their host stars well inside the orbit of Mercury. These planets frequently show evidence for large Hydrogen and Helium envelopes containing several percent of the planet's mass and display a large diversity in mean densities. Here we show that this wide range can be achieved by one or two late giant impacts, which are frequently needed to achieve long-term orbital stability in multiple planet systems once the gas disk has disappeared. We demonstrate using hydrodynamical simulations that a single collision between similarly sized exoplanets can easily reduce the envelope-to-core-mass ratio by a factor of two and show that this leads to a corresponding increase in the observed mean density by factors of 2-3. In addition we investigate how envelope-mass-loss depends on envelope mass, planet radius, semi-major axis, and the mass distribution inside the envelope. We propose that a small number of giant impacts may be responsible for the large observed spread in mean densities, especially for multiple-planet systems containing planets with very different densities and which have not been significantly sculpted by photo evaporation.

preprint2014arXiv

Formation of close in Super-Earths \& Mini-Neptunes: Required Disk Masses \& Their Implications

Recent observations by the {\it Kepler} space telescope have led to the discovery of more than 4000 exoplanet candidates consisting of many systems with Earth- to Neptune-sized objects that reside well inside the orbit of Mercury, around their respective host stars. How and where these close-in planets formed is one of the major unanswered questions in planet formation. Here we calculate the required disk masses for {\it in situ} formation of the {\it Kepler} planets. We find that, if close-in planets formed as {\it isolation masses}, then standard gas-to-dust ratios yield corresponding gas disks that are gravitationally unstable for a significant fraction of systems, ruling out such a scenario. We show that the maximum width of a planet's accretion region in the absence of any migration is $2 v_{esc}/Ω$, where $v_{esc}$ is the escape velocity of the planet and $Ω$ the Keplerian frequency and use it to calculate the required disk masses for {\it in situ} formation with giant impacts. Even with giant impacts, formation without migration requires disk surface densities in solids at semi-major axes less than 0.1~AU of $10^3-10^5 \rm{~g~cm^{-2}}$ implying typical enhancements above the minimum-mass solar nebular (MMSN) by at least a factor of 20. Corresponding gas disks are below, but not far from, the gravitational stability limit. In contrast, formation beyond a few AU is consistent with MMSN disk masses. This suggests that migration of either solids or fully assembled planets is likely to have played a major role in the formation of close-in super-Earths and mini-Neptunes.

preprint2013arXiv

Dynamical and collisional constraints on a stochastic late veneer on the terrestrial planets

Given their tendency to be incorporated into the core during differentiation, the highly-siderophile elements (HSEs) in Earth's mantle are thought to have been accreted as a `late veneer' after the end of the giant impact phase. Bottke et al (2010) proposed that the large Earth-to-Moon HSE abundance ratio can be explained if the late veneer was characterized by large (D = 1000-4000km) impactors. Here we simulate the evolution of the terrestrial planets during a stochastic late veneer phase from the end of accretion until the start of the late heavy bombardment ~500 Myr later. We show that a late veneer population of 0.05 Earth masses dominated by large (D > 1000km) bodies naturally delivers a ~0.01 Earth mass veneer to Earth, consistent with constraints. The eccentricities and inclinations of the terrestrial planets are excited by close encounters with the largest late veneer bodies. We find the best agreement with their post-veneer orbits if either a) the terrestrial planets' pre-veneer angular momentum deficit AMD_0 was less than half of the current one AMD_now, or b) AMD_0 <= AMD_now and the veneer was limited to smaller (D_max <= 2000km) bodies. Veneer impacts on Venus, Earth and Mars were mostly accretionary but on Mercury and the Moon they were mostly erosive. In ~20% of simulations an energetic impact occurred that could have removed >25% of Mercury's mass, thereby increasing its iron mass fraction. We show that, due to the erosive nature of larger impacts, the Moon cannot accrete any material from objects larger than 500-1000km. The large Earth-to-Moon HSE abundance ratio is explained if the late veneer included large impactors (D >= 500-1000km) regardless of their size distribution, as long as most of Earth's veneer came from large bodies. The spin angular momenta imparted by stochastic late veneer impacts was far in excess of the current ones for Mercury and Venus.

preprint2013arXiv

Initial Planetesimal Sizes and the Size Distribution of Small Kuiper Belt Objects

We show, by comparing observations with theoretical models, that the observed Kuiper Belt size distribution is well matched by coagulation models, which start from an initial planetesimal population with radii of about 1km, and subsequent collisional evolution. We find that the observed size distribution for R > 30km has not been modified by collisional evolution over the age of the solar system, and that the size distribution below R ~ 30km has been modified by collisions and that its slope is well matched by collisional evolution models that use published strength laws. We investigate in detail the resulting size distribution of bodies ranging from 0.01km to 30km and find that its slope changes several times as a function of radius before approaching the expected value for an equilibrium collisional cascade of material strength dominated bodies for R < 0.1km. Compared to a single power law size distribution that would span the whole range from 0.01km to 30km, we find a strong deficit of bodies around R ~10km and a strong excess of bodies around 2km. This deficit and excess are caused by an excess mass in small planetesimals in the km size range that was left over from the runaway growth phase and that leaves a signature in the size distribution that is not erased after 4.5 Gyrs of collisional evolution. Observations of the small KBO size distribution can therefore test if large KBOs grew as a result of runaway growth and constrain the initial planetesimal sizes. We find that results from recent KBO occultation surveys and the observed KBO size distribution can be best matched by an initial planetesimal population that contained about equal mass per logarithmic mass bin in bodies ranging from 0.4km to 4km in radius. In addition, we find that we cannot match the observed KBO size distribution if most of the planetesimal mass was contained in bodies that were 10km in radius or larger.

preprint2013arXiv

Overstable Librations can account for the Paucity of Mean Motion Resonances among Exoplanet Pairs

We assess the multi-planet systems discovered by the Kepler satellite in terms of current ideas about orbital migration and eccentricity damping due to planet-disk interactions. Our primary focus is on mean motion resonances. Only a few percent of planet pairs are in close proximity to a resonance. However, predicted migration rates (parameterized by $τ_n=n/{|\dot n|}$) imply that during convergent migration most planets would have been captured into first order resonances. Eccentricity damping (parameterized by $τ_e=e/{|\dot e|}$) offers a plausible resolution. Estimates suggest $τ_e/τ_n\sim (h/a)^2\sim 10^{-2}$, where $h/a$ is the ratio of disk thickness to radius. Together, eccentricity damping and orbital migration give rise to an equilibrium eccentricity, $e_{eq}\sim(τ_e/τ_n)^{1/2}$. Capture is permanent provided $e_{eq}\lesssim μ^{1/3}$, where $μ$ denotes the planet to star mass ratio. But for $e_{eq}\gtrsim μ^{1/3}$, capture is only temporary because librations around equilibrium are overstable and lead to passage through resonance on timescale $τ_e$. Most Kepler planet pairs have $e_{eq}>μ^{1/3}$. Since $τ_n>> τ_e$ is the timescale for migration between neighboring resonances, only a modest percentage of pairs end up trapped in resonances after the disk disappears. Planet pairs close to a mean motion resonance typically exhibit period ratios 1-2% larger than those for exact resonance. The direction of this shift undoubtedly reflects the same asymmetry that requires convergent migration for resonance capture. Permanent resonance capture at these separations from exact resonance would demand $μ(τ_n/τ_e)^{1/2}\gtrsim 0.01$, a value that estimates of $μ$ from transit data and $(τ_e/τ_n)^{1/2}$ from theory are insufficient to match. Plausible alternatives involve eccentricity damping during or after disk dispersal. (Abridged)

preprint2012arXiv

Measuring the Abundance of sub-kilometer sized Kuiper Belt Objects using Stellar Occultations

We present here the analysis of about 19,500 new star hours of low ecliptic latitude observations (|b| < 20 deg) obtained by the Hubble Space Telescope's FGS over a time span of more than nine years; which is an addition to the 12,000 star hours previously analyzed by Schlichting et al. (2009). Our search for stellar occultations by small Kuiper belt objects (KBOs) yielded one new candidate event corresponding to a body with a 530 +/-70m radius at a distance of about 40AU. Using bootstrap simulations, we estimate a probability of approx 5%, that this event is due to random statistical fluctuations within the new data set. Combining this new event with the single KBO occultation reported by Schlichting et al. (2009) we arrive at the following results: 1) The ecliptic latitudes of 6.6 deg and 14.4 deg of the two events are consistent with the observed inclination distribution of larger, 100km-sized KBOs. 2) Assuming that small, sub-km sized KBOs have the same ecliptic latitude distribution as their larger counterparts, we find an ecliptic surface density of KBOs with radii larger than 250m of N(r>250m) = 1.1^{+1.5}_{-0.7} x 10^7 deg^{-2}; if sub-km sized KBOs have instead a uniform ecliptic latitude distribution for -20 deg < b< 20 deg then N(r>250m) = 4.4^{+5.8}_{-2.8} x 10^6 deg^{-2}. This is the best measurement of the surface density of sub-km sized KBOs to date. 3) Assuming the KBO size distribution can be well described by a single power law given by N(>r) \propto r^{1-q}, where N(>r) is the number of KBOs with radii greater than r, and q is the power law index, we find q=3.8+/-0.2 for a KBO ecliptic latitude distribution that follows the observed distribution for larger, 100-km sized KBOs. 4) Regardless of the exact power law, our results suggest that small KBOs are numerous enough to satisfy the required supply rate for the Jupiter family comets. (Abridged)

preprint2012arXiv

The Last Stages of Terrestrial Planet Formation: Dynamical Friction and the Late Veneer

The final stage of terrestrial planet formation consists of the cleanup of residual planetesimals after the giant impact phase. Dynamically, a residual planetesimal population is needed to damp the high eccentricities of the terrestrial planets after the giant impact stage. Geochemically, highly siderophile element (HSE) abundance patterns inferred for the terrestrial planets and the Moon suggest that a total of about 0.01 M_Earth of chondritic material was delivered as `late veneer' by planetesimals to the terrestrial planets after the end of giant impacts. Here we combine these two independent lines of evidence for a leftover population of planetesimals and show that: 1) A residual planetesimal population containing 0.01 M_Earth is able to damp the eccentricities of the terrestrial planets after giant impacts to their observed values. 2) At the same time, this planetesimal population can account for the observed relative amounts of late veneer added to the Earth, Moon and Mars provided that the majority of the late veneer was delivered by small planetesimals with radii <10m. These small planetesimal sizes are required to ensure efficient damping of the planetesimal's velocity dispersion by mutual collisions, which in turn ensures that the planets' accretion cross sections are significantly enhanced by gravitational focusing above their geometric values. Specifically we find, in the limit that the relative velocity between the terrestrial planets and the planetesimals is significantly less than the terrestrial planets' escape velocities, that gravitational focusing yields an accretion ratio Earth/Mars~17, which agrees well with the accretion ratio inferred from HSEs of 12-23. For the Earth-Moon system, we find an accretion ratio of ~200, which is consistent with estimates of 150-700 derived from HSE abundances that include the lunar crust as well as mantle component. (Abridged)

preprint2011arXiv

Self-consistent size and velocity distributions of collisional cascades

The standard theoretical treatment of collisional cascades derives a steady-state size distribution assuming a single constant velocity dispersion for all bodies regardless of size. Here we relax this assumption and solve self-consistently for the bodies' steady-state size and size-dependent velocity distributions. Specifically, we account for viscous stirring, dynamical friction, and collisional damping of the bodies' random velocities in addition to the mass conservation requirement typically applied to find the size distribution in a steady-state cascade. The resulting size distributions are significantly steeper than those derived without velocity evolution. For example, accounting self-consistently for the velocities can change the standard q=3.5 power-law index of the Dohnanyi (1969) differential size spectrum to an index as large as q=4. Similarly, for bodies held together by their own gravity, the corresponding power-law index range 2.88<q<3.14 of Pan & Sari (2005) can steepen to values as large as q=3.26. Our velocity results allow quantitative predictions of the bodies' scale heights as a function of size. Together with our predictions, observations of the scale heights for different sized bodies for the Kuiper belt, the asteroid belt, and extrasolar debris disks may constrain the mass and number of large bodies stirring the cascade as well as the colliding bodies' internal strengths.

preprint2011arXiv

Using Kuiper Belt Binaries to Constrain Neptune's Migration History

Approximately 10-20% of all Kuiper belt objects (KBOs) occupy mean-motion resonances with Neptune. This dynamical configuration likely resulted from resonance capture as Neptune migrated outward during the late stages of planet formation. The details of Neptune's planetesimal-driven migration, including its radial extent and the concurrent eccentricity evolution of the planet, are the subject of considerable debate. Two qualitatively different proposals for resonance capture have been proposed--migration-induced capture driven by smooth outward evolution of Neptune's orbit and chaotic capture driven by damping of the planet's eccentricity near its current semi-major axis. We demonstrate that the distribution of comparable-mass, wide-separation binaries occupying resonant orbits can differentiate between these two scenarios. If migration-induced capture occurred, this fraction records information about the formation locations of different populations of KBOs. Chaotic capture, in contrast, randomizes the orbits of bodies as they are placed in resonance. In particular, migration-induced capture produces the following signatures. The 2:1 resonance should contain a dynamically cold component, with inclinations less than 5-10 degrees, having a binary fraction comparable to that among cold classical KBOs. If the 3:2 resonance also hosts a cold component, its binary fraction should be 20-30% lower than in the cold classical belt. Among cold 2:1 (and if present 3:2) KBOs, objects with eccentricities e < 0.2 should have a binary fraction ~20% larger than those with e > 0.2. Searches for cold components in the binary fractions of resonant KBOs are currently practical. The additional migration-generated trends described here may be distinguished with objects discovered by LSST. (Abstract abridged.)

preprint2011arXiv

Warm Saturns: On the Nature of Rings around Extrasolar Planets that Reside Inside the Ice Line

We discuss the nature of rings that may exist around extrasolar planets. Taking the general properties of rings around the gas giants in the Solar System, we infer the likely properties of rings around exoplanets that reside inside the ice line. Due to their proximity to their host star, rings around such exoplanets must primarily consist of rocky materials. However, we find that despite the higher densities of rock compared to ice, most of the observed extrasolar planets with reliable radii measurements have sufficiently large Roche radii to support rings. For the currently known transiting extrasolar planets, Poynting-Robertson drag is not effective in significantly altering the dynamics of individual ring particles over a time span of $10^8$ years provided that they exceed about 1 m in size. In addition, we show that significantly smaller ring particles can exist in optically thick rings, for which we find typical ring lifetimes ranging from a few times $10^6$ to a few times $10^9$ years. Most interestingly, we find that many of the rings could have nontrivial Laplacian planes due to the increased effects of the orbital quadrupole caused by the exoplanets' proximity to their host star, allowing a constraint on the $J_2$ of extrasolar planets from ring observations. This is particular exciting, since a planet's $J_2$ reveals information about its interior structure. Furthermore, measurements of an exoplanet's oblateness and of its $J_2$, from warped rings, would together place limits on its spin period. Based on the constraints that we have derived for extrasolar rings, we anticipate that the best candidates for ring detections will come from transit observations by the Kepler spacecraft of extrasolar planets with semi-major axes $\sim 0.1$ AU and larger.

preprint2010arXiv

Runaway Growth During Planet Formation: Explaining the Size Distribution of Large Kuiper Belt Objects

Runway growth is an important stage in planet formation during which large protoplanets form, while most of the initial mass remains in small planetesimals. The amount of mass converted into large protoplanets and their resulting size distribution are not well understood. Here, we use analytic work, that we confirm by coagulation simulations, to describe runaway growth and the corresponding evolution of the velocity dispersion. We find that runaway growth proceeds as follows: Initially all the mass resides in small planetesimals, with mass surface density σ, and large protoplanets start to form by accreting small planetesimals. This growth continues until growth by merging large protoplanets becomes comparable to growth by planetesimal accretion. This condition sets in when Σ/σ~α^{3/4} ~ 10^{-3}, where Σis the mass surface density in protoplanets in a given logarithmic mass interval and αis the ratio of the size of a body to its Hill radius. From then on, protoplanetary growth and the evolution of the velocity dispersion become self-similar and Σremains roughly constant, since an increase in Σby accretion of small planetesimals is balanced by a decrease due to merging with large protoplanets. We show that this growth leads to a protoplanet size distribution given by N(>R) \propto R^{-3} where N(>R) is the number of objects with radii greater than R (i.e., a differential power-law index of 4). We apply our results to the Kuiper Belt, which is a relic of runaway growth. Our results successfully match the observed Kuiper belt size distribution, they illuminate the physical processes that shaped it and explain the total mass that is present in large Kuiper belt objects (KBOs) today. This work suggests that the current mass in large KBOs is primordial and that it has not been significantly depleted. Abridged

Hilke E. Schlichting

What is connected

Connect this record

See the researcher in context

Building this map preview

21 published item(s)

The Kepler Giant Planet Search. I: A Decade of Kepler Planet-host Radial Velocities from W. M. Keck Observatory

Chemical equilibrium between Cores, Mantles, and Atmospheres of Super-Earths and Sub-Neptunes, and Implications for their Compositions, Interiors and Evolution

Properties of the radius valley around low mass stars: Predictions from the core-powered mass-loss mechanism

The importance of silicate vapor in determining the structure, radii, and envelope mass fractions of sub-Neptunes

Atmosphere Loss by Aerial Bursts

Signatures of the Core-Powered Mass-Loss Mechanism in the Exoplanet Population: Dependence on Stellar Properties and Observational Predictions

Where are the Extrasolar Mercuries?

Super-Earth Atmospheres: Self-Consistent Gas Accretion and Retention

The Formation of Super-Earths and Mini-Neptunes with Giant Impacts

A Search for Ringed Exoplanets using Kepler Photometry

Stealing the Gas: Giant Impacts and the Large Diversity in Exoplanet Densities

Formation of close in Super-Earths \& Mini-Neptunes: Required Disk Masses \& Their Implications

Dynamical and collisional constraints on a stochastic late veneer on the terrestrial planets

Initial Planetesimal Sizes and the Size Distribution of Small Kuiper Belt Objects

Overstable Librations can account for the Paucity of Mean Motion Resonances among Exoplanet Pairs

Measuring the Abundance of sub-kilometer sized Kuiper Belt Objects using Stellar Occultations

The Last Stages of Terrestrial Planet Formation: Dynamical Friction and the Late Veneer

Self-consistent size and velocity distributions of collisional cascades

Using Kuiper Belt Binaries to Constrain Neptune's Migration History

Warm Saturns: On the Nature of Rings around Extrasolar Planets that Reside Inside the Ice Line

Runaway Growth During Planet Formation: Explaining the Size Distribution of Large Kuiper Belt Objects