Search

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

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

Molecular mapping is a supervised method exploiting the spectral diversity of integral field spectrographs to detect and characterize resolved exoplanets blurred into the stellar halo. We present an evolution of the method to remove the stellar halo and the nuisance of telluric features in the datacubes and access a continuum-subtracted spectra of the planets at R$\sim$4000. We derive planet atmosphere properties from a direct analysis of the planet telluric-corrected absorption spectrum. We applied our methods to the SINFONI observation of the planet $\beta$ Pictoris b. We recover the CO and H$_2$O detections in the atmosphere of $\beta$ Pic b using molecular mapping. We further determine some basic properties of its atmosphere, with $T_\text{eq}=1748^{+3}_{-4}$ K, a sub-solar [Fe/H]=$-0.235^{+0.015}_{-0.013}$ dex, and a solar C/O=$0.551 \pm 0.002$ in contrast with values measured for the same exoplanet with other infrared instruments. We confirm a low projected equatorial velocity of 25$^{+5}_{-6}$ km s$^{-1}$. We are also able to measure, for the first time with a medium-resolution spectrograph, the radial velocity of $\beta$ Pic b relative to the central star at MJD=56910.38 with a km/s precision of $-11.3 \pm 1.1$ km s$^{-1}$, compatible with ephemerides based on the current knowledge of the $\beta$ Pic system., Comment: 20 pages, 25 figures, 5 tables, accepted for publication in A&A

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. 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

For decades, the spectral variations of Beta Pictoris have been modelled as the result of the evaporation of exocomets close to the star, termed falling evaporating bodies (FEBs). Resonant perturbations by a giant planet have been proposed to explain the dynamical origin of these stargrazers. The disk is now known to harbour two giant planets, Beta Pic b and c, orbiting the star at 9.9 au and 2.7 au. While the former almost matches the planet formerly suspected, the discovery of the latter complicates the picture. We first question the stability of the two-planet system. Then we investigate the dynamics of a disk of planetesimals orbiting the star with both planets to check the validity of the FEB generation mechanism. Symplectic N-body simulations are used to determine which regions of the planetesimal disk are dynamically stable. Then we focus on regions where disk particles are able to reach high eccentricities thanks to resonant mechanisms. The first result is that the system is dynamically stable. Both planets may temporarily fall in 7:1 mean motion resonance (MMR). Then, simulations reveal that the whole region extending between ~1.5 au and ~25 au is unstable to planetary perturbations. However, a disk below 1.5 au survives, which appears to constitute an active source of FEBs via high-order MMRs with Beta Pic c. Beta Pic b acts as a distant perturber that helps sustain the whole process. These simulations rule out the preceding FEB generation mechanism model, which placed their origin at around 4-5 au. Conversely, FEBs are likely to originate from a region much further in and related to MMRs with Beta Pic c. That mechanism also appears to last longer, as new planetesimals are able to continuously enter the MMRs and evolve towards the FEB state. Subsequently, the physical nature of the FEBs may differ from that previously thought, and presumably may not be icy., Comment: 14 pages, 10 figures. To appear in Astronomy & Astrophysics

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

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

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

Since 2019, the direct imaging B-star Exoplanet Abundance Study (BEAST) at SPHERE@VLT has been scanning the surroundings of young B-type stars in order to ascertain the ultimate frontiers of giant planet formation. Recently, the $17^{+3}_{-4}$ Myr HIP 81208 was found to host a close-in (~50 au) brown dwarf and a wider (~230 au) late M star around the central 2.6Msun primary. Alongside the continuation of the survey, we are undertaking a complete reanalysis of archival data aimed at improving detection performances so as to uncover additional low-mass companions. We present here a new reduction of the observations of HIP 81208 using PACO ASDI, a recent and powerful algorithm dedicated to processing high-contrast imaging datasets, as well as more classical algorithms and a dedicated PSF-subtraction approach. The combination of different techniques allowed for a reliable extraction of astrometric and photometric parameters. A previously undetected source was recovered at a short separation from the C component of the system. Proper motion analysis provided robust evidence for the gravitational bond of the object to HIP 81208 C. Orbiting C at a distance of ~20 au, this 15Mjup brown dwarf becomes the fourth object of the hierarchical HIP 81208 system. Among the several BEAST stars which are being found to host substellar companions, HIP 81208 stands out as a particularly striking system. As the first stellar binary system with substellar companions around each component ever found by direct imaging, it yields exquisite opportunities for thorough formation and dynamical follow-up studies., Comment: 12 pages, 9 figures, 5 tables. Accepted for publication as a Letter in Astronomy and Astrophysics, section 1. Letters to the Editor

[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

Orbital motion in binary and planetary systems is the main source of precise stellar and planetary mass measurements, and joint analysis of data from multiple observational methods can both lift degeneracies and improve precision. We set out to measure the masses of individual stars in binary systems using all the information brought by the Hipparcos and Gaia absolute astrometric missions. We present BINARYS, a tool which uses the Hipparcos and Gaia absolute astrometric data and combines it with relative astrometry and/or radial velocity measurements to determine the orbit of a binary system. It rigorously combines the Hipparcos and Gaia data (here EDR3), and it can use the Hipparcos Transit Data as needed for binaries where Hipparcos detect significant flux from the secondary component. It also support the case where Gaia resolved the system, giving an astrometric solution for both components. We determine model-independent individual masses for the first time for three systems: the two mature binaries Gl~494 ($M_1=0.584 \pm 0.003 M_{\odot}$ and $M_2=87 \pm 1 M_{\textrm{Jup}}$) and HIP~88745 ($M_1=0.96 \pm 0.02 M_{\odot}$ and $M_2= 0.60^{+ 0.02 }_{- 0.01 } M_{\odot}$), and the younger AB Dor member GJ~2060 ($M_1=0.60 ^{+ 0.06}_{- 0.05} M_{\odot}$ and $M_2=0.45 ^{+ 0.06}_{- 0.05}M_{\odot}$). The latter provides a rare test of evolutionary model predictions at young ages in the low stellar-mass range and sets a lower age limit of 100~Myr for the moving group.

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

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

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

Aims. This work aims at constraining the masses and separations of potential substellar companions to five accelerating stars (HIP 1481, HIP 88399, HIP 96334, HIP 30314 and HIP 116063) using multiple data sets acquired with different techniques. Methods. Our targets were originally observed as part of the SPHERE/SHINE survey, and radial velocity (RV) archive data were also available for four of the five objects. No companions were originally detected in any of these data sets, but the presence of significant proper motion anomalies (PMa) for all the stars strongly suggested the presence of a companion. Combining the information from the PMa with the limits derived from the RV and SPHERE data, we were able to put constraints on the characteristics of the unseen companions. Results. Our analysis led to relatively strong constraints for both HIP 1481 and HIP 88399, narrowing down the companion masses to 2-5 M_Jup and 3-5 M_Jup and separations within 2-15 au and 3-9 au, respectively. Because of the large age uncertainties for HIP 96334, the poor observing conditions for the SPHERE epochs of HIP 30314 and the lack of RV data for HIP 116063, the results for these targets were not as well defined, but we were still able to constrain the properties of the putative companions within a reasonable confidence level. Conclusions. For all five targets, our analysis has revealed that the companions responsible for the PMa signal would be well within reach for future instruments planned for the ELT (e.g., MICADO), which would easily achieve the required contrast and angular resolution. Our results therefore represent yet another confirmation of the power of multi-technique approaches for both the discovery and characterisation of planetary systems., Comment: 11 pages, 7 figures, accepted for publications on A&A

Context. HD 95086 is a young nearby Solar System analog hosting a giant exoplanet orbiting at 57 au from the star between an inner and outer debris belt. The existence of additional planets has been suggested as the mechanism that maintains the broad cavity between the two belts. Aims. We present a dedicated monitoring of HD 95086 with the VLT/SPHERE instrument to refine the orbital and atmospheric properties of HD 95086 b, and to search for additional planets in this system. Methods. SPHERE observations, spread over ten epochs from 2015 to 2019 and including five new datasets, were used. Combined with archival observations, from VLT/NaCo (2012-2013) and Gemini/GPI (2013-2016), the extended set of astrometric measurements allowed us to refine the orbital properties of HD 95086 b. We also investigated the spectral properties and the presence of a circumplanetary disk around HD 95086 b by using the special fitting tool exploring the diversity of several atmospheric models. In addition, we improved our detection limits in order to search for a putative planet c via the K-Stacker algorithm. Results. We extracted for the first time the JH low-resolution spectrum of HD 95086 b by stacking the six best epochs, and confirm its very red spectral energy distribution. Combined with additional datasets from GPI and NaCo, our analysis indicates that this very red color can be explained by the presence of a circumplanetary disk around planet b, with a range of high-temperature solutions (1400-1600 K) and significant extinction (Av > 10 mag), or by a super-solar metallicity atmosphere with lower temperatures (800-1300 K), and small to medium amount of extinction (Av < 10 mag). We do not find any robust candidates for planet c, but give updated constraints on its potential mass and location., Comment: 29 pages, 20 figures, A&A, accepted

Virtually all known exoplanets reside around stars with $M<2.3~M_\odot$; to clarify if the dearth of planets around more massive stars is real, we launched the direct-imaging B-star Exoplanet Abundance STudy (BEAST) survey targeting B stars ($M>2.4~M_\odot$) in the young (5-20 Myr) Scorpius-Centaurus association (Sco-Cen). Here we present the case of a massive ($M \sim 9~M_\odot$) BEAST target, $\mu^2$ Sco. Based on kinematic information, we found that $\mu^2$ Sco is a member of a small group which we label Eastern Lower Scorpius, refining in turn the precision on stellar parameters. Around this star we identified a robustly detected substellar companion ($14.4\pm 0.8 M_J$) at a projected separation of $290\pm 10$ au, and a probable second object ($18.5\pm 1.5 M_J$) at $21\pm 1$ au. The planet-to-star mass ratios of these objects are similar to that of Jupiter to the Sun, and their irradiation is similar to those of Jupiter and Mercury, respectively. The two companions of $\mu^2$ Sco are naturally added to the giant planet b Cen b recently discovered by BEAST; although slightly more massive than the deuterium burning limit, their properties resemble those of giant planets around less massive stars and they are better reproduced by a formation under a planet-like, rather than a star-like scenario. Irrespective of the (needed) confirmation of the inner companion, $\mu^2$ Sco is the first star that would end its life as a supernova that hosts such a system. The tentative high frequency of BEAST discoveries shows that giant planets or small-mass brown dwarfs can form around B stars. When putting this finding in the context of core accretion and gravitational instability, we conclude that the current modeling of both mechanisms is not able to produce this kind of companion. BEAST will pave the way for the first time to an extension of these models to intermediate and massive stars. (abridged), Comment: 26 pages, 19 figures. Accepted for publication in A&A

We present near-infrared Large Binocular Telescope Interferometer LMIRCam imagery of the disk around the Herbig Ae/Be star AB Aurigae. A comparison of surface brightness at Ks (2.16 ${\mu}$m), H2O narrowband (3.08 ${\mu}$m), and L' (3.7 ${\mu}$m) allows us to probe the presence of icy grains in this (pre)transitional disk environment. By applying Reference Differential Imaging PSF subtraction, we detect the disk at high signal to noise in all three bands. We find strong morphological differences between bands, including asymmetries consistent with observed spiral arms within 100 AU in L'. An apparent deficit of scattered light at 3.08 ${\mu}$m relative to bracketing wavelengths (Ks and L') is evocative of ice absorption at the disk surface layer. However, the ${\Delta}$(Ks-H2O) color is consistent with grains with little to no ice (0-5% by mass). The ${\Delta}$(H2O-L') color, conversely, suggests grains with a much higher ice mass fraction (~0.68), and the two colors cannot be reconciled under a single grain population model. Additionally, we find the extremely red ${\Delta}$(Ks-L') disk color cannot be reproduced under conventional scattered light modeling with any combination of grain parameters or reasonable local extinction values. We hypothesize that the scattering surfaces at the three wavelengths are not co-located, and optical depth effects result in each wavelength probing the grain population at different disk surface depths. The morphological similarity between Ks and H2O suggests their scattering surfaces are near one another, lending credence to the ${\Delta}$(Ks-H2O) disk color constraint of < 5% ice mass fraction for the outermost scattering disk layer., Comment: Accepted for publication in AJ

To understand how the multitude of planetary systems that have been discovered come to be, we need to study systems at different evolutionary stages, with different central stars but also in different environments. The most challenging environment for planet formation may be the harsh UV radiation field of nearby massive stars which quickly erodes disks by external photo-evaporation. We have observed the AT Pyx system, located in the head of a cometary globule in the Gum Nebula, to search for signs of ongoing planet formation. We used the extreme adaptive optics imager VLT/SPHERE to observe AT Pyx in polarized light as well as total intensity in the J, H and K-band. Additionally we employed VLT/NACO to observe the system in the L-band. We resolve the disk around AT Pyx in scattered light across multiple wavelengths. We find an extended (>126 au) disk, with an intermediate inclination between 35 deg and 42 deg. The disk shows complex sub-structure and we identify 2 and possibly 3 spiral-like features. Depending on the precise geometry of the disk (which we can not unambiguously infer from our data) the disk may be eccentric with an eccentricity of ~0.16 or partially self-shadowed. The spiral features and possible eccentricity are both consistent with signatures of an embedded gas giant planet equal in mass to Jupiter. Our own observations can rule out brown dwarf companions embedded in the resolved disk, but are not sensitive enough to detect gas giants. AT Pyx is the first disk in a cometray globule in the Gum Nebula which is spatially resolved. By comparison with disks in the Orion Nebula Cluster we note that the extension of the disk may be exceptional for this environment if the external UV radiation field is comparable to other cometary globules in the region. The signposts of ongoing planet formation are intriguing and need to be followed up with higher sensitivity., Comment: 11 pages, 9 figures, accepted for publication in A&A

We initiated a deep-imaging survey of Scorpius-Centaurus A-F stars with predicted warm inner and cold outer belts of debris reminiscent of the architecture of emblematic systems such as HR 8799. We present resolved SPHERE images of a narrow ring of debris around the F5-type star HD 141011 that was observed as part of our survey in 2015, 2016, and 2019. The ring extends up to ~1.1" (~141 au) from the star in the IRDIS and IFS data obtained in 2016 and 2019. The disk is not detected in the 2015 data which are of poorer quality. The disks is best reproduced by models of a noneccentric ring centered on the star with an inclination of $69.1\pm0.9^{\circ}$, a position angle of $-24.6 \pm 1.7^{\circ}$, and a semimajor axis of $127.5\pm3.8$ au. The combination of radial velocity and imaging data excludes brown-dwarf (M>13.6 MJup) companions coplanar with the disk from 0.1 to 0.9 au and from 20 au up to 500 au (90% probability). HD 141011 adds to the growing list of debris disks that are resolved in Sco-Cen. It is one of the faintest disks that are resolved from the ground and has a radial extent and fractional width ($\sim$12.5%) reminiscent of Fomalhaut. Its moderate inclination and large semimajor axis make it a good target for the James Webb Space Telescope and should allow a deeper search for putative companions shaping the dust distribution., Comment: 13 pages, 10 figures. Paper II of our Sco-Cen imaging survey. Accepted in A&A

Context. Protoplanetary disks around young stars often contain substructures like rings, gaps, and spirals that could be caused by interactions between the disk and forming planets. Aims. We aim to study the young (1-3 Myr) star DR Tau in the near-infrared and characterize its disk, which was previously resolved through sub-millimeter interferometry with ALMA, and to search for possible sub-stellar companions embedded into it. Methods. We observed DR Tau with VLT/SPHERE both in polarized light (H broad band) and total intensity (in Y, J, H, and K spectral bands). We also performed L' band observations with LBTI/LMIRCam on the Large Binocular Telescope (LBT). Results. We found two previously undetected spirals extending north-east and south of the star, respectively. We further detected an arc-like structure north of the star. Finally a bright, compact and elongated structure was detected at separation of 303 +/- 10 mas and position angle 21.2 +/- 3.7 degrees, just at the root of the north-east spiral arm. Since this feature is visible both in polarized light and in total intensity and has a flat spectrum it is likely caused by stellar light scattered by dust. Conclusions. The two spiral arms are at different separation from the star, have very different pitch angles, and are separated by an apparent discontinuity, suggesting they might have a different origin. The very open southern spiral arm might be caused by infalling material from late encounters with cloudlets into the formation environment of the star itself. The compact feature could be caused by interaction with a planet in formation still embedded in its dust envelope and it could be responsible for launching the north-east spiral. We estimate a mass of the putative embedded object of the order of few M_Jup ., Comment: 13 pages, 11 Figures. Accepted for publications on A&A

We aim to demonstrate that the presence and mass of an exoplanet can now be effectively derived from the astrometry of another exoplanet. We combined previous astrometry of $\beta$ Pictoris b with a new set of observations from the GRAVITY interferometer. The orbital motion of $\beta$ Pictoris b is fit using Markov chain Monte Carlo simulations in Jacobi coordinates. The inner planet, $\beta$ Pictoris c, was also reobserved at a separation of 96\,mas, confirming the previous orbital estimations. From the astrometry of planet b only, we can (i) detect the presence of $\beta$ Pictoris c and (ii) constrain its mass to $10.04^{+4.53}_{-3.10}\,M_{\rm Jup}$. If one adds the astrometry of $\beta$ Pictoris c, the mass is narrowed down to $9.15^{+1.08}_{-1.06}\,M_{\rm Jup}$. The inclusion of radial velocity measurements does not affect the orbital parameters significantly, but it does slightly decrease the mass estimate to $8.89^{+0.75}_{-0.75}\,M_{\rm Jup}$. With a semimajor axis of $2.68\pm0.02$\,au, a period of $1221\pm15$ days, and an eccentricity of $0.32\pm0.02$, the orbital parameters of $\beta$ Pictoris c are now constrained as precisely as those of $\beta$ Pictoris b. The orbital configuration is compatible with a high-order mean-motion resonance (7:1). The impact of the resonance on the planets' dynamics would then be negligible with respect to the secular perturbations, which might have played an important role in the eccentricity excitation of the outer planet., Comment: Accepted in A&A Letter

Context. Spiral arms in protoplanetary disks could be shown to be the manifestation of density waves launched by protoplanets and propagating in the gaseous component of the disk. At least two point sources have been identified in the L band in the MWC 758 system as planetary mass object candidates. Aims. We used VLT/SPHERE to search for counterparts of these candidates in the H and K bands, and to characterize the morphology of the spiral arms . Methods. The data were processed with now-standard techniques in high-contrast imaging to determine the limits of detection, and to compare them to the luminosity derived from L band observations. Results. In considering the evolutionary, atmospheric, and opacity models we were not able to confirm the two former detections of point sources performed in the L band. In addition, the analysis of the spiral arms from a dynamical point of view does not support the hypothesis that these candidates comprise the origin of the spirals. Conclusions. Deeper observations and longer timescales will be required to identify the actual source of the spiral arms in MWC 758., Comment: Accepted for publication in Astronomy and Astrophysics

We aim to reveal the nature of the reddest known substellar companion HD 206893 B by studying its near-infrared colors and spectral morphology and by investigating its orbital motion. We fit atmospheric models for giant planets and brown dwarfs and perform spectral retrievals with petitRADTRANS and ATMO on the observed GRAVITY, SPHERE, and GPI spectra of HD 206893 B. To recover its unusual spectral features, we include additional extinction by high-altitude dust clouds made of enstatite grains in the atmospheric model fits. We also infer the orbital parameters of HD 206893 B by combining the $\sim 100~\mu\text{as}$ precision astrometry from GRAVITY with data from the literature and constrain the mass and position of HD 206893 C based on the Gaia proper motion anomaly of the system. The extremely red color and the very shallow $1.4~\mu\text{m}$ water absorption feature of HD 206893 B can be fit well with the adapted atmospheric models and spectral retrievals. Altogether, our analysis suggests an age of $\sim 3$-$300~\text{Myr}$ and a mass of $\sim 5$-$30~\text{M}_\text{Jup}$ for HD 206893 B, which is consistent with previous estimates but extends the parameter space to younger and lower-mass objects. The GRAVITY astrometry points to an eccentric orbit ($e = 0.29^{+0.06}_{-0.11}$) with a mutual inclination of $< 34.4~\text{deg}$ with respect to the debris disk of the system. While HD 206893 B could in principle be a planetary-mass companion, this possibility hinges on the unknown influence of the inner companion on the mass estimate of $10^{+5}_{-4}~\text{M}_\text{Jup}$ from radial velocity and Gaia as well as a relatively small but significant Argus moving group membership probability of $\sim 61\%$. However, we find that if the mass of HD 206893 B is $< 30~\text{M}_\text{Jup}$, then the inner companion HD 206893 C should have a mass between $\sim 8$-$15~\text{M}_\text{Jup}$., Comment: 28 pages, 19 figures, accepted for publication in A&A

We present the multiple stellar systems observed within the SpHere INfrared survey for Exoplanet (SHINE). SHINE searched for substellar companions to young stars using high contrast imaging. Although stars with known stellar companions within SPHERE field of view (<5.5 arcsec) were removed from the original target list, we detected additional stellar companions to 78 of the 463 SHINE targets observed so far. 27% of the systems have three or more components. Given the heterogeneity of the sample in terms of observing conditions and strategy, tailored routines were used for data reduction and analysis, some of which were specifically designed for these data sets. We then combined SPHERE data with literature and archival ones, TESS light curves and Gaia parallaxes and proper motions, to characterise these systems as completely as possible. Combining all data, we were able to constrain the orbits of 25 systems. We carefully assessed the completeness of our sample for the separation range 50-500 mas (period range a few years - a few tens of years), taking into account the initial selection biases and recovering part of the systems excluded from the original list due to their multiplicity. This allowed us to compare the binary frequency for our sample with previous studies and highlight some interesting trends in the mass ratio and period distribution. We also found that, for the few objects for which such estimate was possible, the values of the masses derived from dynamical arguments were in good agreement with the model predictions. Stellar and orbital spins appear fairly well aligned for the 12 stars having enough data, which favour a disk fragmentation origin. Our results highlight the importance of combining different techniques when tackling complex problems such as the formation of binaries and show how large samples can be useful for more than one purpose., Comment: 59 pages, 11 tables, 18 figures

Large surveys with new-generation high-contrast imaging instruments are needed to derive the frequency and properties of exoplanet populations with separations from $\sim$5 to 300 AU. A careful assessment of the stellar properties is crucial for a proper understanding of when, where, and how frequently planets form, and how they evolve. The sensitivity of detection limits to stellar age makes this a key parameter for direct imaging surveys. We describe the SpHere INfrared survey for Exoplanets (SHINE), the largest direct imaging planet-search campaign initiated at the VLT in 2015 in the context of the SPHERE Guaranteed Time Observations of the SPHERE consortium. In this first paper we present the selection and the properties of the complete sample of stars surveyed with SHINE, focusing on the targets observed during the first phase of the survey (from February 2015 to February 2017). This early sample composed of 150 stars is used to perform a preliminary statistical analysis of the SHINE data, deferred to two companion papers presenting the survey performance, main discoveries, and the preliminary statistical constraints set by SHINE. Based on a large database collecting the stellar properties of all young nearby stars in the solar vicinity (including kinematics, membership to moving groups, isochrones, lithium abundance, rotation, and activity), we selected the original sample of 800 stars that were ranked in order of priority according to their sensitivity for planet detection in direct imaging with SPHERE. The properties of the stars that are part of the early statistical sample were revisited, including for instance measurements from the GAIA Data Release 2.

Over the past decades, direct imaging has confirmed the existence of substellar companions (exoplanets or brown dwarfs) on wide orbits (>10 au) from their host stars. To understand their formation and evolution mechanisms, we have initiated in 2015 the SPHERE infrared survey for exoplanets (SHINE), a systematic direct imaging survey of young, nearby stars to explore their demographics.} {We aim to detect and characterize the population of giant planets and brown dwarfs beyond the snow line around young, nearby stars. Combined with the survey completeness, our observations offer the opportunity to constrain the statistical properties (occurrence, mass and orbital distributions, dependency on the stellar mass) of these young giant planets.} {In this study, we present the observing and data analysis strategy, the ranking process of the detected candidates, and the survey performances for a subsample of 150 stars, which are representative of the full SHINE sample. The observations were conducted in an homogeneous way from February 2015 to February 2017 with the dedicated ground-based VLT/SPHERE instrument equipped with the IFS integral field spectrograph and the IRDIS dual-band imager covering a spectral range between 0.9 and 2.3 $\mu$m. We used coronographic, angular and spectral differential imaging techniques to reach the best detection performances for this study down to the planetary mass regime.}

Recent observations of resolved cold debris disks at tens of au have revealed that gaps could be a common feature in these Kuiper belt analogues. Such gaps could be evidence for the presence of planets within the gaps or closer-in near the edges of the disk. We present SPHERE observations of HD 92945 and HD 107146, two systems with detected gaps. We constrained the mass of possible companions responsible for the gap to 1-2 M Jup for planets located inside the gap and to less than 5 M Jup for separations down to 20 au from the host star. These limits allow us to exclude some of the possible configurations of the planetary systems proposed to explain the shape of the disks around these two stars. In order to put tighter limits on the mass at very short separations from the star, where direct imaging data are less effective, we also combined our data with astrometric measurements from Hipparcos and Gaia and radial velocity measurements. We were able to limit the separation and the mass of the companion potentially responsible for the proper motion anomaly of HD 107146 to values of 2-7 au and 2-5 M Jup , respectively., Comment: 17 pages, 14 Figures, accepted for publication on MNRAS

The planets HR8799bc display nearly identical colours and spectra as variable young exoplanet analogues such as VHS 1256-1257ABb and PSO J318.5-22, and are likely to be similarly variable. Here we present results from a 5-epoch SPHERE IRDIS broadband-$H$ search for variability in these two planets. HR 8799b aperture photometry and HR 8799bc negative simulated planet photometry share similar trends within uncertainties. Satellite spot lightcurves share the same trends as the planet lightcurves in the August 2018 epochs, but diverge in the October 2017 epochs. We consider $\Delta(mag)_{b} - \Delta(mag)_{c}$ to trace non-shared variations between the two planets, and rule out non-shared variability in $\Delta(mag)_{b} - \Delta(mag)_{c}$ to the 10-20$\%$ level over 4-5 hours. To quantify our sensitivity to variability, we simulate variable lightcurves by inserting and retrieving a suite of simulated planets at similar radii from the star as HR 8799bc, but offset in position angle. For HR 8799b, for periods $<$10 hours, we are sensitive to variability with amplitude $>5\%$. For HR 8799c, our sensitivity is limited to variability $>25\%$ for similar periods., Comment: 41 pages, 24 figures, accepted to MNRAS

Combining adaptive optics and interferometric observations results in a considerable contrast gain compared to single-telescope, extreme AO systems. Taking advantage of this, the ExoGRAVITY project is a survey of known young giant exoplanets located in the range of 0.1'' to 2'' from their stars. The observations provide astrometric data of unprecedented accuracy, being crucial for refining the orbital parameters of planets and illuminating their dynamical histories. Furthermore, GRAVITY will measure non-Keplerian perturbations due to planet-planet interactions in multi-planet systems and measure dynamical masses. Over time, repetitive observations of the exoplanets at medium resolution ($R=500$) will provide a catalogue of K-band spectra of unprecedented quality, for a number of exoplanets. The K-band has the unique properties that it contains many molecular signatures (CO, H$_2$O, CH$_4$, CO$_2$). This allows constraining precisely surface gravity, metallicity, and temperature, if used in conjunction with self-consistent models like Exo-REM. Further, we will use the parameter-retrieval algorithm petitRADTRANS to constrain the C/O ratio of the planets. Ultimately, we plan to produce the first C/O survey of exoplanets, kick-starting the difficult process of linking planetary formation with measured atomic abundances., Comment: SPIE 2020, invited talk

We present K-band interferometric observations of the PDS 70 protoplanets along with their host star using VLTI/GRAVITY. We obtained K-band spectra and 100 $\mu$as precision astrometry of both PDS 70 b and c in two epochs, as well as spatially resolving the hot inner disk around the star. Rejecting unstable orbits, we found a nonzero eccentricity for PDS 70 b of $0.17 \pm 0.06$, a near-circular orbit for PDS 70 c, and an orbital configuration that is consistent with the planets migrating into a 2:1 mean motion resonance. Enforcing dynamical stability, we obtained a 95% upper limit on the mass of PDS 70 b of 10 $M_\textrm{Jup}$, while the mass of PDS 70 c was unconstrained. The GRAVITY K-band spectra rules out pure blackbody models for the photospheres of both planets. Instead, the models with the most support from the data are planetary atmospheres that are dusty, but the nature of the dust is unclear. Any circumplanetary dust around these planets is not well constrained by the planets' 1-5 $\mu$m spectral energy distributions (SEDs) and requires longer wavelength data to probe with SED analysis. However with VLTI/GRAVITY, we made the first observations of a circumplanetary environment with sub-au spatial resolution, placing an upper limit of 0.3~au on the size of a bright disk around PDS 70 b., Comment: 30 pages, 9 figures, Accepted for publication in AJ

Young giant planets and brown dwarf companions emit near-infrared radiation that can be linearly polarized up to several percent. This polarization can reveal the presence of a circumsubstellar accretion disk, rotation-induced oblateness of the atmosphere, or an inhomogeneous distribution of atmospheric dust clouds. We measured the near-infrared linear polarization of 20 known directly imaged exoplanets and brown dwarf companions with the high-contrast imager SPHERE-IRDIS at the VLT. We reduced the data using the IRDAP pipeline to correct for the instrumental polarization and crosstalk with an absolute polarimetric accuracy <0.1% in the degree of polarization. We report the first detection of polarization originating from substellar companions, with a polarization of several tenths of a percent for DH Tau B and GSC 6214-210 B in H-band. By comparing the measured polarization with that of nearby stars, we find that the polarization is unlikely to be caused by interstellar dust. Because the companions have previously measured hydrogen emission lines and red colors, the polarization most likely originates from circumsubstellar disks. Through radiative transfer modeling, we constrain the position angles of the disks and find that the disks must have high inclinations. The presence of these disks as well as the misalignment of the disk of DH Tau B with the disk around its primary star suggest in situ formation of the companions. For the 18 other companions, we do not detect significant polarization and place subpercent upper limits on their degree of polarization. These non-detections may indicate the absence of circumsubstellar disks, a slow rotation rate of young companions, the upper atmospheres containing primarily submicron-sized dust grains, and/or limited cloud inhomogeneity. Finally, we present images of the circumstellar disks of DH Tau, GQ Lup, PDS 70, Beta Pic, and HD 106906., Comment: Accepted for publication in A&A. Shortened abstract. 29 pages, 22 figures

Sirius-like systems are wide binaries composed of a white dwarf (WD) and a companion of a spectral type earlier than M0. The WD progenitor evolves in isolation, but its wind during the AGB phase pollutes the companion surface and transfers some angular momentum. Within SHINE survey that uses SPHERE at the VLT, we acquired images of HD2133, HD114174, and CD-567708 and combined this data with high resolution spectra of the primaries, TESS, and literature data. We performed accurate abundance analyses for the MS. We found brighter J and K magnitudes for HD114174B than obtained previously and extended the photometry down to 0.95 micron. Our new data indicate a higher temperature and then shorter cooling age (5.57+/-0.02 Gyr) and larger mass (0.75+/-0.03 Mo) for this WD than previously assumed. This solved the discrepancy previously found with the age of the MS star. The two other WDs are less massive, indicating progenitors of ~1.3 Mo and 1.5-1.8 Mo for HD2133B and CD-56 7708B, respectively. We were able to derive constraints on the orbit for HD114174 and CD-56 7708. The composition of the MS stars agrees fairly well with expectations from pollution by the AGB progenitors of the WDs: HD2133A has a small enrichment of n-capture elements, which is as expected for pollution by an AGB star with a mass <1.5 Mo; CD-56 7708A is a previously unrecognized mild Ba-star, which is expected due to pollution by an AGB star with a mass in the range of 1.5-3.0 Mo; and HD114174 has a very moderate excess of n-capture elements, which is in agreement with the expectation for a massive AGB star to have a mass >3.0 Mo. On the other hand, none of these stars show the excesses of C that are expected to go along with those of n-capture elements. This might be related to the fact that these stars are at the edges of the mass range where we expect nucleosynthesis related to thermal pulses., Comment: Accepted on Astronomy and Astrophysics. 19 pages, 15 figures

Discovering wide companions of stellar systems allows us to constrain the dynamical environment and age of the latter. We studied four probable wide companions of four different stellar systems. The central stars are V4046 Sgr, HIP 74865, HIP 65426, and HIP 73145, and their probable wide companions are 2MASS J18152222-3249329, 2MASS J15174874-3028484, 2MASS J13242119-5129503, and 2MASS J14571503-3543505 respectively. V4046 Sgr is a member of $\beta$-Pictoris Moving Group while the rest of the stellar systems are acknowledged as members of the Scorpius-Centaurus association. The selected stellar systems are particularly interesting because all of them are already known to possess a low-mass companion and/or a spatially resolved disk. Identifying wider companions of these stars can improve their eligibility as benchmarks for understanding the formation channels of various triple systems, and can help us to determine the orbits of their possibly undiscovered inner, wider companions in case of higher multiplicity. We find that 2MASS J15174874-3028484 (0.11 $M_{\odot}$, 7.4$\pm$0.5 Myr), an already recognized pre-main sequence (PMS) member of Scorpius-Centaurus association, is a highly probable wide companion of HIP 74865. 2MASS J13242119-5129503 (0.04 $M_{\odot}$, 16$\pm$2.2 Myr) is ruled out as a plausible wide companion of HIP 65426, but confirmed to be a new sub-stellar member of the Scorpius-Centaurus association. 2MASS J14571503-3543505 (0.02 $M_{\odot}$, 17.75$\pm$4.15 Myr) is a probable sub-stellar member of the same association, but we cannot confirm whether or not it is gravitationally bound to HIP 73145. 2MASS J18152222-3249329 (0.3 $M_{\odot}$, older than 150 Myr) is determined to be a mildly active main sequence (MS) star, much older than members of $\beta$-Pictoris Moving Group, and unbound to V4046 Sgr despite their similar kinematic features., Comment: 13 pages, 9 figures, 8 tables

Methods used to detect giant exoplanets can be broadly divided into two categories: indirect and direct. Indirect methods are more sensitive to planets with a small orbital period, whereas direct detection is more sensitive to planets orbiting at a large distance from their host star. %, and thus on long orbital period. This dichotomy makes it difficult to combine the two techniques on a single target at once. Simultaneous measurements made by direct and indirect techniques offer the possibility of determining the mass and luminosity of planets and a method of testing formation models. Here, we aim to show how long-baseline interferometric observations guided by radial-velocity can be used in such a way. We observed the recently-discovered giant planet $\beta$ Pictoris c with GRAVITY, mounted on the Very Large Telescope Interferometer (VLTI). This study constitutes the first direct confirmation of a planet discovered through radial velocity. We find that the planet has a temperature of $T = 1250\pm50$\,K and a dynamical mass of $M = 8.2\pm0.8\,M_{\rm Jup}$. At $18.5\pm2.5$\,Myr, this puts $\beta$ Pic c close to a 'hot start' track, which is usually associated with formation via disk instability. Conversely, the planet orbits at a distance of 2.7\,au, which is too close for disk instability to occur. The low apparent magnitude ($M_{\rm K} = 14.3 \pm 0.1$) favours a core accretion scenario. We suggest that this apparent contradiction is a sign of hot core accretion, for example, due to the mass of the planetary core or the existence of a high-temperature accretion shock during formation., Comment: 10 pages

The rotational spectral modulation (spectro-photometric variability) of brown dwarfs is usually interpreted as a sign of the presence of inhomogeneous cloud covers in the atmosphere. This paper aims at exploring the role of temperature fluctuations in these spectral modulations. These fluctuations could naturally arise in a convective atmosphere impacted by diabatic processes such as complex chemistry, i.e. the recently proposed mechanism to explain the L/T transition: CO/CH4 radiative convection. We use the 1D radiative/convective code ATMO with ad-hoc modifications of the temperature gradient to model the rotational spectral modulation of 2MASS 1821, 2MASS 0136, and PSO 318.5-22. Modeling the spectral bright-to-faint ratio of the modulation of 2MASS 1821, 2MASS 0136, and PSO 318.5-22 shows that most spectral characteristics can be reproduced by temperature variations alone. Furthermore, the approximately anti-correlated variability between different wavelengths can be easily interpreted as a change in the temperature gradient in the atmosphere which is the consequence we expect from CO/CH4 radiative convection to explain the L/T transition. The deviation from an exact anti-correlation could then be interpreted as a phase shift similar to the hot-spot shift a different bandpasses in the atmosphere of hot Jupiters. Our results suggest that the rotational spectral modulation from cloud-opacity and temperature variations are degenerate. The detection of direct cloud spectral signatures, e.g. the silicate absorption feature at 10 um, would help to confirm the presence of clouds and their contribution to spectral modulations. Future studies looking at the differences in the spectral modulation of objects with and without the silicate absorption feature may give us some insight on how to distinguish cloud-opacity fluctuations from temperature fluctuations., Comment: accepted in A&A

In recent decades, thousands of substellar companions have been discovered with both indirect and direct methods of detection. In this paper, we focus our attention on substellar companions detected with the direct imaging technique, with the primary goal of investigating their close surroundings and looking for additional companions and satellites, as well as disks and rings. Any such discovery would shed light on many unresolved questions, particularly with regard to their possible formation mechanisms. To reveal bound features of directly imaged companions we need to suppress the contribution from the source itself. Therefore, we developed a method based on the negative fake companion (NEGFC) technique that first estimates the position in the field of view (FoV) and the flux of the imaged companion, then subtracts a rescaled model point spread function (PSF) from the imaged companion. Next it performs techniques, such as angular differential imaging (ADI), to further remove quasi-static patterns of the star. We applied the method to the sample of substellar objects observed with SPHERE during the SHINE GTO survey. Among the 27 planets and brown dwarfs we analyzed, we detected a possible point source close to DH Tau B. This candidate companion was detected in four different SPHERE observations, with an estimated mass of $\sim 1$ M\textsubscript{Jup}, and a mass ratio with respect to the brown dwarf of $1/10$. This binary system, if confirmed, would be the first of its kind, opening up interesting questions for the formation mechanism, evolution, and frequency of such pairs. In order to address the latter, the residuals and contrasts reached for 25 companions in the sample of substellar objects observed with SPHERE were derived. If the DH Tau Bb companion is real, the binary fraction obtained is $\sim 7\%$, which is in good agreement with the results obtained for field brown dwarfs., Comment: 19 pages, 17 figures

The SHINE project is a 500-star survey performed with SPHERE on the VLT for the purpose of directly detecting new substellar companions and understanding their formation and early evolution. Here we present an initial statistical analysis for a subsample of 150 stars that are representative of the full SHINE sample. Our goal is to constrain the frequency of substellar companions with masses between 1 and 75 MJup and semimajor axes between 5 and 300 au. We adopt detection limits as a function of angular separation from the survey data for all stars converted into mass and projected orbital separation using the BEX-COND-hot evolutionary tracks and known distance to each system. Based on the results obtained for each star and on the 13 detections in the sample, we use a MCMC tool to compare our observations to two different types of models. The first is a parametric model based on observational constraints, and the second type are numerical models that combine advanced core accretion and gravitational instability planet population synthesis. Using the parametric model, we show that the frequencies of systems with at least one substellar companion are $23.0_{-9.7}^{+13.5}\%$, $5.8_{-2.8}^{+4.7}\%$, and $12.6_{-7.1}^{+12.9}\%$ for BA, FGK, and M stars, respectively. We also demonstrate that a planet-like formation pathway probably dominates the mass range from 1-75 MJup for companions around BA stars, while for M dwarfs, brown dwarf binaries dominate detections. In contrast, a combination of binary star-like and planet-like formation is required to best fit the observations for FGK stars. Using our population model and restricting our sample to FGK stars, we derive a frequency of $5.7_{-2.8}^{+3.8}\%$, consistent with predictions from the parametric model. More generally, the frequency values that we derive are in excellent agreement with values obtained in previous studies., Comment: 24 pages, 14 figures, 3 tables. Accepted for publication in A&A

The determination of the fundamental properties (mass, separation, age, gravity and atmospheric properties) of brown dwarf companions allows us to infer crucial informations on their formation and evolution mechanisms. Spectroscopy of substellar companions is available to date only for a limited number of objects (and mostly at very low resolution, R<50) because of technical limitations, i.e., contrast and angular resolution. We present medium resolution (R=350), coronagraphic long-slit spectroscopic observations with SPHERE of two substellar companions, HD 1160 B and HD 19467 B. We found that HD 1160 B has a peculiar spectrum that cannot be fitted by spectra in current spectral libraries. A good fit is possible only considering separately the Y+J and the H spectral band. The spectral type is between M5 and M7. We also estimated a T_eff of 2800-2900 K and a log(g) of 3.5-4.0 dex. The low surface gravity seems to favour young age (10-20 Myr) and low mass (~20 M Jup ) for this object. HD 19467 B is instead a fully evolved object with a T_eff of ~1000 K and log g of ~5.0 dex. Its spectral type is T6+/-1., Comment: 14 pages, 14 Figures. Accepted for publication on MNRAS

Context. Detecting and characterizing substellar companions for which the luminosity, mass, and age can be determined independently is of utter importance to test and calibrate the evolutionary models due to uncertainties in their formation mechanisms. HD 19467 is a bright and nearby star hosting a cool brown dwarf companion detected with RV and imaging, making it a valuable object for such studies. Aims. We aim to further characterize the orbital, spectral, and physical properties of the HD 19467 system. Methods. We present new high-contrast imaging data with the SPHERE and NaCo instruments. We also analyze archival data from HARPS, NaCo, HIRES, UVES, and ASAS. We also use proper motion data of the star from Hipparcos and Gaia. Results. We refine the properties of the host star and derive an age of 8.0$^{+2.0}_{-1.0}$ Gyr based on isochrones, gyrochronology, and chemical and kinematic arguments. This estimate is slightly younger than previous estimates of ~9-11 Gyr. No orbital curvature is seen in the current imaging, RV, and astrometric data. From a joint fit of the data, we refine the orbital parameters for HD 19467B: period 398$^{+95}_{-93}$ yr, inclination 129.8$^{+8.1}_{-5.1}$ deg, eccentricity 0.56$\pm$0.09, longitude of the ascending node 134.8$\pm$4.5 deg, and argument of the periastron 64.2$^{+5.5}_{-6.3}$ deg. We assess a dynamical mass of 74$^{+12}_{-9}$ MJ. The fit with atmospheric models of the spectrophotometric data of HD 19467B indicates an atmosphere without clouds or with very thin clouds, an effective temperature of 1042$^{+77}_{-71}$ K, and a large surface gravity of 5.34$^{+0.08}_{-0.09}$ dex. The comparison to model predictions of the bolometric luminosity and dynamical mass of HD 19467B, assuming our system age estimate, indicates a better agreement with the Burrows et al. models; whereas the other evolutionary models used tend to underestimate its cooling rate., Comment: Accepted for publication in A&A, 25 pages, 23 figures, 9 tables. Abstract slightly abridged to match arXiv requirements. Updated after language editing

Dusty debris disks around pre- and main-sequence stars are potential signposts for the existence of planetesimals and exoplanets. Giant planet formation is therefore expected to play a key role in the evolution of the disk. This is indirectly confirmed by extant sub-millimeter near-infrared images of young protoplanetary and cool dusty debris disks around main sequence stars usually showing substantial spatial structures. A majority of recent discoveries of imaged giant planets have been obtained around young, early-type stars hosting a circumstellar disk. In this context, we have carried out a direct imaging program designed to maximize our chances of giant planet discovery and targeting twenty-two young, early-type stars. About half of them show indication of multi-belt architectures. Using the IRDIS dual-band imager and the IFS integral field spectrograph of SPHERE to acquire high-constrast coronagraphic differential near-infrared images, we have conducted a systematic search in the close environment of these young, dusty and early-type stars. We confirmed that companions detected around HIP 34276, HIP 101800 and HIP 117452 are stationary background sources and binary companions. The companion candidates around HIP 8832, HIP 16095 and HIP 95619 are determined as background contamination. For stars for which we infer the presence of debris belts, a theoretical minimum mass for planets required to clear the debris gaps can be calculated . The dynamical mass limit is at least $0.1 M_J$ and can exceed $1 M_J$. Direct imaging data is typically sensitive to planets down to $\sim 3.6 M_J$ at 1 $''$, and $1.7 M_J$ in the best case. These two limits tightly constrain the possible planetary systems present around each target. These systems will be probably detectable with the next generation of planet imagers., Comment: Accepted for publication in Astronomy & Astrophysics. 12 pages, 9 figures