Search

The GPU-based High-order adaptive OpticS Testbench (GHOST) at the European Southern Observatory (ESO) is a new 2-stage extreme adaptive optics (XAO) testbench at ESO. The GHOST is designed to investigate and evaluate new control methods (machine learning, predictive control) for XAO which will be required for instruments such as the Planetary Camera and Spectrograph of ESOs Extremely Large Telescope. The first stage corrections are performed in simulation, with the residual wavefront error at each iteration saved. The residual wavefront errors from the first stage are then injected into the GHOST using a spatial light modulator. The second stage correction is made with a Boston Michromachines Corporation 492 actuator deformable mirror and a pyramid wavefront sensor. The flexibility of the bench also opens it up to other applications, one such application is investigating the flip-flop modulation method for the pyramid wavefront sensor., Comment: spie astronomical telescopes + instrumentation 2022

We present the database of potential targets for the Large Interferometer For Exoplanets (LIFE), a space-based mid-infrared nulling interferometer mission proposed for the Voyage 2050 science program of the European Space Agency (ESA). The database features stars, their planets and disks, main astrophysical parameters, and ancillary observations. It allows users to create target lists based on various criteria to predict, for instance, exoplanet detection yields for the LIFE mission. As such, it enables mission design trade-offs, provides context for the analysis of data obtained by LIFE, and flags critical missing data. Work on the database is in progress, but given its relevance to LIFE and other space missions, including the Habitable Worlds Observatory (HWO), we present its main features here. A preliminary version of the LIFE database is publicly available on the German Astrophysical Virtual Observatory (GAVO)., Comment: RNAAS published, 3 pages, 1 figure

Inferring atmospheric properties of exoplanets from observed spectra is key to understanding their formation, evolution, and habitability. Since traditional Bayesian approaches to atmospheric retrieval (e.g., nested sampling) are computationally expensive, a growing number of machine learning (ML) methods such as neural posterior estimation (NPE) have been proposed. We seek to make ML-based atmospheric retrieval (1) more reliable and accurate with verified results, and (2) more flexible with respect to the underlying neural networks and the choice of the assumed noise models. First, we adopt flow matching posterior estimation (FMPE) as a new ML approach to atmospheric retrieval. FMPE maintains many advantages of NPE, but provides greater architectural flexibility and scalability. Second, we use importance sampling (IS) to verify and correct ML results, and to compute an estimate of the Bayesian evidence. Third, we condition our ML models on the assumed noise level of a spectrum (i.e., error bars), thus making them adaptable to different noise models. Both our noise level-conditional FMPE and NPE models perform on par with nested sampling across a range of noise levels when tested on simulated data. FMPE trains about 3 times faster than NPE and yields higher IS efficiencies. IS successfully corrects inaccurate ML results, identifies model failures via low efficiencies, and provides accurate estimates of the Bayesian evidence. FMPE is a powerful alternative to NPE for fast, amortized, and parallelizable atmospheric retrieval. IS can verify results, thus helping to build confidence in ML-based approaches, while also facilitating model comparison via the evidence ratio. Noise level conditioning allows design studies for future instruments to be scaled up, for example, in terms of the range of signal-to-noise ratios., Comment: Accepted for publication in Astronomy & Astrophysics

The increased brightness temperature of young rocky protoplanets during their magma ocean epoch makes them potentially amenable to atmospheric characterization to distances from the solar system far greater than thermally equilibrated terrestrial exoplanets, offering observational opportunities for unique insights into the origin of secondary atmospheres and the near surface conditions of prebiotic environments. The Large Interferometer For Exoplanets (LIFE) mission will employ a space-based mid-infrared nulling interferometer to directly measure the thermal emission of terrestrial exoplanets. Here, we seek to assess the capabilities of various instrumental design choices of the LIFE mission concept for the detection of cooling protoplanets with transient high-temperature magma ocean atmospheres, in young stellar associations in particular. Using the LIFE mission instrument simulator (LIFEsim) we assess how specific instrumental parameters and design choices, such as wavelength coverage, aperture diameter, and photon throughput, facilitate or disadvantage the detection of protoplanets. We focus on the observational sensitivities of distance to the observed planetary system, protoplanet brightness temperature using a blackbody assumption, and orbital distance of the potential protoplanets around both G- and M-dwarf stars. Our simulations suggest that LIFE will be able to detect (S/N $\geq$ 7) hot protoplanets in young stellar associations up to distances of $\approx$100 pc from the solar system for reasonable integration times (up to $\sim$hours). Detection of an Earth-sized protoplanet orbiting a solar-sized host star at 1 AU requires less than 30 minutes of integration time. M-dwarfs generally need shorter integration times. The contribution from wavelength regions $<$6 $\mu$m is important for decreasing the detection threshold and discriminating emission temperatures., Comment: 18 pages, 19 figures; accepted for publication in A&A

Atmospheric retrievals are widely used to constrain exoplanet properties from observed spectra. We investigate how the common nonphysical retrieval assumptions of vertically constant molecule abundances and cloud-free atmospheres affect our characterization of an exo-Earth (an Earth-twin orbiting a Sun-like star). Specifically, we use a state-of-the-art retrieval framework to explore how assumptions for the $\mathrm{H_2O}$ profile and clouds affect retrievals. In a first step, we validate different retrieval models on a low-noise simulated 1D mid-infrared (MIR) spectrum of Earth. Thereafter, we study how these assumptions affect the characterization of Earth with the Large Interferometer For Exoplanets (LIFE). We run retrievals on LIFE mock observations based on real disk-integrated MIR Earth spectra. The performance of different retrieval models is benchmarked against ground truths derived from remote sensing data. We show that assumptions for the $\mathrm{H_2O}$ abundance and clouds directly affect our characterization. Overall, retrievals that use physically motivated models for the $\mathrm{H_2O}$ profile and clouds perform better on the empirical Earth data. For observations of Earth with LIFE, they yield accurate estimates for the radius, pressure-temperature structure, and the abundances of $\mathrm{CO_2}$, $\mathrm{H_2O}$, and $\mathrm{O_3}$. Further, at $R=100$, a reliable and bias-free detection of the biosignature $\mathrm{CH_4}$ becomes feasible. We conclude that the community must use a diverse range of models for temperate exoplanet atmospheres to build an understanding of how different retrieval assumptions can affect the interpretation of exoplanet spectra. This will enable the characterization of distant habitable worlds and the search for life with future space-based instruments., Comment: Published in ApJ; 17 pages (main text incl. 9 figures and 5 tables) + appendix; comments are welcome

The Enhanced Resolution Imager and Spectrograph (ERIS) is the new near-infrared instrument at the VLT-UT4. ERIS replaces and extends the observational capabilities formerly provided by SINFONI and NACO: integral field spectroscopy at 1 - 2.5 $\mu$m, imaging at 1 - 5 $\mu$m with several options for high-contrast imaging, and long-slit spectroscopy. In particular, a vortex coronagraph is now available for high contrast observations at L and M band. It is implemented using annular groove (or vortex) phase masks (one for each of the L and M bands) in a focal plane, and a Lyot stop in a downstream pupil plane. The vortex coronagraph has a discovery space starting already at $\sim$1$\lambda/D$, and works well in broadbands. However, to reach its optimal performance, it is critical to correct for slow pointing errors onto the vortex phase mask, which mandates a dedicated pointing control strategy. To do so, a control loop based on the QACITS algorithm has been developed and commissioned for ERIS. Good pointing stability is now regularly achieved with errors between 0.01 and 0.02 $\lambda/D$ and a correction rate of 0.2 Hz. In this contribution, we first review the design of the ERIS vortex coronagraph. We then detail the implementation of the QACITS algorithm describing the entire observing sequence, including the calibration steps, the initial centering, and the stabilization during the observing template. We then discuss performance based on commissioning data in terms of pointing accuracy and stability. Finally, we present post-processed contrast curves obtained during commissioning and compare them with NACO vortex data, showing a significant improvement of about 1 mag at all separations., Comment: 15 pages, 13 figures, paper presented at SPIE Astronomical Telescopes + Instrumentation 2024

The Mid-infrared ELT Imager and Spectrograph (METIS) is one of the first-generation scientific instruments for the ELT, built under the supervision of ESO by a consortium of research institutes across and beyond Europe. Designed to cover the 3 to 13 $\mu$m wavelength range, METIS had its final design reviewed in Fall 2022, and has then entered in earnest its manufacture, assembly, integration, and test (MAIT) phase. Here, we present the final design of the METIS high-contrast imaging (HCI) modes. We detail the implementation of the two main coronagraphic solutions selected for METIS, namely the vortex coronagraph and the apodizing phase plate, including their combination with the high-resolution integral field spectrograph of METIS, and briefly describe their respective backup plans (Lyot coronagraph and shaped pupil plate). We then describe the status of the MAIT phase for HCI modes, including a review of the final design of individual components such as the vortex phase masks, the grayscale ring apodizer, and the apodizing phase plates, as well as a description of their on-going performance tests and of our plans for system-level integration and tests. Using end-to-end simulations, we predict the performance that will be reached on sky by the METIS HCI modes in presence of various environmental and instrumental disturbances, including non-common path aberrations and water vapor seeing, and discuss our strategy to mitigate these various effects. We finally illustrate with mock observations and data processing that METIS should be capable of directly imaging temperate rocky planets around the nearest stars., Comment: 16 pages, paper presented at SPIE Astronomical Telescopes + Instrumentation 2024

Elemental abundances of Sun-like stars are crucial for understanding the detailed properties of their planets. However, measuring elemental abundances in M stars is challenging due to their faintness and pervasive molecular features in optical spectra. To address this, elemental abundances of Sun-like stars have been proposed to constrain those of M stars by scaling [X/H] with measured [Fe/H]. This study tests the robustness of this practice using M- and GK-dwarf stellar abundances and rigorous statistical methods. We compile elemental abundances for 43 M dwarfs for 10 major rock-forming elements (Fe, C, O, Mg, Si, Al, Ca, Na, Ni, and Ti) from high-resolution near-infrared stellar surveys. We perform bootstrap-based linear regressions on the M dwarfs to determine the trends of [X/H] vs. [Fe/H] and compare them with GK dwarfs. A 2-sample, multivariate Mahalanobis Distance test is applied to assess the significance of differences in [X/H]--[Fe/H] trends for individual elemental pairs between M and GK dwarfs. The null hypothesis of no significant difference in chemical trends between M and GK dwarfs is strongly rejected for all elements except Si, for which rejection is marginal, and Na and Ni, for which results are inconclusive. This suggests that assuming no difference may lead to biased results and inaccurate constraints on rocky planets around M dwarfs. Therefore, it is crucial for both the stellar and exoplanet communities to recognise these differences. To better understand these differences, we advocate for dedicated modelling techniques for M dwarf atmospheres and more homogeneous abundance analyses. Our statistically constrained trends of [X/H]--[Fe/H] for M dwarfs offer a new constraint on estimating M-dwarf elemental abundances given measured [Fe/H], aiding in characterising the properties of M dwarf-hosted rocky worlds., Comment: 8 pages, 10 figures, 1 table in the main body; 6 appendices, with 13 figures and an open-access code. Accepted by A&A

The main challenge of exoplanet high-contrast imaging (HCI) is to separate the signal of exoplanets from their host stars, which are many orders of magnitude brighter. HCI for ground-based observations is further exacerbated by speckle noise originating from perturbations in the Earth's atmosphere and imperfections in the telescope optics. Various data post-processing techniques are used to remove this speckle noise and reveal the faint planet signal. Often, however, a significant part of the planet signal is accidentally subtracted together with the noise. In the present work, we use explainable machine learning to investigate the reason for the loss of the planet signal for one of the most used post-processing methods: Principal Component Analysis (PCA). We find that PCA learns the shape of the telescope point spread function for high numbers of PCA components. This representation of the noise captures not only the speckle noise, but also the characteristic shape of the planet signal. Building upon these insights, we develop a new post-processing method (4S) that constrains the noise model to minimize this signal loss. We apply our model to 11 archival HCI datasets from the VLT-NACO instrument in the L'-band and find that our model consistently outperforms PCA. The improvement is largest at close separations to the star ($\leq 4 \lambda /D$) providing up to 1.5 magnitudes deeper contrast. This enhancement enables us to detect the exoplanet AF Lep b in data from 2011, 11 years before its subsequent discovery. We present updated orbital parameters for this object., Comment: Submitted to AJ, 26 pages, 15 figures, comments welcome, code available on ReadtheDocs: https://fours.readthedocs.io/en/latest/

The new generation of observatories and instruments (VLT/ERIS, JWST, ELT) motivate the development of robust methods to detect and characterise faint and close-in exoplanets. Molecular mapping and cross-correlation for spectroscopy use molecular templates to isolate a planet's spectrum from its host star. However, reliance on signal-to-noise ratio (S/N) metrics can lead to missed discoveries, due to strong assumptions of Gaussian independent and identically distributed noise. We introduce machine learning for cross-correlation spectroscopy (MLCCS); the method aims to leverage weak assumptions on exoplanet characterisation, such as the presence of specific molecules in atmospheres, to improve detection sensitivity for exoplanets. MLCCS methods, including a perceptron and unidimensional convolutional neural networks, operate in the cross-correlated spectral dimension, in which patterns from molecules can be identified. We test on mock datasets of synthetic planets inserted into real noise from SINFONI at K-band. The results from MLCCS show outstanding improvements. The outcome on a grid of faint synthetic gas giants shows that for a false discovery rate up to 5%, a perceptron can detect about 26 times the amount of planets compared to an S/N metric. This factor increases up to 77 times with convolutional neural networks, with a statistical sensitivity shift from 0.7% to 55.5%. In addition, MLCCS methods show a drastic improvement in detection confidence and conspicuity on imaging spectroscopy. Once trained, MLCCS methods offer sensitive and rapid detection of exoplanets and their molecular species in the spectral dimension. They handle systematic noise and challenging seeing conditions, can adapt to many spectroscopic instruments and modes, and are versatile regarding atmospheric characteristics, which can enable identification of various planets in archival and future data., Comment: 27 pages, 24 figures. Submitted for publication in A&A January 2, 2024. After first iteration with the referee, resubmitted May 17, 2024

No matter how fascinating and exotic other terrestrial planets are revealed to be, nothing generates more excitement than announcements regarding their habitability. From the observation of Mars to present-day efforts toward Venus and the characterization of exoplanets, the search for life, or at least environments that could accommodate life, has been a major drive for space exploration. So far, we have found no other unquestionably habitable world besides Earth. The conditions of the habitability of terrestrial planets have proved elusive, as surface conditions depend on the complex interplay of many processes throughout the evolution of a planet. Here, we review how the interior of a rocky planet can drive the evolution of surface conditions and the atmosphere. Instead of listing criteria assumed to be critical for life, we discuss how the bulk-silicate planet can affect the onset, continuation and cessation of habitability. We then consider how it can be observed and current efforts towards this end.

Direct imaging of Earth-like exoplanets is one of the most prominent scientific drivers of the next generation of ground-based telescopes. Typically, Earth-like exoplanets are located at small angular separations from their host stars, making their detection difficult. Consequently, the adaptive optics (AO) system's control algorithm must be carefully designed to distinguish the exoplanet from the residual light produced by the host star. A new promising avenue of research to improve AO control builds on data-driven control methods such as Reinforcement Learning (RL). RL is an active branch of the machine learning research field, where control of a system is learned through interaction with the environment. Thus, RL can be seen as an automated approach to AO control, where its usage is entirely a turnkey operation. In particular, model-based reinforcement learning (MBRL) has been shown to cope with both temporal and misregistration errors. Similarly, it has been demonstrated to adapt to non-linear wavefront sensing while being efficient in training and execution. In this work, we implement and adapt an RL method called Policy Optimization for AO (PO4AO) to the GHOST test bench at ESO headquarters, where we demonstrate a strong performance of the method in a laboratory environment. Our implementation allows the training to be performed parallel to inference, which is crucial for on-sky operation. In particular, we study the predictive and self-calibrating aspects of the method. The new implementation on GHOST running PyTorch introduces only around 700 microseconds in addition to hardware, pipeline, and Python interface latency. We open-source well-documented code for the implementation and specify the requirements for the RTC pipeline. We also discuss the important hyperparameters of the method, the source of the latency, and the possible paths for a lower latency implementation., Comment: Accepted for publication in JATIS

Atmospheric retrievals (AR) characterize exoplanets by estimating atmospheric parameters from observed light spectra, typically by framing the task as a Bayesian inference problem. However, traditional approaches such as nested sampling are computationally expensive, thus sparking an interest in solutions based on machine learning (ML). In this ongoing work, we first explore flow matching posterior estimation (FMPE) as a new ML-based method for AR and find that, in our case, it is more accurate than neural posterior estimation (NPE), but less accurate than nested sampling. We then combine both FMPE and NPE with importance sampling, in which case both methods outperform nested sampling in terms of accuracy and simulation efficiency. Going forward, our analysis suggests that simulation-based inference with likelihood-based importance sampling provides a framework for accurate and efficient AR that may become a valuable tool not only for the analysis of observational data from existing telescopes, but also for the development of new missions and instruments., Comment: Accepted at the "AI to Accelerate Science and Engineering (AI2ASE)" workshop at AAAI 2024

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

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

In this study, we treat Earth as an exoplanet and investigate our home planet by means of a potential future mid-infrared (MIR) space mission called the Large Interferometer For Exoplanets (LIFE). We combine thermal spectra from an empirical dataset of disk-integrated Earth observations with a noise model for LIFE to create mock observations. We apply a state-of-the-art atmospheric retrieval framework to characterize the planet, assess the potential for detecting the known bioindicators, and investigate the impact of viewing geometry, seasonality, and patchy clouds on the characterization. Key findings include: (1) we are observing a temperate habitable planet with significant abundances of $\mathrm{CO_2}$, $\mathrm{H_2O}$, $\mathrm{O_3}$, and $\mathrm{CH_4}$; (2) seasonal variations in the surface and equilibrium temperature, and in the Bond albedo are detectable; (3) the viewing geometry and the spatially and temporally unresolved nature of our observations only have a minor impact on the characterization; (4) Earth's variable H2O profile and patchy cloud coverage lead to biased retrieval results for the atmospheric structure and trace gas abundances; (5) the limited extent of Earth's seasonal variations in biosignature abundances makes the direct detection of its biosphere through atmospheric seasonality unlikely. Our results suggest that LIFE could correctly identify Earth as a planet where life could thrive, with detectable levels of bioindicators, a temperate climate, and surface conditions allowing liquid surface water. Even if atmospheric seasonality is not easily observed, our study demonstrates that next generation, optimized space missions can assess whether nearby temperate terrestrial exoplanets are habitable or even inhabited., Comment: Published in ApJ. Update: Correction of typo in title. 15 pages (main text incl. 7 figures and 3 tables) + appendix. This paper is part of a series: arXiv:2010.02589, arXiv:2210.05414

The investigation of the atmospheres of closely separated, directly imaged gas giant exoplanets is challenging due to the presence of stellar speckles that pollute their spectrum. To remedy this, the analysis of medium- to high-resolution spectroscopic data via cross-correlation with spectral templates (cross-correlation spectroscopy) is emerging as a leading technique. We aim to define a robust Bayesian framework combining, for the first time, three widespread direct-imaging techniques, namely photometry, low-resolution spectroscopy, and medium-resolution cross-correlation spectroscopy in order to derive the atmospheric properties of close-in directly imaged exoplanets. Our framework CROCODILE (cross-correlation retrievals of directly imaged self-luminous exoplanets) naturally combines the three techniques by adopting adequate likelihood functions. To validate our routine, we simulated observations of gas giants similar to the well-studied $\beta$~Pictoris~b planet and we explored the parameter space of their atmospheres to search for potential biases. We obtain more accurate measurements of atmospheric properties when combining photometry, low- and medium-resolution spectroscopy into atmospheric retrievals than when using the techniques separately as is usually done in the literature. We find that medium-resolution ($R \approx 4000$) K-band cross-correlation spectroscopy alone is not suitable to constrain the atmospheric properties of our synthetic datasets; however, this problem disappears when simultaneously fitting photometry and low-resolution ($R \approx 60$) spectroscopy between the Y and M bands. Our framework allows the atmospheric characterisation of directly imaged exoplanets using the high-quality spectral data that will be provided by the new generation of instruments such as VLT/ERIS, JWST/MIRI, and ELT/METIS.

We present model constraints on the atmospheric structure of HD 106906 b, a planetary-mass companion orbiting at a ~700 AU projected separation around a 15 Myr-old stellar binary, using the APOLLO retrieval code on spectral data spanning 1.1-2.5 $\mu$m. C/O ratios can provide evidence for companion formation pathways, as such pathways are ambiguous both at wide separations and at star-to-companion mass ratios in the overlap between the distributions of planets and brown dwarfs. We benchmark our code against an existing retrieval of the field L dwarf 2M2224-0158, returning a C/O ratio consistent with previous fits to the same JHKs data, but disagreeing in the thermal structure, cloud properties, and atmospheric scale height. For HD 106906 b, we retrieve C/O $=0.53^{+0.15}_{-0.25}$, consistent with the C/O ratios expected for HD 106906's stellar association and therefore consistent with a stellar-like formation for the companion. We find abundances of H$_2$O and CO near chemical equilibrium values for a solar metallicity, but a surface gravity lower than expected, as well as a thermal profile with sharp transitions in the temperature gradient. Despite high signal-to-noise and spectral resolution, more accurate constraints necessitate data across a broader wavelength range. This work serves as preparation for subsequent retrievals in the era of JWST, as JWST's spectral range provides a promising opportunity to resolve difficulties in fitting low-gravity L dwarfs, and also underscores the need for simultaneous comparative retrievals on L dwarf companions with multiple retrieval codes., Comment: 39 pages, 21 figures. Accepted for publication in the Astronomical Journal

Atmospheric retrievals (AR) of exoplanets typically rely on a combination of a Bayesian inference technique and a forward simulator to estimate atmospheric properties from an observed spectrum. A key component in simulating spectra is the pressure-temperature (PT) profile, which describes the thermal structure of the atmosphere. Current AR pipelines commonly use ad hoc fitting functions here that limit the retrieved PT profiles to simple approximations, but still use a relatively large number of parameters. In this work, we introduce a conceptually new, data-driven parameterization scheme for physically consistent PT profiles that does not require explicit assumptions about the functional form of the PT profiles and uses fewer parameters than existing methods. Our approach consists of a latent variable model (based on a neural network) that learns a distribution over functions (PT profiles). Each profile is represented by a low-dimensional vector that can be used to condition a decoder network that maps $P$ to $T$. When training and evaluating our method on two publicly available datasets of self-consistent PT profiles, we find that our method achieves, on average, better fit quality than existing baseline methods, despite using fewer parameters. In an AR based on existing literature, our model (using two parameters) produces a tighter, more accurate posterior for the PT profile than the five-parameter polynomial baseline, while also speeding up the retrieval by more than a factor of three. By providing parametric access to physically consistent PT profiles, and by reducing the number of parameters required to describe a PT profile (thereby reducing computational cost or freeing resources for additional parameters of interest), our method can help improve AR and thus our understanding of exoplanet atmospheres and their habitability., Comment: Accepted for publication in Astronomy & Astrophysics

The next generation of space-based observatories will characterize the atmospheres of low-mass, temperate exoplanets with the direct-imaging technique. This will be a major step forward in our understanding of exoplanet diversity and the prevalence of potentially habitable conditions beyond the Earth. We compute a list of currently known exoplanets detectable with the mid-infrared Large Interferometer For Exoplanets (LIFE) in thermal emission. We also compute the list of known exoplanets accessible to a notional design of the Habitable Worlds Observatory (HWO), observing in reflected starlight. With a pre-existing method, we processed the NASA Exoplanet Archive and computed orbital realizations for each known exoplanet. We derived their mass, radius, equilibrium temperature, and planet-star angular separation. We used the LIFEsim simulator to compute the integration time ($t_{int}$) required to detect each planet with LIFE. A planet is considered detectable if a broadband signal-to-noise ratio $S/N$=7 is achieved over the spectral range $4-18.5\mu$m in $t_{int}\leq$100 hours. We tested whether the planet is accessible to HWO in reflected starlight based on its notional inner and outer working angles, and minimum planet-to-star contrast. LIFE's reference configuration (four 2-m telescopes with 5% throughput and a nulling baseline between 10-100 m) can detect 212 known planets within 20 pc. Of these, 55 are also accessible to HWO in reflected starlight, offering a unique opportunity for synergies in atmospheric characterization. LIFE can also detect 32 known transiting exoplanets. Furthermore, 38 LIFE-detectable planets orbit in the habitable zone, of which 13 with $M_p<5M_\oplus$ and 8 with $5M_\oplus

A mid-infrared nulling-space interferometer is a promising way to characterize thermal light from habitable planet candidates around Sun-like stars. However, one of the main challenges for achieving this ambitious goal is a high-precision stability of the optical path difference (OPD) and amplitude over a few days for planet detection and up to a few weeks for in-depth characterization. Here we propose a new method called phase-space synthesis decomposition (PSSD) to shorten the stability requirement to minutes, significantly relaxing the technological challenges of the mission. Focusing on what exactly modulates the planet signal in the presence of the stellar leak and systematic error, PSSD prioritizes the modulation of the signals along the wavelength domain rather than baseline rotation. Modulation along the wavelength domain allows us to extract source positions in parallel to the baseline vector for each exposure. The sum of the one-dimensional data converts into two-dimensional information. Based on the reconstructed image, we construct a continuous equation and extract the spectra through the singular value decomposition (SVD) while efficiently separating them from a long-term systematic stellar leak. We performed numerical simulations to investigate the feasibility of PSSD for the LIFE mission concept. We confirm that multiple terrestrial planets in the habitable zone around a Sun-like star at 10 pc can be detected and characterized despite high levels and long durations of systematic noise. We also find that PSSD is more robust against a sparse sampling of the array rotation compared to purely rotation-based signal extraction. Using PSSD as signal extraction method significantly relaxes the technical requirements on signal stability and further increases the feasibility of the LIFE mission., Comment: Accepted for publication by A&A - 14 pages main text including 10 Figures

Over the past decade, hundreds of nights have been spent on the worlds largest telescopes to search for and directly detect new exoplanets using high-contrast imaging (HCI). Thereby, two scientific goals are of central interest: First, to study the characteristics of the underlying planet population and distinguish between different planet formation and evolution theories. Second, to find and characterize planets in our immediate Solar neighborhood. Both goals heavily rely on the metric used to quantify planet detections and non-detections. Current standards often rely on several explicit or implicit assumptions about the noise. For example, it is often assumed that the residual noise after data post-processing is Gaussian. While being an inseparable part of the metric, these assumptions are rarely verified. This is problematic as any violation of these assumptions can lead to systematic biases. This makes it hard, if not impossible, to compare results across datasets or instruments with different noise characteristics. We revisit the fundamental question of how to quantify detection limits in HCI. We focus our analysis on the error budget resulting from violated assumptions. To this end, we propose a new metric based on bootstrapping that generalizes current standards to non-Gaussian noise. We apply our method to archival HCI data from the NACO-VLT instrument and derive detection limits for different types of noise. Our analysis shows that current standards tend to give detection limit that are about one magnitude too optimistic in the speckle-dominated regime. That is, HCI surveys may have excluded planets that can still exist., Comment: After first iteration with the referee, resubmitted to AJ. Comments welcome!

ERIS (Enhanced Resolution Imager and Spectrograph) is a new adaptive optics instrument installed at the Cassegrain focus of the VLT-UT4 telescope at the Paranal Observatory in Chile. ERIS consists of two near-infrared instruments: SPIFFIER, an integral field unit (IFU) spectrograph covering J to K bands, and NIX, an imager covering J to M bands. ERIS has an adaptive optics system able to work with both LGS and NGS. The Assembly Integration Verification (AIV) phase of ERIS at the Paranal Observatory was carried out starting in December 2021, followed by several commissioning runs in 2022. This contribution will describe the first preliminary results of the on-sky performance of ERIS during its commissioning and the future perspectives based on the preliminary scientific results., Comment: Proceeding of SPIE Astronomical Telescopes + Instrumentation 2022

We present the statistical analysis of a subsample of 45 young stars surrounded by protoplanetary disks (PPDs). This is the largest imaging survey uniquely focused on PPDs to date. Our goal is to search for young forming companions embedded in the disk material and to constrain their occurrence rate in relation to the formation mechanism. We used principal component analysis based point spread function subtraction techniques to reveal young companions forming in the disks. We calculated detection limits for our datasets and adopted a black-body model to derive temperature upper limits of potential forming planets. We then used Monte Carlo simulations to constrain the population of forming gas giant companions and compare our results to different types of formation scenarios. Our data revealed a new binary system (HD38120) and a recently identified triple system with a brown dwarf companion orbiting a binary system (HD101412), in addition to 12 known companions. Furthermore, we detected signals from 17 disks, two of which (HD72106 and TCrA) were imaged for the first time. We reached median detection limits of L =15.4 mag at 2.0 arcsec, which were used to investigate the temperature of potentially embedded forming companions. We can constrain the occurrence of forming planets with semi-major axis a in [20 - 500] au and Teff in [600 - 3000] K, in line with the statistical results obtained for more evolved systems from other direct imaging surveys. The NaCo-ISPY data confirm that massive bright planets accreting at high rates are rare. More powerful instruments with better sensitivity in the near- to mid-infrared are likely required to unveil the wealth of forming planets sculpting the observed disk substructures., Comment: 25 pages, 16 figures, 3 tables, accepted for publication in A&A

Context. Planets around stars in the solar neighbourhood will be prime targets for characterisation with upcoming large space- and ground-based facilities. Since large-scale exoplanet searches will not be feasible with such telescopes, it is crucial to use currently available data and instruments to find possible target planets before next generation facilities come online. Aims. We aim at detecting new extrasolar planets around stars in the solar neighbourhood by blind radial velocity (RV) search with ESPRESSO. Our target sample consist of nearby stars (d < 11 pc) with little (< 10) or no previous RV measurements. Methods. We use 31 radial velocity measurements obtained with ESPRESSO at the VLT between December 2020 and February 2022 of the nearby M dwarf star (M_star = 0.21 M_sun, d = 10.23 pc) L 363-38 to derive the orbital parameters of the newly discovered planet. In addition, we use TESS photometry and archival VLT/NaCo high contrast imaging data to put further constraints on the orbit inclination and the possible planetary system architecture around L 363-38. Results. We present the detection of a new extrasolar planet orbiting the nearby M dwarf star L 363-38. L 363-38 b is a planet with minimum mass mp sin(i) = 4.67+/-0.43 M_Earth orbiting its star with a period P = 8.781+/-0.007 d, corresponding to a semi-major axis a = 0.048+/-0.006 AU, which is well inside the inner edge of the habitable zone. We further estimate a minimum radius rp sin(i) = 1.55 - 2.75 R_Earth and an equilibrium temperature Teq = 330K., Comment: 12 pages, 11 figures

We assess the dependence of Earth's disk-integrated mid-infrared thermal emission spectrum on observation geometries and investigate which and how spectral features are impacted by seasonality on Earth. We compiled an exclusive dataset containing 2690 disk-integrated thermal emission spectra for four different full-disk observing geometries (North & South Pole centered and Africa & Pacific centred equatorial views) over four consecutive years. The spectra were derived from 2378 spectral channels in the wavelength range from 3.75 to 15.4 micron (nominal resolution $\approx$ 1200) and were recorded by the Atmospheric Infrared Sounder aboard the Aqua satellite. We learned that there is significant seasonal variability in Earth's thermal emission spectrum, and the strength of spectral features of bio-indicators, such as N2O, CH4, O3 and CO2 depends strongly on both season and viewing geometry. In addition, we found a strong spectral degeneracy with respect to the latter two indicating that multi-epoch measurements and time-dependent signals may be required in order to fully characterize planetary environments. Even for Earth and especially for equatorial views, the variations in flux and strength of absorption features in the disk-integrated data are small and typically $\le$ 10%. Disentangling these variations from the noise in future exoplanet observations will be a challenge. However, irrespectively of when the planet will be measured (i.e., day or night or season) the results from mid-infrared observations will remain the same to the zeroth order which is an advantage over reflected light observations., Comment: 21 pages, 15 Figures, 3 Tables

While previous studies have shown a strong preference for a future mid-infrared nulling interferometer space mission to detect planets within the HZ around M dwarfs, we here focus on a more conservative approach toward the concept of habitability and present yield estimates for two stellar samples consisting of nearby (d<20 pc) Sun-like stars (4800-6300 K) and nearby FGK-type stars (3940-7220 K) accessible to such a mission. Our yield estimates are based on recently derived occurrence rates of rocky planets from the Kepler mission and our LIFE exoplanet observation simulation tool LIFEsim, which includes all main astrophysical noise sources, but no instrumental noise sources as yet. Depending on a pessimistic or optimistic extrapolation of the Kepler results, we find that during a 2.5-year search phase, LIFE could detect between ~10-16 (average) or ~5-34 (including 1$\sigma$ uncertainties) rocky planets (0.5-1.5 R${}_\rm{Earth}$) within the optimistic HZ of Sun-like stars and between ~4-6 (average) or ~1-13 (including 1$\sigma$ uncertainties) exo-Earth candidates (EECs) assuming four collector spacecraft equipped with 2 m mirrors and a conservative instrument throughput of 5%. With D=3.5 m or 1 m mirrors, the yield $Y$ changes strongly, following approximately $Y \propto D^{3/2}$. With the larger sample of FGK-type stars, the yield increases to ~16-22 (average) rocky planets within the optimistic HZ and ~5-8 (average) EECs. Furthermore, we find that in addition to the mirror diameter, the yield depends strongly on the total throughput, but only weakly on the exozodiacal dust level and the accessible wavelength range of the mission. When the focus lies entirely on Sun-like stars, larger mirrors (~3 m with 5% total throughput) or a better total throughput (~20% with 2 m mirrors) are required to detect a statistically relevant sample of ~30 rocky planets within the optimistic HZ., Comment: 19 pages, 13 figures, 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

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

Gravitational interactions in star-forming regions are capable of disrupting and destroying planetary systems, as well as creating new ones. In particular, a planet can be stolen, where it is directly exchanged between passing stars during an interaction; or captured, where a planet is first ejected from its birth system and is free-floating for a period of time, before being captured by a passing star. We perform sets of direct N-body simulations of young, substructured star-forming regions, and follow their evolution for 10 Myr in order to determine how many planets are stolen and captured, and their respective orbital properties. We show that in high density star-forming regions, stolen and captured planets have distinct properties. The semimajor axis distribution of captured planets is significantly skewed to wider orbits compared to the semimajor axis distribution of stolen planets and planets that are still orbiting their parent star (preserved planets). However, the eccentricity and inclination distributions of captured and stolen planets are similar, but in turn very different to the inclination and eccentricity distributions of preserved planets. In low-density star-forming regions these differences are not as distinct but could still, in principle, be used to determine whether observed exoplanets have likely formed in situ or have been stolen or captured. We find that the initial degree of spatial and kinematic substructure in a star-forming region is as important a factor as the stellar density in determining whether a planetary system will be altered, disrupted, captured or stolen., Comment: 15 pages, 10 figures, 1 appendix, accepted for publication in MNRAS

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

A major goal in the discovery and characterisation of exoplanets is to identify terrestrial-type worlds that are similar to (or otherwise distinct from) our Earth. Recent results have highlighted the importance of applying devolatilisation -- i.e. depletion of volatiles -- to the chemical composition of planet-hosting stars to constrain bulk composition and interiors of terrestrial-type exoplanets. In this work, we apply such an approach to a selected sample of 13 planet-hosting Sun-like stars, for which high-precision photospheric abundances have been determined in the first paper of the series. With the resultant devolatilised stellar composition (i.e. the model planetary bulk composition) as well as other constraints including mass and radius, we model the detailed mineralogy and interior structure of hypothetical, habitable-zone terrestrial planets ("exo-Earths") around these stars. Model output shows that most of these exo-Earths are expected to have broadly Earth-like composition and interior structure, consistent with conclusions derived independently from analysis of polluted white dwarfs. The exceptions are the Kepler-10 and Kepler-37 exo-Earths, which we predict are strongly oxidised and thus would develop metallic cores much smaller than Earth. Investigating our devolatilisation model at its extremes as well as varying planetary mass and radius (within the terrestrial regime) reveals potential diversities in the interiors of terrestrial planets. By considering (i) high-precision stellar abundances, (ii) devolatilisation, and (iii) planetary mass and radius holistically, this work represents essential steps to explore the detailed mineralogy and interior structure of terrestrial-type exoplanets, which in turn are fundamental for our understanding of planetary dynamics and long-term evolution., Comment: Accepted for publication in MNRAS; 19 pages, 14 figures, 2 tables; a related talk is available at https://www.youtube.com/watch?v=n0TPQ2IrAyA

High-contrast imaging of exoplanets hinges on powerful post-processing methods to denoise the data and separate the signal of a companion from its host star, which is typically orders of magnitude brighter. Existing post-processing algorithms do not use all prior domain knowledge that is available about the problem. We propose a new method that builds on our understanding of the systematic noise and the causal structure of the data-generating process. Our algorithm is based on a modified version of half-sibling regression (HSR), a flexible denoising framework that combines ideas from the fields of machine learning and causality. We adapt the method to address the specific requirements of high-contrast exoplanet imaging data obtained in pupil tracking mode. The key idea is to estimate the systematic noise in a pixel by regressing the time series of this pixel onto a set of causally independent, signal-free predictor pixels. We use regularized linear models in this work; however, other (non-linear) models are also possible. In a second step, we demonstrate how the HSR framework allows us to incorporate observing conditions such as wind speed or air temperature as additional predictors. When we apply our method to four data sets from the VLT/NACO instrument, our algorithm provides a better false-positive fraction than PCA-based PSF subtraction, a popular baseline method in the field. Additionally, we find that the HSR-based method provides direct and accurate estimates for the contrast of the exoplanets without the need to insert artificial companions for calibration in the data sets. Finally, we present first evidence that using the observing conditions as additional predictors can improve the results. Our HSR-based method provides an alternative, flexible and promising approach to the challenge of modeling and subtracting the stellar PSF and systematic noise in exoplanet imaging data., Comment: Accepted for publication in Astronomy & Astrophysics

The Large Interferometer For Exoplanets (LIFE) initiative is developing the science and a technology roadmap for an ambitious space mission featuring a space-based mid-infrared (MIR) nulling interferometer in order to detect the thermal emission of hundreds of exoplanets and characterize their atmospheres. In order to quantify the science potential of such a mission, in particular in the context of technical trade-offs, an instrument simulator is required. In addition, signal extraction algorithms are needed to verify that exoplanet properties (e.g., angular separation, spectral flux) contained in simulated exoplanet datasets can be accurately retrieved. We present LIFEsim, a software tool developed for simulating observations of exoplanetary systems with an MIR space-based nulling interferometer. It includes astrophysical noise sources (i.e., stellar leakage and thermal emission from local zodiacal and exo-zodiacal dust) and offers the flexibility to include instrumental noise terms in the future. LIFEsim provides an accessible way for predicting the expected SNR of future observations as a function of various key instrument and target parameters. The SNRs of the extracted spectra are photon-noise dominated, as expected from our current simulations. From single epoch observations in our mock survey of small ($R < 1.5 R_\mathrm{Earth}$) planets orbiting within the habitable zones of their stars, we find that typical uncertainties in the estimated effective temperature of the exoplanets are $\lesssim$10%, for the exoplanet radius $\lesssim$20%, and for the separation from the host star $\lesssim$2%. SNR values obtained in the signal extraction process deviate less than 10% from purely photon-counting statistics based SNRs. (abridged), Comment: Accepted for publication in A&A; 15 pages (main text incl. 13 figures and 1 table) + appendix; comments are welcome

The Large Binocular Telescope (LBT) has two 8.4-m primary mirrors that produce beams that can be combined coherently in a "Fizeau" interferometric mode. In principle, the Fizeau PSF enables the probing of structure at a resolution up to three times better than that of the adaptive-optics-corrected PSF of a single 8.4-m telescope. In this work, we examined the nearby star Altair (5.13 pc, type A7V, $\sim$100s Myr to $\approx$1.4 Gyr) in the Fizeau mode with the LBT at Br-$\alpha$ (4.05 $\mu$m) and carried out angular differential imaging to search for companions. This work presents the first filled-aperture LBT Fizeau science dataset to benefit from a correcting mirror which provides active phase control. In the analysis of the $\lambda/D$ angular regime, the sensitivity of the dataset is down to $\approx$0.5 $M_{\odot}$ at 1" for a 1.0 Gyr system. This sensitivity remains limited by the small amount of integration time, which is in turn limited by the instability of the Fizeau PSF. However, in the Fizeau fringe regime we attain sensitivities of $\Delta m \approx 5$ at 0.2" and put constraints to companions of 1.3 $M_{\odot}$ down to an inner angle of $\approx$0.15", closer than any previously published direct imaging of Altair. This analysis is a pathfinder for future datasets of this type, and represents some of the first steps to unlocking the potential of the first ELT. Fizeau observations will be able to reach dimmer targets with upgrades to the instrument, in particular the phase detector., Comment: Published in AJ

Planet formation occurs around a wide range of stellar masses and stellar system architectures. An improved understanding of the formation process can be achieved by studying it across the full parameter space, particularly toward the extremes. Earlier studies of planets in close-in orbits around high-mass stars have revealed an increase in giant planet frequency with increasing stellar mass until a turnover point at 1.9 solar masses, above which the frequency rapidly decreases. This could potentially imply that planet formation is impeded around more massive stars, and that giant planets around stars exceeding 3 solar masses may be rare or non-existent. However, the methods used to detect planets in small orbits are insensitive to planets in wide orbits. Here we demonstrate the existence of a planet at 560 times the Sun-Earth distance from the 6-10 solar mass binary b Centauri through direct imaging. The planet-to-star mass ratio of 0.10-0.17% is similar to the Jupiter-Sun ratio, but the separation of the detected planet is ~100 times wider than that of Jupiter. Our results show that planets can reside in much more massive stellar systems than what would be expected from extrapolation of previous results. The planet is unlikely to have formed in-situ through the conventional core accretion mechanism, but might have formed elsewhere and arrived to its present location through dynamical interactions, or might have formed via gravitational instability., Comment: Manuscript version. Published in Nature 9 December 2021

The bulk chemical composition and interior structure of rocky exoplanets are of fundamental importance to understanding their long-term evolution and potential habitability. Observations of the chemical compositions of the solar system rocky bodies and of other planetary systems have increasingly shown a concordant picture that the chemical composition of rocky planets reflects that of their host stars for refractory elements, whereas this expression breaks down for volatiles. This behavior is explained by devolatilization during planetary formation and early evolution. Here, we apply a devolatilization model calibrated with solar system bodies to the chemical composition of our nearest Sun-like stars -- $\alpha$ Centauri A and B -- to estimate the bulk composition of any habitable-zone rocky planet in this binary system ("$\alpha$-Cen-Earth"). Through further modeling of likely planetary interiors and early atmospheres, we find that compared to Earth, such a planet is expected to have (i) a reduced (primitive) mantle that is similarly dominated by silicates albeit enriched in carbon-bearing species (graphite/diamond); (ii) a slightly larger iron core, with a core mass fraction of $38.4_{-5.1}^{+4.7}$ wt% (cf. Earth's 32.5 $\pm$ 0.3 wt%); (iii) an equivalent water-storage capacity; and (iv) a CO$_2$-CH$_4$-H$_2$O-dominated early atmosphere that resembles that of Archean Earth. Further taking into account its $\sim$ 25% lower intrinsic radiogenic heating from long-lived radionuclides, an ancient $\alpha$-Cen-Earth ($\sim$ 1.5-2.5 Gyr older than Earth) is expected to have less efficient mantle convection and planetary resurfacing, with a potentially prolonged history of stagnant-lid regimes., Comment: Accepted to ApJ; 27 pages, 9 figures, 5 tables

GQ Lup B is a young and accreting, substellar companion that appears to drive a spiral arm in the circumstellar disk of its host star. We report high-contrast imaging observations of GQ Lup B with VLT/NACO at 4-5 $\mu$m and medium-resolution integral field spectroscopy with VLT/MUSE. The optical spectrum is consistent with an M9 spectral type, shows characteristics of a low-gravity atmosphere, and exhibits strong H$\alpha$ emission. The $H-M'$ color is $\gtrsim$1 mag redder than field dwarfs with similar spectral types and a detailed analysis of the spectral energy distribution (SED) from optical to mid-infrared wavelengths reveals excess emission in the $L'$, NB4.05, and $M'$ bands. The excess flux is well described by a blackbody component with $T_\mathrm{disk} \approx 460$ K and $R_\mathrm{disk} \approx 65\,R_\mathrm{J}$ and is expected to trace continuum emission from small grains in a protolunar disk. We derive an extinction of $A_V \approx 2.3$ mag from the broadband SED with a suspected origin in the vicinity of the companion. We also combine 15 yr of astrometric measurements and constrain the mutual inclination with the circumstellar disk to $84 \pm 9$ deg, indicating a tumultuous dynamical evolution or a stellar-like formation pathway. From the measured H$\alpha$ flux and the estimated companion mass, $M_\mathrm{p} \approx 30\,M_\mathrm{J}$, we derive an accretion rate of $\dot{M} \approx 10^{-6.5}\,M_\mathrm{J}\,\mathrm{yr}^{-1}$. We speculate that the disk is in a transitional stage in which the assembly of satellites from a pebble reservoir has opened a central cavity while GQ Lup B is in the final stages of its formation., Comment: 28 pages, 11 figures, accepted for publication in AJ

Direct detection and characterization of Earth-like planets around Sun-like stars is a core task for evaluating the prevalence of habitability and life in the Universe. Here, we discuss a promising option for achieving this goal, which is based on placing an occulter in orbit and having it project its shadow onto the E- ELT at the surface of Earth, thus providing a sufficient contrast for imaging and taking spectra of Earth-like planets in the habitable zones of Sun-like stars. Doing so at a sensible fuel budget will require tailored orbits, an occulter with a high area-to-mass ratio, and appropriate instrumentation at the E-ELT. In this White Paper, submitted in response to the ESA Voyage 2050 Call, we outline the fundamental aspects of the concept, and the most important technical developments that will be required to develop a full mission., Comment: Published in special issue of Experimental Astronomy together with other White Papers submitted to ESA in 2019 in the context of the Voyage 2050 long-term planning of the Science Programme. This is the manuscript version; the published journal version is available with open access at: http://link.springer.com/article/10.1007/s10686-021-09792-y

The past decade has seen increasing efforts in detecting and characterising exoplanets by high contrast imaging in the near/mid-infrared, which is the optimal wavelength domain for studying old, cold planets. In this work, we present deep AO imaging observations of the nearby Sun-like star $\epsilon$ Ind A with NaCo ($L^{\prime}$) and NEAR (10-12.5 microns) instruments at VLT, in an attempt to directly detect its planetary companion whose presence has been indicated from radial velocity (RV) and astrometric trends. We derive brightness limits from the non-detection of the companion with both instruments, and interpret the corresponding sensitivity in mass based on both cloudy and cloud-free atmospheric and evolutionary models. For an assumed age of 5 Gyr for the system, we get detectable mass limits as low as 4.4 $M_{\rm J}$ in NaCo $L^{\prime}$ and 8.2 $M_{\rm J}$ in NEAR bands at 1.5$\arcsec$ from the central star. If the age assumed is 1 Gyr, we reach even lower mass limits of 1.7 $M_{\rm J}$ in NaCo $L^{\prime}$ and 3.5 $M_{\rm J}$ in NEAR bands, at the same separation. However, based on the dynamical mass estimate (3.25 $M_{\rm J}$) and ephemerides from astrometry and RV, we find that the non-detection of the planet in these observations puts a constraint of 2 Gyr on the lower age limit of the system. NaCo offers the highest sensitivity to the planetary companion in these observations, but the combination with the NEAR wavelength range adds a considerable degree of robustness against uncertainties in the atmospheric models. This underlines the benefits of including a broad set of wavelengths for detection and characterisation of exoplanets in direct imaging studies., Comment: 13 pages, 9 figures, 2 tables Accepted for publication in section 10. Planets and planetary systems of Astronomy and Astrophysics

The Mid-infrared ELT Imager and Spectrograph (METIS) will provide the Extremely Large Telescope (ELT) with a unique window to the thermal- and mid-infrared (3 - 13 microns). Its single-conjugate adaptive optics (SCAO) system will enable high contrast imaging and integral field unit (IFU) spectroscopy (R~100,000) at the diffraction limit of the ELT. This article describes the science drivers, conceptual design, observing modes, and expected performance of METIS., Comment: 5 pages, 3 figures, published in the ESO Messenger

The Planetary Camera and Spectrograph (PCS) for the Extremely Large Telescope (ELT) will be dedicated to detecting and characterising nearby exoplanets with sizes from sub-Neptune to Earth-size in the neighbourhood of the Sun. This goal is achieved by a combination of eXtreme Adaptive Optics (XAO), coronagraphy and spectroscopy. PCS will allow us not only to take images, but also to look for biosignatures such as molecular oxygen in the exoplanets' atmospheres. This article describes the PCS primary science goals, the instrument concept and the research and development activities that will be carried out over the coming years., Comment: 6 pages, 6 figures

While the occurrence rate of wide giant planets appears to increase with stellar mass at least up through the A-type regime, B-type stars have not been systematically studied in large-scale surveys so far. It therefore remains unclear up to what stellar mass this occurrence trend continues. The B-star Exoplanet Abundance Study (BEAST) is a direct imaging survey with the extreme adaptive optics instrument SPHERE, targeting 85 B-type stars in the young Scorpius-Centaurus (Sco-Cen) region with the aim to detect giant planets at wide separations and constrain their occurrence rate and physical properties. The statistical outcome of the survey will help determine if and where an upper stellar mass limit for planet formation occurs. In this work, we describe the selection and characterization of the BEAST target sample. Particular emphasis is placed on the age of each system, which is a central parameter in interpreting direct imaging observations. We implement a novel scheme for age dating based on kinematic sub-structures within Sco-Cen, which complements and expands upon previous age determinations in the literature. We also present initial results from the first epoch observations, including the detections of ten stellar companions, of which six were previously unknown. All planetary candidates in the survey will need follow up in second epoch observations, which are part of the allocated observational programme and will be executed in the near future., Comment: 22 pages, 15 figures, accepted for publication in A&A