Showing 1–46 of 46 results for author: Lee, E J

Probing Dust and Gas Properties Using Ringed Disks

Authors: Eve J. Lee

Abstract: How rapidly a planet grows in mass and how far they may park from the host star depend sensitively on two non-dimensional parameters: Stokes number St and turbulent $α$. Yet, these parameters remain highly uncertain being difficult or impossible to measure directly. Here, we demonstrate how the ringed disks can be leveraged to obtain St and $α$ separately by constructing a simple toy model that co… ▽ More How rapidly a planet grows in mass and how far they may park from the host star depend sensitively on two non-dimensional parameters: Stokes number St and turbulent $α$. Yet, these parameters remain highly uncertain being difficult or impossible to measure directly. Here, we demonstrate how the ringed disks can be leveraged to obtain St and $α$ separately by constructing a simple toy model that combines dust radial equation of motion under aerodynamic drag and coupling to gas motion with the measured distribution of dust masses in Class 0/I disks. Focusing on known systems with well-resolved dust rings, we find that the range of St and $α$ that are consistent with the measured properties of the rings are small: $10^{-4} \lesssim {\rm St} \lesssim 10^{-2}$ and $10^{-5} \lesssim α\lesssim 10^{-3}$. These low St and $α$ ensure the observed rings are stable against clumping. Even in one marginal case where the formation of bound clumps is possible, further mass growth by pebble accretion is inhibited. Furthermore, the derived low $α$ is consistent with the nearly inviscid regime where Type I migration can be prematurely halted. Our analysis predicts minimal planet population beyond $\sim$10s of au where we observe dust rings and significantly more vigorous planet formation inside $\sim$10 AU, consistent with current exo-giant statistics. We close with discussions on the implications of our results on small planet statistics at large orbital distances. △ Less

Submitted 2 July, 2024; originally announced July 2024.

Comments: Accepted to ApJL

Giant Outer Transiting Exoplanet Mass (GOT `EM) Survey. V. Two Giant Planets in Kepler-511 but Only One Ran Away

Authors: Yayaati Chachan, Paul A. Dalba, Daniel P. Thorngren, Stephen R. Kane, Eve J. Lee, Edward W. Schwieterman, Howard Isaacson, Andrew W. Howard, Matthew J. Payne

Abstract: Systems hosting multiple giant planets are important laboratories for understanding planetary formation and migration processes. We present a nearly decade-long Doppler spectroscopy campaign from the Keck-I telescope to characterize the two transiting giant planets orbiting Kepler-511 on orbits of 27 days and 297 days. The radial velocity measurements yield precise masses for both planets, which w… ▽ More Systems hosting multiple giant planets are important laboratories for understanding planetary formation and migration processes. We present a nearly decade-long Doppler spectroscopy campaign from the Keck-I telescope to characterize the two transiting giant planets orbiting Kepler-511 on orbits of 27 days and 297 days. The radial velocity measurements yield precise masses for both planets, which we use to infer their bulk heavy element content. Both planets contain approximately 30 Earth masses of heavy elements, but their bulk metallicities (i.e., the ratio between metal mass and total mass) are drastically different ($0.86 \pm 0.04$ and $0.22 \pm 0.04$ respectively). Envelope mass loss cannot account for this difference due to the relatively large orbital distance and mass of the inner planet. We conclude that the outer planet underwent runaway gas accretion while the inner planet did not. This bifurcation in accretion histories is likely a result of the accretion of gas with very different metallicities by the two planets or the late formation of the inner planet from a merger of sub-Neptunes. Kepler-511 uniquely demonstrates how giant planet formation can produce strikingly different outcomes even for planets in the same system. △ Less

Submitted 17 May, 2024; originally announced May 2024.

Comments: submitted to AAS journals, comments welcome

Friends not Foes: Strong Correlation between Inner Super-Earths and Outer Gas Giants

Authors: Marta L. Bryan, Eve J. Lee

Abstract: The connection between outer gas giants and inner super-Earths reflects their formation and evolutionary histories. Past work exploring this link has suggested a tentative positive correlation between these two populations, but these studies have been limited by small sample sizes and in some cases sample biases. Here we take a new look at this connection with a sample of 184 super-Earth systems w… ▽ More The connection between outer gas giants and inner super-Earths reflects their formation and evolutionary histories. Past work exploring this link has suggested a tentative positive correlation between these two populations, but these studies have been limited by small sample sizes and in some cases sample biases. Here we take a new look at this connection with a sample of 184 super-Earth systems with publicly available radial velocity data and resolved outer gas giants. We calculate the frequency of outer gas giants (GG) in super-Earth (SE) systems, dividing our sample into metal-rich ([Fe/H] $>$ 0) and metal-poor ([Fe/H]$\leq$0) hosts. We find P(GG$|$SE, [Fe/H]$>$0) = 28.0$^{+4.9}_{-4.6}\%$ and P(GG$|$SE, [Fe/H]$\leq$0) = 4.5$^{+2.6}_{-1.9}\%$. Comparing these conditional occurrence rates to field giant occurrence rates from Rosenthal et al 2021, we show that there is a distinct positive correlation between inner super-Earths and outer gas giants for metal-rich host stars at the 2.7$σ$ level, but this correlation disappears for metal-poor systems. We further find that, around metal-rich stars, the GG/SE correlation enhances slightly for systems with giants that are more distant (beyond 3 AU), more eccentric ($e > 0.2$), and/or in multi-gas giant systems. Such trends disappear around metal-poor stars with the exception of systems of multiple giants in which we observe a tentative anti-correlation. Our findings highlight the critical role metallicity (disk solid budget) plays in shaping the overall planetary architecture. △ Less

Submitted 21 May, 2024; v1 submitted 13 March, 2024; originally announced March 2024.

Comments: accepted to ApJL

A Uniform Analysis of Debris Disks with the Gemini Planet Imager I: An Empirical Search for Perturbations from Planetary Companions in Polarized Light Images

Authors: Katie A. Crotts, Brenda C. Matthews, Gaspard Duchêne, Thomas M. Esposito, Ruobing Dong, Justin Hom, Rebecca Oppenheimer, Malena Rice, Schuyler G. Wolff, Christine H. Chen, Clarissa R. Do Ó, Paul Kalas, Briley L. Lewis, Alycia J. Weinberger, David J. Wilner, Mark Ammons, Pauline Arriaga, Robert J. De Rosa, John H. Debes, Michael P. Fitzgerald, Eileen C. Gonzales, Dean C. Hines, Sasha Hinkley, A. Meredith Hughes, Ludmilla Kolokolova , et al. (15 additional authors not shown)

Abstract: The Gemini Planet Imager (GPI) has excelled in imaging debris disks in the near-infrared. The GPI Exoplanet Survey (GPIES) imaged twenty-four debris disks in polarized $H$-band light, while other programs observed half of these disks in polarized $J$- and/or $K1$-bands. Using these data, we present a uniform analysis of the morphology of each disk to find asymmetries suggestive of perturbations, p… ▽ More The Gemini Planet Imager (GPI) has excelled in imaging debris disks in the near-infrared. The GPI Exoplanet Survey (GPIES) imaged twenty-four debris disks in polarized $H$-band light, while other programs observed half of these disks in polarized $J$- and/or $K1$-bands. Using these data, we present a uniform analysis of the morphology of each disk to find asymmetries suggestive of perturbations, particularly those due to planet-disk interactions. The multi-wavelength surface brightness, the disk color and geometry permit identification of any asymmetries such as warps or disk offsets from the central star. We find that nineteen of the disks in this sample exhibit asymmetries in surface brightness, disk color, disk geometry, or a combination of the three, suggesting that for this sample, perturbations, as seen in scattered light, are common. The relationship between these perturbations and potential planets in the system are discussed. We also explore correlations among stellar temperatures, ages, disk properties, and observed perturbations. We find significant trends between the vertical aspect ratio and the stellar temperature, disk radial extent, and the dust grain size distribution power-law, $q$. We also confirm a trend between the disk color and stellar effective temperature, where the disk becomes increasingly red/neutral with increasing temperature. Such results have important implications on the evolution of debris disk systems around stars of various spectral types. △ Less

Submitted 24 November, 2023; originally announced November 2023.

Comments: 46 pages, 20 figures, 6 tables, accepted for publication in ApJ

The Not-So Dramatic Effect of Advective Flows on Gas Accretion

Authors: Vincent Savignac, Eve J. Lee

Abstract: Super-Earths and mini-Neptunes are the most common types of exoplanets discovered, yet the physics of their formation are still debated. Standard core accretion models in gas-rich environment find that typical mini-Neptune mass planets would blow up into Jupiters before the underlying disk gas dissipates away. The injection of entropy from the protoplanetary disk into forming gaseous envelopes has… ▽ More Super-Earths and mini-Neptunes are the most common types of exoplanets discovered, yet the physics of their formation are still debated. Standard core accretion models in gas-rich environment find that typical mini-Neptune mass planets would blow up into Jupiters before the underlying disk gas dissipates away. The injection of entropy from the protoplanetary disk into forming gaseous envelopes has recently been put forward as a mechanism to delay this runaway accretion, specifically at short orbital distances. Here, we reevaluate this line of reasoning by incorporating recycling flows of gas into a numerical one-dimensional thermodynamic model with more realistic equation of state and opacities and the thermal state of the advective flow. At 0.1 AU, we find that advective flows are only able to produce mini-Neptunes if they can penetrate below ~0.25 of the planet's gravitational sphere of influence. Otherwise, the gas-to-core mass ratio (GCR) reaches above ~10% which is too large to explain the measured properties of mini-Neptunes, necessitating other gas-limiting processes such as late-time core assembly. The effect of entropy advection on gas accretion weakens even further beyond 0.1 AU. We present an updated scaling relation between GCR and the penetration depth of the advective flows which varies non-trivially with orbital distances, core masses and dusty vs. dust-free opacity. We further demonstrate how measurements of planet mass distribution beyond ~1 AU using future instruments such as the Nancy Grace Roman Space Telescope could be used to disambiguate between different formation conditions of gas-poor planets. △ Less

Submitted 9 July, 2024; v1 submitted 9 October, 2023; originally announced October 2023.

Comments: 18 pages, 9 figures. Accepted to ApJ

Small Planets Around Cool Dwarfs: Enhanced Formation Efficiency of Super-Earths around M dwarfs

Authors: Yayaati Chachan, Eve J. Lee

Abstract: Current measurements of planet population as a function of stellar mass show three seemingly contradictory signatures: close-in super-Earths are more prevalent around M dwarfs than FGK dwarfs; inner super-Earths are correlated with outer giants; and outer giants are less common around M dwarfs than FGK dwarfs. Here, we build a simple framework that combines the theory of pebble accretion with the… ▽ More Current measurements of planet population as a function of stellar mass show three seemingly contradictory signatures: close-in super-Earths are more prevalent around M dwarfs than FGK dwarfs; inner super-Earths are correlated with outer giants; and outer giants are less common around M dwarfs than FGK dwarfs. Here, we build a simple framework that combines the theory of pebble accretion with the measurements of dust masses in protoplanetary disks to reconcile all three observations. First, we show that cooler stars are more efficient at converting pebbles into planetary cores at short orbital periods. Second, when disks are massive enough to nucleate a heavy core at 5 AU, more than enough dust can drift in to assemble inner planets, establishing the correlation between inner planets and outer giants. Finally, while stars of varying masses are similarly capable of converting pebbles into cores at long orbital periods, hotter stars are much more likely to harbor more massive dust disks so that the giant planet occurrence rate rises around hotter stars. Our results are valid over a wide range of parameter space for a disk accretion rate that follows $\dot{M}_\star \sim 10^{-8}\,M_\odot\,{\rm yr}^{-1}(M_\star/M_\odot)^2$. We predict a decline in mini-Neptune population (but not necessarily terrestrial planets) around stars lighter than $\sim 0.3-0.5 M_\odot$. Cold giants ($\gtrsim$5 AU), if they exist, should remain correlated with inner planets even around lower mass stars. △ Less

Submitted 28 June, 2023; v1 submitted 1 May, 2023; originally announced May 2023.

Comments: accepted to ApJL

Removal of Hot Saturns in Mass-Radius Plane by Runaway Mass Loss

Authors: Daniel P. Thorngren, Eve J. Lee, Eric D. Lopez

Abstract: The hot Saturn population exhibits a boundary in mass-radius space, such that no planets are observed at a density less than $\sim$0.1 g cm$^{-3}$. Yet, planet interior structure models can readily construct such objects as the natural result of radius inflation. Here, we investigate the role XUV-driven mass-loss plays in sculpting the density boundary by constructing interior structure models tha… ▽ More The hot Saturn population exhibits a boundary in mass-radius space, such that no planets are observed at a density less than $\sim$0.1 g cm$^{-3}$. Yet, planet interior structure models can readily construct such objects as the natural result of radius inflation. Here, we investigate the role XUV-driven mass-loss plays in sculpting the density boundary by constructing interior structure models that include radius inflation, photoevaporative mass loss and a simple prescription of Roche lobe overflow. We demonstrate that planets puffier than $\sim$0.1 g cm$^{-3}$ experience a runaway mass loss caused by adiabatic radius expansion as the gas layer is stripped away, providing a good explanation of the observed edge in mass-radius space. The process is also visible in the radius-period and mass-period spaces, though smaller, high-bulk-metallicity planets can still survive at short periods, preserving a partial record of the population distribution at formation. △ Less

Submitted 8 March, 2023; v1 submitted 21 November, 2022; originally announced November 2022.

Comments: 11 pages, 5 figures, Accepted to ApJ Letters. This revision includes the effect of orbital evolution resulting from Roche lobe overflow; the overall conclusions are the same

Establishing Dust Rings and Forming Planets Within Them

Authors: Eve J. Lee, J. R. Fuentes, Philip F. Hopkins

Abstract: Radio images of protoplanetary disks demonstrate that dust grains tend to organize themselves into rings. These rings may be a consequence of dust trapping within gas pressure maxima wherein the local high dust-to-gas ratio is expected to trigger the formation of planetesimals and eventually planets. We revisit the behavior of dust near gas pressure perturbations enforced by a planet in two-dimens… ▽ More Radio images of protoplanetary disks demonstrate that dust grains tend to organize themselves into rings. These rings may be a consequence of dust trapping within gas pressure maxima wherein the local high dust-to-gas ratio is expected to trigger the formation of planetesimals and eventually planets. We revisit the behavior of dust near gas pressure perturbations enforced by a planet in two-dimensional, shearing box simulations. While dust grains collect into generally long-lived rings, particles with small Stokes parameter $τ_s < 0.1$ tend to advect out of the ring within a few drift timescales. Scaled to the properties of ALMA disks, we find that rings composed of larger particles ($τ_s \geq 0.1$) can nucleate a dust clump massive enough to trigger pebble accretion which proceeds to ingest the entire dust ring well within $\sim$1 Myr. To ensure the survival of the dust rings, we favor a non-planetary origin and typical grain size $τ_s \lesssim 0.05$--0.1. Planet-driven rings may still be possible but if so we would expect the orbital distance of the dust rings to be larger for older systems. △ Less

Submitted 24 August, 2022; v1 submitted 2 June, 2022; originally announced June 2022.

Comments: To appear in ApJ. Paper refocused to forcing by a planet

Creating the Radius Gap without Mass Loss

Authors: Eve J. Lee, Amalia Karalis, Daniel P. Thorngren

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

Submitted 23 October, 2022; v1 submitted 24 January, 2022; originally announced January 2022.

Comments: ApJ in press

Sculpting the sub-Saturn Occurrence Rate via Atmospheric Mass Loss

Authors: Tim Hallatt, Eve J. Lee

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

Submitted 21 October, 2021; v1 submitted 26 May, 2021; originally announced May 2021.

Comments: Accepted to ApJ

Radial Gradients in Dust-to-Gas Ratio Lead to Preferred Region for Giant Planet Formation

Authors: Yayaati Chachan, Eve J. Lee, Heather A. Knutson

Abstract: The Rosseland mean opacity of dust in protoplanetary disks is often calculated assuming the interstellar medium (ISM) size distribution and a constant dust-to-gas ratio. However, the dust size distribution and the dust-to-gas ratio in protoplanetary disks are distinct from those of the ISM. Here, we use simple dust evolution models that incorporate grain growth and transport to calculate the time… ▽ More The Rosseland mean opacity of dust in protoplanetary disks is often calculated assuming the interstellar medium (ISM) size distribution and a constant dust-to-gas ratio. However, the dust size distribution and the dust-to-gas ratio in protoplanetary disks are distinct from those of the ISM. Here, we use simple dust evolution models that incorporate grain growth and transport to calculate the time evolution of mean opacity of dust grains as a function of distance from the star. Dust dynamics and size distribution are sensitive to the assumed value of the turbulence strength $α_{\rm t}$ and the velocity at which grains fragment $v_{\rm frag}$. For moderate-to-low turbulence strengths of $α_{\rm t} \lesssim 10^{-3}$ and substantial differences in $v_{\rm frag}$ for icy and ice-free grains, we find a spatially non-uniform dust-to-gas ratio and grain size distribution that deviate significantly from the ISM values, in agreement with previous studies. The effect of a non-uniform dust-to-gas ratio on the Rosseland mean opacity dominates over that of the size distribution. This spatially varying -- that is non-monotonic -- dust-to-gas ratio creates a region in the protoplanetary disk that is optimal for producing hydrogen-rich planets, potentially explaining the apparent peak in gas giant planet occurrence rate at intermediate distances. The enhanced dust-to-gas ratio within the ice line also suppresses gas accretion rates onto sub-Neptune cores, thus stifling their tendency to undergo runaway gas accretion within disk lifetimes. Finally, our work corroborates the idea that low mass cores with large primordial gaseous envelopes (`super-puffs') originate beyond the ice line. △ Less

Submitted 30 June, 2021; v1 submitted 25 January, 2021; originally announced January 2021.

Comments: accepted to ApJ, 21 pages, 14 figures; added runaway gas accretion prescription and a discussion on accretion heating

Journal ref: ApJ, 919, 63 (2021)

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

Authors: Caroline Piaulet, Björn Benneke, Ryan A. Rubenzahl, Andrew W. Howard, Eve J. Lee, Daniel Thorngren, Ruth Angus, Merrin Peterson, Joshua E. Schlieder, Michael Werner, Laura Kreidberg, Tareq Jaouni, Ian J. M. Crossfield, David R. Ciardi, Erik A. Petigura, John Livingston, Courtney D. Dressing, Benjamin J. Fulton, Charles Beichman, Jessie L. Christiansen, Varoujan Gorjian, Kevin K. Hardegree-Ullman, Jessica Krick, Evan Sinukoff

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

Submitted 26 November, 2020; originally announced November 2020.

Comments: Accepted for publication in AJ, 16 pages, 7 figures

Can Large-Scale Migration Explain the Giant Planet Occurrence Rate?

Authors: Tim Hallatt, Eve J Lee

Abstract: The giant planet occurrence rate rises with orbital period out to at least $\sim$300 days. Large-scale planetary migration through the disk has long been suspected to be the physical origin of this feature, as the timescale of standard Type I migration in a standard solar nebula is longer farther from the star. These calculations also find that typical Jupiter-bearing cores shuttle towards the dis… ▽ More The giant planet occurrence rate rises with orbital period out to at least $\sim$300 days. Large-scale planetary migration through the disk has long been suspected to be the physical origin of this feature, as the timescale of standard Type I migration in a standard solar nebula is longer farther from the star. These calculations also find that typical Jupiter-bearing cores shuttle towards the disk inner edge on timescales orders of magnitude shorter than the gas disk lifetime. The presence of gas giants at myriad distances requires mechanisms to slow large-scale migration. We revisit the migration paradigm by deriving model occurrence rate profiles from migration of cores, mass growth by gas accretion, and planetary gap opening. We show explicitly that the former two processes occur in tandem. Radial transport of planets can slow down significantly once deep gaps are carved out by their interaction with disk gas. Disks are more easily perturbed closer to the star, so accounting for gap opening flattens the final orbital period distribution. To recover the observed rise in occurrence rate, gas giants need to be more massive farther out, which is naturally achieved if their envelopes are dust-free. The range of mass gradients we find to reconcile the observed occurrence rate of Jupiters is too abrupt to account for the mass-period distribution of low-eccentricity giant planets, challenging disk migration as the dominant origin channel of hot and warm Jupiters. Future efforts in characterizing the unbiased mass distribution will place stronger constraints on predictions from migration theory. △ Less

Submitted 15 October, 2020; originally announced October 2020.

Comments: Accepted to ApJ

A Featureless Infrared Transmission Spectrum for the Super-Puff Planet Kepler-79d

Authors: Yayaati Chachan, Daniel Jontof-Hutter, Heather A. Knutson, Danica Adams, Peter Gao, Björn Benneke, Zachory Berta-Thompson, Fei Dai, Drake Deming, Eric Ford, Eve J. Lee, Jessica E. Libby-Roberts, Nikku Madhusudhan, Hannah R. Wakeford, Ian Wong

Abstract: Extremely low density planets ('super-puffs') are a small but intriguing subset of the transiting planet population. With masses in the super-Earth range ($1-10$ M$_{\oplus}$) and radii akin to those of giant planets ($>4$ R$_{\oplus}$), their large envelopes may have been accreted beyond the water snow line and many appear to be susceptible to catastrophic mass loss. Both the presence of water an… ▽ More Extremely low density planets ('super-puffs') are a small but intriguing subset of the transiting planet population. With masses in the super-Earth range ($1-10$ M$_{\oplus}$) and radii akin to those of giant planets ($>4$ R$_{\oplus}$), their large envelopes may have been accreted beyond the water snow line and many appear to be susceptible to catastrophic mass loss. Both the presence of water and the importance of mass loss can be explored using transmission spectroscopy. Here, we present new HST WFC3 spectroscopy and updated Kepler transit depth measurements for the super-puff Kepler-79d. We do not detect any molecular absorption features in the $1.1-1.7$ $μ$m WFC3 bandpass and the combination of Kepler and WFC3 data are consistent with a flat line model, indicating the presence of aerosols in the atmosphere. We compare the shape of Kepler-79d's transmission spectrum to predictions from a microphysical haze model that incorporates an outward particle flux due to ongoing mass loss. We find that photochemical hazes offer an attractive explanation for the observed properties of super-puffs like Kepler-79d, as they simultaneously render the near-infrared spectrum featureless and reduce the inferred envelope mass loss rate by moving the measured radius (optical depth unity surface during transit) to lower pressures. We revisit the broader question of mass loss rates for super-puffs and find that the age estimates and mass loss rates for the majority of super-puffs can be reconciled if hazes move the photosphere from the typically assumed pressure of $\sim 10$ mbar to $\sim 10 \; μ$bar. △ Less

Submitted 19 August, 2020; v1 submitted 12 August, 2020; originally announced August 2020.

Comments: Awaiting publication in AJ. Small updates in Table 6 and Fig 11. Table 3 and 6 will be provided in MRT format upon publication

Journal ref: AJ, 160, 201 (2020)

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

Authors: Eve J. Lee, Nicholas J. Connors

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

Submitted 22 December, 2020; v1 submitted 3 August, 2020; originally announced August 2020.

Comments: Accepted to ApJ. New Figures 5 & 6 compare directly the radius gap carved out purely from late-time accretion with that from photoevaporation

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

Authors: Eve J. Lee, Philip F. Hopkins

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

Submitted 26 March, 2020; v1 submitted 6 March, 2020; originally announced March 2020.

Comments: MNRAS Letters in press

Forming Diverse Super-Earth Systems in Situ

Authors: Mariah G. MacDonald, Rebekah I. Dawson, Sarah J. Morrison, Eve J. Lee, Arjun Khandelwal

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

Submitted 17 January, 2020; originally announced January 2020.

Comments: 18 pages, 9 figures. Submitted to ApJ; revised in response to referee report

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

Authors: Joseph E. Rodriguez, Andrew Vanderburg, Sebastian Zieba, Laura Kreidberg, Caroline V. Morley, Jason D. Eastman, Stephen R. Kane, Alton Spencer, Samuel N. Quinn, Ryan Cloutier, Chelsea X. Huang, Karen A. Collins, Andrew W. Mann, Emily Gilbert, Joshua E. Schlieder, Elisa V. Quintana, Thomas Barclay, Gabrielle Suissa, Ravi kumar Kopparapu, Courtney D. Dressing, George R. Ricker, Roland K. Vanderspek, David W. Latham, Sara Seager, Joshua N. Winn , et al. (21 additional authors not shown)

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

Submitted 10 June, 2020; v1 submitted 3 January, 2020; originally announced January 2020.

Comments: 14 Pages, 5 Figures, 2 Tables, Accepted to AJ

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

Authors: Emily A. Gilbert, Thomas Barclay, Joshua E. Schlieder, Elisa V. Quintana, Benjamin J. Hord, Veselin B. Kostov, Eric D. Lopez, Jason F. Rowe, Kelsey Hoffman, Lucianne M. Walkowicz, Michele L. Silverstein, Joseph E. Rodriguez, Andrew Vanderburg, Gabrielle Suissa, Vladimir S. Airapetian, Matthew S. Clement, Sean N. Raymond, Andrew W. Mann, Ethan Kruse, Jack J. Lissauer, Knicole D. Colón, Ravi kumar Kopparapu, Laura Kreidberg, Sebastian Zieba, Karen A. Collins , et al. (70 additional authors not shown)

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

Submitted 10 July, 2020; v1 submitted 3 January, 2020; originally announced January 2020.

Comments: Accepted for publication in AJ (30 pages, 13 figures, 4 tables, 2 appendices)

Eruptions and Flaring Activity in Emerging Quadrupolar Regions

Authors: P. Syntelis, E. J. Lee, C. W. Fairbairn, V. Archontis, A. W. Hood

Abstract: Context. Some of the most dynamic active regions are associated with complex photospheric magnetic configurations such as quadrupolar regions, and especially ones with a $δ$-spot configuration and a strong Polarity Inversion Line (PIL). Aims. We study the formation and eruption of magnetic flux ropes in quadrupolar regions. Methods. 3D MHD simulations of the partial emergence of a highly twisted f… ▽ More Context. Some of the most dynamic active regions are associated with complex photospheric magnetic configurations such as quadrupolar regions, and especially ones with a $δ$-spot configuration and a strong Polarity Inversion Line (PIL). Aims. We study the formation and eruption of magnetic flux ropes in quadrupolar regions. Methods. 3D MHD simulations of the partial emergence of a highly twisted flux tube from the solar interior into a non-magnetized stratified atmosphere. Results. We emerge two $Ω$-shaped loops forming a quadrupolar region. The emerging flux forms two initially separated bipoles, that later come in contact creating a $δ$-spot central region. Above the two bipoles, two magnetic lobes expand and interact through a series of current sheets at the interface between them. Two recurrent confined eruptions are produced. In both cases, reconnection between sheared low-lying field lines forms a flux rope. Reconnection between the two lobes high in the atmosphere forms retracting field lines that push against the flux rope, creating a current sheet between them. It also forms a third magnetic lobe between the two emerged ones, that later acts as a strapping field. The flux rope eruptions are triggered when the reconnection between the flux ropes and the field above them becomes efficient enough to remove the tension of the overlying field. These reconnection events occur internally in the system. The first erupting flux rope almost fully reconnects with the overlying field. The second eruption is confined by the overlying field. During the confined eruption, the flux rope is enhanced in size, flux and twist, similar to confined-flare-to-flux-rope observations. Proxies of the emission reveal the two erupting filaments channels. A flare arcade is formed only in the second eruption due to the longer-lasting and more efficient reconnection below the flux rope. △ Less

Submitted 3 September, 2019; originally announced September 2019.

Comments: 15 pages, 19 figures, to appear in Astronomy and Astrophysics

Journal ref: A&A 630, A134 (2019)

A multi-wavelength study of the debris disc around 49 Cet

Authors: Nicole Pawellek, Attila Moór, Julien Milli, Ágnes Kóspál, Johan Olofsson, Péter Ábrahám, Miriam Keppler, Quentin Kral, Adriana Pohl, Jean-Charles Augereau, Anthony Boccaletti, Gaël Chauvin, Élodie Choquet, Natalia Engler, Thomas Henning, Maud Langlois, Eve J. Lee, François Ménard, Philippe Thébault, Alice Zurlo

Abstract: In a multi-wavelength study of thermal emission and scattered light images we analyse the dust properties and structure of the debris disc around the A1-type main sequence star 49~Cet. As a basis for this study, we present new scattered light images of the debris disc known to possess both a high amount of dust and gas. The outer region of the disc is revealed in former coronagraphic H-band and ou… ▽ More In a multi-wavelength study of thermal emission and scattered light images we analyse the dust properties and structure of the debris disc around the A1-type main sequence star 49~Cet. As a basis for this study, we present new scattered light images of the debris disc known to possess both a high amount of dust and gas. The outer region of the disc is revealed in former coronagraphic H-band and our new Y-band images from the Very Large Telescope SPHERE instrument. We use the knowledge of the disc's radial extent inferred from ALMA observations and the grain size distribution found by SED fitting to generate semi-dynamical dust models of the disc. We compare the models to scattered light and thermal emission data and find that a disc with a maximum of the surface density at 110~au and shallow edges can describe both thermal emission and scattered light observations. This suggests that grains close to the blow-out limit and large grains stem from the same planetesimal population and are mainly influenced by radiation pressure. The influence of inwards transport processes could not be analysed in this study. △ Less

Submitted 15 July, 2019; originally announced July 2019.

Comments: Accepted by MNRAS

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

Authors: Shreyas Vissapragada, Daniel Jontof-Hutter, Avi Shporer, Heather A. Knutson, Leo Liu, Daniel Thorngren, Eve J. Lee, Yayaati Chachan, Dimitri Mawet, Maxwell A. Millar-Blanchaer, Ricky Nilsson, Samaporn Tinyanont, Gautam Vasisht, Jason Wright

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

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

Comments: 33 pages, 18 figures, accepted by AJ

The Boundary between Gas-rich and Gas-poor Planets

Authors: Eve J. Lee

Abstract: Sub-Saturns straddle the boundary between gas-rich Jupiters and gas-poor super-Earths/sub-Neptunes. Their large radii (4--8$R_\oplus$) suggest that their gas-to-core mass ratios range $\sim$0.1--1.0. With their envelopes as massive as their cores, sub-Saturns are just on the verge of runaway gas accretion; they are expected to be significantly less populous than gas giants. Yet, the observed occur… ▽ More Sub-Saturns straddle the boundary between gas-rich Jupiters and gas-poor super-Earths/sub-Neptunes. Their large radii (4--8$R_\oplus$) suggest that their gas-to-core mass ratios range $\sim$0.1--1.0. With their envelopes as massive as their cores, sub-Saturns are just on the verge of runaway gas accretion; they are expected to be significantly less populous than gas giants. Yet, the observed occurrence rates of sub-Saturns and Jupiters are comparable within $\sim$100 days. We show that in these inner regions of planetary systems, the growth of sub-Saturns/Jupiters is ultimately limited by local and global hydrodynamic flows---runaway accretion terminates and the formation of gas giants is suppressed. We derive a simple analytic formula for the local hydrodynamic accretion rate---an expression that has been previously reported only as an empirical fit to numerical simulations. Evolving simultaneously the background disk gas and the gas accretion onto planetary cores, we find that both the ubiquity of super-Earths/sub-Neptunes and the rarity of gas-rich planets are best explained if an underlying core mass distribution is peaked at $\sim$4.3$M_\oplus$. Within a finite disk lifetime $\sim$10 Myrs, massive cores ($\gtrsim 10M_\oplus$) can become either gas-poor or gas-rich depending on when they assemble but smaller cores ($\lesssim 10M_\oplus$) can only become gas-poor. This wider range of possible outcomes afforded by more massive cores may explain why metal-rich stars harbor a more diverse set of planets. △ Less

Submitted 3 June, 2019; v1 submitted 23 April, 2019; originally announced April 2019.

Comments: ApJ in press. Figure 6 now compares model gas-to-core mass ratio vs. core mass to that inferred from observed exoplanets

The Gemini Planet Imager Exoplanet Survey: Giant Planet and Brown Dwarf Demographics From 10-100 AU

Authors: Eric L. Nielsen, Robert J. De Rosa, Bruce Macintosh, Jason J. Wang, Jean-Baptiste Ruffio, Eugene Chiang, Mark S. Marley, Didier Saumon, Dmitry Savransky, S. Mark Ammons, Vanessa P. Bailey, Travis Barman, Celia Blain, Joanna Bulger, Jeffrey Chilcote, Tara Cotten, Ian Czekala, Rene Doyon, Gaspard Duchene, Thomas M. Esposito, Daniel Fabrycky, Michael P. Fitzgerald, Katherine B. Follette, Jonathan J. Fortney, Benjamin L. Gerard , et al. (40 additional authors not shown)

Abstract: We present a statistical analysis of the first 300 stars observed by the Gemini Planet Imager Exoplanet Survey (GPIES). This subsample includes six detected planets and three brown dwarfs; from these detections and our contrast curves we infer the underlying distributions of substellar companions with respect to their mass, semi-major axis, and host stellar mass. We uncover a strong correlation be… ▽ More We present a statistical analysis of the first 300 stars observed by the Gemini Planet Imager Exoplanet Survey (GPIES). This subsample includes six detected planets and three brown dwarfs; from these detections and our contrast curves we infer the underlying distributions of substellar companions with respect to their mass, semi-major axis, and host stellar mass. We uncover a strong correlation between planet occurrence rate and host star mass, with stars M $>$ 1.5 $M_\odot$ more likely to host planets with masses between 2-13 M$_{\rm Jup}$ and semi-major axes of 3-100 au at 99.92% confidence. We fit a double power-law model in planet mass (m) and semi-major axis (a) for planet populations around high-mass stars (M $>$ 1.5M$_\odot$) of the form $\frac{d^2 N}{dm da} \propto m^αa^β$, finding $α$ = -2.4 $\pm$ 0.8 and $β$ = -2.0 $\pm$ 0.5, and an integrated occurrence rate of $9^{+5}_{-4}$% between 5-13 M$_{\rm Jup}$ and 10-100 au. A significantly lower occurrence rate is obtained for brown dwarfs around all stars, with 0.8$^{+0.8}_{-0.5}$% of stars hosting a brown dwarf companion between 13-80 M$_{\rm Jup}$ and 10-100 au. Brown dwarfs also appear to be distributed differently in mass and semi-major axis compared to giant planets; whereas giant planets follow a bottom-heavy mass distribution and favor smaller semi-major axes, brown dwarfs exhibit just the opposite behaviors. Comparing to studies of short-period giant planets from the RV method, our results are consistent with a peak in occurrence of giant planets between ~1-10 au. We discuss how these trends, including the preference of giant planets for high-mass host stars, point to formation of giant planets by core/pebble accretion, and formation of brown dwarfs by gravitational instability. △ Less

Submitted 10 April, 2019; originally announced April 2019.

Comments: 52 pages, 18 figures. AJ in press

Tracing the Origins and Evolution of Small Planets using Their Orbital Obliquities

Authors: Marshall C. Johnson, George Zhou, Brett C. Addison, David R. Ciardi, Diana Dragomir, Yasuhiro Hasegawa, Eve J. Lee, Songhu Wang, Lauren Weiss

Abstract: We recommend an intensive effort to survey and understand the obliquity distribution of small close-in extrasolar planets over the coming decade. The orbital obliquities of exoplanets--i.e., the relative orientation between the planetary orbit and the stellar rotation--is a key tracer of how planets form and migrate. While the orbital obliquities of smaller planets are poorly explored today, a new… ▽ More We recommend an intensive effort to survey and understand the obliquity distribution of small close-in extrasolar planets over the coming decade. The orbital obliquities of exoplanets--i.e., the relative orientation between the planetary orbit and the stellar rotation--is a key tracer of how planets form and migrate. While the orbital obliquities of smaller planets are poorly explored today, a new generation of facilities coming online over the next decade will make such observations possible en masse. Transit spectroscopic observations with the extremely large telescopes will enable us to measure the orbital obliquities of planets as small as $\sim2R_{\oplus}$ around a wide variety of stars, opening a window into the orbital properties of the most common types of planets. This effort will directly contribute to understanding the formation and evolution of planetary systems, a key objective of the National Academy of Sciences' Exoplanet Science Strategies report. △ Less

Submitted 11 March, 2019; originally announced March 2019.

Comments: Submitted to the Astro2020 call for science white papers. 7 pages, 2 figures

Deep exploration of $ε$ Eridani with Keck Ms-band vortex coronagraphy and radial velocities: mass and orbital parameters of the giant exoplanet

Authors: Dimitri Mawet, Lea Hirsch, Eve J. Lee, Jean-Baptiste Ruffio, Michael Bottom, Benjamin J. Fulton, Olivier Absil, Charles Beichman, Brendan Bowler, Marta Bryan, Elodie Choquet, David Ciardi, Valentin Christiaens, Denis Defrère, Carlos Alberto Gomez Gonzalez, Andrew W. Howard, Elsa Huby, Howard Isaacson, Rebecca Jensen-Clem, Molly Kosiarek, Geoff Marcy, Tiffany Meshkat, Erik Petigura, Maddalena Reggiani, Garreth Ruane , et al. (5 additional authors not shown)

Abstract: We present the most sensitive direct imaging and radial velocity (RV) exploration of $ε$ Eridani to date. $ε$ Eridani is an adolescent planetary system, reminiscent of the early Solar system. It is surrounded by a prominent and complex debris disk which is likely stirred by one or several gas giant exoplanets. The discovery of the RV signature of a giant exoplanet was announced 15 years ago, but h… ▽ More We present the most sensitive direct imaging and radial velocity (RV) exploration of $ε$ Eridani to date. $ε$ Eridani is an adolescent planetary system, reminiscent of the early Solar system. It is surrounded by a prominent and complex debris disk which is likely stirred by one or several gas giant exoplanets. The discovery of the RV signature of a giant exoplanet was announced 15 years ago, but has met with scrutiny due to possible confusion with stellar noise. We confirm the planet with a new compilation and analysis of precise RV data spanning 30 years, and combine it with upper limits from our direct imaging search, the most sensitive ever performed. The deep images were taken in the Ms band (4.7$μ$m) with the vortex coronagraph recently installed in W.M. Keck Observatory's infrared camera NIRC2, which opens a sensitive window for planet searches around nearby adolescent systems. The RV data and direct imaging upper limit maps were combined in an innovative joint Bayesian analysis, providing new constraints on the mass and orbital parameters of the elusive planet. $ε$ Eridani b has a mass of $0.78^{+0.38}_{-0.12}$ $M_{Jup}$ and is orbiting $ε$ Eridani at about $3.48\pm 0.02$ AU with a period of $7.37 \pm 0.07$ years. The eccentricity of $ε$ Eridani b's orbit is $0.07^{+0.06}_{-0.05}$, an order of magnitude smaller than early estimates and consistent with a circular orbit. We discuss our findings from the standpoint of planet-disk interactions and prospects for future detection and characterization with the James Webb Space Telescope. △ Less

Submitted 29 October, 2018; v1 submitted 8 October, 2018; originally announced October 2018.

Comments: Accepted to AJ

On The Nature of Variations in the Measured Star Formation Efficiency of Molecular Clouds

Authors: Michael Y. Grudić, Philip F. Hopkins, Eve J. Lee, Norman Murray, Claude-André Faucher-Giguère, L. Clifton Johnson

Abstract: Measurements of the star formation efficiency (SFE) of giant molecular clouds (GMCs) in the Milky Way generally show a large scatter, which could be intrinsic or observational. We use magnetohydrodynamic simulations of GMCs (including feedback) to forward-model the relationship between the true GMC SFE and observational proxies. We show that individual GMCs trace broad ranges of observed SFE throu… ▽ More Measurements of the star formation efficiency (SFE) of giant molecular clouds (GMCs) in the Milky Way generally show a large scatter, which could be intrinsic or observational. We use magnetohydrodynamic simulations of GMCs (including feedback) to forward-model the relationship between the true GMC SFE and observational proxies. We show that individual GMCs trace broad ranges of observed SFE throughout collapse, star formation, and disruption. Low measured SFEs (<<1%) are "real" but correspond to early stages, the true "per-freefall" SFE where most stars actually form can be much larger. Very high (>>10%) values are often artificially enhanced by rapid gas dispersal. Simulations including stellar feedback reproduce observed GMC-scale SFEs, but simulations without feedback produce 20x larger SFEs. Radiative feedback dominates among mechanisms simulated. An anticorrelation of SFE with cloud mass is shown to be an observational artifact. We also explore individual dense "clumps" within GMCs and show that (with feedback) their bulk properties agree well with observations. Predicted SFEs within the dense clumps are ~2x larger than observed, possibly indicating physics other than feedback from massive (main sequence) stars is needed to regulate their collapse. △ Less

Submitted 25 September, 2018; v1 submitted 21 September, 2018; originally announced September 2018.

Comments: Fixed typo in the arXiv abstract

An Excess of Jupiter Analogs in Super-Earth Systems

Authors: Marta L. Bryan, Heather A. Knutson, Eve J. Lee, BJ Fulton, Konstantin Batygin, Henry Ngo, Tiffany Meshkat

Abstract: We use radial velocity observations to search for long-period gas giant companions in systems hosting inner super-Earth (1-4 R_Earth, 1-10 M_Earth) planets to constrain formation and migration scenarios for this population. We consistently re-fit published RV datasets for 65 stars and find 9 systems with statistically significant trends indicating the presence of an outer companion. We combine the… ▽ More We use radial velocity observations to search for long-period gas giant companions in systems hosting inner super-Earth (1-4 R_Earth, 1-10 M_Earth) planets to constrain formation and migration scenarios for this population. We consistently re-fit published RV datasets for 65 stars and find 9 systems with statistically significant trends indicating the presence of an outer companion. We combine these RV data with AO images to constrain the masses and semi-major axes of these companions. We quantify our sensitivity to the presence of long-period companions by fitting the sample with a power law distribution and find an occurrence rate of 39+/-7% for companions 0.5-20 M_Jup and 1-20 AU. Half of our systems were discovered by the transit method and half were discovered by the RV method. While differences in RV baselines and number of data points between the two samples lead to different sensitivities to distant companions, we find that occurrence rates of gas giant companions in each sample are consistent at the 0.5$σ$ level. We compare the frequency of Jupiter analogs in these systems to the equivalent rate from field star surveys and find that Jupiter analogs are more common around stars hosting super-Earths. We conclude that the presence of outer gas giants does not suppress the formation of inner super-Earths, and that these two populations of planets instead appear to be correlated. We also find that the stellar metallicities of systems with gas giant companions are higher than those without companions, in agreement with the well-established metallicity correlation from RV surveys of field stars. △ Less

Submitted 23 January, 2019; v1 submitted 22 June, 2018; originally announced June 2018.

Comments: published in AJ

A Balanced Budget View on Forming Giant Planets by Pebble Accretion

Authors: Jonathan W. Lin, Eve J. Lee, Eugene Chiang

Abstract: Pebble accretion refers to the assembly of rocky planet cores from particles whose velocity dispersions are damped by drag from circumstellar disc gas. Accretion cross-sections can approach maximal Hill-sphere scales for particles whose Stokes numbers approach unity. While fast, pebble accretion is also lossy. Gas drag brings pebbles to protocores but also sweeps them past; those particles with th… ▽ More Pebble accretion refers to the assembly of rocky planet cores from particles whose velocity dispersions are damped by drag from circumstellar disc gas. Accretion cross-sections can approach maximal Hill-sphere scales for particles whose Stokes numbers approach unity. While fast, pebble accretion is also lossy. Gas drag brings pebbles to protocores but also sweeps them past; those particles with the largest accretion cross-sections also have the fastest radial drift speeds and are the most easily drained out of discs. We present a global model of planet formation by pebble accretion that keeps track of the disc's mass budget. Cores, each initialized with a lunar mass, grow from discs whose finite stores of mm-cm sized pebbles drift inward across all radii in viscously accreting gas. For every 1 $M_\oplus$ netted by a core, at least 10 $M_\oplus$ and possibly much more are lost to radial drift. Core growth rates are typically exponentially sensitive to particle Stokes number, turbulent Mach number, and solid surface density. This exponential sensitivity, when combined with disc migration, tends to generate binary outcomes from 0.1-30 AU: either sub-Earth cores remain sub-Earth, or explode into Jupiters, with the latter migrating inward to varying degrees. When Jupiter-breeding cores assemble from mm-cm sized pebbles, they do so in discs where such particles drain out in $\sim$10$^5$ yr or less; such fast-draining discs do not fit mm-wave observations. △ Less

Submitted 6 August, 2018; v1 submitted 1 June, 2018; originally announced June 2018.

Comments: Accepted to MNRAS

Optically Thin Core Accretion: How Planets Get Their Gas in Nearly Gas-Free Disks

Authors: Eve J. Lee, Eugene Chiang, Jason W. Ferguson

Abstract: Models of core accretion assume that in the radiative zones of accreting gas envelopes, radiation diffuses. But super-Earths/sub-Neptunes (1-4$R_\oplus$, 2-20$M_\oplus$) point to formation conditions that are optically thin: their modest gas masses are accreted from short-lived and gas-poor nebulae reminiscent of the transparent cavities of transitional disks. Planetary atmospheres born in such en… ▽ More Models of core accretion assume that in the radiative zones of accreting gas envelopes, radiation diffuses. But super-Earths/sub-Neptunes (1-4$R_\oplus$, 2-20$M_\oplus$) point to formation conditions that are optically thin: their modest gas masses are accreted from short-lived and gas-poor nebulae reminiscent of the transparent cavities of transitional disks. Planetary atmospheres born in such environments can be optically thin to both incident starlight and internally generated thermal radiation. We construct time-dependent models of such atmospheres, showing that super-Earths/sub-Neptunes can accrete their $\sim$1%-by-mass gas envelopes, and super-puffs/sub-Saturns their $\sim$20%-by-mass envelopes, over a wide range of nebular depletion histories requiring no fine tuning. Although nascent atmospheres can exhibit stratospheric temperature inversions effected by atomic Fe and various oxides that absorb strongly at visible wavelengths, the rate of gas accretion remains controlled by the radiative-convective boundary (rcb) at much greater pressures. For dusty envelopes, the temperature at the rcb $T_{\rm rcb} \simeq 2500$ K is still set by ${\rm H}_2$ dissociation; for dust-depleted envelopes, $T_{\rm rcb}$ tracks the temperature of the visible or thermal photosphere, whichever is deeper, out to at least $\sim$5 AU. The rate of envelope growth remains largely unchanged between the old radiative diffusion models and the new optically thin models, reinforcing how robustly super-Earths form as part of the endgame chapter in disk evolution. △ Less

Submitted 6 February, 2018; v1 submitted 6 October, 2017; originally announced October 2017.

Comments: accepted to MNRAS, new section 4.2 connects our formation scenario of super-Earths to atmospheric mass loss

Magnetospheric Truncation, Tidal Inspiral, and the Creation of Short and Ultra-Short Period Planets

Authors: Eve J. Lee, Eugene Chiang

Abstract: Sub-Neptunes around FGKM dwarfs are evenly distributed in log orbital period down to $\sim$10 days, but dwindle in number at shorter periods. Both the break at $\sim$10 days and the slope of the occurrence rate down to $\sim$1 day can be attributed to the truncation of protoplanetary disks by their host star magnetospheres at co-rotation. We demonstrate this by deriving planet occurrence rate prof… ▽ More Sub-Neptunes around FGKM dwarfs are evenly distributed in log orbital period down to $\sim$10 days, but dwindle in number at shorter periods. Both the break at $\sim$10 days and the slope of the occurrence rate down to $\sim$1 day can be attributed to the truncation of protoplanetary disks by their host star magnetospheres at co-rotation. We demonstrate this by deriving planet occurrence rate profiles from empirical distributions of pre-main-sequence stellar rotation periods. Observed profiles are better reproduced when planets are distributed randomly in disks---as might be expected if planets formed in situ---rather than piled up near disk edges, as would be the case if they migrated in by disk torques. Planets can be brought from disk edges to ultra-short ($< 1$ day) periods by asynchronous equilibrium tides raised on their stars. Tidal migration can account for how ultra-short period planets (USPs) are more widely spaced than their longer period counterparts. Our picture provides a starting point for understanding why metal-rich stars tend to harbor more USPs, and why the sub-Neptune population drops at $\sim$10 days regardless of whether the host star is of type FGK or early M. We predict planet occurrence rates around A stars to also break at short periods, but at $\sim$1 day instead of $\sim$10 days because A stars rotate faster than lower mass stars (this prediction presumes that the planetesimal building blocks of planets can drift inside the dust sublimation radius). △ Less

Submitted 24 April, 2017; v1 submitted 27 February, 2017; originally announced February 2017.

Comments: ApJ in press

Mass and Eccentricity Constraints in the WASP-47 Planetary System from a Simultaneous Analysis of Radial Velocities & Transit Timing Variations

Authors: Lauren M. Weiss, Katherine Deck, Evan Sinukoff, Erik A. Petigura, Eric Agol, Eve J. Lee, Juliette C. Becker, Andrew W. Howard, Howard Isaacson, Ian J. M. Crossfield, Benjamin J. Fulton, Lea Hirsch

Abstract: Measuring precise planet masses, densities, and orbital dynamics in individual planetary systems is an important pathway toward understanding planet formation. The WASP-47 system has an unusual architecture that motivates a complex formation theory. The system includes a hot Jupiter ("b") neighbored by interior ("e") and exterior ("d") sub-Neptunes, and a long-period eccentric giant planet ("c").… ▽ More Measuring precise planet masses, densities, and orbital dynamics in individual planetary systems is an important pathway toward understanding planet formation. The WASP-47 system has an unusual architecture that motivates a complex formation theory. The system includes a hot Jupiter ("b") neighbored by interior ("e") and exterior ("d") sub-Neptunes, and a long-period eccentric giant planet ("c"). We simultaneously modeled transit times from the Kepler K2 Mission and 118 radial velocities to determine precise masses, densities, and Keplerian orbital elements of the WASP-47 planets. Combining RVs and TTVs provides a better estimate of the mass of planet d (13.6\pm2.0~M_\odot) than obtained with only RVs ($12.75\pm2.70~M_\odot$) or TTVs ($16.1\pm3.8~M_\odot$). Planets e and d have high densities for their size, consistent with a history of photo-evaporation and/or formation in a volatile-poor environment. Through our RV and TTV analysis, we find that the planetary orbits have eccentricities similar to the solar system planets. The WASP-47 system has three similarities to our own solar system: (1) the planetary orbits are nearly circular and coplanar, (2) the planets are not trapped in mean motion resonances, and (3) the planets have diverse compositions. None of the current single-process exoplanet formation theories adequately reproduce these three characteristics of the WASP-47 system (or our solar system). We propose that WASP-47, like the solar system, formed in two stages: first, the giant planets formed in a gas-rich disk and migrated to their present locations, and second, the high-density sub-Neptunes formed in situ in a gas-poor environment. △ Less

Submitted 25 March, 2017; v1 submitted 14 December, 2016; originally announced December 2016.

Comments: Accepted for publication in the Astronomical Journal. 18 pages, 10 figures

Physical properties of molecular clouds for the entire Milky Way disk

Authors: Marc-Antoine Miville-Deschênes, Norman Murray, Eve J. Lee

Abstract: This study presents a catalog of 8107 molecular clouds that covers the entire Galactic plane and includes 98% of the $^{12}$CO emission observed within $b\pm5^\circ$. The catalog was produced using a hierarchical cluster identification method applied to the result of a Gaussian decomposition of the Dame et al. data. The total H$_2$ mass in the catalog is $1.2\times10^9 M_\odot$, in agreement with… ▽ More This study presents a catalog of 8107 molecular clouds that covers the entire Galactic plane and includes 98% of the $^{12}$CO emission observed within $b\pm5^\circ$. The catalog was produced using a hierarchical cluster identification method applied to the result of a Gaussian decomposition of the Dame et al. data. The total H$_2$ mass in the catalog is $1.2\times10^9 M_\odot$, in agreement with previous estimates. We find that 30% of the sight lines intersect only a single cloud, with another 25% intersecting only two clouds. The most probable cloud size is $R\sim30$ pc. We find that $M\propto R^{2.2\pm0.2}$, with no correlation between the cloud surface density, $Σ$, and $R$. In contrast with the general idea, we find a rather large range of values of $Σ$, from 2 to $300 M_\odot$ pc$^{-2}$, and a systematic decrease with increasing Galactic radius, $R_{\rm gal}$. The cloud velocity dispersion and the normalization $σ_0=σ_v/R^{1/2}$ both decrease systematically with $R_{\rm gal}$. When studied over the whole Galactic disk, there is a large dispersion in the line width-size relation, and a significantly better correlation between $σ_v$ and $Σ\,R$. The normalization of this correlation is constant to better than a factor of two for $R_{\rm gal}<20$ kpc. This relation is used to disentangle the ambiguity between near and far kinematic distances. We report a strong variation of the turbulent energy injection rate. In the outer Galaxy it may be maintained by accretion through the disk and/or onto the clouds, but neither source can drive the 100 times higher cloud-averaged injection rate in the inner Galaxy. △ Less

Submitted 5 December, 2016; v1 submitted 19 October, 2016; originally announced October 2016.

Comments: 31 pages, accepted by ApJ

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

Authors: Eve J. Lee, Marc-Antoine Miville-Deschenes, Norman Murray

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

Submitted 9 September, 2016; v1 submitted 18 August, 2016; originally announced August 2016.

Comments: accepted to ApJ; we added a discussion of multi-freefall models of star formation

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

Authors: Thomas M. Esposito, Michael P. Fitzgerald, James R. Graham, Paul Kalas, Eve J. Lee, Eugene Chiang, Gaspard Duchene, Jason Wang, Maxwell A. Millar-Blanchaer, Eric Nielsen, S. Mark Ammons, Sebastian Bruzzone, Robert J. De Rosa, Zachary H. Draper, Bruce Macintosh, Franck Marchis, Stanimir A. Metchev, Marshall Perrin, Laurent Pueyo, Abhijith Rajan, Fredrik T. Rantakyro, David Vega, Schuyler Wolff

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

Submitted 16 September, 2016; v1 submitted 19 May, 2016; originally announced May 2016.

Comments: Accepted to AJ; added Figure 5 and minor text edits

Journal ref: The Astronomical Journal 152 (2016) 85

A Primer on Unifying Debris Disk Morphologies

Authors: Eve J. Lee, Eugene Chiang

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

Submitted 9 June, 2016; v1 submitted 19 May, 2016; originally announced May 2016.

Comments: Accepted to ApJ; minor edits made

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

Authors: Rebekah I. Dawson, Eve J. Lee, Eugene Chiang

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

Submitted 22 March, 2016; v1 submitted 15 December, 2015; originally announced December 2015.

Comments: 19 pages, 14 figures. Accepted by ApJ (in press)

Two Transiting Low Density Sub-Saturns from K2

Authors: Erik A. Petigura, Andrew W. Howard, Eric D. Lopez, Katherine M. Deck, Benjamin J. Fulton, Ian J. M. Crossfield, David R. Ciardi, Eugene Chiang, Eve J. Lee, Howard Isaacson, Charles A. Beichman, Brad M. S. Hansen, Joshua E. Schlieder, Evan Sinukoff

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

Submitted 9 June, 2016; v1 submitted 13 November, 2015; originally announced November 2015.

Comments: 12 pages, 6 figures, 3 tables, published in ApJ. v2: text updated to ApJ version, replaced EPIC-2037 with K2-24

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

Authors: Eve J. Lee, Eugene Chiang

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

Submitted 24 November, 2015; v1 submitted 29 October, 2015; originally announced October 2015.

Comments: Accepted to ApJ

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

Authors: Eve J. Lee, Eugene Chiang

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

Submitted 20 August, 2015; originally announced August 2015.

Comments: 9 pages, 6 figures, accepted to ApJ

A Metallicity Recipe for Rocky Planets

Authors: Rebekah I. Dawson, Eugene Chiang, Eve J. Lee

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

Submitted 18 August, 2015; v1 submitted 23 June, 2015; originally announced June 2015.

Comments: Submitted to MNRAS on April 22, 2015. Accepted July 15, 2015. In press

Helical Blowout Jets in the Sun: Untwisting and Propagation of Waves

Authors: E. J. Lee, V. Archontis, A. W. Hood

Abstract: We report on numerical experiment of the recurrent onset of helical "blowout" jets in an emerging flux region (EFR). We find that these jets are running with velocities of $\sim 100-250\,\mathrm{kms^{-1}}$ and they transfer a vast amount of heavy plasma into the outer solar atmosphere. During their emission, they undergo an untwisting motion as a result of reconnection between the twisted emerging… ▽ More We report on numerical experiment of the recurrent onset of helical "blowout" jets in an emerging flux region (EFR). We find that these jets are running with velocities of $\sim 100-250\,\mathrm{kms^{-1}}$ and they transfer a vast amount of heavy plasma into the outer solar atmosphere. During their emission, they undergo an untwisting motion as a result of reconnection between the twisted emerging and the non-twisted pre-existing magnetic field in the solar atmosphere. For the first time in the context of blowout jets, we provide a direct evidence that their untwisting motion is associated with the propagation of torsional Alfvén waves in the corona. △ Less

Submitted 15 December, 2014; originally announced December 2014.

Comments: 6 pages, 5 figures, Accepted for publication in the Astrophysical Journal Letters

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

Authors: Eve J. Lee, Eugene Chiang, Chris W. Ormel

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

Submitted 8 December, 2014; v1 submitted 11 September, 2014; originally announced September 2014.

Comments: 18 pages, 10 figures, Final proofed version, 2014 ApJ, 797, 95. Includes conceptual corrections to Section 4.2 on the amount of nebular depletion required for scenario (b)

Time Varying Dynamical Star Formation Rate

Authors: Eve J. Lee, Philip Chang, Norman Murray

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

Submitted 3 February, 2015; v1 submitted 16 June, 2014; originally announced June 2014.

Comments: 14 pages, 14 figures, to appear in ApJ. Added purely compressive turbulence driving for FLASH runs

Milky Way Star Forming Complexes and the Turbulent Motion of the Galaxy's Molecular Gas

Authors: E. J. Lee, N. Murray, M. Rahman

Abstract: We analyze Spitzer GLIMPSE, MSX, and WMAP images of the Milky Way to identify 8 micron and free-free sources in the Galaxy. Seventy-two of the eighty-eight WMAP sources have coverage in the GLIMPSE and MSX surveys suitable for identifying massive star forming complexes (SFC). We measure the ionizing luminosity functions of the SFCs and study their role in the turbulent motion of the Galaxy's molec… ▽ More We analyze Spitzer GLIMPSE, MSX, and WMAP images of the Milky Way to identify 8 micron and free-free sources in the Galaxy. Seventy-two of the eighty-eight WMAP sources have coverage in the GLIMPSE and MSX surveys suitable for identifying massive star forming complexes (SFC). We measure the ionizing luminosity functions of the SFCs and study their role in the turbulent motion of the Galaxy's molecular gas. We find a total Galactic free-free flux f_ν = 46177.6 Jy; the 72 WMAP sources with full 8 micron coverage account for 34263.5 Jy (~75%), with both measurements made at ν=94GHz (W band). We find a total of 280 SFCs, of which 168 have unique kinematic distances and free-free luminosities. We use a simple model for the radial distribution of star formation to estimate the free-free and ionizing luminosity for the sources lacking distance determinations. The total dust-corrected ionizing luminosity is Q = 2.9 \pm 0.5 x 10^53 photons s^-1, which implies a galactic star formation rate of 1.2 \pm 0.2 M_{\sun} yr^-1. We present the (ionizing) luminosity function of the SFCs, and show that 24 sources emit half the ionizing luminosity of the Galaxy. The SFCs appear as bubbles in GLIMPSE or MSX images; the radial velocities associated with the bubble walls allow us to infer the expansion velocity of the bubbles. We calculate the kinetic luminosity of the bubble expansion and compare it to the turbulent luminosity of the inner molecular disk. SFCs emitting 80% of the total galactic free-free luminosity produce a kinetic luminosity equal to 65% of the turbulent luminosity in the inner molecular disk. This suggests that the expansion of the bubbles is a major driver of the turbulent motion of the inner Milky Way molecular gas. △ Less

Submitted 18 April, 2012; originally announced April 2012.

Comments: 40 pages, 12 figures, Accepted for publication in ApJ

Radiation-Hydrodynamic Simulations of Protostellar Outflows: Synthetic Observations and Data Comparisons

Authors: S. S. R. Offner, E. J. Lee, A. A. Goodman, H. Arce

Abstract: We present results from three-dimensional, self-gravitating, radiation-hydrodynamic simulations of low-mass protostellar outflows. We construct synthetic observations in 12CO in order to compare with observed outflows and evaluate the effects of beam resolution and outflow orientation on inferred outflow properties. To facilitate the comparison, we develop a quantitative prescription for measuring… ▽ More We present results from three-dimensional, self-gravitating, radiation-hydrodynamic simulations of low-mass protostellar outflows. We construct synthetic observations in 12CO in order to compare with observed outflows and evaluate the effects of beam resolution and outflow orientation on inferred outflow properties. To facilitate the comparison, we develop a quantitative prescription for measuring outflow opening angles. Using this prescription, we demonstrate that, in both simulations and synthetic observations, outflow opening angles broaden with time similarly to observed outflows. However, the interaction between the outflowing gas and the turbulent core envelope produces significant asymmetry between the red and blue shifted outflow lobes. We find that applying a velocity cutoff may result in outflow masses that are underestimated by a factor 5 or more, and masses derived from optically thick CO emission further underpredict the mass of the high-velocity gas by a factor of 5-10. Derived excitation temperatures indicate that outflowing gas is hotter than the ambient gas with temperature rising over time, which is in agreement with the simulation gas temperatures. However, excitation temperatures are otherwise not well correlated with the actual gas temperature. △ Less

Submitted 25 October, 2011; originally announced October 2011.

Comments: 13 pages, accepted to ApJ