Abstract
Upcoming James Webb Space Telescope 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 1D and 3D 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 1Dl atmospheric radiative-convective equilibrium (RCE) model. This model has heritage from the Fortran-based code, 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 RCE 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).
Abstract Recent JWST eclipse spectra of the high-density hot Saturn HD 149026b between 2.35 and 5.08 μ m have allowed for in-depth study of its atmosphere. To understand its atmospheric properties, ...we have created a grid of 1D radiative-convective–thermochemical–equilibrium atmosphere models and spectra with PICASO 3.0 . In agreement with previous work, we find that the presence of gaseous TiO creates a thermal inversion, which is inconsistent with the data. The presence of gaseous VO, however, which condenses at temperatures 200 K cooler, does not cause such inversions but alters the temperature–pressure profile of the atmosphere. We estimate an atmospheric metallicity of 14 − 8 + 12 × solar without VO and 20 − 8 + 11 × solar with VO, a factor of ∼10 times smaller than previous work from Bean et al., who relied on atmosphere retrievals. We attribute this significant difference in metallicity to a larger temperature gradient at low pressures in radiative equilibrium models. Such models with lower metallicities readily fit the strong CO 2 feature at 4.3 μ m. Our lower estimated metallicity makes HD 149026b more consistent with the mass–metallicity relationship for other giant planets. We find a C/O ratio of 0.67 − 0.27 + 0.06 with and without VO. The best-fit heat redistribution factor without VO is 1.17, a very high value suggesting very little dayside energy transport and no energy transport to the nightside. The heat redistribution factor shrinks to a more plausible value of 0.91 − 0.05 + 0.05 , with VO, which we regard as circumstantial evidence for the molecule in the atmosphere of HD 149026b.
Abstract
We reanalyze the near-infrared spectra of the young extrasolar giant planet 51 Eridani b, which was originally presented in Macintosh et al. and Rajan et al. using modern atmospheric models, ...including a self-consistent treatment of disequilibrium chemistry due to turbulent vertical mixing. In addition, we investigate the possibility that significant opacity from micrometeors or other impactors in the planet’s atmosphere may be responsible for shaping the observed spectral energy distribution (SED). We find that disequilibrium chemistry is useful for describing the mid-infrared colors of the planet’s spectra, especially in regard to photometric data at the
M
band around 4.5
μ
m, which is the result of superequilibrium abundances of carbon monoxide, while the micrometeors are unlikely to play a pivotal role in shaping the SED. The best-fitting, micrometeoroid dust–free, disequilibrium chemistry, patchy cloud model has the following parameters: effective temperature
T
eff
= 681 K with clouds (or without clouds, i.e., the grid temperature
T
grid
= 900 K), surface gravity
g
= 1000 m s
−2
, sedimentation efficiency
f
sed
= 10, vertical eddy diffusion coefficient
K
zz
= 10
3
cm
2
s
−1
, cloud hole fraction
f
hole
= 0.2, and planet radius
R
planet
= 1.0
R
Jup
.
Abstract
We present the first JWST spectral energy distribution of a Y dwarf. This spectral energy distribution of the Y0 dwarf WISE J035934.06−540154.6 consists of low-resolution (
λ
/Δ
λ
∼100) ...spectroscopy from 1–12
μ
m and three photometric points at 15, 18, and 21
μ
m. The spectrum exhibits numerous fundamental, overtone, and combination rotational–vibrational bands of H
2
O, CH
4
, CO, CO
2
, and NH
3
, including the previously unidentified
ν
3
band of NH
3
at 3
μ
m. Using a Rayleigh–Jeans tail to account for the flux emerging at wavelengths greater than 21
μ
m, we measure a bolometric luminosity of 1.523 ± 0.090 × 10
20
W. We determine a semiempirical effective temperature estimate of
467
−
18
+
16
K using the bolometric luminosity and evolutionary models to estimate a radius. Finally, we compare the spectrum and photometry to a grid of atmospheric models and find reasonably good agreement with a model having
T
eff
= 450 K, log
g
= 3.25 cm s
−2
, and M/H = −0.3. However, the low surface gravity implies an extremely low mass of 1
M
Jup
and a very young age of 20 Myr, the latter of which is inconsistent with simulations of volume-limited samples of cool brown dwarfs.
Abstract Disequilibrium chemistry due to vertical mixing in the atmospheres of many brown dwarfs and giant exoplanets is well established. Atmosphere models for these objects typically parameterize ...mixing with the highly uncertain K zz diffusion parameter. The role of mixing in altering the abundances of C-N-O-bearing molecules has mostly been explored for atmospheres with a solar composition. However, atmospheric metallicity and the C/O ratio also impact atmospheric chemistry. Therefore, we present the Sonora Elf Owl grid of self-consistent cloud-free 1D radiative-convective equilibrium model atmospheres for JWST observations, which includes a variation in K zz across several orders of magnitude and also encompasses subsolar to supersolar metallicities and C/O ratios. We find that the impact of K zz on the T ( P ) profile and spectra is a strong function of both T eff and metallicity. For metal-poor objects, K zz has large impacts on the atmosphere at significantly higher T eff than in metal-rich atmospheres, where the impact of K zz is seen to occur at lower T eff . We identify significant spectral degeneracies between varying K zz and metallicity in multiple wavelength windows, in particular, at 3–5 μ m. We use the Sonora Elf Owl atmospheric grid to fit the observed spectra of a sample of nine early to late T-type objects from T eff = 550–1150 K. We find evidence for very inefficient vertical mixing in these objects, with inferred K zz values lying in the range between ∼10 1 and 10 4 cm 2 s −1 . Using self-consistent models, we find that this slow vertical mixing is due to the observations, which probe mixing in the deep detached radiative zone in these atmospheres.
Abstract WASP-69 b is a hot, inflated, Saturn-mass planet (0.26 M Jup ) with a zero-albedo equilibrium temperature of 963 K. Here, we report the JWST 2–12 μ m emission spectrum of the planet ...consisting of two eclipses observed with NIRCam grism time series and one eclipse observed with the MIRI low-resolution spectrometer (LRS). The emission spectrum shows absorption features of water vapor, carbon dioxide, and carbon monoxide, but no strong evidence for methane. WASP-69 b’s emission spectrum is poorly fit by cloud-free homogeneous models. We find three possible model scenarios for the planet: (1) a scattering model that raises the brightness at short wavelengths with a free geometric albedo parameter; (2) a cloud-layer model that includes high-altitude silicate aerosols to moderate long-wavelength emission; and (3) a two-region model that includes significant dayside inhomogeneity and cloud opacity with two different temperature–pressure profiles. In all cases, aerosols are needed to fit the spectrum of the planet. The scattering model requires an unexpectedly high geometric albedo of 0.64. Our atmospheric retrievals indicate inefficient redistribution of heat and an inhomogeneous dayside distribution, which is tentatively supported by MIRI LRS broadband eclipse maps that show a central concentration of brightness. Our more plausible models (2 and 3) retrieve chemical abundances enriched in heavy elements relative to solar composition by 6× to 14× solar and a C/O ratio of 0.65–0.94, whereas the less plausible highly reflective scenario (1) retrieves a slightly lower metallicity and lower C/O ratio.
Abstract
We observed HD 19467 B with JWST’s NIRCam in six filters spanning 2.5–4.6
μ
m with the long-wavelength bar coronagraph. The brown dwarf HD 19467 B was initially identified through a ...long-period trend in the radial velocity of the G3V star HD 19467. HD 19467 B was subsequently detected via coronagraphic imaging and spectroscopy, and characterized as a late-T type brown dwarf with an approximate temperature ∼1000 K. We observed HD 19467 B as a part of the NIRCam GTO science program, demonstrating the first use of the NIRCam Long Wavelength Bar coronagraphic mask. The object was detected in all six filters (contrast levels of 2 × 10
−4
to 2 × 10
−5
) at a separation of 1.″6 using angular differential imaging and synthetic reference differential imaging. Due to a guide star failure during the acquisition of a preselected reference star, no reference star data were available for post-processing. However, reference differential imaging was successfully applied using synthetic point-spread functions developed from contemporaneous maps of the telescope’s optical configuration. Additional radial velocity data (from Keck/HIRES) are used to constrain the orbit of HD 19467 B. Photometric data from TESS are used to constrain the properties of the host star, particularly its age. NIRCam photometry, spectra, and photometry from the literature, and improved stellar parameters are used in conjunction with recent spectral and evolutionary substellar models to derive the physical properties of HD 19467 B. Using an age of 9.4 ± 0.9 Gyr inferred from spectroscopy, Gaia astrometry, and TESS asteroseismology, we obtain a model-derived mass of 62 ± 1
M
J
, which is consistent within 2
σ
with the dynamically derived mass of
81
−
12
+
14
M
J
.
Abstract We report observations of the atmospheric transmission spectrum of the sub-Neptune exoplanet GJ 3470 b taken using the Near-Infrared Camera on JWST. Combined with two archival Hubble Space ...Telescope/Wide-Field Camera 3 transit observations and 15 archival Spitzer transit observations, we detect water, methane, sulfur dioxide, and carbon dioxide in the atmosphere of GJ 3470 b, each with a significance of >3 σ . GJ 3470 b is the lowest-mass—and coldest—exoplanet known to show a substantial sulfur dioxide feature in its spectrum, at M p = 11.2 M ⊕ and T eq = 600 K. This indicates that disequilibrium photochemistry drives sulfur dioxide production in exoplanet atmospheres over a wider range of masses and temperatures than has been reported or expected. The water, carbon dioxide, and sulfur dioxide abundances we measure indicate an atmospheric metallicity of approximately 100× solar. We see further evidence for disequilibrium chemistry in our inferred methane abundance, which is significantly lower than expected from equilibrium models consistent with our measured water and carbon dioxide abundances.
Abstract
Wolf 359 (CN Leo, GJ 406, Gaia DR3 3864972938605115520) is a low-mass star in the fifth-closest neighboring system (2.41 pc). Because of its relative youth and proximity, Wolf 359 offers a ...unique opportunity to study substellar companions around M stars using infrared high-contrast imaging and radial velocity monitoring. We present the results of
Ms
-band (4.67
μ
m) vector vortex coronagraphic imaging using Keck-NIRC2 and add 12 Keck-HIRES and 68 MAROON-X velocities to the radial velocity baseline. Our analysis incorporates these data alongside literature radial velocities from CARMENES, the High Accuracy Radial velocity Planet Searcher, and Keck-HIRES to rule out the existence of a close (
a
< 10 au) stellar or brown dwarf companion and the majority of large gas giant companions. Our survey does not refute or confirm the long-period radial velocity candidate, Wolf 359 b (
P
∼ 2900 days), but rules out the candidate's existence as a large gas giant (>4
M
Jup
) assuming an age of younger than 1 Gyr. We discuss the performance of our high-contrast imaging survey to aid future observers using Keck-NIRC2 in conjunction with the vortex coronagraph in the
Ms
band and conclude by exploring the direct imaging capabilities with JWST to observe Jupiter- and Neptune-mass planets around Wolf 359.