Exoplanet atmospheres

Abstract:

In Titan’s atmosphere, the chemistry of simple hydrocarbons (e.g., CH4 and C2H2) and nitrogen bearing species (e.g., N2 and CN) represents an important link between molecular species and the ubiquitous organic haze that gives Titan its characteristic orange hue. Here we present a new search for two previously undetected molecules, triacetylene (C6H2) and the gas phase dicyanoacetylene (C4N2), using the Echelon-Cross-Echelle Spectrograph instrument on board the Stratospheric Observatory for Infrared Astronomy aircraft. We do not detect these two molecules but determine upper limits for their mixing ratios and column abundances. We find the 3σ upper limits on the uniform volume mixing ratio (VMR) above 100 km for C6H2 to be 4.3 × 10−11, which is lower than the photochemical model predictions. This new upper limit suggests that the growth of linear molecules is inhibited. We also put a strict upper limit on the uniform VMR for gas phase C4N2 above 125 km to be 1.0 × 10−10. This upper limit is well below the saturation mixing ratio at this altitude for C4N2 and greatly limits the feasibility of C4N2 forming ice from condensation.

Abstract:

Atmospheric convection behaves differently in hydrogen-rich atmospheres compared to higher mean molecular weight atmospheres due to compositional gradients of tracers. Previous 1D studies predict that when a condensable tracer exceeds a critical mixing ratio in H2-rich atmospheres, convection is inhibited, leading to the formation of radiative layers where the temperature decreases faster with height than in convective profiles. We use 3D convection-resolving simulations to test whether convection is inhibited in H2-rich atmospheres when the tracer mixing ratio exceeds the critical threshold, while including processes neglected in 1D, e.g., turbulent mixing and evaporation. We run two sets of simulations. First, we perform simulations initialized on saturated isothermal states and find that compositional gradients can destabilize isothermal atmospheres. Second, we perform simulations initialized on adiabatic profiles, which show distinct, stable inhibition layers form when the condensable tracer exceeds the critical threshold. Within the inhibition layer, only a small amount of energy is carried by latent heat flux, and turbulent mixing transports a small amount of tracer upward, but both are generally too weak to sustain substantial tracer or heat transport. The thermal profile gradually relaxes to a steep radiative state, but radiative relaxation timescales are long. Our results suggest stable layers driven by condensation-induced convective inhibition form in H2-rich atmospheres, including those of sub-Neptune exoplanets.

Abstract:

The formation and evolution of giant planets remain incompletely understood, with mounting evidence that many close-in giants may have migrated from their birth locations. The detection of helium escaping the atmosphere of exoplanets has provided a powerful new tracer of atmospheric escape and exoplanetary evolution. Here, using high-precision spectroscopic observations from the James Webb Space Telescope (JWST) Near Infrared Imager and Slitless Spectrograph (NIRISS) in single-object slitless spectroscopy mode (SOSS) mode, we report the detection of substantial helium absorption during the pre-transit phase of WASP-107 b (17σ), as well as in the transit and post-transit phases. This unique continuous helium absorption begins approximately 1.5 h before the planet’s ingress and reveals the presence of an extended thermosphere. The observations show a maximum transit depth of 2.395 ± 0.01% near the helium triplet (36σ; at the NIRISS-SOSS resolution of ~700). Our ellipsoidal model of the planetary thermosphere matches the measured light curve well, suggesting an outflow extending to tens of planetary radii. Furthermore, we confidently detect water absorption (log10H2O = −2.5 ± 0.6), superimposed with a short-wavelength slope that we attribute to a prominent signature from unocculted stellar spots (5.2σ), rather than a small-particle haze slope. We place an upper limit on the abundance of K (log10K < −4.86, or K/H < 75× stellar) at 2σ, which is consistent with the O/H supersolar metallicity estimate. Together with the supersolar water abundance and the evidence for vigorous atmospheric escape, these findings suggest that WASP-107 b has undergone inward migration in its recent past, probably accompanied by strong tidal heating that continues to sustain its inflated atmosphere and mass loss. This investigation underscores the transformative potential of JWST for investigating planetary evolution.

Abstract:

Context. Characterizations of giant exoplanets such as β Pictoris b (hereafter β Pic b) are now routinely performed with multiple spectrographs and imagers exploring different spectral bandwidths and resolutions, allowing for atmospheric retrieval of spectra with or without the conservation of the planet spectral continuum. The accounting of data multimodality in the analysis could provide a more comprehensive determination of the planets physical and chemical properties and inform on their formation history. Aims. We present the first VLTI observations at R λ ∼4000 of β Pic b obtained for an exoplanet with GRAVITY at such a high resolution. We upgraded the forward modelling code ForMoSA to account for the data multimodality, including low-, medium-, and high-resolution spectroscopy based on both a direct model-data comparison and an analysis of cross-correlation signals. We used the ForMoSA code to refine the constraints on the atmospheric properties of the exoplanet and evaluated the sensitivity of the retrieved values to the input dataset. Methods. We obtained four high-signal-to-noise (S/N ∼ 20) spectra of β Pic b in the K band with GRAVITY at R λ ∼4000 conserving both the pseudo-continuum and the pattern of molecular absorptions. We used ForMoSA with four grids of self-consistent forward models (Exo-REM, ATMO, BT-Settl, and Sonora) to explore different T e ff , log(g), metallicity, C/O, and 12 CO/ 13 CO ratio values. We then combined the GRAVITY spectra with published 1–5 µm photometry (NaCo, VisAO, NICI, and SPHERE), low-to-mediumresolution ( R λ ≤ 700 broadband, 0.9–7 µm) spectra, and echelle spectra covering narrower bandwidths ( R λ ∼ 100 000, 2.1–5.2 µm). Results. Sonora and Exo-REM are statistically preferred among all four models, regardless of the dataset used. Exo-REM predicts T eff  = 1607.45 −6.20 +4.85 K and log(g) = 4.46 −0.04 +0.02 dex when using only the GRAVITY epochs, whereas we have T eff  = 1502.74 −2.14 +2.32 K log(g) = 4.00 ± 0.01 dex when incorporating all available datasets. The inclusion of archival data significantly affects all retrieved posteriors. When using all datasets, C/O mostly remains solar (0.552 −0.002 +0.003 ), while [M/H] reaches super-solar values (0.50 ± 0.01). We report the first tentative constraint on the isotopic ratio log( 12 CO/ 13 CO) = 1.12 −0.08 +0.11 in β Pic b’s atmosphere; however, we note that this detection remains inconclusive due to telluric residuals affecting both the GRAVITY and SINFONI data. Additionally, we estimated the bolometric luminosity as log(L/L ⊙ ) = −4.01 −0.05 +0.04 dex. Using a system age of 23 ± 3 Myr, along with this bolometric luminosity and the constraints on the dynamical mass of β Pic b, we were able to constrain the maximum of heavy element content of the planet to be on the order of 5% (20–80 M Earth ). Conclusions. The joint access to the pseudo-continuum and molecular lines in the K band provided by GRAVITY have a significant impact on the retrieved metallicity, possibly owing to the collision-induced absorption driving the continuum shape of the K band. The echelle spectra do not dominate the final fit with respect to lower resolution data covering a broader portion of the spectral energy distribution and the latter keeps encapsulating more robust information on T eff . Future multimodal frameworks should include a weighting scheme to account for the bandwidth and central wavelength of the observations.

Abstract:

Although TRAPPIST-1’s temperate planets have the highest transmission signals of any known system, flares contaminate 50%–70% of transits at the 1000 ppm level, far above 100 ppm secondary atmospheric signals. Efforts to mitigate flare contamination and assess impacts on radiation environments are each hampered by a lack of empirical spectral analysis and physics-based modeling. We present spectrotemporal analysis and radiative-hydrodynamic modeling of 5.5 hr of NIRISS and NIRSpec observations of six TRAPPIST-1 flares of 2.2–8.7 × 1030 erg. The flare lines and continua are characterized using grid searches of RADYN beam-heating models spanning 104 times in electron beam parameters. Best-fit models indicate these flares result from moderate-intensity beams with emergent electron fluxes of Fe = 1012 erg s−1 cm−2 and energies ≤37 keV, although all models overpredict the Paschen jump. These models predict X-ray and extreme UV (XUV), far-UV, and near-UV counterparts to the IR peak fluxes of 8.9–28.9 × 1027, 4.3–13.9 × 1026, and 3.4–11.4 × 1027 erg s−1, respectively. Scaling the flare rate into the XUV suggests flaring contributes 1.35 −0.15+2.0× quiescence yr−1. We bin integrations of similar flare effective temperature to construct fiducial flare spectra from 2000 to 4500 K, in order to develop separate empirical and RADYN-based mitigation pipelines. Both pipelines are applied to all 5.5 hr of R = 10 data, resulting in maximum residuals from 1 to 2.8 μm of 100–140 ppm and typical residuals of 54 ± 14 and 65 ± 17 ppm for the empirical and RADYN-based pipelines, respectively. Injection testing supports a 3σ detection capability for CO2 atmospheres with features of 150–250 ppm, with weak evidence (Bayes factor ≈ 3) still obtained at 130 ppm. Our results motivate multiwavelength observations to improve model fidelity and test high-energy predictions.