Search

Exoplanets form from circumstellar protoplanetary discs whose fundamental properties (notably their extent, composition, mass, temperature and lifetime) depend on the host star properties, such as their mass and luminosity. B-stars are among the most massive stars and their protoplanetary discs test extreme conditions for exoplanet formation. This paper investigates the frequency of giant planet companions around young B-stars (median age of 16 Myr) in the Scorpius-Centaurus association, the closest association containing a large population of B-stars. We systematically search for massive exoplanets with the high-contrast direct imaging instrument SPHERE using the data from the BEAST survey, that targets an homogeneous sample of young B-stars from the wide Sco-Cen association. We derive accurate detection limits in case of non-detections. We found evidence in previous papers for two substellar companions around 42 stars. The masses of these companions are straddling the ~13 Jupiter mass deuterium burning limit but their mass ratio with respect to their host star is close to that of Jupiter. We derive a frequency of such massive planetary mass companions around B stars of 11-5+7%, accounting for the survey sensitivity. The discoveries of substellar companions bcen b and mu2sco B happened after only few stars in the survey had been observed, raising the possibility that massive Jovian planets might be common around B-stars. However our statistical analysis show that the occurrence rate of such planets is similar around B-stars and around solar-type stars of similar age, while B-star companions exhibit low mass ratios and larger semi-major axis., Comment: Accepted in A&A

The Roman Space Telescope will be a critical mission to demonstrate high-contrast imaging technologies allowing for the characterisation of exoplanets in reflected light. It will demonstrate $10^{-7}$ contrast limits or better at 3--9 $\lambda / D$ separations with active wavefront control for the first time in space. The detection limits for the Coronagraph Instrument are expected to be set by wavefront variations between the science target and the reference star observations. We are investigating methods to use the deformablel mirrors to methodically probe the impact of such variations on the coronagraphic PSF, generating a PSF library during observations of the reference star to optimise the starlight subtraction at post-processing. We are collaborating with STScI to test and validate these methods in lab using the HiCAT tested, a high-contrast imaging lab platform dedicated to system-level developments for future space missions. In this paper, we will present the first applications of these methods on HiCAT.

The Roman Coronagraph Instrument will be the first space facility equipped with deformable mirrors (DMs). These will lead to reach a contrast of $10^{-8}$ or better in a dark hole between $3-9 \lambda/D$. Post-processing techniques play an important role in increasing the contrast limits. Our work investigates how DMs can be used to calibrate the instrument response to controlled wavefront error maps and to improve the post-processing performance. To this goal, we are developing a simulation pipeline, CAPyBARA, that includes both a propagation model of the Coronagraph and a post-processing module and produces starlight subtracted images of a science target. This pipeline will allow us to investigate alternative observing strategies and test their performance for the Roman Coronagraph. Here we present the first version of the simulator: it currently reproduces the optical propagation, which consists in the hybrid Lyot coronagraph optical structure and dark-hole digging technique (Electric Field Conjugation coupled with $\beta$-bumping), the environment (quasi-static aberration) and the post-processing. With it, we mimic a Coronagraph Instrument observing sequence, which consists in first acquiring reference star data before slewing to the scientific target, and we investigate how the evolution of quasi-static aberrations deteriorate the contrast limit in the dark hole. We simulate a science target with planets at high contrast with their star and we perform a first post-processing analysis with classical subtraction techniques. Here we present the CAPyBARA simulator, as well as some first results. The next step will be to generate PSF libraries by injecting pre-calibrated probes on the DMs (in open loop) during the reference star acquisition and compute a PCA model. Later, we will compare the performance gain obtained with the modulated-DM reference library over standard approaches (RDI).

Context: HR2562B is a planetary-mass companion located 0.56arcsec (19au) from its host star. It is one of a few L/T transitional objects orbiting a young star. This companion provides insight into the evolution of young objects in the L/T transition. However, its key physical properties, such as Teff and mass, remain poorly constrained, with large uncertainties (34% for Teff, 22% for log(g)) based on near-infrared observations alone. Aims: We aim to refine these uncertainties, especially for Teff (1200-1700K) and log(g) (4-5), using new MIR data from the JWST/MIRI filters (10.65, 11.40, and 15.50 microns), and better understand the companion's chemical composition and its role in the L/T transition. Methods: MIRI data were processed using reference star differential imaging, revealing HR2562B at high S/N (16) in all 3 filters. We used 2 atmospheric models, ATMO and ExoREM, to fit the SED, combining MIR and NIR datasets. Additionally, we used CMD with brown dwarfs to explore the chemical composition of HR2562B's atmosphere and compare it to another L/T transition object, VHS1256b. Results: Our analysis improved the temperature precision (Teff=1255+-15K) by 6x compared to previous estimates. We also narrowed its luminosity to -4.69+-0.01 dex. Surface gravity remains uncertain (4.4-4.8), and its mass is estimated between 8 and 18.5Mj, depending on modeling and astrometry. Sensitivity analysis revealed the ability to detect objects between 2-5Mj at 100au. Conclusions: HR2562B likely has a near cloud-free atmosphere, with the ATMO model fitting better than ExoREM. Silicate absorption features are weak, requiring further spectroscopic observations. While HR2562B and VHS1256b share similarities, they are in different evolutionary stages, making HR2562B key to understanding young objects in the L/T transition. It is likely a planetary-mass companion, suggesting a reclassification as HR2562b.

In this study, we explored the fundamental contrast limit of NIRCam coronagraphy observations, representing the achievable performance with post-processing techniques. This limit is influenced by photon noise and readout noise, with complex noise propagation through post-processing methods like principal component analysis. We employed two approaches: developing a formula based on simplified scenarios and using Markov Chain Monte Carlo (MCMC) methods, assuming Gaussian noise properties and uncorrelated pixel noise. Tested on datasets HIP\,65426, AF\,Lep, and HD\,114174, the MCMC method provided accurate but computationally intensive estimates. The analytical approach offered quick, reliable estimates closely matching MCMC results in simpler scenarios. Our findings showed the fundamental contrast curve is significantly deeper than the current achievable contrast limit obtained with post-processing techniques at shorter separations, being 10 times deeper at $0.5''$ and 4 times deeper at $1''$. At greater separations, increased exposure time improves sensitivity, with the transition between photon and readout noise dominance occurring between $2''$ and $3''$. We conclude the analytical approach is a reliable estimate of the fundamental contrast limit, offering a faster alternative to MCMC. These results emphasize the potential for greater sensitivity at shorter separations, highlighting the need for improved or new post-processing techniques to enhance JWST NIRCam sensitivity or contrast curve., Comment: Proceeding of the SPIE July 2024

The search for biosignatures in potentially habitable exoplanets is one of the major astrophysics' drivers for the coming decades, and the prime science goal of the HWO NASA mission, a large UV-Optical-IR space telescope to be launched in the 2040s. To reach this goal, it will be equipped with state-of-the-art high-contrast spectro-imaging capabilities enabling the detection of exoplanets 10^10 times fainter than their host stars, a formidable challenge given today's best detection limits at ~10^-6 contrast levels. This goal puts stringent constraints on the entire observatory, and demands the optimization at the system level to leverage the performance of individual sub-systems. However, while image processing techniques are a key asset to reach the ultimate performance, the science and technological definition of the mission concepts mostly rely on the coronagraph and wavefront control to reject the starlight, assuming a conservative gain of ~10 in sensitivity from image processing, extrapolated from performance obtained with classical techniques on Hubble observations. In the ESCAPE project, we investigate integrated solutions for optimizing the observing methods and data processing techniques with future space telescopes, making use of their wavefront sensors and deformable mirrors. The Roman Space Telescope, scheduled for launch in 2026, will be a critical milestone to demonstrate key technologies ahead of HWO with the Coronagraph instrument, and is thus a unique opportunity to also test and validate innovative image processing techniques. Here we present the rational, methodology, and timeline of the ESCAPE project., Comment: Proceedings of the 2024 SPIE Astronomical Telescopes + Instrumentation conference

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

PLATO (PLAnetary Transits and Oscillations of stars) is ESA's M3 mission designed to detect and characterise extrasolar planets and perform asteroseismic monitoring of a large number of stars. PLATO will detect small planets (down to <2 R_(Earth)) around bright stars (<11 mag), including terrestrial planets in the habitable zone of solar-like stars. With the complement of radial velocity observations from the ground, planets will be characterised for their radius, mass, and age with high accuracy (5 %, 10 %, 10 % for an Earth-Sun combination respectively). PLATO will provide us with a large-scale catalogue of well-characterised small planets up to intermediate orbital periods, relevant for a meaningful comparison to planet formation theories and to better understand planet evolution. It will make possible comparative exoplanetology to place our Solar System planets in a broader context. In parallel, PLATO will study (host) stars using asteroseismology, allowing us to determine the stellar properties with high accuracy, substantially enhancing our knowledge of stellar structure and evolution. The payload instrument consists of 26 cameras with 12cm aperture each. For at least four years, the mission will perform high-precision photometric measurements. Here we review the science objectives, present PLATO's target samples and fields, provide an overview of expected core science performance as well as a description of the instrument and the mission profile at the beginning of the serial production of the flight cameras. PLATO is scheduled for a launch date end 2026. This overview therefore provides a summary of the mission to the community in preparation of the upcoming operational phases.

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

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

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

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

Dual-field interferometric observations with VLTI/GRAVITY sometimes require the use of a "binary calibrator", a binary star whose individual components remain unresolved by the interferometer, with a separation between 400 and 2000 mas for observations with the Units Telescopes (UTs), or 1200 to 3000 mas for the Auxiliary Telescopes (ATs). The separation vector also needs to be predictable to within 10 mas for proper pointing of the instrument. Up until now, no list of properly vetted calibrators was available for dual-field observations with VLTI/GRAVITY on the UTs. Our objective is to compile such a list, and make it available to the community. We identify a list of candidates from the Washington Double Star (WDS) catalogue, all with appropriate separations and brightness, scattered over the Southern sky. We observe them as part of a dedicated calibration programme, and determine whether these objects are true binaries (excluding higher multiplicities resolved interferometrically but unseen by imaging), and extract measurements of the separation vectors. We combine these new measurements with those available in the WDS to determine updated orbital parameters for all our vetted calibrators. We compile a list of 13 vetted binary calibrators for observations with VLTI/GRAVITY on the UTs, and provide orbital estimates and astrometric predictions for each of them. We show that our list guarantees that there are always at least two binary calibrators at airmass < 2 in the sky over the Paranal observatory, at any point in time. Any Principal Investigator wishing to use the dual-field mode of VLTI/GRAVITY with the UTs can now refer to this list to select an appropriate calibrator. We encourage the use of "whereistheplanet" to predict the astrometry of these calibrators, which seamlessly integrates with "p2Gravity" for VLTI/GRAVITY dual-field observing material preparation., Comment: 16 pages, 6 figures. Submitted to A&A

Aims. Measuring dynamical masses of substellar companions is a powerful tool to test models of mass-luminosity-age relations, as well as determining observational features that constrain the boundary between stellar and substellar companions. In order to dynamically constrain the mass of such companions, we use multiple exoplanet measurement techniques to remove degeneracies in the orbital fits of these objects and place tight constraints on their model-independent masses. Methods. We combine long-period radial-velocity data from the CORALIE survey with relative astrometry from direct imaging with VLT/SPHERE, along with astrometric accelerations from Hipparcos-Gaia eDR3 to perform a combined orbital fit and measure precise dynamical masses of two newly discovered benchmark brown dwarfs. Results. We report the discovery of HD112863B and HD206505B, which are two new benchmark likely brown dwarfs that sit at the substellar-stellar boundary, with precise dynamical masses. We perform an orbital fit which yields dynamical masses for HD112863B and HD206505B to be $77.1^{+2.9}_{-2.8}~M_{\rm{Jup}}$ and $79.8\pm1.8~M_{\rm{Jup}}$ respectively. The orbital period for HD112863B is determined to be $21.59\pm0.05$ years and the orbital period of HD206505B is determined to be ${50.9}_{-1.5}^{+1.7}$ years. From the $H$ and $K$ band photometry from IRDIS data taken with VLT/SPHERE, we estimate the spectral types of both HD112863B and HD206505B to be early-mid L-types., Comment: Accepted for publication in A&A, minor language edits added

Aims. We aim to expand the atmospheric exploration of AF Lep b by modeling all available observations obtained with SPHERE at VLT (between 0.95-1.65, at 2.105, and 2.253 $\mu$m, and NIRC2 at Keck (at 3.8 $\mu$m) with self-consistent atmospheric models. Methods. To understand the physical properties of this exoplanet, we used ForMoSA. This forward-modeling code compares observations with grids of pre-computed synthetic atmospheric spectra using Bayesian inference methods. We used Exo-REM, an atmospheric radiative-convective equilibrium model, including the effects of non-equilibrium processes and clouds. Results. From the atmospheric modeling we derive solutions at a low effective temperature of ~750 K. Our analysis also favors a metal-rich atmosphere (>0.4) and solar to super-solar carbon-to-oxygen ratio (~0.6). We tested the robustness of the estimated values for each parameter by cross-validating our models using the leave-one-out strategy, where all points are used iteratively as validation points. Our results indicate that the photometry point at 3.8 $\mu$m strongly drives the metal-rich and super-solar carbon-to-oxygen solutions. Conclusions. Our atmospheric forward-modeling analysis strongly supports the planetary nature of AF Lep b. Its spectral energy distribution is consistent with that of a young, cold, early-T super-Jovian planet. We recover physically consistent solutions for the surface gravity and radius, which allows us to reconcile atmospheric forward modeling with evolutionary models, in agreement with the previously published complementary analysis done by retrievals. Finally, we identified that future data at longer wavelengths are mandatory before concluding about the metal-rich nature of AF Lep b., Comment: 7 pages, 7 figures, 6 tables. Accepted for publication on the 9th of January 2024 in A&A

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

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

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

We present aperture masking interferometry (AMI) observations of the star HIP 65426 at $3.8\,\rm{\mu m}$ as a part of the \textit{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 \textit{JWST} coronagraphs. When combined with \textit{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 contrast of $\Delta m_{F380M}{\sim }7.8$\,mag relative to the host star at a separation of ${\sim}0.07\arcsec$ but detect no additional companions interior to the known companion HIP\,65426\,b. Our observations thus rule out companions more massive than $10{-}12\,\rm{M\textsubscript{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 \textit{JWST} is sensitive to planetary mass companions orbiting at the water frost line, even for more distant stars at $\sim$100\,pc. This result will allow the planning and successful execution of future observations to probe the inner regions of nearby stellar systems, opening essentially unexplored parameter space., Comment: 15 pages, 9 figures, submitted to ApJ Letters

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

The discovery of planets orbiting at less than 1 au from their host star and less massive than Saturn in various exoplanetary systems revolutionized our theories of planetary formation. The fundamental question is whether these close-in low-mass planets could have formed in the inner disk interior to 1 au, or whether they formed further out in the planet-forming disk and migrated inward. Exploring the role of additional giant planets in these systems may help us to pinpoint their global formation and evolution. We searched for additional substellar companions by using direct imaging in systems known to host close-in small planets. The use of direct imaging complemented by radial velocity and astrometric detection limits enabled us to explore the giant planet and brown dwarf demographics around these hosts to investigate the potential connection between both populations. We carried out a direct imaging survey with VLT/SPHERE to look for outer giant planets and brown dwarf companions in 27 systems hosting close-in low-mass planets discovered by radial velocity. Our sample is composed of very nearby (<20pc) planetary systems, orbiting G-, K-, and M-type mature (0.5-10Gyr) stellar hosts. We performed homogeneous direct imaging data reduction and analysis to search for and characterize point sources, and derived robust statistical detection limits. Of 337 point-source detections, we do not find any new bound companions. We recovered the emblematic very cool T-type brown dwarf GJ229B. Our typical sensitivities in direct imaging range from 5 to 30 MJup beyond 2 au. The non-detection of massive companions is consistent with predictions based on models of planet formation by core accretion. Our pilot study opens the way to a multi-technique approach for the exploration of very nearby exoplanetary systems with future ground-based and space observatories., Comment: 49 pages including 31 pages of appendices and references, 31 figures, A&A, accepted

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

A major endeavor of this decade is the direct characterization of young giant exoplanets at high spectral resolution to determine the composition of their atmosphere and infer their formation processes and evolution. Such a goal represents a major challenge owing to their small angular separation and luminosity contrast with respect to their parent stars. Instead of designing and implementing completely new facilities, it has been proposed to leverage the capabilities of existing instruments that offer either high contrast imaging or high dispersion spectroscopy, by coupling them using optical fibers. In this work we present the implementation and first on-sky results of the HiRISE instrument at the very large telescope (VLT), which combines the exoplanet imager SPHERE with the recently upgraded high resolution spectrograph CRIRES using single-mode fibers. The goal of HiRISE is to enable the characterization of known companions in the $H$ band, at a spectral resolution of the order of $R = \lambda/\Delta\lambda = 100\,000$, in a few hours of observing time. We present the main design choices and the technical implementation of the system, which is constituted of three major parts: the fiber injection module inside of SPHERE, the fiber bundle around the telescope, and the fiber extraction module at the entrance of CRIRES. We also detail the specific calibrations required for HiRISE and the operations of the instrument for science observations. Finally, we detail the performance of the system in terms of astrometry, temporal stability, optical aberrations, and transmission, for which we report a peak value of $\sim$3.9% based on sky measurements in median observing conditions. Finally, we report on the first astrophysical detection of HiRISE to illustrate its potential., Comment: 17 pages, 15 figures, 3 tables. Accepted in A&A on 24 October 2023

We propose to explore a cascade extreme Adaptive optics (ExAO) approach with a second stage based on a Zernike wavefront sensor (ZWFS) for exoplanet imaging and spectroscopy. Most exoplanet imagers currently use a single-stage ExAO to correct for the effects of atmospheric turbulence and produce high-Strehl images of observed stars in the near-infrared. While such systems enable the observation of warm gaseous companions around nearby stars, adding a second-stage AO enables to push the wavefront correction further and possibly observe colder or smaller planets. This approach is currently investigated in different exoplanet imagers (VLT/SPHERE, Mag-AOX, Subaru/SCExAO) by considering a Pyramid wavefront sensor (PWFS) in the second arm to measure the residual atmospheric turbulence left from the first stage. Since these aberrations are expected to be very small (a few tens of nm in the near-infrared domain), we propose to investigate an alternative approach based on the ZWFS. This sensor is a promising concept with a small capture range to estimate residual wavefront errors thanks to its large sensitivity, simple phase reconstruction and easiness of implementation. In this contribution, we perform preliminary tests on the GHOST testbed at ESO to validate this approach experimentally. Additional experiments with petalling effects are also showed, giving promising wavefront correction results. Finally, we briefly discuss a first comparison between PWFS-based and ZWFS-based second-stage AO to draw preliminary conclusions on the interests of both schemes for exoplanet imaging and spectroscopy with the upgrade of the current exoplanet imagers and the envisioned ExAO instruments for ELTs., Comment: 17 pages, 10 figures, pre-print of the proceeding of the AO4ELT7 conference held in June 2023 in Avignon, France

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

An edge-on debris disk was detected in 2015 around the young, nearby A0V star HD 110058. The disk showed features resembling those seen in the disk of beta Pictoris that could indicate the presence of a perturbing planetary-mass companion in the system. We investigated new and archival scattered light images of the disk in order to characterise its morphology and spectrum. In particular, we analysed the disk's warp to constrain the properties of possible planetary perturbers. Our work uses data from two VLT/SPHERE observations and archival data from HST/STIS. We measured the morphology of the disk by analysing vertical profiles along the length of the disk to extract the centroid spine position and vertical height. We extracted the surface brightness and reflectance spectrum of the disk. We detect the disk between 20 au (with SPHERE) and 150 au (with STIS), at a position angle of 159.6$^\circ\pm$0.6$^\circ$. Analysis of the spine shows an asymmetry between the two sides of the disk, with a 3.4$^\circ\pm$0.9$^\circ$ warp between ~20 au and 60 au. The disk is marginally vertically resolved in scattered light, with a vertical aspect ratio of 9.3$\pm$0.7% at 45 au. The extracted reflectance spectrum is featureless, flat between 0.95 micron and 1.1 micron, and red from 1.1 micron to 1.65 micron. The outer parts of the disk are also asymmetric with a tilt between the two sides compatible with a disk made of forward-scattering particles and seen not perfectly edge-on, suggesting an inclination of <84$^\circ$. The presence of an undetected planetary-mass companion on an inclined orbit with respect to the disk could explain the warp. The misalignment of the inner parts of the disk with respect to the outer disk suggests a warp that has not yet propagated to the outer parts of the disk, favouring the scenario of an inner perturber as the origin of the warp., Comment: 17 pages, 15 figures, 3 tables; accepted for publication in A&A

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

Radio and near-infrared observations have observed dozens of protoplanetary disks that host spiral arm features. Numerical simulations have shown that companions may excite spiral density waves in protoplanetary disks via companion-disk interaction. However, the lack of direct observational evidence for spiral-driving companions poses challenges to current theories of companion-disk interaction. Here we report multi-epoch observations of the binary system HD 100453 with the Spectro-Polarimetric High-contrast Exoplanet REsearch (SPHERE) facility at the Very Large Telescope. By recovering the spiral features via robustly removing starlight contamination, we measure spiral motion across 4 yr to perform dynamical motion analyses. The spiral pattern motion is consistent with the orbital motion of the eccentric companion. With this first observational evidence of a companion driving a spiral arm among protoplanetary disks, we directly and dynamically confirm the long-standing theory on the origin of spiral features in protoplanetary disks. With the pattern motion of companion-driven spirals being independent of companion mass, here we establish a feasible way of searching for hidden spiral-arm-driving planets that are beyond the detection of existing ground-based high-contrast imagers., Comment: Accepted for publication in Astronomy and Astrophysics; 12 pages, 9 figures

Recent observations from B-star Exoplanet Abundance Study (BEAST) have illustrated the existence of sub-stellar companions around very massive stars. In this paper, we present the detection of two lower mass companions to a relatively nearby ($148.7^{+1.5}_{-1.3}$ pc), young ($17^{+3}_{-4}$ Myr), bright (V=$6.632\pm0.006$ mag), $2.58\pm0.06~ M_{\odot}$ B9V star HIP 81208 residing in the Sco-Cen association, using the Spectro-Polarimetric High-contrast Exoplanet REsearch (SPHERE) instrument at the Very Large Telescope (VLT) in Chile. Analysis of the photometry obtained gives mass estimates of $67^{+6}_{-7}~M_J$ for the inner companion and $0.135^{+0.010}_{-0.013}~M_{\odot}$ for the outer companion, indicating the former to be most likely a brown dwarf and the latter to be a low-mass star. The system is compact but unusual, as the orbital planes of the two companions are likely close to orthogonal. The preliminary orbital solutions we derived for the system indicate that the star and the two companions are likely in a Kozai resonance, rendering the system dynamically very interesting for future studies., Comment: 18 pages, 14 figures, 5 tables Accepted for publication in the 10. Planets and planetary systems section of A&A

Since 2005 ESO has been working with its community and industry to develop an extremely large optical/infrared telescope. ESO's Extremely Large Telescope, or ELT for short, is a revolutionary ground-based telescope that will have a 39-metre main mirror and will be the largest visible and infrared light telescope in the world. To address specific topics that are needed for the science operations and calibrations of the telescope, thirteen specific working groups were created to coordinate the effort between ESO, the instrument consortia, and the wider community. We describe here the goals of these working groups as well as their achievements so far., Comment: 8 pages, 4 figures

Wavefront sensing and control (WFSC) will play a key role in improving the stability of future large segmented space telescopes while relaxing the thermo-mechanical constraints on the observatory structure. Coupled with a coronagraph to reject the light of an observed bright star, WFSC enables the generation and stabilisation of a dark hole (DH) in the star image to perform planet observations. While WFSC traditionally relies on a single wavefront sensor (WFS) input to measure wavefront errors, the next generation of instruments will require several WFSs to address aberrations with different sets of spatial and temporal frequency contents. The multiple measurements produced in such a way will then have to be combined and converted to commands for deformable mirrors (DMs) to modify the wavefront subsequently. We asynchronously operate a loop controlling the high-order modes digging a DH and a control loop that uses the rejected light by a Lyot coronagraph with a Zernike wavefront sensor to stabilize the low-order aberrations. Using the HiCAT testbed with a segmented telescope aperture, we implement concurrent operations and quantify the expected cross-talk between the two controllers. We then present experiments that alternate high-order and low-order control loops to identify and estimate their respective contributions. We show an efficient combination of the high-order and low-order control loops, keeping a DH contrast better than 5 x 10-8 over a 30 min experiment and stability improvement by a factor of 1.5. In particular, we show a contrast gain of 1.5 at separations close to the DH inner working angle, thanks to the low-order controller contribution. Concurrently digging a DH and using the light rejected by a Lyot coronagraph to stabilize the wavefront is a promising path towards exoplanet imaging and spectroscopy with future large space observatories.

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

We present new optical and near-IR images of debris disk around the F-type star HD 114082. We obtained direct imaging observations and analysed the TESS photometric time series data of this target with a goal to search for planetary companions and to characterise the morphology of the debris disk and the scattering properties of dust particles. HD 114082 was observed with the VLT/SPHERE instrument: the IRDIS camera in the K band together with the IFS in the Y, J and H band using the ADI technique as well as IRDIS in the H band and ZIMPOL in the I_PRIME band using the PDI technique. The scattered light images were fitted with a 3D model for single scattering in an optically thin dust disk. We performed aperture photometry in order to derive the scattering and polarized phase functions, polarization fraction and spectral scattering albedo for the dust particles in the disk. This method was also used to obtain the reflectance spectrum of the disk to retrieve the disk color and study the dust reflectivity in comparison to the debris disk HD 117214. We also performed the modeling of the HD 114082 light curve measured by TESS using the models for planet transit and stellar activity to put constraints on radius of the detected planet and its orbit. The debris disk appears as an axisymmetric debris belt with a radius of ~0.37$"$ (35 au), inclination of ~83$^\circ$ and a wide inner cavity. Dust particles in HD 114082 have a maximum polarization fraction of ~17% and a high reflectivity which results in a spectral scattering albedo of 0.65. The analysis of TESS photometric data reveals a transiting planetary companion to HD 114082 with a radius of $\sim$1~$\rm R_{J}$ on an orbit with a semi-major axis of $0.7 \pm 0.4$ au. Combining different data, we reach deep sensitivity limits in terms of companion masses down to ~5$M_{\rm Jup}$ at 50 au, and ~10 $M_{\rm Jup}$ at 30 au from the central star., Comment: 27 pages

We aim to revisit the system AB Pic which has a known companion at the exoplanet/ brown-dwarf boundary. We based this study on a rich set of observations to investigate the companion's orbit and atmosphere. We composed a spectrum of AB Pic b merging archival VLT/SINFONI K-band data, with published spectra at J and H-band (SINFONI) and Lp-band (Magellan-AO), and photometric measurements (HST and Spitzer). We modeled the spectrum with ForMoSA, based on two atmospheric models: ExoREM and BT-SETTL13. We determined the orbital properties of b fitting the astrometric measurements from NaCo (2003 and 2004) and SPHERE (2015). The orbital solutions favor a semi-major axis of $\sim$190au viewed edge-on. With Exo-REM, we derive a T$_{eff}$ of 1700$\pm$50K and surface gravity of 4.5$\pm$0.3dex, consistent with previous works, and we report for the first time a C/O ratio of 0.58$\pm$0.08 ($\sim$solar). The posteriors are sensitive to the wavelength interval and the family of models used. Given the 2.1hr rotation period and our vsin(i) of $\sim$73km/s, we estimate for the first time the true obliquity to be $\sim$45 or $\sim$135deg, indicating a significant misalignment between the planet's spin and orbit orientations. Finally, a proper motion anomaly between the Hipparcos and Gaia eDR3 compared to our SPHERE detection limits and adapted radial velocity limits indicate the existence of a $\sim$6M$_{Jup}$ inner planet orbiting from 2 to 10au (40-200mas). The possible existence of an inner companion and the likely miss-alignment of the spin axis orientation strongly favor a formation path by gravitational instability or core accretion within a disk closer inside followed by dynamical interactions. Confirmation and characterization of planet c and access to a broader wavelength coverage for planet b will be essential to probe the uncertainties associated with the parameters., Comment: 17 pages, 13 Figures, 6 Tables

We present a new open-source data-reduction pipeline to reconstruct spectral data cubes from raw SPHERE integral-field spectrograph (IFS) data. The pipeline is written in Python and based on the pipeline that was developed for the CHARIS IFS. It introduces several improvements to SPHERE data analysis that ultimately produce significant improvements in postprocessing sensitivity. We first used new data to measure SPHERE lenslet point spread functions (PSFs) at the four laser calibration wavelengths. These lenslet PSFs enabled us to forward-model SPHERE data, to extract spectra using a least-squares fit, and to remove spectral crosstalk using the measured lenslet PSFs. Our approach also reduces the number of required interpolations, both spectral and spatial, and can preserve the original hexagonal lenslet geometry in the SPHERE IFS. In the case of least-squares extraction, no interpolation of the data is performed. We demonstrate this new pipeline on the directly imaged exoplanet 51 Eri b and on observations of the hot white dwarf companion to HD 2133. The extracted spectrum of HD 2133B matches theoretical models, demonstrating spectrophotometric calibration that is good to a few percent. Postprocessing on two 51 Eri b data sets demonstrates a median improvement in sensitivity of 80% and 30% for the 2015 and 2017 data, respectively, compared to the use of cubes reconstructed by the SPHERE Data Center. The largest improvements are seen for poorer observing conditions. The new SPHERE pipeline takes less than three minutes to produce a data cube on a modern laptop, making it practical to reprocess all SPHERE IFS data., Comment: 17 pages, 11 figures. Software available at: https://github.com/PrincetonUniversity/charis-dep

SPHERE+ is a proposed upgrade of the SPHERE instrument at the VLT, which is intended to boost the current performances of detection and characterization for exoplanets and disks. SPHERE+ will also serve as a demonstrator for the future planet finder (PCS) of the European ELT. The main science drivers for SPHERE+ are 1/ to access the bulk of the young giant planet population down to the snow line ($3-10$ au), to bridge the gap with complementary techniques (radial velocity, astrometry); 2/ to observe fainter and redder targets in the youngest ($1-10$\,Myr) associations compared to those observed with SPHERE to directly study the formation of giant planets in their birth environment; 3/ to improve the level of characterization of exoplanetary atmospheres by increasing the spectral resolution in order to break degeneracies in giant planet atmosphere models. Achieving these objectives requires to increase the bandwidth of the xAO system (from $\sim$1 to 3\,kHz) as well as the sensitivity in the infrared (2 to 3\,mag). These features will be brought by a second stage AO system optimized in the infrared with a pyramid wavefront sensor. As a new science instrument, a medium resolution integral field spectrograph will provide a spectral resolution from 1000 to 5000 in the J and H bands. This paper gives an overview of the science drivers, requirements and key instrumental trade-off that were done for SPHERE+ to reach the final selected baseline concept., Comment: To appear in the Proceedings of the SPIE Astronomical Telescopes + Instrumentation (2022), 13 pages, 6 figure

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

Young M-type binaries are particularly useful for precise isochronal dating by taking advantage of their extended pre-main sequence evolution. Orbital monitoring of these low-mass objects becomes essential in constraining their fundamental properties, as dynamical masses can be extracted from their Keplerian motion. Here, we present the combined efforts of the AstraLux Large Multiplicity Survey, together with a filler sub-programme from the SpHere INfrared Exoplanet (SHINE) project and previously unpublished data from the FastCam lucky imaging camera at the Nordical Optical Telescope (NOT) and the NaCo instrument at the Very Large Telescope (VLT). Building on previous work, we use archival and new astrometric data to constrain orbital parameters for 20 M-type binaries. We identify that eight of the binaries have strong Bayesian probabilities and belong to known young moving groups (YMGs). We provide a first attempt at constraining orbital parameters for 14 of the binaries in our sample, with the remaining six having previously fitted orbits for which we provide additional astrometric data and updated Gaia parallaxes. The substantial orbital information built up here for four of the binaries allows for direct comparison between individual dynamical masses and theoretical masses from stellar evolutionary model isochrones, with an additional three binary systems with tentative individual dynamical mass estimates likely to be improved in the near future. We attained an overall agreement between the dynamical masses and the theoretical masses from the isochrones based on the assumed YMG age of the respective binary pair. The two systems with the best orbital constrains for which we obtained individual dynamical masses, J0728 and J2317, display higher dynamical masses than predicted by evolutionary models., Comment: Accepted for publication in A&A. 38 pages, 29 figures, 11 tables

Reference-star differential imaging (RDI) is a promising technique in high-contrast imaging that is thought to be more sensitive to exoplanets and disks than angular differential imaging (ADI) at short angular separations (i.e., <0.3"). However, it is unknown whether the performance of RDI on ground-based instruments can be improved by using all the archival data to optimize the subtraction of stellar contributions. We characterize the performance of RDI on SPHERE/IRDIS data in direct imaging of exoplanets and disks. We made use of all the archival data in H23 obtained by SPHERE/IRDIS in the past five years to build a master reference library and perform RDI. In the point-source detection, RDI can outperform ADI at small angular separations (<0.4") if the observing conditions are around the median conditions of our master reference library. On average, RDI has a gain of ~0.8 mag over ADI at 0.15" separation for observations under median conditions. We demonstrate that including more reference targets in the master reference library can indeed help to improve the performance of RDI. In disk imaging, RDI can reveal more disk features and provide a more robust recovery of the disk morphology. We resolve 33 disks in total intensity (19 planet-forming disks and 14 debris disks), and 4 of them can only be detected with RDI. Two disks are resolved in scattered light for the first time. Three disks are detected in total intensity for the first time. The master reference library we built in this work can be easily implemented into legacy or future SPHERE surveys to perform RDI, achieving better performance than that of ADI. To obtain optimal RDI gains over ADI, we recommend future observations be carried out under seeing conditions of 0.6"-0.8"., Comment: Accepted for publication in Astronomy and Astrophysics; 19 pages, 14 figures, 3 tables

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

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

HR\,8799 is a young planetary system composed of 4 planets and a double debris belt. Being the first multi-planetary system discovered with the direct imaging technique, it has been observed extensively since 1998. This wide baseline of astrometric measurements, counting over 50 observations in 20 years, permits a detailed orbital and dynamical analysis of the system. To explore the orbital parameters of the planets, their dynamical history, and the planet-to-disk interaction, we made follow-up observations of the system during the VLT/SPHERE GTO program. We obtained 21 observations, most of them in favorable conditions. In addition, we observed HR\,8799 with the instrument LBT/LUCI. All the observations were reduced with state-of-the-art algorithms implemented to apply the spectral and angular differential imaging method. We re-reduced the SPHERE data obtained during the commissioning of the instrument and in 3 open-time programs to have homogeneous astrometry. The precise position of the 4 planets with respect to the host star was calculated by exploiting the fake negative companions method. To improve the orbital fitting, we also took into account all of the astrometric data available in the literature. From the photometric measurements obtained in different wavelengths, we estimated the planets' masses following the evolutionary models. We obtained updated parameters for the orbits with the assumption of coplanarity, relatively small eccentricities, and periods very close to the 2:1 resonance. We also refined the dynamical mass of each planet and the parallax of the system (24.49 $\pm$ 0.07 mas). We also conducted detailed $N$-body simulations indicating possible positions of a~putative fifth innermost planet with a mass below the present detection limits of $\simeq 3$~\MJup., Comment: 27 pages, 15 figures, A&A in press