Showing 1–50 of 60 results for author: Schlichting, H E

The miscibility of hydrogen and water in planetary atmospheres and interiors

Authors: Akash Gupta, Lars Stixrude, Hilke E. Schlichting

Abstract: Many planets in the solar system and across the galaxy have hydrogen-rich atmospheres overlying more heavy element-rich interiors with which they interact for billions of years. Atmosphere-interior interactions are thus crucial to understanding the formation and evolution of these bodies. However, this understanding is still lacking in part because the relevant pressure-temperature conditions are… ▽ More Many planets in the solar system and across the galaxy have hydrogen-rich atmospheres overlying more heavy element-rich interiors with which they interact for billions of years. Atmosphere-interior interactions are thus crucial to understanding the formation and evolution of these bodies. However, this understanding is still lacking in part because the relevant pressure-temperature conditions are extreme. We conduct molecular dynamics simulations based on Density Functional Theory to investigate how hydrogen and water interact over a wide range of pressure and temperature, encompassing the interiors of Neptune-sized and smaller planets. We determine the critical curve at which a single homogeneous phase exsolves into two separate, hydrogen-rich and water-rich phases, finding good agreement with existing experimental data. We find that the temperature along the critical curve increases with increasing pressure and shows the influence of a change in fluid structure from molecular to atomic near 30 GPa and 3000 K, which may impact magnetic field generation. The internal temperatures of many exoplanets, including TOI-270 d and K2-18 b may lie entirely above the critical curve: the envelope is expected to consist of a single homogeneous hydrogen-water fluid, that is much less susceptible to atmospheric loss as compared with a pure hydrogen envelope. As planets cool, they cross the critical curve, leading to rainout of water-rich fluid and an increase in internal luminosity. Compositions of the resulting outer, hydrogen-rich, and inner, water-rich envelopes depend on age and instellation and are governed by thermodynamics. Rainout of water may be occurring in Uranus and Neptune at present. △ Less

Submitted 5 July, 2024; originally announced July 2024.

Comments: 17 pages, 4 figures

Blowin' in the non-isothermal wind: core-powered mass loss with hydrodynamic radiative transfer

Authors: William Misener, Matthäus Schulik, Hilke E. Schlichting, James E. Owen

Abstract: The mass loss rates of planets undergoing core-powered escape are usually modeled using an isothermal Parker-type wind at the equilibrium temperature, $T_\mathrm{eq}$. However, the upper atmospheres of sub-Neptunes may not be isothermal if there are significant differences between the opacity to incident visible and outgoing infrared radiation. We model bolometrically-driven escape using aiolos, a… ▽ More The mass loss rates of planets undergoing core-powered escape are usually modeled using an isothermal Parker-type wind at the equilibrium temperature, $T_\mathrm{eq}$. However, the upper atmospheres of sub-Neptunes may not be isothermal if there are significant differences between the opacity to incident visible and outgoing infrared radiation. We model bolometrically-driven escape using aiolos, a hydrodynamic radiative-transfer code that incorporates double-gray opacities, to investigate the process's dependence on the visible-to-infrared opacity ratio, $γ$. For a value of $γ\approx 1$, we find that the resulting mass loss rates are well-approximated by a Parker-type wind with an isothermal temperature $T = T_\mathrm{eq}/2^{1/4}$. However, we show that over a range of physically plausible values of $γ$, the mass loss rates can vary by orders of magnitude, ranging from $10^{-5} \times$ the isothermal rate for low $γ$ to $10^5 \times$ the isothermal rate for high $γ$. The differences in mass loss rates are largest for small planet radii, while for large planet radii, mass loss rates become nearly independent of $γ$ and approach the isothermal approximation. We incorporate these opacity-dependent mass loss rates into a self-consistent planetary mass and energy evolution model and show that lower/higher $γ$ values lead to more/less hydrogen being retained after core-powered mass loss. In some cases, the choice of opacities determines whether or not a planet can retain a significant primordial hydrogen atmosphere. The dependence of escape rate on the opacity ratio may allow atmospheric escape observations to directly constrain a planet's opacities and therefore its atmospheric composition. △ Less

Submitted 24 May, 2024; originally announced May 2024.

Comments: 24 pages, 10 figures. Submitted to ApJ

JWST Reveals CH$_4$, CO$_2$, and H$_2$O in a Metal-rich Miscible Atmosphere on a Two-Earth-Radius Exoplanet

Authors: Björn Benneke, Pierre-Alexis Roy, Louis-Philippe Coulombe, Michael Radica, Caroline Piaulet, Eva-Maria Ahrer, Raymond Pierrehumbert, Joshua Krissansen-Totton, Hilke E. Schlichting, Renyu Hu, Jeehyun Yang, Duncan Christie, Daniel Thorngren, Edward D. Young, Stefan Pelletier, Heather A. Knutson, Yamila Miguel, Thomas M. Evans-Soma, Caroline Dorn, Anna Gagnebin, Jonathan J. Fortney, Thaddeus Komacek, Ryan MacDonald, Eshan Raul, Ryan Cloutier , et al. (6 additional authors not shown)

Abstract: Even though sub-Neptunes likely represent the most common outcome of planet formation, their natures remain poorly understood. In particular, planets near 1.5-2.5$\,R_\oplus$ often have bulk densities that can be explained equally well with widely different compositions and interior structures, resulting in grossly divergent implications for their formation. Here, we present the full 0.6-5.2… ▽ More Even though sub-Neptunes likely represent the most common outcome of planet formation, their natures remain poorly understood. In particular, planets near 1.5-2.5$\,R_\oplus$ often have bulk densities that can be explained equally well with widely different compositions and interior structures, resulting in grossly divergent implications for their formation. Here, we present the full 0.6-5.2$\,μ\mathrm{m}$ JWST NIRISS/SOSS+NIRSpec/G395H transmission spectrum of the 2.2$\,R_\oplus$ TOI-270d ($4.78\,M_\oplus$, $T_\mathrm{eq}$=350-380 K), delivering unprecedented sensitivity for atmospheric characterization in the sub-Neptune regime. We detect five vibrational bands of CH$_4$ at 1.15, 1.4, 1.7, 2.3, and 3.3$\,μ$m (9.4$σ$), the signature of CO$_2$ at 4.3$\,μ$m (4.8$σ$), water vapor (2.5$σ$), and potential signatures of SO$_2$ at 4.0$\,μ\mathrm{m}$ and CS$_2$ at 4.6$\,μ\mathrm{m}$. Intriguingly, we find an overall highly metal-rich atmosphere, with a mean molecular weight of $5.47_{-1.14}^{+1.25}$. We infer an atmospheric metal mass fraction of $58_{-12}^{+8}\%$ and a C/O of $0.47_{-0.19}^{+0.16}$, indicating that approximately half the mass of the outer envelope is in high-molecular-weight volatiles (H$_2$O, CH$_4$, CO, CO$_2$) rather than H$_2$/He. We introduce a sub-Neptune classification scheme and identify TOI-270d as a "miscible-envelope sub-Neptune" in which H$_2$/He is well-mixed with the high-molecular-weight volatiles in a miscible supercritical metal-rich envelope. For a fully miscible envelope, we conclude that TOI-270d's interior is $90_{-4}^{+3}\,$wt$\,\%$ rock/iron, indicating that it formed as a rocky planet that accreted a few wt % of H$_2$/He, with the overall envelope metal content explained by magma-ocean/envelope reactions without the need for significant ice accretion. TOI-270d may well be an archetype of the overall population of sub-Neptunes. △ Less

Submitted 5 March, 2024; originally announced March 2024.

Comments: 25 pages, 12 figures

Fleeting but not Forgotten: the Imprint of Escaping Hydrogen Atmospheres on Super-Earth Interiors

Authors: James G. Rogers, Hilke E. Schlichting, Edward D. Young

Abstract: Small, close-in exoplanets are divided into two sub-populations: super-Earths and sub-Neptunes. Most super-Earths are thought to have lost their primordially accreted hydrogen-dominated atmospheres via thermally driven winds. We consider the global chemical equilibrium of super-Earths and the lasting impacts of their fleeting hydrogen atmospheres. We find that hydrogen is efficiently sequestered i… ▽ More Small, close-in exoplanets are divided into two sub-populations: super-Earths and sub-Neptunes. Most super-Earths are thought to have lost their primordially accreted hydrogen-dominated atmospheres via thermally driven winds. We consider the global chemical equilibrium of super-Earths and the lasting impacts of their fleeting hydrogen atmospheres. We find that hydrogen is efficiently sequestered into the interior, oxidising iron and endogenously producing $\sim0.5-1.0\%$ water by mass. As the atmospheres of super-Earths are continuously sculpted by mass loss and chemical equilibration, they remain hydrogen-dominated by mole (number) fraction but become steam-dominated by mass, which may be observable with JWST for planets transitioning across the radius valley. One of the main effects of efficient sequestration of hydrogen into the interior is to produce an under-dense bulk interior compared to that of Earth. We predict bulk densities of super-Earths to be $\sim 5.0 \text{ g cm}^{-3}$ for a $1M_\oplus$ planet, which is consistent with high-precision mass measurements and also population-level inference analyses from atmospheric escape models. △ Less

Submitted 30 May, 2024; v1 submitted 21 February, 2024; originally announced February 2024.

Comments: 16 pages, 4 figures. Accepted for publication in ApJ. Comments welcome

A Tale of Two Peas-In-A-Pod: The Kepler-323 and Kepler-104 Systems

Authors: C. Alexander Thomas, Lauren M. Weiss, Howard Isaacson, Hilke E. Schlichting, Corey Beard, Casey L. Brinkman, Ashley Chontos, Paul Dalba, Fei Dai, Steven Giacalone, Jack Lubin, Judah Van Zandt, Malena Rice

Abstract: In order to understand the relationship between planet multiplicity, mass, and composition, we present newly measured masses of five planets in two planetary systems: Kepler-323 and Kepler-104. We used the HIRES instrument at the W.M. Keck Observatory to collect 79 new radial velocity measurements (RVs) for Kepler-323, which we combined with 48 literature RVs from TNG/HARPS-N. We also conducted a… ▽ More In order to understand the relationship between planet multiplicity, mass, and composition, we present newly measured masses of five planets in two planetary systems: Kepler-323 and Kepler-104. We used the HIRES instrument at the W.M. Keck Observatory to collect 79 new radial velocity measurements (RVs) for Kepler-323, which we combined with 48 literature RVs from TNG/HARPS-N. We also conducted a reanalysis of the Kepler-104 system, using 44 previously published RV measurements. Kepler-323 b and c have masses of $2.0^{+1.2}_{-1.1}$ M$_\oplus$ and 6.5$\pm1.6$ M$_\oplus$, respectively, whereas the three Kepler-104 planets are more massive (10.0$\pm2.8$ M$_\oplus$, $7.1^{+3.8}_{-3.5}$ M$_\oplus$, and $5.5^{+4.6}_{-3.5}$ M$_\oplus$ for planets b, c, and d, respectively). The Kepler-104 planets have densities consistent with rocky cores overlaid with gaseous envelopes ($4.1^{+1.2}_{-1.1}$ g/cc, $2.9^{+1.7}_{-1.5}$ g/cc, and $1.6^{+1.5}_{-1.1}$ g/cc respectively), whereas the Kepler-323 planets are consistent with having rocky compositions ($4.5^{+2.8}_{-2.4}$ g/cc and $9.9^{+2.7}_{-2.5}$ g/cc). The Kepler-104 system has among the lowest values for gap complexity ($\mathcal{C}$ = 0.004) and mass partitioning ($\mathcal{Q}$ = 0.03); whereas, the Kepler-323 planets have a mass partitioning similar to that of the Inner Solar System ($\mathcal{Q}$ = 0.28 and $\mathcal{Q}$ = 0.24, respectively). For both exoplanet systems, the uncertainty in the mass partitioning is affected equally by (1) individual mass errors of the planets and (2) the possible existence of undetected low-mass planets, meaning that both improved mass characterization and improved sensitivity to low-mass planets in these systems would better elucidate the mass distribution among the planets. △ Less

Submitted 20 February, 2024; originally announced February 2024.

Comments: Accepted for publication in Astronomical Journal

Shallower radius valley around low-mass hosts: Evidence for icy planets, collisions or high-energy radiation scatter

Authors: Cynthia S. K. Ho, James G. Rogers, Vincent Van Eylen, James E. Owen, Hilke E. Schlichting

Abstract: The radius valley, i.e., a dearth of planets with radii between 1.5 and 2 Earth radii, provides insights into planetary formation and evolution. Using homogenously revised planetary parameters from Kepler 1-minute short cadence light curves, we remodel transits of 72 small planets mostly orbiting low-mass stars, improving the precision and accuracy of planet parameters. By combining this sample wi… ▽ More The radius valley, i.e., a dearth of planets with radii between 1.5 and 2 Earth radii, provides insights into planetary formation and evolution. Using homogenously revised planetary parameters from Kepler 1-minute short cadence light curves, we remodel transits of 72 small planets mostly orbiting low-mass stars, improving the precision and accuracy of planet parameters. By combining this sample with a similar sample of planets around higher-mass stars, we determine the depth of the radius valley as a function of stellar mass. We find that the radius valley is shallower for low-mass stars compared to their higher mass counterparts. Upon comparison, we find that theoretical models of photoevaporation under-predict the number of planets observed inside the radius valley for low-mass stars: with decreasing stellar mass, the predicted fraction of planets inside the valley remains approximately constant whereas the observed fraction increases. We argue that this provides evidence for the presence of icy planets around low-mass stars. Alternatively, planets orbiting low-mass stars undergo more frequent collisions and scatter in the stars' high-energy output may also cause planets to fill the valley. We predict that more precise mass measurements for planets orbiting low mass stars would be able to distinguish between these scenarios. △ Less

Submitted 11 June, 2024; v1 submitted 22 January, 2024; originally announced January 2024.

Comments: 19 pages, accepted for publication in MNRAS

Under the light of a new star: evolution of planetary atmospheres through protoplanetary disc dispersal and boil-off

Authors: James G. Rogers, James E. Owen, Hilke E. Schlichting

Abstract: The atmospheres of small, close-in exoplanets are vulnerable to rapid mass-loss during protoplanetary disc dispersal via a process referred to as `boil-off', in which confining pressure from the local gas disc reduces, inducing atmospheric loss and contraction. We construct self-consistent models of planet evolution during gaseous core accretion and boil-off. As the surrounding disc gas dissipates… ▽ More The atmospheres of small, close-in exoplanets are vulnerable to rapid mass-loss during protoplanetary disc dispersal via a process referred to as `boil-off', in which confining pressure from the local gas disc reduces, inducing atmospheric loss and contraction. We construct self-consistent models of planet evolution during gaseous core accretion and boil-off. As the surrounding disc gas dissipates, we find that planets lose mass via subsonic breeze outflows which allow causal contact to exist between disc and planet. Planets initially accrete of order $\sim 10\%$ in atmospheric mass, however, boil-off can remove $\gtrsim 90\%$ of this mass during disc dispersal. We show that a planet's final atmospheric mass fraction is strongly dictated by the ratio of cooling timescale to disc dispersal timescale, as well as the planet's core mass and equilibrium temperature. With contributions from core cooling and radioactivity, we show that core luminosity eventually leads to the transition from boil-off to core-powered mass-loss. We find that smaller mass planets closest to their host star may have their atmospheres completely stripped through a combination of boil-off and core-powered mass-loss during disc dispersal, implying the existence of a population-level radius gap emerging as the disc disperses. We additionally consider the transition from boil-off/core-powered mass-loss to X-ray/EUV (XUV) photoevaporation by considering the penetration of stellar XUV photons below the planet's sonic surface. Finally, we show that planets may open gaps in their protoplanetary discs during the late stages of boil-off, which may enhance mass-loss rates. △ Less

Submitted 20 February, 2024; v1 submitted 20 November, 2023; originally announced November 2023.

Comments: 18 pages, 8 figures. Accepted for publication in MNRAS

CHIMERA Occultation Constraints on the Abundance of Kilometer-scale Kuiper Belt Objects

Authors: Qicheng Zhang, Gregg W. Hallinan, Navtej S. Saini, Hilke E. Schlichting, Leon K. Harding, Jennifer W. Milburn

Abstract: Occultations provide indirect sensitivity to the number density of small Kuiper Belt objects (KBOs) too faint to directly detect telescopically. We present results from the Caltech HI-speed Multicolor camERA (CHIMERA) survey with the Palomar Hale Telescope, which monitored stars over the central 5'x5' of the M22 globular cluster along the ecliptic plane for serendipitous occultations by kilometer-… ▽ More Occultations provide indirect sensitivity to the number density of small Kuiper Belt objects (KBOs) too faint to directly detect telescopically. We present results from the Caltech HI-speed Multicolor camERA (CHIMERA) survey with the Palomar Hale Telescope, which monitored stars over the central 5'x5' of the M22 globular cluster along the ecliptic plane for serendipitous occultations by kilometer-scale KBOs over 63 hr across 24 nights at a 33 Hz frame rate simultaneously in i' and g'. We adapted dense-field photometry and occultation template fitting techniques to this dataset, finding a 95% confidence upper limit on the occultation rate corresponding to an ecliptic sky density of <10^7 deg^-2 of >1 km diameter classical KBOs. We discuss a few of the occultation-like light curve signatures at the edge of the sensitivity limit responsible for setting the upper bounds, and their likely nonviability as true occultations. △ Less

Submitted 7 November, 2023; originally announced November 2023.

Comments: 13 pages, 7 figures, 1 table; AJ, in press

Journal ref: Astron.J. 166 (2023) 242

Investigating the Atmospheric Mass Loss of the Kepler-105 Planets Straddling the Radius Gap

Authors: Aaron Householder, Lauren M. Weiss, James E. Owen, Howard Isaacson, Andrew W. Howard, Daniel Fabrycky, Leslie A. Rogers, Hilke E. Schlichting, Benjamin J. Fulton, Erik A. Petigura, Steven Giacalone, Joseph M. Akana Murphy, Corey Beard, Ashley Chontos, Fei Dai, Judah Van Zandt, Jack Lubin, Malena Rice, Alex S. Polanski, Paul Dalba, Sarah Blunt, Emma V. Turtelboom, Ryan Rubenzahl, Casey Brinkman

Abstract: An intriguing pattern among exoplanets is the lack of detected planets between approximately $1.5$ R$_\oplus$ and $2.0$ R$_\oplus$. One proposed explanation for this "radius gap" is the photoevaporation of planetary atmospheres, a theory that can be tested by studying individual planetary systems. Kepler-105 is an ideal system for such testing due to the ordering and sizes of its planets. Kepler-1… ▽ More An intriguing pattern among exoplanets is the lack of detected planets between approximately $1.5$ R$_\oplus$ and $2.0$ R$_\oplus$. One proposed explanation for this "radius gap" is the photoevaporation of planetary atmospheres, a theory that can be tested by studying individual planetary systems. Kepler-105 is an ideal system for such testing due to the ordering and sizes of its planets. Kepler-105 is a sun-like star that hosts two planets straddling the radius gap in a rare architecture with the larger planet closer to the host star ($R_b = 2.53\pm0.07$ R$_\oplus$, $P_b = 5.41$ days, $R_c = 1.44\pm0.04$ R$_\oplus$, $P_c = 7.13$ days). If photoevaporation sculpted the atmospheres of these planets, then Kepler-105b would need to be much more massive than Kepler-105c to retain its atmosphere, given its closer proximity to the host star. To test this hypothesis, we simultaneously analyzed radial velocities (RVs) and transit timing variations (TTVs) of the Kepler-105 system, measuring disparate masses of $M_b = 10.8\pm2.3$ M$_\oplus$ ($ ρ_b = 0.97\pm0.22$ g cm$^{-3}$) and $M_c = 5.6\pm1.2$ M$_\oplus $ ($ρ_c = 2.64\pm0.61$ g cm$^{-3}$). Based on these masses, the difference in gas envelope content of the Kepler-105 planets could be entirely due to photoevaporation (in 76\% of scenarios), although other mechanisms like core-powered mass loss could have played a role for some planet albedos. △ Less

Submitted 5 December, 2023; v1 submitted 20 September, 2023; originally announced September 2023.

Comments: 14 pages, 3 figures, 2 tables, accepted for publication in The Astronomical Journal

The DECam Ecliptic Exploration Project (DEEP): V. The Absolute Magnitude Distribution of the Cold Classical Kuiper Belt

Authors: Kevin J. Napier, Hsing-Wen Lin, David W. Gerdes, Fred C. Adams, Anna M. Simpson, Matthew W. Porter, Katherine G. Weber, Larissa Markwardt, Gabriel Gowman, Hayden Smotherman, Pedro H. Bernardinelli, Mario Jurić, Andrew J. Connolly, J. Bryce Kalmbach, Stephen K. N. Portillo, David E. Trilling, Ryder Strauss, William J. Oldroyd, Chadwick A. Trujillo, Colin Orion Chandler, Matthew J. Holman, Hilke E. Schlichting, Andrew McNeill, the DEEP Collaboration

Abstract: The DECam Ecliptic Exploration Project (DEEP) is a deep survey of the trans-Neptunian solar system being carried out on the 4-meter Blanco telescope at Cerro Tololo Inter-American Observatory in Chile using the Dark Energy Camera (DECam). By using a shift-and-stack technique to achieve a mean limiting magnitude of $r \sim 26.2$, DEEP achieves an unprecedented combination of survey area and depth,… ▽ More The DECam Ecliptic Exploration Project (DEEP) is a deep survey of the trans-Neptunian solar system being carried out on the 4-meter Blanco telescope at Cerro Tololo Inter-American Observatory in Chile using the Dark Energy Camera (DECam). By using a shift-and-stack technique to achieve a mean limiting magnitude of $r \sim 26.2$, DEEP achieves an unprecedented combination of survey area and depth, enabling quantitative leaps forward in our understanding of the Kuiper Belt populations. This work reports results from an analysis of twenty 3 sq.\ deg.\ DECam fields along the invariable plane. We characterize the efficiency and false-positive rates for our moving-object detection pipeline, and use this information to construct a Bayesian signal probability for each detected source. This procedure allows us to treat all of our Kuiper Belt Object (KBO) detections statistically, simultaneously accounting for efficiency and false positives. We detect approximately 2300 candidate sources with KBO-like motion at S/N $>6.5$. We use a subset of these objects to compute the luminosity function of the Kuiper Belt as a whole, as well as the Cold Classical (CC) population. We also investigate the absolute magnitude ($H$) distribution of the CCs, and find consistency with both an exponentially tapered power-law, which is predicted by streaming instability models of planetesimal formation, and a rolling power law. Finally, we provide an updated mass estimate for the Cold Classical Kuiper Belt of $M_{CC}(H_r < 12) = 0.0017^{+0.0010}_{-0.0004} M_{\oplus}$, assuming albedo $p = 0.15$ and density $ρ= 1$ g cm$^{-3}$. △ Less

Submitted 18 September, 2023; originally announced September 2023.

Comments: Accepted by PSJ

The DECam Ecliptic Exploration Project (DEEP) IV: Constraints on the shape distribution of bright TNOs

Authors: R. Strauss, D. E. Trilling, P. H. Bernardinelli, C. Beach, W. J. Oldroyd, S. S. Sheppard, H. E. Schlichting, D. W. Gerdes, F. C. Adams, C. O. Chandler, C. Fuentes, M. J. Holman, M. Jurić, H. W. Lin, L. Markwardt, A. McNeill, M. Mommert, K. J. Napier, M. J. Payne, D. Ragozzine, A. S. Rivkin, H. Smotherman, C. A. Trujillo

Abstract: We present the methods and results from the discovery and photometric measurement of 26 bright (VR $>$ 24 trans-Neptunian objects (TNOs) during the first year (2019-20) of the DECam Ecliptic Exploration Project (DEEP). The DEEP survey is an observational TNO survey with wide sky coverage, high sensitivity, and a fast photometric cadence. We apply a computer vision technique known as a progressive… ▽ More We present the methods and results from the discovery and photometric measurement of 26 bright (VR $>$ 24 trans-Neptunian objects (TNOs) during the first year (2019-20) of the DECam Ecliptic Exploration Project (DEEP). The DEEP survey is an observational TNO survey with wide sky coverage, high sensitivity, and a fast photometric cadence. We apply a computer vision technique known as a progressive probabilistic Hough transform to identify linearly-moving transient sources within DEEP photometric catalogs. After subsequent visual vetting, we provide a photometric and astrometric catalog of our TNOs. By modeling the partial lightcurve amplitude distribution of the DEEP TNOs using Monte Carlo techniques, we find our data to be most consistent with an average TNO axis ratio b/a $<$ 0.5, implying a population dominated by non-spherical objects. Based on ellipsoidal gravitational stability arguments, we find our data to be consistent with a TNO population containing a high fraction of contact binaries or other extremely non-spherical objects. We also discuss our data as evidence that the expected binarity fraction of TNOs may be size-dependent. △ Less

Submitted 7 September, 2023; originally announced September 2023.

Mapping out the parameter space for photoevaporation and core-powered mass-loss

Authors: James E. Owen, Hilke E. Schlichting

Abstract: Understanding atmospheric escape in close-in exoplanets is critical to interpreting their evolution. We map out the parameter space over which photoevaporation and core-powered mass loss dominate atmospheric escape. Generally, the transition between the two regimes is determined by the location of the Bondi radius (i.e. the sonic point of core-powered outflow) relative to the penetration depth of… ▽ More Understanding atmospheric escape in close-in exoplanets is critical to interpreting their evolution. We map out the parameter space over which photoevaporation and core-powered mass loss dominate atmospheric escape. Generally, the transition between the two regimes is determined by the location of the Bondi radius (i.e. the sonic point of core-powered outflow) relative to the penetration depth of XUV photons. Photoevaporation dominates the loss when the XUV penetration depth lies inside the Bondi radius ($R_{XUV}<R_B$) and core-powered mass-loss when XUV radiation is absorbed higher up in the flow ($R_B<R_{XUV}$). The transition between the two regimes occurs at a roughly constant ratio of the planet's radius to its Bondi radius, with the exact value depending logarithmically on planetary and stellar properties. In general, core-powered mass-loss dominates for lower-gravity planets with higher equilibrium temperatures, and photoevaporation dominates for higher-gravity planets with lower equilibrium temperatures. However, planets can transition between these two mass-loss regimes during their evolution, and core-powered mass loss can ``enhance'' photo-evaporation over a significant region of parameter space. Interestingly, a planet that is ultimately stripped by core-powered mass-loss has likely only ever experienced core-powered mass-loss. In contrast a planet that is ultimately stripped by photoevaporation could have experienced an early phase of core-powered mass-loss. Applying our results to the observed super-Earth population suggests that it contains significant fractions of planets where each mechanism controlled the final removal of the H/He envelope, although photoevaporation appears to be responsible for the final carving of the exoplanet radius-valley. △ Less

Submitted 21 December, 2023; v1 submitted 31 July, 2023; originally announced August 2023.

Comments: Accepted for publication in MNRAS

Earth shaped by primordial H$_2$ atmospheres

Authors: Edward D. Young, Anat Shahar, Hilke E. Schlichting

Abstract: Earth's water, intrinsic oxidation state, and metal core density are fundamental chemical features of our planet. Studies of exoplanets provide a useful context for elucidating the source of these chemical traits. Planet formation and evolution models demonstrate that rocky exoplanets commonly formed with hydrogen-rich envelopes that were lost over time. These findings suggest that Earth may also… ▽ More Earth's water, intrinsic oxidation state, and metal core density are fundamental chemical features of our planet. Studies of exoplanets provide a useful context for elucidating the source of these chemical traits. Planet formation and evolution models demonstrate that rocky exoplanets commonly formed with hydrogen-rich envelopes that were lost over time. These findings suggest that Earth may also have formed from bodies with H$_2$-rich primary atmospheres. Here we use a self-consistent thermodynamic model to show that Earth's water, core density, and overall oxidation state can all be sourced to equilibrium between H$_2$-rich primary atmospheres and underlying magma oceans in its progenitor planetary embryos. Water is produced from dry starting materials resembling enstatite chondrites as oxygen from magma oceans reacts with hydrogen. Hydrogen derived from the atmosphere enters the magma ocean and eventually the metal core at equilibrium, causing metal density deficits matching that of Earth. Oxidation of the silicate rocks from solar-like to Earth-like oxygen fugacities also ensues as Si, along with H and O, alloys with Fe in the cores. Reaction with hydrogen atmospheres and metal-silicate equilibrium thus provides a simple explanation for fundamental features of Earth's geochemistry that is consistent with rocky planet formation across the galaxy. △ Less

Submitted 16 April, 2023; originally announced April 2023.

Comments: 3 main figures, 5 auxiliary figures

Journal ref: Nature, v. 616 (7956), 306-311 (2023)

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

Authors: Lauren M. Weiss, Howard Isaacson, Andrew W. Howard, Benjamin J. Fulton, Daniel Fabrycky, Daniel Jontof-Hutter, Jason H. Steffen, Hilke E. Schlichting, Corey Beard, Casey L. Brinkman, Ashley Chontos, Steven Giacalone, Michelle Hill, Molly R. Kosiarek, Mason G. MacDougall, Teo Močnik, Alex S. Polanski, Emma V. Turtelboom, Dakotah Tyler, Judah Van Zandt

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

Submitted 5 January, 2024; v1 submitted 31 March, 2023; originally announced April 2023.

Comments: Published in the Astrophysical Journal Supplement Series

Journal ref: ApJS 270 8 (2024)

Atmospheres as windows into sub-Neptune interiors: coupled chemistry and structure of hydrogen-silane-water envelopes

Authors: William Misener, Hilke E. Schlichting, Edward D. Young

Abstract: Sub-Neptune exoplanets are commonly hypothesized to consist of a silicate-rich magma ocean topped by a hydrogen-rich atmosphere. Previous work studying the outgassing of silicate material has demonstrated that such atmosphere-interior interactions can affect the atmosphere's overall structure and extent. But these models only considered SiO in an atmosphere of hydrogen gas, without considering che… ▽ More Sub-Neptune exoplanets are commonly hypothesized to consist of a silicate-rich magma ocean topped by a hydrogen-rich atmosphere. Previous work studying the outgassing of silicate material has demonstrated that such atmosphere-interior interactions can affect the atmosphere's overall structure and extent. But these models only considered SiO in an atmosphere of hydrogen gas, without considering chemical reactions between them. Here we couple calculations of the chemical equilibrium between H, Si, and O species with an atmospheric structure model. We find that substantial amounts of silane, SiH$_4$, and water, H$_2$O, are produced by the interaction between the silicate-rich interior and hydrogen-rich atmosphere. These species extend high into the atmosphere, though their abundance is greatest at the hottest, deepest regions. For example, for a 4 $M_\oplus$ planet with an equilibrium temperature of 1000 K, a base temperature of 5000 K, and a 0.1 $M_\oplus$ hydrogen envelope, silicon species and water can comprise 30 percent of the atmosphere by number at the bottom of the atmosphere. Due to this abundance enhancement, we find that convection is inhibited at temperatures $\gtrsim 2500$ K. This temperature is lower, implying that the resultant non-convective region is thicker, than was found in previous models which did not account for atmospheric chemistry. Our findings show that significant endogenous water is produced by magma-hydrogen interactions alone, without the need to accrete ice-rich material. We discuss the observability of the signatures of atmosphere-interior interaction and directions for future work, including condensate lofting and more complex chemical networks. △ Less

Submitted 21 June, 2023; v1 submitted 16 March, 2023; originally announced March 2023.

Comments: 13 pages, 12 figures, accepted for publication in MNRAS

Conclusive evidence for a population of water-worlds around M-dwarfs remains elusive

Authors: James G. Rogers, Hilke E. Schlichting, James E. Owen

Abstract: The population of small, close-in exoplanets is bifurcated into super-Earths and sub-Neptunes. We calculate physically motivated mass-radius relations for sub-Neptunes, with rocky cores and H/He dominated atmospheres, accounting for their thermal evolution, irradiation and mass-loss. For planets $\lesssim 10~$M$_\oplus$, we find that sub-Neptunes retain atmospheric mass fractions that scale with p… ▽ More The population of small, close-in exoplanets is bifurcated into super-Earths and sub-Neptunes. We calculate physically motivated mass-radius relations for sub-Neptunes, with rocky cores and H/He dominated atmospheres, accounting for their thermal evolution, irradiation and mass-loss. For planets $\lesssim 10~$M$_\oplus$, we find that sub-Neptunes retain atmospheric mass fractions that scale with planet mass and show that the resulting mass-radius relations are degenerate with results for `water-worlds' consisting of a 1:1 silicate-to-ice composition ratio. We further demonstrate that our derived mass-radius relation is in excellent agreement with the observed exoplanet population orbiting M-dwarfs and that planet mass and radii alone are insufficient to determine the composition of some sub-Neptunes. Finally, we highlight that current exoplanet demographics show an increase in the ratio of super-Earths to sub-Neptunes with both stellar mass (and therefore luminosity) and age, which are both indicative of thermally driven atmospheric escape processes. Therefore, such processes should not be ignored when making compositional inferences in the mass-radius diagram. △ Less

Submitted 3 April, 2023; v1 submitted 11 January, 2023; originally announced January 2023.

Comments: 16 pages, 2 figures. Accepted for publication in ApJL

The Exoplanet Radius Valley from Gas-driven Planet Migration and Breaking of Resonant Chains

Authors: Andre Izidoro, Hilke E. Schlichting, Andrea Isella, Rajdeep Dasgupta, Christian Zimmermann, Bertram Bitsch

Abstract: The size frequency distribution of exoplanet radii between 1 and 4$R_{\oplus}$ is bimodal with peaks at $\sim$1.4 $R_{\oplus}$ and $\sim$2.4 $R_{\oplus}$, and a valley at $\sim$1.8$R_{\oplus}$. This radius valley separates two classes of planets -- usually referred to as "super-Earths" and "mini-Neptunes" -- and its origin remains debated. One model proposes that super-Earths are the outcome of ph… ▽ More The size frequency distribution of exoplanet radii between 1 and 4$R_{\oplus}$ is bimodal with peaks at $\sim$1.4 $R_{\oplus}$ and $\sim$2.4 $R_{\oplus}$, and a valley at $\sim$1.8$R_{\oplus}$. This radius valley separates two classes of planets -- usually referred to as "super-Earths" and "mini-Neptunes" -- and its origin remains debated. One model proposes that super-Earths are the outcome of photo-evaporation or core-powered mass-loss stripping the primordial atmospheres of the mini-Neptunes. A contrasting model interprets the radius valley as a dichotomy in the bulk compositions, where super-Earths are rocky planets and mini-Neptunes are water-ice rich worlds. In this work, we test whether the migration model is consistent with the radius valley and how it distinguishes these views. In the migration model, planets migrate towards the disk inner edge forming a chain of planets locked in resonant configurations. After the gas disk dispersal, orbital instabilities "break the chains" and promote late collisions. This model broadly matches the period-ratio and planet-multiplicity distributions of Kepler planets, and accounts for resonant chains such as TRAPPIST-1, Kepler-223, and TOI-178. Here, by combining the outcome of planet formation simulations with compositional mass-radius relationships, and assuming complete loss of primordial H-rich atmospheres in late giant-impacts, we show that the migration model accounts for the exoplanet radius valley and the intra-system uniformity ("peas-in-a-pod") of Kepler planets. Our results suggest that planets with sizes of $\sim$1.4 $R_{\oplus}$ are mostly rocky, whereas those with sizes of $\sim$2.4 $R_{\oplus}$ are mostly water-ice rich worlds. Our results do not support an exclusively rocky composition for the cores of mini-Neptunes. △ Less

Submitted 11 October, 2022; originally announced October 2022.

Comments: Accepted for publicationn in ApJL

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

Authors: Akash Gupta, Lorraine Nicholson, Hilke E. Schlichting

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

Submitted 31 August, 2022; v1 submitted 27 May, 2022; originally announced May 2022.

Comments: 9 pages, 5 figures and 1 table. Accepted for publication in MNRAS

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

Authors: William Misener, Hilke E. Schlichting

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

Submitted 20 June, 2022; v1 submitted 12 January, 2022; originally announced January 2022.

Comments: 13 pages, 7 figures, accepted for publication in MNRAS after revision

Atmosphere Loss by Aerial Bursts

Authors: Isabella L. Trierweiler, Hilke E. Schlichting

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

Submitted 14 December, 2021; originally announced December 2021.

Comments: 10 pages, 9 figures

The fundamentals of Lyman-alpha exoplanet transits

Authors: James E. Owen, Ruth A. Murray-Clay, Ethan Schreyer, Hilke E. Schlichting, David Ardila, Akash Gupta, R. O. Parke Loyd, Evgenya L. Shkolnik, David K. Sing, Mark R. Swain

Abstract: Lyman-$α$ transits have been detected from several nearby exoplanets and are one of our best insights into the atmospheric escape process. However, due to ISM absorption, we typically only observe the transit signature in the blue-wing, making them challenging to interpret. This challenge has been recently highlighted by non-detections from planets thought to be undergoing vigorous escape. Pioneer… ▽ More Lyman-$α$ transits have been detected from several nearby exoplanets and are one of our best insights into the atmospheric escape process. However, due to ISM absorption, we typically only observe the transit signature in the blue-wing, making them challenging to interpret. This challenge has been recently highlighted by non-detections from planets thought to be undergoing vigorous escape. Pioneering 3D simulations have shown that escaping hydrogen is shaped into a cometary tail receding from the planet. Motivated by this work, we develop a simple model to interpret Lyman-$α$ transits. Using this framework, we show that the Lyman-$α$ transit depth is primarily controlled by the properties of the stellar tidal field rather than details of the escape process. Instead, the transit duration provides a direct measurement of the velocity of the planetary outflow. This result arises because the underlying physics is the distance a neutral hydrogen atom can travel before it is photoionized in the outflow. Thus, higher irradiation levels, expected to drive more powerful outflows, produce weaker, shorter Lyman-$α$ transits because the outflowing gas is ionized more quickly. Our framework suggests that the generation of energetic neutral atoms may dominate the transit signature early, but the acceleration of planetary material produces long tails. Thus, Lyman-$α$ transits do not primarily probe the mass-loss rates. Instead, they inform us about the velocity at which the escape mechanism is ejecting material from the planet, providing a clean test of predictions from atmospheric escape models. △ Less

Submitted 18 November, 2022; v1 submitted 11 November, 2021; originally announced November 2021.

Comments: Accepted for publication in MNRAS

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

Authors: Hilke E. Schlichting, Edward D. Young

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

Submitted 8 March, 2022; v1 submitted 21 July, 2021; originally announced July 2021.

Comments: Submitted to PSJ, 25 pages, 18 figures

Photoevaporation vs. core-powered mass-loss: model comparison with the 3D radius gap

Authors: James G. Rogers, Akash Gupta, James E. Owen, Hilke E. Schlichting

Abstract: The EUV/X-ray photoevaporation and core-powered mass-loss models are both capable of reproducing the bimodality in the sizes of small, close-in exoplanets observed by the \textit{Kepler} space mission, often referred to as the "radius gap". However, it is unclear which of these two mechanisms dominates the atmospheric mass-loss which is likely sculpting the radius gap. In this work, we propose a n… ▽ More The EUV/X-ray photoevaporation and core-powered mass-loss models are both capable of reproducing the bimodality in the sizes of small, close-in exoplanets observed by the \textit{Kepler} space mission, often referred to as the "radius gap". However, it is unclear which of these two mechanisms dominates the atmospheric mass-loss which is likely sculpting the radius gap. In this work, we propose a new method of differentiating between the two models, which relies on analysing the radius gap in 3D parameter space. Using models for both mechanisms, and by performing synthetic transit surveys we predict the size and characteristics of a survey capable of discriminating between the two models. We find that a survey of $\gtrsim 5000$ planets, with a wide range in stellar mass and measurement uncertainties at a $\lesssim 5\%$ level is sufficient. Our methodology is robust against moderate false positive contamination of $\lesssim 10\%$. We perform our analysis on two surveys (which do not satisfy our requirements): the \textit{California Kepler Survey} and the \textit{Gaia-Kepler Survey} and find, unsurprisingly, that both data-sets are consistent with either model. We propose a hypothesis test to be performed on future surveys which can robustly ascertain which of the two mechanisms formed the radius gap, provided one dominates over the other. △ Less

Submitted 7 October, 2021; v1 submitted 7 May, 2021; originally announced May 2021.

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

To cool is to keep: Residual H/He atmospheres of super-Earths and sub-Neptunes

Authors: William Misener, Hilke E. Schlichting

Abstract: Super-Earths and sub-Neptunes are commonly thought to have accreted hydrogen/helium envelopes, consisting of a few to ten percent of their total mass, from the primordial gas disk. Subsequently, hydrodynamic escape driven by core-powered mass-loss and/or photo-evaporation likely stripped much of these primordial envelopes from the lower-mass and closer-in planets to form the super-Earth population… ▽ More Super-Earths and sub-Neptunes are commonly thought to have accreted hydrogen/helium envelopes, consisting of a few to ten percent of their total mass, from the primordial gas disk. Subsequently, hydrodynamic escape driven by core-powered mass-loss and/or photo-evaporation likely stripped much of these primordial envelopes from the lower-mass and closer-in planets to form the super-Earth population. In this work we show that after undergoing core-powered mass-loss, some super-Earths can retain small residual H/He envelopes. This retention is possible because, for significantly depleted atmospheres, the density at the radiative-convective boundary drops sufficiently such that thhe cooling time-scale becomes less than the mass-loss time-scale. The residual envelope is therefore able to contract, terminating further mass loss. Using analytic calculations and numerical simulations, we show that the mass of primordial H/He envelope retained as a fraction of the planet's total mass, $f_{ret}$, increases with increasing planet mass, $M_{c}$, and decreases with increasing equilibrium temperature, $T_{eq}$, scaling as $f_{ret} \propto M_{c}^{3/2} T_{eq}^{-1/2} \exp{[M_{c}^{3/4} T_{eq}^{-1}]}$. $f_{ret}$ varies from $< 10^{-8}$ to about $10^{-3}$ for typical super-Earth parameters. To first order, the exact amount of left-over H/He depends on the initial envelope mass, the planet mass, its equilibrium temperature, and the envelope's opacity. These residual hydrogen envelopes reduce the atmosphere's mean molecular weight compared to a purely secondary atmosphere, a signature observable by current and future facilities. These remnant atmospheres may, however, in many cases be vulnerable to long-term erosion by photo-evaporation. Any residual hydrogen envelope likely plays an important role in the long-term physical evolution of super-Earths, including their geology and geochemistry. △ Less

Submitted 16 March, 2021; originally announced March 2021.

Comments: 17 pages, 11 figures. Resubmitted to MNRAS after revision

Caught in the Act: Core-powered Mass-loss Predictions for Observing Atmospheric Escape

Authors: Akash Gupta, Hilke E. Schlichting

Abstract: Past studies have demonstrated that atmospheric escape by the core-powered mass-loss mechanism can explain a multitude of observations associated with the radius valley that separates the super-Earth and sub-Neptune planet populations. Complementing such studies, in this work, we present a shortlist of planets that could be losing their atmospheres today if their evolution is indeed primarily dict… ▽ More Past studies have demonstrated that atmospheric escape by the core-powered mass-loss mechanism can explain a multitude of observations associated with the radius valley that separates the super-Earth and sub-Neptune planet populations. Complementing such studies, in this work, we present a shortlist of planets that could be losing their atmospheres today if their evolution is indeed primarily dictated by core-powered mass-loss. We use Bayesian inference analysis on our planet evolution and mass-loss model to estimate the posteriors of the parameters that encapsulate the current state of a given planet, given their published masses, radii and host star properties. Our models predict that the following planets could be losing their atmospheres today at a rate $\gtrsim 10^7$ g/s at 50\% confidence level: pi Men c, Kepler-60 d, Kepler-60 b, HD 86226 c, EPIC 249893012 b, Kepler-107 c, HD 219134 b, Kepler-80 e, Kepler-138 d and GJ 9827 d. As a by-product of our Bayesian inference analysis, we were also able to identify planets that most-likely harbor either secondary atmospheres abundant with high mean-molecular weight species, low-density interiors abundant with ices, or both. The planets belonging to this second category are WASP-47 e, Kepler-78 b, Kepler-10 b, CoRoT-7 b, HD 80653 b, 55 Cnc e and Kepler-36 b. While the aforementioned lists are by no means exhaustive, we believe that candidates presented here can serve as useful input for target selection for future surveys and for testing the importance of core-powered mass-loss in individual planetary systems. △ Less

Submitted 16 April, 2021; v1 submitted 15 March, 2021; originally announced March 2021.

Comments: 16 pages, 7 figures, and 3 tables. Accepted to MNRAS

Where are the Extrasolar Mercuries?

Authors: Alexandra E. Doyle, Beth Klein, Hilke E. Schlichting, Edward D. Young

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

Submitted 13 August, 2020; originally announced August 2020.

Comments: 22 pages, 13 figures; accepted for publication in ApJ

Journal ref: ApJ 901 10 (2020)

Losing Oceans: The Effects of Composition on the Thermal Component of Impact-driven Atmospheric Loss

Authors: John B. Biersteker, Hilke E. Schlichting

Abstract: The formation of the solar system's terrestrial planets concluded with a period of giant impacts. Previous works examining the volatile loss caused by the impact shock in the moon-forming impact find atmospheric losses of at most 20-30 per cent and essentially no loss of oceans. However, giant impacts also result in thermal heating, which can lead to significant atmospheric escape via a Parker-typ… ▽ More The formation of the solar system's terrestrial planets concluded with a period of giant impacts. Previous works examining the volatile loss caused by the impact shock in the moon-forming impact find atmospheric losses of at most 20-30 per cent and essentially no loss of oceans. However, giant impacts also result in thermal heating, which can lead to significant atmospheric escape via a Parker-type wind. Here we show that H2O and other high-mean molecular weight outgassed species can be efficiently lost through this thermal wind if present in a hydrogen-dominated atmosphere, substantially altering the final volatile inventory of terrestrial planets. Specifically, we demonstrate that a giant impact of a Mars-sized embryo with a proto-Earth can remove several Earth oceans' worth of H2O, and other heavier volatile species, together with a primordial hydrogen-dominated atmosphere. These results may offer an explanation for the observed depletion in Earth's light noble gas budget and for its depleted xenon inventory, which suggest that Earth underwent significant atmospheric loss by the end of its accretion. Because planetary embryos are massive enough to accrete primordial hydrogen envelopes and because giant impacts are stochastic and occur concurrently with other early atmospheric evolutionary processes, our results suggest a wide diversity in terrestrial planet volatile budgets. △ Less

Submitted 2 July, 2020; originally announced July 2020.

Comments: submitted to MNRAS, 9 pages, 4 figures

Near-equilibrium isotope fractionation during planetesimal evaporation

Authors: E. D. Young, A. Shahar, F. Nimmo, H. E. Schlichting, E. A. Schauble, H. Tang, J. Labidi

Abstract: Silicon and Mg in differentiated rocky bodies exhibit heavy isotope enrichments that have been attributed to evaporation of partially or entirely molten planetesimals. We evaluate the mechanisms of planetesimal evaporation in the early solar system and the conditions that controlled attendant isotope fractionations. Energy balance at the surface of a body accreted within ~1 Myr of CAI formation an… ▽ More Silicon and Mg in differentiated rocky bodies exhibit heavy isotope enrichments that have been attributed to evaporation of partially or entirely molten planetesimals. We evaluate the mechanisms of planetesimal evaporation in the early solar system and the conditions that controlled attendant isotope fractionations. Energy balance at the surface of a body accreted within ~1 Myr of CAI formation and heated from within by 26Al decay results in internal temperatures exceeding the silicate solidus, producing a transient magma ocean with a thin surface boundary layer of order < 1 meter that would be subject to foundering. Bodies that are massive enough to form magma oceans by radioisotope decay (ge 0.1%) can retain hot rock vapor even in the absence of ambient nebular gas. We find that a steady-state rock vapor forms within minutes to hours and results from a balance between rates of magma evaporation and atmospheric escape. Vapor pressure buildup adjacent to the surfaces of the evaporating magmas would have inevitably led to an approach to equilibrium isotope partitioning between the vapor phase and the silicate melt. Numerical simulations of this near-equilibrium evaporation process for a body with a radius of ~ 700 km yield a steady-state far-field vapor pressure corresponding to 95% saturation. Approaches to equilibrium isotope fractionation between vapor and melt should have been the norm during planet formation due to the formation of steady-state rock vapor atmospheres and/or the presence of protostellar gas. We model the Si and Mg isotopic composition of bulk Earth and show that the best fit is for a carbonaceous chondrite-like source material with about 12% loss of Mg and 15% loss of Si resulting from near-equilibrium evaporation into the solar protostellar disk of hydrogen gas on timescales of 10,000 to 100,000 years. △ Less

Submitted 12 September, 2019; originally announced September 2019.

Comments: 35 pages, 15 figures

Journal ref: Icarus 323, 1-15 (2019)

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

Authors: Akash Gupta, Hilke E. Schlichting

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

Submitted 9 January, 2020; v1 submitted 8 July, 2019; originally announced July 2019.

Comments: 16 pages and 11 figures. Submitted to MNRAS. Revised version

New Constraints From Dust Lines On The Surface Densities Of Protoplanetary Disks

Authors: Diana Powell, Ruth Murray-Clay, Laura M. Pérez, Hilke E. Schlichting, Mickey Rosenthal

Abstract: We present new determinations of disk surface density, independent of an assumed dust opacity, for a sample of 7 bright, diverse protoplanetary disks using measurements of disk dust lines. We develop a robust method for determining the location of dust lines by modeling disk interferometric visibilities at multiple wavelengths. The disks in our sample have newly derived masses that are 9-27% of th… ▽ More We present new determinations of disk surface density, independent of an assumed dust opacity, for a sample of 7 bright, diverse protoplanetary disks using measurements of disk dust lines. We develop a robust method for determining the location of dust lines by modeling disk interferometric visibilities at multiple wavelengths. The disks in our sample have newly derived masses that are 9-27% of their host stellar mass, substantially larger than the minimum mass solar nebula. All are stable to gravitational collapse except for one which approaches the limit of Toomre-Q stability. Our mass estimates are 2-15 times larger than estimates from integrated optically thin dust emission. We derive depleted dust-to-gas ratios with typical values of ~$10^{-3}$ in the outer disk. Using coagulation models we derive dust surface density profiles that are consistent with millimeter dust observations. In these models, the disks formed with an initial dust mass that is a factor of ~10 greater than is presently observed. Of the three disks in our sample with resolved CO line emission, the masses of HD 163296, AS 209, and TW Hya are roughly 3, 115, and 40 times more massive than estimates from CO respectively. This range indicates that CO depletion is not uniform across different disks and that dust is a more robust tracer of total disk mass. Our method of determining surface density using dust lines is robust even if particles form as aggregates and is useful even in the presence of dust substructure caused by pressure traps. The low Toomre-Q values observed in this sample indicate that at least some disks do not accrete efficiently. △ Less

Submitted 8 May, 2019; originally announced May 2019.

Comments: 29 pages, 12 figures, submitted to ApJ, revised following comments from referee

Atmospheric Mass Loss from High Velocity Giant Impacts

Authors: Almog Yalinewich, Hilke E. Schlichting

Abstract: Using moving mesh hydrodynamic simulations, we determine the shock propagation and resulting ground velocities for a planet hit by a high velocity impactor. We use our results to determine the atmospheric mass loss caused by the resulting ground motion due to the impact shock wave. We find that there are two distinct shock propagation regimes: In the limit in which the impactor is significantly sm… ▽ More Using moving mesh hydrodynamic simulations, we determine the shock propagation and resulting ground velocities for a planet hit by a high velocity impactor. We use our results to determine the atmospheric mass loss caused by the resulting ground motion due to the impact shock wave. We find that there are two distinct shock propagation regimes: In the limit in which the impactor is significantly smaller than the target ($R_i<< R_t$), the solutions are self-similar and the shock velocity at a fixed point on the target scale as $m_i^{2/3}$, where $m_i$ is the mass of the impactor. In addition, the ground velocities follow a universal profile given by $v_g/v_i=(14.2x^2-25.3x+11.3)/(x^2-2.5x+1.9) +2\ln{R_i/R_t}$, where $x=\sin\left(θ/2\right)$, $θ$ is the latitude on the target measured from the impact site, and $v_g$ and $v_i$ are the ground velocity and impact velocity, respectively. In contrast, in the limit in which the impactor is comparable to the size of the target ($R_i \sim R_t$), we find that shock velocities decline with the mass of the impactor significantly more weakly than $m_i^{2/3}$. We use the resulting surface velocity profiles to calculate the atmospheric mass loss for a large range of impactor masses and impact velocities and apply them to the Kepler-36 system and the Moon forming impact. Finally, we present and generalise our results in terms of the $v_g/v_i$ and the impactor to target size ratio ($R_i/R_t$) such that they can easily be applied to other collision scenarios. △ Less

Submitted 16 April, 2019; v1 submitted 28 November, 2018; originally announced November 2018.

Sculpting the Valley in the Radius Distribution of Small Exoplanets as a by-product of Planet Formation: The Core-Powered Mass-Loss Mechanism

Authors: Akash Gupta, Hilke E. Schlichting

Abstract: Recent observations revealed a bimodal radius distribution of small, short-period exoplanets with a paucity in their occurrence, a radius `valley', around $1.5-2.0$ $R_\oplus$. In this work, we investigate the effect of a planet's own cooling luminosity on its thermal evolution and atmospheric mass-loss (core-powered mass-loss) and determine its observational consequences for the radius distributi… ▽ More Recent observations revealed a bimodal radius distribution of small, short-period exoplanets with a paucity in their occurrence, a radius `valley', around $1.5-2.0$ $R_\oplus$. In this work, we investigate the effect of a planet's own cooling luminosity on its thermal evolution and atmospheric mass-loss (core-powered mass-loss) and determine its observational consequences for the radius distribution of small, close-in exoplanets. Using simple analytical descriptions and numerical simulations, we demonstrate that planetary evolution based on the core-powered mass-loss mechanism alone (i.e., without any photoevaporation) can produce the observed valley in the radius distribution. Our results match the valley's location, shape and slope in planet radius-orbital period parameter space, and the relative magnitudes of the planet occurrence rate above and below the valley. We find that the slope of the valley is, to first order, dictated by the atmospheric mass-loss timescale at the Bondi radius and given by $\text{d log} R_p/ \text{d log} P \simeq 1/(3(1-β))$ which evaluates to $ -0.11$ for $ β\simeq 4$, where $M_c/M_\oplus = (R_c/R_{\oplus})^β (ρ_{c*}/ρ_{\oplus})^{β/3}$ is the mass-radius relation of the core. This choice for $β$ yields good agreement with observations and attests to the significance of internal compression for massive planetary cores. We further find that the location of the valley scales as $ρ_{c*}^{-4/9}$ and that the observed planet population must have predominantly rocky cores with typical water-ice fractions of less than $\sim 20\%$. Furthermore, we show that the relative magnitude of the planet occurrence rate above and below the valley is sensitive to the details of the planet-mass distribution but that the location of the valley is not. △ Less

Submitted 18 April, 2019; v1 submitted 7 November, 2018; originally announced November 2018.

Comments: 11 pages and 6 figures. Revised version. Submitted to MNRAS

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 487, Issue 1, July 2019, Pages 24-33

Atmospheric mass loss due to giant impacts: the importance of the thermal component for hydrogen-helium envelopes

Authors: John B. Biersteker, Hilke E. Schlichting

Abstract: Systems of close-in super-Earths display striking diversity in planetary bulk density and composition. Giant impacts are expected to play a role in the formation of many of these worlds. Previous works, focused on the mechanical shock caused by a giant impact, have shown that these impacts can eject large fractions of the planetary envelope, offering a partial explanation for the observed spread i… ▽ More Systems of close-in super-Earths display striking diversity in planetary bulk density and composition. Giant impacts are expected to play a role in the formation of many of these worlds. Previous works, focused on the mechanical shock caused by a giant impact, have shown that these impacts can eject large fractions of the planetary envelope, offering a partial explanation for the observed spread in exoplanet compositions. Here, we examine the thermal consequences of giant impacts, and show that the atmospheric loss caused by these effects can significantly exceed that caused by mechanical shocks for hydrogen-helium (H/He) envelopes. When a giant impact occurs, part of the impact energy is converted into thermal energy, heating the rocky core and the envelope. We find that the ensuing thermal expansion of the envelope can lead to a period of sustained, rapid mass loss through a Parker wind, resulting in the partial or complete erosion of the H/He envelope. The fraction of the envelope lost depends on the planet's orbital distance from its host star and its initial thermal state, and hence age. Planets closer to their host stars are more susceptible to thermal atmospheric loss triggered by impacts than ones on wider orbits. Similarly, younger planets, with rocky cores which are still hot and molten from formation, suffer greater atmospheric loss. This is especially interesting because giant impacts are expected to occur $10{-}100~\mathrm{Myr}$ after formation. For planets where the thermal energy of the core is much greater than the envelope energy, the impactor mass required for significant atmospheric removal is $M_\mathrm{imp} / M_p \sim μ/ μ_c \sim 0.1$, approximately the ratio of the heat capacities of the envelope and core. When the envelope energy dominates the total energy budget, complete loss can occur when the impactor mass is comparable to the envelope mass. △ Less

Submitted 18 September, 2018; originally announced September 2018.

Comments: 10 pages, 9 figures

Formation of Super-Earths

Authors: Hilke E Schlichting

Abstract: Super-Earths are the most abundant planets known to date and are characterized by having sizes between that of Earth and Neptune, typical orbital periods of less than 100 days and gaseous envelopes that are often massive enough to significantly contribute to the planet's overall radius. Furthermore, super-Earths regularly appear in tightly-packed multiple-planet systems, but resonant configuration… ▽ More Super-Earths are the most abundant planets known to date and are characterized by having sizes between that of Earth and Neptune, typical orbital periods of less than 100 days and gaseous envelopes that are often massive enough to significantly contribute to the planet's overall radius. Furthermore, super-Earths regularly appear in tightly-packed multiple-planet systems, but resonant configurations in such systems are rare. This chapters summarizes current super-Earth formation theories. It starts from the formation of rocky cores and subsequent accretion of gaseous envelopes. We follow the thermal evolution of newly formed super-Earths and discuss their atmospheric mass loss due to disk dispersal, photoevaporation, core-cooling and collisions. We conclude with a comparison of observations and theoretical predictions, highlighting that even super-Earths that appear as barren rocky cores today likely formed with primordial hydrogen and helium envelopes and discuss some paths forward for the future. △ Less

Submitted 8 February, 2018; originally announced February 2018.

Comments: Invited review accepted for publication in the 'Handbook of Exoplanets,' Planet Formation section, Springer Reference Works, Juan Antonio Belmonte and Hans Deeg, Eds

Core-powered mass loss and the radius distribution of small exoplanets

Authors: Sivan Ginzburg, Hilke E. Schlichting, Re'em Sari

Abstract: Recent observations identify a valley in the radius distribution of small exoplanets, with planets in the range $1.5-2.0\,{\rm R}_{\oplus}$ significantly less common than somewhat smaller or larger planets. This valley may suggest a bimodal population of rocky planets that are either engulfed by massive gas envelopes that significantly enlarge their radius, or do not have detectable atmospheres at… ▽ More Recent observations identify a valley in the radius distribution of small exoplanets, with planets in the range $1.5-2.0\,{\rm R}_{\oplus}$ significantly less common than somewhat smaller or larger planets. This valley may suggest a bimodal population of rocky planets that are either engulfed by massive gas envelopes that significantly enlarge their radius, or do not have detectable atmospheres at all. One explanation of such a bimodal distribution is atmospheric erosion by high-energy stellar photons. We investigate an alternative mechanism: the luminosity of the cooling rocky core, which can completely erode light envelopes while preserving heavy ones, produces a deficit of intermediate sized planets. We evolve planetary populations that are derived from observations using a simple analytical prescription, accounting self-consistently for envelope accretion, cooling and mass loss, and demonstrate that core-powered mass loss naturally reproduces the observed radius distribution, regardless of the high-energy incident flux. Observations of planets around different stellar types may distinguish between photoevaporation, which is powered by the high-energy tail of the stellar radiation, and core-powered mass loss, which depends on the bolometric flux through the planet's equilibrium temperature that sets both its cooling and mass-loss rates. △ Less

Submitted 6 March, 2018; v1 submitted 4 August, 2017; originally announced August 2017.

Comments: Matched to published version

Journal ref: MNRAS 476, 759-765 (2018)

Constraints on the Growth and Spin of the Supermassive Black Hole in M32 From High Cadence Visible Light Observations

Authors: R. Chary, G. Hallinan, L. K. Harding, N. S. Saini, H. E. Schlichting

Abstract: We present 1-second cadence observations of M32 (NGC221) with the CHIMERA instrument at the Hale 200-inch telescope of the Palomar Observatory. Using field stars as a baseline for relative photometry, we are able to construct a light curve of the nucleus in the g-prime and r-prime band with 1sigma=36 milli-mag photometric stability. We derive a temporal power spectrum for the nucleus and find no e… ▽ More We present 1-second cadence observations of M32 (NGC221) with the CHIMERA instrument at the Hale 200-inch telescope of the Palomar Observatory. Using field stars as a baseline for relative photometry, we are able to construct a light curve of the nucleus in the g-prime and r-prime band with 1sigma=36 milli-mag photometric stability. We derive a temporal power spectrum for the nucleus and find no evidence for a time-variable signal above the noise as would be expected if the nuclear black hole were accreting gas. Thus, we are unable to constrain the spin of the black hole although future work will use this powerful instrument to target more actively accreting black holes. Given the black hole mass of (2.5+/-0.5)*10^6 Msun inferred from stellar kinematics, the absence of a contribution from a nuclear time-variable signal places an upper limit on the accretion rate which is 4.6*10^{-8} of the Eddington rate, a factor of two more stringent than past upper limits from HST. The low mass of the black hole despite the high stellar density suggests that the gas liberated by stellar interactions was primarily at early cosmic times when the low-mass black hole had a small Eddington luminosity. This is at least partly driven by a top-heavy stellar initial mass function at early cosmic times which is an efficient producer of stellar mass black holes. The implication is that supermassive black holes likely arise from seeds formed through the coalescence of 3-100 Msun mass black holes that then accrete gas produced through stellar interaction processes. △ Less

Submitted 1 May, 2017; originally announced May 2017.

Comments: 8 pages, 3 figures, submitted to the Astrophysical Journal, comments welcome

Avoiding resonance capture in multi-planet extrasolar systems

Authors: Margaret Pan, Hilke E. Schlichting

Abstract: A commonly noted feature of the population of multi-planet extrasolar systems is the rarity of planet pairs in low-order mean-motion resonances. We revisit the physics of resonance capture via convergent disk-driven migration. We point out that for planet spacings typical of stable configurations for Kepler systems, the planets can routinely maintain a small but nonzero eccentricity due to gravita… ▽ More A commonly noted feature of the population of multi-planet extrasolar systems is the rarity of planet pairs in low-order mean-motion resonances. We revisit the physics of resonance capture via convergent disk-driven migration. We point out that for planet spacings typical of stable configurations for Kepler systems, the planets can routinely maintain a small but nonzero eccentricity due to gravitational perturbations from their neighbors. Together with the upper limit on the migration rate needed for capture, the finite eccentricity can make resonance capture difficult or impossible in Sun-like systems for planets smaller than ~Neptune-sized. This mass limit on efficient capture is broadly consistent with observed exoplanet pairs that have mass determinations: of pairs with the heavier planet exterior to the lighter planet -- which would have been undergoing convergent migration in their disks -- those in or nearly in resonance are much more likely to have total mass greater than two Neptune masses than to have smaller masses. The agreement suggests that the observed paucity of resonant pairs around sun-like stars may simply arise from a small resonance capture probability for lower-mass planets. Planet pairs that thereby avoid resonance capture are much less likely to collide in an eventual close approach than to simply migrate past one another to become a divergently migrating pair with the lighter planet exterior. For systems around M stars we expect resonant pairs to be much more common, since there the minimum mass threshhold for efficient capture is about an Earth mass. △ Less

Submitted 28 April, 2017; v1 submitted 25 April, 2017; originally announced April 2017.

Comments: 14 pages, 4 figures; submitted to ApJ; corrected minor typos

Using Ice and Dust Lines to Constrain the Surface Densities of Protoplanetary Disks

Authors: Diana Powell, Ruth Murray-Clay, Hilke E. Schlichting

Abstract: We present a novel method for determining the surface density of protoplanetary disks through consideration of disk 'dust lines' which indicate the observed disk radial scale at different observational wavelengths. This method relies on the assumption that the processes of particle growth and drift control the radial scale of the disk at late stages of disk evolution such that the lifetime of the… ▽ More We present a novel method for determining the surface density of protoplanetary disks through consideration of disk 'dust lines' which indicate the observed disk radial scale at different observational wavelengths. This method relies on the assumption that the processes of particle growth and drift control the radial scale of the disk at late stages of disk evolution such that the lifetime of the disk is equal to both the drift timescale and growth timescale of the maximum particle size at a given dust line. We provide an initial proof of concept of our model through an application to the disk TW Hya and are able to estimate the disk dust-to-gas ratio, CO abundance, and accretion rate in addition to the total disk surface density. We find that our derived surface density profile and dust-to-gas ratio are consistent with the lower limits found through measurements of HD gas. The CO ice line also depends on surface density through grain adsorption rates and drift and we find that our theoretical CO ice line estimates have clear observational analogues. We further apply our model to a large parameter space of theoretical disks and find three observational diagnostics that may be used to test its validity. First we predict that the dust lines of disks other than TW Hya will be consistent with the normalized CO surface density profile shape for those disks. Second, surface density profiles that we derive from disk ice lines should match those derived from disk dust lines. Finally, we predict that disk dust and ice lines will scale oppositely, as a function of surface density, across a large sample of disks. △ Less

Submitted 16 August, 2019; v1 submitted 15 April, 2017; originally announced April 2017.

Comments: Published in ApJ, May 2017

CHIMERA: a wide-field, multi-color, high-speed photometer at the prime focus of the Hale telescope

Authors: Leon K. Harding, Gregg Hallinan, Jennifer Milburn, Paul Gardner, Nick Konidaris, Navtej Singh, Michael Shao, Jagmit Sandhu, Gillian Kyne, Hilke E. Schlichting

Abstract: The Caltech HIgh-speed Multi-color camERA (CHIMERA) is a new instrument that has been developed for use at the prime focus of the Hale 200-inch telescope. Simultaneous optical imaging in two bands is enabled by a dichroic beam splitter centered at 567 nm, with Sloan u' and g' bands available on the blue arm and Sloan r', i' and z_s' bands available on the red arm. Additional narrow-band filters wi… ▽ More The Caltech HIgh-speed Multi-color camERA (CHIMERA) is a new instrument that has been developed for use at the prime focus of the Hale 200-inch telescope. Simultaneous optical imaging in two bands is enabled by a dichroic beam splitter centered at 567 nm, with Sloan u' and g' bands available on the blue arm and Sloan r', i' and z_s' bands available on the red arm. Additional narrow-band filters will also become available as required. An Electron Multiplying CCD (EMCCD) detector is employed for both optical channels, each capable of simultaneously delivering sub-electron effective read noise under multiplication gain and frame rates of up to 26 fps full frame (several 1000 fps windowed), over a fully corrected 5 x 5 arcmin field of view. CHIMERA was primarily developed to enable the characterization of the size distribution of sub-km Kuiper Belt Objects via stellar occultation, a science case that motivates the frame-rate, the simultaneous multi-color imaging and the wide field of view of the instrument. In addition, it also has unique capability in the detection of faint near-Earth asteroids and will be used for the monitoring of short duration transient and periodic sources, particularly those discovered by the intermediate Palomar Transient Factory (iPTF), and the upcoming Zwicky Transient Facility (ZTF). △ Less

Submitted 12 January, 2016; originally announced January 2016.

Comments: Accepted to Mon. Not. R. Astron. Soc. (MNRAS), 15 pages, 9 figures, 3 tables

Journal ref: 2016, MNRAS, 457 (3), p. 3036-3049

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

Authors: Sivan Ginzburg, Hilke E. Schlichting, Re'em Sari

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

Submitted 16 April, 2016; v1 submitted 24 December, 2015; originally announced December 2015.

Comments: Updated (refereed) version

A Search for Ringed Exoplanets using Kepler Photometry

Authors: Matthew Z. Heising, Geoffrey W. Marcy, Hilke E. Schlichting

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

Submitted 16 November, 2015; v1 submitted 3 November, 2015; originally announced November 2015.

Comments: 41 pages, 11 figures, Accepted for publication in ApJ; corrected typos, added discussion of photometric signal of exomoons in section 2.3

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

Authors: Niraj K. Inamdar, Hilke E. Schlichting

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

Submitted 10 December, 2015; v1 submitted 7 October, 2015; originally announced October 2015.

Comments: 14 pages, 4 figures. Updated figures and with referee's comments addressed and incorporated

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

Authors: Niraj K. Inamdar, Hilke E. Schlichting

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

Submitted 13 November, 2016; v1 submitted 14 December, 2014; originally announced December 2014.

Comments: 10 pages, 10 figures. Abstract abridged for submission. Fixed typographical error in Eq. 6 exponent

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

Authors: Hilke E. Schlichting

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

Submitted 4 October, 2014; originally announced October 2014.

Comments: Accepted for publication in ApJL

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

Authors: Peter Goldreich, Hilke E. Schlichting

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

Submitted 21 November, 2013; v1 submitted 21 August, 2013; originally announced August 2013.

Comments: Accepted for publication in AJ. Corrected typos, updated references, added Appendix

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

Authors: Sean N. Raymond, Hilke E. Schlichting, Franck Hersant, Franck Selsis

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

Submitted 18 June, 2013; originally announced June 2013.

Comments: Icarus, in press. 33 pages, 5 figures

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

Authors: Hilke E. Schlichting, Cesar I. Fuentes, David E. Trilling

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

Submitted 1 April, 2013; v1 submitted 30 January, 2013; originally announced January 2013.

Comments: Updated Figures to include latest TAOS results and updated references

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

Authors: Hilke E. Schlichting, Eran O. Ofek, Re'em Sari, Edmund P. Nelan, Avishay Gal-Yam, Michael Wenz, Philip Muirhead, Nikta Javanfar, Mario Livio

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

Submitted 30 October, 2012; originally announced October 2012.

Comments: Accepted for publication in the Astrophysical Journal, 32 pages, 15 figures. arXiv admin note: text overlap with arXiv:0912.2996

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

Authors: Hilke E. Schlichting, Paul H. Warren, Qing-Zhu Yin

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

Submitted 22 April, 2012; v1 submitted 28 February, 2012; originally announced February 2012.

Comments: accepted for publication in ApJ, 9 pages, 4 figures; minor corrections, additional references added

Self-consistent size and velocity distributions of collisional cascades

Authors: Margaret Pan, Hilke E. Schlichting

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

Submitted 2 November, 2011; originally announced November 2011.

Comments: 23 pages, 3 figures, 1 table; submitted to ApJ