Source author record

Eve J. Lee

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

19works
4topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

19 published item(s)

preprint2022arXiv

Creating the Radius Gap without Mass Loss

The observed exoplanet population features a gap in the radius distribution that separates the smaller super-Earths ($\lesssim$1.7 Earth radii) from the larger sub-Neptunes ($\sim$1.7--4 Earth radii). While mass loss theories can explain many of the observed features of this radius valley, it is difficult to reconcile them with potentially rising population of terrestrials beyond orbital periods of $\sim$30 days. We investigate the ability of gas accretion during the gas-poor phase of disk evolution to reproduce both the location of the observed radius gap and the existence of long-period terrestrial planets. Updating the analytic scaling relations of gas accretion rate accounting for the shrinking of the bound radius by hydrodynamic effects and deriving a more realistic disk temperature profile, we find that the late-stage gas accretion alone is able to carve out the observed radius gap, with slopes $R_{\rm gap} \propto P^{-0.096}$ and $R_{\rm gap} \propto M_\star^{0.15}$ for top-heavy; and $R_{\rm gap} \propto P^{-0.089}$ and $R_{\rm gap} \propto M_\star^{0.22}$ for bottom-heavy core mass distributions, in good agreement with observations. The general morphology of the primordial radius gap is stable against a range of disk gas density and disk accretion rate with the latter affecting mostly the population of large planets ($\gtrsim$3--4$R_\oplus$). The peaks and valleys in the radius distribution were likely set in place primordially while post-formation mass loss further tunes the exoplanetary population. We provide potential observational tests that may be possible with TESS, PLATO and Roman Space Telescope.

preprint2021arXiv

Sculpting the sub-Saturn Occurrence Rate via Atmospheric Mass Loss

The sub-Saturn ($\sim$4--8$R_{\oplus}$) occurrence rate rises with orbital period out to at least $\sim$300 days. In this work we adopt and test the hypothesis that the decrease in their occurrence towards the star is a result of atmospheric mass loss, which can transform sub-Saturns into sub-Neptunes ($\lesssim$4$R_{\oplus}$) more efficiently at shorter periods. We show that under the mass loss hypothesis, the sub-Saturn occurrence rate can be leveraged to infer their underlying core mass function, and by extension that of gas giants. We determine that lognormal core mass functions peaked near $\sim$10--20$M_{\oplus}$ are compatible with the sub-Saturn period distribution, the distribution of observationally-inferred sub-Saturn cores, and gas accretion theories. Our theory predicts that close-in sub-Saturns should be $\sim$50\% less common and $\sim$30\% more massive around rapidly rotating stars; this should be directly testable for stars younger than $\lesssim$500 Myr. We also predict that the sub-Jovian desert becomes less pronounced and opens up at smaller orbital periods around M stars compared to solar-type stars ($\sim$0.7 days vs.~$\sim$3 days). We demonstrate that exceptionally low-density sub-Saturns, "Super-Puffs", can survive intense hydrodynamic escape to the present day if they are born with even larger atmospheres than they currently harbor; in this picture, Kepler 223 d began with an envelope $\sim$1.5$\times$ the mass of its core and is currently losing its envelope at a rate $\sim$2$\times 10^{-3}M_{\oplus}~\mathrm{Myr}^{-1}$. If the predictions from our theory are confirmed by observations, the core mass function we predict can also serve to constrain core formation theories of gas-rich planets.

preprint2020arXiv

Diffuser-Assisted Infrared Transit Photometry for Four Dynamically Interacting \textit{Kepler} Systems

We present ground-based infrared transit observations for four dynamically interacting \textit{Kepler} planets, including Kepler-29b, Kepler-36c, KOI-1783.01, and Kepler-177c, obtained using the Wide-field Infrared Camera on the Hale 200" telescope at Palomar Observatory. By utilizing an engineered diffuser and custom guiding software, we mitigate time-correlated telluric and instrumental noise sources in these observations. We achieve an infrared photometric precision comparable to or better than that of space-based observatories such as the \textit{Spitzer Space Telescope}, and detect transits with greater than 3$σ$ significance for all planets. For Kepler-177c ($J=13.9$) our measurement uncertainties are only $1.2\times$ the photon noise limit and 1.9 times better than the predicted photometric precision for \textit{Spitzer} IRAC photometry of this same target. We find that a single transit observation obtained $4-5$ years after the end of the original \textit{Kepler} mission can reduce dynamical mass uncertainties by as much as a factor of three for these systems. Additionally, we combine our new observations of KOI-1783.01 with information from the literature to confirm the planetary nature of this system. We discuss the implications of our new mass and radius constraints in the context of known exoplanets with low incident fluxes, and we note that Kepler-177c may be a more massive analog to the currently known super-puffs given its core mass (3.8$\pm0.9M_\Earth$) and large gas-to-core ratio (2.8$\pm0.7$). Our demonstrated infrared photometric performance opens up new avenues for ground-based observations of transiting exoplanets previously thought to be restricted to space-based investigation.

preprint2020arXiv

Forming Diverse Super-Earth Systems in Situ

Super-Earths and mini-Neptunes exhibit great diversity in their compositional and orbital properties. Their bulk densities span a large range, from those dense enough to be purely rocky to those needing a substantial contribution from volatiles to their volumes. Their orbital configurations range from compact, circular multi-transiting systems like Kepler-11 to systems like our Solar System's terrestrial planets with wider spacings and modest but significant eccentricities and mutual inclinations. Here we investigate whether a continuum of formation conditions resulting from variation in the amount of solids available in the inner disk can account for the diversity of orbital and compositional properties observed for super Earths, including the apparent dichotomy between single transiting and multiple transiting system. We simulate in situ formation of super-Earths via giant impacts and compare to the observed Kepler sample. We find that intrinsic variations among disks in the amount of solids available for in situ formation can account for the orbital and compositional diversity observed among Kepler's transiting planets. Our simulations can account for the planets' distributions of orbital period ratios, transit duration ratios, and transit multiplicity; higher eccentricities for single than multi transiting planets; smaller eccentricities for larger planets; scatter in the mass-radius relation, including lower densities for planets with masses measured with TTVs than RVs; and similarity in planets' sizes and spacings within each system. Our findings support the theory that variation among super-Earth and mini-Neptune properties is primarily locked in by different in situ formation conditions, rather than arising stochastically through subsequent evolution.

preprint2020arXiv

Most Stars (and Planets?) Are Born in Intense Radiation Fields

Protostars and young stars are strongly spatially "clustered" or "correlated" within their natal giant molecular clouds (GMCs). We demonstrate that such clustering leads to the conclusion that the incident bolometric radiative flux upon a random young star/disc is enhanced (relative to volume-averaged fluxes) by a factor which increases with the total stellar mass of the complex. Because the Galactic cloud mass function is top-heavy, the typical star in our Galaxy experienced a much stronger radiative environment than those forming in well-observed nearby (but relatively small) clouds, exceeding fluxes in the Orion Nebular Cluster by factors of $\gtrsim$30. Heating of the circumstellar disc around a median young star is dominated by this external radiation beyond $\sim 50\,$AU. And if discs are not well-shielded by ambient dust, external UV irradiation can dominate over the host star down to sub-AU scales. Another consequence of stellar clustering is an extremely broad Galaxy-wide distribution of incident flux (spanning $>10$ decades), with half the Galactic star formation in a substantial "tail" towards even more intense background radiation. We also show that the strength of external irradiation is amplified super-linearly in high-density environments such as the Galactic centre, starbursts, or high-redshift galaxies.

preprint2020arXiv

Primordial Radius Gap and Potentially Broad Core Mass Distributions of Super-Earths and Sub-Neptunes

The observed radii distribution of {\it Kepler} exoplanets reveals two distinct populations: those that are more likely to be terrestrials ($\lesssim1.7R_\oplus$) and those that are more likely to be gas-enveloped ($\gtrsim2R_\oplus$). There exists a clear gap in the distribution of radii that separates these two kinds of planets. Mass loss processes like photoevaporation by high energy photons from the host star have been proposed as natural mechanisms to carve out this radius valley. These models favor underlying core mass function of sub-Neptunes that is sharply peaked at $\sim$4--8$M_\oplus$ but the radial-velocity follow-up of these small planets hint at a more bottom-heavy mass function. By taking into account the initial gas accretion in gas-poor (but not gas-empty) nebula, we demonstrate that 1) the observed radius valley is a robust feature that is initially carved out at formation during late-time gas accretion; and 2) that it can be reconciled with core mass functions that are broad extending well into sub-Earth regime. The maximally cooled isothermal limit prohibits cores lighter than $\sim$1--2$M_\oplus$ from accreting enough mass to appear gas-enveloped. The rocky-to-enveloped transition established at formation produces a gap in the radius distribution that shifts to smaller radii farther from the star, similar to that observed. For the best agreement with the data, our late-time gas accretion model favors dust-free accretion in hotter disks with cores slightly less dense than the Earth ($\sim$0.8$ρ_\oplus$) drawn from a mass function that is as broad as $dN/dM_{\rm core} \propto M_{\rm core}^{-0.7}$.

preprint2020arXiv

The First Habitable Zone Earth-Sized Planet From TESS II: $Spitzer$ Confirms TOI-700 d

We present $Spitzer$ 4.5$μ$m observations of the transit of TOI-700 d, a habitable zone Earth-sized planet in a multiplanet system transiting a nearby M-dwarf star (TIC 150428135, 2MASS J06282325-6534456). TOI-700 d has a radius of $1.144^{+0.062}_{-0.061}R_\oplus$ and orbits within its host star's conservative habitable zone with a period of 37.42 days ($T_\mathrm{eq} \sim 269$K). TOI-700 also hosts two small inner planets (R$_b$=$1.037^{+0.065}_{-0.064}R_\oplus$ & R$_c$=$2.65^{+0.16}_{-0.15}R_\oplus$) with periods of 9.98 and 16.05 days, respectively. Our $Spitzer$ observations confirm the TESS detection of TOI-700 d and remove any remaining doubt that it is a genuine planet. We analyze the $Spitzer$ light curve combined with the 11 sectors of TESS observations and a transit of TOI-700 c from the LCOGT network to determine the full system parameters. Although studying the atmosphere of TOI-700 d is not likely feasible with upcoming facilities, it may be possible to measure the mass of TOI-700 d using state-of-the-art radial velocity instruments (expected RV semi-amplitude of $\sim$70 cm/s).

preprint2020arXiv

The First Habitable Zone Earth-sized Planet from TESS. I: Validation of the TOI-700 System

We present the discovery and validation of a three-planet system orbiting the nearby (31.1 pc) M2 dwarf star TOI-700 (TIC 150428135). TOI-700 lies in the TESS continuous viewing zone in the Southern Ecliptic Hemisphere; observations spanning 11 sectors reveal three planets with radii ranging from 1 R$_\oplus$ to 2.6 R$_\oplus$ and orbital periods ranging from 9.98 to 37.43 days. Ground-based follow-up combined with diagnostic vetting and validation tests enable us to rule out common astrophysical false-positive scenarios and validate the system of planets. The outermost planet, TOI-700 d, has a radius of $1.19\pm0.11$ R$_\oplus$ and resides in the conservative habitable zone of its host star, where it receives a flux from its star that is approximately 86% of the Earth's insolation. In contrast to some other low-mass stars that host Earth-sized planets in their habitable zones, TOI-700 exhibits low levels of stellar activity, presenting a valuable opportunity to study potentially-rocky planets over a wide range of conditions affecting atmospheric escape. While atmospheric characterization of TOI-700 d with the James Webb Space Telescope (JWST) will be challenging, the larger sub-Neptune, TOI-700 c (R = 2.63 R$_\oplus$), will be an excellent target for JWST and beyond. TESS is scheduled to return to the Southern Hemisphere and observe TOI-700 for an additional 11 sectors in its extended mission, which should provide further constraints on the known planet parameters and searches for additional planets and transit timing variations in the system.

preprint2020arXiv

WASP-107b's density is even lower: a case study for the physics of planetary gas envelope accretion and orbital migration

With a mass in the Neptune regime and a radius of Jupiter, WASP-107b presents a challenge to planet formation theories. Meanwhile, the planet's low surface gravity and the star's brightness also make it one of the most favorable targets for atmospheric characterization. Here, we present the results of an extensive 4-year Keck/HIRES radial-velocity (RV) follow-up program of the WASP-107 system and provide a detailed study of the physics governing the accretion of its gas envelope. We reveal that WASP-107b's mass is only 1.8 Neptune masses ($M_b = 30.5 \pm 1.7$ $M_\oplus$). The resulting extraordinarily low density suggests that WASP-107b has a H/He envelope mass fraction of $> 85$% unless it is substantially inflated. The corresponding core mass of $<4.6$ $M_\oplus$ at 3$σ$ is significantly lower than what is traditionally assumed to be necessary to trigger massive gas envelope accretion. We demonstrate that this large gas-to-core mass ratio most plausibly results from the onset of accretion at $\gtrsim 1$ AU onto a low-opacity, dust-free atmosphere and subsequent migration to the present-day $a_b = 0.0566 \pm 0.0017$ AU. Beyond WASP-107b, we also detect a second more massive planet ($M_c \sin i = 0.36 \pm 0.04$ $M_{J}$) on a wide eccentric orbit ($e_c = 0.28 \pm 0.07$) which may have influenced the orbital migration and spin-orbit misalignment of WASP-107b. Overall, our new RV observations and envelope accretion modeling provide crucial insights into the intriguing nature of WASP-107b and the system's formation history. Looking ahead, WASP-107b will be a keystone planet to understand the physics of gas envelope accretion.

preprint2016arXiv

A Primer on Unifying Debris Disk Morphologies

A "minimum model" for debris disks consists of a narrow ring of parent bodies, secularly forced by a single planet on a possibly eccentric orbit, colliding to produce dust grains that are perturbed by stellar radiation pressure. We demonstrate how this minimum model can reproduce a wide variety of disk morphologies imaged in scattered starlight. Five broad categories of disk shape can be captured: "rings," "needles," "ships-and-wakes," "bars," and "moths (a.k.a. fans)," depending on the viewing geometry. Moths can also sport "double wings." We explain the origin of morphological features from first principles, exploring the dependence on planet eccentricity, disk inclination dispersion, and the parent body orbital phases at which dust grains are born. A key determinant in disk appearance is the degree to which dust grain orbits are apsidally aligned. Our study of a simple steady-state (secularly relaxed) disk should serve as a reference for more detailed models tailored to individual systems. We use the intuition gained from our guidebook of disk morphologies to interpret, informally, the images of a number of real-world debris disks. These interpretations suggest that the farthest reaches of planetary systems are perturbed by eccentric planets, possibly just a few Earth masses each.

preprint2016arXiv

Bringing "The Moth" to Light: A Planet-Sculpting Scenario for the HD 61005 Debris Disk

The HD 61005 debris disk ("The Moth") stands out from the growing collection of spatially resolved circumstellar disks by virtue of its unusual swept-back morphology, brightness asymmetries, and dust ring offset. Despite several suggestions for the physical mechanisms creating these features, no definitive answer has been found. In this work, we demonstrate the plausibility of a scenario in which the disk material is shaped dynamically by an eccentric, inclined planet. We present new Keck NIRC2 scattered-light angular differential imaging of the disk at 1.2-2.3 microns that further constrains its outer morphology (projected separations of 27-135 AU). We also present complementary Gemini Planet Imager 1.6 micron total intensity and polarized light detections that probe down to projected separations less than 10 AU. To test our planet-sculpting hypothesis, we employed secular perturbation theory to construct parent body and dust distributions that informed scattered-light models. We found that this method produced models with morphological and photometric features similar to those seen in the data, supporting the premise of a planet-perturbed disk. Briefly, our results indicate a disk parent body population with a semimajor axis of 40-52 AU and an interior planet with an eccentricity of at least 0.2. Many permutations of planet mass and semimajor axis are allowed, ranging from an Earth mass at 35 AU to a Jupiter mass at 5 AU.

preprint2016arXiv

Correlations between compositions and orbits established by the giant impact era of planet formation

The giant impact phase of terrestrial planet formation establishes connections between super-Earths' orbital properties (semimajor axis spacings, eccentricities, mutual inclinations) and interior compositions (the presence or absence of gaseous envelopes). Using N-body simulations and analytic arguments, we show that spacings derive not only from eccentricities, but also from inclinations. Flatter systems attain tighter spacings, a consequence of an eccentricity equilibrium between gravitational scatterings, which increase eccentricities, and mergers, which damp them. Dynamical friction by residual disk gas plays a critical role in regulating mergers and in damping inclinations and eccentricities. Systems with moderate gas damping and high solid surface density spawn gas-enveloped super-Earths with tight spacings, small eccentricities, and small inclinations. Systems in which super-Earths coagulate without as much ambient gas, in disks with low solid surface density, produce rocky planets with wider spacings, larger eccentricities, and larger mutual inclinations. A combination of both populations can reproduce the observed distributions of spacings, period ratios, transiting planet multiplicities, and transit duration ratios exhibited by Kepler super-Earths. The two populations, both formed in situ, also help to explain observed trends of eccentricity vs. planet size, and bulk density vs. method of mass measurement (radial velocities vs. transit timing variations). Simplifications made in this study --- including the limited timespan of the simulations, and the approximate treatments of gas dynamical friction and gas depletion history --- should be improved upon in future work to enable a detailed quantitative comparison to the observations.

preprint2016arXiv

Observational Evidence of Dynamic Star Formation Rate in Milky Way Giant Molecular Clouds

Star formation on galactic scales is known to be a slow process, but whether it is slow on smaller scales is uncertain. We cross-correlate 5469 giant molecular clouds (GMCs) from a new all-sky catalog with 256 star forming complexes (SFCs) to build a sample of 191 SFC-GMC complexes---collections of multiple clouds each matched to 191 SFCs. The total mass in stars harbored by these clouds is inferred from WMAP free-free fluxes. We measure the GMC mass, the virial parameter, the star formation efficiency $ε$ and the star formation rate per free-fall time $ε_{\rm ff}$. Both $ε$ and $ε_{\rm ff}$ range over 3--4 orders of magnitude. We find that 68.3% of the clouds fall within $σ_{\logε}=0.79\pm0.22\,{\rm dex}$ and $σ_{\logε_{\rm ff}}=0.91\pm0.22\,{\rm dex}$ about the median. Compared to these observed scatters, a simple model with a time independent $ε_{\rm ff}$ that depends on the host GMC properties predicts $σ_{\logε_{\rm ff}}=$0.12-0.24. Allowing for a time-variable $ε_{\rm ff}$, we can recover the large dispersion in the rate of star formation. This strongly suggests that star formation in the Milky Way is a dynamic process on GMC scales. We also show that the surface star formation rate profile of the Milky Way correlates well with the molecular gas surface density profile.

preprint2016arXiv

Two Transiting Low Density Sub-Saturns from K2

We report the discovery and confirmation of two sub-Saturn planets orbiting a bright (V = 11.3), metal-rich ([Fe/H] = 0.42 $\pm$ 0.04 dex) G3 dwarf in the K2 Campaign 2 field. The planets are 5.68 $\pm$ 0.56 Earth-radii and 7.82 $\pm$ 0.72 Earth-radii and have orbital periods of 20.8851 $\pm$ 0.0003 d and 42.3633$\pm$0.0006 d, near to the 2:1 mean-motion resonance. We obtained 32 radial velocities (RVs) with Keck/HIRES and detected the reflex motion due to EPIC-203771098b and c. These planets have masses of 21.0 $\pm$ 5.4 Earth-masses and 27.0 $\pm$ 6.9 Earth-masses, respectively. With low densities of 0.63 $\pm$ 0.25 g/cc and 0.31 $\pm$ 0.12 g/cc, respectively, the planets require thick envelopes of H/He to explain their large sizes and low masses. Interior structure models predict that the planets have fairly massive cores of 17.6 $\pm$ 4.3 Earth-masses and 16.1 $\pm$ 4.2 Earth-masses, respectively. They may have formed exterior to their present locations, accreted their H/He envelopes at large orbital distances, and migrated in as a resonant pair. The proximity to resonance, large transit depths, and host star brightness offer rich opportunities for TTV follow-up. Finally, the low surface gravities of the EPIC-203771098 planets make them favorable targets for transmission spectroscopy by HST, Spitzer, and JWST.

preprint2015arXiv

A Metallicity Recipe for Rocky Planets

Planets with sizes between those of Earth and Neptune divide into two populations: purely rocky bodies whose atmospheres contribute negligibly to their sizes, and larger gas-enveloped planets possessing voluminous and optically thick atmospheres. We show that whether a planet forms rocky or gas-enveloped depends on the solid surface density of its parent disk. Assembly times for rocky cores are sensitive to disk solid surface density. Lower surface densities spawn smaller planetary embryos; to assemble a core of given mass, smaller embryos require more mergers between bodies farther apart and therefore exponentially longer formation times. Gas accretion simulations yield a rule of thumb that a rocky core must be at least 2$M_\oplus$ before it can acquire a volumetrically significant atmosphere from its parent nebula. In disks of low solid surface density, cores of such mass appear only after the gas disk has dissipated, and so remain purely rocky. Higher surface density disks breed massive cores more quickly, within the gas disk lifetime, and so produce gas-enveloped planets. We test model predictions against observations, using planet radius as an observational proxy for gas-to-rock content and host star metallicity as a proxy for disk solid surface density. Theory can explain the observation that metal-rich stars host predominantly gas-enveloped planets.

preprint2015arXiv

Breeding Super-Earths and Birthing Super-Puffs in Transitional Disks

The riddle posed by super-Earths (1-4$R_\oplus$, 2-20$M_\oplus$) is that they are not Jupiters: their core masses are large enough to trigger runaway gas accretion, yet somehow super-Earths accreted atmospheres that weigh only a few percent of their total mass. We show that this puzzle is solved if super-Earths formed late, as the last vestiges of their parent gas disks were about to clear. This scenario would seem to present fine-tuning problems, but we show that there are none. Ambient gas densities can span many (up to 9) orders of magnitude, and super-Earths can still robustly emerge after $\sim$0.1-1 Myr with percent-by-weight atmospheres. Super-Earth cores are naturally bred in gas-poor environments where gas dynamical friction has weakened sufficiently to allow constituent protocores to merge. So little gas is present at the time of core assembly that cores hardly migrate by disk torques: formation of super-Earths can be in situ. The picture --- that close-in super-Earths form in a gas-poor (but not gas-empty) inner disk, fed continuously by gas that bleeds inward from a more massive outer disk --- recalls the largely evacuated but still accreting inner cavities of transitional protoplanetary disks. We also address the inverse problem presented by super-puffs: an uncommon class of short-period planets seemingly too voluminous for their small masses (4-10$R_\oplus$, 2-6$M_\oplus$). Super-puffs easily acquire their thick atmospheres as dust-free, rapidly cooling worlds outside $\sim$1 AU where nebular gas is colder, less dense, and therefore less opaque. Unlike super-Earths which can form in situ, super-puffs migrated in to their current orbits; they are expected to form the outer links of mean-motion resonant chains, and to exhibit greater water content. We close by confronting observations and itemizing remaining questions.

preprint2015arXiv

Time Varying Dynamical Star Formation Rate

We present numerical evidence of dynamic star formation in which the accreted stellar mass grows superlinearly with time, roughly as $t^2$. We perform simulations of star formation in self-gravitating hydrodynamic and magneto-hydrodynamic turbulence that is continuously driven. By turning the self-gravity of the gas in the simulations on or off, we demonstrate that self-gravity is the dominant physical effect setting the mass accretion rate at early times before feedback effects take over, contrary to theories of turbulence-regulated star formation. We find that gravitational collapse steepens the density profile around stars, generating the power-law tail on what is otherwise a lognormal density probability distribution function. Furthermore, we find turbulent velocity profiles to flatten inside collapsing regions, altering the size-linewidth relation. This local flattening reflects enhancements of turbulent velocity on small scales, as verified by changes to the velocity power spectra. Our results indicate that gas self-gravity dynamically alters both density and velocity structures in clouds, giving rise to a time-varying star formation rate. We find that a substantial fraction of the gas that forms stars arrives via low density flows, as opposed to accreting through high density filaments.

preprint2015arXiv

To Cool is to Accrete: Analytic Scalings for Nebular Accretion of Planetary Atmospheres

Planets acquire atmospheres from their parent circumstellar disks. We derive a general analytic expression for how the atmospheric mass grows with time $t$, as a function of the underlying core mass $M_{\rm core}$ and nebular conditions, including the gas metallicity $Z$. Planets accrete as much gas as can cool: an atmosphere's doubling time is given by its Kelvin-Helmholtz time. Dusty atmospheres behave differently from atmospheres made dust-free by grain growth and sedimentation. The gas-to-core mass ratio (GCR) of a dusty atmosphere scales as GCR $\propto t^{0.4} M_{\rm core}^{1.7} Z^{-0.4} μ_{\rm rcb}^{3.4}$, where $μ_{\rm rcb} \propto 1/(1-Z)$ (for $Z$ not too close to 1) is the mean molecular weight at the innermost radiative-convective boundary. This scaling applies across all orbital distances and nebular conditions for dusty atmospheres; their radiative-convective boundaries, which regulate cooling, are not set by the external environment, but rather by the internal microphysics of dust sublimation, H$_2$ dissociation, and the formation of H$^-$. By contrast, dust-free atmospheres have their radiative boundaries at temperatures $T_{\rm rcb}$ close to nebular temperatures $T_{\rm out}$, and grow faster at larger orbital distances where cooler temperatures, and by extension lower opacities, prevail. At 0.1 AU in a gas-poor nebula, GCR $\propto t^{0.4} T_{\rm rcb}^{-1.9} M_{\rm core}^{1.6} Z^{-0.4} μ_{\rm rcb}^{3.3}$, while beyond 1 AU in a gas-rich nebula, GCR $\propto t^{0.4} T_{\rm rcb}^{-1.5} M_{\rm core}^1 Z^{-0.4}μ_{\rm rcb}^{2.2}$. We confirm our analytic scalings against detailed numerical models for objects ranging in mass from Mars (0.1 $M_\oplus$) to the most extreme super-Earths (10-20 $M_\oplus$), and explain why heating from planetesimal accretion cannot prevent the latter from undergoing runaway gas accretion.

preprint2014arXiv

Make Super-Earths, Not Jupiters: Accreting Nebular Gas onto Solid Cores at 0.1 AU and Beyond

Close-in super-Earths having radii 1--4 $R_\oplus$ may possess hydrogen atmospheres comprising a few percent by mass of their rocky cores. We determine the conditions under which such atmospheres can be accreted by cores from their parent circumstellar disks. Accretion from the nebula is problematic because it is too efficient: we find that 10-$M_\oplus$ cores embedded in solar metallicity disks tend to undergo runaway gas accretion and explode into Jupiters, irrespective of orbital location. The threat of runaway is especially dire at $\sim$0.1 AU, where solids may coagulate on timescales orders of magnitude shorter than gas clearing times; thus nascent atmospheres on close-in orbits are unlikely to be supported against collapse by planetesimal accretion. The time to runaway accretion is well approximated by the cooling time of the atmosphere's innermost convective zone, whose extent is controlled by where H$_2$ dissociates. Insofar as the temperatures characterizing H$_2$ dissociation are universal, timescales for core instability tend not to vary with orbital distance --- and to be alarmingly short for 10-$M_\oplus$ cores. Nevertheless, in the thicket of parameter space, we identify two scenarios, not mutually exclusive, that can reproduce the preponderance of percent-by-mass atmospheres for super-Earths at $\sim$0.1 AU, while still ensuring the formation of Jupiters at $\gtrsim 1$ AU. Scenario (a): planets form in disks with dust-to-gas ratios that range from $\sim$20$\times$ solar at 0.1 AU to $\sim$2$\times$ solar at 5 AU. Scenario (b): the final assembly of super-Earth cores from mergers of proto-cores --- a process that completes quickly at $\sim$0.1 AU once begun --- is delayed by gas dynamical friction until just before disk gas dissipates completely.