Recently, we introduced PLanetary Atmospheric Tool for Observer Noobs (PLATON), a Python package that calculates model transmission spectra for exoplanets and retrieves atmospheric characteristics ...based on observed spectra. We now expand its capabilities to include the ability to compute secondary eclipse depths. We have also added the option to calculate models using the correlated-k method for radiative transfer, which improves accuracy without sacrificing speed. Additionally, we update the opacities in PLATON-many of which were generated using old or proprietary line lists-using the most recent and complete public line lists. These opacities are made available at R = 1000 and R = 10,000 over the 0.3-30 m range, and at R = 375,000 in select near-IR bands, making it possible to utilize PLATON for ground-based high-resolution cross-correlation studies. To demonstrate PLATON's new capabilities, we perform a retrieval on published Hubble Space Telescope (HST) and Spitzer transmission and emission spectra of the archetypal hot Jupiter HD 189733b. This is the first joint transit and secondary eclipse retrieval for this planet in the literature, as well as the most comprehensive set of both transit and eclipse data assembled for a retrieval to date. We find that these high signal-to-noise data are well matched by atmosphere models with a C/O ratio of and a metallicity of times solar where the terminator is dominated by extended nanometer-sized haze particles at optical wavelengths. These are among the smallest uncertainties reported to date for an exoplanet, demonstrating both the power and the limitations of HST and Spitzer exoplanet observations.
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 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
M
̇
⋆
∼
10
−
8
M
⊙
yr
−
1
(
M
⋆
/
M
⊙
)
2
. We predict a decline in mini-Neptune population (but not necessarily terrestrial planets) around stars lighter than ∼0.3–0.5
M
⊙
. Cold giants (≳5 au), if they exist, should remain correlated with inner planets even around lower-mass stars.
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 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 the 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
α
t
and the velocity at which grains fragment
v
frag
. For moderate-to-low turbulence strengths of
α
t
≲ 10
−3
and substantial differences in
v
frag
for icy and ice-free grains, we find a spatially nonuniform dust-to-gas ratio and grain size distribution that deviate significantly from the ISM values, in agreement with previous studies. The effect of a nonuniform 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 the 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.
Abstract
The metastable helium line at 1083 nm can be used to probe the extended upper atmospheres of close-in exoplanets and thus provide insight into their atmospheric mass loss, which is likely to ...be significant in sculpting their population. We used an ultra-narrow band filter centered on this line to observe two transits of the low-density gas giant HAT-P-18b, using the 200″ Hale Telescope at Palomar Observatory, and report the detection of its extended upper atmosphere. We constrain the excess absorption to be 0.46% ± 0.12% in our 0.635 nm bandpass, exceeding the transit depth from the Transiting Exoplanet Survey Satellite (TESS) by 3.9
σ
. If we fit this signal with a 1D Parker wind model, we find that it corresponds to an atmospheric mass loss rate between
8.3
−
1.9
+
2.8
×
10
−
5
M
J
Gyr
−1
and
2.63
−
0.64
+
0.46
×
10
−
3
M
J
Gyr
−1
for thermosphere temperatures ranging from 4000 K to 13,000 K, respectively. With a
J
magnitude of 10.8, this is the faintest system for which such a measurement has been made to date, demonstrating the effectiveness of this approach for surveying mass loss on a diverse sample of close-in gas giant planets.
Abstract
Relating planet formation to atmospheric composition has been a long-standing goal of the planetary science community. So far, most modeling studies have focused on predicting the enrichment ...of heavy elements and the C/O ratio in giant planet atmospheres. Although this framework provides useful constraints on the potential formation locations of gas giant exoplanets, carbon and oxygen measurements alone are not enough to determine where a given gas giant planet originated. Here, we show that characterizing the abundances of refractory elements (e.g., silicon and iron) can break these degeneracies. Refractory elements are present in the solid phase throughout most of the disk, and their atmospheric abundances therefore reflect the solid-to-gas accretion ratio during formation. We introduce a new framework that parameterizes the atmospheric abundances of gas giant exoplanets in the form of three ratios: Si/H, O/Si, and C/Si. Si/H traces the solid-to-gas accretion ratio of a planet and is loosely equivalent to earlier notions of “metallicity.” For O/Si and C/Si, we present a global picture of their variation with distance and time based on what we know from the solar system meteorites and an updated understanding of the variations of thermal processing within protoplanetary disks. We show that ultrahot Jupiters are ideal targets for atmospheric characterization studies using this framework as we can measure the abundances of refractories, oxygen, and carbon in the gas phase. Finally, we propose that hot Jupiters with silicate clouds and low water abundances might have accreted their envelopes between the soot line and the water snow line.
Abstract
Kepler-289 is a three-planet system containing two sub-Neptunes and one cool giant planet orbiting a young, Sun-like star. All three planets exhibit transit timing variations (TTVs), with ...both adjacent planet pairs having orbital periods close to the 2:1 orbital resonance. We observe two transits of Kepler-289c with the Wide-field InfraRed Camera on the 200″ Hale Telescope at Palomar Observatory, using diffuser-assisted photometry to achieve space-like photometric precision from the ground. These new transit observations extend the original four-year Kepler TTV baseline by an additional 7.5 yr. We rereduce the archival Kepler data with an improved stellar activity correction and carry out a joint fit with the Palomar data to constrain the transit shapes and derive updated transit times. We then model the TTVs to determine the masses of the three planets and constrain their densities and bulk compositions. Our new analysis improves on previous mass and density constraints by a factor of two or more for all three planets, with the innermost planet showing the largest improvement. Our updated atmospheric mass fractions for the inner two planets indicate that they have hydrogen-rich envelopes, consistent with their location on the upper side of the radius valley. We also constrain the heavy element composition of the outer Saturn-mass planet, Kepler-289c, for the first time, finding that it contains 30.5 ± 6.9
M
⊕
of metals. We use dust evolution models to show that Kepler-289c must have formed beyond 1 au, and likely beyond 3 au, and then migrated inward.
Abstract
Early in their lives, planets endure extreme amounts of ionizing radiation from their host stars. For planets with primordial hydrogen and helium-rich envelopes, this can lead to substantial ...mass loss. Direct observations of atmospheric escape in young planetary systems can help elucidate this critical stage of planetary evolution. In this work, we search for metastable helium absorption—a tracer of tenuous gas in escaping atmospheres—during transits of three planets orbiting the young solar analog V1298 Tau. We characterize the stellar helium line using HET/HPF, and find that it evolves substantially on timescales of days to months. The line is stable on hour-long timescales except for one set of spectra taken during the decay phase of a stellar flare, where absoprtion increased with time. Utilizing a beam-shaping diffuser and a narrowband filter centered on the helium feature, we observe four transits with Palomar/WIRC: two partial transits of planet d (
P
= 12.4 days), one partial transit of planet b (
P
= 24.1 days), and one full transit of planet c (
P
= 8.2 days). We do not detect the transit of planet c, and we find no evidence of excess absorption for planet b, with Δ
R
b
/
R
⋆
< 0.019 in our bandpass. We find a tentative absorption signal for planet d with Δ
R
d
/
R
⋆
= 0.0205 ± 0.054, but the best-fit model requires a substantial (−100 ± 14 minutes) transit-timing offset on a two-month timescale. Nevertheless, our data suggest that V1298 Tau d may have a high present-day mass-loss rate, making it a priority target for follow-up observations.
We present the transmission spectrum of HAT-P-12b through a joint analysis of data obtained from the Hubble Space Telescope Space Telescope Imaging Spectrograph and Wide Field Camera 3 and Spitzer, ...covering the wavelength range 0.3-5.0 m. We detect a muted water vapor absorption feature at 1.4 m attenuated by clouds, as well as a Rayleigh scattering slope in the optical indicative of small particles. We interpret the transmission spectrum using both the state-of-the-art atmospheric retrieval code SCARLET and the aerosol microphysics model CARMA. These models indicate that the atmosphere of HAT-P-12b is consistent with a broad range of metallicities between several tens to a few hundred times solar, a roughly solar C/O ratio, and moderately efficient vertical mixing. Cloud models that include condensate clouds do not readily generate the submicron particles necessary to reproduce the observed Rayleigh scattering slope, while models that incorporate photochemical hazes composed of soot or tholins are able to match the full transmission spectrum. From a complementary analysis of secondary eclipses by Spitzer, we obtain measured depths of 0.042% 0.013% and 0.045% 0.018% at 3.6 and 4.5 m, respectively, which are consistent with a blackbody temperature of 890+60−70 K and indicate efficient day-night heat recirculation. HAT-P-12b joins the growing number of well-characterized warm planets that underscore the importance of clouds and hazes in our understanding of exoplanet atmospheres.
Abstract Previous attempts have been made to characterize the atmospheres of directly imaged planets at low resolution ( R ∼ 10–100 s), but the presence of clouds has often led to degeneracies in the ...retrieved atmospheric abundances with cloud opacity and temperature structure that bias retrieved compositions. In this study, we perform retrievals on the ultrayoung (≲5 Myr) directly imaged planet ROXs 42B b with both a downsampled low-resolution JHK -band spectrum from Gemini/NIFS and Keck/OSIRIS, and a high-resolution K -band spectrum from pre-upgrade Keck/NIRSPAO. Using the atmospheric retrieval framework of petitRADTRANS , we analyze both data sets individually and combined. We additionally fit for the stellar abundances and other physical properties of the host stars, a young M spectral type binary, using the SPHINX model grid. We find that the measured C/O, 0.50 ± 0.05, and metallicity, Fe/H = −0.67 ± 0.35, for ROXs 42B b from our high-resolution spectrum agree with those of its host stars within 1 σ . The retrieved parameters from the high-resolution spectrum are also independent of our choice of cloud model. In contrast, the retrieved parameters from the low-resolution spectrum show strong degeneracies between the clouds and the retrieved metallicity and temperature structure. When we retrieve both data sets together, we find that these degeneracies are reduced but not eliminated, and the final results remain highly sensitive to cloud modeling choices. We conclude that high-resolution spectroscopy offers the most promising path for reliably determining atmospheric compositions of directly imaged companions independent of their cloud properties.
Low-mass low-density planets discovered by Kepler in the super-Earth mass regime typically have large radii for their inferred masses, implying the presence of H2-He atmospheres. These planets are ...vulnerable to atmospheric mass loss due to heating by the parent star's XUV flux. Models coupling atmospheric mass loss with thermal evolution predicted a bimodal distribution of planetary radii, which has gained observational support. However, a key component that has been ignored in previous studies is the dissolution of these gases into the molten core of rock and iron that constitute most of their mass. Such planets have high temperatures (>2000 K) and pressures (∼kbars) at the core-envelope boundary, ensuring a molten surface and a subsurface reservoir of hydrogen that can be 5-10 times larger than the atmosphere. This study bridges this gap by coupling the thermal evolution of the planet and the mass loss of the atmosphere with the thermodynamic equilibrium between the dissolved H2 and the atmospheric H2 (Henry's law). Dissolution in the interior allows a planet to build a larger hydrogen repository during the planet formation stage. We show that the dissolved hydrogen outgasses to buffer atmospheric mass loss. The slow cooling of the planet also leads to outgassing because solubility decreases with decreasing temperature. Dissolution of hydrogen in the interior therefore increases the atmosphere retention ability of super-Earths. The study highlights the importance of including the temperature- and pressure-dependent solubility of gases in magma oceans and coupling outgassing to planetary evolution models.