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

The Mid-infrared ELT Imager and Spectrograph (METIS) is a first-generation instrument for the Extremely Large Telescope (ELT), Europe's next-generation 39 m ground-based telescope for optical and infrared wavelengths. METIS will offer diffraction-limited imaging, low- and medium-resolution slit spectroscopy, and coronagraphy for high-contrast imaging between 3 and 13 microns, as well as high-resolution integral field spectroscopy between 3 and 5 microns. The main METIS science goals are the detection and characterisation of exoplanets, the investigation of proto-planetary disks, and the formation of planets. The Single-Conjugate Adaptive Optics (SCAO) system corrects atmospheric distortions and is thus essential for diffraction-limited observations with METIS. Numerous challenging aspects of an ELT Adaptive Optics (AO) system are addressed in the mature designs for the SCAO control system and the SCAO hardware module: the complex interaction with the telescope entities that participate in the AO control, wavefront reconstruction with a fragmented and moving pupil, secondary control tasks to deal with differential image motion, non-common path aberrations and mis-registration. A K-band pyramid wavefront sensor and a GPU-based Real-Time Computer (RTC), tailored to the needs of METIS at the ELT, are core components. This current paper serves as a natural sequel to our previous work presented in Hippler et al. (2018). It includes updated performance estimations in terms of several key performance indicators, including achieved contrast curves. We outline all important design decisions that were taken, and present the major challenges we faced and the main analyses carried out to arrive at these decisions and eventually the final design. We also elaborate on our testing and verification strategy, and, last not least, comprehensively present the full design, hardware and software., Comment: 51 pages, 29 figures, accepted to Experimental Astronomy

Ring resonators are crucial in silicon photonics for various applications, but conventional designs face performance trade-offs. Here a third-order polynomial interconnected circular (TOPIC) bend is proposed to revolutionize the ring designs fundamentally. The TOPIC bend has a unique feature of continuous curvature and curvature derivative, which is theoretically derived to be essential for waveguide loss optimization. With the TOPIC bend, the silicon ring resonators demonstrated here have achieved three records to the best of our knowledge: the smallest radius (0.7 $\mathrm{\mu m}$) for silicon rings resonating with single guided mode, the lowest thermal tuning power (5.85 mW/$\pi$) for silicon rings with FSR $\geq$3.2 THz, and the first silicon ring-based WDM 32$\times$100 GHz filter. The filter has doubled the channel amount compared to the state of the art, and meanwhile achieved low insertion loss (1.91 $\pm$ 0.28 dB) and low tuning power (283 GHz/mW). Moreover, the TOPIC bend is not limited to ring applications, it can also be used to create bends with an arbitrary angle, with the advantages of ultra-compact radius and heater integration, which are expected to replace all circular bends in integrated photonics, greatly reducing system size and power consumption.

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

This paper considers the problem of minimizing a differentiable function with locally Lipschitz continuous gradient on the algebraic variety of all $m$-by-$n$ real matrices of rank at most $r$. Several definitions of stationarity exist for this nonconvex problem. Among them, Bouligand stationarity is the strongest necessary condition for local optimality. Only a handful of algorithms generate a sequence in the variety whose accumulation points are provably Bouligand stationary. Among them, the most parsimonious with (truncated) singular value decompositions (SVDs) or eigenvalue decompositions can still require a truncated SVD of a matrix whose rank can be as large as $\min\{m, n\}-r+1$ if the gradient does not have low rank, which is computationally prohibitive in the typical case where $r \ll \min\{m, n\}$. This paper proposes a first-order algorithm that generates a sequence in the variety whose accumulation points are Bouligand stationary while requiring SVDs of matrices whose smaller dimension is always at most $r$. A standard measure of Bouligand stationarity converges to zero along the bounded subsequences at a rate at least $O(1/\sqrt{i+1})$, where $i$ is the iteration counter. Furthermore, a rank-increasing scheme based on the proposed algorithm is presented, which can be of interest if the parameter $r$ is potentially overestimated.

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

We give bounds on geodesic distances on the Stiefel manifold, derived from new geometric insights. The considered geodesic distances are induced by the one-parameter family of Riemannian metrics introduced by H\"uper et al. (2021), which contains the well-known Euclidean and canonical metrics. First, we give the best Lipschitz constants between the distances induced by any two members of the family of metrics. Then, we give a lower and an upper bound on the geodesic distance by the easily computable Frobenius distance. We give explicit families of pairs of matrices that depend on the parameter of the metric and the dimensions of the manifold, where the lower and the upper bound are attained. These bounds aim at improving the theoretical guarantees and performance of minimal geodesic computation algorithms by reducing the initial velocity search space. In addition, these findings contribute to advancing the understanding of geodesic distances on the Stiefel manifold and their applications., Comment: 27 pages, 11 figures

To leverage the angular resolution of interferometry at high contrast, one must employ specialized beam-combiners called interferometric nullers. Nullers discard part of the astrophysical information to optimize the recording of light present in the dark fringe of the central source. Asgard/NOTT will deploy a beam-combination scheme offering good instrumental noise rejection when phased appropriately, but for which information is degenerate on the outputs, prompting a dedicated tuning strategy using the science detector. The dispersive effect of water vapor can be corrected with prisms forming a variable thickness of glass. But observations in the L band suffer from an additional and important chromatic effect due to longitudinal atmospheric dispersion coming from a resonance of CO2 at 4.3 micron. To compensate for this effect efficiently, a novel type of compensation device will be deployed leveraging a gas cell of variable length at ambient pressure. After reviewing the impact of water vapor and CO2, we present the design of this atmospheric dispersion compensation device and describe a strategy to maintain this tuning on-sky., Comment: 7 pages, 5 figures, Proc of SPIE Astronomical Telescopes + Instrumentation 2024 Yokohama, Japan, Optical and Infrared Interferometry and Imaging VI

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

High contrast imaging (HCI) is fundamentally limited by wavefront aberrations, and the ability to perform wavefront sensing from focal plane images is key to reach the full potential of ground and space-based instruments. Vortex focal plane mask coupled with downstream pupil (Lyot) stop stands as one of the best small-angle coronagraphs, but is also sensitive to low-order aberrations. Here, we revisit the behavior of the vortex phase mask, from entrance pupil down to the final detector plane, with Zernike polynomials as input phase aberrations. In particular we develop a second-order expansion that allows us to analyze the phase retrieval properties in a more intuitive and accurate way than previously proposed. With this formalism, we show how the azimuthal vortex modulation modifies the phase retrieval properties compared to normal imaging. In particular, our results suggest that images obtained with a scalar vortex coronagraph can be used for unambiguous focal-plane wavefront sensing in any practical situation. We compare our results with numerical simulations and discuss practical implementation in coronagraphic instruments., Comment: 7 pages, 7 figures, paper presented at SPIE Astronomical Telescopes + Instrumentation 2024

Non-common path quasi-static and differential aberrations are one of the big hurdles of direct imaging for current and future high-contrast imaging instruments. They increase speckle and photon noise thus reducing the achievable contrast and lead to a significant hit in HCI performance. The Mid-infrared ELT Imager and Spectrograph (METIS) will provide high-contrast imaging, including vortex coronagraphy in L, M and N bands, with the ultimate goal of directly imaging temperate rocky planets around the nearest stars. Ground-based mid-infrared observations are however also impacted by water vapor inhomogeneities in the atmosphere, which generate additional chromatic turbulence not corrected by the near-infrared adaptive optics. This additional source of wavefront error (WFE) significantly impacts HCI performance, and even dominates the WFE budget in N band. Instantaneous focal plane wavefront sensing is thus required to mitigate its impact. In this context, we propose to implement a novel wavefront sensing approach for the vortex coronagraph using an asymmetric Lyot stop and machine learning. The asymmetric pupil stop allows for the problem to become solvable, lifting the ambiguity on the sign of even Zernike modes. Choosing the Lyot plane instead of the entrance pupil for this mask is also not arbitrary: it preserves the rejection efficiency of the coronagraph and minimizes the impact of the asymmetry on the throughput. Last but not least, machine learning allows us to solve this inversion problem which is non-linear and lacks an analytical solution. In this contribution, we present our concept, our simulation framework, our results and a first laboratory demonstration of the technique., Comment: 8 pages, 4 figures, paper presented at SPIE Astronomical Telescopes + Instrumentation 2024

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

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

Nulling interferometry is a powerful observing technique to reach exoplanets and circumstellar dust at separations too small for direct imaging with single-dish telescopes and too large for indirect methods. With near-future instrumentation, it bears the potential to detect young, hot planets near the snow lines of their host stars. A future space mission could detect and characterize a large number of rocky, habitable-zone planets around nearby stars at thermal-infrared wavelengths. The null self-calibration is a method aiming at modelling the statistical distribution of the nulled signal. It has proven to be more sensitive and accurate than average-based data reduction methods in nulling interferometry. This statistical approach opens the possibility of designing a GPU-based Python package to reduce the data from any of these instruments, by simply providing the data and a simulator of the instrument. GRIP is a toolbox to reduce nulling and interferometric data based on the statistical self-calibration method. In this article, we present the main features of GRIP as well as applications on real data.

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

SAXO+ is a planned enhancement of the existing SAXO, the VLT/ SPHERE adaptive optics system, deployed on ESO's Very Large Telescope. This upgrade is designed to significantly enhance the instrument's capacity to detect and analyze young Jupiter-like planets. The pivotal addition in SAXO+ is a second-stage adaptive optics system featuring a dedicated near-infrared pyramid wavefront sensor and a second deformable mirror. This secondary stage is strategically integrated to address any residual wavefront errors persisting after the initial correction performed by the current primary AO loop, SAXO. However, several recent studies clearly showed that in good conditions, even in the current system SAXO, non-common path aberrations (NCPAs) are the limiting factor of the final normalized intensity in focal plane, which is the final metric for ground-based high-contrast instruments. This is likely to be even more so the case with the new AO system, with which the AO residuals will be minimized. Several techniques have already been extensively tested on SPHERE in internal source and/or on-sky and will be presented in this paper. However, the use of a new type of sensor for the second stage, a pyramid wavefront sensor, will likely complicate the correction of these aberrations. Using an end-to-end AO simulation tool, we conducted simulations to gauge the effect of measured SPHERE NCPAs in the coronagraphic image on the second loop system and their correction using focal plane wavefront sensing systems. We finally analyzed how the chosen position of SAXO+ in the beam will impact the evolution of the NCPAs in the new instrument., Comment: 16 pages, 10 figures, submitted to the proceedings of SPIE Astronomical Telescopes + Instrumentation 2024, 13096-357

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

Optimization over the set of matrices $X$ that satisfy $X^\top B X = I_p$, referred to as the generalized Stiefel manifold, appears in many applications involving sampled covariance matrices such as the canonical correlation analysis (CCA), independent component analysis (ICA), and the generalized eigenvalue problem (GEVP). Solving these problems is typically done by iterative methods that require a fully formed $B$. We propose a cheap stochastic iterative method that solves the optimization problem while having access only to random estimates of $B$. Our method does not enforce the constraint in every iteration; instead, it produces iterations that converge to critical points on the generalized Stiefel manifold defined in expectation. The method has lower per-iteration cost, requires only matrix multiplications, and has the same convergence rates as its Riemannian optimization counterparts that require the full matrix $B$. Experiments demonstrate its effectiveness in various machine learning applications involving generalized orthogonality constraints, including CCA, ICA, and the GEVP., Comment: This v3 is a corrected version of the ICML 2024 paper (PMLR 235:49226-49248); see the errata at the end

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

We demonstrate a design for a high-performance $2 \times 2$ splitter meeting the essential requirements of broadband coupling, support for arbitrary coupling ratio, ultra low-loss, high fabrication tolerance, and a compact footprint. This is achieved based on a rigorous coupled mode theory analysis of the broadband response of the bent directional coupler (DC) and by demonstrating a full coupling model, with measured broadband values of 0.4, 0.5, 0.6, and 0.7. As a benchmark, we demonstrate a 0.5:0.5 splitter that significantly reduces coupling variation from 0.391 in the traditional DC to just 0.051 over an 80 nm wavelength span. This represents a remarkable 7.67 times reduction in coupling variation. Further, newly-invented low-loss bends were used in the proposed design leading to an ultra low-loss design with negligible excess loss ($\mathrm{0.003 \pm 0.013 \ dB}$). The proposed 0.5:0.5 silicon strip waveguide-based design is tolerant and shows consistently low coupling variation over a full 300 mm wafer showcasing a maximum cross coupling variation of 0.112 over 80 nm wavelength range, at the extreme edge of the wafer. Futhermore, we augmented the wafer mapping with a waveguide width fabrication tolerance study, confirming the tolerance of the device with a mere 0.061 maximum coupling variation with a waveguide width deviation of $\pm 20$ nm over 80 nm wavelength range. These specs make the proposed splitter an attractive component for practical applications with mass production.

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

We exhibit conjugate points on the Stiefel manifold endowed with any member of the family of Riemannian metrics introduced by H\"uper et al. (2021). This family contains the well-known canonical and Euclidean metrics. An upper bound on the injectivity radius of the Stiefel manifold in the considered metric is then obtained as the minimum between the length of the geodesic along which the points are conjugate and the length of certain geodesic loops. Numerical experiments support the conjecture that the obtained upper bound is in fact equal to the injectivity radius., Comment: v2: fixed MSC codes; fixed typo in the penultimate sentence of the proof of Theorem 6.1; added section 9 (Concluding remarks)

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

To the best of our knowledge, this paper has achieved the lowest thermal tuning power (5.85 mW/pi) for silicon rings with FSR>=3.2 THz, and the first silicon ring-based WDM-32x100 GHz filter., Comment: 3 pages submitted for conference

We report on the realization of a quantum sensor based on trapped atom interferometry in an optical lattice for the measurement of atom-surface interactions, with sub-micrometer-level control of the mean atom-surface separation distance. The force sensor reaches a short-term sensitivity of 3.4 x 10 --28 N at 1 s and a long-term stability of 4 qN (4 x 10 --30 N). We perform force measurements in the 0-300 $\mu$m range, and despite significant stray forces caused by adsorbed atoms on the surface, we obtain evidence of the Casimir-Polder force. Short-range forces are one of the many frontiers of modern physics [1, 2]. In the submillimeter scales, quantum electrodynamics (QED) interactions are dominant, and give rise in the case of atom-surface interactions to the Casimir-Polder force [3]. Since the first highlight of this force [4], several different methods [5] have been able to bring out Casimir-Polder forces, notably by measuring the transmission of an atomic beam through a micronsized cavity [6], diffracting matter waves on a surface [7] or performing spectroscopy in vapor cells [8, 9]. However these approaches have struggled to achieve the high measurement sensitivity required to detect the very weak forces involved all while maintaining a good understanding of the setup geometry, particularly the distance separating atoms from the surface. Few experiments have achieved measuring Casimir-Polder forces while controlling directly the atom-surface distance. In the range from tens to hundreds of nanometers, the Casimir-Polder potential was measured directly by reflecting the atoms on an evanescent field [10, 11]. In

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

We propose a numerical method to obtain an adequate value for the upper bound on the rank for the tensor completion problem on the variety of third-order tensors of bounded tensor-train rank. The method is inspired by the parametrization of the tangent cone derived by Kutschan (2018). A proof of the adequacy of the upper bound for a related low-rank tensor approximation problem is given and an estimated rank is defined to extend the result to the low-rank tensor completion problem. Some experiments on synthetic data illustrate the approach and show that the method is very robust, e.g., to noise on the data.

Circumstellar disk images have highlighted a wide variety of morphological features. Recovering disk images from high-contrast angular differential imaging (ADI) sequences are however generally affected by geometrical biases, leading to unreliable inference of the morphology of extended disk features. Recently, two types of approaches have been proposed to recover more robust disk images from ADI sequences: iterative principal component analysis, and inverse problem approaches. We introduce MUSTARD, a new IP-based algorithm designed to address the problem of the flux invariant to the rotation in ADI sequences; a limitation inherent to the ADI observing strategy, and discuss the advantages of IP approaches with respect to PCA-based algorithms. The MUSTARD model relies on the addition of morphological priors on the disk and speckle field to a standard IP approach to tackle rotation-invariant signal in circumstellar disk images. We compare the performance of MUSTARD, I-PCA, and standard PCA on a sample of high-contrast imaging data sets acquired in different observing conditions, after injecting a variety of synthetic disk models at different contrast levels. MUSTARD significantly improves the recovery of rotation-invariant signal in disk images, especially for data sets obtained in good observing conditions. However, the MUSTARD model is shown to inadequately handle unstable ADI data sets, and to provide shallower detection limits than PCA-based approaches. MUSTARD has the potential to deliver more robust disk images by introducing a prior to address the inherent ambiguity of ADI observations. However, the effectiveness of the prior is partly hindered by our limited knowledge of the morphological and temporal properties of the stellar speckle halo. In light of this limitation, we suggest that the algorithm could be improved by enforcing a prior based on a library of reference stars, Comment: 35 pages, 38 figures, Accepted for publication in Astronomy & Astrophysics journal the 25 September 2023

The detection and characterization of Earth-like exoplanets around Sun-like stars is a primary science motivation for the Habitable Worlds Observatory. However, the current best technology is not yet advanced enough to reach the 10^-10 contrasts at close angular separations and at the same time remain insensitive to low-order aberrations, as would be required to achieve high-contrast imaging of exo-Earths. Photonic technologies could fill this gap, potentially doubling exo-Earth yield. We review current work on photonic coronagraphs and investigate the potential of hybridized designs which combine both classical coronagraph designs and photonic technologies into a single optical system. We present two possible systems. First, a hybrid solution which splits the field of view spatially such that the photonics handle light within the inner working angle and a conventional coronagraph that suppresses starlight outside it. Second, a hybrid solution where the conventional coronagraph and photonics operate in series, complementing each other and thereby loosening requirements on each subsystem. As photonic technologies continue to advance, a hybrid or fully photonic coronagraph holds great potential for future exoplanet imaging from space., Comment: Conference Proceedings of SPIE: Techniques and Instrumentation for Detection of Exoplanets XI, vol. 12680 (2023)

Metasurfaces offer a flexible framework for the manipulation of light properties in the realm of thin film optics. Specifically, the polarization of light can be effectively controlled through the use of thin phase plates. This study aims to introduce a surrogate optimization framework for these devices. The framework is applied to develop two kinds of vortex phase masks (VPMs) tailored for application in astronomical high-contrast imaging. Computational intelligence techniques are exploited to optimize the geometric features of these devices. The large design space and computational limitations necessitate the use of surrogate models like partial least squares Kriging, radial basis functions, or neural networks. However, we demonstrate the inadequacy of these methods in modeling the performance of VPMs. To address the shortcomings of these methods, a data-efficient evolutionary optimization setup using a deep neural network as a highly accurate and efficient surrogate model is proposed. The optimization process in this study employs a robust particle swarm evolutionary optimization scheme, which operates on explicit geometric parameters of the photonic device. Through this approach, optimal designs are developed for two design candidates. In the most complex case, evolutionary optimization enables optimization of the design that would otherwise be impractical (requiring too much simulations). In both cases, the surrogate model improves the reliability and efficiency of the procedure, effectively reducing the required number of simulations by up to 75% compared to conventional optimization techniques.

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

Recent observations have shown that protoplanetary disks around young stars can embed a wide variety of features. Raw disk images produced by high-contrast imaging instruments are corrupted by slowly varying residual stellar light in the form of quasi-static speckles. Hence, image processing is required to remove speckles from images and to recover circumstellar signals. Current algorithms that rely on the mainstream angular differential imaging (ADI) observing technique are however limited by geometrical biases, and therefore face a major challenge to reliably infer the morphology of extended disk features. In the last two years, four algorithms have been developed for this task, with three of them based on inverse problem (IP) approaches: REXPACO, MAYONNAISE, and \MUSTARD. In this presentation, we will (i) present the new MUSTARD algorithm and (ii) discuss the advantages of IP compared to others methods based on systematic tests., Comment: Presented in ISCS23

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