Search

Effective image post-processing algorithms are vital for the successful direct imaging of exoplanets. Standard PSF subtraction methods use techniques based on a low-rank approximation to separate the rotating planet signal from the quasi-static speckles, and rely on signal-to-noise ratio maps to detect the planet. These steps do not interact or feed each other, leading to potential limitations in the accuracy and efficiency of exoplanet detection. We aim to develop a novel approach that iteratively finds the flux of the planet and the low-rank approximation of quasi-static signals, in an attempt to improve upon current PSF subtraction techniques. In this study, we extend the standard L2 norm minimization paradigm to an L1 norm minimization framework to better account for noise statistics in the high contrast images. Then, we propose a new method, referred to as Alternating Minimization Algorithm with Trajectory, that makes a more advanced use of estimating the low-rank approximation of the speckle field and the planet flux by alternating between them and utilizing both L1 and L2 norms. For the L1 norm minimization, we propose using L1 norm low-rank approximation, a low-rank approximation computed using an exact block-cyclic coordinate descent method, while we use randomized singular value decomposition for the L2 norm minimization. Additionally, we enhance the visibility of the planet signal using a likelihood ratio as a postprocessing step. Numerical experiments performed on a VLT/SPHERE-IRDIS dataset show the potential of AMAT to improve upon the existing approaches in terms of higher S/N, sensitivity limits, and ROC curves. Moreover, for a systematic comparison, we used datasets from the exoplanet data challenge to compare our algorithm to other algorithms in the challenge, and AMAT with likelihood ratio map performs better than most algorithms tested on the exoplanet data challenge., Comment: Submitted to A&A, 20 pages, 20 figures

In this communication, we report on the results of the second phase of the Exoplanet Imaging Data Challenge started in 2019. This second phase focuses on the characterization of point sources (exoplanet signals) within multispectral high-contrast images from ground-based telescopes. We collected eight data sets from two high-contrast integral field spectrographs (namely Gemini-S/GPI and VLT/SPHERE-IFS) that we calibrated homogeneously, and in which we injected a handful of synthetic planetary signals (ground truth) to be characterized by the data challenge participants. The tasks of the participants consist of (1) extracting the precise astrometry of each injected planetary signals, and (2) extracting the precise spectro-photometry of each injected planetary signal. Additionally, the participants may provide the 1-sigma uncertainties on their estimation for further analyses. When available, the participants can also provide the posterior distribution used to estimate the position/spectrum and uncertainties. The data are permanently available on a Zenodo repository and the participants can submit their results through the EvalAI platform. The EvalAI submission platform opened on April 2022 and closed on the 31st of May 2024. In total, we received 4 valid submissions for the astrometry estimation and 4 valid submissions for the spectrophotometry (each submission, corresponding to one pipeline, has been submitted by a unique participant). In this communication, we present an analysis and interpretation of the results., Comment: Submitted to SPIE Astronomical Telescopes + Instrumentation 2024, Adaptive Optics Systems IX, Paper 13097-40

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

Mid-infrared imaging of exoplanets and disks is now possible with the coronagraphs of the MIRI on the JWST. This wavelength range unveils new features of young directly imaged systems and allows us to obtain new constraints for characterizing the atmosphere of young giant exoplanets and associated disks. These observations aim to characterize the atmosphere of the planet HD 95086 b by adding mid-infrared information so that the various hypotheses about its atmospheric parameters values can be unraveled. Improved images of circumstellar disks are provided. We present the MIRI coronagraphic imaging of the system HD 95086 obtained with the F1065C, F1140, and F2300C filters at central wavelengths of 10.575, 11.3, and 23 microns, respectively. We explored the method for subtracting the stellar diffraction pattern in the particular case when bright dust emitting at short separation is present. Furthermore, we compared different methods for extracting the photometry of the planet. Using the atmospheric models Exo-REM and ATMO, we measured the atmospheric parameters of HD 95086 b. The planet HD 95086 b and the contribution from the inner disk are detected at the two shortest MIRI wavelengths F1065C and F1140C. The outer colder belt is imaged at 23 microns. The mid-infrared photometry provides better constraints on the atmospheric parameters. We evaluate a temperature of 800-1050 K, consistent with one previous hypothesis that only used NIR data. The radius measurement of 1.0-1.14 RJup is better aligned with evolutionary models, but still smaller than predicted. These observations allow us to refute the hypothesis of a warm circumplanetary disk. HD 95086 is one of the first exoplanetary systems to be revealed at mid-infrared wavelengths. This highlights the interests and challenges of observations at these wavelengths., Comment: version published in A&A

We present Super-RDI, a unique framework for the application of reference star differential imaging (RDI) to Keck/NIRC2 high-contrast imaging observations with the vortex coronagraph. Super-RDI combines frame selection and signal-to-noise ratio (S/N) optimization techniques with a large multi-year reference point spread function (PSF) library to achieve optimal PSF subtraction at small angular separations. We compile a $\sim$7000 frame reference PSF library based on a set of 288 new Keck/NIRC2 $L'$ sequences of 237 unique targets acquired between 2015 and 2019 as part of two planet-search programs, one focusing on nearby young M dwarfs and the other targeting members of the Taurus star-forming region. For our dataset, synthetic companion injection-recovery tests reveal that frame selection with the mean-squared error (MSE) metric combined with KLIP-based PSF subtraction using 1000-3000 frames and $<$500 principal components yields the highest average S/N for injected synthetic companions. We uniformly reduce targets in the young M-star survey with both Super-RDI and angular differential imaging (ADI). For the typical parallactic angle rotation of our dataset ($\sim$10$^\circ$), Super-RDI performs better than a widely used implementation of ADI at separations $\lesssim$0.4" ($\approx$5 $\lambda$/$D$) gaining an average of 0.25 mag in contrast at 0.25" and 0.4 mag in contrast at 0.15". This represents a performance improvement in separation space over RDI with single-night reference star observations ($\sim$100 frame PSF libraries) applied to a similar Keck/NIRC2 dataset in previous work. We recover two known brown dwarf companions and provide detection limits for 155 targets in the young M-star survey. Our results demonstrate that increasing the PSF library size with careful selection of reference frames can improve the performance of RDI with the Keck/NIRC2 vortex coronagraph in $L'$., Comment: Accepted to AJ. 41 pages, 13 figures, 4 tables. Data presented in this work: https://doi.org/10.5281/zenodo.12747613

Recent Atacama Large Millimeter/submillimeter Array (ALMA) observations of protoplanetary disks in the millimeter continuum have shown a variety of radial gaps, cavities, and spiral features. These substructures may be signposts for ongoing planet formation, and therefore these systems are promising targets for direct imaging planet searches in the near-infrared. To this end, we present results from a deep imaging survey in the $L'$-band (3.8 $\mu$m) with the Keck/NIRC2 vortex coronagraph to search for young planets in 43 disks with resolved features in the millimeter continuum or evidence for gaps/central cavities from their spectral energy distributions. Although we do not detect any new point sources, using the vortex coronagraph allows for high sensitivity to faint sources at small angular separations (down to ${\sim}$0$^{\prime\prime}$.1), allowing us to place strong upper limits on the masses of potential gas giant planets. We compare our mass sensitivities to the masses of planets derived using ALMA observations, and while we are sensitive to $\sim$1 M$_{Jup}$ planets in the gaps in some of our systems, we are generally not sensitive to planets of the masses expected from the ALMA observations. In addition to placing upper limits on the masses of gas giant planets that could be interacting with the dust in the disks to form the observed millimeter substructures, we are also able to map the micron-sized dust as seen in scattered light for 8 of these systems. Our large sample of systems also allows us to investigate limits on planetary accretion rates and disk viscosities., Comment: 23 pages, 14 figures, 3 tables, accepted for publication in AJ

With the advent of JWST, we acquire unprecedented insights into the physical and chemical structure of the inner regions of planet-forming disks where terrestrial planet formation occurs. The very low-mass stars (VLMS) are known to have a high occurrence rate of the terrestrial planets around them. Exploring the chemical composition of the gas in these inner regions of the disks can aid a better understanding of the connection between planet-forming disks and planets. The MIRI mid-Infrared Disk Survey (MINDS) project is a large JWST Guaranteed Time program to characterize the chemistry and physical state of planet-forming and debris disks. We use the JWST-MIRI/MRS spectrum to investigate the gas and dust composition of the planet-forming disk around the very low-mass star Sz28 (M5.5, 0.12\,M$_{\odot}$). We use the dust-fitting tool (DuCK) to determine the dust continuum and to get constraints on the dust composition and grain sizes. We use 0D slab models to identify and fit the molecular spectral features, yielding estimates on the temperature, column density and the emitting area. To test our understanding of the chemistry in the disks around VLMS, we employ the thermo-chemical disk model {P{\tiny RO}D{\tiny I}M{\tiny O}} and investigate the reservoirs of the detected hydrocarbons. We explore how the C/O ratio affects the inner disk chemistry. JWST reveals a plethora of hydrocarbons, including \ce{CH3}, \ce{CH4}, \ce{C2H2}, \ce{^{13}CCH2}, \ce{C2H6}, \ce{C3H4}, \ce{C4H2} and \ce{C6H6} suggesting a disk with a gaseous C/O\,>\,1. Additionally, we detect \ce{CO2}, \ce{^{13}CO2}, \ce{HCN}, and \ce{HC3N}. \ce{H2O} and OH are absent in the spectrum. We do not detect PAHs. Photospheric stellar absorption lines of \ce{H2O} and \ce{CO} are identified. Notably, our radiation thermo-chemical disk models are able to produce these detected hydrocarbons in the surface layers of the disk when the ..., Comment: accepted for publication in A&A

NOTT (formerly Hi-5) is the L'-band (3.5-4.0~microns) nulling interferometer of Asgard, an instrument suite in preparation for the VLTI visitor focus. The primary scientific objectives of NOTT include characterizing (i) young planetary systems near the snow line, a critical region for giant planet formation, and (ii) nearby main-sequence stars close to the habitable zone, with a focus on detecting exozodiacal dust that could obscure Earth-like planets. In 2023-2024, the final warm optics have been procured and assembled in a new laboratory at KU Leuven. First fringes and null measurements were obtained using a Gallium Lanthanum Sulfide (GLS) photonic chip that was also tested at cryogenic temperatures. In this paper, we present an overall update of the NOTT project with a particular focus on the cold mechanical design, the first results in the laboratory with the final NOTT warm optics, and the ongoing Asgard integration activities. We also report on other ongoing activities such as the characterization of the photonic chip (GLS, LiNbO3, SiO), the development of the exoplanet science case, the design of the dispersion control module, and the progress with the self-calibration data reduction software., Comment: 11 pages (incl. 5 figures); Proc. SPIE Astronomical Telescopes + Instrumentation 2024 (Yokohama; Japan), Optical and Infrared Interferometry and Imaging VI

High-contrast imaging (HCI) is a technique designed to observe faint signals near bright sources, such as exoplanets and circumstellar disks. The primary challenge in revealing the faint circumstellar signal near a star is the presence of quasi-static speckles, which can produce patterns on the science images that are as bright, or even brighter, than the signal of interest. Strategies such as angular differential imaging (ADI) or reference-star differential imaging (RDI) aim to provide a means of removing the quasi-static speckles in post-processing. In this paper, we present and discuss the adaptation of state-of-the-art algorithms, initially designed for ADI, to jointly leverage angular and reference-star differential imaging (ARDI) for direct high-contrast imaging of circumstellar disks. Using a collection of high-contrast imaging data sets, we assess the performance of ARDI in comparison to ADI and RDI based on iterative principal component analysis (IPCA). These diverse data sets are acquired under various observing conditions and include the injection of synthetic disk models at various contrast levels. Our results demonstrate that ARDI with IPCA improves the quality of recovered disk images and the sensitivity to planets embedded in disks, compared to ADI or RDI individually. This enhancement is particularly pronounced when dealing with extended sources exhibiting highly ambiguous structures that cannot be accurately retrieved using ADI alone, and when the quality of the reference frames is suboptimal, leading to an underperformance of RDI. We finally apply our method to a sample of real observations of protoplanetary disks taken in star-hopping mode, and propose to revisit the protoplanetary claims associated with these disks., Comment: 21 pages, 18 figures, Accepted to A&A

Very low-mass stars (those <0.3 solar masses) host orbiting terrestrial planets more frequently than other types of stars, but the compositions of those planets are largely unknown. We use mid-infrared spectroscopy with the James Webb Space Telescope to investigate the chemical composition of the planet-forming disk around ISO-ChaI 147, a 0.11 solar-mass star. The inner disk has a carbon-rich chemistry: we identify emission from 13 carbon-bearing molecules including ethane and benzene. We derive large column densities of hydrocarbons indicating that we probe deep into the disk. The high carbon to oxygen ratio we infer indicates radial transport of material within the disk, which we predict would affect the bulk composition of any planets forming in the disk., Comment: Published, 36 pages, 8 figures

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 advent of high-contrast imaging instruments combined with medium-resolution spectrographs allows spectral and temporal dimensions to be combined with spatial dimensions to detect and potentially characterize exoplanets with higher sensitivity. We develop a new method to effectively leverage the spectral and spatial dimensions in integral-field spectroscopy (IFS) datasets using a supervised deep-learning algorithm to improve the detection sensitivity to high-contrast exoplanets. We begin by applying a data transform whereby the IFS datasets are replaced by cross-correlation coefficient tensors obtained by cross-correlating our data with young gas giant spectral template spectra. This transformed data is then used to train machine learning (ML) algorithms. We train a 2D CNN and 3D LSTM with our data. We compare the ML models with a non-ML algorithm, based on the STIM map of arXiv:1810.06895. We test our algorithms on simulated young gas giants in a dataset that contains no known exoplanet, and explore the sensitivity of algorithms to detect these exoplanets at contrasts ranging from 1e-3 to 1e-4 at different radial separations. We quantify the sensitivity using modified receiver operating characteristic curves (mROC). We discover that the ML algorithms produce fewer false positives and have a higher true positive rate than the STIM-based algorithm, and the true positive rate of ML algorithms is less impacted by changing radial separation. We discover that the velocity dimension is an important differentiating factor. Through this paper, we demonstrate that ML techniques have the potential to improve the detection limits and reduce false positives for directly imaged planets in IFS datasets, after transforming the spectral dimension into a radial velocity dimension through a cross-correlation operation., Comment: Accepted for publication in A&A on 23/04/2024. Total 15 pages of text, 7 figures

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

This work aims to measure the mass accretion rate, the accretion luminosity, and more generally the physical conditions of the warm emitting gas in the inner disk of the very low-mass star 2MASS-J16053215-1933159. We investigate the source mid-infrared spectrum for atomic and molecular hydrogen line emission. We present the full James Webb Space Telescope (JWST) Mid-InfraRed Instrument (MIRI) Medium Resolution Spectrometer (MRS) spectrum of the protoplanetary disk around the very low-mass star 2MASS-J16053215-1933159 from the MINDS GTO program, previously shown to be abundant in hydrocarbon molecules. We analyzed the atomic and molecular hydrogen lines in this source by fitting one or multiple Gaussian profiles. We then built a rotational diagram for the H2 lines to constrain the rotational temperature and column density of the gas. Finally, we compared the observed atomic line fluxes to predictions from two standard emission models. We identify five molecular hydrogen pure rotational lines and 16 atomic hydrogen recombination lines. The spectrum indicates optically thin emission for both species. We use the molecular hydrogen lines to constrain the mass and temperature of the warm emitting gas. The HI (7-6) recombination line is used to measure the mass accretion rate and luminosity onto the central source. HI recombination lines can also be used to derive the physical properties of the gas using atomic recombination models. The JWST-MIRI MRS observations for the very low-mass star 2MASS-J16053215-1933159 reveal a large number of emission lines, many originating from atomic and molecular hydrogen because we are able to look into the disk warm molecular layer. Their analysis constrains the physical properties of the emitting gas and showcases the potential of JWST to deepen our understanding of the physical and chemical structure of protoplanetary disks, Comment: 11 pages, 9 figures, 2 tables

The MRS mode of the JWST-MIRI instrument has been shown to be a powerful tool to characterise the molecular gas emission of the inner region of planet-forming disks. Here, we analyse the spectrum of the compact T-Tauri disk DR Tau, which is complemented by high spectral resolution (R~60000-90000) CO ro-vibrational observations. Various molecular species, including CO, CO$_2$, HCN, and C$_2$H$_2$ are detected in the JWST-MIRI spectrum, for which excitation temperatures of T~325-900 K are retrieved using LTE slab models. The high-resolution CO observations allow for a full treatment of the line profiles, which show evidence for two components of the main isotopologue, $^{12}$CO: a broad component tracing the Keplerian disk and a narrow component tracing a slow disk wind. Rotational diagrams yield excitation temperatures of T>725 K for CO, with consistently lower temperatures found for the narrow components, suggesting that the disk wind is launched from a larger distance. The inferred excitation temperatures for all molecules suggest that CO originates from the highest atmospheric layers close to the host star, followed by HCN and C$_2$H$_2$, which emit, together with $^{13}$CO, from slightly deeper layers, whereas the CO$_2$ originates from even deeper inside or further out in the disk. Additional analysis of the $^{12}$CO line wings hint at a misalignment between the inner (i~20 degrees) and outer disk (i~5 degrees). Finally, we emphasise the need for complementary high-resolution CO observations, as in combination with the JWST-MIRI observations they can be used to characterise the CO kinematics and the physical and chemical conditions of the other observed molecules with respect to CO., Comment: Accepted for publication in Astronomy & Astrophysics on 20/03/2024

The study of protoplanetary disks has become increasingly important with the Kepler satellite finding that exoplanets are ubiquitous around stars in our galaxy and the discovery of enormous diversity in planetary system architectures and planet properties. High-resolution near-IR and ALMA images show strong evidence for ongoing planet formation in young disks. The JWST MIRI mid-INfrared Disk Survey (MINDS) aims to (1) investigate the chemical inventory in the terrestrial planet-forming zone across stellar spectral type, (2) follow the gas evolution into the disk dispersal stage, and (3) study the structure of protoplanetary and debris disks in the thermal mid-IR. The MINDS survey will thus build a bridge between the chemical inventory of disks and the properties of exoplanets. The survey comprises 52 targets (Herbig Ae stars, T Tauri stars, very low-mass stars and young debris disks). We primarily obtain MIRI/MRS spectra with high S/N (~100-500) covering the complete wavelength range from 4.9 to 27.9 {\mu}m. For a handful of selected targets we also obtain NIRSpec IFU high resolution spectroscopy (2.87-5.27 {\mu}m). We will search for signposts of planet formation in thermal emission of micron-sized dust - information complementary to near-IR scattered light emission from small dust grains and emission from large dust in the submillimeter wavelength domain. We will also study the spatial structure of disks in three key systems that have shown signposts for planet formation, TW Hya and HD 169142 using the MIRI coronagraph at 15.5 {\mu}m and 10.65 {\mu}m respectively and PDS70 using NIRCam imaging in the 1.87 {\mu}m narrow and the 4.8 {\mu}m medium band filter. ..., Comment: accepted for publication in PASP

Context. Two protoplanets have recently been discovered within the PDS 70 protoplanetary disk. JWST/NIRCam offers a unique opportunity to characterize them and their birth environment at wavelengths difficult to access from the ground. Aims. We aim to image the circumstellar environment of PDS 70 at 1.87 $\mu$m and 4.83 $\mu$m, assess the presence of Pa-$\alpha$ emission due to accretion onto the protoplanets, and probe any IR excess indicative of heated circumplanetary material. Methods. We obtain non-coronagraphic JWST/NIRCam images of PDS 70 within the MINDS (MIRI mid-INfrared Disk Survey) program. We leverage the Vortex Image Processing (VIP) package for data reduction, and develop dedicated routines for optimal stellar PSF subtraction, unbiased imaging of the disk, and protoplanet flux measurement in this type of dataset. A radiative transfer model of the disk is used to disentangle the contributions from the disk and the protoplanets. Results. We re-detect both protoplanets and identify extended emission after subtracting a disk model, including a large-scale spiral-like feature. We interpret its signal in the direct vicinity of planet c as tracing the accretion stream feeding its circumplanetary disk, while the outer part of the feature may rather reflect asymmetric illumination of the outer disk. We also report a bright signal consistent with a previously proposed protoplanet candidate enshrouded in dust, near the 1:2:4 mean-motion resonance with planets b and c. The 1.87 $\mu$m flux of planet b is consistent with atmospheric model predictions, but not that of planet c. We discuss potential origins for this discrepancy, including significant Pa-$\alpha$ line emission. The 4.83 $\mu$m fluxes of planets b and c suggest enshrouding dust or heated CO emission from their circumplanetary environment., Comment: 6+11 pages, 3+10 figures (text+appendix). Accepted for publication in A&A

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

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

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

Looking to the future of exo-Earth imaging from the ground, core technology developments are required in visible extreme adaptive optics (ExAO) to enable the observation of atmospheric features such as oxygen on rocky planets in visible light. UNDERGROUND (Ultra-fast AO techNology Determination for Exoplanet imageRs from the GROUND), a collaboration built in Feb. 2023 at the Optimal Exoplanet Imagers Lorentz Workshop, aims to (1) motivate oxygen detection in Proxima Centauri b and analogs as an informative science case for high-contrast imaging and direct spectroscopy, (2) overview the state of the field with respect to visible exoplanet imagers, and (3) set the instrumental requirements to achieve this goal and identify what key technologies require further development., Comment: SPIE Proceeding: 2023 / 12680-67

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

The understanding of planet formation has changed recently, embracing the new idea of pebble accretion. This means that the influx of pebbles from the outer regions of planet-forming disks to their inner zones could determine the composition of planets and their atmospheres. The solid and molecular components delivered to the planet-forming region can be best characterized by mid-infrared spectroscopy. With Spitzer low-resolution (R=100, 600) spectroscopy, this approach was limited to the detection of abundant molecules such as H2O, C2H2, HCN and CO2. This contribution will present first results of the MINDS (MIRI mid-IR Disk Survey, PI: Th. Henning) project. Due do the sensitivity and spectral resolution (R~1500-3500) provided by JWST we now have a unique tool to obtain the full inventory of chemistry in the inner disks of solar-types stars and brown dwarfs, including also less abundant hydrocarbons and isotopologues. The Integral Field Unit (IFU) capabilities enable at the same time spatial studies of the continuum and line emission in extended sources such as debris disks, the flying saucer and also the search for mid-IR signatures of forming planets in systems such as PDS70. These JWST observations are complementary to ALMA and NOEMA observations of the outer disk chemistry; together these datasets provide an integral view of the processes occurring during the planet formation phase., Comment: 14 pages, 8 figures, authors version of the manuscript submitted 22.1.2023 for the Faraday Discussion "Astrochemistry at high resolution", accepted 21.3.2023

NASA is engaged in planning for a Habitable Worlds Observatory (HabWorlds), a coronagraphic space mission to detect rocky planets in habitable zones and establish their habitability. Surface liquid water is central to the definition of planetary habitability. Photometric and polarimetric phase curves of starlight reflected by an exoplanet can reveal ocean glint, rainbows and other phenomena caused by scattering by clouds or atmospheric gas. Direct imaging missions are optimised for planets near quadrature, but HabWorlds' coronagraph may obscure the phase angles where such optical features are strongest. The range of accessible phase angles for a given exoplanet will depend on the planet's orbital inclination and/or the coronagraph's inner working angle (IWA). We use a recently-created catalog relevant to HabWorlds of 164 stars to estimate the number of exo-Earths that could be searched for ocean glint, rainbows, and polarization effects due to Rayleigh scattering. We find that the polarimetric Rayleigh scattering peak is accessible in most of the exo-Earth planetary systems. The rainbow due to water clouds at phase angles of ${\sim}20-60^\circ$ would be accessible with HabWorlds for a planet with an Earth equivalent instellation in ${\sim}{46}$ systems, while the ocean glint signature at phase angles of ${\sim}130-170^\circ$ would be accessible in ${\sim}{16}$ systems, assuming an IWA${=}62$ mas ($3\lambda/D$). Improving the IWA${=}41$ mas ($2\lambda/D$) increases accessibility to rainbows and glints by factors of approximately 2 and 3, respectively. By observing these scattering features, HabWorlds could detect a surface ocean and water cycle, key indicators of habitability., Comment: MNRAS accepted, 9 pages, 8 figures, 3 tables

Terrestrial and sub-Neptune planets are expected to form in the inner ($<10~$AU) regions of protoplanetary disks. Water plays a key role in their formation, although it is yet unclear whether water molecules are formed in-situ or transported from the outer disk. So far Spitzer Space Telescope observations have only provided water luminosity upper limits for dust-depleted inner disks, similar to PDS 70, the first system with direct confirmation of protoplanet presence. Here we report JWST observations of PDS 70, a benchmark target to search for water in a disk hosting a large ($\sim54~$AU) planet-carved gap separating an inner and outer disk. Our findings show water in the inner disk of PDS 70. This implies that potential terrestrial planets forming therein have access to a water reservoir. The column densities of water vapour suggest in-situ formation via a reaction sequence involving O, H$_2$, and/or OH, and survival through water self-shielding. This is also supported by the presence of CO$_2$ emission, another molecule sensitive to UV photodissociation. Dust shielding, and replenishment of both gas and small dust from the outer disk, may also play a role in sustaining the water reservoir. Our observations also reveal a strong variability of the mid-infrared spectral energy distribution, pointing to a change of inner disk geometry., Comment: To appear in Nature on 24 July 2023. 21 pages, 10 figures; includes extended data. Part of the JWST MINDS Guaranteed Time Observations program's science enabling products. Spectra downloadable on Zenodo at https://zenodo.org/record/7991022

MIRI/MRS on board the JWST allows us to probe the inner regions of protoplanetary disks. Here we examine the disk around the classical T Tauri star Sz 98, which has an unusually large dust disk in the millimetre with a compact core. We focus on the H$_2$O emission through both its ro-vibrational and pure rotational emission. Furthermore, we compare our chemical findings with those obtained for the outer disk from Atacama Large Millimeter/submillimeter Array (ALMA) observations. In order to model the molecular features in the spectrum, the continuum was subtracted and LTE slab models were fitted. The spectrum was divided into different wavelength regions corresponding to H$_2$O lines of different excitation conditions, and the slab model fits were performed individually per region. We confidently detect CO, H$_2$O, OH, CO$_2$, and HCN in the emitting layers. The isotopologue H$^{18}_2$O is not detected. Additionally, no other organics, including C$_2$H$_2$, are detected. This indicates that the C/O ratio could be substantially below unity, in contrast with the outer disk. The H$_2$O emission traces a large radial disk surface region, as evidenced by the gradually changing excitation temperatures and emitting radii. The OH and CO$_2$ emission are relatively weak. It is likely that H$_2$O is not significantly photodissociated; either due to self-shielding against the stellar irradiation, or UV-shielding from small dust particles. The relative emitting strength of the different identified molecular features point towards UV-shielding of H$_2$O in the inner disk of Sz 98, with a thin layer of OH on top. The majority of the organic molecules are either hidden below the dust continuum, or not present. In general, the inferred composition points to a sub-solar C/O ratio (<0.5) in the inner disk, in contrast with the larger than unity C/O ratio in the gas in the outer disk found with ALMA., Comment: Updated to accepted version. Accepted by A&A on September 23 2023. 18 pages, 11 figures

Carbon is an essential element for life but how much can be delivered to young planets is still an open question. The chemical characterization of planet-forming disks is a crucial step in our understanding of the diversity and habitability of exoplanets. Very low-mass stars ($<0.2~M_{\odot}$) are interesting targets because they host a rich population of terrestrial planets. Here we present the JWST detection of abundant hydrocarbons in the disk of a very low-mass star obtained as part of the MIRI mid-INfrared Disk Survey (MINDS). In addition to very strong and broad emission from C$_2$H$_2$ and its $^{13}$C$^{12}$CH$_2$ isotopologue, C$_4$H$_2$, benzene, and possibly CH$_4$ are identified, but water, PAH and silicate features are weak or absent. The lack of small silicate grains implies that we can look deep down into this disk. These detections testify to an active warm hydrocarbon chemistry with a high C/O ratio in the inner 0.1 au of this disk, perhaps due to destruction of carbonaceous grains. The exceptionally high C$_2$H$_2$/CO$_2$ and C$_2$H$_2$/H$_2$O column density ratios suggest that oxygen is locked up in icy pebbles and planetesimals outside the water iceline. This, in turn, will have significant consequences for the composition of forming exoplanets., Comment: version submitted to Nature Astronomy

We propose to use low-rank matrix approximation using the component-wise L1-norm for direct imaging of exoplanets. Exoplanet detection by direct imaging is a challenging task for three main reasons: (1) the host star is several orders of magnitude brighter than exoplanets, (2) the angular distance between exoplanets and star is usually very small, and (3) the images are affected by the noises called speckles that are very similar to the exoplanet signal both in shape and intensity. We first empirically examine the statistical noise assumptions of the L1 and L2 models, and then we evaluate the performance of the proposed L1 low-rank approximation (L1-LRA) algorithm based on visual comparisons and receiver operating characteristic (ROC) curves. We compare the results of the L1-LRA with the widely used truncated singular value decomposition (SVD) based on the L2 norm in two different annuli, one close to the star and one far away., Comment: 13 pages, 4 figures, BNAIC/BeNeLearn 2023

Direct imaging of exoplanets is a challenging task due to the small angular distance and high contrast relative to their host star, and the presence of quasi-static noise. We propose a new statistical method for direct imaging of exoplanets based on a likelihood ratio detection map, which assumes that the noise after the background subtraction step obeys a Laplacian distribution. We compare the method with two detection approaches based on signal-to-noise ratio (SNR) map after performing the background subtraction by the widely used Annular Principal Component Analysis (AnnPCA). The experimental results on the Beta Pictoris data set show the method outperforms SNR maps in terms of achieving the highest true positive rate (TPR) at zero false positive rate (FPR).

[abridged] We analyzed mid-infrared high-contrast coronagraphic images of the beta Pictoris system, taking advantage of the NEAR experiment using the VLT/VISIR instrument. The goal of our analysis is to investigate both the detection of the planet beta Pictoris b and of the disk features at mid-IR wavelengths. In addition, by combining several epochs of observation, we expect to constrain the position of the known clumps and improve our knowledge on the dynamics of the disk. To evaluate the planet b flux contribution, we extracted the photometry and compared it to the flux published in the literature. In addition, we used previous data from T-ReCS and VISIR, to study the evolution of the position of the southwest clump that was initially observed in the planetary disk back in 2003. While we did not detect the planet b, we were able to put constraints on the presence of circumplanetary material, ruling out the equivalent of a Saturn-like planetary ring around the planet. The disk presents several noticeable structures, including the known southwest clump. Using a 16-year baseline, sampled with five epochs of observations, we were able to examine the evolution of the clump: the clump orbits in a Keplerian motion with an sma of 56.1+-0.4 au. In addition to the known clump, the images clearly show the presence of a second clump on the northeast side of the disk and fainter and closer structures that are yet to be confirmed. We found correlations between the CO clumps detected with ALMA and the mid-IR images. If the circumplanetary material were located at the Roche radius, the maximum amount of dust determined from the flux upper limit around beta Pictoris b would correspond to the mass of an asteroid of 5 km in diameter. Finally, the Keplerian motion of the southwestern clump is possibly indicative of a yet-to-be-detected planet or signals the presence of a vortex., Comment: Accepted in Astronomy and Astrophysics

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!

The mid-infrared (IR) regime is well suited to directly detect the thermal signatures of exoplanets in our solar neighborhood. The NEAR experiment: demonstration of high-contrast imaging (HCI) capability at ten microns, can reach sub-mJy detection sensitivity in a few hours of observation time, which is sufficient to detect a few Jupiter mass planets in nearby systems. One of the big limitations for HCI in the mid-IR is thermal sky-background. In this work, we show that precipitate water vapor (PWV) is the principal contributor to thermal sky background and science PSF quality. In the presence of high PWV, the HCI performance is significantly degraded in the background limited regime., Comment: 8 pages, 4 figures, conference proceedings (SPIE Astronomical telescopes and instrumentation 2022)

Supervised deep learning was recently introduced in high-contrast imaging (HCI) through the SODINN algorithm, a convolutional neural network designed for exoplanet detection in angular differential imaging (ADI) datasets. The benchmarking of HCI algorithms within the Exoplanet Imaging Data Challenge (EIDC) showed that (i) SODINN can produce a high number of false positives in the final detection maps, and (ii) algorithms processing images in a more local manner perform better. This work aims to improve the SODINN detection performance by introducing new local processing approaches and adapting its learning process accordingly. We propose NA-SODINN, a new deep learning binary classifier based on a convolutional neural network (CNN) that better captures image noise correlations in ADI-processed frames by identifying noise regimes. Our new approach was tested against its predecessor, as well as two SODINN-based hybrid models and a more standard annular-PCA approach, through local receiving operating characteristics (ROC) analysis of ADI sequences from the VLT/SPHERE and Keck/NIRC-2 instruments. Results show that NA-SODINN enhances SODINN in both sensitivity and specificity, especially in the speckle-dominated noise regime. NA-SODINN is also benchmarked against the complete set of submitted detection algorithms in EIDC, in which we show that its final detection score matches or outperforms the most powerful detection algorithms.Throughout the supervised machine learning case, this study illustrates and reinforces the importance of adapting the task of detection to the local content of processed images., Comment: A&A in press

We propose a direct imaging method for the detection of exoplanets based on a combined low-rank plus structured sparse model. For this task, we develop a dictionary of possible effective circular trajectories a planet can take during the observation time, elements of which can be efficiently computed using rotation and convolution operation. We design a simple alternating iterative hard-thresholding algorithm that jointly promotes a low-rank background and a sparse exoplanet foreground, to solve the non-convex optimisation problem. The experimental comparison on the $\beta$-Pictoris exoplanet benchmark dataset shows that our method has the potential to outperform the widely used Annular PCA for specific planet light intensities in terms of the Receiver operating characteristic (ROC) curves., Comment: 5 pages, 4 figures, ICASSP 2023

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

We present JWST-MIRI MRS spectra of the protoplanetary disk around the low-mass T Tauri star GW Lup from the MIRI mid-INfrared Disk Survey (MINDS) GTO program. Emission from $^{12}$CO$_{2}$, $^{13}$CO$_{2}$, H$_{2}$O, HCN, C$_{2}$H$_{2}$, and OH is identified with $^{13}$CO$_{2}$ being detected for the first time in a protoplanetary disk. We characterize the chemical and physical conditions in the inner few au of the GW Lup disk using these molecules as probes. The spectral resolution of JWST-MIRI MRS paired with high signal-to-noise data is essential to identify these species and determine their column densities and temperatures. The $Q$-branches of these molecules, including those of hot-bands, are particularly sensitive to temperature and column density. We find that the $^{12}$CO$_{2}$ emission in the GW Lup disk is coming from optically thick emission at a temperature of $\sim$400 K. $^{13}$CO$_{2}$ is optically thinner and based on a lower temperature of $\sim$325 K, may be tracing deeper into the disk and/or a larger emitting radius than $^{12}$CO$_{2}$. The derived $N_{\rm{CO_{2}}}$/$N_{\rm{H_{2}O}}$ ratio is orders of magnitude higher than previously derived for GW Lup and other targets based on \textit{Spitzer}-IRS data. This high column density ratio may be due to an inner cavity with a radius in between the H$_{2}$O and CO$_{2}$ snowlines and/or an overall lower disk temperature. This paper demonstrates the unique ability of JWST to probe inner disk structures and chemistry through weak, previously unseen molecular features., Comment: 15 pages, 10 figures. Accepted to ApJL

Observing dynamical interactions between planets and disks is key to understanding their formation and evolution. Two protoplanets have recently been discovered within PDS 70's protoplanetary disk, along with an arm-like structure towards the north-west of the star. Our aim is to constrain the morphology and origin of this arm-like structure, and to assess whether it could trace a spiral density wave caused by the dynamical interaction between the planet PDS 70c and the disk. We analyze polarized and angular differential imaging (PDI and ADI) data taken with VLT/SPHERE, spanning six years of observations. PDI data sets are reduced using the IRDAP polarimetric data reduction pipeline, while ADI data sets are processed using MUSTARD, a novel inverse problem algorithm to tackle the geometrical biases spoiling the images previously used for the analysis of this disk. We confirm the presence of the arm-like structure in all PDI and ADI datasets. We do not observe a south-east symmetric arm with respect to the disk minor axis, which seems to disfavor the previous hypothesis that the arm is the footprint of a double-ring structure. If the structure traces a spiral density wave following the motion of PDS 70c, we would expect $11\overset{\circ}{.}28^{+2\overset{\circ}{.}20}_{-0\overset{\circ}{.}86}$ rotation for the spiral in six years. However, we do not measure any significant movement of the structure. If the arm-like structure is a planet-driven spiral arm, the observed lack of rotation would suggest that the assumption of rigid-body rotation may be inappropriate for spirals induced by planets. We suggest that the arm-like structure may rather trace a vortex appearing as a one-armed spiral in scattered light due to projection effects. The vortex hypothesis accounts for both the lack of observed rotation and the presence of a nearby sub-mm continuum asymmetry detected with ALMA., Comment: Accepted for publication in Astronomy&Astrophysics. 11 pages, 10 figures

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

In the past decade, HCI surveys provided new insights about the frequency and properties of substellar companions at separation larger than 5 au. In this context, our study aims to detect and characterise potential exoplanets and brown dwarfs within debris disks, by considering the SHARDDS survey, which gathers 55 Main Sequence stars with known bright debris disk. We rely on the AutoRSM framework to perform an in-depth analysis of the targets, via the computation of detection maps and contrast curves. A clustering approach is used to divide the set of targets in multiple subsets, in order to reduce the computation time by estimating a single optimal parametrisation for each considered subset. The use of Auto-RSM allows to reach high contrast at short separations, with a median contrast of 10-5 at 300 mas, for a completeness level of 95%. Detection maps generated with different approaches are used along with contrast curves, to identify potential planetary companions. A new planetary characterisation algorithm, based on the RSM framework, is developed and tested successfully, showing a higher astrometric and photometric precision for faint sources compared to standard approaches. Apart from the already known companion of HD206893 and two point-like sources around HD114082 which are most likely background stars, we did not detect any new companion around other stars. A correlation study between achievable contrasts and parameters characterising HCI sequences highlights the importance of the strehl, wind speed and wind driven halo to define the quality of high contrast images. Finally, planet detection and occurrence frequency maps are generated and show, for the SHARDDS survey, a high detection rate between 10 and 100 au for substellar companions with mass >10MJ.

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

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

Most of the high-contrast imaging (HCI) data-processing techniques used over the last 15 years have relied on the angular differential imaging (ADI) observing strategy, along with subtraction of a reference point spread function (PSF) to generate exoplanet detection maps. Recently, a new algorithm called regime switching model (RSM) map has been proposed to take advantage of these numerous PSF-subtraction techniques; RSM uses several of these techniques to generate a single probability map. Selection of the optimal parameters for these PSF-subtraction techniques as well as for the RSM map is not straightforward, is time consuming, and can be biased by assumptions made as to the underlying data set. We propose a novel optimisation procedure that can be applied to each of the PSF-subtraction techniques alone, or to the entire RSM framework. The optimisation procedure consists of three main steps: (i) definition of the optimal set of parameters for the PSF-subtraction techniques using the contrast as performance metric, (ii) optimisation of the RSM algorithm, and (iii) selection of the optimal set of PSF-subtraction techniques and ADI sequences used to generate the final RSM probability map. The optimisation procedure is applied to the data sets of the exoplanet imaging data challenge (EIDC), which provides tools to compare the performance of HCI data-processing techniques. The data sets consist of ADI sequences obtained with three state-of-the-art HCI instruments: SPHERE, NIRC2, and LMIRCam. The results of our analysis demonstrate the interest of the proposed optimisation procedure, with better performance metrics compared to the earlier version of RSM, as well as to other HCI data-processing techniques.

$\epsilon$~Eridani is a young planetary system hosting a complex multi-belt debris disk and a confirmed Jupiter-like planet orbiting at 3.48 AU from its host star. Its age and architecture are thus reminiscent of the early Solar System. The most recent study of Mawet et al. 2019, which combined radial velocity (RV) data and Ms-band direct imaging upper limits, started to constrain the planet's orbital parameters and mass, but are still affected by large error bars and degeneracies. Here we make use of the most recent data compilation from three different techniques to further refine $\epsilon$~Eridani~b's properties: RVs, absolute astrometry measurements from the Hipparcos~and Gaia~missions, and new Keck/NIRC2 Ms-band vortex coronagraph images. We combine this data in a Bayesian framework. We find a new mass, $M_b$ = $0.66_{-0.09}^{+0.12}$~M$_{Jup}$, and inclination, $i$ = $77.95_{-21.06}^{\circ+28.50}$, with at least a factor 2 improvement over previous uncertainties. We also report updated constraints on the longitude of the ascending node, the argument of the periastron, and the time of periastron passage. With these updated parameters, we can better predict the position of the planet at any past and future epoch, which can greatly help define the strategy and planning of future observations and with subsequent data analysis. In particular, these results can assist the search for a direct detection with JWST and the Nancy Grace Roman Space Telescope's coronagraph instrument (CGI)., Comment: Published in AJ

Giant exoplanets on 10-100 au orbits have been directly imaged around young stars. The peak of the thermal emission from these warm young planets is in the near-infrared (~1-5 microns), whereas mature, temperate exoplanets (i.e., those within their stars' habitable zones) radiate primarily in the mid-infrared (mid-IR: ~10 microns). If the background noise in the mid-IR can be mitigated, then exoplanets with low masses--including rocky exoplanets--can potentially be imaged in very deep exposures. Here, we review the recent results of the Breakthrough Watch/New Earths in the Alpha Centauri Region (NEAR) program on the Very Large Telescope (VLT) in Chile. NEAR pioneered a ground-based mid-IR observing approach designed to push the capabilities for exoplanet imaging with a specific focus on the closest stellar system, Alpha Centauri. NEAR combined several new optical technologies--including a mid-IR optimized coronagraph, adaptive optics system, and rapid chopping strategy to mitigate noise from the central star and thermal background within the habitable zone. We focus on the lessons of the VLT/NEAR campaign to improve future instrumentation--specifically on strategies to improve noise mitigation through chopping. We also present the design and commissioning of the Large Binocular Telescope's Exploratory Survey for Super-Earths Orbiting Nearby Stars (LESSONS), an experiment in the Northern hemisphere that is building on what was learned from NEAR to further push the sensitivity of mid-IR imaging. Finally, we briefly discuss some of the possibilities that mid-IR imaging will enable for exoplanet science., Comment: Submitted proceedings to SPIE Optical Engineering + Applications 2021, Techniques and Instrumentation for Detection of Exoplanets X

Direct imaging is a powerful exoplanet discovery technique that is complementary to other techniques and offers great promise in the era of 30 meter class telescopes. Space-based transit surveys have revolutionized our understanding of the frequency of planets at small orbital radii around Sun-like stars. The next generation of extremely large ground-based telescopes will have the angular resolution and sensitivity to directly image planets with $R < 4R_\oplus$ around the very nearest stars. Here, we predict yields from a direct imaging survey of a volume-limited sample of Sun-like stars with the Mid-Infrared ELT Imager and Spectrograph (METIS) instrument, planned for the 39 m European Southern Observatory (ESO) Extremely Large Telescope (ELT) that is expected to be operational towards the end of the decade. Using Kepler occurrence rates, a sample of stars with spectral types A-K within 6.5 pc, and simulated contrast curves based on an advanced model of what is achievable from coronagraphic imaging with adaptive optics, we estimated the expected yield from METIS using Monte Carlo simulations. We find the METIS expected yield of planets in the N2 band (10.10 - 12.40 $\mu$m) is 1.14 planets, which is greater than comparable observations in the L (3.70 - 3.95 $\mu$m) and M (4.70 - 4.90 $\mu$m) bands. We also determined a 24.6% chance of detecting at least one Jovian planet in the background limited regime assuming a 1 hour integration. We calculated the yield per star and estimate optimal observing revisit times to increase the yield. We also analyzed a northern hemisphere version of this survey and found there are additional targets worth considering. In conclusion, we present an observing strategy aimed to maximize the possible yield for limited telescope time, resulting in 1.48 expected planets in the N2 band., Comment: 11 pages, 6 Figures, 4 Tables, A&A 2021 in press. Revision: Minor clarifications implemented during publisher review

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