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

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

Overstability of the 2:1 mean motion resonance: Exploring disc parameters with hydrodynamic simulations

Authors: Zahra Afkanpour, Sareh Ataiee, Alexandros Ziampras, Anna B. T. Penzlin, Rafael Sfair, Christoph Schäfer, Wilhelm Kley, Hilke Schlichting

Abstract: Resonant planetary migration in protoplanetary discs can lead to an interplay between the resonant interaction of planets and their disc torques called overstability. While theoretical predictions and N-body simulations hinted at its existence, there was no conclusive evidence until hydrodynamical simulations were performed. Our primary purpose is to find a hydrodynamic setup that induces overst… ▽ More Resonant planetary migration in protoplanetary discs can lead to an interplay between the resonant interaction of planets and their disc torques called overstability. While theoretical predictions and N-body simulations hinted at its existence, there was no conclusive evidence until hydrodynamical simulations were performed. Our primary purpose is to find a hydrodynamic setup that induces overstability in a planetary system with two moderate-mass planets in a first-order 2:1 mean motion resonance. We also aim to analyse the impact of key disc parameters, namely the viscosity, surface density, and aspect ratio, on the occurrence of overstability in this planetary system when the masses of the planets are kept constant. We performed 2D locally isothermal hydrodynamical simulations of two planets, with masses of 5 and 10 $M_{\oplus}$, in a 2:1 resonance. Upon identifying the fiducial model in which the system exhibits overstability, we performed simulations with different disc parameters to explore the effects of the disc on the overstability of the system. We observe an overstable planetary system in our hydrodynamic simulations. In the parameter study, we note that overstability occurs in discs characterised by low surface density and low viscosity. Increasing the surface density reduces the probability of overstability within the system. A limit cycle was observed in a specific viscous model with $αあるふぁ_νにゅー = 10^{-3}$. In almost all our models, planets create partial gaps in the disc, which affects both the migration timescale and structure of the planetary system. We demonstrate the existence of overstability using hydrodynamic simulations but find deviations from the analytic approximation and show that the main contribution to this deviation can be attributed to dynamic gap opening. △ Less

Submitted 8 April, 2024; originally announced April 2024.

Comments: 16 pages (including appendices), 27 figures, 3 tables, accepted for publication by A&A

Journal ref: A&A 686, A277 (2024)

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

The DECam Ecliptic Exploration Project (DEEP) II. Observational Strategy and Design

Authors: Chadwick A. Trujillo, Cesar Fuentes, David W. Gerdes, Larissa Markwardt, Scott S. Sheppard, Ryder Strauss, Colin Orion Chandler, William J. Oldroyd, David E. Trilling, Hsing Wen Lin, Fred C. Adams, Pedro H. Bernardinelli, Matthew J. Holman, Mario Juric, Andrew McNeill, Michael Mommert, Kevin J. Napier, Matthew J. Payne, Darin Ragozzine, Andrew S. Rivkin, Hilke Schlichting, Hayden Smotherman

Abstract: We present the DECam Ecliptic Exploration Project (DEEP) survey strategy including observing cadence for orbit determination, exposure times, field pointings and filter choices. The overall goal of the survey is to discover and characterize the orbits of a few thousand Trans-Neptunian Objects (TNOs) using the Dark Energy Camera (DECam) on the Cerro Tololo Inter-American Observatory (CTIO) Blanco 4… ▽ More We present the DECam Ecliptic Exploration Project (DEEP) survey strategy including observing cadence for orbit determination, exposure times, field pointings and filter choices. The overall goal of the survey is to discover and characterize the orbits of a few thousand Trans-Neptunian Objects (TNOs) using the Dark Energy Camera (DECam) on the Cerro Tololo Inter-American Observatory (CTIO) Blanco 4 meter telescope. The experiment is designed to collect a very deep series of exposures totaling a few hours on sky for each of several 2.7 square degree DECam fields-of-view to achieve a magnitude of about 26.2 using a wide VR filter which encompasses both the V and R bandpasses. In the first year, several nights were combined to achieve a sky area of about 34 square degrees. In subsequent years, the fields have been re-visited to allow TNOs to be tracked for orbit determination. When complete, DEEP will be the largest survey of the outer solar system ever undertaken in terms of newly discovered object numbers, and the most prolific at producing multi-year orbital information for the population of minor planets beyond Neptune at 30 au. △ Less

Submitted 30 October, 2023; originally announced October 2023.

Comments: 29 pages, 4 figures and 4 tables

The DECam Ecliptic Exploration Project (DEEP) VI: first multi-year observations of trans-Neptunian objects

Authors: Hayden Smotherman, Pedro H. Bernardinelli, Stephen K. N. Portillo, Andrew J. Connolly, J. Bryce Kalmbach, Steven Stetzler, Mario Juric, Dino Bektesvic, Zachary Langford, Fred C. Adams, William J. Oldroyd, Matthew J. Holman, Colin Orion Chandler, Cesar Fuentes, David W. Gerdes, Hsing Wen Lin, Larissa Markwardt, Andrew McNeill, Michael Mommert, Kevin J. Napier, Matthew J. Payne, Darin Ragozzine, Andrew S. Rivkin, Hilke Schlichting, Scott S. Sheppard , et al. (3 additional authors not shown)

Abstract: We present the first set of trans-Neptunian objects (TNOs) observed on multiple nights in data taken from the DECam Ecliptic Exploration Project (DEEP). Of these 110 TNOs, 105 do not coincide with previously known TNOs and appear to be new discoveries. Each individual detection for our objects resulted from a digital tracking search at TNO rates of motion, using two to four hour exposure sets, and… ▽ More We present the first set of trans-Neptunian objects (TNOs) observed on multiple nights in data taken from the DECam Ecliptic Exploration Project (DEEP). Of these 110 TNOs, 105 do not coincide with previously known TNOs and appear to be new discoveries. Each individual detection for our objects resulted from a digital tracking search at TNO rates of motion, using two to four hour exposure sets, and the detections were subsequently linked across multiple observing seasons. This procedure allows us to find objects with magnitudes $m_{VR} \approx 26$. The object discovery processing also included a comprehensive population of objects injected into the images, with a recovery and linking rate of at least $94\%$. The final orbits were obtained using a specialized orbit fitting procedure that accounts for the positional errors derived from the digital tracking procedure. Our results include robust orbits and magnitudes for classical TNOs with absolute magnitudes $H \sim 10$, as well as a dynamically detached object found at 76 au (semi-major axis $a\approx 77 \, \mathrm{au}$). We find a disagreement between our population of classical TNOs and the CFEPS-L7 three component model for the Kuiper belt. △ Less

Submitted 5 October, 2023; originally announced October 2023.

Comments: Accepted to AJ, companion paper do DEEP III. Objects will be released in the journal version (or contacting the authors)

The DECam Ecliptic Exploration Project (DEEP) III: Survey characterization and simulation methods

Authors: Pedro H. Bernardinelli, Hayden Smotherman, Zachary Langford, Stephen K. N. Portillo, Andrew J. Connolly, J. Bryce Kalmbach, Steven Stetzler, Mario Juric, William J. Oldroyd, Hsing Wen Lin, Fred C. Adams, Colin Orion Chandler, Cesar Fuentes, David W. Gerdes, Matthew J. Holman, Larissa Markwardt, Andrew McNeill, Michael Mommert, Kevin J. Napier, Matthew J. Payne, Darin Ragozzine, Andrew S. Rivkin, Hilke Schlichting, Scott S. Sheppard, Ryder Strauss , et al. (2 additional authors not shown)

Abstract: We present a detailed study of the observational biases of the DECam Ecliptic Exploration Project's (DEEP) B1 data release and survey simulation software that enables direct statistical comparisons between models and our data. We inject a synthetic population of objects into the images, and then subsequently recover them in the same processing as our real detections. This enables us to characteriz… ▽ More We present a detailed study of the observational biases of the DECam Ecliptic Exploration Project's (DEEP) B1 data release and survey simulation software that enables direct statistical comparisons between models and our data. We inject a synthetic population of objects into the images, and then subsequently recover them in the same processing as our real detections. This enables us to characterize the survey's completeness as a function of apparent magnitudes and on-sky rates of motion. We study the statistically optimal functional form for the magnitude, and develop a methodology that can estimate the magnitude and rate efficiencies for all survey's pointing groups simultaneously. We have determined that our peak completeness is on average 80\% in each pointing group, and our magnitude drops to $25\%$ of this value at $m_{25} = 26.22$. We describe the freely available survey simulation software and its methodology. We conclude by using it to infer that our effective search area for objects at 40 au is $14.8°^2$, and that our lack of dynamically cold distant objects means that there at most $8\times 10^3$ objects with $60 < a < 80$ au and absolute magnitudes $H \leq 8$. △ Less

Submitted 5 October, 2023; originally announced October 2023.

Comments: Accepted to AJ, companion paper to DEEP VI

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.

The DECam Ecliptic Exploration Project (DEEP): I. Survey description, science questions, and technical demonstration

Authors: David E. Trilling, David W. Gerdes, Mario Juric, Chadwick A. Trujillo, Pedro H. Bernardinelli, Kevin J. Napier, Hayden Smotherman, Ryder Strauss, Cesar Fuentes, Matthew J. Holman, Hsing Wen Lin, Larissa Markwardt, Andrew McNeill, Michael Mommert, William J. Oldroyd, Matthew J. Payne, Darin Ragozzine, Andrew S. Rivkin, Hilke Schlichting, Scott S. Sheppard, Fred C. Adams, Colin Orion Chandler

Abstract: We present here the DECam Ecliptic Exploration Project (DEEP), a three year NOAO/NOIRLab Survey that was allocated 46.5 nights to discover and measure the properties of thousands of trans-Neptunian objects (TNOs) to magnitudes as faint as VR~27, corresponding to sizes as small as 20 km diameter. In this paper we present the science goals of this project, the experimental design of our survey, and… ▽ More We present here the DECam Ecliptic Exploration Project (DEEP), a three year NOAO/NOIRLab Survey that was allocated 46.5 nights to discover and measure the properties of thousands of trans-Neptunian objects (TNOs) to magnitudes as faint as VR~27, corresponding to sizes as small as 20 km diameter. In this paper we present the science goals of this project, the experimental design of our survey, and a technical demonstration of our approach. The core of our project is "digital tracking," in which all collected images are combined at a range of motion vectors to detect unknown TNOs that are fainter than the single exposure depth of VR~23 mag. Through this approach we reach a depth that is approximately 2.5 magnitudes fainter than the standard LSST "wide fast deep" nominal survey depth of 24.5 mag. DEEP will more than double the number of known TNOs with observational arcs of 24 hours or more, and increase by a factor of 10 or more the number of known small (<50 km) TNOs. We also describe our ancillary science goals, including measuring the mean shape distribution of very small main belt asteroids, and briefly outline a set of forthcoming papers that present further aspects of and preliminary results from the DEEP program. △ Less

Submitted 6 September, 2023; originally announced September 2023.

Comments: AJ, in press. First in a series of papers

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

The UV-SCOPE Mission: Ultraviolet Spectroscopic Characterization Of Planets and their Environments

Authors: David R. Ardila, Evgenya Shkolnik, John Ziemer, Mark Swain, James E. Owen, Michael Line, R. O. Parke Loyd, R. Glenn Sellar, Travis Barman, Courtney Dressing, William Frazier, April D. Jewell, Robert J. Kinsey, Carl C. Liebe, Joshua D. Lothringer, Luz Maria Martinez-Sierra, James McGuire, Victoria Meadows, Ruth Murray-Clay, Shouleh Nikzad, Sarah Peacock, Hilke Schlichting, David Sing, Kevin Stevenson, Yen-Hung Wu

Abstract: UV-SCOPE is a mission concept to determine the causes of atmospheric mass loss in exoplanets, investigate the mechanisms driving aerosol formation in hot Jupiters, and study the influence of the stellar environment on atmospheric evolution and habitability. As part of these investigations, the mission will generate a broad-purpose legacy database of time-domain ultraviolet (UV) spectra for nearly… ▽ More UV-SCOPE is a mission concept to determine the causes of atmospheric mass loss in exoplanets, investigate the mechanisms driving aerosol formation in hot Jupiters, and study the influence of the stellar environment on atmospheric evolution and habitability. As part of these investigations, the mission will generate a broad-purpose legacy database of time-domain ultraviolet (UV) spectra for nearly 200 stars and planets. The observatory consists of a 60 cm, f/10 telescope paired to a long-slit spectrograph, yielding simultaneous, almost continuous coverage between 1203 Å and 4000 Å, with resolutions ranging from 6000 to 240. The efficient instrument provides throughputs > 4% (far-UV; FUV) and > 15% (near-UV; NUV), comparable to HST/COS and much better than HST/STIS, over the same spectral range. A key design feature is the LiF prism, which serves as a dispersive element and provides high throughput even after accounting for radiation degradation. The use of two delta-doped Electron-Multiplying CCD detectors with UV-optimized, single-layer anti-reflection coatings provides high quantum efficiency and low detector noise. From the Earth-Sun second Lagrangian point, UV-SCOPE will continuously observe planetary transits and stellar variability in the full FUV-to-NUV range, with negligible astrophysical background. All these features make UV-SCOPE the ideal instrument to study exoplanetary atmospheres and the impact of host stars on their planets. UV-SCOPE was proposed to NASA as a Medium Explorer (MidEx) mission for the 2021 Announcement of Opportunity. If approved, the observatory will be developed over a 5-year period. Its primary science mission takes 34 months to complete. The spacecraft carries enough fuel for 6 years of operations. △ Less

Submitted 19 August, 2022; originally announced August 2022.

Comments: 12 pages, 7 figures, 1 table. Conference presentation, 17 July 2022, SPIE Astronomical Telescopes and Instrumentation, Montreal, Canada

Multiwavelength Vertical Structure in the AUえーゆー Mic Debris Disk: Characterizing the Collisional Cascade

Authors: David Vizgan, A. Meredith Hughes, Evan S. Carter, Kevin M. Flaherty, Margaret Pan, Eugene Chiang, Hilke Schlichting, David J. Wilner, Sean M. Andrews, John M. Carpenter, Atilla Moór, Meredith A. MacGregor

Abstract: Debris disks are scaled-up analogs of the Kuiper Belt in which dust is generated by collisions between planetesimals. In the "collisional cascade" model of debris disks, dust lost to radiation pressure and winds is constantly replenished by grinding collisions between planetesimals. The model assumes that collisions are destructive and involve large velocities; this assumption has not been tested… ▽ More Debris disks are scaled-up analogs of the Kuiper Belt in which dust is generated by collisions between planetesimals. In the "collisional cascade" model of debris disks, dust lost to radiation pressure and winds is constantly replenished by grinding collisions between planetesimals. The model assumes that collisions are destructive and involve large velocities; this assumption has not been tested beyond our Solar System. We present 0"25 ($\approx$2.4 au) resolution observations of the $λらむだ$ = 450 $μみゅー$m dust continuum emission from the debris disk around the nearby M dwarf AUえーゆー Microscopii with the Atacama Large Millimeter/submillimeter Array. We use parametric models to describe the disk structure, and an MCMC algorithm to explore the posterior distributions of the model parameters; we fit the structure of the disk to both our data and archival $λらむだ= 1.3$ mm data (Daley et al. 2019), from which we obtain two aspect ratio measurements at 1.3 mm ($h_{1300}$ = 0.025$^{+0.008}_{-0.002}$) and at 450 $μみゅー$m ($h_{450}$ = 0.019$^{+0.006}_{-0.001}$), as well as the grain size distribution index $q =$ 3.03 $\pm$ 0.02. Contextualizing our aspect ratio measurements within the modeling framework laid out in Pan & Schlichting (2012), we derive a power law index of velocity dispersion as a function of grain size $p = 0.28 \pm 0.06$ for the AUえーゆー Mic debris disk. This result implies that smaller bodies are more easily disrupted than larger bodies by collisions, which is inconsistent with the strength regime usually assumed for such small bodies. Possible explanations for this discrepancy are discussed. △ Less

Submitted 11 July, 2022; originally announced July 2022.

Comments: 17 pages, 6 figures, 3 tables. Accepted for publication in ApJ

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

AUえーゆー Microscopii in the FUV: Observations in Quiescence, During Flares, and Implications for AUえーゆー Mic b and c

Authors: Adina D. Feinstein, Kevin France, Allison Youngblood, Girish M. Duvvuri, DJ Teal, P. Wilson Cauley, Darryl Z. Seligman, Eric Gaidos, Eliza M. R. Kempton, Jacob L. Bean, Hannah Diamond-Lowe, Elisabeth Newton, Sivan Ginzburg, Peter Plavchan, Peter Gao, Hilke Schlichting

Abstract: High energy X-ray and ultraviolet (UV) radiation from young stars impacts planetary atmospheric chemistry and mass loss. The active $\sim 22$ Myr M dwarf AUえーゆー Mic hosts two exoplanets orbiting interior to its debris disk. Therefore, this system provides a unique opportunity to quantify the effects of stellar XUV irradiation on planetary atmospheres as a function of both age and orbital separation. I… ▽ More High energy X-ray and ultraviolet (UV) radiation from young stars impacts planetary atmospheric chemistry and mass loss. The active $\sim 22$ Myr M dwarf AUえーゆー Mic hosts two exoplanets orbiting interior to its debris disk. Therefore, this system provides a unique opportunity to quantify the effects of stellar XUV irradiation on planetary atmospheres as a function of both age and orbital separation. In this paper we present over 5 hours of Far-UV (FUV) observations of AUえーゆー Mic taken with the Cosmic Origins Spectrograph (COS; 1070-1360 Angstrom) on the Hubble Space Telescope (HST). We provide an itemization of $120$ emission features in the HST/COS FUV spectrum and quantify the flux contributions from formation temperatures ranging from $10^4-10^7$ K. We detect 13flares in the FUV white-light curve with energies ranging from $10^{29} - 10^{31}$ ergs. The majority of the energy in each of these flares is released from the transition region between the chromosphere and the corona. There is a 100$\times$ increase in flux at continuum wavelengths $λらむだ< 1100$ Angstrom in each flare which may be caused by thermal Bremsstrahlung emission. We calculate that the baseline atmospheric mass-loss rate for AUえーゆー Mic b is $\sim 10^8$ g s$^{-1}$, although this rate can be as high as $\sim 10^{14}$ g s$^{-1}$ during flares with $L_\textrm{flare} \simeq 10^{33}$ erg s$^{-1}$. Finally, we model the transmission spectra for AUえーゆー Mic b and c with a new panchromatic spectrum of AUえーゆー Mic c and motivate future JWST observations of these planets. △ Less

Submitted 12 August, 2022; v1 submitted 19 May, 2022; originally announced May 2022.

Comments: 29 pages, 12 figures, 3 tables, Accepted to AJ

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

Oxygen fugacities of extrasolar rocks: Evidence for an Earth-like geochemistry of exoplanets

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

Abstract: Oxygen fugacity is a measure of rock oxidation that influences planetary structure and evolution. Most rocky bodies in the Solar System formed at oxygen fugacities approximately five orders of magnitude higher than a hydrogen-rich gas of solar composition. It is unclear whether this oxidation of rocks in the Solar System is typical among other planetary systems. We exploit the elemental abundances… ▽ More Oxygen fugacity is a measure of rock oxidation that influences planetary structure and evolution. Most rocky bodies in the Solar System formed at oxygen fugacities approximately five orders of magnitude higher than a hydrogen-rich gas of solar composition. It is unclear whether this oxidation of rocks in the Solar System is typical among other planetary systems. We exploit the elemental abundances observed in six white dwarfs polluted by the accretion of rocky bodies to determine the fraction of oxidized iron in those extrasolar rocky bodies and therefore their oxygen fugacities. The results are consistent with the oxygen fugacities of Earth, Mars, and typical asteroids in the Solar System, suggesting that at least some rocky exoplanets are geophysically and geochemically similar to Earth. △ Less

Submitted 28 October, 2019; originally announced October 2019.

Comments: 12 pages of main text, 3 figures in the main text

Journal ref: Doyle et al. (2019) Science 366 (6463), 356-359

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

The Mass of Stirring Bodies in the AUえーゆー Mic Debris Disk Inferred from Resolved Vertical Structure

Authors: Cail Daley, A. Meredith Hughes, Evan S. Carter, Kevin Flaherty, Zachary Lambros, Margaret Pan, Hilke Schlichting, Eugene Chiang, Mark Wyatt, David Wilner, Sean Andrews, John Carpenter

Abstract: The vertical distribution of dust in debris disks is sensitive to the number and size of large planetesimals dynamically stirring the disk, and is therefore well-suited for constraining the prevalence of otherwise unobservable Uranus and Neptune analogs. Information regarding stirring bodies has previously been inferred from infrared and optical observations of debris disk vertical structure, but… ▽ More The vertical distribution of dust in debris disks is sensitive to the number and size of large planetesimals dynamically stirring the disk, and is therefore well-suited for constraining the prevalence of otherwise unobservable Uranus and Neptune analogs. Information regarding stirring bodies has previously been inferred from infrared and optical observations of debris disk vertical structure, but theoretical works predict that the small particles traced by short-wavelength observations will be `puffed up' by radiation pressure, yielding only upper limits. The large grains that dominate the disk emission at millimeter wavelengths are much less sensitive to the effects of stellar radiation or stellar winds, and therefore trace the underlying mass distribution more directly. Here we present ALMA 1.3 mm dust continuum observations of the debris disk around the nearby M star AUえーゆー Mic. The 3 au spatial resolution of the observations, combined with the favorable edge-on geometry of the system, allows us to measure the vertical thickness of the disk. We report a scale height-to-radius aspect ratio of $h = 0.031_{-0.004}^{+0.005}$ between radii of $\sim 23$ au and $\sim 41$ au. Comparing this aspect ratio to a theoretical model of size-dependent velocity distributions in the collisional cascade, we find that the perturbing bodies embedded in the local disk must be larger than about 400 km, and the largest perturbing body must be smaller than roughly 1.8 M$_\odot$. These measurements rule out the presence of a gas giant or Neptune analog near the $\sim 40$ au outer edge of the debris ring, but are suggestive of large planetesimals or an Earth-sized planet stirring the dust distribution. △ Less

Submitted 8 April, 2019; v1 submitted 29 March, 2019; originally announced April 2019.

Comments: 21 pages, 8 figures, 5 tables. Accepted for publication in ApJ

The Super-Earth Opportunity - Search for Habitable Exoplanets in the 2020s

Authors: Renyu Hu, Charles A. Beichman, David Brain, Pin Chen, Mario Damiano, Rebekah Dawson, A. James Friedson, Yasuhiro Hasagawa, Andrew Howard, Robert Johnson, Tiffany Kataria, Richard Kidd, Edwin Kite, Heather Knutson, Wladimir Lyra, Michael Mischna, Noah Planavsky, Chris Reinhard, Hilke Schlichting, Sara Seager, Christophe Sotin, Mark Swain, Neal Turner, Robert West, Yuk Yung , et al. (1 additional authors not shown)

Abstract: The recent discovery of a staggering diversity of planets beyond the Solar System has brought with it a greatly expanded search space for habitable worlds. The Kepler exoplanet survey has revealed that most planets in our interstellar neighborhood are larger than Earth and smaller than Neptune. Collectively termed super-Earths and mini-Neptunes, some of these planets may have the conditions to sup… ▽ More The recent discovery of a staggering diversity of planets beyond the Solar System has brought with it a greatly expanded search space for habitable worlds. The Kepler exoplanet survey has revealed that most planets in our interstellar neighborhood are larger than Earth and smaller than Neptune. Collectively termed super-Earths and mini-Neptunes, some of these planets may have the conditions to support liquid water oceans, and thus Earth-like biology, despite differing in many ways from our own planet. In addition to their quantitative abundance, super-Earths are relatively large and are thus more easily detected than true Earth twins. As a result, super-Earths represent a uniquely powerful opportunity to discover and explore a panoply of fascinating and potentially habitable planets in 2020 - 2030 and beyond. △ Less

Submitted 12 March, 2019; originally announced March 2019.

Comments: Science white paper submitted to the 2020 Astronomy and Astrophysics Decadal Survey

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

Determining Exoplanetary Oblateness Using Transit Depth Variations

Authors: John B. Biersteker, Hilke Schlichting

Abstract: The measurement of an exoplanet's oblateness and obliquity provides insights into the planet's internal structure and formation history. Previous work using small differences in the shape of the transit light curve has been moderately successful, but was hampered by the small signal and extreme photometric precision required. The measurement of changes in transit depth, caused by the spin precessi… ▽ More The measurement of an exoplanet's oblateness and obliquity provides insights into the planet's internal structure and formation history. Previous work using small differences in the shape of the transit light curve has been moderately successful, but was hampered by the small signal and extreme photometric precision required. The measurement of changes in transit depth, caused by the spin precession of an oblate planet, was proposed as an alternative method. Here, we present the first attempt to measure these changes. Using Kepler photometry, we examined the brown dwarf Kepler-39b and the warm Saturn Kepler-427b. We could not reliably constrain the oblateness of Kepler-39b. We find transit depth variations for Kepler-427b at $90.1\%$ significance ($1.65σしぐま$) consistent with a precession period of $P_\mathrm{prec} = 5.45^{+0.46}_{-0.37}~\mathrm{years}$ and an oblateness, $f~=~0.19^{+0.32}_{-0.16}$. This oblateness is comparable to Solar System gas giants, and would raise questions about the dynamics and tidal synchronization of Kepler-427b. △ Less

Submitted 29 August, 2017; originally announced August 2017.

Comments: 20 pages, 14 figures, accepted for publication in AJ

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