Exoplanets and Stellar Physics

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:

We present a method for computing upper and lower limits to the expected duration of planetary transits given a range for the parameters of the host star, while explicitly accounting for non-zero impact parameter and eccentricity, and placing a basic constraint on the orbital stability through a minimum planet–star separation at periastron. We find that, especially at longer periods, the transit can be considerably shorter or longer than previous searches have assumed. No transits are known with such short or long transit durations, but it is unclear whether this is a real feature of the planet population or a combination of transit probability, observational bias, and detection bias.

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.

Abstract:

<jats:title>Abstract</jats:title> <jats:p>We have carried out a detailed analysis of the 3.4 μm spectral feature arising from Polycyclic Aromatic Hydrocarbons (PAH), using JWST archival data. For the first time in an external galaxy (NGC 6240), we have identified two distinct spectral components of the PAH 3.4 μm feature: a shorter wavelength component at 3.395 μm, which we attribute to short aliphatic chains tightly attached to the aromatic rings of the PAH molecules; and a longer wavelength feature at 3.405 μm that arises from longer, more fragile, aliphatic chains that are weakly attached to the parent PAH molecule. These longer chains are more easily destroyed by far-ultraviolet photons (&amp;gt;5eV) and PAH thermal emission only occurs where PAH molecules are shielded from more energetic photons by dense molecular gas. We see a very strong correlation in the morphology of the PAH 3.395 μm feature with the PAH 3.3 μm emission, the latter arising from robust aromatic PAH molecules. We also see an equally strong correlation between the PAH 3.405 μm morphology and the warm molecular gas, as traced by H2 vibrational lines. We show that the flux ratio PAH 3.395/PAH 3.405 &amp;lt; 0.3 corresponds strongly to regions where the PAH molecules are shielded by dense molecular gas, so that only modestly energetic UV photons penetrate to excite the PAHs. Our work shows that PAH 3.405 μm and PAH 3.395 μm emission features can provide robust diagnostics of the physical conditions of the interstellar medium in external galaxies, and can be used to quantify the energies of the photon field penetrating molecular clouds.</jats:p>