-
Two-Dimensional Eclipse Mapping of the Hot Jupiter WASP-43b with JWST MIRI/LRS
Authors:
Mark Hammond,
Taylor J. Bell,
Ryan C. Challener,
Neil T. Lewis,
Megan Weiner Mansfield,
Isaac Malsky,
Emily Rauscher,
Jacob L. Bean,
Ludmila Carone,
João M. Mendonça,
Lucas Teinturier,
Xianyu Tan,
Nicolas Crouzet,
Laura Kreidberg,
Giuseppe Morello,
Vivien Parmentier,
Jasmina Blecic,
Jean-Michel Désert,
Christiane Helling,
Pierre-Olivier Lagage,
Karan Molaverdikhani,
Matthew C. Nixon,
Benjamin V. Rackham,
Jingxuan Yang
Abstract:
We present eclipse maps of the two-dimensional thermal emission from the dayside of the hot Jupiter WASP-43b, derived from an observation of a phase curve with the JWST MIRI/LRS instrument. The observed eclipse shapes deviate significantly from those expected for a planet emitting uniformly over its surface. We fit a map to this deviation, constructed from spherical harmonics up to order…
▽ More
We present eclipse maps of the two-dimensional thermal emission from the dayside of the hot Jupiter WASP-43b, derived from an observation of a phase curve with the JWST MIRI/LRS instrument. The observed eclipse shapes deviate significantly from those expected for a planet emitting uniformly over its surface. We fit a map to this deviation, constructed from spherical harmonics up to order $\ell_{\rm max}=2$, alongside the planetary, orbital, stellar, and systematic parameters. This yields a map with a meridionally-averaged eastward hot-spot shift of $(7.75 \pm 0.36)^{\circ}$, with no significant degeneracy between the map and the additional parameters. We show the latitudinal and longitudinal contributions of the day-side emission structure to the eclipse shape, finding a latitudinal signal of $\sim$200 ppm and a longitudinal signal of $\sim$250 ppm. To investigate the sensitivity of the map to the method, we fix the non-mapping parameters and derive an "eigenmap" fitted with an optimised number of orthogonal phase curves, which yields a similar map to the $\ell_{\rm max}=2$ map. We also fit a map up to $\ell_{\rm max}=3$, which shows a smaller hot-spot shift, with a larger uncertainty. These maps are similar to those produced by atmospheric simulations. We conclude that there is a significant mapping signal which constrains the spherical harmonic components of our model up to $\ell_{\rm max}=2$. Alternative mapping models may derive different structures with smaller-scale features; we suggest that further observations of WASP-43b and other planets will drive the development of more robust methods and more accurate maps.
△ Less
Submitted 25 April, 2024;
originally announced April 2024.
-
The XUV-driven escape of the planets around TOI-431 & $ν^2$ Lupi
Authors:
George W. King,
Lía R. Corrales,
Jorge Fernández Fernández,
Peter J. Wheatley,
Isaac Malsky,
Ares Osborn,
David Armstrong
Abstract:
One of the leading mechanisms invoked to explain the existence of the radius valley is atmospheric mass loss driven by X-ray and extreme-ultraviolet irradiation, with this process stripping the primordial envelopes of young, small planets to produce the observed bimodal distribution. We present an investigation into the TOI-431 and $ν^2$ Lupi planetary systems, both of which host planets either si…
▽ More
One of the leading mechanisms invoked to explain the existence of the radius valley is atmospheric mass loss driven by X-ray and extreme-ultraviolet irradiation, with this process stripping the primordial envelopes of young, small planets to produce the observed bimodal distribution. We present an investigation into the TOI-431 and $ν^2$ Lupi planetary systems, both of which host planets either side of the radius valley, to determine if their architectures are consistent with evolution by the XUV mechanism. With $\textit{XMM-Newton}$, we measure the current X-ray flux of each star, and see evidence for a stellar flare in the TOI-431 observations. We then simulate the evolution of all of the transiting planets across the two systems in response to the high-energy irradiation over their lifetimes. We use the measured X-ray fluxes as an anchor point for the XUV time evolution in our simulations, and employ several different models of estimating mass loss rates. While the simulations for TOI-431b encountered a problem with the initial calculated radii, we estimate a likely short ($\sim$ Myr) timespan for primordial envelope removal using reasonable assumptions for the initial planet. $ν^2$ Lupi b is likely harder to strip, but is achieved in a moderate fraction of our simulations. None of our simulations stripped any of the lower density planets of their envelope, in line with prediction. We conclude that both systems are consistent with expectations for generation of the radius valley through XUV photoevaporation.
△ Less
Submitted 16 April, 2024;
originally announced April 2024.
-
Nightside clouds and disequilibrium chemistry on the hot Jupiter WASP-43b
Authors:
Taylor J. Bell,
Nicolas Crouzet,
Patricio E. Cubillos,
Laura Kreidberg,
Anjali A. A. Piette,
Michael T. Roman,
Joanna K. Barstow,
Jasmina Blecic,
Ludmila Carone,
Louis-Philippe Coulombe,
Elsa Ducrot,
Mark Hammond,
João M. Mendonça,
Julianne I. Moses,
Vivien Parmentier,
Kevin B. Stevenson,
Lucas Teinturier,
Michael Zhang,
Natalie M. Batalha,
Jacob L. Bean,
Björn Benneke,
Benjamin Charnay,
Katy L. Chubb,
Brice-Olivier Demory,
Peter Gao
, et al. (58 additional authors not shown)
Abstract:
Hot Jupiters are among the best-studied exoplanets, but it is still poorly understood how their chemical composition and cloud properties vary with longitude. Theoretical models predict that clouds may condense on the nightside and that molecular abundances can be driven out of equilibrium by zonal winds. Here we report a phase-resolved emission spectrum of the hot Jupiter WASP-43b measured from 5…
▽ More
Hot Jupiters are among the best-studied exoplanets, but it is still poorly understood how their chemical composition and cloud properties vary with longitude. Theoretical models predict that clouds may condense on the nightside and that molecular abundances can be driven out of equilibrium by zonal winds. Here we report a phase-resolved emission spectrum of the hot Jupiter WASP-43b measured from 5-12 $μ$m with JWST's Mid-Infrared Instrument (MIRI). The spectra reveal a large day-night temperature contrast (with average brightness temperatures of 1524$\pm$35 and 863$\pm$23 Kelvin, respectively) and evidence for water absorption at all orbital phases. Comparisons with three-dimensional atmospheric models show that both the phase curve shape and emission spectra strongly suggest the presence of nightside clouds which become optically thick to thermal emission at pressures greater than ~100 mbar. The dayside is consistent with a cloudless atmosphere above the mid-infrared photosphere. Contrary to expectations from equilibrium chemistry but consistent with disequilibrium kinetics models, methane is not detected on the nightside (2$σ$ upper limit of 1-6 parts per million, depending on model assumptions).
△ Less
Submitted 23 January, 2024;
originally announced January 2024.
-
A Direct Comparison between the use of Double Gray and Multiwavelength Radiative Transfer in a General Circulation Model with and without Radiatively Active Clouds
Authors:
Isaac Malsky,
Emily Rauscher,
Michael T. Roman,
Elspeth K. H. Lee,
Hayley Beltz,
Arjun Savel,
Eliza M. R. Kempton,
L. Cinque
Abstract:
Inhomogeneous cloud formation and wavelength-dependent phenomena are expected to shape hot Jupiter atmospheres. We present a General Circulation Model (GCM) with multiwavelength "picket fence" radiative transfer and radiatively active, temperature dependent clouds, and compare the results to a double gray routine. The double gray method inherently fails to model polychromatic effects in hot Jupite…
▽ More
Inhomogeneous cloud formation and wavelength-dependent phenomena are expected to shape hot Jupiter atmospheres. We present a General Circulation Model (GCM) with multiwavelength "picket fence" radiative transfer and radiatively active, temperature dependent clouds, and compare the results to a double gray routine. The double gray method inherently fails to model polychromatic effects in hot Jupiter atmospheres, while picket fence captures these non-gray aspects and performs well compared to fully wavelength-dependent methods. We compare both methods with radiatively active clouds and cloud-free models, assessing the limitations of the double gray method. Although there are broad similarities, the picket fence models have larger day-night side temperature differences, non-isothermal upper atmospheres, and multiwavelength effects in the presence of radiatively active clouds. We model the well-known hot Jupiters HD 189733 b and HD 209458 b. For the hotter HD 209458 b, the picket fence method prevents clouds from thermostating dayside temperatures, resulting in hotter upper atmospheres and the dissipation of dayside clouds. Differences in the temperature structures are then associated with nuanced differences in the circulation patterns and clouds. Models of the cooler HD 189733 b have global cloud coverage, regardless of radiative transfer scheme, whereas there are larger differences in the models of HD 209458 b, particularly in the extent of the partial cloud coverage on its dayside. This results in minor changes to the thermal and reflected light phase curves of HD 189733 b, but more significant differences for the picket fence and double gray versions of HD 209458 b.
△ Less
Submitted 26 February, 2024; v1 submitted 2 November, 2023;
originally announced November 2023.
-
Helium Enhanced Planets Along the Upper Edge of the Radius Valley
Authors:
Isaac Malsky,
Leslie Rogers,
Eliza M. R. Kempton,
Nadejda Marounina
Abstract:
The low mean densities of sub-Neptunes imply that they formed within a few million years and accreted primordial envelopes. Because these planets receive a total X-ray and extreme ultra-violet flux that is comparable to the gravitational binding energy of their envelopes, their primordial hydrogen-helium atmospheres are susceptible to mass loss. Models of photoevaporating sub-Neptunes have so far…
▽ More
The low mean densities of sub-Neptunes imply that they formed within a few million years and accreted primordial envelopes. Because these planets receive a total X-ray and extreme ultra-violet flux that is comparable to the gravitational binding energy of their envelopes, their primordial hydrogen-helium atmospheres are susceptible to mass loss. Models of photoevaporating sub-Neptunes have so far assumed that envelope compositions remain constant over time. However, preferential loss of atmospheric hydrogen has the potential to change their compositions. Here, by modeling the thermal and compositional evolution of sub-Neptunes undergoing atmospheric escape with diffusive separation between hydrogen and helium, we show that planets with radii between 1.6 and 2.5 that of Earth can become helium-enhanced from billions of years of photoevaporation, obtaining helium mass fractions in excess of 40%. Atmospheric helium enhancement can be detected through transmission spectra, providing a novel observational test for whether atmospheric escape creates the radius valley.
△ Less
Submitted 27 September, 2023;
originally announced September 2023.
-
A reflective, metal-rich atmosphere for GJ 1214b from its JWST phase curve
Authors:
Eliza M. -R. Kempton,
Michael Zhang,
Jacob L. Bean,
Maria E. Steinrueck,
Anjali A. A. Piette,
Vivien Parmentier,
Isaac Malsky,
Michael T. Roman,
Emily Rauscher,
Peter Gao,
Taylor J. Bell,
Qiao Xue,
Jake Taylor,
Arjun B. Savel,
Kenneth E. Arnold,
Matthew C. Nixon,
Kevin B. Stevenson,
Megan Mansfield,
Sarah Kendrew,
Sebastian Zieba,
Elsa Ducrot,
Achrène Dyrek,
Pierre-Olivier Lagage,
Keivan G. Stassun,
Gregory W. Henry
, et al. (8 additional authors not shown)
Abstract:
There are no planets intermediate in size between Earth and Neptune in our Solar System, yet these objects are found around a substantial fraction of other stars. Population statistics show that close-in planets in this size range bifurcate into two classes based on their radii. It is hypothesized that the group with larger radii (referred to as "sub-Neptunes") is distinguished by having hydrogen-…
▽ More
There are no planets intermediate in size between Earth and Neptune in our Solar System, yet these objects are found around a substantial fraction of other stars. Population statistics show that close-in planets in this size range bifurcate into two classes based on their radii. It is hypothesized that the group with larger radii (referred to as "sub-Neptunes") is distinguished by having hydrogen-dominated atmospheres that are a few percent of the total mass of the planets. GJ 1214b is an archetype sub-Neptune that has been observed extensively using transmission spectroscopy to test this hypothesis. However, the measured spectra are featureless, and thus inconclusive, due to the presence of high-altitude aerosols in the planet's atmosphere. Here we report a spectroscopic thermal phase curve of GJ 1214b obtained with JWST in the mid-infrared. The dayside and nightside spectra (average brightness temperatures of 553 $\pm$ 9 and 437 $\pm$ 19 K, respectively) each show >3$σ$ evidence of absorption features, with H$_2$O as the most likely cause in both. The measured global thermal emission implies that GJ 1214b's Bond albedo is 0.51 $\pm$ 0.06. Comparison between the spectroscopic phase curve data and three-dimensional models of GJ 1214b reveal a planet with a high metallicity atmosphere blanketed by a thick and highly reflective layer of clouds or haze.
△ Less
Submitted 10 May, 2023;
originally announced May 2023.
-
The Hazy and Metal-Rich Atmosphere of GJ 1214 b Constrained by Near and Mid-Infrared Transmission Spectroscopy
Authors:
Peter Gao,
Anjali A. A. Piette,
Maria E. Steinrueck,
Matthew C. Nixon,
Michael Zhang,
Eliza M. R. Kempton,
Jacob L. Bean,
Emily Rauscher,
Vivien Parmentier,
Natasha E. Batalha,
Arjun B. Savel,
Kenneth E. Arnold,
Michael T. Roman,
Isaac Malsky,
Jake Taylor
Abstract:
The near-infrared transmission spectrum of the warm sub-Neptune exoplanet GJ 1214 b has been observed to be flat and featureless, implying a high metallicity atmosphere with abundant aerosols. Recent JWST MIRI LRS observations of a phase curve of GJ 1214 b showed that its transmission spectrum is flat out into the mid-infrared. In this paper, we use the combined near- and mid-infrared transmission…
▽ More
The near-infrared transmission spectrum of the warm sub-Neptune exoplanet GJ 1214 b has been observed to be flat and featureless, implying a high metallicity atmosphere with abundant aerosols. Recent JWST MIRI LRS observations of a phase curve of GJ 1214 b showed that its transmission spectrum is flat out into the mid-infrared. In this paper, we use the combined near- and mid-infrared transmission spectrum of GJ 1214 b to constrain its atmospheric composition and aerosol properties. We generate a grid of photochemical haze models using an aerosol microphysics code for a number of background atmospheres spanning metallicities from 100 to 1000 $\times$ solar, as well as a steam atmosphere scenario. The flatness of the combined data set largely rules out atmospheric metallicities $\leq$300 $\times$ solar due to their large corresponding molecular feature amplitudes, preferring values $\geq$1000 $\times$ solar and column haze production rates $\geq$10$^{-10}$ g cm$^{-2}$ s$^{-1}$. The steam atmosphere scenario with similarly high haze production rates also exhibit sufficiently small molecular features to be consistent with the transmission spectrum. These compositions imply that atmospheric mean molecular weights $\geq$15 g mol$^{-1}$ are needed to fit the data. Our results suggest that haze production is highly efficient on GJ 1214 b and could involve non-hydrocarbon, non-nitrogen haze precursors. Further characterization of GJ 1214 b's atmosphere would likely require multiple transits and eclipses using JWST across the near and mid-infrared, potentially complemented by groundbased high resolution transmission spectroscopy.
△ Less
Submitted 9 May, 2023;
originally announced May 2023.
-
Magnetic Effects and 3D Structure in Theoretical High-Resolution Transmission Spectra of Ultrahot Jupiters: the Case of WASP-76b
Authors:
Hayley Beltz,
Emily Rauscher,
Eliza Kempton,
Isaac Malsky,
Arjun Savel
Abstract:
High resolution spectroscopy has allowed for unprecedented levels of atmospheric characterization, especially for the hottest gas giant exoplanets known as ultrahot Jupiters (UHJs). High-resolution spectra are sensitive to 3D effects, making complex 3D atmospheric models important for interpreting data. Moreover, these planets are expected to host magnetic fields that will shape their resulting at…
▽ More
High resolution spectroscopy has allowed for unprecedented levels of atmospheric characterization, especially for the hottest gas giant exoplanets known as ultrahot Jupiters (UHJs). High-resolution spectra are sensitive to 3D effects, making complex 3D atmospheric models important for interpreting data. Moreover, these planets are expected to host magnetic fields that will shape their resulting atmospheric circulation patterns, but little modeling work has been done to investigate these effects. In this paper, we generate high-resolution transmission spectra from General Circulation Models for the canonical UHJ WASP-76b with three different magnetic treatments in order to understand the influence of magnetic forces on the circulation. In general, spectra from all models have increasingly blueshifted net Doppler shifts as transit progresses, but we find that the differing temperature and wind fields in the upper atmospheres of these models result in measurable differences. We find that magnetic effects may be contributing to the unusual trends previously seen in transmission for this planet. Our $B=3$ Gauss active drag model in particular shows unique trends not found in the models with simpler or no magnetic effects. The net Doppler shifts are additionally influenced by the dominant opacity sources in each wavelength range considered, as each species probes different regions of the atmosphere and are sensitive to spatial differences in the circulation. This work highlights the ongoing need for models of planets in this temperature regime to consider both 3D and magnetic effects when interpreting high resolution transmission spectra.
△ Less
Submitted 22 May, 2023; v1 submitted 27 February, 2023;
originally announced February 2023.
-
Diagnosing limb asymmetries in hot and ultra-hot Jupiters with high-resolution transmission spectroscopy
Authors:
Arjun B. Savel,
Eliza M. -R. Kempton,
Emily Rauscher,
Thaddeus D. Komacek,
Jacob L. Bean,
Matej Malik,
Isaac Malsky
Abstract:
Due to their likely tidally synchronized nature, (ultra)hot Jupiter atmospheres should experience strongly spatially heterogeneous instellation. The large irradiation contrast and resulting atmospheric circulation induce temperature and chemical gradients that can produce asymmetries across the eastern and western limbs of these atmospheres during transit. By observing an (ultra)hot Jupiter's tran…
▽ More
Due to their likely tidally synchronized nature, (ultra)hot Jupiter atmospheres should experience strongly spatially heterogeneous instellation. The large irradiation contrast and resulting atmospheric circulation induce temperature and chemical gradients that can produce asymmetries across the eastern and western limbs of these atmospheres during transit. By observing an (ultra)hot Jupiter's transmission spectrum at high spectral resolution, these asymmetries can be recovered -- namely through net Doppler shifts originating from the exoplanet's atmosphere yielded by cross-correlation analysis. Given the range of mechanisms at play, identifying the underlying cause of observed asymmetry is nontrivial. In this work, we explore sources and diagnostics of asymmetries in high-resolution cross-correlation spectroscopy of hot and ultra-hot Jupiters using both parameterized and self-consistent atmospheric models. If an asymmetry is observed, we find that it can be difficult to attribute it to equilibrium chemistry gradients because many other processes can produce asymmetries. Identifying a molecule that is chemically stable over the temperature range of a planetary atmosphere can help establish a ``baseline'' to disentangle the various potential causes of limb asymmetries observed in other species. We identify CO as an ideal molecule, given its stability over nearly the entirety of the ultra-hot Jupiter temperature range. Furthermore, we find that if limb asymmetry is due to morning terminator clouds, blueshifts for a number of species should decrease during transit. Finally, by comparing our forward models to Kesseli et al. (2022), we demonstrate that binning high-resolution spectra into two phase bins provides a desirable trade-off between maintaining signal to noise and resolving asymmetries.
△ Less
Submitted 4 January, 2023;
originally announced January 2023.
-
A Lack of Variability Between Repeated Spitzer Phase Curves of WASP-43b
Authors:
Matthew M. Murphy,
Thomas G. Beatty,
Michael T. Roman,
Isaac Malsky,
Alex Wingate,
Grace Ochs,
L. Cinque,
Hayley Beltz,
Emily Rauscher,
Eliza M. -R. Kempton,
Kevin B. Stevenson
Abstract:
Though the global atmospheres of hot Jupiters have been extensively studied using phase curve observations, the level of time variability in these data is not well constrained. To investigate possible time variability in a planetary phase curve, we observed two full-orbit phase curves of the hot Jupiter WASP-43b at 4.5 microns using the Spitzer Space Telescope, and reanalyzed a previous 4.5 micron…
▽ More
Though the global atmospheres of hot Jupiters have been extensively studied using phase curve observations, the level of time variability in these data is not well constrained. To investigate possible time variability in a planetary phase curve, we observed two full-orbit phase curves of the hot Jupiter WASP-43b at 4.5 microns using the Spitzer Space Telescope, and reanalyzed a previous 4.5 micron phase curve from Stevenson et al. (2017). We find no significant time variability between these three phase curves, which span timescales of weeks to years. The three observations are best fit by a single phase curve with an eclipse depth of 3907 +- 85 ppm, a dayside-integrated brightness temperature of 1479 +- 13 K, a nightside-integrated brightness temperature of 755 +- 46 K, and an eastward-shifted peak of 10.4 +- 1.8 degrees. To model our observations, we performed 3D general circulation model simulations of WASP-43b with simple cloud models of various vertical extents. In comparing these simulations to our observations, we find that WASP-43b likely has a cloudy nightside that transitions to a relatively cloud-free dayside. We estimate that any change in WASP-43bs vertical cloud thickness of more than three pressure scale heights is inconsistent with our observed upper limit on variation. These observations, therefore, indicate that WASP-43bs clouds are stable in their vertical and spatial extent over timescales up to several years. These results strongly suggest that atmospheric properties derived from previous, single Spitzer phase curve observations of hot Jupiters likely show us the equilibrium properties of these atmospheres.
△ Less
Submitted 15 February, 2023; v1 submitted 6 December, 2022;
originally announced December 2022.
-
CAMEMBERT: A Mini-Neptunes GCM Intercomparison, Protocol Version 1.0. A CUISINES Model Intercomparison Project
Authors:
Duncan A. Christie,
Elspeth K. H. Lee,
Hamish Innes,
Pascal A. Noti,
Benjamin Charnay,
Thomas J. Fauchez,
Nathan J. Mayne,
Russell Deitrick,
Feng Ding,
Jennifer J. Greco,
Mark Hammond,
Isaac Malsky,
Avi Mandell,
Emily Rauscher,
Michael T. Roman,
Denis E. Sergeev,
Linda Sohl,
Maria E. Steinrueck,
Martin Turbet,
Eric T. Wolf,
Maria Zamyatina,
Ludmila Carone
Abstract:
With an increased focus on the observing and modelling of mini-Neptunes, there comes a need to better understand the tools we use to model their atmospheres. In this paper, we present the protocol for the CAMEMBERT (Comparing Atmospheric Models of Extrasolar Mini-neptunes Building and Envisioning Retrievals and Transits) project, an intercomparison of general circulation models (GCMs) used by the…
▽ More
With an increased focus on the observing and modelling of mini-Neptunes, there comes a need to better understand the tools we use to model their atmospheres. In this paper, we present the protocol for the CAMEMBERT (Comparing Atmospheric Models of Extrasolar Mini-neptunes Building and Envisioning Retrievals and Transits) project, an intercomparison of general circulation models (GCMs) used by the exoplanetary science community to simulate the atmospheres of mini-Neptunes. We focus on two targets well studied both observationally and theoretically with planned JWST Cycle 1 observations: the warm GJ~1214b and the cooler K2-18b. For each target, we consider a temperature-forced case, a clear sky dual-grey radiative transfer case, and a clear sky multi band radiative transfer case, covering a range of complexities and configurations where we know differences exist between GCMs in the literature. This paper presents all the details necessary to participate in the intercomparison, with the intention of presenting the results in future papers. Currently, there are eight GCMs participating (ExoCAM, Exo-FMS, FMS PCM, Generic PCM, MITgcm, RM-GCM, THOR, and the UM), and membership in the project remains open. Those interested in participating are invited to contact the authors.
△ Less
Submitted 8 November, 2022;
originally announced November 2022.
-
Magnetic Drag and 3-D Effects in Theoretical High-Resolution Emission Spectra of Ultrahot Jupiters: the Case of WASP-76b
Authors:
Hayley Beltz,
Emily Rauscher,
Eliza M. -R Kempton,
Isaac Malsky,
Grace Ochs,
Mireya Arora,
Arjun Savel
Abstract:
Ultrahot Jupiters are ideal candidates to explore with high-resolution emission spectra. Detailed theoretical studies are necessary to investigate the range of spectra we can expect to see from these objects throughout their orbit, because of the extreme temperature and chemical longitudinal gradients that exist across day and nightside regions. Using previously published 3D GCM models of WASP-76b…
▽ More
Ultrahot Jupiters are ideal candidates to explore with high-resolution emission spectra. Detailed theoretical studies are necessary to investigate the range of spectra we can expect to see from these objects throughout their orbit, because of the extreme temperature and chemical longitudinal gradients that exist across day and nightside regions. Using previously published 3D GCM models of WASP-76b with different treatments of magnetic drag, we post-process the 3D atmospheres to generate high-resolution emission spectra for two wavelength ranges and throughout the planet's orbit. We find that the high-resolution emission spectra vary strongly as a function of phase, at times showing emission features, absorption features, or both, which are a direct result of the 3D structure of the planet. At phases exhibiting both emission and absorption features, the Doppler shift differs in direction between the two spectral features, making them differentiable instead of canceling each other out. Through the use of cross-correlation, we find different patterns in net Doppler shift for models with different treatments of drag: the nightside spectra show opposite signs in their Doppler shift, while the dayside phases have a reversal in the trend of net shift with phase. Finally, we caution researchers from using a single spectral template throughout the planet's orbit; this can bias the corresponding net Doppler shift returned, as it can pick up on a bright region on the edge of the planet disk that is highly red- or blue-shifted.
△ Less
Submitted 19 September, 2022; v1 submitted 27 April, 2022;
originally announced April 2022.
-
A Multi-Planet System's Sole Super-Puff: Exploring Allowable Physical Parameters for the Cold Super-Puff HIP 41378 f
Authors:
Michelle Belkovski,
Juliette Becker,
Alex Howe,
Isaac Malsky,
Konstantin Batygin
Abstract:
The census of known exoplanets exhibits a variety of physical parameters, including densities that are measured to span the range from less dense than styrofoam to more dense than iron. These densities represent a large diversity of interior structures. Despite this staggering diversity, recent analyses have shown that the densities of planets that orbit a common star exhibit remarkable uniformity…
▽ More
The census of known exoplanets exhibits a variety of physical parameters, including densities that are measured to span the range from less dense than styrofoam to more dense than iron. These densities represent a large diversity of interior structures. Despite this staggering diversity, recent analyses have shown that the densities of planets that orbit a common star exhibit remarkable uniformity. A fascinating exception to this is the system HIP 41378 (also known as K2-93), which contains a super-puff planet, HIP 41378 f, as well as several planets with more typical bulk densities. The range of densities in this system begs the question of what physical processes are responsible for the disparate planetary structures in this system. In this paper, we consider how the densities of the planets in the HIP 41378 system would have changed over time as the host star evolved and the planets' atmospheres were subsequently affected by the evolving insolation level. We also present a range of allowable core masses for HIP 41378 f based on the measured planet parameters, and comment on the feasibility of the proposed existence of planetary rings around HIP 41378 f as an explanation for its current low density.
△ Less
Submitted 31 March, 2022;
originally announced March 2022.
-
Irradiation-driven escape of primordial planetary atmospheres II. Evaporation efficiency of sub-Neptunes through hot Jupiters
Authors:
Andrea Caldiroli,
Francesco Haardt,
Elena Gallo,
Riccardo Spinelli,
Isaac Malsky,
Emily Rauscher
Abstract:
Making use of the publicly available 1D photoionization hydrodynamics code ATES we set out to investigate the combined effects of planetary gravitational potential energy ($φ_p\equiv GM_p/R_p$) and stellar X-ray and Extreme Ultraviolet (XUV) irradiation ($F_{\rm XUV}$) on the evaporation efficiency ($η$) of moderately-to-highly irradiated gaseous planets, from sub-Neptunes through hot Jupiters. We…
▽ More
Making use of the publicly available 1D photoionization hydrodynamics code ATES we set out to investigate the combined effects of planetary gravitational potential energy ($φ_p\equiv GM_p/R_p$) and stellar X-ray and Extreme Ultraviolet (XUV) irradiation ($F_{\rm XUV}$) on the evaporation efficiency ($η$) of moderately-to-highly irradiated gaseous planets, from sub-Neptunes through hot Jupiters. We show that the (known) existence of a threshold potential above which energy-limited escape (i.e., $η\simeq 1$) is unattainable can be inferred analytically. For $\log φ_p\gtrsim \log φ_p^{\rm thr}\approx [12.9-13.2]$ (in cgs units), most of the energy absorption occurs where the average kinetic energy acquired by the ions through photo-electron collisions is insufficient for escape. This causes the evaporation efficiency to plummet with increasing $φ_p$,. Whether or not planets with $φ_p\lesssim φ_p^{\rm thr}$ exhibit energy-limited outflows is regulated primarily by the stellar irradiation level. Specifically, for low-gravity planets, above $F_{\rm XUV}\simeq 10^{4-5}$ erg cm$^{-2}$s$^{-1}$ Ly$α$ losses overtake adiabatic and advective cooling and the evaporation efficiency of low-gravity planets drops below the energy-limited approximation, albeit remaining largely independent of $φ_p$Further, we show that whereas $η$ increases as $F_{\rm XUV}$ increases for planets above $φ^{\rm thr}_p$, the opposite is true for low-gravity planets. This behavior can be understood by examining the relative fractional contributions of advective and radiative losses as a function of atmospheric temperature. This novel framework enables a reliable, physically motivated prediction of the expected evaporation efficiency for a given planetary system; an analytical approximation of the best-fitting $η$ is given in the appendix.
△ Less
Submitted 6 April, 2022; v1 submitted 1 December, 2021;
originally announced December 2021.
-
Modeling the high-resolution emission spectra of clear and cloudy non-transiting hot Jupiters
Authors:
Isaac Malsky,
Emily Rauscher,
Eliza M. R. Kempton,
Michael Roman,
Deryl Long,
Caleb K. Harada
Abstract:
The advent of high-resolution spectroscopy as a method for exoplanet atmospheric characterization has expanded our capability to study non-transiting planets, increasing the number of planets accessible for observation. Many of the most favorable targets for atmospheric characterization are hot Jupiters, where we expect large spatial variation in physical conditions such as temperature, wind speed…
▽ More
The advent of high-resolution spectroscopy as a method for exoplanet atmospheric characterization has expanded our capability to study non-transiting planets, increasing the number of planets accessible for observation. Many of the most favorable targets for atmospheric characterization are hot Jupiters, where we expect large spatial variation in physical conditions such as temperature, wind speed, and cloud coverage, making viewing geometry important. Three-dimensional models have generally simulated observational properties of hot Jupiters assuming edge-on viewing, which neglects planets without near edge-on orbits. As the first investigation of how orbital inclination manifests in high-resolution emission spectra, we use a General Circulation Model to simulate the atmospheric structure of Upsilon Andromedae b, a non-transiting hot Jupiter. In order to accurately capture scattering from clouds, we implement a generalized two-stream radiative transfer routine for inhomogeneous multiple scattering atmospheres. We compare models with and without clouds, as cloud coverage intensifies spatial variations. Cloud coverage increases both the net Doppler shifts and the variation of the continuum flux amplitude over the course of the planet's orbit. As orbital inclination decreases, four key features also decrease in both the clear and cloudy models: 1) the average continuum flux level, 2) the amplitude of the variation in continuum with orbital phase, 3) net Doppler shifts of spectral lines, and 4) Doppler broadening in the spectra. Models capable of treating inhomogeneous cloud coverage and different viewing geometries are critical in understanding high-resolution emission spectra, enabling an additional avenue to investigate these extreme atmospheres.
△ Less
Submitted 11 October, 2021;
originally announced October 2021.
-
Irradiation-driven escape of primordial planetary atmospheres I. The ATES photoionization hydrodynamics code
Authors:
Andrea Caldiroli,
Francesco Haardt,
Elena Gallo,
Riccardo Spinelli,
Isaac Malsky,
Emily Rauscher
Abstract:
Intense X-ray and ultraviolet stellar irradiation can heat and inflate the atmospheres of closely orbiting exoplanets, driving mass outflows that may be significant enough to evaporate a sizable fraction of the planet atmosphere over the system lifetime. The recent surge in the number of known exoplanets, together with the imminent deployment of new ground and space-based facilities for exoplanet…
▽ More
Intense X-ray and ultraviolet stellar irradiation can heat and inflate the atmospheres of closely orbiting exoplanets, driving mass outflows that may be significant enough to evaporate a sizable fraction of the planet atmosphere over the system lifetime. The recent surge in the number of known exoplanets, together with the imminent deployment of new ground and space-based facilities for exoplanet discovery and characterization, requires a prompt and efficient assessment of the most promising targets for intensive spectroscopic follow-ups. To this purpose, we developed ATES (ATmospheric EScape); a new hydrodynamics code that is specifically designed to compute the temperature, density, velocity and ionization fraction profiles of highly irradiated planetary atmospheres, along with the current, steady-state mass loss rate. ATES solves the one-dimensional Euler, mass and energy conservation equations in radial coordinates through a finite-volume scheme. The hydrodynamics module is paired with a photoionization equilibrium solver that includes cooling via bremsstrahlung, recombination and collisional excitation/ionization for the case of a primordial atmosphere entirely composed of atomic hydrogen and helium, whilst also accounting for advection of the different ion species. Compared against the results of 14 moderately-to-highly irradiated planets simulated with The PLUTO-CLOUDY Interface (TPCI), ATES yields remarkably good agreement at a significantly smaller fraction of the computational time. A convergence study shows that ATES recovers stable, steady-state hydrodynamic solutions for systems with $\log(-φ_p) \lesssim 12.9 + 0.17\log F_{\rm XUV}$. Incidentally, atmospheres of systems above this threshold are generally thought to be undergoing Jeans escape. The code, which also features a user-friendly graphic interface, is available publicly as an online repository.
△ Less
Submitted 16 August, 2021; v1 submitted 18 June, 2021;
originally announced June 2021.
-
Non-detection of Helium in the Upper Atmospheres of Three Sub-Neptune Exoplanets
Authors:
David Kasper,
Jacob L. Bean,
Antonija Oklopčić,
Isaac Malsky,
Eliza M. -R. Kempton,
Jean-Michel Désert,
Leslie A. Rogers,
Megan Mansfield
Abstract:
We present a search for helium in the upper atmospheres of three sub-Neptune size planets to investigate the origins of these ubiquitous objects. The detection of helium for a low density planet would be strong evidence for the presence of a primary atmosphere accreted from the protoplanetary nebula because large amounts of helium are not expected in the secondary atmospheres of rocky planets. We…
▽ More
We present a search for helium in the upper atmospheres of three sub-Neptune size planets to investigate the origins of these ubiquitous objects. The detection of helium for a low density planet would be strong evidence for the presence of a primary atmosphere accreted from the protoplanetary nebula because large amounts of helium are not expected in the secondary atmospheres of rocky planets. We used Keck+NIRSPEC to obtain high-resolution transit spectroscopy of the planets GJ1214b, GJ9827d, and HD97658b around the 10,833 Ang He triplet feature. We did not detect helium absorption for any of the planets despite achieving a high level of sensitivity. We used the non-detections to set limits on the planets' thermosphere temperatures and atmospheric loss rates by comparing grids of 1D models to the data. We also performed coupled interior structure and atmospheric loss calculations, which suggest that the bulk atmospheres (winds) of the planets would be at most modestly enhanced (depleted) in helium relative to their primordial composition. Our lack of detections of the helium triplet for GJ1214b and GJ9827d are highly inconsistent with the predictions of models for the present day mass loss on these planets. Higher signal-to-noise data would be needed to detect the helium feature predicted for HD97658b. We identify uncertainties in the EUV fluxes of the host stars and the lack of detailed mass loss models specifically for cool and metal-enhanced atmospheres as the main limitations to the interpretation of our results. Ultimately, our results suggest that the upper atmospheres of sub-Neptune planets are fundamentally different than those of gas giant planets.
△ Less
Submitted 5 October, 2020; v1 submitted 25 July, 2020;
originally announced July 2020.
-
Coupled Thermal and Compositional Evolution of Photo Evaporating Planet Envelopes
Authors:
Isaac Malsky,
Leslie A. Rogers
Abstract:
Photo-evaporative mass loss sculpts the atmospheric evolution of tightly-orbiting sub-Neptune-mass exoplanets. To date, models of the mass loss from warm Neptunes have assumed that the atmospheric abundances remain constant throughout the planet's evolution. However, the cumulative effects of billions of years of escape modulated by diffusive separation and preferential loss of hydrogen can lead t…
▽ More
Photo-evaporative mass loss sculpts the atmospheric evolution of tightly-orbiting sub-Neptune-mass exoplanets. To date, models of the mass loss from warm Neptunes have assumed that the atmospheric abundances remain constant throughout the planet's evolution. However, the cumulative effects of billions of years of escape modulated by diffusive separation and preferential loss of hydrogen can lead to planetary envelopes that are enhanced in helium and metals relative to hydrogen (Hu et al. 2015). We have performed the first self-consistent calculations of the coupled thermal, mass-loss, and compositional evolution of hydrogen-helium envelopes surrounding sub-Neptune mass planets. We extended the MESA (Modules for Experiments in Stellar Astrophysics) stellar evolution code to model the evolving envelope abundances of photo-evaporating planets. We find that GJ 436b, the planet that originally inspired Hu et al. (2015) to propose the formation of helium enhanced planetary atmospheres, requires a primordial envelope that is too massive to become helium enhanced. Nonetheless, we show that helium enhancement is possible for planets with masses similar to GJ 436b after only several Gyr of mass loss. These planets have $R_p\lesssim 3.00~R_\oplus$, initial $f_{\rm env} < 0.5\%$, irradiation flux $\sim$10$^1$-10$^3$ times that of Earth, and obtain final helium fractions in excess of Y=0.40 in our models. The results of preferential envelope loss may have observable consequences on mass-radius relations and atmospheric spectra for sub-Neptune populations.
△ Less
Submitted 17 February, 2020; v1 submitted 15 February, 2020;
originally announced February 2020.
-
Signatures of Clouds in Hot Jupiter Atmospheres: Modeled High Resolution Emission Spectra from 3D General Circulation Models
Authors:
Caleb K. Harada,
Eliza M. -R. Kempton,
Emily Rauscher,
Michael Roman,
Isaac Malsky,
Marah Brinjikji,
Victoria diTomasso
Abstract:
Observations of scattered light and thermal emission from hot Jupiter exoplanets have suggested the presence of inhomogeneous aerosols in their atmospheres. 3D general circulation models (GCMs) that attempt to model the effects of aerosols have been developed to understand the physical processes that underlie their dynamical structures. In this work, we investigate how different approaches to aero…
▽ More
Observations of scattered light and thermal emission from hot Jupiter exoplanets have suggested the presence of inhomogeneous aerosols in their atmospheres. 3D general circulation models (GCMs) that attempt to model the effects of aerosols have been developed to understand the physical processes that underlie their dynamical structures. In this work, we investigate how different approaches to aerosol modeling in GCMs of hot Jupiters affect high-resolution thermal emission spectra throughout the duration of the planet's orbit. Using results from a GCM with temperature-dependent cloud formation, we calculate spectra of a representative hot Jupiter with different assumptions regarding the vertical extent and thickness of clouds. We then compare these spectra to models in which clouds are absent or simply post-processed (i.e., added subsequently to the completed clear model). We show that the temperature-dependent treatment of clouds in the GCM produces high-resolution emission spectra that are markedly different from the clear and post-processed cases -- both in the continuum flux levels and line profiles -- and that increasing the vertical extent and thickness of clouds leads to bigger changes in these features. We evaluate the net Doppler shifts of the spectra induced by global winds and the planet's rotation and show that they are strongly phase-dependent, especially for models with thicker and more extended clouds. This work further demonstrates the importance of radiative feedback in cloudy atmospheric models of hot Jupiters, as this can have a significant impact on interpreting spectroscopic observations of exoplanet atmospheres.
△ Less
Submitted 13 January, 2021; v1 submitted 4 December, 2019;
originally announced December 2019.