Search

2MASS 1207 b, the first directly imaged planetary-mass companion, has been instrumental in advancing our understanding of exoplanets and brown dwarfs over the past 20 years. We have performed extensive atmospheric retrieval analyses of 2MASS 1207 b's JWST/NIRSpec spectrum using petitRADTRANS and a new atmospheric inhomogeneity framework, which characterizes homogeneous atmospheres, patchy clouds, cloud-free hot spots, or the combination of patchy clouds and spots. Among 24 retrieval runs with various assumptions, the most statistically preferred model corresponds to the patchy cloud scheme, with Teff$=1174^{+4}_{-3}$ K, log(g)=$3.62^{+0.03}_{-0.02}$ dex, and R$=1.399^{+0.008}_{-0.010}$ R$_{\rm Jup}$, along with near-solar atmospheric compositions of [M/H]$=-0.05\pm0.03$ dex and C/O$=0.440\pm0.012$. This model suggests ~9% of 2MASS 1207 b's atmosphere is covered by thin iron clouds, producing L-dwarf-like spectra, while the remaining 91% consists of thick silicate and iron clouds, emitting blackbody-like spectra. These thin-cloud and thick-cloud regions resemble Jupiter's belts and zones, respectively; this scenario is consistently supported by other retrieval runs incorporating inhomogeneous atmospheres. We demonstrate that 2MASS 1207 b's weak CO absorption can be explained by the veiling effects of patchy thick clouds; the absence of 3.3um CH4 absorption is attributed to its hot thermal structure, which naturally leads to a CO-dominant, CH4-deficient atmosphere. The retrieved atmospheric models also match the observed variability amplitudes of 2MASS 1207 b. Our analysis reveals that the inferred atmospheric properties show significant scatter in less statistically preferred retrieval runs but converge to consistent values among the preferred ones. This underscores the importance of exploring diverse assumptions in retrievals to avoid biased interpretations of atmospheres and formation., Comment: 51 pages, 22 figures, and 8 tables. Submitted to AJ and fully revised following the referee report. This work is dedicated to the 20th anniversary of the discovery of 2MASS 1207 b, the first directly imaged planetary-mass companion. A review of 20 years of literature on the 2MASS 1207 A+b system is presented in Section 2

T-type brown dwarfs present an opportunity to explore atmospheres teeming with molecules such as H2O, CH4 and NH3, which exhibit a wealth of absorption features in the mid-infrared. With JWST, we can finally explore this chemistry in detail, including for the coldest brown dwarfs that were not yet discovered in the Spitzer era. This allows precise derivations of the molecular abundances, which in turn informs our understanding of vertical transport in these atmospheres and can provide clues about the formation of cold brown dwarfs and exoplanets. This study presents the first JWST/MRS mid-IR spectrum (R ~ 1500-3000) of a T-dwarf: the T8.5+T9 brown dwarf binary WISE J045853.90+643451.9. We fit the spectrum using a parameterized P-T profile and free molecular abundances (i.e., a retrieval analysis), treating the binary as unresolved. We find a good fit with a cloud-free atmosphere and identify H2O, CH4 and NH3 features. Moreover, we make the first detections of HCN and C2H2 (at 13.4$\sigma$ and 9.5$\sigma$ respectively) in any brown dwarf atmosphere. The detection of HCN suggests intense vertical mixing ($K_{zz}\sim10^{11}$cm$^2$s$^{-1}$), challenging previous literature derivations of $K_{zz}$ values for T-type brown dwarfs. Even more surprising is the C2H2 detection, which cannot be explained with existing atmospheric models for isolated objects. This result challenges model assumptions about vertical mixing, and/or our understanding of the C2H2 chemical network, or might hint towards a more complex atmospheric processes such as magnetic fields driving aurorae, or lightning driving ionization. These findings open a new frontier in studying carbon chemistry within brown dwarf atmospheres., Comment: Accepted for publication in The Astrophysical Journal Letters. 13 pages (5 figures) + appendices

The beta Pictoris system, with its two directly imaged planets beta Pic b and beta Pic c and its well characterised debris disk, is a prime target for detailed characterisation of young planetary systems. Here, we present high-resolution and high-contrast LM band spectroscopy with CRIRES+ of the system, primarily for the purpose of atmospheric characterisation of beta Pic b. We developed methods for determining slit geometry and wavelength calibration based on telluric absorption and emission lines, as well as methods for PSF modelling and subtraction, and artificial planet injection, in order to extract and characterise planet spectra at a high S/N and spectral fidelity. Through cross-correlation with model spectra, we detected H2O absorption for planet b in each of the 13 individual observations spanning four different spectral settings. This provides a clear confirmation of previously detected water absorption, and allowed us to derive an exquisite precision on the rotational velocity of beta Pic b, v_rot = 20.36 +/- 0.31 km/s, which is consistent within error bars with previous determinations. We also observed a tentative H2O cross-correlation peak at the expected position and velocity of planet c; the feature is however not at a statistically significant level. Despite a higher sensitivity to SiO than earlier studies, we do not confirm a tentative SiO feature previously reported for planet b. When combining data from different epochs and different observing modes for the strong H2O feature of planet b, we find that the S/N grows considerably faster when sets of different spectral settings are combined, compared to when multiple data sets of the same spectral setting are combined. This implies that maximising spectral coverage is often more important than maximising integration depth when investigating exoplanetary atmospheres using cross-correlation techniques., Comment: 18 pages, 21 figure, accepted for publication in A&A

Direct imaging observations are biased towards wide-separation, massive companions that have degenerate formation histories. Although the majority of exoplanets are expected to form via core accretion, most directly imaged exoplanets have not been convincingly demonstrated to follow this formation pathway. We obtained new interferometric observations of the directly imaged giant planet AF Lep b with the VLTI/GRAVITY instrument. We present three epochs of 50$\mu$as relative astrometry and the K-band spectrum of the planet for the first time at a resolution of R=500. Using only these measurements, spanning less than two months, and the Hipparcos-Gaia Catalogue of Accelerations, we are able to significantly constrain the planet's orbit; this bodes well for interferometric observations of planets discovered by Gaia DR4. Including all available measurements of the planet, we infer an effectively circular orbit ($e<0.02, 0.07, 0.13$ at $1, 2, 3 \sigma$) in spin-orbit alignment with the host, and a measure a dynamical mass of $M_\mathrm{p}=3.75\pm0.5\,M_\mathrm{Jup}$. Models of the spectrum of the planet show that it is metal rich ([M/H]$=0.75\pm0.25$), with a C/O ratio encompassing the solar value. This ensemble of results show that the planet is consistent with core accretion formation., Comment: Accepted to the Astronomical Journal. 12 figures, 4 tables

For modeling the spectra of exoplanets one must know their atmospheric composition. This is necessary because the abundance of molecules, atoms, ions and condensates is needed to construct the total cross-section for the interaction between electro-magnetic radiation and matter. In addition, when solving for the temperature structure of an atmosphere the so-called adiabatic temperature gradient must be known, which describes the pressure-temperature dependence in convectively unstable regions well. Depending on the planetary properties, the composition and adiabatic gradients may be well described by equilibrium chemistry, which means that chemical reactions occur faster than any other relevant processes in the atmosphere, such as mixing. What is more, the equilibrium assumption often serves as a useful starting point for non-equilibrium calculations. Efficient and easy-to-use codes for determining equilibrium abundances are therefore needed. Here we report on our easyCHEM Python package that calculates atmospheric compositions and adiabatic temperature gradients in chemical equilibrium for any user-specified elemental composition., Comment: Submitted for review to the Journal of Open Source Software (JOSS)

We present a transmission spectrum of the misaligned hot Jupiter WASP-15b from 2.8-5.2 microns observed with JWST's NIRSpec/G395H grating. Our high signal to noise data, which has negligible red noise, reveals significant absorption by H$_2$O ($4.2\sigma$) and CO$_2$ ($8.9\sigma$). From independent data reduction and atmospheric retrieval approaches, we infer that WASP-15b's atmospheric metallicity is super-solar ($\gtrsim 15\times$ solar) and its carbon-to-oxygen ratio is consistent with solar, that together imply planetesimal accretion. Our general circulation model simulations for WASP-15b suggest that the carbon-to-oxygen we measure at the limb is likely representative of the entire photosphere due to the mostly uniform spatial distribution of H$_2$O, CO$_2$ and CO. We additionally see evidence for absorption by SO$_2$ and absorption at 4.9$\mu$m, for which the current leading candidate is OCS, albeit with several caveats. If confirmed, this would be the first detection of OCS in an exoplanet atmosphere and point towards complex photochemistry of sulphur-bearing species in the upper atmosphere. These are the first observations from the BOWIE-ALIGN survey which is using JWST's NIRSpec/G395H instrument to compare the atmospheric compositions of aligned/low-obliquity and misaligned/high-obliquity hot Jupiters around F stars above the Kraft break. The goal of our survey is to determine whether the atmospheric composition differs across two populations of planets that have likely undergone different migration histories (disc versus disc-free) as evidenced by their obliquities (aligned versus misaligned)., Comment: 27 pages, 25 figures, 6 tables. Accepted to MNRAS on 30th January 2025

High-resolution spectroscopic characterization of young super-Jovian planets enables precise constraints on elemental and isotopic abundances of their atmospheres. As part of the ESO SupJup Survey, we present high-resolution spectral observations of two wide-orbit super-Jupiters in YSES 1 (or TYC 8998-760-1) using the upgraded VLT/CRIRES+ (R~100,000) in K-band. We carry out free atmospheric retrieval analyses to constrain chemical and isotopic abundances, temperature structures, rotation velocities, and radial velocities. We confirm the previous detection of 13CO in YSES 1 b at a higher significance of 12.6{\sigma}, but point to a higher 12CO/13CO ratio of 88+/-13 (1{\sigma} confidence interval), consistent with the primary's isotope ratio 66+/-5. We retrieve a solar-like composition in YSES 1 b with a C/O=0.57+/-0.01, indicating a formation via gravitational instability or core accretion beyond the CO iceline. Additionally, the observations lead to detections of H2O and CO in the outer planet YSES 1 c at 7.3{\sigma} and 5.7{\sigma}, respectively. We constrain the atmospheric C/O ratio of YSES 1 c to be either solar or subsolar (C/O=0.36+/-0.15), indicating the accretion of oxygen-rich solids. The two companions have distinct vsini, 5.34+/-0.14 km/s for YSES 1 b and 11.3+/-2.1 km/s for YSES 1 c, despite their similar natal environments. This may indicate different spin axis inclinations or effective magnetic braking by the long-lived circumplanetary disk around YSES 1 b. YSES 1 represents an intriguing system for comparative studies of super-Jovian companions and linking present atmospheres to formation histories., Comment: 25 pages, 11 figure, accepted for publication in AJ. The extracted CRIRES+ spectra of the YSES-1 system can be found at https://doi.org/10.5281/zenodo.13664032

We propose a comprehensive survey of giant planets ranging from close-in highly irradiated hot Jupiters to young, wide-orbit directly imaged planets. The combination of two established techniques for probing planetary atmospheric compositions (time-series transit observations and high-contrast spectroscopy) will provide an unprecedented window into gaseous planet compositions across a range of equilibrium temperatures (100-2000 K), orbital separations (0.1-100 au), and system ages (10 Myr-1 Gyr). To-date, compositional measurements of both transiting and directly imaged planets suggest two distinct formation pathways: star-like formation for directly imaged planets and planet-like formation for transiting planets. By leveraging the combined technical and theoretical expertise of the transiting and direct imaging communities, we can obtain a holistic view of giant planet formation, migration, accretion, and evolution. From the results of this comprehensive program, we will begin to answer one of the fundamental outstanding questions in our understanding of giant planets: where and how in the protoplanetary disk do giant planets form?, Comment: 4 pages, 2 figures. This white paper was submitted to the "Working Group on Strategic Exoplanet Initiatives with HST and JWST" call (https://sites.google.com/view/exoplanet-strategy-wg) and referenced in the final working group report (10.48550/arXiv.2404.02932)

We report results from 8 hours of JWST/MIRI LRS spectroscopic monitoring directly followed by 7 hours of JWST/NIRSpec prism spectroscopic monitoring of the benchmark binary brown dwarf WISE 1049AB, the closest, brightest brown dwarfs known. We find water, methane, and CO absorption features in both components, including the 3.3 $\mu$m methane absorption feature and a tentative detection of small grain ($<$ 1$\mu$m) silicate absorption at $>$8.5 $\mu$m in WISE 1049A. Both components vary significantly ($>$1$\%$), with WISE 1049B displaying larger variations than WISE 1049A. Using K-means clustering, we find three main transition points in wavelength for both components of the binary: 1) change in behavior at $\sim$2.3 $\mu$m coincident with a CO absorption bandhead, 2) change in behavior at 4.2 $\mu$m, close to the CO fundamental band at $\lambda >$ 4.4 $\mu$m, and 3) change in behavior at 8.3-8.5 $\mu$m, potentially corresponding to silicate absorption. We interpret the lightcurves observed with both NIRSpec and MIRI as likely stemming from 1) a deep pressure level driving the double-peaked variability seen in WISE 1049B at wavelengths $<$2.3 $\mu$m and $>$8.5 $\mu$m, 2) an intermediate pressure level shaping the lightcurve morphology between 2.3 and 4.2 $\mu$m, and 3) a higher-altitude pressure level producing single-peaked and plateaued lightcurve behavior between 4.2 and 8.5 $\mu$m., Comment: 28 pages, 27 figures, accepted to MNRAS

The metal mass fractions of gas giants are a powerful tool to constrain their formation mechanisms and evolution. The metal content is inferred by comparing mass and radius measurements with interior structure and evolution models. In the midst of the JWST, CHEOPS, TESS, and the forthcoming PLATO era, we are at the brink of obtaining unprecedented precision in radius, age and atmospheric metallicity measurements. To prepare for this wealth of data, we present the GAS gianT modeL for Interiors (GASTLI), an easy-to-use, publicly available Python package. The code is optimized to rapidly calculate mass-radius relations, and radius and luminosity thermal evolution curves for a variety of envelope compositions and core mass fractions. Its applicability spans planets with masses $17 \ M_{\oplus} < M < 6 \ M_{Jup}$, and equilibrium temperatures $T_{eq} < 1000$ K. The interior model is stratified in a core composed of water and rock, and an envelope constituted by H/He and metals (water). The interior is coupled to a grid of self-consistent, cloud-free atmospheric models to determine the atmospheric and boundary interior temperature, as well as the contribution of the atmosphere to the total radius. We successfully validate GASTLI by comparing it to previous work and data of the Solar System's gas giants and Neptune. We also test GASTLI on the Neptune-mass exoplanet HAT-P-26 b, finding a bulk metal mass fraction between 0.60-0.78 and a core mass of 8.5-14.4 $M_{\oplus}$. Finally, we explore the impact of different equations of state and assumptions, such as C/O ratio and transit pressure, in the estimation of bulk metal mass fraction. These differences between interior models entail a change in radius of up to 2.5% for Jupiter-mass planets, but more than 10\% for Neptune-mass. These are equivalent to variations in core mass fraction of 0.07, or 0.10 in envelope metal mass fraction., Comment: 18 pages, 9 figures. In review in Astronomy & Astrophysics

Using Keck Planet Imager and Characterizer (KPIC) high-resolution ($R$~35000) spectroscopy from 2.29-2.49 $\mu$m, we present uniform atmospheric retrievals for eight young substellar companions with masses of ~10-30 $M_\textrm{Jup}$, orbital separations spanning ~50-360 au, and $T_\textrm{eff}$ between ~1500-2600 K. We find that all companions have solar C/O ratios, and metallicities, to within the 1-2$\sigma$ level, with the measurements clustered around solar composition. Stars in the same stellar associations as our systems have near-solar abundances, so these results indicate that this population of companions is consistent with formation via direct gravitational collapse. Alternatively, core accretion outside the CO snowline would be compatible with our measurements, though the high mass ratios of most systems would require rapid core assembly and gas accretion in massive disks. On a population level, our findings can be contrasted with abundance measurements for directly imaged planets with m<10 $M_\textrm{Jup}$, which show tentative atmospheric metal enrichment. In addition, the atmospheric compositions of our sample of companions are distinct from those of hot Jupiters, which most likely form via core accretion. For two companions with $T_\textrm{eff}$~1700-2000 K (kap And b and GSC 6214-210 b), our best-fit models prefer a non-gray cloud model with >3$\sigma$ significance. The cloudy models yield 2-3$\sigma$ lower $T_\textrm{eff}$ for these companions, though the C/O and [C/H] still agree between cloudy and clear models at the $1\sigma$ level. Finally, we constrain 12CO/13CO for three companions with the highest S/N data (GQ Lup b, HIP 79098 b, and DH Tau b), and report $v$sin($i$) and radial velocities for all companions., Comment: Accepted to ApJ, 52 pages, 20 figures

Ground-based high-resolution spectra provide a powerful tool for characterising exoplanet atmospheres. However, they are greatly hampered by the dominating telluric and stellar lines, which need to be removed prior to any analysis. Such removal techniques ("preparing pipelines") deform the spectrum, hence a key point is to account for this process in the forward models used in retrievals. We develop a formal derivation on how to prepare froward models for retrievals, in the case where the telluric and instrumental deformations can be represented as a matrix multiplied element-wise with the data. We also introduce the notion of "Bias Pipeline Metric" (BPM), that can be used to compare the bias potential of preparing pipelines. We use the resulting framework to retrieve simulated observations of 1-D and 3-D exoplanet atmospheres and to re-analyse high-resolution ($\mathcal{R} \approx 80\,400$) near infrared (0.96--1.71 $\mu$m) CARMENES transit data of HD~189733~b. We compare these results with those obtained from a CCF analysis. With our fiducial retrieval, we find a blueshift of the absorption features of $-5.51^{+0.66}_{-0.53}$ km$\cdot$s$^{-1}$. In addition, we retrieve a H$_2$O $\log_{10}$(VMR) of $-2.39^{+0.12}_{-0.16}$ and a temperature of $660^{+6}_{-11}$ K. We are also able to put upper limits for the abundances of CH$_4$, CO, H$_2$S, HCN and NH$_3$, consistent with a sub-solar metallicity atmosphere enriched in H$_2$O. We retrieve a broadened line shape, consistent with rotation- and wind-induced line broadening. Finally, we find a lower limit for the pressure of an opaque cloud consistent with a clear atmosphere, and find no evidence for hazes., Comment: 40 pages, 25 figures, accepted in The Astronomical Journal

The unprecedented medium-resolution (R~1500-3500) near- and mid-infrared (1-18um) spectrum provided by JWST for the young (140+/-20Myr) low-mass (12-20MJup) L-T transition (L7) companion VHS1256b gives access to a catalogue of molecular absorptions. In this study, we present a comprehensive analysis of this dataset utilizing a forward modelling approach, applying our Bayesian framework, ForMoSA. We explore five distinct atmospheric models to assess their performance in estimating key atmospheric parameters: Teff, log(g), [M/H], C/O, gamma, fsed, and R. Our findings reveal that each parameter's estimate is significantly influenced by factors such as the wavelength range considered and the model chosen for the fit. This is attributed to systematic errors in the models and their challenges in accurately replicating the complex atmospheric structure of VHS1256b, notably the complexity of its clouds and dust distribution. To propagate the impact of these systematic uncertainties on our atmospheric property estimates, we introduce innovative fitting methodologies based on independent fits performed on different spectral windows. We finally derived a Teff consistent with the spectral type of the target, considering its young age, which is confirmed by our estimate of log(g). Despite the exceptional data quality, attaining robust estimates for chemical abundances [M/H] and C/O, often employed as indicators of formation history, remains challenging. Nevertheless, the pioneering case of JWST's data for VHS1256b has paved the way for future acquisitions of substellar spectra that will be systematically analyzed to directly compare the properties of these objects and correct the systematics in the models., Comment: 32 pages, 16 figures, 6 tables, 2 appendices

Context: High-resolution spectrographs fed by adaptive optics (AO) provide a unique opportunity to characterize directly imaged exoplanets. Observations with such instruments allow us to probe the atmospheric composition, spin rotation, and radial velocity of the planet, thereby helping to reveal information on its formation and migration history. The recent upgrade of the Cryogenic High-Resolution Infrared Echelle Spectrograph (CRIRES+) at the VLT makes it a highly suitable instrument for characterizing directly imaged exoplanets. Aims: In this work, we report on observations of $\beta$ Pictoris b with CRIRES+ and use them to constrain the planets atmospheric properties and update the estimation of its spin rotation. Methods: The data were reduced using the open-source \textit{pycrires} package. We subsequently forward-modeled the stellar, planetary, and systematic contribution to the data to detect molecules in the planet's atmosphere. We also used atmospheric retrievals to provide new constraints on its atmosphere. Results: We confidently detected water and carbon monoxide in the atmosphere of $\beta$ Pictoris b and retrieved a slightly sub-solar carbon-to-oxygen ratio, which is in agreement with previous results. The interpretation is hampered by our limited knowledge of the C/O ratio of the host star. We also obtained a much improved constraint on its spin rotation of $19.9 \pm 1.0$ km/s, which gives a rotation period of $8.7 \pm 0.8$ hours, assuming no obliquity. We find that there is a degeneracy between the metallicity and clouds, but this has minimal impact on the retrieved C/O, $v\sin{i}$, and radial velocity. Our results show that CRIRES+ is performing well and stands as a highly useful instrument for characterizing directly imaged planets., Comment: Accepted for publication in A&A

WASP-107b is a warm ($\sim$740 K) transiting planet with a Neptune-like mass of $\sim$30.5 $M_{\oplus}$ and Jupiter-like radius of $\sim$0.94 $R_{\rm J}$, whose extended atmosphere is eroding. Previous observations showed evidence for water vapour and a thick high-altitude condensate layer in WASP-107b's atmosphere. Recently, photochemically produced sulphur dioxide (SO$_2$) was detected in the atmosphere of a hot ($\sim$1,200 K) Saturn-mass planet from transmission spectroscopy near 4.05 $\mu$m, but for temperatures below $\sim$1,000 K sulphur is predicted to preferably form sulphur allotropes instead of SO$_2$. Here we report the 9$\sigma$-detection of two fundamental vibration bands of SO$_2$, at 7.35 $\mu$m and 8.69 $\mu$m, in the transmission spectrum of WASP-107b using the Mid-Infrared Instrument (MIRI) of the JWST. This discovery establishes WASP-107b as the second irradiated exoplanet with confirmed photochemistry, extending the temperature range of exoplanets exhibiting detected photochemistry from $\sim$1,200 K down to $\sim$740 K. Additionally, our spectral analysis reveals the presence of silicate clouds, which are strongly favoured ($\sim$7$\sigma$) over simpler cloud setups. Furthermore, water is detected ($\sim$12$\sigma$), but methane is not. These findings provide evidence of disequilibrium chemistry and indicate a dynamically active atmosphere with a super-solar metallicity.

Brown dwarfs serve as ideal laboratories for studying the atmospheres of giant exoplanets on wide orbits as the governing physical and chemical processes in them are nearly identical. Understanding the formation of gas giant planets is challenging, often involving the endeavour to link atmospheric abundance ratios, such as the carbon-to-oxygen (C/O) ratio, to formation scenarios. However, the complexity of planet formation requires additional tracers, as the unambiguous interpretation of the measured C/O ratio is fraught with complexity. Isotope ratios, such as deuterium-to-hydrogen and 14N/15N, offer a promising avenue to gain further insight into this formation process, mirroring their utility within the solar system. For exoplanets only a handful of constraints on 12C/13C exist, pointing to the accretion of 13C-rich ice from beyond the disks' CO iceline. Here we report on the mid-infrared detection of the 14NH3 and 15NH3 isotopologues in the atmosphere of a cool brown dwarf with an effective temperature of 380 K in a spectrum taken with the Mid-InfraRed Instrument of the James Webb Space Telescope. As expected, our results reveal a 14N/15N value consistent with star-like formation by gravitational collapse, demonstrating that this ratio can be accurately constrained. Since young stars and their planets should be more strongly enriched in the 15N isotope, we expect that 15NH3 will be detectable in a number of cold, wide-separation exoplanets., Comment: Accepted by Nature. 28 pages, 7 figures, uses nature3.cls

To understand high precision observations of exoplanets and brown dwarfs, we need detailed and complex general circulation models (GCMs) that incorporate hydrodynamics, chemistry, and radiation. For this study, we specifically examined the coupling between chemistry and radiation in GCMs and compared different methods for the mixing of opacities of different chemical species in the correlated-k assumption, when equilibrium chemistry cannot be assumed. We propose a fast machine learning method based on DeepSets (DS), which effectively combines individual correlated-k opacities (k-tables). We evaluated the DS method alongside other published methods such as adaptive equivalent extinction (AEE) and random overlap with rebinning and resorting (RORR). We integrated these mixing methods into our GCM (expeRT/MITgcm) and assessed their accuracy and performance for the example of the hot Jupiter HD~209458 b. Our findings indicate that the DS method is both accurate and efficient for GCM usage, whereas RORR is too slow. Additionally, we observed that the accuracy of AEE depends on its specific implementation and may introduce numerical issues in achieving radiative transfer solution convergence. We then applied the DS mixing method in a simplified chemical disequilibrium situation, where we modeled the rainout of TiO and VO, and confirmed that the rainout of TiO and VO would hinder the formation of a stratosphere. To further expedite the development of consistent disequilibrium chemistry calculations in GCMs, we provide documentation and code for coupling the DS mixing method with correlated-k radiative transfer solvers. The DS method has been extensively tested to be accurate enough for GCMs; however, other methods might be needed for accelerating atmospheric retrievals., Comment: Accepted for publication in A&A, language edited version

We present aperture masking interferometry (AMI) observations of the star HIP 65426 at $3.8\,\rm{\mu m}$ as a part of the JWST Direct Imaging Early Release Science (ERS) program obtained using the Near Infrared Imager and Slitless Spectrograph (NIRISS) instrument. This mode provides access to very small inner working angles (even separations slightly below the Michelson limit of $0.5\lambda/D$ for an interferometer), which are inaccessible with the classical inner working angles of the JWST coronagraphs. When combined with JWST's unprecedented infrared sensitivity, this mode has the potential to probe a new portion of parameter space across a wide array of astronomical observations. Using this mode, we are able to achieve a $5\sigma$ contrast of $\Delta m{\sim}7.62{\pm}0.13$ mag relative to the host star at separations ${\gtrsim}0.07{"}$, and the contrast deteriorates steeply at separations ${\lesssim}0.07{"}$. However, we detect no additional companions interior to the known companion HIP 65426 b (at separation ${\sim}0.82{"}$ or, $87^{+108}_{-31}\,\rm{au}$). Our observations thus rule out companions more massive than $10{-}12\,\rm{M_{Jup}}$ at separations ${\sim}10{-}20\,\rm{au}$ from HIP 65426, a region out of reach of ground or space-based coronagraphic imaging. These observations confirm that the AMI mode on JWST is sensitive to planetary mass companions at close-in separations (${\gtrsim}0.07{"}$), even for thousands of more distant stars at $\sim$100 pc, in addition to the stars in the nearby young moving groups as stated in previous works. This result will allow the planning and successful execution of future observations to probe the inner regions of nearby stellar systems, opening an essentially unexplored parameter space., Comment: 23 pages, 14 figures, accepted for publication in ApJL

We present a performance analysis for the aperture masking interferometry (AMI) mode on board the James Webb Space Telescope Near Infrared Imager and Slitless Spectrograph (JWST/NIRISS). Thanks to self-calibrating observables, AMI accesses inner working angles down to and even within the classical diffraction limit. The scientific potential of this mode has recently been demonstrated by the Early Release Science (ERS) 1386 program with a deep search for close-in companions in the HIP 65426 exoplanetary system. As part of ERS 1386, we use the same data set to explore the random, static, and calibration errors of NIRISS AMI observables. We compare the observed noise properties and achievable contrast to theoretical predictions. We explore possible sources of calibration errors and show that differences in charge migration between the observations of HIP 65426 and point-spread function calibration stars can account for the achieved contrast curves. Lastly, we use self-calibration tests to demonstrate that with adequate calibration NIRISS F380M AMI can reach contrast levels of $\sim9-10$ mag at $\gtrsim \lambda/D$. These tests lead us to observation planning recommendations and strongly motivate future studies aimed at producing sophisticated calibration strategies taking these systematic effects into account. This will unlock the unprecedented capabilities of JWST/NIRISS AMI, with sensitivity to significantly colder, lower-mass exoplanets than lower-contrast ground-based AMI setups, at orbital separations inaccessible to JWST coronagraphy., Comment: 20 pages, 12 figures, accepted to Astrophysical Journal Letters

petitRADTRANS (pRT) is a fast radiative transfer code used for computing emission and transmission spectra of exoplanet atmospheres, combining a FORTRAN back end with a Python based user interface. It is widely used in the exoplanet community with 222 references in the literature to date, and has been benchmarked against numerous similar tools. The spectra calculated with pRT can be used as a forward model for fitting spectroscopic data using Monte Carlo techniques, commonly referred to as an atmospheric retrieval. The new retrieval module combines fast forward modelling with nested sampling codes, allowing for atmospheric retrievals on a large range of different types of exoplanet data. Thus it is now possible to use pRT to easily and quickly infer the atmospheric properties of exoplanets in both transmission and thermal emission., Comment: 5 pages, 1 figure, published in the Journal of Open Source Software

AF Lep A+b is a remarkable planetary system hosting a gas-giant planet that has the lowest dynamical mass among directly imaged exoplanets. We present an in-depth analysis of the atmospheric composition of the star and planet to probe the planet's formation pathway. Based on new high-resolution spectroscopy of AF Lep A, we measure a uniform set of stellar parameters and elemental abundances (e.g., [Fe/H] = $-0.27 \pm 0.31$ dex). The planet's dynamical mass ($2.8^{+0.6}_{-0.5}$ M$_{\rm Jup}$) and orbit are also refined using published radial velocities, relative astrometry, and absolute astrometry. We use petitRADTRANS to perform chemically-consistent atmospheric retrievals for AF Lep b. The radiative-convective equilibrium temperature profiles are incorporated as parameterized priors on the planet's thermal structure, leading to a robust characterization for cloudy self-luminous atmospheres. This novel approach is enabled by constraining the temperature-pressure profiles via the temperature gradient $(d\ln{T}/d\ln{P})$, a departure from previous studies that solely modeled the temperature. Through multiple retrievals performed on different portions of the $0.9-4.2$ $\mu$m spectrophotometry, along with different priors on the planet's mass and radius, we infer that AF Lep b likely possesses a metal-enriched atmosphere ([Fe/H] $> 1.0$ dex). AF Lep b's potential metal enrichment may be due to planetesimal accretion, giant impacts, and/or core erosion. The first process coincides with the debris disk in the system, which could be dynamically excited by AF Lep b and lead to planetesimal bombardment. Our analysis also determines $T_{\rm eff} \approx 800$ K, $\log{(g)} \approx 3.7$ dex, and the presence of silicate clouds and dis-equilibrium chemistry in the atmosphere. Straddling the L/T transition, AF Lep b is thus far the coldest exoplanet with suggested evidence of silicate clouds., Comment: AJ, in press. Main text: Pages 1-32, Figures 1-15, Tables 1-6. All figures and tables after References belong to the Appendix (Pages 32-58, Figures 16-20, Table 7). For supplementary materials, please refer to the Zenodo repository https://doi.org/10.5281/zenodo.8267466

Many post-processing algorithms have been developed in order to better separate the signal of a companion from the bright light of the host star, but the effect of such algorithms on the shape of exoplanet spectra extracted from integral field spectrograph data is poorly understood. The resulting spectra are affected by noise that is correlated in wavelength space due to both optical and data processing effects. Within the framework of Bayesian atmospheric retrievals, we aim to understand how these correlations and other systematic effects impact the inferred physical parameters. We consider three algorithms (KLIP, PynPoint and ANDROMEDA), optimizing the choice of algorithmic parameters using a series of injection tests into archival SPHERE and GPI data of the HR 8799 system. The wavelength-dependent covariance matrix is calculated to provide a measure of instrumental and algorithmic systematics. We perform atmospheric retrievals using petitRADTRANS on optimally extracted spectra to measure how these data processing systematics influence the retrieved parameter distributions. The choice of data processing algorithm and parameters significantly impact the accuracy of retrieval results, with the mean posterior parameter bias ranging from 1 to 3 $\sigma$ from the true input parameters. Including the full covariance matrix in the likelihood improves the accuracy of inferred parameters, and cannot be accounted for using ad hoc scaling parameters in the retrieval framework. Using the Bayesian information criterion and other statistical measures as a heuristic goodness-of-fit metrics, the retrievals including the full covariance matrix are favoured when compared to using only the diagonal elements., Comment: 22 pages, 13 figures, accepted to Astronomy & Astrophysics

Retrieving the physical parameters from spectroscopic observations of exoplanets is key to understanding their atmospheric properties. Exoplanetary atmospheric retrievals are usually based on approximate Bayesian inference and rely on sampling-based approaches to compute parameter posterior distributions. Accurate or repeated retrievals, however, can result in very long computation times due to the sequential nature of sampling-based algorithms. We aim to amortize exoplanetary atmospheric retrieval using neural posterior estimation (NPE), a simulation-based inference algorithm based on variational inference and normalizing flows. In this way, we aim (i) to strongly reduce inference time, (ii) to scale inference to complex simulation models with many nuisance parameters or intractable likelihood functions, and (iii) to enable the statistical validation of the inference results. We evaluate NPE on a radiative transfer model for exoplanet spectra petitRADTRANS, including the effects of scattering and clouds. We train a neural autoregressive flow to quickly estimate posteriors and compare against retrievals computed with MultiNest. NPE produces accurate posterior approximations while reducing inference time down to a few seconds. We demonstrate the computational faithfulness of our posterior approximations using inference diagnostics including posterior predictive checks and coverage, taking advantage of the quasi-instantaneous inference time of NPE. Our analysis confirms the reliability of the approximate posteriors produced by NPE. The accuracy and reliability of the inference results produced by NPE establishes it as a promising approach for atmospheric retrievals. Amortization of the posterior inference makes repeated inference on several observations computationally inexpensive since it does not require on-the-fly simulations, making the retrieval efficient, scalable, and testable., Comment: The paper has been submitted to AandA after a final revision

Transmission spectroscopy is a prime technique to study the chemical composition and structure of exoplanetary atmospheres. Strong excess absorption signals have been detected in the optical Na I D1, 2 Fraunhofer lines during transits of hot Jupiters, which are attributed to the planetary atmospheres and allow us to constrain its structure. We study the atmosphere of WASP-7 b by means of high-resolution transit spectroscopy in the sodium lines. We analyzed a spectral transit time series of 89 high-resolution spectra of the hot Jupiter WASP-7 b that was observed using the Ultraviolet and Visual Echelle Spectrograph (UVES). We used the telluric lines for an accurate alignment of the spectra and carry out a telluric correction with molecfit. Stellar magnetic activity was monitored by investigating chromospheric lines such as the Ca II H and K and hydrogen H$\alpha$ lines. Finally, we obtained transmission spectra and light curves for various lines. The star shows no identifiable flares and, if any, marginal changes in activity during our observing run. The sodium transmission spectra and corresponding light curves clearly show signs of the Rossiter-McLaughlin effect (RM) and the stellar center-to-limb variation (CLV) that we modeled using synthetic spectra. A statistically significant, narrow absorption feature with a line contrast of 0.50$\pm$0.06% (at $\sim 8.3\sigma$ level) and a full width at half maximum (FWHM) of 0.13$\pm$0.02 A is detected at the location of the Na I D$_2$ line. For the Na I D$_1$ line signal, we derived a line contrast of 0.13$\pm$0.04% (at $\sim 3.2\sigma$ level), which we consider a tentative detection. In addition, we provide upper limits for absorption by the hydrogen Balmer lines (H$\alpha$, H$\beta$, and H$\gamma$), K I $\lambda$7699 A, Ca II H and K, and infra-red triplet (IRT) lines., Comment: Accepted for publication in A&A

Atmospheric retrieval studies are essential to determine the science requirements for future generation missions, such as the Large Interferometer for Exoplanets (LIFE). The use of heterogeneous absorption cross-sections might be the cause of systematic effects in retrievals, which could bias a correct characterization of the atmosphere. In this contribution we quantified the impact of differences in line list provenance, broadening coefficients, and line wing cut-offs in the retrieval of an Earth twin exoplanet orbiting a Sun-like star at 10 pc from the observer, as it would be observed with LIFE. We ran four different retrievals on the same input spectrum, by varying the opacity tables that the Bayesian retrieval framework was allowed to use. We found that the systematics introduced by the opacity tables could bias the correct estimation of the atmospheric pressure at the surface level, as well as an accurate retrieval of the abundance of some species in the atmosphere (such as CO$_2$ and N$_2$O). We argue that differences in the line wing cut-off might be the major source of errors. We highlight the need for more laboratory and modeling efforts, as well as inter-model comparisons of the main radiative transfer models and Bayesian retrieval frameworks. This is especially relevant in the context of LIFE and future generation missions, to identify issues and critical points for the community to jointly work together to prepare for the analysis of the upcoming observations., Comment: 24 pages, 12 figures. Proceedings SPIE Volume 12180, Space Telescopes and Instrumentation 2022: Optical, Infrared, and Millimeter Wave; 121803L (2022)

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

Direct imaging studies have mainly used low-resolution spectroscopy ($R\sim20-100$) to study the atmospheres of giant exoplanets and brown dwarf companions, but the presence of clouds has often led to degeneracies in the retrieved atmospheric abundances (e.g. C/O, metallicity). This precludes clear insights into the formation mechanisms of these companions. The Keck Planet Imager and Characterizer (KPIC) uses adaptive optics and single-mode fibers to transport light into NIRSPEC ($R\sim35,000$ in $K$ band), and aims to address these challenges with high-resolution spectroscopy. Using an atmospheric retrieval framework based on petitRADTRANS, we analyze KPIC high-resolution spectrum ($2.29-2.49~\mu$m) and archival low-resolution spectrum ($1-2.2~\mu$m) of the benchmark brown dwarf HD 4747 B ($m=67.2\pm1.8~M_{\rm{Jup}}$, $a=10.0\pm0.2$ au, $T_{\rm eff}\approx1400$ K). We find that our measured C/O and metallicity for the companion from the KPIC high-resolution spectrum agree with that of its host star within $1-2\sigma$. The retrieved parameters from the $K$ band high-resolution spectrum are also independent of our choice of cloud model. In contrast, the retrieved parameters from the low-resolution spectrum are highly sensitive to our chosen cloud model. Finally, we detect CO, H$_2$O, and CH$_4$ (volume mixing ratio of log(CH$_4$)=$-4.82\pm0.23$) in this L/T transition companion with the KPIC data. The relative molecular abundances allow us to constrain the degree of chemical disequilibrium in the atmosphere of HD 4747 B, and infer a vertical diffusion coefficient that is at the upper limit predicted from mixing length theory., Comment: 33 pages, 16 figures, Accepted to ApJ

We present simultaneous 0.65-2.5 microns medium resolution (3300 < R < 8100) VLT/X-Shooter spectra of the young low-mass (19+/-5MJup) L-T transition object VHS 1256-1257 b, a known spectroscopic analogue of HR8799d. The companion is a prime target for the JWST Early Release Science (ERS) and one of the highest-amplitude variable brown-dwarf known to date. We compare the spectrum to the custom grids of cloudless ATMO models exploring different atmospheric composition with the Bayesian inference tool ForMoSA. We also re-analyze low-resolution HST/WFC3 1.10-1.67 microns spectra at minimum and maximum variability to contextualize the X-Shooter data interpretation. The models reproduce the slope and most molecular absorption from 1.10 to 2.48 microns self-consistently but fail to provide a radius consistent with evolutionary model predictions. They do not reproduce consistently the optical spectrum and the depth of the K I doublets in the J-band. We derive Teff = 1380+/-54 K, log(g) = 3.97+/-0.48 dex, [M/H] = 0.21+/-0.29, and C/O > 0.63. Our inversion of the HST/WFC3 spectra suggests a relative change of 27+6-5 K of the disk-integrated Teff correlated with the near-infrared brightness. Our data anchor the characterization of that object in the near-infrared and could be used jointly to the ERS mid-infrared data to provide the most detailed characterization of an ultracool dwarf to date., Comment: 9 pages, 4 figures, 4 tables

Constraining planet formation based on the atmospheric composition of exoplanets is a fundamental goal of the exoplanet community. Existing studies commonly try to constrain atmospheric abundances, or to analyze what abundance patterns a given description of planet formation predicts. However, there is also a pressing need to develop methodologies that investigate how to transform atmospheric compositions into planetary formation inferences. In this study we summarize the complexities and uncertainties of state-of-the-art planet formation models and how they influence planetary atmospheric compositions. We introduce a methodology that explores the effect of different formation model assumptions when interpreting atmospheric compositions. We apply this framework to the directly imaged planet HR 8799e. Based on its atmospheric composition, this planet may have migrated significantly during its formation. We show that including the chemical evolution of the protoplanetary disk leads to a reduced need for migration. Moreover, we find that pebble accretion can reproduce the planet's composition, but some of our tested setups lead to too low atmospheric metallicities, even when considering that evaporating pebbles may enrich the disk gas. We conclude that the definitive inversion from atmospheric abundances to planet formation for a given planet may be challenging, but a qualitative understanding of the effects of different formation models is possible, opening up pathways for new investigations., Comment: Accepted for publication in ApJ, base formation inversion code is available at: https://gitlab.com/mauricemolli/formation-inversion

An important future goal in exoplanetology is to detect and characterize potentially habitable planets. Using nulling interferometry, LIFE will allow us to constrain the radius and effective temperature of (terrestrial) exoplanets, as well as provide unique information about their atmospheric structure and composition. We explore the potential of LIFE in characterizing emission spectra of Earth at various stages of its evolution. We perform Bayesian retrievals on simulated spectra of 8 different scenarios, which correspond to cloud-free and cloudy spectra of four different epochs of the evolution of the Earth. Assuming a distance of 10 pc and a Sun-like host star, we simulate observations obtained with LIFE using its simulator LIFEsim, considering all major astrophysical noise sources. With the nominal spectral resolution (R=50) and signal-to-noise ratio (assumed to be S/N=10 at 11.2 $\mu$m), we can identify the main spectral features of all the analyzed scenarios (most notably CO$_2$, H$_2$O, O$_3$, CH$_4$). This allows us to distinguish between inhabited and lifeless scenarios. Results suggest that particularly O$_3$ and CH$_4$ yield an improved abundance estimate by doubling the S/N from 10 to 20. We conclude that the baseline requirements for R and S/N are sufficient for LIFE to detect O$_3$ and CH$_4$ in the atmosphere of an Earth-like planet with an abundance of O$_2$ of around 2% in volume mixing ratio. This information is relevant in terms of the LIFE mission planning. We also conclude that cloud-free retrievals of cloudy planets can be used to characterize the atmospheric composition of terrestrial habitable planets, but not the thermal structure of the atmosphere. From the inter-model comparison performed, we deduce that differences in the opacity tables (caused by e.g. a different line wing treatment) may be an important source of systematic errors., Comment: 18 pages (main text, incl. 11 figures) + appendix; submitted to A&A; comments are very welcome! Fifth paper of LIFE telescope series. First: arXiv:2101.07500v4, Second: arXiv:2203.00471, Third: arXiv:2112.02054, Sixth: arXiv:2201.04891

Simulations with a 3D general circulation model (GCM) suggest that one potential driver behind the observed radius inflation in hot Jupiters may be the downward advection of energy from the highly irradiated photosphere into the deeper layers. Here, we compare dynamical heat transport within the non-inflated hot Jupiter WASP-43b and the canonical inflated hot Jupiter HD 209458b, with similar effective temperatures. We investigate to what extent the radiatively driven heating and cooling in the photosphere (at pressures smaller than 1 bar) influence the deeper temperature profile (at pressures between 1 to 700 bar). Our simulations with the new non-gray 3D radiation-hydrodynamical model expeRT/MITgcm show that the deep temperature profile of WASP-43b is associated with a relatively cold adiabat. The deep layers of HD 209458b, however, do not converge and remain nearly unchanged regardless of whether a cold or a hot initial state is used. Furthermore, we show that different flow structures in the deep atmospheric layers arise. There, we find that WASP-43b exhibits a deep equatorial jet, driven by the relatively fast tidally locked rotation of this planet (0.81 days), as compared to HD 209458b (3.47 days). However, by comparing simulations with different rotation periods, we find that the resulting flow structures only marginally influence the temperature evolution in the deep atmosphere, which is almost completely dominated by radiative heating and cooling. Furthermore, we find that the evolution of deeper layers can influence the 3D temperature structure in the photosphere of WASP-43b. Thus, dayside emission spectra of WASP-43b may shed more light onto the dynamical processes occurring at greater depths., Comment: accepted for publication in A&A

Ultra-hot Jupiters have dayside temperatures at which most molecules are expected to thermally dissociate. The dissociation of water vapour results in the production of the hydroxyl radical (OH). While OH absorption is easily observed in near-infrared spectra of M dwarfs, which have similar effective temperatures as ultra-hot Jupiters, it is often not considered when studying the atmospheres of ultra-hot Jupiters. We use high-resolution spectroscopic observations of a transit of WASP-76b obtained using CARMENES to study the presence of OH. After validating the OH line list, we generate model transit spectra of WASP-76b with petitRADTRANS. The data are corrected for stellar and telluric contamination and cross-correlated with the model spectra. After combining all cross-correlation functions from the transit, a detection map is constructed. Constraints on the planet properties from the OH absorption are obtained from a Markov chain Monte Carlo analysis. OH is detected in the atmosphere of WASP-76b with a peak signal-to-noise ratio of 6.1. From the retrieval we obtain $K_p=232 \pm 12$ km/s and a blueshift of $-13.2 \pm 1.6$ km/s, which are offset from the expected velocities. Considering the fast spin rotation of the planet, the blueshift is best explained with the signal predominantly originating from the evening terminator and the presence of a strong dayside-to-nightside wind. The increased $K_p$ over its expected value (196.5 km/s) is, however, a bit puzzling. The signal is found to be broad, with a full width at half maximum of $16.8^{+4.6}_{-4.0}$ km/s. The retrieval results in a weak constraint on the mean temperature of 2700-3700 K at the pressure range of the OH signal. We show that OH is readily observable in the transit spectra of ultra-hot Jupiters. Studying this molecule can provide insights into the molecular dissociation processes in the atmospheres of such planets., Comment: Accepted for publication in A&A

[abridged] Measurements of relative isotope abundances can provide unique insights into the formation and evolution histories of celestial bodies. The five stable isotopes of titanium are used to study the early history of the solar system and constrain Galactic chemical models. The minor isotopes of titanium are relatively abundant compared to those of other elements, making them more accessible for challenging observations. We assessed the feasibility of performing titanium isotope measurements in exoplanet atmospheres, and in particular, whether processing techniques used for high-resolution spectroscopy affect the derived isotope ratios. We used an archival high-dispersion CARMENES spectrum of the M-dwarf GJ 1002 as a proxy for an exoplanet observed at very high signal-to-noise. Spectral retrievals using petitRADTRANS models were performed on both narrow (7045-7090 {\AA}) and wide (7045-7500 {\AA}) wavelength regions, resulting in isotope ratios and uncertainties. These retrievals were repeated on the spectrum with its continuum removed to mimic typical high-dispersion exoplanet observations. The relative abundances of all minor Ti isotopes are found to be slightly enhanced compared to terrestrial values. Loss of continuum information from broadband spectral filtering has little effect on the isotope ratios. The CARMENES spectrum was subsequently degraded by adding varying levels of Gaussian noise to estimate the signal-to-noise requirements for future exoplanet atmospheric observations. For the wide wavelength range, a spectrum with signal-to-noise of 5 is required to determine the isotope ratios with relative errors $\lesssim$10%. Super Jupiters at large angular separations from their host star are the most accessible exoplanets, requiring about an hour of observing time on 8-meter-class telescopes, and less than a minute of observing time with the future Extremely Large Telescope., Comment: Accepted to A&A; 13 pages, 6 figures

Linking atmospheric characteristics of planets to their formation pathways is a central theme in the study of extrasolar planets. Although the 12C/13C isotope ratio shows little variation in the solar system, the atmosphere of a super-Jupiter was recently shown to be rich in 13CO, possibly as a result of dominant ice accretion beyond the CO snowline during its formation. In this paper we aim to measure the 12CO/13CO isotopologue ratio of a young, isolated brown dwarf. While the general atmospheric characteristics of young, low mass brown dwarfs are expected to be very similar to those of super-Jupiters, their formation pathways may be different, leading to distinct isotopologue ratios. We analyse archival K-band spectra of the L dwarf 2MASS J03552337+1133437 taken with NIRSPEC at the Keck telescope. A free retrieval analysis is applied to the data to determine the isotopologue ratio 12CO/13CO in its atmosphere. The isotopologue 13CO, is detected in the atmosphere through the cross-correlation method at a signal-to-noise of ~8.4. The detection significance is determined to be ~9.5 sigma using Bayesian model comparison between two retrieval models. We retrieve an isotopologue ratio 12CO/13CO of 97(+25/-18), marginally higher than the local interstellar standard. Its C/O ratio of ~0.56 is consistent with the solar value. Although only one super-Jupiter and one brown dwarf have now a measured 12CO/13CO ratio, it is intriguing that they are different, possibly hinting to distinct formation pathways. Albeit spectroscopic similarities, isolated brown dwarfs may experience a top-down formation via gravitational collapse, which resembles star formation, while giant exoplanets favorably form through core accretion, which potentially alters isotopologue ratios in their atmospheres depending on the material they accrete from protoplanetary disks. (abridged), Comment: Accepted to A&A; 9 pages, 4 figures

Isotope abundance ratios play an important role in astronomy and planetary sciences, providing insights in the origin and evolution of the Solar System, interstellar chemistry, and stellar nucleosynthesis. In contrast to deuterium/hydrogen ratios, carbon isotope ratios are found to be roughly constant (~89) in the Solar System, but do vary on galactic scales with 12C/13C~68 in the current local interstellar medium. In molecular clouds and protoplanetary disks, 12CO/13CO isotopologue ratios can be altered by ice and gas partitioning, low-temperature isotopic ion exchange reactions, and isotope-selective photodissociation. Here we report on the detection of 13CO in the atmosphere of the young, accreting giant planet TYC 8998-760-1 b at a statistical significance of >6 sigma. Marginalizing over the planet's atmospheric temperature structure, chemical composition, and spectral calibration uncertainties, suggests a 12CO/13CO ratio of 31 [+17,-10] (90% confidence), a significant enrichment in 13C with respect to the terrestrial standard and the local interstellar value. Since the current location of TYC 8998 b at >160 au is far beyond the CO snowline, we postulate that it accreted a significant fraction of its carbon from ices enriched in 13C through fractionation. Future isotopologue measurements in exoplanet atmospheres can provide unique constraints on where, when and how planets are formed., Comment: Published in Nature, July 14, 2021. A short video presentation can be found on exoplanet-talks.org at https://exoplanet-talks.org/talk/377

We report the discovery of a transiting, temperate, Neptune-sized exoplanet orbiting the nearby ($d$ = 27.5 pc), M3V star TOI-1231 (NLTT 24399, L 248-27, 2MASS J10265947-5228099). The planet was detected using photometric data from the Transiting Exoplanet Survey Satellite and followed up with observations from the Las Cumbres Observatory and the Antarctica Search for Transiting ExoPlanets program. Combining the photometric data sets, we find that the newly discovered planet has a radius of 3.65$^{+0.16}_{-0.15}$ R$_{\oplus}$, and an orbital period of 24.246 days. Radial velocity measurements obtained with the Planet Finder Spectrograph on the Magellan Clay telescope confirm the existence of the planet and lead to a mass measurement of 15.5$\pm$3.3 M$_{\oplus}$. With an equilibrium temperature of just 330K TOI-1231 b is one of the coolest small planets accessible for atmospheric studies thus far, and its host star's bright NIR brightness (J=8.88, K$_{s}$=8.07) make it an exciting target for HST and JWST. Future atmospheric observations would enable the first comparative planetology efforts in the 250-350 K temperature regime via comparisons with K2-18 b. Furthermore, TOI-1231's high systemic radial velocity (70.5 k\ms) may allow for the detection of low-velocity hydrogen atoms escaping the planet by Doppler shifting the H I Ly-alpha stellar emission away from the geocoronal and ISM absorption features., Comment: 20 pages, 9 figures. Accepted for publication in The Astronomical Journal

The atmospheres of synchronously rotating exoplanets are intrinsically three-dimensional, and fast vertical and horizontal winds are expected to mix the atmosphere, driving the chemical composition out of equilibrium. Due to the longer computation times associated with multi-dimensional forward models, horizontal mixing has only been investigated for a few case studies. In this paper, we aim to generalize the impact of horizontal and vertical mixing on the chemistry of exoplanet atmospheres over a large parameter space. We do this by applying a sequence of post-processed forward models for a large grid of synchronously rotating gaseous exoplanets, where we vary the effective temperature (between 400 K and 2600 K), surface gravity, and rotation rate. We find that there is a dichotomy in the horizontal homogeneity of the chemical abundances. Planets with effective temperatures below 1400 K tend to have horizontally homogeneous, vertically quenched chemical compositions, while planets hotter than 1400 K exhibit large compositional day-night differences for molecules such as methane. Furthermore, we find that the planet's rotation rate impacts the planetary climate, and thus also the molecular abundances and transmission spectrum. By employing a hierarchical modelling approach, we assess the relative importance of disequilibrium chemistry on the exoplanet transmission spectrum, and conclude that the temperature has the most profound impact. Temperature differences are also the main cause of limb asymmetries, which we estimate could be observable with the James Webb Space Telescope. This work highlights the value of applying a consistent modelling setup to a broad parameter space in exploratory theoretical research., Comment: 51 pages, 37 figures, accepted for publication in MNRAS. Erratum: Figure 17 corrected

Ariel, the Atmospheric Remote-sensing Infrared Exoplanet Large-survey, was adopted as the fourth medium-class mission in ESA's Cosmic Vision programme to be launched in 2029. During its 4-year mission, Ariel will study what exoplanets are made of, how they formed and how they evolve, by surveying a diverse sample of about 1000 extrasolar planets, simultaneously in visible and infrared wavelengths. It is the first mission dedicated to measuring the chemical composition and thermal structures of hundreds of transiting exoplanets, enabling planetary science far beyond the boundaries of the Solar System. The payload consists of an off-axis Cassegrain telescope (primary mirror 1100 mm x 730 mm ellipse) and two separate instruments (FGS and AIRS) covering simultaneously 0.5-7.8 micron spectral range. The satellite is best placed into an L2 orbit to maximise the thermal stability and the field of regard. The payload module is passively cooled via a series of V-Groove radiators; the detectors for the AIRS are the only items that require active cooling via an active Ne JT cooler. The Ariel payload is developed by a consortium of more than 50 institutes from 16 ESA countries, which include the UK, France, Italy, Belgium, Poland, Spain, Austria, Denmark, Ireland, Portugal, Czech Republic, Hungary, the Netherlands, Sweden, Norway, Estonia, and a NASA contribution., Comment: Ariel Definition Study Report, 147 pages. Reviewed by ESA Science Advisory Structure in November 2020. Original document available at: https://www.cosmos.esa.int/documents/1783156/3267291/Ariel_RedBook_Nov2020.pdf/

[abridged] Efficient absorption of stellar UV and visible radiation by TiO and VO is predicted to drive temperature inversions in the upper atmospheres of hot Jupiters. However, few inversions or detections of TiO or VO have been reported, and results are often contradictory. Using the improved ExoMol Toto line list, we searched for TiO emission in the dayside spectrum of WASP-33b using the same data in which the molecule was previously detected with an older line list at 4.8$\sigma$. We intended to confirm the molecular detection and quantify the signal improvement offered by the ExoMol Toto line list. Data from the High Dispersion Spectrograph on the Subaru Telescope was extracted and reduced in an identical manner to the previous study. Stellar and telluric contamination were then removed. High-resolution TiO emission models of WASP-33b were created using the radiative transfer code petitRADTRANS and cross-correlated with the data. We measure a 4.3$\sigma$ TiO emission signature using the ExoMol Toto models, corresponding to a WASP-33b orbital velocity semi-amplitude of $K_\mathrm{p}=252.9^{+5.0}_{-5.3}\ \mathrm{km\ s^{-1}}$ and a system velocity of $v_\mathrm{sys}=-23.0^{+4.7}_{-4.6}\ \mathrm{km\ s^{-1}}$. Injection-recovery tests using models based on the new and earlier line lists indicate that if the new models provide a perfect match to the planet spectrum, the significance of the TiO detection should have increased by a factor of $\sim$2. Although the TiO signal we find is statistically significant, comparison with previous works makes our result too ambiguous to claim a clear-cut detection. Unexpectedly, the new ExoMol Toto models provide a weaker signal than that found previously, which is offset in $K_\mathrm{p}$-$v_\mathrm{sys}$ space. This sheds some doubt on both detections, especially in light of a recently published TiO non-detection using a different dataset., Comment: 17 pages, 15 figures, Accepted for publication in A&A

We present a transmission spectrum for the warm (500-600K) sub-Neptune HD3167c obtained using the Hubble Space Telescope Wide Field Camera 3 infrared spectrograph. We combine these data, which span the 1.125-1.643 micron wavelength range, with broadband transit measurements made using Kepler/K2 (0.6-0.9 micron) and Spitzer/IRAC (4-5 micron). We find evidence for absorption by at least one of H2O, HCN, CO2, and CH4 (Bayes factor 7.4; 2.5-sigma significance), although the data precision does not allow us to unambiguously discriminate between these molecules. The transmission spectrum rules out cloud-free hydrogen-dominated atmospheres with metallicities <100x solar at >5.8-sigma confidence. In contrast, good agreement with the data is obtained for cloud-free models assuming metallicities >700x solar. However, for retrieval analyses that include the effect of clouds, a much broader range of metallicities (including subsolar) is consistent with the data, due to the degeneracy with cloud-top pressure. Self-consistent chemistry models that account for photochemistry and vertical mixing are presented for the atmosphere of HD3167c. The predictions of these models are broadly consistent with our abundance constraints, although this is primarily due to the large uncertainties on the latter. Interior structure models suggest the core mass fraction is >40%, independent of a rock or water core composition, and independent of atmospheric envelope metallicity up to 1000x solar. We also report abundance measurements for fifteen elements in the host star, showing that it has a very nearly solar composition., Comment: Accepted for publication in the Astronomical Journal

Non-rocky sub-jovian exoplanets in high irradiation environments are rare. LTT 9979b, also known as TESS Object of Interest (TOI) 193.01, is one of the few such planets discovered to date, and the first example of an ultra-hot Neptune. The planet's bulk density indicates that it has a substantial atmosphere, so to investigate its atmospheric composition and shed further light on its origin, we obtained {\it Spitzer} IRAC secondary eclipse observations of LTT 9979b at 3.6 and 4.5 $\mu$m. We combined the {\it Spitzer} observations with a measurement of the secondary eclipse in the {\it TESS} bandpass. The resulting secondary eclipse spectrum strongly prefers a model that includes CO absorption over a blackbody spectrum, incidentally making LTT 9979b the first {\it TESS} exoplanet (and the first ultra-hot Neptune) with evidence of a spectral feature in its atmosphere. We did not find evidence of a thermal inversion, at odds with expectations based on the atmospheres of similarly-irradiated hot Jupiters. We also report a nominal dayside brightness temperature of 2305 $\pm$ 141 K (based on the 3.6 $\mu$m secondary eclipse measurement), and we constrained the planet's orbital eccentricity to $e < 0.01$ at the 99.7 \% confidence level. Together with our analysis of LTT 9979b's thermal phase curves reported in a companion paper, our results set the stage for similar investigations of a larger sample of exoplanets discovered in the hot Neptune desert, investigations which are key to uncovering the origin of this population., Comment: 12 pages, 5 figures; accepted to ApJ Letters

The orbital parameters of warm Jupiters serve as a record of their formation history, providing constraints on formation scenarios for giant planets on close and intermediate orbits. Here, we report the discovery of TIC 237913194b, detected in full frame images from Sectors 1 and 2 of TESS, ground-based photometry (CHAT, LCOGT), and FEROS radial velocity time series. We constrain its mass to $M_\mathrm{P} = 1.942_{-0.091}^{+0.091}\,{\rm M_{J}} $, and its radius to $R_\mathrm{P} = 1.117_{-0.047}^{+0.054}\,{\rm R_J}$, implying a bulk density similar to Neptune's. It orbits a G-type star (${\rm M}_{\star} = 1.026_{-0.055}^{+0.057}\,{\rm M}_{\odot}$, $V = 12.1$ mag) with a period of $15.17\,$d on one of the most eccentric orbits of all known warm giants ($e \approx 0.58$). This extreme dynamical state points to a past interaction with an additional, undetected massive companion. A tidal evolution analysis showed a large tidal dissipation timescale, suggesting that the planet is not a progenitor for a hot Jupiter caught during its high-eccentricity migration. TIC 237913194b further represents an attractive opportunity to study the energy deposition and redistribution in the atmosphere of a warm Jupiter with high eccentricity., Comment: 20 pages, 14 figures. Accepted for publication in AJ

The circumstellar disk of PDS 70 hosts two forming planets, which are actively accreting gas from their environment. In this work, we report the first detection of PDS 70 b in the Br$\alpha$ and $M'$ filters with VLT/NACO, a tentative detection of PDS 70 c in Br$\alpha$, and a reanalysis of archival NACO $L'$ and SPHERE $H23$ and $K12$ imaging data. The near side of the disk is also resolved with the Br$\alpha$ and $M'$ filters, indicating that scattered light is non-negligible at these wavelengths. The spectral energy distribution of PDS 70 b is well described by blackbody emission, for which we constrain the photospheric temperature and photospheric radius to $T_\mathrm{eff}=1193 \pm 20$ K and $R=3.0 \pm 0.2$ $R_\mathrm{J}$. The relatively low bolometric luminosity, $\log(L/L_\odot) = -3.79 \pm 0.02$, in combination with the large radius, is not compatible with standard structure models of fully convective objects. With predictions from such models, and adopting a recent estimate of the accretion rate, we derive a planetary mass and radius in the range of $M_\mathrm{p}\approx 0.5-1.5$ $M_\mathrm{J}$ and $R_\mathrm{p}\approx 1-2.5$ $R_\mathrm{J}$, independently of the age and post-formation entropy of the planet. The blackbody emission, large photospheric radius, and the discrepancy between the photospheric and planetary radius suggests that infrared observations probe an extended, dusty environment around the planet, which obscures the view on its molecular composition. Finally, we derive a rough upper limit on the temperature and radius of potential excess emission from a circumplanetary disk, $T_\mathrm{eff}\lesssim256$ K and $R\lesssim245$ $R_\mathrm{J}$, but we do find weak evidence that the current data favors a model with a single blackbody component., Comment: 19 pages, 7 figures, accepted for publication in A&A

The radii of hot Jupiters are still not fully understood and all of the proposed explanations are based on the idea that these close-in giant planets possess hot interiors. We approach the radius anomaly problem by adopting a statistical approach. We infer the internal luminosity for the sample of hot Jupiters, study its effect on the interior structure, and put constraints on which mechanism is the dominant one. We develop a flexible and robust hierarchical Bayesian model that couples the interior structure of exoplanets to their observed properties. We apply the model to 314 hot Jupiters and infer the internal luminosity distribution for each planet and study at the population level ({\it i}) the mass-luminosity-radius distribution and as a function of equilibrium temperature the distributions of the ({\it ii}) heating efficiency, ({\it iii}) internal temperature, and the ({\it iv}) pressure of the radiative-convective-boundary (RCB). We find that hot Jupiters tend to have high internal luminosity leading to hot interiors. This has important consequences on the cooling rate and we find that the RCB is located at low pressures. Assuming that the ultimate source of the extra heating is the irradiation from the host star, we illustrate that the heating efficiency follows a Gaussian distribution, in agreement with previous results. We discuss our findings in the context of the proposed heating mechanisms and illustrate that ohmic dissipation, advection of potential temperature, and thermal tides are in agreement with certain trends inferred from our analysis and thus all three models can explain aspects of the observations. We provide new insights on the interior structure of hot Jupiters and show that with our current knowledge it is still challenging to firmly identify the universal mechanism driving the inflated radii., Comment: 27 pages and 12 figures. Accepted in A&A. Source code can be found at https://github.com/psarkis/bloatedHJs and data at https://www.space.unibe.ch/research/research_groups/planets_in_time/numerical_data/index_eng.html

A publicly available database of opacities for molecules of astrophysical interest, ExoMolOP, has been compiled for over 80 species, based on the latest line list data from the ExoMol, HITEMP and MoLLIST databases. These data are generally suitable for characterising high temperature exoplanet or cool stellar/substellar atmospheres, and have been computed at a variety of pressures and temperatures, with a few molecules included at room-temperature only from the HITRAN database. The data are formatted in different ways for four different exoplanet atmosphere retrieval codes; ARCiS, TauREx, NEMESIS and petitRADTRANS, and include both cross-sections (at R~=~$\frac{\lambda}{\Delta \lambda}$~=~15,000) and k-tables (at R~=~$\frac{\lambda}{\Delta \lambda}$~=~1000) for the 0.3~-~50$\mu$m wavelength region. Opacity files can be downloaded and used directly for these codes. Atomic data for alkali metals Na and K are also included, using data from the NIST database and the latest line shapes for the resonance lines. Broadening parameters have been taken from the literature where available, or from those for a known molecule with similar molecular properties where no broadening data are available. The data are available from www.exomol.com., Comment: 25 pages, 9 figures. Accepted for publication in A&A

Recent observations of the potentially habitable planets TRAPPIST-1 e, f, and g suggest that they possess large water mass fractions of possibly several tens of wt% of water, even though the host star's activity should drive rapid atmospheric escape. These processes can photolyze water, generating free oxygen and possibly desiccating the planet. After the planets formed, their mantles were likely completely molten with volatiles dissolving and exsolving from the melt. In order to understand these planets and prepare for future observations, the magma ocean phase of these worlds must be understood. To simulate these planets, we have combined existing models of stellar evolution, atmospheric escape, tidal heating, radiogenic heating, magma ocean cooling, planetary radiation, and water-oxygen-iron geochemistry. We present MagmOc, a versatile magma ocean evolution model, validated against the rocky Super-Earth GJ 1132b and early Earth. We simulate the coupled magma ocean-atmospheric evolution of TRAPPIST-1 e, f, and g for a range of tidal and radiogenic heating rates, as well as initial water contents between 1 and 100 Earth oceans. We also reanalyze the structures of these planets and find they have water mass fractions of 0-0.23, 0.01-0.21, and 0.11-0.24 for planets e, f, and g, respectively. Our model does not make a strong prediction about the water and oxygen content of the atmosphere of TRAPPIST-1 e at the time of mantle solidification. In contrast, the model predicts that TRAPPIST-1 f and g would have a thick steam atmosphere with a small amount of oxygen at that stage. For all planets that we investigated, we find that only 3-5% of the initial water will be locked in the mantle after the magma ocean solidified., Comment: Accepted for publication in Astrobiology, 24 pages, 9 figures, 7 tables

Observations suggest an abundance of water and paucity of methane in the majority of observed exoplanetary atmospheres. We isolate the effect of atmospheric processes to investigate possible causes. Previously, we studied the effect of effective temperature, surface gravity, metallicity, carbon-to-oxygen ratio, and stellar type assuming cloud-free thermochemical equilibrium and disequilibrium chemistry. However, under these assumptions, methane remains a persisting spectral feature in the transmission spectra of exoplanets over a certain parameter space, the Methane Valley. In this work we investigate the role of clouds on this domain and we find that clouds change the spectral appearance of methane in two direct ways: 1) by heating-up the photosphere of colder planets, and 2) by obscuring molecular features. The presence of clouds also affects methane features indirectly: 1) cloud heating results in more evaporation of condensates and hence releases additional oxygen, causing water dominated spectra of colder carbon-poor exoplanets, and 2) HCN/CO production results in a suppression of depleted methane features by these molecules. The presence of HCN/CO and a lack of methane could be an indication of cloud formation on hot exoplanets. Cloud heating can also deplete ammonia. Therefore, a simultaneous depletion of methane and ammonia is not unique to photochemical processes. We propose that the best targets for methane detection are likely to be massive but smaller planets with a temperature around 1450 K orbiting colder stars. We also construct Spitzer synthetic color-maps and find that clouds can explain some of the high contrast observations by IRAC's channel 1 and 2., Comment: 24 pages, 19 figures, accepted in ApJ (relevant papers: arXiv:1809.09629 and arXiv:1908.09847)

Knowing the bulk elemental abundances of an exoplanetary atmosphere is not an easy task, but it is crucial to understand the formation history of planets. The purpose of this work is to show that the observability of TiO features at optical wavelengths in the transmission spectra of hot Jupiter atmospheres is sensitive to the bulk chemical properties of the atmosphere. For that, we run a grid of chemical models which include TiO formation and destruction, for the ultra-hot Jupiter WASP-19b and an ultra-hot version of HD~209458b. We take into account non-equilibrium chemistry and changes in the temperature and pressure structure of these atmospheres caused by different C/O ratios. We calculate synthetic transmission spectra for these models, and study the relative strengths of TiO and \ce{H2O} features quantitatively. To compare with observations, we use a model independent metric for molecular abundances, $\Delta Z_{\rm TiO-H_2O}/H_{\rm eq}$ that has been previously used in observational studies of exoplanetary atmospheres. We find that with this metric we can differentiate between different chemical models and place constraints on the atmosphere's bulk carbon and oxygen abundance. From chemical considerations we expected that the TiO abundance would depend on the bulk nitrogen, however found that any change to N/H did not result in changes to the resulting TiO. We apply our method to a set of known exoplanets that have been observed in the relevant optical wavelengths and find good agreement between low resolution observations and our model for WASP-121b, marginally good agreement with WASP-79b, WASP-76b, and WASP-19b, and poorer agreement with HD 209458b. Our method can be particularly helpful for indirect studies of the bulk abundances of carbon and oxygen., Comment: 14 pages, 13 figures, accepted for publication in A&A