Search

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)

With a temperature of $\sim 285$ K WISE0855 is the coldest brown dwarf observed so far. Using the James Webb Space Telescope (JWST) we obtained observations that allow us to characterize WISE0855s atmosphere focusing on vertical variation in the water steam abundance, measuring trace gas abundances and receiving bulk parameters for this cold object. We observed the ultra cool dwarf WISE0855 using the Mid-Infrared Instrument Medium Resolution Spectrometer (MIRI/MRS) onboard JWST at a spectral resolution of up to 3750. We combined the observation with published data from the Near Infrared Spectrograph (NIRSpec) G395M and PRISM modes yielding a spectrum ranging from 0.8 to 22 um. We apply atmospheric retrievals using petitRADTRANS to measure atmospheric abundances, the pressure-temperature structure, radius and gravity of the brown dwarf. We also employ publicly available clear and cloudy self-consistent grid models to estimate bulk properties of the atmosphere such as the effective temperature, radius, gravity and metallicity. Atmospheric retrievals constrain a variable water abundance profile in the atmosphere, as predicted by equilibrium chemistry. We detect the 15NH3 isotopologue and infer a ratio of mass fraction of 14NH3/15NH3 = 332+63-43 for the clear retrieval. We measure the bolometric luminosity by integrating the presented spectrum and obtain a value of log(L/L$_{\odot}$) = -7.291+/-0.008. The detected water depletion indicates that water condenses out in the upper atmosphere due to the very low effective temperature of WISE0855. The height in the atmosphere where this occurs is covered by the MIRI/MRS data, and thus demonstrates the potential of MIRI to characterize cold gas giants atmospheres. Comparing the data to retrievals and self-consistent grid models, we do not detect signs for water ice clouds, although their spectral features have been predicted in previous studies., Comment: Submitted to A&A, 29 pages, 21 figures

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 C/O is consistent with solar, that together imply planetesimal accretion. Our GCM simulations for WASP-15b suggest that the C/O 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: 24 pages, 23 figures, 6 tables. Submitted to MNRAS

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)

HD 209458 b is one of the most studied exoplanets to date. Despite this, atmospheric characterisation studies yielded inconsistent species detections and abundances. Values reported for the C/O ratio range from 0.1 to 1.0. Of particular interest is the simultaneous detection of H2O and HCN reported by some studies using high-resolution ground-based observations, which would require the atmospheric C/O ratio to be fine-tuned to a narrow interval around 1. HCN has however not been detected from recent space-based observations. We aim to provide an independent study of HD 209458 b's atmosphere with high-resolution observations, in order to infer the presence of several species, including H2O and HCN. We observed four primary transits of HD 209458 b at a high resolution (R=92000) with CRIRES+ in the near infrared (band H, 1.4--1.8 um). After reducing the data with pycrires, we prepared the data using the SysRem algorithm and performed a cross-correlation (CCF) analysis of the transmission spectra. We also compared the results with those obtained from simulated datasets constructed by combining the Exo-REM self-consistent model with the petitRADTRANS package. Combining the four transits, we detect H2O with a signal-to-noise CCF metric of 8.7. This corresponds to a signal emitted at $K_p=151.3^{+31.1}_{-23.4}$ km/s and blueshifted by $-6^{+1}_{-2}$ km/s, consistent with what is expected for HD 209458 b. We do not detect any other species among C2H2, CH4, CO, CO2, H2S, HCN, and NH3. Comparing this with our simulated datasets, this result is consistent with a C/O ratio of 0.1 and an opaque cloud top pressure of 50 Pa, at a 3 times solar metallicity. This would also be consistent with recent JWST observations. However, none of the simulated results obtained with a bulk C/O ratio of 0.8, a value suggested by previous studies using GIANO-B and CRIRES, are consistent with our observations., Comment: 15 pages, 13 figures, 5 tables, accepted for publication in Astronomy & Astrophysics

The first generation of ELT instruments includes an optical-infrared high-resolution spectrograph, indicated as ELT-HIRES and recently christened ANDES (ArmazoNes high Dispersion Echelle Spectrograph). ANDES consists of three fibre-fed spectrographs ([U]BV, RIZ, YJH) providing a spectral resolution of $\sim$100,000 with a minimum simultaneous wavelength coverage of 0.4-1.8 $\mu$m with the goal of extending it to 0.35-2.4 $\mu$m with the addition of a U arm to the BV spectrograph and a separate K band spectrograph. It operates both in seeing- and diffraction-limited conditions and the fibre feeding allows several, interchangeable observing modes including a single conjugated adaptive optics module and a small diffraction-limited integral field unit in the NIR. Modularity and fibre-feeding allow ANDES to be placed partly on the ELT Nasmyth platform and partly in the Coud\'e room. ANDES has a wide range of groundbreaking science cases spanning nearly all areas of research in astrophysics and even fundamental physics. Among the top science cases, there are the detection of biosignatures from exoplanet atmospheres, finding the fingerprints of the first generation of stars, tests on the stability of Nature's fundamental couplings, and the direct detection of the cosmic acceleration. The ANDES project is carried forward by a large international consortium, composed of 35 Institutes from 13 countries, forming a team of almost 300 scientists and engineers which include the majority of the scientific and technical expertise in the field that can be found in ESO member states., Comment: SPIE astronomical telescope and instrumentation 2024, in press

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

Young brown dwarfs exhibit atmospheric characteristics similar to those of super-Jupiters, providing a unique opportunity to study planetary atmospheres. The ESO SupJup Survey, utilizing CRIRES$^+$ on the Very Large Telescope, aims to assess the role of $^{12}$C/$^{13}$C as a formation tracer. We present observations of three young brown dwarfs: 2MASS J12003792-7845082, TWA 28, and 2MASS J08561384-1342242, with the goal of constraining their chemical compositions, thermal profiles, surface gravities, spin rotations, and $^{12}$C/$^{13}$C. Atmospheric retrievals of CRIRES$^+$ K-band spectra were conducted using the radiative transfer code petitRADTRANS coupled with the Bayesian inference algorithm MultiNest, resulting in a detailed characterization of the atmospheres of these objects. We report the volume mixing ratios of main molecular and atomic species, including the novel detection of hydrogen fluoride (HF) in a brown dwarf's atmosphere, and determine $^{12}$C/$^{13}$C values of $81^{+28}_{-19}$ and $79^{+20}_{-14}$ in the atmospheres of TWA 28 and J0856, respectively, with strong significance ($>3\sigma$). Tentative evidence ($\sim 2\sigma$) of $^{13}$C in J1200 was found, with $^{12}$C/$^{13}$C = $114^{+69}_{-33}$, along with $^{18}$O detected at moderate significance in J0856 (3.3$\sigma$) and TWA 28 (2.1$\sigma$). The retrieved thermal profiles indicate hot atmospheres (2300-2600 K) with low surface gravities and slow spins, consistent with young objects. The consistent carbon isotope ratios among the three objects, showing no significant deviation from the local ISM, suggest a fragmentation-based formation mechanism similar to star formation. The tentative detection of $^{18}$O in two objects highlights the potential of high-resolution spectroscopy to probe additional isotope ratios, such as $^{16}$O/$^{18}$O, in the atmospheres of brown dwarfs and super-Jupiters., Comment: Published in A&A Volume 689, September 2024

Following the recommendations to NASA and ESA, the search for life on exoplanets will be a priority in the next decades. Two direct imaging space mission concepts are being developed: the Habitable Worlds Observatory (HWO) and the Large Interferometer for Exoplanets (LIFE). HWO focuses on reflected light spectra in the ultraviolet/visible/near-infrared (UV/VIS/NIR), while LIFE captures the mid-infrared (MIR) emission of temperate exoplanets. We assess the potential of HWO and LIFE in characterizing a cloud-free Earth twin orbiting a Sun-like star at 10 pc, both separately and synergistically, aiming to quantify the increase in information from joint atmospheric retrievals on a habitable planet. We perform Bayesian retrievals on simulated data from an HWO-like and a LIFE-like mission separately, then jointly, considering the baseline spectral resolutions currently assumed for these concepts and using two increasingly complex noise simulations. HWO would constrain H$_2$O, O$_2$, and O$_3$, in the atmosphere, with ~ 100 K uncertainty on the temperature profile. LIFE would constrain CO$_2$, H$_2$O, O$_3$ and provide constraints on the thermal atmospheric structure and surface temperature (~ 10 K uncertainty). Both missions would provide an upper limit on CH$_4$. Joint retrievals on HWO and LIFE data would accurately define the atmospheric thermal profile and planetary parameters, decisively constrain CO$_2$, H$_2$O, O$_2$, and O$_3$, and weakly constrain CO and CH$_4$. The detection significance is greater or equal to single-instrument retrievals. Both missions provide specific information to characterize a terrestrial habitable exoplanet, but the scientific yield is maximized with synergistic UV/VIS/NIR+MIR observations. Using HWO and LIFE together will provide stronger constraints on biosignatures and life indicators, potentially transforming the search for life in the universe., Comment: 16 pages (main text, incl. 12 figures) + appendix; accepted for publication in A&A (current version: post 1st revision). Thirteenth paper of LIFE telescope series

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

Since 2019, GRAVITY has provided direct observations of giant planets and brown dwarfs at separations of down to 95 mas from the host star. Some of these observations have provided the first direct confirmation of companions previously detected by indirect techniques (astrometry and radial velocities). We want to improve the observing strategy and data reduction in order to lower the inner working angle of GRAVITY in dual-field on-axis mode. We also want to determine the current limitations of the instrument when observing faint companions with separations in the 30-150 mas range. To improve the inner working angle, we propose a fiber off-pointing strategy during the observations to maximize the ratio of companion-light-to-star-light coupling in the science fiber. We also tested a lower-order model for speckles to decouple the companion light from the star light. We then evaluated the detection limits of GRAVITY using planet injection and retrieval in representative archival data. We compare our results to theoretical expectations. We validate our observing and data-reduction strategy with on-sky observations; first in the context of brown dwarf follow-up on the auxiliary telescopes with HD 984 B, and second with the first confirmation of a substellar candidate around the star Gaia DR3 2728129004119806464. With synthetic companion injection, we demonstrate that the instrument can detect companions down to a contrast of $8\times 10^{-4}$ ($\Delta \mathrm{K}= 7.7$ mag) at a separation of 35 mas, and a contrast of $3\times 10^{-5}$ ($\Delta \mathrm{K}= 11$ mag) at 100 mas from a bright primary (K<6.5), for 30 min exposure time. With its inner working angle and astrometric precision, GRAVITY has a unique reach in direct observation parameter space. This study demonstrates the promising synergies between GRAVITY and Gaia for the confirmation and characterization of substellar companions., Comment: 16 pages, 14 figures. Submitted to A&A

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

It has been proposed that the distinct formation and evolution of exoplanets and brown dwarfs may affect the chemical and isotopic content of their atmospheres. Recent work has indeed shown differences in the $^{12}$C/$^{13}$C isotope ratio, provisionally attributed to the top-down formation of brown dwarfs and the core accretion pathway of super-Jupiters. The ESO SupJup Survey aims to disentangle the formation pathways of isolated brown dwarfs and planetary-mass companions using chemical and isotopic tracers. The survey uses high-resolution spectroscopy with the recently upgraded VLT/CRIRES$^+$ spectrograph, covering a total of 49 targets. Here, we present the first results: an atmospheric characterisation of DENIS J0255-4700, an isolated brown dwarf near the L-T transition. We analyse its K-band spectrum using a retrieval framework where the radiative transfer code petitRADTRANS is coupled to PyMultiNest. Gaussian Processes are employed to model inter-pixel correlations and we adopt an updated parameterisation of the PT-profile. Abundances of CO, H$_2$O, CH$_4$, and NH$_3$ are retrieved for this fast-rotating L-dwarf. The ExoMol H$_2$O line list provides a significantly better fit than that of HITEMP. A free-chemistry retrieval is strongly favoured over equilibrium chemistry, caused by an under-abundance of CH$_4$. The free-chemistry retrieval constrains a super-solar C/O-ratio of $\sim0.68$ and a solar metallicity. We find tentative evidence ($\sim3\sigma$) for the presence of $^{13}$CO, with a constraint on the isotope ratio of $\mathrm{^{12}C/^{13}C}=184^{+61}_{-40}$ and a lower limit of $\gtrsim97$, suggesting a depletion of $^{13}$C compared to the interstellar medium ($\sim68$). High-resolution, high signal-to-noise K-band spectra provide an excellent means to constrain the chemistry and isotopic content of sub-stellar objects, as is the main objective of the ESO SupJup Survey., Comment: Submitted to A&A

Context. High-resolution spectroscopy has the potential to drive a better understanding of the atmospheric composition, physics, and dynamics of young exoplanets and brown dwarfs, bringing clear insights into the formation channel of individual objects. Aims. Using the Keck Planet Imager and Characterizer (KPIC; R = 35,000), we aim to characterize a young brown dwarf HD 984 B. By measuring its C/O and 12CO/13CO ratios, we expect to gain new knowledge about its origin by confirming the difference in the formation pathways between brown dwarfs and super-Jupiters. Methods. We analysed the KPIC high-resolution spectrum (2.29-2.49 {\mu}m) of HD 984 B using an atmospheric retrieval framework based on nested sampling and petitRADTRANS, using both clear and cloudy models. Results. Using our best-fit model, we find C/O = 0.50+0.01-0.01 (0.01 is the statistical error) for HD 984 B which agrees with that of its host star within 1{\sigma} (0.40+0.20-0.20). We also retrieve an isotopolog 12CO/13CO ratio of 98+20-25 in its atmosphere, which is similar to that of the Sun. In addition, HD 984 B has a substellar metallicity with [Fe/H] = -0.62+0.02-0.02. Finally, we find that most of the retrieved parameters are independent of our choice of retrieval model. Conclusions. From our measured C/O and 12CO/13CO, the favored formation mechanism of HD 984 B seems to be via gravitational collapse or disk instability and not core accretion, which is a favored formation mechanism for giant exoplanets with m < 13 MJup and semimajor axis between 10 and 100 au. However, with only a few brown dwarfs with a measured 12CO/13CO ratio, similar analyses using high-resolution spectroscopy will become essential in order to determine planet formation processes more precisely., Comment: 14 pages, 8 figures

With four companions at separations from 16 to 71 au, HR 8799 is a unique target for direct imaging, presenting an opportunity for the comparative study of exoplanets with a shared formation history. Combining new VLTI/GRAVITY observations obtained within the ExoGRAVITY program with archival data, we perform a systematic atmospheric characterisation of all four planets. We explore different levels of model flexibility to understand the temperature structure, chemistry and clouds of each planet using both petitRADTRANS atmospheric retrievals and fits to self-consistent radiative-convective equilibrium models. Using Bayesian Model Averaging to combine multiple retrievals, we find that the HR 8799 planets are highly enriched in metals, with [M/H] $\gtrsim$1, and have stellar to super-stellar C/O ratios. The C/O ratio increases with increasing separation from $0.55^{+0.12}_{-0.10}$ for d to $0.78^{+0.03}_{-0.04}$ for b, with the exception of the innermost planet which has a C/O ratio of $0.87\pm0.03$. By retrieving a quench pressure and using a disequilibrium chemistry model we derive vertical mixing strengths compatible with predictions for high-metallicity, self-luminous atmospheres. Bayesian evidence comparisons strongly favour the presence of HCN in HR 8799 c and e, as well as CH$_{4}$ in HR 8799 c, with detections at $>5\sigma$ confidence. All of the planets are cloudy, with no evidence for patchiness. The clouds of c, d and e are best fit by silicate clouds lying above a deep iron cloud layer, while the clouds of the cooler HR 8799 b are more likely composed of Na$_{2}$S. With well defined atmospheric properties, future exploration of this system is well positioned to unveil further detail in these planets, extending our understanding of the composition, structure, and formation history of these siblings., Comment: 45 pages, 25 figures, accepted by Astronomy & Astrophysics

Precise mass constraints are vital for the characterisation of brown dwarfs and exoplanets. Here we present how the combination of data obtained by Gaia and GRAVITY can help enlarge the sample of substellar companions with measured dynamical masses. We show how the Non-Single-Star (NSS) two-body orbit catalogue contained in Gaia DR3 can be used to inform high-angular-resolution follow-up observations with GRAVITY. Applying the method presented in this work to eight Gaia candidate systems, we detect all eight predicted companions, seven of which were previously unknown and five are of a substellar nature. Among the sample is Gaia DR3 2728129004119806464 B, which - detected at an angular separation of (34.01 $\pm$ 0.15) mas from the host - is the closest substellar companion ever imaged. This translates to a semi-major axis of (0.938 $\pm$ 0.023) AU. WT 766 B, detected at a greater angular separation, was confirmed to be on an orbit exhibiting an even smaller semi-major axis of (0.676 $\pm$ 0.008) AU. The GRAVITY data were then used to break the host-companion mass degeneracy inherent to the Gaia NSS orbit solutions as well as to constrain the orbital solutions of the respective target systems. Knowledge of the companion masses enabled us to further characterise them in terms of their ages, effective temperatures, and radii via the application of evolutionary models. The inferred ages exhibit a distinct bias towards values younger than what is to be expected based on the literature. The results serve as an independent validation of the orbital solutions published in the NSS two-body orbit catalogue and show that the combination of astrometric survey missions and high-angular-resolution direct imaging holds great promise for efficiently increasing the sample of directly imaged companions in the future, especially in the light of Gaia's upcoming DR4 and the advent of GRAVITY+., Comment: Accepted for publication in Astronomy and Astrophysics

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

Young, low-mass Brown Dwarfs orbiting early-type stars, with low mass ratios ($q\lesssim0.01$), appear intrinsically rare and present a formation dilemma: could a handful of these objects be the highest mass outcomes of ``planetary" formation channels (bottom up within a protoplanetary disk), or are they more representative of the lowest mass ``failed binaries" (formed via disk fragmentation, or core fragmentation)? Additionally, their orbits can yield model-independent dynamical masses, and when paired with wide wavelength coverage and accurate system age estimates, can constrain evolutionary models in a regime where the models have a wide dispersion depending on initial conditions. We present new interferometric observations of the $16\,\mathrm{Myr}$ substellar companion HD~136164~Ab (HIP~75056~Ab) with VLTI/GRAVITY and an updated orbit fit including proper motion measurements from the Hipparcos-Gaia Catalogue of Accelerations. We estimate a dynamical mass of $35\pm10\,\mathrm{M_J}$ ($q\sim0.02$), making HD~136164~Ab the youngest substellar companion with a dynamical mass estimate. The new mass and newly constrained orbital eccentricity ($e=0.44\pm0.03$) and separation ($22.5\pm1\,\mathrm{au}$) could indicate that the companion formed via the low-mass tail of the Initial Mass Function. Our atmospheric fit to the \texttt{SPHINX} M-dwarf model grid suggests a sub-solar C/O ratio of $0.45$, and $3\times$ solar metallicity, which could indicate formation in the circumstellar disk via disk fragmentation. Either way, the revised mass estimate likely excludes ``bottom-up" formation via core accretion in the circumstellar disk. HD~136164~Ab joins a select group of young substellar objects with dynamical mass estimates; epoch astrometry from future \textit{Gaia} data releases will constrain the dynamical mass of this crucial object further., Comment: Accepted for publication in The Astronomical Journal. 9 figures, 3 tables

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

The multi planet system HR8799 is the first target observed with MIRI's coronagraphs as part of the MIRI-EC Guaranteed Time Observations exoplanets programme in Nov. 2022. We obtained deep observations in three coronagraphic filters from 10 to 15mic (F1065C, F1140C, F1550C), and one standard imaging filter at 20 mic (F2100W), with the goal to extract the photometry of the four planets, as well as to detect and investigate the distribution of circumstellar dust. Using dedicated observations of a reference star, we tested several algorithms to subtract the stellar diffraction pattern while preserving the fluxes of planets, which can be significantly affected by over-subtraction. Measuring correctly the planet's flux values requires accounting for the attenuation by the coronagraphs as a function of their position, and to estimate the normalisation with respect to the central star. We tested several procedures to derive averaged photometric values and error bars. These observations have enabled us to obtain two main results. First of all, the four planets in the system are well recovered, and their mid-IR fluxes, combined with near-IR flux values from the literature, are compared to two exoplanet atmosphere models, ATMO and Exo-REM. As a main outcome, the MIRI photometric data points imply larger radii (0.86 or 1.07 RJ for planet b) and cooler temperatures (950 or 1100 K for planet b), especially for planet b, in better agreement with evolutionary models. Second of all, these JWST/MIRI coronagraphic data also deliver the first spatially resolved detection of the inner warm debris disk, the radius of which is constrained to about 15 au, with flux densities comparable, but lower than former unresolved spectroscopic measurements with Spitzer. abridged..., Comment: submitted on Sep; 8th, 2023

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: 20 pages, 12 figures, submitted to 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

Giant exoplanets have been directly imaged over orders of magnitude of orbital separations, prompting theoretical and observational investigations of their formation pathways. In this paper, we present new VLTI/GRAVITY astrometric data of HIP 65426 b, a cold, giant exoplanet which is a particular challenge for most formation theories at a projected separation of 92 au from its primary. Leveraging GRAVITY's astrometric precision, we present an updated eccentricity posterior that disfavors large eccentricities. The eccentricity posterior is still prior-dependent, and we extensively interpret and discuss the limits of the posterior constraints presented here. We also perform updated spectral comparisons with self-consistent forward-modeled spectra, finding a best fit ExoREM model with solar metallicity and C/O=0.6. An important caveat is that it is difficult to estimate robust errors on these values, which are subject to interpolation errors as well as potentially missing model physics. Taken together, the orbital and atmospheric constraints paint a preliminary picture of formation inconsistent with scattering after disk dispersal. Further work is needed to validate this interpretation. Analysis code used to perform this work is available at https://github.com/sblunt/hip65426., Comment: 28 pages, 14 figures. Revised and resubmitted to AJ

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

Tension remains between the observed and modeled properties of substellar objects, but objects in binary orbits, with known dynamical masses can provide a way forward. HD 72946 B is a recently imaged brown dwarf companion to the nearby, solar type star. We achieve $\sim100~\mu\mathrm{as}$ relative astrometry of HD 72946 B in the K-band using VLTI/GRAVITY, unprecedented for a benchmark brown dwarf. We fit an ensemble of measurements of the orbit using orbitize! and derive a strong dynamical mass constraint $\mathrm{M_B}=69.5\pm0.5~\mathrm{M_{Jup}}$ assuming a strong prior on the host star mass $\mathrm{M_A}=0.97\pm0.01~\mathrm{M_\odot}$ from an updated stellar analysis. We fit the spectrum of the companion to a grid of self-consistent BT-Settl-CIFIST model atmospheres, and perform atmospheric retrievals using petitRADTRANS. A dynamical mass prior only marginally influences the sampled distribution on effective temperature, but has a large influence on the surface gravity and radius, as expected. The dynamical mass alone does not strongly influence retrieved pressure-temperature or cloud parameters within our current retrieval setup. Independent of cloud prescription and prior assumptions, we find agreement within $\pm2\,\sigma$ between the C/O ratio of the host ($0.52\pm0.05)$ and brown dwarf ($0.43$ to $0.63$), as expected from a molecular cloud collapse formation scenario, but our retrieved metallicities are implausibly high ($0.6-0.8$) in light of an excellent agreement of the data with the solar abundance model grid. Future work on our retrieval framework will seek to resolve this tension. Additional study of low surface-gravity objects is necessary to assess the influence of a dynamical mass prior on atmospheric analysis., Comment: Accepted for publication in the Astrophysical Journal. 32 pages, 14 figures, 7 tables. v2&3 correct errors in co-author's affiliations, figure rendering, and some grant acknowledgements

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

The Large Interferometer For Exoplanets (LIFE) initiative aims to develop a space based mid-infrared (MIR) nulling interferometer to measure the thermal emission spectra of temperate terrestrial exoplanets. We investigate how well LIFE could characterize a cloudy Venus-twin exoplanet to: (1) test our retrieval routine on a realistic non-Earth-like MIR spectrum of a known planet, (2) investigate how clouds impact retrievals, (3) refine the LIFE requirements derived in previous Earth-centered studies. We run retrievals for simulated LIFE observations of a Venus-twin exoplanet orbiting a Sun-like star located 10 pc from the observer. By assuming different models (cloudy and cloud-free) we analyze the performance as a function of the quality of the LIFE observation. This allows us to determine how well atmosphere and clouds are characterizable depending on the quality of the spectrum. Our study shows that the current minimal resolution ($R=50$) and signal-to-noise ($S/N=10$ at $11.2\mu$m) requirements for LIFE suffice to characterize the structure and composition of a Venus-like atmosphere above the cloud deck if an adequate model is chosen. However, we cannot infer cloud properties. The accuracy of the retrieved planet radius ($R_{pl}$), equilibrium temperature ($T_{eq}$), and Bond albedo ($A_B$) depend on the choice of model. Generally, a cloud-free model performs best and thus the presence of clouds cannot be inferred. This model dependence of retrieval results emphasizes the importance of developing a community-wide best-practice for atmospheric retrieval studies. If we consider higher quality spectra (especially $S/N=20$), we can infer the presence of clouds and pose first constraints on their structure., Comment: Accepted for publication in A&A; Updates: LIFE-series number changed from 10 to 9, additional references added to the discussion, language editing; 15 pages (main text incl. 8 figures and 6 tables) + appendix; comments are welcome. This paper is part of a series on the LIFE telescope. Related series papers: arXiv:2101.07500, arXiv:2203.00471, arXiv:2112.02054, arXiv:2204.10041

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

[Full abstract in the paper] We aim to better constrain the atmospheric properties of the directly imaged exoplanet 51~Eri~b by using a retrieval approach on higher signal-to-noise data than previously reported. In this context, we also compare the results of using the atmospheric retrieval code \texttt{petitRADTRANS} vs a self-consistent model to fit atmospheric parameters. We present a higher signal-to-noise $YH$ spectrum of the planet and revised $K1K2$ photometry (M$_{K1} = 15.11 \pm 0.04$ mag, M$_{K2} = 17.11 \pm 0.38$ mag). The best-fit parameters obtained using an atmospheric retrieval differ from previous results using self-consistent models. In general, we find that our solutions tend towards cloud-free atmospheres (e.g. log $\tau_{\rm clouds} = -5.20 \pm 1.44$). For our ``nominal'' model with new data, we find a lower metallicity ([Fe/H] $= 0.26\pm$0.30 dex) and C/O ratio ($0.38\pm0.09$), and a slightly higher effective temperature (T$_{\rm{eff}} = 807\pm$45 K) than previous studies. The surface gravity (log $g = 4.05\pm0.37$) is in agreement with the reported values in the literature within uncertainties. We estimate the mass of the planet to be between 2 and 4 M$_{\rm{Jup}}$. When comparing with self-consistent models, we encounter a known correlation between the presence of clouds and the shape of the $P-T$ profiles. Our findings support the idea that results from atmospheric retrievals should not be discussed in isolation, but rather along with self-consistent temperature structures obtained using the retrieval's best-fit parameters., Comment: Accepted for publication in A&A. 21 pages, 7 figures in the main text and 9 figures in the Appendix

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)

Ultra-hot Jupiters are highly irradiated gas giant exoplanets on close-in orbits around their host stars. We analyzed high-resolution spectra from CARMENES, HARPS-N, and ESPaDOnS taken over eight observation nights to study the emission spectrum of WASP-33b and draw conclusions about its atmosphere. By applying the cross-correlation technique, we detected the spectral signatures of Ti I, V I, and a tentative signal of Ti II for the first time via emission spectroscopy. These detections are an important finding because of the fundamental role of Ti- and V-bearing species in the planetary energy balance. Moreover, we assessed and confirm the presence of OH, Fe I, and Si I from previous studies. The spectral lines are all detected in emission, which unambiguously proves the presence of an inverted temperature profile in the planetary atmosphere. By performing retrievals on the emission lines of all the detected species, we determined a relatively weak atmospheric thermal inversion extending from approximately 3400 K to 4000 K. We infer a supersolar metallicity close to 1.5 dex in the planetary atmosphere, and find that its emission signature undergoes significant line broadening with a Gaussian FWHM of about 4.5 km/s. Also, we find that the atmospheric temperature profile retrieved at orbital phases far from the secondary eclipse is about 300 K to 700 K cooler than that measured close to the secondary eclipse, which is consistent with different day- and nightside temperatures. Moreover, retrievals performed on the emission lines of the individual chemical species lead to consistent results, which gives additional confidence to our retrieval method. Increasing the number of species included in the retrieval and expanding the set of retrieved atmospheric parameters will further advance our understanding of exoplanet atmospheres., Comment: Accepted for publication in A&A

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

Long term precise radial velocity (RV) monitoring of the nearby star HD206893, as well as anomalies in the system proper motion, have suggested the presence of an additional, inner companion in the system. Here we describe the results of a multi-epoch search for the companion responsible for this RV drift and proper motion anomaly using the VLTI/GRAVITY instrument. Utilizing information from ongoing precision RV measurements with the HARPS spectrograph, as well as Gaia host star astrometry, we report a high significance detection of the companion HD206893c over three epochs, with clear evidence for Keplerian orbital motion. Our astrometry with $\sim$50-100 $\mu$arcsec precision afforded by GRAVITY allows us to derive a dynamical mass of 12.7$^{+1.2}_{-1.0}$ M$_{\rm Jup}$ and an orbital separation of 3.53$^{+0.08}_{-0.06}$ au for HD206893c. Our fits to the orbits of both companions in the system utilize both Gaia astrometry and RVs to also provide a precise dynamical estimate of the previously uncertain mass of the B component, and therefore derive an age of $155\pm15$ Myr. We find that theoretical atmospheric/evolutionary models incorporating deuterium burning for HD206893c, parameterized by cloudy atmospheres provide a good simultaneous fit to the luminosity of both HD206893B and c. In addition to utilizing long-term RV information, this effort is an early example of a direct imaging discovery of a bona fide exoplanet that was guided in part with Gaia astrometry. Utilizing Gaia astrometry is expected to be one of the primary techniques going forward to identify and characterize additional directly imaged planets. Lastly, this discovery is another example of the power of optical interferometry to directly detect and characterize extrasolar planets where they form at ice-line orbital separations of 2-4\,au., Comment: Accepted to A&A

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

Ultra-hot Jupiters (UHJs) are expected to possess temperature inversion layers in their dayside atmospheres. Recent thermal emission observations have discovered several atomic and molecular species along with temperature inversions in UHJs. We observed the thermal emission spectra of two UHJs (WASP-33b and WASP-189b) with the GIANO-B high-resolution near-infrared spectrograph. Using the cross-correlation technique, we detected carbon monoxide (CO) in the dayside atmospheres of both planets. The detected CO lines are in emission, which agrees with previous discoveries of iron emission lines and temperature inversions in the two planets. This is the first detection of CO lines in emission with high-resolution spectroscopy. Further retrieval work combining the CO lines with other spectral features will enable a comprehensive understanding of the atmospheric properties such as temperature structures and C/O ratios. The detected CO and iron emission lines of WASP-189b have redshifted radial velocities of several km/s, which likely originate from a dayside to nightside wind in its atmosphere. Such a redshifted velocity has not been detected for the emission lines of WASP-33b, suggesting that the atmospheric circulation patterns of the two UHJs may be different., Comment: 8 pages, 9 figures; accepted for publication in A&A Letter

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