Abstract
We present a forward-modeling framework using the Bayesian inference tool Starfish and cloudless Sonora-Bobcat model atmospheres to analyze low-resolution (
R
≈ 80–250) near-infrared ...(1.0–2.5
μ
m) spectra of T dwarfs. Our approach infers effective temperatures, surface gravities, metallicities, radii, and masses, and by accounting for uncertainties from model interpolation and correlated residuals due to instrumental effects and modeling systematics, produces more realistic parameter posteriors than traditional (
χ
2
-based) spectral-fitting analyses. We validate our framework by fitting the model atmospheres themselves and finding negligible offsets between derived and input parameters. We apply our methodology to three well-known benchmark late-T dwarfs, HD 3651B, GJ 570D, and Ross 458C, using both solar and nonsolar metallicity atmospheric models. We also derive these benchmarks’ physical properties using their bolometric luminosities, their primary stars’ ages and metallicities, and Sonora-Bobcat evolutionary models. Assuming the evolutionary-based parameters are more robust, we find our atmospheric-based, forward-modeling analysis produces two outcomes. For HD 3615B and GJ 570D, spectral fits provide accurate
T
eff
and
R
but underestimated
log
g
(by ≈1.2 dex) and
Z
(by ≈0.35 dex), likely due to the systematics from modeling the potassium line profiles. For Ross 458C, spectral fits provide accurate
log
g
and
Z
but overestimated
T
eff
(by ≈120 K) and underestimated
R
(by ≈1.6×), likely because our model atmospheres lack clouds, reduced vertical temperature gradients, or disequilibrium processes. Finally, the spectroscopically inferred masses of these benchmarks are all considerably underestimated.
The Kepler mission revealed a class of planets known as "super-puffs," with masses only a few times larger than Earth'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 reanalyzed 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 & Kreidberg, 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.
An important focus of exoplanet research is the determination of the atmospheric temperature structure of strongly irradiated gas giant planets, or hot Jupiters. HD 209458b is the prototypical ...exoplanet for atmospheric thermal inversions, but this assertion does not take into account recently obtained data or newer data reduction techniques. We reexamine this claim by investigating all publicly available Spitzer Space Telescope secondary-eclipse photometric data of HD 209458b and performing a self-consistent analysis. We employ data reduction techniques that minimize stellar centroid variations, apply sophisticated models to known Spitzer systematics, and account for time-correlated noise in the data. We derive new secondary-eclipse depths of 0.119% + or - 0.007%, 0.123% + or - 0.006%, 0.134% + or - 0.035%, and 0.215% + or - 0.008% in the 3.6, 4.5, 5.8, and 8.0 mum bandpasses, respectively. We feed these results into a Bayesian atmospheric retrieval analysis and determine that it is unnecessary to invoke a thermal inversion to explain our secondary-eclipse depths. The data are well fitted by a temperature model that decreases monotonically between pressure levels of 1 and 0.01 bars. We conclude that there is no evidence for a thermal inversion in the atmosphere of HD 209458b.
Phase curve observations provide an opportunity to study the energy budgets of exoplanets by quantifying the amount of heat redistributed from their daysides to their nightsides. Theories 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 ∼2000 K. In order to study the impact of H2 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 and a nightside temperature of . We observe a phase curve amplitude of 0.609 0.020 and an offset of . 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 H2 dissociation and recombination provide a better match to the data. The GCMs, however, predict a maximum phase offset of 5°, which disagrees with our observations at >5 confidence. This discrepancy may be due to magnetic effects in the planet's highly ionized atmosphere.
ABSTRACT
The exoplanet WASP-12b is the prototype for the emerging class of ultrahot, Jupiter-mass exoplanets. Past models have predicted – and near-ultraviolet observations have shown – that this ...planet is losing mass. We present an analysis of two sets of 3.6 and 4.5 $\mu \rm{m}$Spitzer phase curve observations of the system which show clear evidence of infrared radiation from gas stripped from the planet, and the gas appears to be flowing directly toward or away from the host star. This accretion signature is only seen at 4.5 $\mu \rm{m}$, not at 3.6 $\mu \rm{m}$, which is indicative either of CO emission at the longer wavelength or blackbody emission from cool, ≲600 K gas. It is unclear why WASP-12b is the only ultrahot Jupiter to exhibit this mass-loss signature, but perhaps WASP-12b’s orbit is decaying as some have claimed, while the orbits of other exoplanets may be more stable; alternatively, the high-energy irradiation from WASP-12A may be stronger than the other host stars. We also find evidence for phase offset variability at the level of 6.4σ (46.2°) at 3.6 $\mu \rm{m}$.
Exoplanet host star activity, in the form of unocculted starspots or faculae, alters the observed transmission and emission spectra of the exoplanet. This effect can be exacerbated when combining ...data from different epochs if the stellar photosphere varies between observations due to activity. Here, we present a method to characterize and correct for relative changes due to stellar activity by exploiting multi-epoch ( visits/transits) observations to place them in a consistent reference frame. Using measurements from portions of the planet's orbit where negligible planet transmission or emission can be assumed, we determine changes to the stellar spectral amplitude. With the analytical methods described here, we predict the impact of stellar variability on transit observations. Supplementing these forecasts with Kepler-measured stellar variabilities for F-, G-, K-, and M-dwarfs, and predicted transit precisions by the James Webb Space Telescope's (JWST) NIRISS, NIRCam, and MIRI, we conclude that stellar activity does not impact infrared transiting exoplanet observations of most presently known or predicted TESS targets by current or near-future platforms, such as JWST, as activity-induced spectral changes are below the measurement precision.
The Neptune-mass GJ 436b is one of the most studied transiting exoplanets with repeated measurements of its thermal emission and transmission spectra. We build on previous studies to answer ...outstanding questions about this planet, including its potentially high metallicity and tidal heating of its interior. We present new observations of GJ 436b's thermal emission at 3.6 and 4.5 m, which reduce uncertainties in estimates of GJ 436b's flux at those wavelengths and demonstrate consistency between Spitzer observations spanning more than 7 yr. We analyze the Spitzer thermal emission photometry and Hubble WFC3 transmission spectrum. We use a dual-pronged modeling approach of both self-consistent and retrieval models. We vary the metallicity, intrinsic luminosity from tidal heating, disequilibrium chemistry, and heat redistribution. We also study clouds and photochemical hazes, but do not find strong evidence for either. The self-consistent and retrieval models combine to suggest that GJ 436b has a high atmospheric metallicity, with best fits at or above several hundred times solar metallicity, tidal heating warming its interior with best-fit intrinsic effective temperatures around 300-350 K, and disequilibrium chemistry. High metal enrichments (>600× solar) occur from the accretion of rocky, rather than icy, material. Assuming the interior temperature Tint ∼ 300-350 K, we find a dissipation factor Q′ ∼ 2 × 105-106, larger than Neptune's Q′, implying a long tidal circularization timescale for the orbit. We suggest that Neptune-mass planets may be more diverse than imagined, with metal enhancements spanning several orders of magnitude, to perhaps over 1000× solar metallicity. High-fidelity observations with instruments like the James Webb Space Telescope will be critical for characterizing this diversity.
ABSTRACT
Ultra-hot Jupiters are perfect targets for transmission spectroscopy. However, their atmospheres feature strong spatial variations in temperature, chemistry, dynamics, cloud coverage, and ...scale height. This makes transit observations at high spectral resolution challenging to interpret. In this work, we model the cross-correlation signal of five chemical species – Fe, CO, H2O, OH, and TiO – on WASP-76b, a benchmark ultra-hot Jupiter. We compute phase-dependent high-resolution transmission spectra of three-dimensional (3D) SPARC/MITgcm models. The spectra are obtained with gCMCRT, a 3D Monte-Carlo radiative-transfer code. We find that, on top of atmospheric dynamics, the phase-dependent Doppler shift of the absorption lines in the planetary rest frame is shaped by the combined effect of planetary rotation and the unique 3D spatial distribution of chemical species. For species probing the dayside (e.g. refractories or molecules like CO and OH), the two effects act in tandem, leading to increasing blueshifts with orbital phase. For species that are depleted on the dayside (e.g. H2O and TiO), the two effects act in an opposite manner, and could lead to increasing redshifts during the transit. This behaviour yields species-dependent offsets from a planet’s expected Kp value that can be much larger than planetary wind speeds. The offsets are usually negative for refractory species. We provide an analytical formula to estimate the size of a planet’s Kp offsets, which can serve as a prior for atmospheric retrievals. We conclude that observing the phase-resolved absorption signal of multiple species is key to constraining the 3D thermochemical structure and dynamics of ultra-hot Jupiters.
ABSTRACT
We have developed an open-source pipeline for the analysis of Spitzer/IRAC channel 1 and 2 time-series photometry, incorporating some of the most popular decorrelation methods. We applied ...this pipeline to new phase curve observations of ultra-hot Jupiters MASCARA-1b and KELT-16b, and we performed the first comprehensive reanalysis of 15 phase curves. We find that MASCARA-1b and KELT-16b have phase offsets of $6^{+11}_{-11}~^{\circ }$W and $38^{+16}_{-15}~^{\circ }$W, dayside temperatures of $2952^{+100}_{-97}$ and $3070^{+160}_{-150}$ K, and nightside temperatures of $1300^{+340}_{-340}$ and $1900^{+430}_{-440}$ K, respectively. We confirm a strong correlation between dayside and irradiation temperatures with a shallower dependence for nightside temperature. We also find evidence that the normalized phase curve amplitude (peak-to-trough divided by eclipse depth) is correlated with stellar effective temperature. In addition, while our different models often retrieve similar parameters, significant differences occasionally arise between them, as well as between our preferred model and the literature values. Nevertheless, our preferred models are consistent with published phase offsets to within −8 ± 21 degrees (−1.6 ± 3.2 sigma), and normalized phase curve amplitudes are on average reproduced to within −0.01 ± 0.24 (−0.1 ± 1.6 sigma). Finally, we find that BLISS performs best in most cases, but not all; we therefore recommend future analyses consider numerous detector models to ensure an optimal fit and to assess model dependencies.
ABSTRACT Upcoming space-based coronagraphic instruments in the next decade will perform reflected light spectroscopy and photometry of cool directly imaged extrasolar giant planets. We are developing ...a new atmospheric retrieval methodology to help assess the science return and inform the instrument design for such future missions, and ultimately interpret the resulting observations. Our retrieval technique employs a geometric albedo model coupled with both a Markov chain Monte Carlo Ensemble Sampler (emcee) and a multimodal nested sampling algorithm (MultiNest) to map the posterior distribution. This combination makes the global evidence calculation more robust for any given model and highlights possible discrepancies in the likelihood maps. As a proof of concept, our current atmospheric model contains one or two cloud layers, methane as a major absorber, and a H2-He background gas. This 6-to-9 parameter model is appropriate for Jupiter-like planets and can be easily expanded in the future. In addition to deriving the marginal likelihood distribution and confidence intervals for the model parameters, we perform model selection to determine the significance of methane and cloud detection as a function of expected signal-to-noise ratio in the presence of spectral noise correlations. After internal validation, the method is applied to realistic spectra of Jupiter, Saturn, and HD 99492c, a model observing target. We find that the presence or absence of clouds and methane can be determined with high confidence, while parameter uncertainties are model dependent and correlated. Such general methods will also be applicable to the interpretation of direct imaging spectra of cloudy terrestrial planets.