Search

Owing to their similarities with giant exoplanets, brown dwarf companions of stars provide insights into the fundamental processes of planet formation and evolution. From their orbits, several brown dwarf companions are found to be more massive than theoretical predictions given their luminosities and the ages of their host stars (e.g. Brandt et al. 2021, Cheetham et al. 2018, Li et al. 2023). Either the theory is incomplete or these objects are not single entities. For example, they could be two brown dwarfs each with a lower mass and intrinsic luminosity (Brandt et al. 2021, Howe et al. 2024). The most problematic example is Gliese 229 B (Nakajima et al. 1995, Oppenheimer et al. 1995), which is at least 2-6 times less luminous than model predictions given its dynamical mass of $71.4\pm0.6$ Jupiter masses ($M_{\rm Jup}$) (Brandt et al. 2021). We observed Gliese 229 B with the GRAVITY interferometer and, separately, the CRIRES+ spectrograph at the Very Large Telescope. Both sets of observations independently resolve Gliese 229 B into two components, Gliese 229 Ba and Bb, settling the conflict between theory and observations. The two objects have a flux ratio of $0.47\pm0.03$ at a wavelength of 2 $\mu$m and masses of $38.1\pm1.0$ and $34.4\pm1.5$ $M_{\rm Jup}$, respectively. They orbit each other every 12.1 days with a semimajor axis of 0.042 astronomical units (AU). The discovery of Gliese 229 BaBb, each only a few times more massive than the most massive planets, and separated by 16 times the Earth-moon distance, raises new questions about the formation and prevalence of tight binary brown dwarfs around stars., Comment: Published in Nature. The Version of Record of this article is located at https://www.nature.com/articles/s41586-024-08064-x

Ground-based thermal infrared observations face substantial challenges in correcting the predominant background emitted as thermal radiation from the atmosphere and the telescope itself. With the upcoming 40\,m class ELTs, unprecedented sensitivities from ground will be reached, underlining the need of even more sophisticated background correction strategies. This study aims to investigate the impact of thermal backgrounds on ground-based observations and identify possible limiting factors in dedicated correction strategies. We evaluate temporal and spatial characteristics of the thermal background in direct imaging data obtained with different telescopes and observation modes. In particular, three distinct datasets, acquired using VLT/NACO and KECK/NIRC2, are analyzed. Our analysis reveals that the observations are not fully photon shot noise limited, but exhibit additional sensitivity losses caused by imperfect background compensation in the different datasets. We identify correlations between background fluctuations and the activity of the adaptive optics system. We hypothesize that the pupil modulation of the adaptive optics mirrors introduces high frequency spatial and temporal fluctuations to the background, which could ultimately constrain the detection limit if they are not compensated adequately., Comment: 11 pages, 9 figures. Accepted by Publications of the Astronomical Society of the Pacific

The dust- and gas-rich protoplanetary disks around young stellar systems play a key role in star and planet formation. While considerable progress has recently been made in probing these disks on large scales of a few tens of astronomical units (au), the central au needs to be more investigated. We aim at unveiling the physical processes at play in the innermost regions of the strongly accreting T Tauri Star S CrA N by means of near-infrared interferometric observations. The K-band continuum emission is well reproduced with an azimuthally-modulated dusty ring. As the star alone cannot explain the size of this sublimation front, we propose that magnetospheric accretion is an important dust-heating mechanism leading to this continuum emission. The differential analysis of the Hydrogen Br$\gamma$ line is in agreement with radiative transfer models combining magnetospheric accretion and disk winds. Our observations support an origin of the Br$\gamma$ line from a combination of (variable) accretion-ejection processes in the inner disk region.

The first generation of ELT instruments includes an optical-infrared high-resolution spectrograph, indicated as ELT-HIRES and recently christened ANDES (ArmazoNes high Dispersion Echelle Spectrograph). ANDES consists of three fibre-fed spectrographs ([U]BV, RIZ, YJH) providing a spectral resolution of $\sim$100,000 with a minimum simultaneous wavelength coverage of 0.4-1.8 $\mu$m with the goal of extending it to 0.35-2.4 $\mu$m with the addition of a U arm to the BV spectrograph and a separate K band spectrograph. It operates both in seeing- and diffraction-limited conditions and the fibre feeding allows several, interchangeable observing modes including a single conjugated adaptive optics module and a small diffraction-limited integral field unit in the NIR. Modularity and fibre-feeding allow ANDES to be placed partly on the ELT Nasmyth platform and partly in the Coud\'e room. ANDES has a wide range of groundbreaking science cases spanning nearly all areas of research in astrophysics and even fundamental physics. Among the top science cases, there are the detection of biosignatures from exoplanet atmospheres, finding the fingerprints of the first generation of stars, tests on the stability of Nature's fundamental couplings, and the direct detection of the cosmic acceleration. The ANDES project is carried forward by a large international consortium, composed of 35 Institutes from 13 countries, forming a team of almost 300 scientists and engineers which include the majority of the scientific and technical expertise in the field that can be found in ESO member states., Comment: SPIE astronomical telescope and instrumentation 2024, in press

We present new VLTI/GRAVITY near-infrared interferometric measurements of the angular size of the innermost hot dust continuum for 14 type 1 AGNs. The angular sizes are resolved on scales of ~0.7 mas and the inferred ring radii range from 0.028 to 1.33 pc, comparable to those reported previously and a factor 10-20 smaller than the mid-infrared sizes in the literature. Combining our new data with previously published values, we compile a sample of 25 AGN with bolometric luminosity ranging from $10^{42}$ to $10^{47} \rm erg~s^{-1}$, with which we study the radius-luminosity (R-L) relation for the hot dust structure. Our interferometric measurements of radius are offset by a factor 2 from the equivalent relation derived through reverberation mapping. Using a simple model to explore the dust structure's geometry, we conclude that this offset can be explained if the 2 um emitting surface has a concave shape. Our data show that the slope of the relation is in line with the canonical $R \propto L^{0.5}$ when using an appropriately non-linear correction for bolometric luminosity. In contrast, using optical luminosity or applying a constant bolometric correction to it results in a significant deviation in the slope, suggesting a potential luminosity dependence on the spectral energy distribution. Over four orders of magnitude in luminosity, the intrinsic scatter around the R-L relation is 0.2 dex, suggesting a tight correlation between innermost hot dust structure size and the AGN luminosity., Comment: 14 pages, 5 figures, 4 tables, accepted for publication in A&A

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

The detection of low-mass planets orbiting the nearest stars is a central stake of exoplanetary science, as they can be directly characterized much more easily than their distant counterparts. Here, we present the results of our long-term astrometric observations of the nearest binary M-dwarf Gliese 65 AB (GJ65), located at a distance of only 2.67 pc. We monitored the relative astrometry of the two components from 2016 to 2023 with the VLTI/GRAVITY interferometric instrument. We derived highly accurate orbital parameters for the stellar system, along with the dynamical masses of the two red dwarfs. The GRAVITY measurements exhibit a mean accuracy per epoch of 50-60 microarcseconds in 1.5h of observing time using the 1.8m Auxiliary Telescopes. The residuals of the two-body orbital fit enable us to search for the presence of companions orbiting one of the two stars (S-type orbit) through the reflex motion they imprint on the differential A-B astrometry. We detected a Neptune-mass candidate companion with an orbital period of p = 156 +/- 1 d and a mass of m = 36 +/- 7 Mearth. The best-fit orbit is within the dynamical stability region of the stellar pair. It has a low eccentricity, e = 0.1 - 0.3, and the planetary orbit plane has a moderate-to-high inclination of i > 30{\deg} with respect to the stellar pair, with further observations required to confirm these values. These observations demonstrate the capability of interferometric astrometry to reach microarcsecond accuracy in the narrow-angle regime for planet detection by reflex motion from the ground. This capability offers new perspectives and potential synergies with Gaia in the pursuit of low-mass exoplanets in the solar neighborhood., Comment: Corresponding authors: G.Bourdarot, P.Kervella, O.Pfuhl. Accepted in A&A Letters

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

Temporal variability in the photometric and spectroscopic properties of protoplanetary disks is common in YSO. However, evidence pointing toward changes in their morphology over short timescales has only been found for a few sources, mainly due to a lack of high cadence observations at mas resolution. We combine GRAVITY multi-epoch observations of HD98922 at mas resolution with PIONIER archival data covering a total time span of 11 years. We interpret the interferometric visibilities and spectral energy distribution with geometrical models and through radiative transfer techniques. We investigated high-spectral-resolution quantities to obtain information on the properties of the HI BrG-line-emitting region. The observations are best fitted by a model of a crescent-like asymmetric dust feature located at 1 au and accounting for 70% of the NIR emission. The feature has an almost constant magnitude and orbits the central star with a possible sub-Keplerian period of 12 months, although a 9 month period is another, albeit less probable, solution. The radiative transfer models show that the emission originates from a small amount of carbon-rich (25%) silicates, or quantum-heated particles located in a low-density region. Among different possible scenarios, we favor hydrodynamical instabilities in the inner disk that can create a large vortex. The high spectral resolution differential phases in the BrG-line show that the hot-gas component is offset from the star and in some cases is located between the star and the crescent feature. The scale of the emission does not favor magnetospheric accretion as a driving mechanism. The scenario of an asymmetric disk wind or a massive accreting substellar or planetary companion is discussed. With this unique observational data set for HD98922, we reveal morphological variability in the innermost 2 au of its disk region., Comment: 45 pages, 20 figures, accepted by and to be published in Astronomy & Astrophysics (A&A)

Tight relationships exist in the local universe between the central stellar properties of galaxies and the mass of their supermassive black hole. These suggest galaxies and black holes co-evolve, with the main regulation mechanism being energetic feedback from accretion onto the black hole during its quasar phase. A crucial question is how the relationship between black holes and galaxies evolves with time; a key epoch to probe this relationship is at the peaks of star formation and black hole growth 8-12 billion years ago (redshifts 1-3). Here we report a dynamical measurement of the mass of the black hole in a luminous quasar at a redshift of 2, with a look back time of 11 billion years, by spatially resolving the broad line region. We detect a 40 micro-arcsecond (0.31 pc) spatial offset between the red and blue photocenters of the H$\alpha$ line that traces the velocity gradient of a rotating broad line region. The flux and differential phase spectra are well reproduced by a thick, moderately inclined disk of gas clouds within the sphere of influence of a central black hole with a mass of 3.2x10$^{8}$ solar masses. Molecular gas data reveal a dynamical mass for the host galaxy of 6x10$^{11}$ solar masses, which indicates an under-massive black hole accreting at a super-Eddington rate. This suggests a host galaxy that grew faster than the supermassive black hole, indicating a delay between galaxy and black hole formation for some systems., Comment: 5 pages Main text, 8 figures, 2 tables, to be published in Nature, under embargo until 29 January 2024 16:00 (London)

The 2nd generation VLTI instrument GRAVITY is a two-field infrared interferometer operating in the K band between 1.97 and 2.43 $\mu$m with either the four 8 m or the four 1.8 m telescopes of the Very Large Telescope (VLT). Beams collected by the telescopes are corrected with adaptive optics systems and the fringes are stabilized with a fringe-tracking system. A metrology system allows the measurement of internal path lengths in order to achieve high-accuracy astrometry. High sensitivity and high interferometric accuracy are achieved thanks to (i) correction of the turbulent phase, (ii) the use of low-noise detectors, and (iii) the optimization of photometric and coherence throughput. Beam combination and most of the beam transport are performed with single-mode waveguides in vacuum and at low temperature. In this paper, we present the functions and performance achieved with weakly birefringent standard single-mode fiber systems in GRAVITY. Fibered differential delay lines (FDDLs) are used to dynamically compensate for up to 6 mm of delay between the science and reference targets. Fibered polarization rotators allow us to align polarizations in the instrument and make the single-mode beam combiner close to polarization neutral. The single-mode fiber system exhibits very low birefringence (less than 23{\deg}), very low attenuation (3.6-7 dB/km across the K band), and optimized differential dispersion (less than 2.04 $\mu$rad cm2 at zero extension of the FDDLs). As a consequence, the typical fringe contrast losses due to the single-mode fibers are 6% to 10% in the lowest-resolution mode and 5% in the medium- and high-resolution modes of the instrument for a photometric throughput of the fiber chain of the order of 90%. There is no equivalent of this fiber system to route and modally filter beams with delay and polarization control in any other K-band beamcombiner., Comment: 14 pages including 2 appendices, 8 figures in the main text plus 2 figures in the appendices, final version published in A&A

We aim to investigate the origin of the HI Br$\gamma$ emission in young stars by using GRAVITY to image the innermost region of circumstellar disks, where important physical processes such as accretion and winds occur. With high spectral and angular resolution, we focus on studying the continuum and the HI Br$\gamma$-emitting area of the Herbig star HD58647. Using VLTI-GRAVITY, we conducted observations of HD58647 with both high spectral and high angular resolution. Thanks to the extensive $uv$ coverage, we were able to obtain detailed images of the circumstellar environment at a sub-au scale, specifically capturing the continuum and the Br$\gamma$-emitting region. Through the analysis of velocity-dispersed images and photocentre shifts, we were able to investigate the kinematics of the HI Br$\gamma$-emitting region. The recovered continuum images show extended emission where the disk major axis is oriented along a position angle of 14\degr. The size of the continuum emission at 5-sigma levels is $\sim$ 1.5 times more extended than the sizes reported from geometrical fitting (3.69 mas $\pm$ 0.02 mas). This result supports the existence of dust particles close to the stellar surface, screened from the stellar radiation by an optically thick gaseous disk. Moreover, for the first time with GRAVITY, the hot gas component of HD58647 traced by the Br$\gamma$ ,has been imaged. This allowed us to constrain the size of the Br$\gamma$-emitting region and study the kinematics of the hot gas; we find its velocity field to be roughly consistent with gas that obeys Keplerian motion. The velocity-dispersed images show that the size of the hot gas emission is from a more compact region than the continuum (2.3 mas $\pm$ 0.2 mas). Finally, the line phases show that the emission is not entirely consistent with Keplerian rotation, hinting at a more complex structure in the hot gaseous disk., Comment: accepted in A&A on 20/11/2023

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 dark compact object at the centre of the Milky Way is well established to be a supermassive black hole with mass $M_{\bullet} \sim 4.3 \cdot 10^6 \, M_{\odot}$, but the nature of its environment is still under debate. In this work, we used astrometric and spectroscopic measurements of the motion of the star S2, one of the closest stars to the massive black hole, to determine an upper limit on an extended mass composed of a massive vector field around Sagittarius A*. For a vector with effective mass $10^{-19} \, \rm eV \lesssim m_s \lesssim 10^{-18} \, \rm eV$, our Markov Chain Monte Carlo analysis shows no evidence for such a cloud, placing an upper bound $M_{\rm cloud} \lesssim 0.1\% M_{\bullet}$ at $3\sigma$ confidence level. We show that dynamical friction exerted by the medium on S2 motion plays no role in the analysis performed in this and previous works, and can be neglected thus., Comment: 9 pages, 5 figures, accepted to MNRAS

The goal of this work is to characterize the polarization effects of the VLTI and GRAVITY. This is needed to calibrate polarimetric observations with GRAVITY for instrumental effects and to understand the systematic error introduced to the astrometry due to birefringence when observing targets with a significant intrinsic polarization. By combining a model of the VLTI light path and its mirrors and dedicated experimental data, we construct a full polarization model of the VLTI UTs and the GRAVITY instrument. We first characterize all telescopes together to construct a UT calibration model for polarized targets. We then expand the model to include the differential birefringence. With this, we can constrain the systematic errors for highly polarized targets. Together with this paper, we publish a standalone Python package to calibrate the instrumental effects on polarimetric observations. This enables the community to use GRAVITY to observe targets in a polarimetric observing mode. We demonstrate the calibration model with the galactic center star IRS 16C. For this source, we can constrain the polarization degree to within 0.4 % and the polarization angle within 5 deg while being consistent with the literature. Furthermore, we show that there is no significant contrast loss, even if the science and fringe-tracker targets have significantly different polarization, and we determine that the phase error in such an observation is smaller than 1 deg, corresponding to an astrometric error of 10 {\mu}as. With this work, we enable the use of the polarimetric mode with GRAVITY/UTs and outline the steps necessary to observe and calibrate polarized targets. We demonstrate that it is possible to measure the intrinsic polarization of astrophysical sources with high precision and that polarization effects do not limit astrometric observations of polarized targets., Comment: Accepted by A&A

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

We present new astrometric and polarimetric observations of flares from Sgr A* obtained with GRAVITY, the near-infrared interferometer at ESO's Very Large Telescope Interferometer (VLTI), bringing the total sample of well-covered astrometric flares to four and polarimetric ones to six, where we have for two flares good coverage in both domains. All astrometric flares show clockwise motion in the plane of the sky with a period of around an hour, and the polarization vector rotates by one full loop in the same time. Given the apparent similarities of the flares, we present a common fit, taking into account the absence of strong Doppler boosting peaks in the light curves and the EHT-measured geometry. Our results are consistent with and significantly strengthen our model from 2018: We find that a) the combination of polarization period and measured flare radius of around nine gravitational radii ($9 R_g \approx 1.5 R_{ISCO}$, innermost stable circular orbit) is consistent with Keplerian orbital motion of hot spots in the innermost accretion zone. The mass inside the flares' radius is consistent with the $4.297 \times 10^6 \; \text{M}_\odot$ measured from stellar orbits at several thousand $R_g$. This finding and the diameter of the millimeter shadow of Sgr A* thus support a single black hole model. Further, b) the magnetic field configuration is predominantly poloidal (vertical), and the flares' orbital plane has a moderate inclination with respect to the plane of the sky, as shown by the non-detection of Doppler-boosting and the fact that we observe one polarization loop per astrometric loop. Moreover, c) both the position angle on sky and the required magnetic field strength suggest that the accretion flow is fueled and controlled by the winds of the massive, young stars of the clockwise stellar disk 1-5 arcsec from Sgr A*, in agreement with recent simulations., Comment: 10 pages, 12 figures. Accepted for publication in A&A

The motion of S2, one of the stars closest to the Galactic Centre, has been measured accurately and used to study the compact object at the centre of the Milky Way. It is commonly accepted that this object is a supermassive black hole but the nature of its environment is open to discussion. Here, we investigate the possibility that dark matter in the form of an ultralight scalar field ``cloud'' clusters around Sgr~A*. We use the available data for S2 to perform a Markov Chain Monte Carlo analysis and find the best-fit estimates for a scalar cloud structure. Our results show no substantial evidence for such structures. When the cloud size is of the order of the size of the orbit of S2, we are able to constrain its mass to be smaller than $0.1\%$ of the central mass, setting a strong bound on the presence of new fields in the galactic centre., Comment: Published on MNRAS. References added

In the Milky Way the central massive black hole, SgrA*, coexists with a compact nuclear star cluster that contains a sub-parsec concentration of fast-moving young stars called S-stars. Their location and age are not easily explained by current star formation models, and in several scenarios the presence of an intermediate-mass black hole (IMBH) has been invoked. We use GRAVITY astrometric and SINFONI, KECK, and GNIRS spectroscopic data of S2 to investigate whether a second massive object could be present deep in the Galactic Centre (GC) in the form of an IMBH binary companion to SgrA*. To solve the three-body problem, we used a post-Newtonian framework and consider two types of settings: (i) a hierarchical set-up where the star S2 orbits the SgrA* - IMBH binary and (ii) a non-hierarchical set-up where the IMBH trajectory lies outside the S2 orbit. In both cases we explore the full 20-dimensional parameter space by employing a Bayesian dynamic nested sampling method. For the hierarchical case we find: IMBH masses > 2000 Msun on orbits with smaller semi-major axes than S2 are largely excluded. For the non-hierarchical case the parameter space contains several pockets of valid IMBH solutions. However, a closer analysis of their impact on the resident stars reveals that IMBHs on semi-major axes larger than S2 tend to disrupt the S-star cluster in less than a million years. This makes the existence of an IMBH among the S-stars highly unlikely. The current S2 data do not formally require the presence of an IMBH. If an IMBH hides in the GC, it has to be either a low-mass IMBH inside the S2 orbit that moves on a short and significantly inclined trajectory or an IMBH with a semi-major axis >1". We provide the parameter maps of valid IMBH solutions in the GC and discuss the general structure of our results. (abridged), Comment: 11 pages, 6 figures, published in A $\&$ A

The mass-loss mechanisms in M-type AGB stars are not well understood, in particular, the formation of dust-driven winds from the innermost gaseous layers around these stars. One way to understand the gas-dust interaction in these regions and its impact on the mass-loss mechanisms is through the analysis of high-resolution observations of the stellar surface and its closest environment. We aim at characterizing the inner circumstellar environment (~3 R*) of the M-type Mira star R Car in the near-infrared at different phases of a pulsation period. We used GRAVITY interferometric observations in the K-band obtained at two different epochs over 2018. Those data were analyzed using parametric models and image reconstruction of both the pseudo-continuum and the CO band-heads observed. The reported data are the highest angular resolution observations on the source in the K-band. We determine sizes of R Car's stellar disk of 16.67 +- 0.05 mas (3.03 au) in January 2018 and 14.84+-0.06 mas (2.70 au) in February 2018, respectively. From our physical model, we determined temperatures and size ranges for the innermost CO layer detected around R Car. We find that magnesium composites, Mg2SiO4 and MgSiO3, have temperatures and condensation distances consistent with the ones obtained for the CO layer model and pure-line reconstructed images, being them the most plausible dust types responsible of wind formation. Our reconstructed images show evidence of asymmetrical and inhomogeneous structures, which might trace a complex and perhaps clumpy structure of the CO molecule distribution. Our work demonstrates that the conditions for dust nucleation and thus for initialising dust-driven winds in M-type AGB stars are met in R Car and we identify Magnesium composites as the most probable candidates. This observational evidence is crucial to constrain the role of convection and pulsation in M-type stars., Comment: Accepted to be published by A&A, 20 pages, 21 figures

[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

Aims: We aim to spatially and spectrally resolve the Br-gamma hydrogen emission line with the methods of interferometry in order to examine the kinematics of the hydrogen gas emission region in the inner accretion disk of a sample of solar-like young stellar objects. The goal is to identify trends and categories among the sources of our sample and to discuss whether or not they can be tied to different origin mechanisms associated with Br-gamma emission in T Tauri stars, chiefly and most prominently magnetospheric accretion. Methods: We observed a sample of seven T Tauri stars for the first time with VLTI GRAVITY, recording spectra and spectrally dispersed interferometric quantities across the Br-gamma line in the NIR K-band. We use them to extract the size of the Br-gamma emission region and the photocenter shifts. To assist in the interpretation, we also make use of radiative transfer models of magnetospheric accretion to establish a baseline of expected interferometric signatures if accretion is the primary driver of Br-gamma emission. Results: From among our sample, we find that five of the seven T~Tauri stars show an emission region with a half-flux radius in the range broadly expected for magnetospheric truncation. Two of the five objects also show Br-gamma emission primarily originating from within the corotation radius, while two other objects exhibit extended emission on a scale beyond 10 R$_*$, one of them even beyond the K~band continuum half-flux radius of 11.3 R$_*$. Conclusions: We find strong evidence to suggest that for the two weakest accretors in the sample, magnetospheric accretion is the primary driver of Br-gamma radiation. The results for the remaining sources imply either partial or strong contributions coming from spatially extended emission components in the form of outflows, such as stellar or disk winds., Comment: Version corresponds to the one in A&A production. The author list has been amended with "GRAVITY Collaboration" as first entry in keeping with the convention established for papers published by the GRAVITY consortium in recent years

This work focuses on active galactic nuclei (AGNs), and the relation between the sizes of the hot dust continuum and the broad-line region (BLR). We find that the continuum size measured using optical/near-infrared interferometry (OI) is roughly twice that measured by reverberation mapping (RM). Both OI and RM continuum sizes show a tight relation with the H$\beta$ BLR size with only an intrinsic scatter of 0.25 dex. The masses of supermassive black holes (BHs) can hence be simply derived from a dust size in combination with a broad line width and virial factor. Since the primary uncertainty of these BH masses comes from the virial factor, the accuracy of the continuum-based BH masses is close to those based on the RM measurement of the broad emission line. Moreover, the necessary continuum measurements can be obtained on a much shorter timescale than those required monitoring for RM, and are also more time efficient than those needed to resolve the BLR with OI. The primary goal of this work is to demonstrate measuring the BH mass based on the dust continuum size with our first calibration of the $R_\mathrm{BLR}$-$R_\mathrm{d}$ relation. The current limitation and caveats are discussed in detail. Future GRAVITY observations are expected to improve the continuum-based method and have the potential to measure BH masses for a large sample of AGNs in the low-redshift Universe., Comment: Accepted for publication in A&A; 11 pages, 5 figures, 3 tables

As part of the GRAVITY$^{+}$ project, the near-infrared beam combiner GRAVITY and the VLTI are currently undergoing a series of significant upgrades to further improve the performance and sky coverage. The instrumental changes will be transformational, and for instance uniquely position GRAVITY to observe the broad line region of hundreds of Active Galactic Nuclei (AGN) at a redshift of two and higher. The increased sky coverage is achieved by enlarging the maximum angular separation between the celestial science object (SC) and the off-axis fringe tracking (FT) star from currently 2 arcseconds (arcsec) up to unprecedented 30 arcsec, limited by the atmospheric conditions. This was successfully demonstrated at the VLTI for the first time., Comment: Proceeding of SPIE Astronomical Telescopes + Instrumentation 2022

With the upgrade from GRAVITY to GRAVITY$^+$ the instrument will evolve into an all-sky interferometer that can observe faint targets, such as high redshift AGN. Observing the faintest targets requires reducing the noise sources in GRAVITY as much as possible. The dominant noise source, especially in the blue part of the spectrum, is the backscattering of the metrology laser light onto the detector. To reduce this noise we introduce two new metrology modes. With a combination of small hardware changes and software adaptations, we can dim the metrology laser during the observation without losing the phase referencing. For single beam targets, we can even turn off the metrology laser for the maximum SNR on the detector. These changes lead to an SNR improvement of over a factor of two averaged over the whole spectrum and up to a factor of eight in the part of the spectrum currently dominated by laser noise., Comment: Proceeding of SPIE Astronomical Telescopes + Instrumentation 2022

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

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

GRAVITY+ is the upgrade of GRAVITY and the Very Large Telescope Interferometer (VLTI) with wide-separation fringe tracking, new adaptive optics, and laser guide stars on all four 8~m Unit Telescopes (UTs), for ever fainter, all-sky, high contrast, milliarcsecond interferometry. Here we present the design and first results of the first phase of GRAVITY+, called GRAVITY Wide. GRAVITY Wide combines the dual-beam capabilities of the VLTI and the GRAVITY instrument to increase the maximum separation between the science target and the reference star from 2 arcseconds with the 8 m UTs up to several 10 arcseconds, limited only by the Earth's turbulent atmosphere. This increases the sky-coverage of GRAVITY by two orders of magnitude, opening up milliarcsecond resolution observations of faint objects, and in particular the extragalactic sky. The first observations in 2019 - 2022 include first infrared interferometry of two redshift $z\sim2$ quasars, interferometric imaging on the binary system HD 105913A, and repeated observations of multiple star systems in the Orion Trapezium Cluster. We find the coherence loss between the science object and fringe-tracking reference star well described by the turbulence of the Earth's atmosphere. We confirm that the larger apertures of the UTs result in higher visibilities for a given separation due to larger overlap of the projected pupils on sky and give predictions for visibility loss as a function of separation to be used for future planning., Comment: 15 pages, 12 figures, 5 tables. Accepted by A&A

Multiplicity is a ubiquitous characteristic of massive stars. Multiple systems offer us a unique observational constraint on the formation of high-mass systems. Herschel 36 A is a massive triple system composed of a close binary (Ab1-Ab2) and an outer component (Aa). We measured the orbital motion of the outer component of Herschel 36 A using infrared interferometry with the AMBER and PIONIER instruments of ESO's Very Large Telescope Interferometer. Our immediate aims are to constrain the masses of all components of this system and to determine if the outer orbit is co-planar with the inner one. Reported spectroscopic data for all three components of this system and our interferometric data allow us to derive full orbital solutions for the outer orbit Aa-Ab and the inner orbit Ab1-Ab2. For the first time, we derive the absolute masses of mAa = 22.3 +/- 1.7 M_sun, mAb1 = 20.5 +/- 1.5 M_sun and mAb2 = 12.5 +/- 0.9 M_sun. Despite not being able to resolve the close binary components, we infer the inclination of their orbit by imposing the same parallax as the outer orbit. Inclinations derived from the inner and outer orbits imply a modest difference of about 22 deg. between the orbital planes. We discuss this result and the formation of Herschel 36 A in the context of Core Accretion and Competitive Accretion models, which make different predictions regarding the statistic of the relative orbital inclinations., Comment: accepted to be published by MNRAS, 8 pages, 5 figures

The stars orbiting the compact radio source Sgr A* in the Galactic Centre are precision probes of the gravitational field around the closest massive black hole. In addition to adaptive optics assisted astrometry (with NACO / VLT) and spectroscopy (with SINFONI / VLT, NIRC2 / Keck and GNIRS / Gemini) over three decades, since 2016/2017 we have obtained 30-100 mu-as astrometry with the four-telescope interferometric beam combiner GRAVITY / VLTI reaching a sensitivity of mK = 20 when combining data from one night. We present the simultaneous detection of several stars within the diffraction limit of a single telescope, illustrating the power of interferometry. The new data for the stars S2, S29, S38 and S55 yield significant accelerations between March and July 2021, as these stars pass the pericenters of their orbits between 2018 and 2023. This allows for a high-precision determination of the gravitational potential around Sgr A*. Our data are in excellent agreement with general relativity orbits around a single central point mass, M = 4.30 x 10^6 M_sun with a precision of about +-0.25%. We improve the significance of our detection of the Schwarzschild precession in the S2 orbit to 7 sigma. Assuming plausible density profiles, an extended mass component inside S2's apocentre (= 0.23" or 2.4 x 10^4 R_S) must be 3000 M_sun (1 sigma), or 0.1% of M. Adding the enclosed mass determinations from 13 stars orbiting Sgr A* at larger radii, the innermost radius at which the excess mass beyond Sgr A* tentatively is seen is r = 2.5" >= 10x the apocentre of S2. This is in full harmony with the stellar mass distribution (including stellar-mass black holes) obtained from the spatially resolved luminosity function., Comment: published in A&A

Stellar orbits at the Galactic Center provide a very clean probe of the gravitational potential of the supermassive black hole. They can be studied with unique precision, beyond the confusion limit of a single telescope, with the near-infrared interferometer GRAVITY. Imaging is essential to search the field for faint, unknown stars on short orbits which potentially could constrain the black hole spin. Furthermore, it provides the starting point for astrometric fitting to derive highly accurate stellar positions. Here, we present $\mathrm{G^R}$, a new imaging tool specifically designed for Galactic Center observations with GRAVITY. The algorithm is based on a Bayesian interpretation of the imaging problem, formulated in the framework of information field theory and building upon existing works in radio-interferometric imaging. Its application to GRAVITY observations from 2021 yields the deepest images to date of the Galactic Center on scales of a few milliarcseconds. The images reveal the complicated source structure within the central $100\,\mathrm{mas}$ around Sgr A*, where we detected the stars S29 and S55 and confirm S62 on its trajectory, slowly approaching Sgr A*. Furthermore, we were able to detect S38, S42, S60, and S63 in a series of exposures for which we offset the fiber from Sgr A*. We provide an update on the orbits of all aforementioned stars. In addition to these known sources, the images also reveal a faint star moving to the west at a high angular velocity. We cannot find any coincidence with any known source and, thus, we refer to the new star as S300. From the flux ratio with S29, we estimate its K-band magnitude as $m_\mathrm{K}\left(\mathrm{S300}\right)\simeq 19.0 - 19.3$. Images obtained with CLEAN confirm the detection., Comment: 24 pages, 14 figures

These protoplanetary disks in T Tauri stars play a central role in star and planet formation. We spatially resolve at sub-au scales the innermost regions of a sample of T Tauri's disks to better understand their morphology and composition. We extended our homogeneous data set of 27 Herbig stars and collected near-IR K-band observations of 17 T Tauri stars, spanning effective temperatures and luminosities in the ranges of ~4000-6000 K and ~0.4-10 Lsun. We focus on the continuum emission and develop semi-physical geometrical models to fit the interferometric data and search for trends between the properties of the disk and the central star. The best-fit models of the disk's inner rim correspond to wide rings. We extend the Radius-luminosity relation toward the smallest luminosities (0.4-10 Lsun) and find the R~L^(1/2) trend is no longer valid, since the K-band sizes measured with GRAVITY are larger than the predicted sizes from sublimation radius computation. No clear correlation between the K-band half-flux radius and the mass accretion rate is seen. Having magnetic truncation radii in agreement with the K-band GRAVITY sizes would require magnetic fields as strong as a few kG, which should have been detected, suggesting that accretion is not the main process governing the location of the half-flux radius of the inner dusty disk. Our measurements agree with models that take into account the scattered light. The N-to-K band size ratio may be a proxy for disentangling disks with silicate features in emission from disks with weak and/or in absorption silicate features. When comparing inclinations and PA of the inner disks to those of the outer disks (ALMA) in nine objects of our sample, we detect misalignments for four objects., Comment: Accepted for publication in A&A; 20 pages, 15 figures, 9 tables

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

The formation and evolution of planetary systems impact the primordial accretion disk. HD141569 is the only known pre-main sequence star characterized by a hybrid disk. Observations probed the outer-disk structure showing a complex system of rings and interferometric observations attempted to characterize its inner 5 au region, but derived limited constraints. The goal of this work was to explore with new high-resolution interferometric observations the properties of the dust and gas in the internal regions of HD141569. We observed HD141569 on mas scales with GRAVITY/VLTI in the near-infrared at low and high spectral resolution. We interpreted the visibilities and spectral energy distribution with geometrical models and radiative transfer techniques to constrain the dust emission. We analyzed the high spectral resolution quantities to investigate the properties of the Br-Gamma line emitting region. Thanks to the combination of three different epochs, GRAVITY resolves the inner dusty disk in the K band. Data modeling shows that an IR excess of about 6% is spatially resolved and that the origin of this emission is confined in a ring of material located at a radius of 1 au from the star with a width smaller than 0.3 au. The MCMax modeling suggests that this emission could originate from a small amount of QHPs, while large silicate grain models cannot reproduce at the same time the observational constraints on the properties of near-IR and mid-IR fluxes. The differential phases in the Br-Gamma line clearly show an S-shape that can be best reproduced witha gas disk in Keplerian rotation, confined within 0.09 au. This is also hinted at by the double-peaked Br-Gamma emission line shape. The modeling of the continuum and gas emission shows that the inclination and position angle of these two components are consistent with a system showing relatively coplanar rings on all scales., Comment: Accepted for publication in A&A; 25 pages, 15 figures, 5 tables

The angular size of the broad line region (BLR) of the nearby active galactic nucleus (AGN) NGC 3783 has been spatially resolved by recent observations with VLTI/GRAVITY. A reverberation mapping (RM) campaign has also recently obtained high quality light curves and measured the linear size of the BLR in a way that is complementary to the GRAVITY measurement. The size and kinematics of the BLR can be better constrained by a joint analysis that combines both GRAVITY and RM data. This, in turn, allows us to obtain the mass of the supermassive black hole in NGC3783 with an accuracy that is about a factor of two better than that inferred from GRAVITY data alone. We derive $M_\mathrm{BH}=2.54_{-0.72}^{+0.90}\times 10^7\,M_\odot$. Finally, and perhaps most notably, we are able to measure a geometric distance to NGC 3783 of $39.9^{+14.5}_{-11.9}$ Mpc. We are able to test the robustness of the BLR-based geometric distance with measurements based on the Tully-Fisher relation and other indirect methods. We find the geometric distance is consistent with other methods within their scatter. We explore the potential of BLR-based geometric distances to directly constrain the Hubble constant, $H_0$, and identify differential phase uncertainties as the current dominant limitation to the $H_0$ measurement precision for individual sources., Comment: 9 pages and 5 figures in main text, Accepted for publication in A&A

[Full abstract in the paper] In recent years, protoplanetary disks with spiral structures have been detected in scattered light, millimeter continuum, and CO gas emission. The mechanisms causing these structures are still under debate. A popular scenario to drive the spiral arms is the one of a planet perturbing the material in the disk. However, if the disk is massive, gravitational instability is usually the favored explanation. Multiwavelength studies could be helpful to distinguish between the two scenarios. So far, only a handful of disks with spiral arms have been observed in both scattered light and millimeter continuum. We aim to perform an in-depth characterization of the protoplanetary disk morphology around WaOph 6 analyzing data obtained at different wavelengths, as well as to investigate the origin of the spiral features in the disk. We present the first near-infrared polarimetric observations of WaOph 6 obtained with SPHERE at the VLT and compare them to archival millimeter continuum ALMA observations. We traced the spiral features in both data sets and estimated the respective pitch angles. We discuss the different scenarios that can give rise to the spiral arms in WaOph 6. We tested the planetary perturber hypothesis by performing hydrodynamical and radiative transfer simulations to compare them with scattered light and millimeter continuum observations., Comment: Accepted for publication in A&A. 17 pages, 9 figures in the main text and 3 figures in the Appendix

Protoplanetary disks drive some of the formation process (e.g., accretion, gas dissipation, formation of structures, etc.) of stars and planets. Understanding such physical processes is one of the main astrophysical questions. HD 163296 is an interesting young stellar object for which infrared and sub-millimeter observations have shown a prominent circumstellar disk with gaps plausibly created by forming planets. This study aims at characterizing the morphology of the inner disk in HD 163296 with multi-epoch near-infrared interferometric observations performed with GRAVITY at the Very Large Telescope Interferometer (VLTI). Our goal is to depict the K-band (lambda_0 ~ 2.2 um) structure of the inner rim with milliarcsecond (sub-au) angular resolution. Our data is complemented with archival PIONIER (H-band; lambda_0 ~ 1.65 um) data of the source. We performed a Gradient Descent parametric model fitting to recover the sub-au morphology of our source. Our analysis shows the existence of an asymmetry in the disk surrounding the central star of HD 163296. We confirm variability of the disk structure in the inner ~2 mas (0.2 au). While variability of the inner disk structure in this source has been suggested by previous interferometric studies, this is the first time that it is confirmed in the H- and K-bands by using a complete analysis of the closure phases and squared visibilities over several epochs. Because of the separation from the star, position changes, and persistence of this asymmetric structure on timescales of several years, we argue that it is a dusty feature (e.g., a vortex or dust clouds), probably, made by a mixing of sillicate and carbon dust and/or refractory grains, inhomogeneously distributed above the mid-plane of the disk., Comment: Accepted to be published in Astronomy and Astrophysics; main-body: 11 pages, 3 figures and 3 tables

We report the time-resolved spectral analysis of a bright near-infrared and moderate X-ray flare of Sgr A*. We obtained light curves in the $M$-, $K$-, and $H$-bands in the mid- and near-infrared and in the $2-8~\mathrm{keV}$ and $2-70~\mathrm{keV}$ bands in the X-ray. The observed spectral slope in the near-infrared band is $\nu L_\nu\propto \nu^{0.5\pm0.2}$; the spectral slope observed in the X-ray band is $\nu L_\nu \propto \nu^{-0.7\pm0.5}$. We tested synchrotron and synchrotron self-Compton (SSC) scenarios. The observed near-infrared brightness and X-ray faintness, together with the observed spectral slopes, pose challenges for all models explored. We rule out a scenario in which the near-infrared emission is synchrotron emission and the X-ray emission is SSC. A one-zone model in which both the near-infrared and X-ray luminosity are produced by SSC and a model in which the luminosity stems from a cooled synchrotron spectrum can explain the flare. In order to describe the mean SED, both models require specific values of the maximum Lorentz factor $\gamma_{max}$, which however differ by roughly two orders of magnitude: the SSC model suggests that electrons are accelerated to $\gamma_{max}\sim 500$, while cooled synchrotron model requires acceleration up to $\gamma_{max}\sim5\times 10^{4}$. The SSC scenario requires electron densities of $10^{10}~\mathrm{cm^{-3}}$ much larger than typical ambient densities in the accretion flow, and thus require in an extraordinary accretion event. In contrast, assuming a source size of $1R_s$, the cooled synchrotron scenario can be realized with densities and magnetic fields comparable with the ambient accretion flow. For both models, the temporal evolution is regulated through the maximum acceleration factor $\gamma_{max}$, implying that sustained particle acceleration is required to explain at least a part of the temporal evolution of the flare., Comment: accepted for publication in Astronomy & Astrophysics; preview abstract shortened due to arXiv requirements

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

GCIRS 7, the brightest star in the Galactic central parsec, formed $6\pm2$ Myr ago together with dozens of massive stars in a disk orbiting the central black-hole. It has been argued that GCIRS 7 is a pulsating body, on the basis of photometric variability. We present the first medium-resolution ($R=500$), K-band spectro-interferometric observations of GCIRS 7, using the GRAVITY instrument with the four auxiliary telescopes of the ESO VLTI. We looked for variations using two epochs, namely 2017 and 2019. We find GCIRS 7 to be moderately resolved with a uniform-disk photospheric diameter of $\theta^*_\text{UD}=1.55 \pm 0.03$ mas ($R^*_\text{UD}=1368 \pm 26$ $R_\odot$) in the K-band continuum. The narrow-band uniform-disk diameter increases above 2.3 $\mu$m, with a clear correlation with the CO band heads in the spectrum. This correlation is aptly modeled by a hot ($T_\text{L}=2368\pm37$ K), geometrically thin molecular shell with a diameter of $\theta_\text{L}=1.74\pm0.03$ mas, as measured in 2017. The shell diameter increased ($\theta_\text{L}=1.89\pm0.03$ mas), while its temperature decreased ($T_\text{L}=2140\pm42$ K) in 2019. In contrast, the photospheric diameter $\theta^*_\text{UD}$ and the extinction up to the photosphere of GCIRS 7 ($A_{\mathrm{K}_\mathrm{S}}=3.18 \pm 0.16$) have the same value within uncertainties at the two epochs. In the context of previous interferometric and photo-spectrometric measurements, the GRAVITY data allow for an interpretation in terms of photospheric pulsations. The photospheric diameter measured in 2017 and 2019 is significantly larger than previously reported using the PIONIER instrument ($\theta_*=1.076 \pm 0.093$ mas in 2013 in the H band). The parameters of the photosphere and molecular shell of GCIRS 7 are comparable to those of other red supergiants that have previously been studied using interferometry., Comment: 12 pages, 11 figures, 3 tables. Accepted for publication in Astronomy and Astrophysics (A&A)

Stars form by accreting material from their surrounding disks. There is a consensus that matter flowing through the disk is channelled onto the stellar surface by the stellar magnetic field. This is thought to be strong enough to truncate the disk close to the so-called corotation radius where the disk rotates at the same rate as the star. Spectro-interferometric studies in young stellar objects show that Hydrogen is mostly emitted in a region of a few milliarcseconds across, usually located within the dust sublimation radius. Its origin is still a matter of debate and it can be interpreted as coming from the stellar magnetosphere, a rotating wind or a disk. In the case of intermediate-mass Herbig AeBe stars, the fact that the Br gamma emission is spatially resolved rules out that most of the emission comes from the magnetosphere. This is due to the weak magnetic fields (some tenths of G) detected in these sources, resulting in very compact magnetospheres. In the case of T Tauri sources, their larger magnetospheres should make them easier to resolve. However, the small angular size of the magnetosphere (a few tenths of milliarcseconds), along with the presence of winds emitting in Hydrogen make the observations interpretation challenging. Here, we present direct evidence of magnetospheric accretion by spatially resolving the inner disk of the 60 pc T Tauri star TW Hydrae through optical long baseline interferometry. We find that the hydrogen near-infrared emission comes from a region approximately 3.5 stellar radii (R*) across. This region is within the continuum dusty disk emitting region (Rcont = 7 R*) and smaller than the corotation radius which is twice as big. This indicates that the hydrogen emission originates at the accretion columns, as expected in magnetospheric accretion models, rather than in a wind emitted at much larger distance (>1au)., Comment: Publish in Nature

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

The young T Tauri star WW Cha was recently proposed to be a close binary object with strong infrared and submillimeter excess associated with circum-system emission. This makes WW Cha a very interesting source for studying the influence of dynamical effects on circumstellar as well as circumbinary material. We derive the relative astrometric positions and flux ratios of the stellar companion in WW Cha from the interferometric model fitting of observations made with the VLTI instruments AMBER, PIONIER, and GRAVITY in the near-infrared from 2011 to 2020. For two epochs, the resulting uv-coverage in spatial frequencies permits us to perform the first image reconstruction of the system in the K band. The positions of nine epochs are used to determine the orbital elements and the total mass of the system. We find the secondary star orbiting the primary with a period of T=206.55 days, a semimajor axis of a=1.01 au, and a relatively high eccentricity of e=0.45. Combining the orbital solution with distance measurements from Gaia DR2 and the analysis of evolutionary tracks, the dynamical mass of Mtot=3.20 Msol can be explained by a mass ratio between ~0.5 and 1. The orbital angular momentum vector is in close alignment with the angular momentum vector of the outer disk as measured by ALMA and SPHERE. The analysis of the relative photometry suggests the presence of infrared excess surviving in the system and likely originating from truncated circumstellar disks. The flux ratio between the two components appears variable, in particular in the K band, and may hint at periods of triggered higher and lower accretion or changes in the disks' structures. The knowledge of the orbital parameters, combined with a relatively short period, makes WW Cha an ideal target for studying the interaction of a close young T Tauri binary with its surrounding material, such as time-dependent accretion phenomena., Comment: Accepted for publication in A&A; 24 pages, 14 figures, 3 tables; affiliations corrected

Using VLTI/GRAVITY and SINFONI data, we investigate the sub-pc gas and dust structure around the nearby type 1 AGN hosted by NGC 3783. The K-band coverage of GRAVITY uniquely allows a simultaneous analysis of the size and kinematics of the broad line region (BLR), the size and structure of the near-IR continuum emitting hot dust, and the size of the coronal line region (CLR). We find the BLR probed through broad Br$\gamma$ emission is well described by a rotating, thick disk with a radial distribution of clouds peaking in the inner region. In our BLR model the physical mean radius of 16 light days is nearly twice the 10 day time lag that would be measured, which matches very well the 10 day time lag that has been measured by reverberation mapping. We measure a hot dust FWHM size of 0.74 mas (0.14 pc) and further reconstruct an image of the hot dust which reveals a faint (5% of the total flux) offset cloud which we interpret as an accreting cloud heated by the central AGN. Finally, we directly measure the FWHM size of the nuclear CLR as traced by the [CaVIII] and narrow Br$\gamma$ line. We find a FWHM size of 2.2 mas (0.4 pc), fully in line with the expectation of the CLR located between the BLR and narrow line region. Combining all of these measurements together with larger scale near-IR integral field unit and mid-IR interferometry data, we are able to comprehensively map the structure and dynamics of gas and dust from 0.01--100 pc., Comment: 20 pages and 18 figures in main text, Accepted for publication in A&A

The GRAVITY instrument on the ESO VLTI pioneers the field of high-precision near-infrared interferometry by providing astrometry at the $10 - 100\,\mu$as level. Measurements at such high precision crucially depend on the control of systematic effects. Here, we investigate how aberrations introduced by small optical imperfections along the path from the telescope to the detector affect the astrometry. We develop an analytical model that describes the impact of such aberrations on the measurement of complex visibilities. Our formalism accounts for pupil-plane and focal-plane aberrations, as well as for the interplay between static and turbulent aberrations, and successfully reproduces calibration measurements of a binary star. The Galactic Center observations with GRAVITY in 2017 and 2018, when both Sgr A* and the star S2 were targeted in a single fiber pointing, are affected by these aberrations at a level of less than 0.5 mas. Removal of these effects brings the measurement in harmony with the dual beam observations of 2019 and 2020, which are not affected by these aberrations. This also resolves the small systematic discrepancies between the derived distance $R_0$ to the Galactic Center reported previously.