Search

Photochemistry is a fundamental process of planetary atmospheres that regulates the atmospheric composition and stability. However, no unambiguous photochemical products have been detected in exoplanet atmospheres to date. Recent observations from the JWST Transiting Exoplanet Early Release Science Program found a spectral absorption feature at 4.05 $\mu$m arising from SO$_2$ in the atmosphere of WASP-39b. WASP-39b is a 1.27-Jupiter-radii, Saturn-mass (0.28 M$_J$) gas giant exoplanet orbiting a Sun-like star with an equilibrium temperature of $\sim$1100 K. The most plausible way of generating SO$_2$ in such an atmosphere is through photochemical processes. Here we show that the SO$_2$ distribution computed by a suite of photochemical models robustly explains the 4.05 $\mu$m spectral feature identified by JWST transmission observations with NIRSpec PRISM (2.7$\sigma$) and G395H (4.5$\sigma$). SO$_2$ is produced by successive oxidation of sulphur radicals freed when hydrogen sulphide (H$_2$S) is destroyed. The sensitivity of the SO$_2$ feature to the enrichment of the atmosphere by heavy elements (metallicity) suggests that it can be used as a tracer of atmospheric properties, with WASP-39b exhibiting an inferred metallicity of $\sim$10$\times$ solar. We further point out that SO$_2$ also shows observable features at ultraviolet and thermal infrared wavelengths not available from the existing observations., Comment: 39 pages, 14 figures, accepted to be published in Nature

Measuring the metallicity and carbon-to-oxygen (C/O) ratio in exoplanet atmospheres is a fundamental step towards constraining the dominant chemical processes at work and, if in equilibrium, revealing planet formation histories. Transmission spectroscopy provides the necessary means by constraining the abundances of oxygen- and carbon-bearing species; however, this requires broad wavelength coverage, moderate spectral resolution, and high precision that, together, are not achievable with previous observatories. Now that JWST has commenced science operations, we are able to observe exoplanets at previously uncharted wavelengths and spectral resolutions. Here we report time-series observations of the transiting exoplanet WASP-39b using JWST's Near InfraRed Camera (NIRCam). The long-wavelength spectroscopic and short-wavelength photometric light curves span 2.0 - 4.0 $\mu$m, exhibit minimal systematics, and reveal well-defined molecular absorption features in the planet's spectrum. Specifically, we detect gaseous H$_2$O in the atmosphere and place an upper limit on the abundance of CH$_4$. The otherwise prominent CO$_2$ feature at 2.8 $\mu$m is largely masked by H$_2$O. The best-fit chemical equilibrium models favour an atmospheric metallicity of 1-100$\times$ solar (i.e., an enrichment of elements heavier than helium relative to the Sun) and a sub-stellar carbon-to-oxygen (C/O) ratio. The inferred high metallicity and low C/O ratio may indicate significant accretion of solid materials during planet formation or disequilibrium processes in the upper atmosphere., Comment: 35 pages, 13 figures, 3 tables, Nature, accepted

Young exoplanets are attractive targets for atmospheric characterization to explore the early phase of planetary evolution and the surrounding environment. Recent observations of the 10 Myr young Neptune-sized exoplanet K2-33b revealed that the planet's transit depth drastically decreases from the optical to near-infrared wavelengths. Thao et al. (2022) suggested that a thick planetary haze and/or stellar spots may be the cause; however, even the best-fit model only barely explains the data. Here, we propose that the peculiar transmission spectrum may indicate that K2-33b possess a circumplanetary dust ring, an analog of Jupiter's dust ring. We demonstrate that the ring could produce a steep slope in the transmission spectrum even if its optical depth is as low as $\sim {10}^{-2}$. We then apply a novel joint atmosphere-ring retrieval to K2-33b and find that the ring scenario could well explain the observed spectrum for various possible ring compositions. Importantly, the dust ring also exhibits prominent absorption features of ring particles around $\sim$10 $\rm{\mu}m$, whose shape and strength depend on the composition of the ring. Thus, future observations by JWST-MIRI would be able to test not only the ring hypothesis but also, if it indeed exists, to constrain the composition of the ring -- providing a unique opportunity to explore the origins of the dust ring around its parent planet, soon after the planetary system's formation., Comment: 8 pages, 4 figures, Accepted for publication in ApJL. Please also see Thao et al. (2022) for the transmission spectrum of K2-33b

We present high-resolution dayside thermal emission observations of the exoplanet WASP-18b using IGRINS on Gemini South. We remove stellar and telluric signatures using standard algorithms, and we extract the planet signal via cross correlation with model spectra. We detect the atmosphere of WASP-18b at a signal-to-noise ratio (SNR) of 5.9 using a full chemistry model, measure H2O (SNR=3.3), CO (SNR=4.0), and OH (SNR=4.8) individually, and confirm previous claims of a thermal inversion layer. The three species are confidently detected (>4$\sigma$) with a Bayesian inference framework, which we also use to retrieve abundance, temperature, and velocity information. For this ultra-hot Jupiter (UHJ), thermal dissociation processes likely play an important role. Retrieving abundances constant with altitude and allowing the temperature-pressure profile to freely adjust results in a moderately super-stellar carbon to oxygen ratio (C/O=0.75^{+0.14}_{-0.17}) and metallicity ([M/H]=1.03^{+0.65}_{-1.01}). Accounting for undetectable oxygen produced by thermal dissociation leads to C/O=0.45^{+0.08}_{-0.10} and [M/H]=1.17^{+0.66}_{-1.01}. A retrieval that assumes radiative-convective-thermochemical-equilibrium and naturally accounts for thermal dissociation constrains C/O<0.34 (2$\sigma$) and [M/H]=0.48^{+0.33}_{-0.29}, in line with the chemistry of the parent star. Looking at the velocity information, we see a tantalising signature of different Doppler shifts at the level of a few km/s for different molecules, which might probe dynamics as a function of altitude and location on the planet disk. Our results demonstrate that ground-based, high-resolution spectroscopy at infrared wavelengths can provide meaningful constraints on the compositions and climate of highly irradiated planets. This work also elucidates potential pitfalls with commonly employed retrieval assumptions when applied to UHJ spectra., Comment: 27 pages, 18 figures, submitted to AAS Journals. Community feedback welcome

We present the highest fidelity spectrum to date of a planetary-mass object. VHS 1256 b is a $<$20 M$_\mathrm{Jup}$ widely separated ($\sim$8\arcsec, a = 150 au), young, planetary-mass companion that shares photometric colors and spectroscopic features with the directly imaged exoplanets HR 8799 c, d, and e. As an L-to-T transition object, VHS 1256 b exists along the region of the color-magnitude diagram where substellar atmospheres transition from cloudy to clear. We observed VHS 1256~b with \textit{JWST}'s NIRSpec IFU and MIRI MRS modes for coverage from 1 $\mu$m to 20 $\mu$m at resolutions of $\sim$1,000 - 3,700. Water, methane, carbon monoxide, carbon dioxide, sodium, and potassium are observed in several portions of the \textit{JWST} spectrum based on comparisons from template brown dwarf spectra, molecular opacities, and atmospheric models. The spectral shape of VHS 1256 b is influenced by disequilibrium chemistry and clouds. We directly detect silicate clouds, the first such detection reported for a planetary-mass companion., Comment: Accepted ApJL. Iterations of spectra reduced by the ERS team are hosted at this link: https://github.com/bemiles/JWST_VHS1256b_Reduction/tree/main/reduced_spectra

We present JWST Early Release Science (ERS) coronagraphic observations of the super-Jupiter exoplanet, HIP 65426 b, with the Near-Infrared Camera (NIRCam) from 2-5 $\mu$m, and with the Mid-Infrared Instrument (MIRI) from 11-16 $\mu$m. At a separation of $\sim$0.82" (86$^{+116}_{-31}$ au), HIP 65426 b is clearly detected in all seven of our observational filters, representing the first images of an exoplanet to be obtained by JWST, and the first ever direct detection of an exoplanet beyond 5 $\mu$m. These observations demonstrate that JWST is exceeding its nominal predicted performance by up to a factor of 10, depending on separation and subtraction method, with measured 5$\sigma$ contrast limits of $\sim$1$\times10^{-5}$ and $\sim$2$\times10^{-4}$ at 1" for NIRCam at 4.4 $\mu$m and MIRI at 11.3 $\mu$m, respectively. These contrast limits provide sensitivity to sub-Jupiter companions with masses as low as 0.3$M_\mathrm{Jup}$ beyond separations of $\sim$100 au. Together with existing ground-based near-infrared data, the JWST photometry are well fit by a BT-SETTL atmospheric model from 1-16 $\mu$m, and span $\sim$97% of HIP 65426 b's luminous range. Independent of the choice of model atmosphere we measure an empirical bolometric luminosity that is tightly constrained between $\mathrm{log}\!\left(L_\mathrm{bol}/L_{\odot}\right)$=-4.31 to $-$4.14, which in turn provides a robust mass constraint of 7.1$\pm$1.2 $M_\mathrm{Jup}$. In totality, these observations confirm that JWST presents a powerful and exciting opportunity to characterise the population of exoplanets amenable to high-contrast imaging in greater detail., Comment: 35 pages, 16 figures, 4 tables, 1 wonderful telescope; Submitted to AAS Journals

Evidence of disequilibrium chemistry due to vertical mixing in the atmospheres of many T and Y-dwarfs has been inferred due to enhanced mixing ratios of CO and reduced NH$_3$. Atmospheric models of planets and brown dwarfs typically parameterize this vertical mixing phenomenon with the vertical eddy diffusion coefficient, $K_{\rm zz}$. While $K_{\rm zz}$ can perhaps be approximated in the convective regions in the atmosphere with mixing length theory, in radiative regions the strength of vertical mixing is uncertain by many orders of magnitude. With a new grid of self-consistent 1D model atmospheres from $T_{\rm eff}$ of 400 - 1000 K, computed with a new radiative-convective equilibrium python code PICASO 3.0, we aim to assess how molecular abundances and corresponding spectra can be used as a probe of depth-dependent $K_{\rm zz}$. At a given surface gravity, we find non-monotonic behavior in the CO abundance as a function of $T_{\rm eff}$, as chemical abundances are sometimes quenched in either of two potential atmospheric convective zones, or quenched in either of two possible radiative zones. The temperature structure and chemical quenching behavior also changes with gravity. We compare our models with available near-infrared and M-band spectroscopy of several T and Y-dwarfs and assess their atmospheric vertical mixing profiles. We also compare to color-magnitude diagrams and make predictions for JWST spectra. This work yields new constraints, and points the way to significant future gains, in determining $K_{\rm zz}$, a fundamental atmospheric parameter in substellar atmospheres, with significant implications for chemistry and cloud modeling., Comment: Accepted for Publication in The Astrophysical Journal, 31 Pages, 22 Figures. The model grid will be released via Zenodo

Carbon dioxide (CO2) is a key chemical species that is found in a wide range of planetary atmospheres. In the context of exoplanets, CO2 is an indicator of the metal enrichment (i.e., elements heavier than helium, also called "metallicity"), and thus formation processes of the primary atmospheres of hot gas giants. It is also one of the most promising species to detect in the secondary atmospheres of terrestrial exoplanets. Previous photometric measurements of transiting planets with the Spitzer Space Telescope have given hints of the presence of CO2 but have not yielded definitive detections due to the lack of unambiguous spectroscopic identification. Here we present the detection of CO2 in the atmosphere of the gas giant exoplanet WASP-39b from transmission spectroscopy observations obtained with JWST as part of the Early Release Science Program (ERS). The data used in this study span 3.0 to 5.5 {\mu}m in wavelength and show a prominent CO2 absorption feature at 4.3 {\mu}m (26{\sigma} significance). The overall spectrum is well matched by one-dimensional, 10x solar metallicity models that assume radiative-convective-thermochemical equilibrium and have moderate cloud opacity. These models predict that the atmosphere should have water, carbon monoxide, and hydrogen sulfide in addition to CO2, but little methane. Furthermore, we also tentatively detect a small absorption feature near 4.0 {\mu}m that is not reproduced by these models., Comment: 27 pages, 6 figures, Accepted for publication in Nature, data and models available at https://doi.10.5281/zenodo.6959427

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)., Comment: Under review in ApJ; posting after incorporating the first referee response; official code release will coincide with paper acceptance; comments welcome; 29 Pages; 15 Figures; 3 Tables

With the advent of direct imaging spectroscopy, the number of spectra from brown dwarfs and extrasolar gas giants is growing rapidly. Many brown dwarfs and extrasolar gas giants exhibit spectroscopic and photometric variability, which is likely the result of weather patterns. However, for the foreseeable future, point-source observations will be the only viable method to extract brown dwarf and exoplanet spectra. Models have been able to reproduce the observed variability, but ground truth observations are required to verify their results. To that end, we provide visual and near-infrared spectra of Jupiter and Saturn obtained from the \emph{Cassini} VIMS instrument. We disk-integrate the VIMS spectral cubes to simulate the spectra of Jupiter and Saturn as if they were directly imaged exoplanets or brown dwarfs. We present six empirical disk-integrated spectra for both Jupiter and Saturn with phase coverage of $1.7^\circ$ to $133.5^\circ$ and $39.6^\circ$ to $110.2^\circ$, respectively. To understand the constituents of these disk-integrated spectra, we also provide end member (single feature) spectra for permutations of illumination and cloud density, as well as for Saturn's rings. In tandem, these disk-integrated and end member spectra provide the ground truth needed to analyze point source spectra from extrasolar gas giants and brown dwarfs. Lastly, we discuss the impact that icy rings, such as Saturn's, have on disk-integrated spectra and consider the feasibility of inferring the presence of rings from direct imaging spectra.

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., Comment: Main text 16 pages, 11 figures

One of the strongest ${\rm Na~I}$ features was observed in WASP-96b. To confirm this novel detection, we provide a new 475-825nm transmission spectrum obtained with Magellan/IMACS, which indeed confirms the presence of a broad sodium absorption feature. We find the same result when reanalyzing the 400-825nm VLT/FORS2 data. We also utilize synthetic data to test the effectiveness of two common detrending techniques: (1) a Gaussian processes (GP) routine, and (2) common-mode correction followed by polynomial correction (CMC+Poly). We find that both methods poorly reproduce the absolute transit depths but maintain their true spectral shape. This emphasizes the importance of fitting for offsets when combining spectra from different sources or epochs. Additionally, we find that for our datasets both methods give consistent results, but CMC+Poly is more accurate and precise. We combine the Magellan/IMACS and VLT/FORS2 spectra with literature 800-1644nm HST/WFC3 spectra, yielding a global spectrum from 400-1644nm. We used the PLATON and Exoretrievals retrieval codes to interpret this spectrum, and find that both yield relatively deeper pressures where the atmosphere is optically thick at log-pressures between $1.3^{+1.0}_{-1.1}$ and 0.29$^{+1.86}_{-2.02}$ bars, respectively. Exoretrievals finds a solar to super-solar ${\rm Na~I}$ and ${\rm H_2O}$ log-mixing ratios of $-5.4^{+2.0}_{-1.9}$ and $-4.5^{+2.0}_{-2.0}$, respectively, while PLATON finds an overall metallicity of $log_{10}(Z/Z_{\odot}) = -0.49^{+1.0}_{-0.37}$dex. Therefore, our findings are in agreement with literature and support the inference that the terminator of WASP-96b has few aerosols obscuring prominent features in the optical to near-infrared (near-IR) spectrum., Comment: ACCEPT by AJ July 5th 2022

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., Comment: Accepted for publication in MNRAS

About 70-80% of stars in our solar and galactic neighborhood are M dwarfs. They span a range of low masses and temperatures relative to solar-type stars, facilitating molecule formation throughout their atmospheres. Standard stellar atmosphere models primarily designed for FGK stars face challenges when characterizing broadband molecular features in spectra of cool stars. Here, we introduce SPHINX--a new 1-D self-consistent radiative-convective thermochemical equilibrium chemistry model grid of atmospheres and spectra for M dwarfs in low-resolution (R~250). We incorporate the latest pre-computed absorption cross-sections with pressure-broadening for key molecules dominant in late-K, early/main-sequence-M stars. We then validate our grid models by acquiring fundamental properties (Teff, log(g), [M/H], radius, and C/O) for 10 benchmark M+G binary stars with known host metallicities and 10 M dwarfs with interferometrically measured angular diameters. Incorporating a Gaussian-process inference tool Starfish, we account for correlated and systematic noise in low-resolution (spectral stitching of SpeX, SNIFS, and STIS) observations and derive robust estimates of fundamental M dwarf atmospheric parameters. Additionally, we assess the influence of photospheric heterogeneity on acquired [M/H] and find that it could explain some deviations from observations. We also probe whether the model-assumed convective mixing-length parameter influences inferred radii, effective temperature, and [M/H] and again find that may explain discrepancies between interferometry observations and model-derived stellar parameters for cooler M dwarfs. Mainly, we show the unique strength in leveraging broadband molecular absorption features occurring in low-resolution M dwarf spectra and demonstrate the ability to improve constraints on fundamental properties of exoplanet hosts and late brown dwarf companions., Comment: Submitted to ApJ

We present Gemini South/IGRINS observations of the 1060 K T6 dwarf 2MASS J08173001$-$6155158 with unprecedented resolution ($R\equiv\lambda/\Delta\lambda=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 $\mu$m, and a molecular hydrogen (H$_2$) feature at 2.1218 $\mu$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., Comment: 45 pages, 63 figures, 5 tables, Accepted to MNRAS

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 $\Delta 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\mu$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., Comment: Accepted for publication in MNRAS

The direct characterization of exoplanetary systems with high contrast imaging is among the highest priorities for the broader exoplanet community. As large space missions will be necessary for detecting and characterizing exo-Earth twins, developing the techniques and technology for direct imaging of exoplanets is a driving focus for the community. For the first time, JWST will directly observe extrasolar planets at mid-infrared wavelengths beyond 5$\mu$m, deliver detailed spectroscopy revealing much more precise chemical abundances and atmospheric conditions, and provide sensitivity to analogs of our solar system ice-giant planets at wide orbital separations, an entirely new class of exoplanet. However, in order to maximise the scientific output over the lifetime of the mission, an exquisite understanding of the instrumental performance of JWST is needed as early in the mission as possible. In this paper, we describe our 55-hour Early Release Science Program that will utilize all four JWST instruments to extend the characterisation of planetary mass companions to $\sim$15$\mu$m as well as image a circumstellar disk in the mid-infrared with unprecedented sensitivity. Our program will also assess the performance of the observatory in the key modes expected to be commonly used for exoplanet direct imaging and spectroscopy, optimize data calibration and processing, and generate representative datasets that will enable a broad user base to effectively plan for general observing programs in future cycles., Comment: 21 pages, 7 figures. Accepted for Publication in PASP

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., Comment: 10 pages, 5 figures, 15 pages Supplementary Information, 3 Supplementary Figures

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., Comment: 31 pages, 16 figures, accepted for publication in PSJ

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$ $\mu$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\sigma$ 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., Comment: 12 pages, 7 figures, submitted to The Astronomical Journal

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., Comment: 16 pages, 5 figures, updated with minor formatting corrections

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., Comment: 10 pages, 8 figures, published to ApJ

We report observations of the recently discovered warm Neptune TOI-674 b (5.25 \rearth{}, 23.6 \mearth{}) with the Hubble Space Telescope's Wide Field Camera 3 instrument. TOI-674 b is in the Neptune desert, an observed paucity of Neptune-size exoplanets at short orbital periods. Planets in the desert are thought to have complex evolutionary histories due to photoevaporative mass loss or orbital migration, making identifying the constituents of their atmospheres critical to understanding their origins. We obtained near-infrared transmission spectroscopy of the planet's atmosphere with the G141 grism. After extracting, detrending, and fitting the spectral lightcurves to measure the planet's transmission spectrum, we used the petitRADTRANS atmospheric spectral synthesis code to perform retrievals on the planet's atmosphere to identify which absorbers are present. These results show moderate evidence for increased absorption at 1.4 $\mu$m due to water vapor at 2.9$\sigma$ (Bayes factor = 15.8), as well as weak evidence for the presence of clouds at 2.2$\sigma$ (Bayes factor = 4.0). TOI-674 b is a strong candidate for further study to refine the water abundance, which is poorly constrained by our data. We also incorporated new TESS short-cadence optical photometry, as well as Spitzer/IRAC data, and re-fit the transit parameters for the planet. We find the planet to have the following transit parameters: $R_p/R_* = 0.1135\pm0.0006$, $T_0 = 2458544.523792\pm0.000452$ BJD, and $P = 1.977198\pm0.00007$ d. These measurements refine the planet radius estimate and improve the orbital ephemerides for future transit spectroscopy observations of this highly intriguing warm Neptune., Comment: Accepted to the Astronomical Journal. 20 pages, 17 figures

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 $\mu$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., Comment: 16 pages, 7 figures, Accepted for publication in ApJ, Please also see a companion paper Alam et al. (2022) for the application of ring transmission spectrum to HIP 41378 f (arXiv:2201.02686)

Traditionally, ground-based spectrophotometric observations probing transiting exoplanet atmospheres have employed a linear map between comparison and target star light curves (e.g. via differential spectrophotometry) to correct for systematics contaminating the transit signal. As an alternative to this conventional method, we introduce a new Gaussian Processes (GP) regression-based method to analyse ground-based spectrophotometric data. Our new method allows for a generalised non-linear mapping between the target transit light curves and the time series used to detrend them. This represents an improvement compared to previous studies because the target and comparison star fluxes are affected by different telluric and instrumental systematics, which are complex and non-linear. We apply our method to six Gemini/GMOS transits of the warm (T$_{\rm eq}$ = 990 K) Neptune HAT-P-26b. We obtain on average $\sim$20 % better transit depth precision and residual scatter on the white light curve compared to the conventional method when using the comparison star light curve as a GP regressor and $\sim$20% worse when explicitly not using the comparison star. Ultimately, with only a cost of 30% precision on the transmission spectra, our method overcomes the necessity of using comparison stars in the instrument field of view, which has been one of the limiting factors for ground-based observations of the atmospheres of exoplanets transiting bright stars. We obtain a flat transmission spectrum for HAT-P-26b in the range of 490-900 nm that can be explained by the presence of a grey opacity cloud deck, and indications of transit timing variations, both of which are consistent with previous measurements., Comment: Accepted for publication in MNRAS

Relatively little is understood about the atmospheric composition of temperate to warm exoplanets (equilibrium temperature $T_{\rm eq}<$ 1000 K), as many of them are found to have uncharacteristically flat transmission spectra. Their flattened spectra are likely due to atmospheric opacity sources such as planet-wide photochemical hazes and condensation clouds. We compile the transmission spectra of 25 warm exoplanets previously observed by the Hubble Space Telescope and quantify the haziness of each exoplanet using a normalized amplitude of the water absorption feature ($A_{\rm H}$). By examining the relationships between $A_{\rm H}$ and various planetary and stellar forcing parameters, we endeavor to find correlations of haziness associated with planetary properties. We adopt new statistical correlation tests that are more suitable for the small, non-normally distributed warm exoplanet sample. Our analysis shows that none of the parameters hold statistically significant correlation with $A_{\rm H}$ ($p \le 0.01$) with the addition of new exoplanet data, including the previously identified linear trends between $A_{\rm H}$ and $T_{\rm{eq}}$ or hydrogen-helium envelope mass fraction (f$_{\rm{HHe}}$). This suggests that haziness in warm exoplanets is not simply controlled by any single planetary/stellar parameter. Among all the parameters we investigated, planet gravity ($g_{\rm p}$), atmospheric scale height ($H$), planet density ($\rho_{\rm p}$), orbital eccentricity ($e$), and age of the star ($t_{\rm age}$) hold tentative correlations with $A_{\rm H}$. Specifically, lower $H$, higher $g_{\rm p}$, $\rho_{\rm p}$, $e$, or $t_{\rm age}$ may lead to clearer atmospheres. We still need more observations and laboratory experiments to fully understand the complex physics and chemistry involved in creating hazy warm exoplanets., Comment: 29 pages, 6 figures, 6 tables, accepted by ApJ

We report \textit{Spitzer} full-orbit phase observations of the eccentric hot Jupiter XO-3b at 3.6 and 4.5 $\mu$m. Our new eclipse depth measurements of $1770 \pm 180$ ppm at 3.6 $\mu$m and $1610 \pm 70$ ppm at 4.5 $\mu$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'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 $\mu$m observations are pristine, the 3.6 $\mu$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 $\mu$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., Comment: 23 pages, 14 figures, accepted for publication in AJ

Venus's past climate evolution is uncertain. General circulation model simulations permit a habitable climate as late as ~0.7 Ga, and there is suggestive-albeit inconclusive-evidence for previous liquid water from surface geomorphology and mineralogy. However, it is unclear whether a habitable past can be reconciled with Venus's inferred atmospheric evolution. In particular, the lack of leftover atmospheric oxygen argues against recent water loss. Here, we apply a fully coupled model of Venus's atmospheric-interior-climate evolution from post-accretion magma ocean to present. The model self-consistently tracks C-, H-, and O-bearing volatiles and surface climate through the entirety of Venus's history. Atmospheric escape, mantle convection, melt production, outgassing, deep water cycling, and carbon cycling are explicitly coupled to climate and redox evolution. Plate tectonic and stagnant lid histories are considered. Using this coupled model, we conclude that both a habitable Venusian past and one where Venus never possessed liquid surface water can be reconciled with known constraints. Specifically, either scenario can reproduce bulk atmospheric composition, inferred surface heat flow, and observed $^{40}$Ar and $^{4}$He. Moreover, the model suggests that Venus could have been habitable with a ~100 m global ocean as late as 1 Ga, without violating any known constraints. In fact, if diffusion-limited water loss is throttled by a cool, CO$_2$-dominated upper atmosphere, then a habitable past is tentatively favored by our model. This escape throttling makes it difficult to simultaneously recover negligible water vapor and ~90 bar CO$_2$ in the modern atmosphere without temporarily sequestering carbon in the interior via silicate weathering to enhance H escape., Comment: Complete text 27 pages (main text 16 pages), including 9 figures and 7 supplementary figures

Measurements of the atmospheric carbon (C) and oxygen (O) relative to hydrogen (H) in hot Jupiters (relative to their host stars) provide insight into their formation location and subsequent orbital migration. Hot Jupiters that form beyond the major volatile (H2O/CO/CO2) ice lines and subsequently migrate post disk-dissipation are predicted have atmospheric carbon-to-oxygen ratios (C/O) near 1 and subsolar metallicities, whereas planets that migrate through the disk before dissipation are predicted to be heavily polluted by infalling O-rich icy planetesimals, resulting in C/O < 0.5 and super-solar metallicities. Previous observations of hot Jupiters have been able to provide bounded constraints on either H2O or CO, but not both for the same planet, leaving uncertain the true elemental C and O inventory and subsequent C/O and metallicity determinations. Here we report spectroscopic observations of a typical transiting hot Jupiter, WASP-77Ab. From these, we determine the atmospheric gas volume mixing ratio constraints on both H2O and CO (9.5$\times 10^{-5}$ - 1.5$\times 10^{-4}$ and 1.2$\times 10^{-4}$ - 2.6$\times 10^{-4}$, respectively). From these bounded constraints, we are able to derive the atmospheric C/H (0.35$^{+0.17}_{-0.10}$ $\times$ Solar) and O/H (0.32 $^{+0.12}_{-0.08}$ $\times$ Solar) abundances and the corresponding atmospheric carbon-to-oxygen ratio (C/O=0.59$\pm$0.08; the solar value is 0.55). The sub-solar (C+O)/H (0.33$^{+0.13}_{-0.09}$ $\times$ Solar) is suggestive of a metal-depleted atmosphere relative to what is expected for Jovian-like planets while the near solar value of C/O rules out the disk-free migration/C-rich atmosphere scenario., Comment: This is the accepted "pre-proof" version. Minor editorial reference/figure/abstract differences from published version. Published version here: https://www.nature.com/articles/s41586-021-03912-6

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., Comment: accepted to ApJ

The emergent spectra of close-in, giant exoplanets ("hot Jupiters") are expected to be distinct from those of self-luminous objects with similar effective temperatures because hot Jupiters are primarily heated from above by their host stars rather than internally from the release of energy from their formation. Theoretical models predict a continuum of dayside spectra for hot Jupiters as a function of irradiation level, with the coolest planets having absorption features in their spectra, intermediate-temperature planets having emission features due to thermal inversions, and the hottest planets having blackbody-like spectra due to molecular dissociation and continuum opacity from the H- ion. Absorption and emission features have been detected in the spectra of a number of individual hot Jupiters, and population-level trends have been observed in photometric measurements. However, there has been no unified, population-level study of the thermal emission spectra of hot Jupiters such as has been done for cooler brown dwarfs and transmission spectra of hot Jupiters. Here we show that hot Jupiter secondary eclipse spectra centered around a water absorption band at 1.4 microns follow a common trend in water feature strength with temperature. The observed trend is broadly consistent with model predictions for how the thermal structures of solar-composition planets vary with irradiation level. Nevertheless, the ensemble of planets exhibits some degree of scatter around the mean trend for solar composition planets. The spread can be accounted for if the planets have modest variations in metallicity and/or elemental abundance ratios, which is expected from planet formation models. (abridged abstract), Comment: 24 pages, 8 figures, published in Nature Astronomy

The detection of disk-integrated polarization from Luhman 16A and B in H-band, and subsequent modeling, has been interpreted in the framework of zonal cloud bands on these bodies. Recently, Tan and Showman (2021) investigated three-dimensional atmospheric circulation and cloud structures of brown dwarfs with general circulation models (GCMs), and their simulations yield complex cloud distributions showing some aspects of zonal jets, but also complex vortices that cannot be captured by a simple model. Here we use these 3D GCMs specific to Luhman 16A and B, along with the three-dimensional Monte Carlo radiative transfer code ARTES, to calculate their polarization signals. We adopt the 3D temperature-pressure and cloud profiles from the GCMs as our input atmospheric structures. Our polarization calculations at 1.6 $\mu$m agree well with the measured degree of linear polarization from both Luhman 16 A and B. Our calculations reproduce the measured polarization for both the objects with cloud particle sizes between 0.5-1 \,$\mu$m for Luhman 16 A and 5 \,$\mu$m for Luhman 16 B. We find that the degree of linear polarization can vary on hour-long timescales over the course of a rotation period. We also show that models with azimuthally symmetric band-like cloud geometries, typically used for interpreting polarimetry observations of brown dwarfs, over-predict the polarization signal if the cloud patterns do not include complex vortices within these bands. This exploratory work shows that GCMs are promising for modeling and interpreting polarization signals of brown dwarfs., Comment: 13 pages, 6 figures, Accepted for publication in The Astrophysical journal

Recent observations of the ultra-hot Jupiter WASP-76b have revealed a diversity of atmospheric species. Here we present new high-resolution transit spectroscopy of WASP-76b with GRACES at the Gemini North Observatory, serving as a baseline for the Large and Long Program "Exploring the Diversity of Exoplanet Atmospheres at High Spectral Resolution" (Exoplanets with Gemini Spectroscopy, or ExoGemS for short). With a broad spectral range of $400 - 1050$ nm, these observations allow us to search for a suite of atomic species. We recover absorption features due to neutral sodium (Na I), and report a new detection of the ionized calcium (Ca II) triplet at ~ $850$ nm in the atmosphere of WASP-76b, complementing a previous detection of the Ca II H & K lines. The triplet has line depths of $0.295 \pm 0.034$% at ~ $849.2$ nm, $0.574 \pm 0.041$% at ~ $854.2$ nm, and $0.454 \pm 0.024$% at ~ $866.2$ nm, corresponding to effective radii close to (but within) the planet's Roche radius. These measured line depths are significantly larger than those predicted by model LTE and NLTE spectra obtained on the basis of a pressure-temperature profile computed assuming radiative equilibrium. The discrepancy suggests that the layers probed by our observations are either significantly hotter than predicted by radiative equilibrium and/or in a hydrodynamic state. Our results shed light on the exotic atmosphere of this ultra-hot world, and will inform future analyses from the ExoGemS survey., Comment: 19 pages, 10 figures, 3 tables. Accepted for publication in Astrophysical Journal Letters

This brief review focuses on methods and applications of modeling exoplanetary atmospheres. We discuss various kinds of state of the art self-consistent and retrieval models in 1D and multi-D with a focus on open questions and short- and long-term goals in the field. Expertise previously developed in modeling cool stellar atmospheres and in modeling solar system planetary atmospheres has proven valuable to the field, and will continue to do so. We described upcoming opportunities for making progress in our understanding of atmospheres, and close with what we see as the field's challenges., Comment: Accepted in Autumn of 2020. To appear as a book chapter in "ExoFrontiers: Big questions in exoplanetary science", Ed. N Madhusudhan (Bristol: IOP Publishing Ltd) AAS-IOP ebooks https://iopscience.iop.org/bookListInfo/aas-iop-astronomy

We present a new generation of substellar atmosphere and evolution models, appropriate for application to studies of L, T, and Y-type brown dwarfs and self-luminous extrasolar planets. The atmosphere models describe the expected temperature-pressure profiles and emergent spectra of atmospheres in radiative-convective equilibrium with effective temperatures and gravities within the ranges $200\le T_{\rm eff}\le2400\,\rm K$ and $2.5\le \log g \le 5.5$. These ranges encompass masses from about 0.5 to 85 Jupiter masses for a set of metallicities ($[{\rm M/H}] = -0.5$ to $+0.5$), C/O ratios (from 0.5 to 1.5 times that of solar), and ages. The evolution tables describe the cooling of these substellar objects through time. These models expand the diversity of model atmospheres currently available, notably to cooler effective temperatures and greater ranges in C/O. Notable improvements from past such models include updated opacities and atmospheric chemistry. Here we describe our modeling approach and present our initial tranche of models for cloudless, chemical equilibrium atmospheres. We compare the modeled spectra, photometry, and evolution to various datasets., Comment: 27 pages, 16 figures, accepted for Astrophysical Journal. Models available at https://doi.org/10.5281/zenodo.5063476

Photochemical hazes are important opacity sources in temperate exoplanet atmospheres, hindering current observations from characterizing exoplanet atmospheric compositions. The haziness of an atmosphere is determined by the balance between haze production and removal. However, the material-dependent removal physics of the haze particles is currently unknown under exoplanetary conditions. Here we provide experimentally-measured surface energies for a grid of temperate exoplanet hazes to characterize haze removal in exoplanetary atmospheres. We found large variations of surface energies for hazes produced under different energy sources, atmospheric compositions, and temperatures. The surface energies of the hazes were found to be the lowest around 400 K for the cold plasma samples, leading to the lowest removal rates. We show a suggestive correlation between haze surface energy and atmospheric haziness with planetary equilibrium temperature. We hypothesize that habitable zone exoplanets could be less hazy, as they would possess high-surface-energy hazes which can be removed efficiently., Comment: 6 figures, 2 tables, 23 pages

The large radii of many hot Jupiters can only be matched by models that have hot interior adiabats, and recent theoretical work has shown that the interior evolution of hot Jupiters has a significant impact on their atmospheric structure. Due to its inflated radius, low gravity, and ultra-hot equilibrium temperature, WASP-76b is an ideal case study for the impact of internal evolution on observable properties. Hot interiors should most strongly affect the non-irradiated side of the planet, and thus full phase curve observations are critical to ascertain the effect of the interior on the atmospheres of hot Jupiters. In this work, we present the first Spitzer phase curve observations of WASP-76b. We find that WASP-76b has an ultra-hot day side and relatively cold nightside with brightness temperatures of $2471 \pm 27~\mathrm{K}$/$1518 \pm 61~\mathrm{K}$ at $3.6~\micron$ and $2699 \pm 32~\mathrm{K}$/$1259 \pm 44~\mathrm{K}$ at $4.5~\micron$, respectively. These results provide evidence for a dayside thermal inversion. Both channels exhibit small phase offsets of $0.68 \pm 0.48^{\circ}$ at $3.6~\micron$ and $0.67 \pm 0.2^{\circ}$ at $4.5~\mu\mathrm{m}$. We compare our observations to a suite of general circulation models that consider two end-members of interior temperature along with a broad range of frictional drag strengths. Strong frictional drag is necessary to match the small phase offsets and cold nightside temperatures observed. From our suite of cloud-free GCMs, we find that only cases with a cold interior can reproduce the cold nightsides and large phase curve amplitude at $4.5~\micron$, hinting that the hot interior adiabat of WASP-76b does not significantly impact its atmospheric dynamics or that clouds blanket its nightside., Comment: 24 pages, 10 Figures, 5 Tables. Accepted to AJ. Co-First Authors

We present ground-based optical transmission spectroscopy of the low-density hot Jupiter WASP-88b covering the wavelength range 4413-8333 {\AA} with the FORS2 spectrograph on the Very Large Telescope. The FORS2 white light curves exhibit a significant time-correlated noise which we model using a Gaussian Process and remove as a wavelength-independent component from the spectroscopic light curves. We analyse complementary photometric observations from the Transiting Exoplanet Survey Satellite and refine the system properties and ephemeris. We find a featureless transmission spectrum with increased absorption towards shorter wavelengths. We perform an atmospheric retrieval analysis with the AURA code, finding tentative evidence for haze in the upper atmospheric layers and a lower likelihood for a dense cloud deck. Whilst our retrieval analysis results point toward clouds and hazes, further evidence is needed to definitively reject a clear-sky scenario., Comment: 19 pages, 12 figures, 3 tables, Accepted for publication in MNRAS

Terrestrial planets with large water inventories are likely ubiquitous and will be among the first Earth-sized planets to be characterized with upcoming telescopes. It has previously been argued that waterworlds-particularly those possessing more than 1% H$_2$O-experience limited melt production and outgassing due to the immense pressure overburden of their overlying oceans, unless subject to high internal heating. But an additional, underappreciated obstacle to outgassing on waterworlds is the high solubility of volatiles in high-pressure melts. Here, we investigate this phenomenon and show that volatile solubilities in melts probably prevent almost all magmatic outgassing from waterworlds. Specifically, for Earth-like gravity and oceanic crust composition, oceans or water ice exceeding 10-100 km in depth (0.1-1 GPa) preclude the exsolution of volatiles from partial melt of silicates. This solubility limit compounds the pressure overburden effect as large surface oceans limit both melt production and degassing from any partial melt that is produced. We apply these calculations to Trappist-1 planets to show that, given current mass and radius constraints and implied surface water inventories, Trappist-1f and -1g are unlikely to experience volcanic degassing. While other mechanisms for interior-surface volatile exchange are not completely excluded, the suppression of magmatic outgassing simplifies the range of possible atmospheric evolution trajectories and has implications for interpretation of ostensible biosignature gases, which we illustrate with a coupled model of planetary interior-climate-atmosphere evolution., Comment: 16 pages, 11 figures

We present a ground-based optical transmission spectrum for the warm Saturn-mass exoplanet WASP-110b from two transit observations made with the FOcal Reducer and Spectrograph (FORS2) on the Very Large Telescope (VLT). The spectrum covers the wavelength range from 4000 to 8333\AA, which is binned in 46 transit depths measured to an averaged precision of 220 parts per million (ppm) over an averaged 80\AA~bin for a Vmag=12.8 star. The measured transit depths are unaffected by a dilution from a close A-type field dwarf, which was fully resolved. The overall main characteristic of the transmission spectrum is an increasing radius with wavelength and a lack of the theoretically predicted pressure-broadened sodium and potassium absorption features for a cloud-free atmosphere. We analyze archival high-resolution optical spectroscopy and find evidence for low to moderate activity of the host star, which we take into account in the atmospheric retrieval analysis. Using the AURA retrieval code, we find that the observed transmission spectrum can be best explained by a combination of unocculted stellar faculae and a cloud deck. Transmission spectra of cloud-free and hazy atmospheres are rejected at a high confidence. With a possible cloud deck at its terminator, WASP-110b joins the increasing population of irradiated hot-Jupiter exoplanets with cloudy atmospheres observed in transmission., Comment: Accepted for publication in The Astronomical Journal

Sub-Neptunes (Rp~1.25-4 REarth) remain the most commonly detected exoplanets to date. However, it remains difficult for observations to tell whether these intermediate-sized exoplanets have surfaces and where their surfaces are located. Here we propose that the abundances of trace species in the visible atmospheres of these sub-Neptunes can be used as proxies for determining the existence of surfaces and approximate surface conditions. As an example, we used a state-of-the-art photochemical model to simulate the atmospheric evolution of K2-18b and investigate its final steady-state composition with surfaces located at different pressures levels (Psurf). We find the surface location has a significant impact on the atmospheric abundances of trace species, making them deviate significantly from their thermochemical equilibrium and "no-surface" conditions. This result arises primarily because the pressure-temperature conditions at the surface determine whether photochemically-produced species can be recycled back to their favored thermochemical-equilibrium forms and transported back to the upper atmosphere. For an assumed H2-rich atmosphere for K2-18b, we identify seven chemical species that are most sensitive to the existence of surfaces: ammonia (NH3), methane (CH4), hydrogen cyanide (HCN), acetylene (C2H2), ethane (C2H6), carbon monoxide (CO), and carbon dioxide (CO2). The ratio between the observed and the no-surface abundances of these species, can help distinguish the existence of a shallow surface (Psurf < 10 bar), an intermediate surface (10 bar < Psurf < 100 bar), and a deep surface (Psurf > 100 bar). This framework can be applied together with future observations to other sub-Neptunes of interest., Comment: 32 pages, 13 figures, accepted in ApJ

Terrestrial exoplanets likely form initial atmospheres through outgassing during and after accretion, although there is currently no first-principles understanding of how to connect a planet's bulk composition to its early atmospheric properties. Important insights into this connection can be gained by assaying meteorites, representative samples of planetary building blocks. We perform laboratory outgassing experiments that use a mass spectrometer to measure the abundances of volatiles released when meteorite samples are heated to 1200 $^{\circ}$C. We find that outgassing from three carbonaceous chondrite samples consistently produce H$_2$O-rich (averaged ~66 %) atmospheres but with significant amounts of CO (~18 %) and CO$_2$ (~15 %) as well as smaller quantities of H$_2$ and H$_2$S (up to 1 %). These results provide experimental constraints on the initial chemical composition in theoretical models of terrestrial planet atmospheres, supplying abundances for principal gas species as a function of temperature., Comment: 42 pages, 9 Figures, 2 Tables, 4 Supplementary Figures, 3 Supplementary Tables, published in Nature Astronomy (2021)

Oxygen is a promising exoplanet biosignature due to the evolutionary advantage conferred by harnessing starlight for photosynthesis, and the apparent low likelihood of maintaining oxygen-rich atmospheres without life. Hypothetical scenarios have been proposed for non-biological oxygen accumulation on planets around late M-dwarfs, where the extended pre-main sequence may favor abiotic O$_2$ accumulation. In contrast, abiotic oxygen accumulation on planets around F, G, and K-type stars is seemingly less likely, provided they possess substantial non-condensable gas inventories. The comparative robustness of oxygen biosignatures around larger stars has motivated plans for next-generation telescopes capable of oxygen detection on planets around sun-like stars. However, the general tendency of terrestrial planets to develop oxygen-rich atmospheres across a broad range of initial conditions and evolutionary scenarios has not been explored. Here, we use a coupled thermal-geochemical-climate model of terrestrial planet evolution to illustrate three scenarios whereby significant abiotic oxygen can accumulate around sun-like stars, even when significant non-condensable gas inventories are present. For Earth-mass planets, we find abiotic oxygen can accumulate to modern levels if (1) the CO$_2$:H$_2$O ratio of the initial volatile inventory is high, (2) the initial water inventory exceeds ~50 Earth oceans, or (3) the initial water inventory is very low. Fortunately, these three abiotic oxygen scenarios could be distinguished from biological oxygen with observations of other atmospheric constituents or characterizing the planetary surface. This highlights the need for broadly capable next-generation telescopes that are equipped to constrain surface water inventories via time-resolved photometry and search for temporal biosignatures or disequilibrium biosignatures to assess whether oxygen is biogenic., Comment: Main text 20 pages, 7 figures. Supporting materials 45 pages

Optical, reflected light eclipse observations provide a direct probe of the exoplanet scattering properties, such as from aerosols. We present here the photometric, reflected light observations of WASP-43b using the HST WFC3/UVIS instrument with the F350LP filter (346-822nm) encompassing the entire optical band. This is the first reflected light, photometric eclipse using UVIS in scanning mode; as such we further detail our scanning extraction and analysis pipeline Arctor. Our HST WFC3/UVIS eclipse light curve for WASP-43 b derived a 3-{\sigma} upper limit of 67 ppm on the eclipse depth, which implies that WASP-43b has a very dark dayside atmosphere. With our atmospheric modeling campaign, we compared our reflected light constraints with predictions from global circulation and cloud models, benchmarked with HST and Spitzer observations of WASP-43b. We infer that we do not detect clouds on the dayside within the pressure levels probed by HST WFC3/UVIS with the F350LP filter (P > 1 bar). This is consistent with the GCM predictions based on previous WASP-43b observations. Dayside emission spectroscopy results from WASP-43b with HST and Spitzer observations are likely to not be significantly affected by contributions from cloud particles., Comment: 29 pages, 22 figures, accepted to AAS/ApJ

We provide a brief review of many aspects of the planetary physics of hot Jupiters. Our aim is to cover most of the major areas of current study while providing the reader with additional references for more detailed follow-up. We first discuss giant planet formation and subsequent orbital evolution via disk-driven torques or dynamical interactions. More than one formation pathway is needed to understand the population. Next, we examine our current understanding of the evolutionary history and current interior structure of the planets, where we focus on bulk composition as well as viable models to explain the inflated radii of the population. Finally we discuss aspects of their atmospheres in the context of observations and 1D and 3D models, including atmospheric structure and escape, spectroscopic signatures, and complex atmospheric circulation. The major opacity sources in these atmospheres, including alkali metals, water vapor, and others, are discussed. We discuss physics that control the 3D atmospheric circulation and day-to-night temperature structures. We conclude by suggesting important future work for still-open questions., Comment: References fixed, Figures revised, text edited slightly to better match final version, open access at: http://dx.doi.org/10.1029/2020JE006629

Constraining Jupiter's internal structure is crucial for understanding its formation and evolution history. Recent interior models of Jupiter that fit Juno's measured gravitational field suggest an inhomogeneous interior and potentially the existence of a diluted core. These models, however, strongly depend on the model assumptions and the equations of state used. A complementary modelling approach is to use empirical structure models. These can later be used to reveal new insights on the planetary interior and be compared to standard models. Here we present empirical structure models of Jupiter where the density profile is constructed by piecewise-polytropic equations. With these models we investigate the relation between the normalized moment of inertia (MoI) and the gravitational moments $J_2$ and $J_4$. Given that only the first few gravitational moments of Jupiter are measured with high precision, we show that an accurate and independent measurement of the MoI value could be used to further constrain Jupiter's interior. An independent measurement of the MoI with an accuracy better than $\sim 0.1\%$ could constrain Jupiter's core region and density discontinuities in its envelope. We find that models with a density discontinuity at $\sim$ 1 Mbar, as would produce a presumed hydrogen-helium separation, correspond to a fuzzy core in Jupiter. We next test the appropriateness of using polytropes, by comparing them with empirical models based on polynomials. We conclude that both representations result in similar density profiles and ranges of values for quantities like core mass and MoI., Comment: 16 pages, 12 figures, 5 tables, accepted for publication in ApJ

The unexpectedly large radii of hot Jupiters are a longstanding mystery whose solution will provide important insights into their interior physics. Many potential solutions have been suggested, which make diverse predictions about the details of inflation. In particular, although any valid model must allow for maintaining large planetary radii, only some allow for radii to increase with time. This reinflation process would potentially occur when the incident flux on the planet is increased. In this work, we examine the observed population of hot Jupiters to see if they grow as their parent stars brighten along the main sequence. We consider the relation between radius and other observables, including mass, incident flux, age, and fractional age (age over main sequence lifetime), and show that main sequence brightening is often sufficient to produce detectable reinflation. We further argue that these provide strong evidence for the relatively rapid reinflation of giant planets, and discuss the implications for proposed heating mechanisms. In our population analysis we also find evidence for a "delayed-cooling effect", wherein planets cool and contract far more slowly than expected. While not capable of explaining the observed radii alone, it may represent an important component of the effect. Finally, we identify a weak negative relationship between stellar metallicity and planet radius which is presumably the result of enhanced planetary bulk metallicity around metal-rich stars and has important implications for planet formation theory., Comment: 10 pages, 3 figures, 1 table. Revision improves clarity, but scientific results are unchanged. Accepted to The Astrophysical Journal Letters