Researcher profile

Jonathan J. Fortney

Jonathan J. Fortney contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
26works
0followers
3topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

26 published item(s)

preprint2022arXiv

A 1.46-2.48 $μ$m Spectroscopic Atlas of a T6 Dwarf (1060 K) Atmosphere with IGRINS: First Detections of H$_2$S and H$_2$, and Verification of H$_2$O, CH$_4$, and NH$_3$ Line Lists

We present Gemini South/IGRINS observations of the 1060 K T6 dwarf 2MASS J08173001$-$6155158 with unprecedented resolution ($R\equivλ/Δλ=45\,000$) and signal-to-noise ratio (SNR > 200) for a late-type T dwarf. We use this benchmark observation to test the reliability of molecular line lists used up-to-date atmospheric models. We determine which spectroscopic regions should be used to estimate the parameters of cold brown dwarfs and, by extension, exoplanets. We present a detailed spectroscopic atlas with molecular identifications across the $H$ and $K$ bands of the near-infrared. We find that water (H$_2$O) line lists are overall reliable. We find the most discrepancies amongst older methane (CH$_4$) line lists, and that the most up-to-date CH$_4$ line lists correct many of these issues. We identify individual ammonia (NH$_3$) lines, a hydrogen sulfide (H$_2$S) feature at 1.5900 $μ$m, and a molecular hydrogen (H$_2$) feature at 2.1218 $μ$m. These are the first unambiguous detections of H$_2$S and H$_2$ absorption features in an extra-solar atmosphere. With the H$_2$ detection, we place an upper limit on the atmospheric dust concentration of this T6 dwarf: at least 500 times less than the interstellar value, implying that the atmosphere is effectively dust-free. We additionally identify several features that do not appear in the model spectra. Our assessment of the line lists is valuable for atmospheric model applications to high-dispersion, low-SNR, high-background spectra, such as an exoplanet around a star. We demonstrate a significant enhancement in the detection of the CH$_4$ absorption signal in this T6 dwarf with the most up-to-date line lists.

preprint2022arXiv

A Framework for Characterizing Transmission Spectra of Exoplanets with Circumplanetary Rings

Recent observations revealed that several extremely low-density exoplanets show featureless transmission spectra. While atmospheric aerosols are a promising explanation for both the low density and featureless spectra, there is another attractive possibility: the presence of circumplanetary rings. Previous studies suggested that rings cause anomalously large transit radii. However, it remains poorly understood how rings affect the transmission spectrum. Here, we provide a framework to characterize the transmission spectra of ringed exoplanets. We develop an analytical prescription to include rings in the transmission spectra for arbitrarily viewing geometries. We also establish a simple post-processing model that can include the ring's effects on precomputed ring-free spectra. The ring flattens the transmission spectrum for a wide range of viewing geometries, consistent with the featureless spectra of extremely low-density exoplanets. Near-future observations by JWST at longer wavelengths would be able to distinguish the aerosol and ring scenarios. We also find that rocky rings might cause a silicate feature at $\sim$10 $μ$m if the ring's optical depth is around unity. Thus, the ring's spectral features, if detected, would provide tight constrains on the physical properties of exoplanetary rings. We also discuss the ring's stability and suggest that thick rings are sustainable only at the equilibrium temperature of $\lesssim$300 K for the ring's age comparable to Kepler planets. This might indicate the intrinsic deficit of thick rings in the Kepler samples, unless rings are much younger than the planets as suggested for Saturn.

preprint2022arXiv

A new method to correct for host star variability in multi-epoch observations of exoplanet transmission spectra

Transmission spectra of exoplanets orbiting active stars suffer from wavelength-dependent effects due to stellar photospheric heterogeneity. WASP-19b, an ultra-hot Jupiter (T$_{eq}$ $\sim$ 2100 K), is one such strongly irradiated gas-giant orbiting an active solar-type star. We present optical (520-900 nm) transmission spectra of WASP-19b obtained across eight epochs using the Gemini Multi-Object Spectrograph (GMOS) on the Gemini-South telescope. We apply our recently developed Gaussian Processes regression based method to model the transit light curve systematics and extract the transmission spectrum at each epoch. We find that WASP-19b's transmission spectrum is affected by stellar variability at individual epochs. We report an observed anticorrelation between the relative slopes and offsets of the spectra across all epochs. This anticorrelation is consistent with the predictions from the forward transmission models, which account for the effect of unocculted stellar spots and faculae measured previously for WASP-19. We introduce a new method to correct for this stellar variability effect at each epoch by using the observed correlation between the transmission spectral slopes and offsets. We compare our stellar variability corrected GMOS transmission spectrum with previous contradicting MOS measurements for WASP-19b and attempt to reconcile them. We also measure the amplitude and timescale of broadband stellar variability of WASP-19 from TESS photometry, which we find to be consistent with the effect observed in GMOS spectroscopy and ground-based broadband photometric long-term monitoring. Our results ultimately caution against combining multi-epoch optical transmission spectra of exoplanets orbiting active stars before correcting each epoch for stellar variability.

preprint2022arXiv

Confirmation of Water Absorption in the Thermal Emission Spectrum of the Hot Jupiter WASP-77Ab with HST/WFC3

Secondary eclipse observations of hot Jupiters can reveal both their compositions and thermal structures. Previous observations have shown a diversity of hot Jupiter eclipse spectra, including absorption features, emission features, and featureless blackbody-like spectra. We present a secondary eclipse spectrum of the hot Jupiter WASP-77Ab observed between $1-5$ $μ$m with the Hubble Space Telescope (HST) and the Spitzer Space Telescope. The HST observations show signs of water absorption indicative of a non-inverted thermal structure. We fit the data with both a one-dimensional free retrieval and a grid of one-dimensional self-consistent forward models to confirm this non-inverted structure. The free retrieval places a $3σ$ lower limit on the atmospheric water abundance of $\log(n_\mathrm{H_2O})>-4.78$ and can not constrain the CO abundance. The grid fit produces a slightly super-stellar metallicity and constrains the carbon-to-oxygen ratio to less than or equal to the solar value. We also compare our data to recent high-resolution observations of WASP-77Ab taken with the Gemini-South/IGRINS spectrograph and find that our observations are consistent with the best-fit model to the high-resolution data. However, the metallicity derived from the IGRINS data is significantly lower than that derived from our self-consistent model fit. We find that this difference may be due to disequilibrium chemistry, and the varying results between the models applied here demonstrate the difficulty of constraining disequilibrium chemistry with low-resolution, low wavelength coverage data alone. Future work to combine observations from IGRINS, HST, and JWST will improve our estimate of the atmospheric composition of WASP-77Ab.

preprint2022arXiv

Microphysics of Water Clouds in the Atmospheres of Y Dwarfs and Temperate Giant Planets

Water clouds are expected to form on Y dwarfs and giant planets with equilibrium temperatures near or below that of Earth, drastically altering their atmospheric compositions and their albedos and thermal emission spectra. Here we use the 1D Community Aerosol and Radiation Model for Atmospheres (CARMA) to investigate the microphysics of water clouds on cool substellar worlds to constrain their typical particle sizes and vertical extent, taking into consideration nucleation and condensation, which have not been considered in detail for water clouds in H/He atmospheres. We compute a small grid of Y dwarf and temperate giant exoplanet atmosphere models with water clouds forming through homogeneous nucleation and heterogeneous nucleation on cloud condensation nuclei composed of meteoritic dust, organic photochemical hazes, and upwelled potassium chloride cloud particles. We present comparisons with the Ackerman & Marley parameterization of cloud physics to extract the optimal sedimentation efficiency parameter (f$_{sed}$) using Virga. We find that no Virga model replicates the CARMA water clouds exactly and that a transition in f$_{sed}$ occurs from the base of the cloud to the cloud top. Furthermore, we generate simulated thermal emission and geometric albedo spectra and find large, wavelength-dependent differences between the CARMA and Virga models, with different gas absorption bands reacting differently to the different cloud distributions and particularly large differences in the M band. Therefore, constraining the vertically-dependent properties of water clouds will be essential to estimating the gas abundances in these atmospheres.

preprint2022arXiv

PICASO 3.0: A One-Dimensional Climate Model for Giant Planets and Brown Dwarfs

Upcoming James Webb Space Telescope (JWST) observations will allow us to study exoplanet and brown dwarf atmospheres in great detail. The physical interpretation of these upcoming high signal-to-noise observations requires precise atmospheric models of exoplanets and brown dwarfs. While several one-dimensional and three-dimensional atmospheric models have been developed in the past three decades, these models have often relied on simplified assumptions like chemical equilibrium and are also often not open-source, which limits their usage and development by the wider community. We present a python-based one-dimensional atmospheric radiative-convective equilibrium model. This model has heritage from the Fortran-based code (Marley et al.,1996} which has been widely used to model the atmospheres of Solar System objects, brown dwarfs, and exoplanets. In short, the basic capability of the original model is to compute the atmospheric state of the object under radiative-convective equilibrium given its effective or internal temperature, gravity, and host--star properties (if relevant). In the new model, which has been included within the well-utilized code-base PICASO, we have added these original features as well as the new capability of self-consistently treating disequilibrium chemistry. This code is widely applicable to Hydrogen-dominated atmospheres (e.g., brown dwarfs and giant planets).

preprint2022arXiv

Predictions for Observable Atmospheres of Trappist-1 Planets from a Fully Coupled Atmosphere-Interior Evolution Model

The Trappist-1 planets provide a unique opportunity to test the current understanding of rocky planet evolution. The James Webb Space Telescope is expected to characterize the atmospheres of these planets, potentially detecting CO$_2$, CO, H$_2$O, CH$_4$, or abiotic O$_2$ from water photodissociation and subsequent hydrogen escape. Here, we apply a coupled atmosphere-interior evolution model to the Trappist-1 planets to anticipate their modern atmospheres. This model, which has previously been validated for Earth and Venus, connects magma ocean crystallization to temperate geochemical cycling. Mantle convection, magmatic outgassing, atmospheric escape, crustal oxidation, a radiative-convective climate model, and deep volatile cycling are explicitly coupled to anticipate bulk atmospheres and planetary redox evolution over 8 Gyr. By adopting a Monte Carlo approach that samples a broad range of initial conditions and unknown parameters, we make some tentative predictions about current Trappist-1 atmospheres. We find that anoxic atmospheres are probable, but not guaranteed, for the outer planets; oxygen produced via hydrogen loss during the pre-main sequence is typically consumed by crustal sinks. In contrast, oxygen accumulation on the inner planets occurs in around half of all models runs. Complete atmospheric erosion is possible but not assured for the inner planets (occurs in 20%-50% of model runs), whereas the outer planets retain significant surface volatiles in virtually all model simulations. For all planets that retain substantial atmospheres, CO$_2$-dominated or CO$_2$-O$_2$ atmospheres are expected; water vapor is unlikely to be a detectable atmospheric constituent in most cases. There are necessarily many caveats to these predictions, but the ways in which they misalign with upcoming observations will highlight gaps in terrestrial planet knowledge.

preprint2022arXiv

Solar-to-supersolar sodium and oxygen absolute abundances for a "hot Saturn" orbiting a metal-rich star

We present new analysis of infrared transmission spectroscopy of the cloud-free hot-Saturn WASP-96b performed with the Hubble and Spitzer Space Telescopes (HST and Spitzer). The WASP-96b spectrum exhibits the absorption feature from water in excellent agreement with synthetic spectra computed assuming a cloud-free atmosphere. The HST-Spitzer spectrum is coupled with Very Large Telescope (VLT) optical transmission spectroscopy which reveals the full pressure-broadened profile of the sodium absorption feature and enables the derivation of absolute abundances. We confirm and correct for a spectral offset of $ΔR_{\rm p}/R_{\ast}=(-4.29^{+0.31}_{-0.37})\,\times10^{-3}$ of the VLT data relative to the HST-Spitzer spectrum. This offset can be explained by the assumed radius for the common-mode correction of the VLT spectra, which is a well-known feature of ground-based transmission spectroscopy. We find evidence for a lack of chromospheric and photometric activity of the host star which, therefore, make a negligible contribution to the offset. We measure abundances for Na and O that are consistent with solar to supersolar, with abundances relative to solar values of $21^{+27}_{-14}$ and $7^{+11}_{-4}$, respectively. We complement the transmission spectrum with new thermal emission constraints from Spitzer observations at 3.6 and $4.5μ$m, which are best explained by the spectrum of an atmosphere with a temperature decreasing with altitude. A fit to the spectrum assuming an isothermal blackbody atmosphere constrains the dayside temperature to be $T_{\rm{p}}$=$1545$$\pm$$90$K.

preprint2022arXiv

The Case and Context for Atmospheric Methane as an Exoplanet Biosignature

Methane has been proposed as an exoplanet biosignature. Imminent observations with the James Webb Space Telescope may enable methane detections on potentially habitable exoplanets, so it is essential to assess in what planetary contexts methane is a compelling biosignature. Methane's short photochemical lifetime in terrestrial planet atmospheres implies that abundant methane requires large replenishment fluxes. While methane can be produced by a variety of abiotic mechanisms such as outgassing, serpentinizing reactions, and impacts, we argue that, in contrast to an Earth-like biosphere, known abiotic processes cannot easily generate atmospheres rich in CH$_4$ and CO$_2$ with limited CO due to the strong redox disequilibrium between CH$_4$ and CO$_2$. Methane is thus more likely to be biogenic for planets with 1) a terrestrial bulk density, high mean-molecular-weight and anoxic atmosphere, and an old host star; 2) an abundance of CH$_4$ that implies surface fluxes exceeding what could be supplied by abiotic processes; and 3) atmospheric CO$_2$ with comparatively little CO.

preprint2022arXiv

The Promise and Limitations of Precision Gravity: Application to the Interior Structure of Uranus and Neptune

We study the constraining power of a high-precision measurement of the gravity field for Uranus and Neptune, as could be delivered by a low periapse orbiter. Our study is practical, assessing the possible deliverables and limitations of such a mission with respect to the structure of the planets. Our study is also academic, assessing in a general way the relative importance of the low order gravity, high order gravity, rotation rate, and moment of inertia (MOI) in constraining planetary structure. We attempt to explore all possible interior density structures of a planet that are consistent with hypothetical gravity data, via MCMC sampling of parameterized density profiles. When the gravity field is poorly known, as it is today, uncertainties in the rotation rate on the order of 10 minutes are unimportant, as they are interchangeable with uncertainties in the gravity coefficients. By the same token, when the gravity field is precisely determined the rotation rate must be known to comparable precision. When gravity and rotation are well known the MOI becomes well-constrained, limiting the usefulness of independent MOI determinations unless they are extraordinarily precise. For Uranus and Neptune, density profiles can be well-constrained. However, the non-uniqueness of the relative roles of H/He, watery volatiles, and rock in the deep interior will still persist with high-precision gravity data. Nevertheless, the locations and magnitudes (in pressure-space) of any large-scale composition gradient regions can likely be identified, offering a crucially better picture of the interiors of Uranus or Neptune.

preprint2022arXiv

Understanding planetary context to enable life detection on exoplanets and test the Copernican principle

The search for life on exoplanets is motivated by the universal ways in which life could modify its planetary environment. Atmospheric gases such as oxygen and methane are promising candidates for such environmental modification due to the evolutionary benefits their production would confer. However, confirming that these gases are produced by life, rather than by geochemical or astrophysical processes, will require a thorough understanding of planetary context, including the expected counterfactual atmospheric evolution for lifeless planets. Here, we evaluate current understanding of planetary context for several candidate biosignatures and their upcoming observability. We review the contextual framework for oxygen and describe how conjectured abiotic oxygen scenarios may be testable. In contrast to oxygen, current understanding of how planetary context controls non-biological methane (CH$_4$) production is limited, even though CH$_4$ biosignatures in anoxic atmospheres may be readily detectable with the James Webb Space Telescope. We assess environmental context for CH$_4$ biosignatures and conclude that abundant atmospheric CH$_4$ coexisting with CO$_2$, and CO:CH$_4$ << 1 is suggestive of biological production, although precise thresholds are dependent on stellar context and sparsely characterized abiotic CH$_4$ scenarios. A planetary context framework is also considered for alternative or agnostic biosignatures. Whatever the distribution of life in the Universe, observations of terrestrial exoplanets in coming decades will provide a quantitative understanding of the atmospheric evolution of lifeless worlds. This knowledge will inform future instrument requirements to either corroborate the presence of life elsewhere or confirm its apparent absence.

preprint2021arXiv

The Sonora Substellar Atmosphere Models. II. Cholla: A Grid of Cloud-free, Solar Metallicity Models in Chemical Disequilibrium for the JWST Era

Exoplanet and brown dwarf atmospheres commonly show signs of disequilibrium chemistry. In the James Webb Space Telescope era high resolution spectra of directly imaged exoplanets will allow the characterization of their atmospheres in more detail, and allow systematic tests for the presence of chemical species that deviate from thermochemical equilibrium in these atmospheres. Constraining the presence of disequilibrium chemistry in these atmospheres as a function of parameters such as their effective temperature and surface gravity will allow us to place better constrains in the physics governing these atmospheres. This paper is part of a series of works presenting the Sonora grid of atmosphere models (Marley et al 2021, Morley et al in prep.). In this paper we present a grid of cloud-free, solar metallicity atmospheres for brown dwarfs and wide separation giant planets with key molecular species such as CH4, H2O, CO and NH3 in disequilibrium. Our grid covers atmospheres with Teff~[500 K,1300 K], logg~[3.0,5.5] (cgs) and an eddy diffusion parameter of logKzz=2, 4 and 7 (cgs). We study the effect of different parameters within the grid on the temperature and composition profiles of our atmospheres. We discuss their effect on the near-infrared colors of our model atmospheres and the detectability of CH4, H2O, CO and NH3 using the JWST. We compare our models against existing MKO and Spitzer observations of brown dwarfs and verify the importance of disequilibrium chemistry for T dwarf atmospheres. Finally, we discuss how our models can help constrain the vertical structure and chemical composition of these atmospheres.

preprint2021arXiv

Thermal Phase Curves of XO-3b: an Eccentric Hot Jupiter at the Deuterium Burning Limit

We report \textit{Spitzer} full-orbit phase observations of the eccentric hot Jupiter XO-3b at 3.6 and 4.5 $μ$m. Our new eclipse depth measurements of $1770 \pm 180$ ppm at 3.6 $μ$m and $1610 \pm 70$ ppm at 4.5 $μ$m show no evidence of the previously reported dayside temperature inversion. We also empirically derive the mass and radius of XO-3b and its host star using Gaia DR3&#39;s parallax measurement and find a planetary mass $M_p=11.79 \pm 0.98 ~M_{\rm{Jup}}$ and radius $R_p=1.295 \pm 0.066 ~R_{\rm{Jup}}$. We compare our \textit{Spitzer} observations with multiple atmospheric models to constrain the radiative and advective properties of XO-3b. While the decorrelated 4.5 $μ$m observations are pristine, the 3.6 $μ$m phase curve remains polluted with detector systematics due to larger amplitude intrapixel sensitivity variations in this channel. We focus our analysis on the more reliable 4.5 $μ$m phase curve and fit an energy balance model with solid body rotation to estimate the zonal wind speed and the pressure of the bottom of the mixed layer. Our energy balance model fit suggests an eastward equatorial wind speed of $3.13 ^{+0.26} _{-0.83}$ km/s, an atmospheric mixed layer down to $2.40 ^{+0.92} _{-0.16}$ bar, and Bond albedo of $0.106 ^{+0.008} _{-0.106}$. We assume that the wind speed and mixed layer depth are constant throughout the orbit. We compare our observations with a 1D planet-averaged model predictions at apoapse and periapse and 3D general circulation model (GCM) predictions for XO-3b. We also investigate the inflated radius of XO-3b and find that it would require an unusually large amount of internal heating to explain the observed planetary radius.

preprint2020arXiv

2D Retrieval Frameworks for Hot Jupiter Phase Curves

Spectroscopic phase curves provide unique access to the three-dimensional properties of transiting exoplanet atmospheres. However, a modeling framework must be developed to deliver accurate inferences of atmospheric properties for these complex data sets. Here, we develop an approach to retrieve temperature structures and molecular abundances from phase curve spectra at any orbital phase. In the context of a representative hot Jupiter with a large day-night temperature contrast, we examine the biases in typical one-dimensional (1D) retrievals as a function of orbital phase/geometry, compared to two-dimensional (2D) models that appropriately capture the disk-integrated phase geometry. We guide our intuition by applying our new framework on a simulated HST+Spitzer phase curve data set in which the &#34;truth&#34; is known, followed by an application to the spectroscopic phase curve of the canonical hot Jupiter, WASP-43b. We also demonstrate the retrieval framework on simulated JWST phase curve observations. We apply our new geometric framework to a joint-fit of all spectroscopic phases, assuming longitudinal molecular abundance homogeneity, resulting in an a factor of 2 improvement in abundances precision when compared to individual phase constraints. With a 1D retrieval model on simulated HST+Spitzer data, we find strongly biased molecular abundances for CH$_4$ and CO$_2$ at most orbital phases. With 2D, the day and night profiles retrieved from WASP-43b remain consistent throughout the orbit. JWST retrievals show that a 2D model is strongly favored at all orbital phases. Based on our new 2D retrieval implementation, we provide recommendations on when 1D models are appropriate and when more complex phase geometries involving multiple TP profiles are required to obtain an unbiased view of tidally locked planetary atmospheres.

preprint2020arXiv

Arid or Cloudy: Characterizing the Atmosphere of the super-Earth 55 Cancri e using High-Resolution Spectroscopy

The nearby super-Earth 55 Cnc e orbits a bright (V = 5.95 mag) star with a period of ~ 18 hours and a mass of ~ 8 Earth masses. Its atmosphere may be water-rich and have a large scale-height, though attempts to characterize it have yielded ambiguous results. Here we present a sensitive search for water and TiO in its atmosphere at high spectral resolution using the Gemini North telescope and the GRACES spectrograph. We combine observations with previous observations from Subaru and CFHT, improving the constraints on the presence of water vapor. We adopt parametric models with an updated planet radius based on recent measurements, and use a cross-correlation technique to maximize sensitivity. Our results are consistent with atmospheres that are cloudy or contain minimal amounts of water and TiO. Using these parametric models, we rule out a water-rich atmosphere (VMR >= 0.1%) with a mean molecular weight of <= 15 g/mol at a 3 sigma confidence level, improving on the previous limit by a significant margin. For TiO, we rule out a mean molecular weight of <= 5 g/mol with a 3 sigma confidence level for a VMR greater than 10^-8; for a VMR of greater than 10^-7, the limit rises to a mean molecular weight of <= 10 g/mol. We can rule out low mean-molecular-weight chemical equilibrium models both including and excluding TiO/VO at very high confidence levels (> 10 sigma). Overall, our results are consistent with an atmosphere with a high mean molecular weight and/or clouds, or no atmosphere.

preprint2020arXiv

Evidence for H$_{2}$ Dissociation and Recombination Heat Transport in the Atmosphere of KELT-9b

Phase curve observations provide an opportunity to study the full energy budgets of exoplanets by quantifying the amount of heat redistributed from their daysides to their nightsides. Theories explaining the properties of phase curves for hot Jupiters have focused on the balance between radiation and dynamics as the primary parameter controlling heat redistribution. However, recent phase curves have shown deviations from the trends that emerge from this theory, which has led to work on additional processes that may affect hot Jupiter energy budgets. One such process, molecular hydrogen dissociation and recombination, can enhance energy redistribution on ultra-hot Jupiters with temperatures above $\sim2000$ K. In order to study the impact of H$_{2}$ dissociation on ultra-hot Jupiters, we present a phase curve of KELT-9b observed with the Spitzer Space Telescope at 4.5 $μ$m. KELT-9b is the hottest known transiting planet, with a 4.5-$μ$m dayside brightness temperature of $4566^{+140}_{-136}$ K and a nightside temperature of $2556^{+101}_{-97}$ K. We observe a phase curve amplitude of $0.609 \pm 0.020$ and a hot spot offset of $18.7^{+2.1}_{-2.3}$ degrees. The observed amplitude is too small to be explained by a simple balance between radiation and advection. General circulation models (GCMs) and an energy balance model that include the effects of H$_{2}$ dissociation and recombination provide a better match to the data. The GCMs, however, predict a maximum hot spot offset of $5$ degrees, which disagrees with our observations at $>5σ$ confidence. This discrepancy may be due to magnetic effects in the planet&#39;s highly ionized atmosphere.

preprint2020arXiv

Observations of Disequilibrium CO Chemistry in the Coldest Brown Dwarfs

Cold brown dwarfs are excellent analogs of widely separated, gas giant exoplanets, and provide insight into the potential atmospheric chemistry and physics we may encounter in objects discovered by future direct imaging surveys. We present a low resolution R $\sim$ 300 $M$-band spectroscopic sequence of seven brown dwarfs with effective temperatures between 750 K and 250 K along with Jupiter. These spectra reveal disequilibrium abundances of carbon monoxide (CO) produced by atmospheric quenching. We use the eddy diffusion coefficient (K$_{zz}$) to estimate the strength of vertical mixing in each object. The K$_{zz}$ values of cooler gaseous objects are close to their theoretical maximum and warmer objects show weaker mixing, likely due to less efficient convective mixing in primarily radiative layers. The CO-derived K$_{zz}$ values imply that disequilibrium phosphine (PH$_{3}$) should be easily observable in all of the brown dwarfs, but none as yet show any evidence for PH$_{3}$ absorption. We find that ammonia is relatively insensitive to atmospheric quenching at these effective temperatures. We are able to improve the fit to WISE 0855&#39;s $M$-band spectrum by including both CO and water clouds in the atmospheric model.

preprint2020arXiv

Smaller than expected bright-spot offsets in Spitzer phase curves of the hot Jupiter Qatar-1b

We present \textit{Spitzer} full-orbit thermal phase curves of the hot Jupiter Qatar-1b, a planet with the same equilibrium temperature---and intermediate surface gravity and orbital period---as the well-studied planets HD 209458b and WASP-43b. We measure secondary eclipse of $0.21 \pm 0.02 \%$ at $3.6~μ$m and $0.30 \pm 0.02 \%$ at $4.5~μ$m, corresponding to dayside brightness temperatures of $1542^{+32}_{-31}$~K and $1557^{+35}_{-36}$~K, respectively, consistent with a vertically isothermal dayside. The respective nightside brightness temperatures are $1117^{+76}_{-71}$~K and $1167^{+69}_{-74}$~K, in line with a trend that hot Jupiters all have similar nightside temperatures. We infer a Bond albedo of $0.12_{-0.16}^{+0.14}$ and a moderate day-night heat recirculation efficiency, similar to HD 209458b. General circulation models for HD 209458b and WASP-43b predict that their bright-spots should be shifted east of the substellar point by tens of degrees, and these predictions were previously confirmed with \textit{Spitzer} full-orbit phase curve observations. The phase curves of Qatar-1b are likewise expected to exhibit eastward offsets. Instead, the observed phase curves are consistent with no offset: $11^{\circ}\pm 7^{\circ}$ at $3.6~μ$m and $-4^{\circ}\pm 7^{\circ}$ at $4.5~μ$m. The discrepancy in circulation patterns between these three otherwise similar planets points to the importance of secondary parameters like rotation rate and surface gravity, and the presence or absence of clouds, in determining atmospheric conditions on hot Jupiters.

preprint2020arXiv

Statistical Characterization of Hot Jupiter Atmospheres using Spitzer&#39;s Secondary Eclipses

We report 78 secondary eclipse depths for a sample of 36 transiting hot Jupiters observed at 3.6- and 4.5 microns using the Spitzer Space Telescope. Our eclipse results for 27 of these planets are new, and include highly irradiated worlds such as KELT-7b, WASP-87b, WASP-76b, and WASP-64b, and important targets for JWST such as WASP-62b. We find that WASP-62b has a slightly eccentric orbit e cos(omega) = 0.00614+/- 0.00064, and we confirm the eccentricity of HAT-P-13b and WASP-14b. The remainder are individually consistent with circular orbits, but we find statistical evidence for eccentricity increasing with orbital period in our range from 1 to 5 days. Our day-side brightness temperatures for the planets yield information on albedo and heat redistribution, following Cowan and Agol (2011). Planets having maximum day side temperatures exceeding ~ 2200K are consistent with zero albedo and distribution of stellar irradiance uniformly over the day-side hemisphere. Our most intriguing result is that we detect a systematic difference between the emergent spectra of these hot Jupiters as compared to blackbodies. The ratio of observed brightness temperatures, Tb(4.5)/Tb(3.6), increases with equilibrium temperature by 100 +/- 24 parts-per-million per Kelvin, over the entire temperature range in our sample (800K to 2500K). No existing model predicts this trend over such a large range of temperature. We suggest that this may be due to a structural difference in the atmospheric temperature profile between the real planetary atmospheres as compared to models.

preprint2020arXiv

The Intrinsic Temperature and Radiative-Convective Boundary Depth in the Atmospheres of Hot Jupiters

In giant planet atmosphere modelling, the intrinsic temperature $T_\mathrm{int}$ and radiative-convective boundary (RCB) are important lower boundary conditions. Often in one-dimensional radiative-convective models and in three-dimensional general circulation models it is assumed that $T_\mathrm{int}$ is similar to that of Jupiter itself, around 100 K, which yields a RCB around 1 kbar for hot Jupiters. In this work, we show that the inflated radii, and hence high specific entropy interiors, of hot Jupiters suggest much higher $T_\mathrm{int}$ values. Assuming the effect is primarily due to current heating (rather than delayed cooling), we derive an equilibrium relation between $T_\mathrm{eq}$ and $T_\mathrm{int}$, showing that the latter can take values as high as 700 K. In response, the RCB moves upward in the atmosphere. Using one-dimensional radiative-convective atmosphere models, we find RCBs of only a few bars, rather than the kilobar typically supposed. This much shallower RCB has important implications for the atmospheric structure, vertical and horizontal circulation, interpretations of phase curves, and the effect of deep cold traps on cloud formation.

preprint2019arXiv

Evaluating Climate Variability Of The Canonical Hot Jupiters Hd 189733b & Hd 209458b Through Multi-epoch Eclipse Observations

Here we present the analysis of multi-epoch secondary eclipse observations of HD 189733b and HD 209458b as a probe of temporal variability in the planetary climate using both Spitzer channels 1 and 2 (3.6 and 4.5 um). Constraining temporal variability will inform models and identify physical processes occurring at either length scales too small to directly observe or at pressure levels that are inaccessible to transit observations. We do not detect statistically significant variability and are able to place useful upper limits on the IR variability amplitudes in these atmospheres. There are very few planets with multi-epoch observations at the required precision to probe variability in dayside emission. The observations considered in this study span several years, providing insight into temporal variability at multiple timescales. In the case of HD 189733b, the best fit eclipse depths for the channel 2 observations exhibit a scatter of 102 ppm about a median depth of 1827 ppm and in channel 1 exhibit a scatter of 88 ppm about a median depth of 1481 ppm. For HD 209458b, the best fit eclipse depths for the channel 2 observations exhibit a scatter of 22 ppm about a median depth of 1406 ppm and in channel 1 exhibit a scatter of 131 ppm about a median depth of 1092 ppm. The precision and scatter in these observations allow us to constrain variability to less than (5.6% and 6.0%) and (12% and 1.6%) for channels (1,2) of HD 189733b and HD 209458b respectively. There is a difference in the best fit eclipse timing compared to the predicted time consistent with an offset hotspot as predicted by GCMs and confirmed in previous phase curve observations.

preprint2019arXiv

Evidence for a Dichotomy in the Interior Structures of Jupiter and Saturn from Helium Phase Separation

We examine the comparative thermal evolution of Jupiter and Saturn applying recent theoretical results for helium&#39;s immiscibility in fluid metallic hydrogen. The redistribution of helium in their interiors proceeds very differently for the two planets. We confirm that based on Jupiter&#39;s atmospheric helium depletion as observed in situ by the Galileo entry probe, Jupiter&#39;s interior helium has differentiated modestly, and we present models reconciling Jupiter&#39;s helium depletion, radius, and heat flow at the solar age. Jupiter&#39;s recently revised Bond albedo implies a lower intrinsic flux for the planet, accommodating more luminosity from helium differentiation such that mildly superadiabatic interiors can satisfy all constraints. The same phase diagram applied to the less massive Saturn predicts dramatic helium differentiation to the degree that Saturn inevitably forms a helium-rich shell or core, an outcome previously proposed by Stevenson & Salpeter and others. The luminosity from Saturn&#39;s helium differentiation is sufficient to extend its cooling time to the solar age, even for adiabatic interiors. This model predicts Saturn&#39;s atmospheric helium to be depleted to Y=0.07+/-0.01, corresponding to a He/H_2 mixing ratio 0.036+/-0.006. We also show that neon differentiation may have contributed to both planets&#39; luminosity in the past. These results demonstrate that Jupiter and Saturn&#39;s thermal evolution can be explained self-consistently with a single physical model, and emphasize that nontrivial helium distributions should be considered in future models for Saturn&#39;s internal structure and dynamo.

preprint2019arXiv

Investigating Trends in Atmospheric Compositions of Cool Gas Giant Planets Using Spitzer Secondary Eclipses

We present new 3.6 and 4.5 micron secondary eclipse measurements for five cool (less than approximately 1000 K) transiting gas giant planets: HAT-P-15b, HAT-P-17b, HAT-P-18b, HAT-P-26b, and WASP-69b. We detect eclipses in at least one bandpass for all planets except HAT-P-15b. We confirm and refine the orbital eccentricity of HAT-P-17b, which is also the only planet in our sample with a known outer companion. We compare our measured eclipse depths in these two bands, which are sensitive to the relative abundances of methane versus carbon monoxide and carbon dioxide, respectively, to predictions from 1D atmosphere models for each planet. For planets with hydrogen-dominated atmospheres and equilibrium temperatures cooler than approximately 1000 K, this ratio should vary as a function of both atmospheric metallicity and the carbon-to-oxygen ratio. For HAT-P-26b, our observations are in good agreement with the low atmospheric metallicity inferred from transmission spectroscopy. We find that all four of the planets with detected eclipses are best matched by models with relatively efficient circulation of energy to the nightside. We see no evidence for a solar-system-like correlation between planet mass and atmospheric metallicity, but instead identify a potential (1.9 sigma) correlation between the inferred methane/(carbon monoxide + carbon dioxide) ratio and stellar metallicity. Our ability to characterize this potential trend is limited by the relatively large uncertainties in the stellar metallicity values. Our observations provide a first look at the brightness of these planets at wavelengths accessible to the James Webb Space Telescope, which will be able to resolve individual methane, carbon monoxide, and carbon dioxide bands and provide much stronger constraints on their atmospheric compositions.

preprint2019arXiv

Saturn&#39;s Probable Interior: An Exploration of Saturn&#39;s Potential Interior Density Structures

The gravity field of a giant planet is typically our best window into its interior structure and composition. Through comparison of a model planet&#39;s calculated gravitational potential with the observed potential, inferences can be made about interior quantities, including possible composition and the existence of a core. Necessarily, a host of assumptions go into such calculations, making every inference about a giant planet&#39;s structure strongly model dependent. In this work we present a more general picture by setting Saturn&#39;s gravity field, as measured during the \emph{Cassini} Grand Finale, as a likelihood function driving a Markov-Chain Monte Carlo exploration of the possible interior density profiles. The result is a posterior distribution of the interior structure that is not tied to assumed composition, thermal state, or material equations of state. Constraints on interior structure derived in this Bayesian framework are necessarily less informative, but are also less biased and more general. These empirical and probabilistic constraints on the density structure are our main data product which we archive for continued analysis. We find that the outer half of Saturn&#39;s radius is relatively well constrained, and we interpret our findings as suggesting a significant metal enrichment, in line with atmospheric abundances from remote sensing. As expected, the inner half of Saturn&#39;s radius is less well-constrained by gravity, but we generally find solutions that include a significant density enhancement, which can be interpreted as a core, although this core is often lower in density and larger in radial extent than typically found by standard models. This is consistent with a dilute core and/or composition gradients.

preprint2019arXiv

The Featureless Transmission Spectra of Two Super-Puff Planets

The Kepler mission revealed a class of planets known as &#39;&#39;super-puffs,&#39;&#39; with masses only a few times larger than Earth&#39;s but radii larger than Neptune, giving them very low mean densities. All three of the known planets orbiting the young solar-type star Kepler 51 are super-puffs. The Kepler 51 system thereby provides an opportunity for a comparative study of the structures and atmospheres of this mysterious class of planets, which may provide clues about their formation and evolution. We observed two transits each of Kepler 51b and 51d with the Wide Field Camera 3 (WFC3) on the Hubble Space Telescope. Combining new WFC3 transit times with re-analyzed Kepler data and updated stellar parameters, we confirmed that all three planets have densities lower than 0.1 g/cm$^{3}$. We measured the WFC3 transmission spectra to be featureless between 1.15 and 1.63 $μ$m, ruling out any variations greater than 0.6 scale heights (assuming a H/He dominated atmosphere), thus showing no significant water absorption features. We interpreted the flat spectra as the result of a high-altitude aerosol layer (pressure $<$3 mbar) on each planet. Adding this new result to the collection of flat spectra that have been observed for other sub-Neptune planets, we find support for one of the two hypotheses introduced by Crossfield and Kreidberg (2017), that planets with cooler equilibrium temperatures have more high-altitude aerosols. We strongly disfavor their other hypothesis that the H/He mass fraction drives the appearance of large amplitude transmission features.

preprint2019arXiv

The Precision of Mass Measurements Required for Robust Atmospheric Characterization of Transiting Exoplanets

Two of TESS&#39;s major science goals are to measure masses for 50 planets smaller than 4 Earth radii and to discover high-quality targets for atmospheric characterization efforts. It is important that these two goals are linked by quantifying what precision of mass constraint is required to yield robust atmospheric properties of planets. Here, we address this by conducting retrievals on simulated JWST transmission spectra under various assumptions for the degree of uncertainty in the planet&#39;s mass for a representative population of seven planets ranging from terrestrials to warm Neptunes to hot Jupiters. Only for the cloud-free, low metallicity gas giants are we able to infer exoplanet mass from transmission spectroscopy alone, to ~10% accuracy. For low metallicity cases (<4xSolar) we are able to accurately constrain atmospheric properties without prior knowledge of the planet&#39;s mass. For all other cases (including terrestrial-like planets), atmospheric properties can only be inferred with a mass precision of better than ~50%. At this level, though, the widths of the posterior distributions of the atmospheric properties are dominated by the uncertainties in mass. With a precision of ~20%, the widths of the posterior distributions are dominated by the spectroscopic data quality. Therefore, as a rule-of-thumb, we recommend: a ~50% mass precision for initial atmospheric characterization and a ~20% mass precision for more detailed atmospheric analyses.