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

Cepheid masses continue to be important tests of evolutionary tracks for intermediate mass stars as well as important predictors of their future fate. For systems where the secondary is a B star, {\it Hubble Space Telescope} ultraviolet spectra have been obtained. From these spectra a temperature can be derived, and from this a mass of the companion M$_2$. Once {\it Gaia} DR4 is available, proper motions can be used to determine the inclination of the orbit. Combining mass of the companion, M$_2$, the mass function from the ground-based orbit of the Cepheid and the inclination produces the mass of the Cepheid, M$_1$. The Cepheid system FN Vel is used here to demonstrate this approach and what limits can be put on the Cepheid mass for inclination between 50 and 130$^o$., Comment: accepted by AJ

Wolf-Rayet stars (WRs) are one of the final evolutionary stages of massive stars and immediate progenitors of stellar-mass black holes. Their multiplicity forms an important anchor point in single and binary population models for predicting gravitational-wave progenitors. Recent spectroscopic campaigns have suggested incompatible multiplicity fractions and period distributions for N- and C-rich Galactic WRs (WNs and WCs) at short as well as long orbital periods, in contradiction with evolutionary model predictions. In this work, we employed infrared interferometry using the $K$-band instrument GRAVITY at the VLTI to investigate the multiplicity of WRs at long periods and explore the nature of their companions. We present a survey of 39 Galactic WRs, including 11 WN, 15 WC and 13 H-rich WN (WNh) stars. We detected wide companions with GRAVITY for only four stars: WR 48, WR 89, WR 93 and WR 115. Combining with spectroscopic studies, we arrived at multiplicity fractions of $f^{\rm WN}_{\rm obs} = 0.55\pm0.15$, $f^{\rm WC}_{\rm obs} = 0.40\pm0.13$ and $f^{\rm WNh}_{\rm obs} = 0.23\pm0.12$. In addition, we also found other features in the GRAVITY dataset such as (i) a diffuse extended component in over half the WR sample; (ii) five known spectroscopic binaries resolved in differential phase data and (iii) spatially resolved winds in four stars: WR 16, WR 31a, WR 78 and WR 110. Our survey reveals a lack of intermediate (few 100s d) and long- (few years to decades) period WR systems. The 200-d peak in the period distributions of WR+OB and BH+OB binaries predicted by Case B mass-transfer binary evolution models is not seen in our data. The rich companionship of their O-type progenitors in this separation range suggest that the WR progenitor stars expand and interact with their companions, most likely through unstable mass-transfer, resulting in either a short-period system or a merger., Comment: Accepted in A&A

We resolve the multiple images of the binary-lens microlensing event ASASSN-22av using the GRAVITY instrument of the Very Large Telescope Interferometer (VLTI). The light curves show weak binary perturbations, complicating the analysis, but the joint modeling with the VLTI data breaks several degeneracies, arriving at a strongly favored solution. Thanks to precise measurements of angular Einstein radius \theta_E = 0.724 +/- 0.002 mas and microlens parallax, we determine that the lens system consists of two M dwarfs with masses of M_1 = 0.258 +/- 0.008 M_sun and M_2 = 0.130 +/- 0.007 M_sun, a projected separation of r_\perp = 6.83 +/- 0.31 AU and a distance of D_L = 2.29 +/- 0.08 kpc. The successful VLTI observations of ASASSN-22av open up a new path for studying intermediate-separation (i.e., a few AUs) stellar-mass binaries, including those containing dark compact objects such as neutron stars and stellar-mass black holes., Comment: Accepted by ApJ. See the ancillary file for animation associated with Fig. 8

Interferometric observations of gravitational microlensing events offer an opportunity for precise, efficient, and direct mass and distance measurements of lensing objects, especially those of isolated neutron stars and black holes. However, such observations were previously possible for only a handful of extremely bright events. The recent development of a dual-field interferometer, GRAVITY Wide, has made it possible to reach out to significantly fainter objects, and increase the pool of microlensing events amenable to interferometric observations by two orders of magnitude. Here, we present the first successful observation of a microlensing event with GRAVITY Wide and the resolution of microlensed images in the event OGLE-2023-BLG-0061/KMT-2023-BLG-0496. We measure the angular Einstein radius of the lens with a sub-percent precision, $\theta_{\rm E} = 1.280 \pm 0.009$ mas. Combined with the microlensing parallax detected from the event light curve, the mass and distance to the lens are found to be $0.472 \pm 0.012 M_{\odot}$ and $1.81 \pm 0.05$ kpc, respectively. We present the procedure for the selection of targets for interferometric observations, and discuss possible systematic effects affecting GRAVITY Wide data. This detection demonstrates the capabilities of the new instrument and it opens up completely new possibilities for the follow-up of microlensing events, and future routine discoveries of isolated neutron stars and black holes., Comment: submitted to AAS Journals

Stars with initial masses larger than 8 solar masses undergo substantial mass loss through mechanisms that remain elusive. Unraveling the origins of this mass loss is important for comprehending the evolutionary path of these stars, the type of supernova explosion and whether they become neutron stars or black hole remnants. In 2022 December, RW Cep experienced the Great Dimming in its visible brightness, presenting a unique opportunity to understand mass loss mechanisms. Our previous observations of RW Cep from the CHARA Array, taken during the dimming phase, show a compelling asymmetry in the star images, with a darker zone on the west side of the star indicating presence of dust in front of the star in our line of sight. Here, we present multi-epoch observations from CHARA while the star re-brightened in 2023. We created images using three image reconstruction methods and an analytical model fit. Comparisons of images acquired during the dimming and re-brightening phases reveal remarkable differences. Specifically, the west side of RW Cep, initially obscured during the dimming phase, reappeared during the subsequent re-brightening phase and the measured angular diameter became larger by 8%. We also observed image changes from epoch to epoch while the star is brightening indicating the time evolution of dust in front of the star. We suggest that the dimming of RW Cep was a result from a recent surface mass ejection event, generating a dust cloud that partially obstructed the stellar photosphere., Comment: Accepted by The Astrophysical Journal Letters (ApJL)

The 30 year orbit of the Cepheid Polaris has been followed with observations by the CHARA Array (Center for High Angular Resolution Astronomy) from 2016 through 2021. An additional measurement has been made with speckle interferometry at the Apache Point Observatory. Detection of the companion is complicated by its comparative faintness--an extreme flux ratio. Angular diameter measurements appear to show some variation with pulsation phase. Astrometric positions of the companion were measured with a custom grid-based model-fitting procedure and confirmed with the CANDID software. These positions were combined with the extensive radial velocities discussed by Torres (2023) to fit an orbit. Because of the imbalance of the sizes of the astrometry and radial velocity datasets, several methods of weighting are discussed. The resulting mass of the Cepheid is 5.13$\pm$ 0.28 $M_\odot$. Because of the comparatively large eccentricity of the orbit (0.63), the mass derived is sensitive to the value found for the eccentricity. The mass combined with the distance shows that the Cepheid is more luminous than predicted for this mass from evolutionary tracks. The identification of surface spots is discussed. This would give credence to the identification of photometric variation with a period of approximately 120 days as a rotation period. Polaris has some unusual properties (rapid period change, a phase jump, variable amplitude, unusual polarization). However, a pulsation scenario involving pulsation mode, orbital periastron passage (Torres 2023), and low pulsation amplitude can explain these characteristics within the framework of pulsation seen in Cepheids., Comment: ApJ in press

The Cepheid AW Per is a component in a multiple system with a long period orbit. The radial velocities of Griffin (2016) cover the 38 year orbit well. An extensive program of interferometry with the CHARA array is reported here, from which the long period orbit is determined. In addition, a {\it Hubble Space Telescope} high resolution spectrum in the ultraviolet demonstrates that the companion is itself a binary with nearly equal mass components. These data combined with a distance from {\it Gaia} provide a mass of the Cepheid (primary) of M$_1$ = 6.79 $\pm$ 0.85 $M_\odot$. The combined mass of the secondary is M$_S$ = 8.79 $\pm$ 0.50 $M_\odot$. The accuracy of the mass will be improved after the fourth Gaia data release expected in approximately two years., Comment: Accepted for ApJ

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

Massive stars (those larger than 8 solar masses at formation) have radiative envelopes that cannot sustain a dynamo, the mechanism that produces magnetic fields in lower-mass stars. Despite this, approximately 7\% of massive stars have observed magnetic fields, the origin of which is debated. We used multi-epoch interferometric and spectroscopic observations to characterize HD 148937, a binary system of two massive stars. We found that only one star is magnetic and that it appears younger than its companion. The system properties and a surrounding bipolar nebula can be reproduced with a model in which two stars merged (in a previous triple system) to produce the magnetic massive star. Our results provide observational evidence that magnetic fields form in at least some massive stars through stellar mergers., Comment: Full paper with supplementary materials. 59 pages, 18 figures

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

An inner companion has recently been discovered orbiting the prototype of classical Cepheids, delta Cep, whose orbital parameters are still not fully constrained. We collected new precise radial velocity measurements of delta Cep in 2019 using the HARPS-N spectrograph mounted at the Telescopio Nazionale Galileo. Using these radial velocity measurements, we aimed to improve the orbital parameters of the system. We considered a template available in the literature as a reference for the radial velocity curve of the pulsation of the star. We then calculated the residuals between our global dataset (composed of the new 2019 observations plus data from the literature) and the template as a function of the pulsation phase and the barycentric Julian date. This provides the orbital velocity of the Cepheid component. Using a Bayesian tool, we derived the orbital parameters of the system. Considering priors based on already published Gaia constraints, we find for the orbital period a maximum a posteriori probability of Porb=9.32+/-0.03 years (uncertainties correspond to the 95% highest density probability interval), and we obtain an eccentricity e=0.71+/-0.02, a semimajor axis a=0.029 +/-0.003 arcsecond, and a center-of-mass velocity V0=-17.28+/-0.08 km/s, among other parameters. In this short analysis we derive the orbital parameters of the delta Cep inner binary system and provide a cleaned radial velocity curve of the pulsation of the star, which will be used to study its Baade-Wesselink projection factor in a future publication., Comment: accepted for publication in A&A Letters

Isolated black holes (BHs) and neutron stars (NSs) are largely undetectable across the electromagnetic spectrum. For this reason, our only real prospect of observing these isolated compact remnants is via microlensing; a feat recently performed for the first time. However, characterisation of the microlensing events caused by BHs and NSs is still in its infancy. In this work, we perform N-body simulations to explore the frequency and physical characteristics of microlensing events across the entire sky. Our simulations find that every year we can expect $88_{-6}^{+6}$ BH, $6.8_{-1.6}^{+1.7}$ NS and $20^{+30}_{-20}$ stellar microlensing events which cause an astrometric shift larger than 2~mas. Similarly, we can expect $21_{-3}^{+3}$ BH, $18_{-3}^{+3}$ NS and $7500_{-500}^{+500}$ stellar microlensing events which cause a bump magnitude larger than 1~mag. Leveraging a more comprehensive dynamical model than prior work, we predict the fraction of microlensing events caused by BHs as a function of Einstein time to be smaller than previously thought. Comparison of our microlensing simulations to events in Gaia finds good agreement. Finally, we predict that in the combination of Gaia and GaiaNIR data there will be $14700_{-900}^{+600}$ BH and $1600_{-200}^{+300}$ NS events creating a centroid shift larger than 1~mas and $330_{-120}^{+100}$ BH and $310_{-100}^{+110}$ NS events causing bump magnitudes $> 1$. Of these, $<10$ BH and $5_{-5}^{+10}$ NS events should be detectable using current analysis techniques. These results inform future astrometric mission design, such as GaiaNIR, as they indicate that, compared to stellar events, there are fewer observable BH events than previously thought., Comment: Accepted for publication in MNRAS. 15 pages

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

Previous studies of the protoplanetary disk HD 163296 revealed that the morphology of its sub-au infrared emission encompasses the terminal sublimation front of dust grains, referred to as the inner rim, but also extends into the (supposedly) dust-free region within it. Here, we present a set of radiative hydrostatic simulations of the inner rim in order to assess how much the rim alone can contribute to the observed interferometric visibilities $V$, half-light radii $R_{\mathrm{hl}}$, and fractional disk fluxes $\mathcal{F}$ in the wavelength range $1.5$--$13\,\mu\mathrm{m}$. In our set of models, we regulate the cooling efficiency of the disk via the boundary condition for radiation diffusion and we also modify the shape of the sublimation front. We find that when the cooling efficiency is reduced, the infrared photosphere at the rim becomes hotter, leading to an increase of $R_{\mathrm{hl}}$ sufficient to match the observations. However, the near-infrared disk flux is typically too low ($\mathcal{F}\simeq0.25$ at $1.5\,\mu\mathrm{m}$), resulting in H-band visibility curves located above the observed data. We show that the match to the H-band observations up to moderate baselines can be improved when a wall-shaped rather than curved sublimation front is considered. Nevertheless, our model visibilities always exhibit a bounce at long baselines, which is not observed, confirming the need for additional emission interior to the rim. In summary, our study illustrates how the temperature structure and geometry of the inner rim needs to change in order to boost the rim's infrared emission., Comment: Accepted for publication in AJ, 20 pages, 10 figures

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)

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

Rapid rotation and nonradial pulsations enable Be stars to build decretion disks, where the characteristic line emission forms. A major but unconstrained fraction of Be stars owe their rapid rotation to mass and angular-momentum transfer in a binary. The faint, stripped companions can be helium-burning subdwarf OB-type stars (sdOBs), white dwarfs (WDs), or neutron stars. We present optical/near-IR CHARA interferometry of 37 Be stars selected for spectroscopic indications of low-mass companions. From multi-epoch $H$- and/or $K$-band interferometry plus radial velocities and parallaxes collected elsewhere, we constructed 3D orbits and derived flux ratios and absolute dynamical masses of both components for six objects, quadrupling the number of anchor points for evolutionary models. In addition, a new wider companion was identified for the known Be + sdO binary 59 Cyg, while auxiliary VLTI/GRAVITY spectrointerferometry confirmed circumstellar matter around the sdO companion to HR 2142. On the other hand, we failed to detect any companion to the six Be stars with $\gamma$ Cas-like X-ray emission, with sdOB and main-sequence companions of the expected spectroscopic mass being ruled out for the X-ray-prototypical stars $\gamma$ Cas and $\pi$ Aqr, leaving the elusive WD companions as the most likely companions, as well as a likely explanation of the X-rays. No low-mass main-sequence close companions were identified in the other stars., Comment: Accepted to ApJ

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 report discovering an exoplanet from following up a microlensing event alerted by Gaia. The event Gaia22dkv is toward a disk source rather than the traditional bulge microlensing fields. Our primary analysis yields a Jovian planet with M_p = 0.59^{+0.15}_{-0.05} M_J at a projected orbital separation r_perp = 1.4^{+0.8}_{-0.3} AU, and the host is a ~1.1 M_sun turnoff star at ~1.3 kpc. At r'~14, the host is far brighter than any previously discovered microlensing planet host, opening up the opportunity of testing the microlensing model with radial velocity (RV) observations. RV data can be used to measure the planet's orbital period and eccentricity, and they also enable searching for inner planets of the microlensing cold Jupiter, as expected from the ''inner-outer correlation'' inferred from Kepler and RV discoveries. Furthermore, we show that Gaia astrometric microlensing will not only allow precise measurements of its angular Einstein radius theta_E, but also directly measure the microlens parallax vector and unambiguously break a geometric light-curve degeneracy, leading to definitive characterization of the lens system., Comment: Accepted by AJ

Some evolved binaries, namely post-asymptotic giant branch binaries, are surrounded by stable and massive circumbinary disks similar to protoplanetary disks found around young stars. Around 10% of these disks are transition disks: they have a large inner cavity in the dust. Previous interferometric measurements and modeling have ruled out the cavity being formed by dust sublimation and suggested that the cavity is due to a massive circumbinary planet that traps the dust in the disk and produces the observed depletion of refractory elements on the surface of the post-AGB star. In this study, we test alternative scenario in which the large cavity could be due to dynamical truncation from the inner binary. We performed near-infrared interferometric observations with the CHARA Array on the archetype of such a transition disk around a post-AGB binary: AC Her. We detect the companion at ten epochs over 4 years and determine the 3-dimensional orbit using these astrometric measurements in combination with the radial velocity time series. This is the first astrometric orbit constructed for a post-AGB binary system. We derive the best-fit orbit with a semi-major axis $2.01 \pm 0.01$ mas ($2.83\pm0.08$ au), inclination $(142.9 \pm 1.1)^\circ$ and longitude of the ascending node $(155.1 \pm 1.8)^\circ$. We find that the theoretical dynamical truncation and dust sublimation radius are at least $\sim3\times$ smaller than the observed inner disk radius ($\sim21.5$ mas or 30 au). This strengthens the hypothesis that the origin of such a cavity is due to the presence of a circumbinary planet., Comment: Accepted to be published in The Astrophysical Journal

Many properties of circumbinary discs around evolved post-asymptotic giant branch (post-AGB) binary systems are similar to those of protoplanetary discs around young stars. The deficits of near-infrared (near-IR) flux in the spectral energy distributions (SEDs) of these systems hints towards large dust-free cavities that are reminiscent of transition discs as are commonly observed around young stars. We aim to assess the size of the inner rim of 6 post-AGB binary systems with lack in the near-IR like this. We used resolved mid-infrared (mid-IR) high-angular resolution observations of VLTI/MATISSE and VLTI/MIDI. We compared these inner rim sizes to 5 systems with available MATISSE data that were identified to host a disc starting at the dust sublimation radius. We used geometric ring models to estimate the inner rim sizes, the relative flux contributions of the star, the ring, and an over-resolved emission, the orientation of the ring, and the spectral dependences of the components. We find that the inner dust rims of the targets with a lack of near-IR excess in their SEDs are 2.5 to 7.5 times larger than the theoretical dust sublimation radii, and inner rim sizes of the systems that do not show this deficit are similar to those of their theoretical dust sublimation radii. The physical radii of the inner rims of these transition discs around post-AGB binaries are 3-25 au, which are larger than the disc sizes inferred for transition discs around young stars with VLTI/MIDI. With mid-IR interferometric data, we directly confirm the transition disc nature of six circumbinary discs around post-AGB binary systems. Future observational and modelling efforts are needed to progress in our understanding of the structure, origin, and evolution of these transition discs, Comment: accepted for publication in A&A. 13 pages, including appendices

We aim to measure very precise and accurate model-independent masses and distances of detached binary stars. Precise masses at the $< 1$% level are necessary to test and calibrate stellar interior and evolution models, while precise and independent orbital parallaxes are essential to check for the next Gaia data releases. We combined RV measurements with interferometric observations to determine orbital and physical parameters of ten double-lined spectroscopic systems. We report new relative astrometry from VLTI/GRAVITY and, for some systems, new VLT/UVES spectra to determine the radial velocities of each component. We measured the distance of ten binary systems and the mass of their components with a precision as high as 0.03% (average level 0.2%). They are combined with other stellar parameters (effective temperatures, radii, flux ratios, etc.) to fit stellar isochrones and determine their evolution stage and age. We also compared our orbital parallaxes with Gaia and showed that half of the stars are beyond $1\sigma$ with our orbital parallaxes; although, their RUWE is below the frequently used cutoff of 1.4 for reliable Gaia astrometry. By fitting the telluric features in the GRAVITY spectra, we also estimated the accuracy of the wavelength calibration to be $\sim 0.02$% in high and medium spectral resolution modes. We demonstrate that combining spectroscopic and interferometric observations of binary stars provides extremely precise and accurate dynamical masses and orbital parallaxes. As they are detached binaries, they can be used as benchmark stars to calibrate stellar evolution models and test the Gaia parallaxes., Comment: Accepted for publication in A&A

Inner regions of protoplanetary disks host many complex physical processes such as star-disk interactions, magnetic fields, planet formation, and the migration of new planets. To directly study this region requires milli-arcsecond angular resolution, beyond the diffraction limit of the world's largest optical telescopes and even too small for the mm-wave interferometer ALMA. However, we can use infrared interferometers to image the inner astronomical unit. Here, we present new results from the CHARA and VLTI arrays for the young and luminous Herbig Be star HD 190073. We detect a sub-AU cavity surrounded by a ring-like structure that we interpret as the dust destruction front. We model the shape with 6 radial profiles, 3 symmetric and 3 asymmetric, and present a model-free image reconstruction. All the models are consistent with a near face-on disk with inclination $\lesssim 20^\circ$, and we measure an average ring radius of 1.4 $\pm 0.2$ mas (1.14 AU). Around $48\%$ of the total flux comes from the disk with ~$15\%$ of that emission appearing to emerge from inside the inner rim. The cause of emission is still unclear, perhaps due to different dust grain compositions or gas emission. The skewed models and the imaging point to an off-center star, possibly due to binarity. Our image shows a sub-AU structure, which seems to move between the two epochs inconsistently with Keplerian motion and we discuss possible explanations for this apparent change., Comment: 20 pages, 9 figures

The study of the multiplicity of massive stars gives hints on their formation processes and their evolutionary paths, which are still not fully understood. Large separation binaries (>50 milliseconds of arc, mas) can be probed by adaptive-optics-assisted direct imaging and sparse aperture masking, while close binaries can be resolved by photometry and spectroscopy. However, optical long baseline interferometry is mandatory to establish the multiplicity of Galactic massive stars at the separation gap between 1 and 50 mas. In this paper, we aim to demonstrate the capability of the new interferometric instrument MIRC-X, located at the CHARA Array, to study the multiplicity of O-type stars and therefore probe the full range of separation for more than 120 massive stars (H<7.5 mag). We initiated a pilot survey of bright O-type stars (H<6.5mag) observable with MIRC-X. We observed 29 O-type stars, including two systems in average atmospheric conditions around a magnitude of H=7.5 mag. We systematically reduced the obtained data with the public reduction pipeline of the instrument. We analyzed the reduced data using the dedicated python software CANDID to detect companions. Out of these 29 systems, we resolved 19 companions in 17 different systems with angular separations between ~0.5 and 50 mas. This results in a multiplicity fraction fm=17/29=0.59+/-0.09, and an average number of companions fc=19/29=0.66+/-0.13. Those results are in agreement with the results of the SMASH+ survey in the Southern Hemisphere. Thirteen of these companions have been resolved for the first time, including the companion responsible for the nonthermal emission in Cyg OB2-5 A and the confirmation of the candidate companion of HD 47129 suggested by SMASH+. A large survey on more than 120 northern O-type stars (H<7.5) is possible with MIRC-X and will be fruitful., Comment: 15 pages, 9 figures, 5 tables, accepted in A&A

The Baade-Wesselink (BW) method of distance determination of Cepheids is used to calibrate the distance scale. Various versions of this method are mainly based on interferometry and/or the surface-brightness color relation (SBCR). We quantify the impact of the SBCR, its slope, and its zeropoint on the projection factor. This quantity is used to convert the pulsation velocity into the radial velocity in the BW method. We also study the impact of extinction and of a potential circumstellar environment on the projection factor. We analyzed HARPS-N spectra of eta Aql to derive its radial velocity curve using different methods. We then applied the inverse BW method using various SBCRs in the literature in order to derive the BW projection factor. We find that the choice of the SBCR is critical: a scatter of about 8% is found in the projection factor for different SBCRs in the literature. The uncertainty on the coefficients of the SBCR affects the statistical precision of the projection factor only little (1-2\%). Confirming previous studies, we find that the method with which the radial velocity curve is derived is also critical, with a potential difference on the projection factor of 9%. An increase of 0.1 in E(B-V) translates into a decrease in the projection factor of 3%. A 0.1 magnitude effect of a circumstellar envelope (CSE) in the visible domain is rather small on the projection factor, about 1.5%. However, we find that a 0.1 mag infrared excess in the K band due to a CSE can increase the projection factor by about 6%. The impact of the surface-brightness color relation on the BW projection factor is found to be critical. Efforts should be devoted in the future to improve the SBCR of Cepheids empirically, but also theoretically, taking their CSE into account as well., Comment: Accepted for publication in Astronomy & Astrophysics

The GRAVITY instrument has been revolutionary for near-infrared interferometry by pushing sensitivity and precision to previously unknown limits. With the upgrade of GRAVITY and the Very Large Telescope Interferometer (VLTI) in GRAVITY+, these limits will be pushed even further, with vastly improved sky coverage, as well as faint-science and high-contrast capabilities. This upgrade includes the implementation of wide-field off-axis fringe-tracking, new adaptive optics systems on all Unit Telescopes, and laser guide stars in an upgraded facility. GRAVITY+ will open up the sky to the measurement of black hole masses across cosmic time in hundreds of active galactic nuclei, use the faint stars in the Galactic centre to probe General Relativity, and enable the characterisation of dozens of young exoplanets to study their formation, bearing the promise of another scientific revolution to come at the VLTI., Comment: Published in the ESO Messenger

Kernel phase imaging (KPI) enables the direct detection of substellar companions and circumstellar dust close to and below the classical (Rayleigh) diffraction limit. We present a kernel phase analysis of JWST NIRISS full pupil images taken during the instrument commissioning and compare the performance to closely related NIRISS aperture masking interferometry (AMI) observations. For this purpose, we develop and make publicly available the custom "Kpi3Pipeline" enabling the extraction of kernel phase observables from JWST images. The extracted observables are saved into a new and versatile kernel phase FITS file (KPFITS) data exchange format. Furthermore, we present our new and publicly available "fouriever" toolkit which can be used to search for companions and derive detection limits from KPI, AMI, and long-baseline interferometry observations while accounting for correlated uncertainties in the model fitting process. Among the four KPI targets that were observed during NIRISS instrument commissioning, we discover a low-contrast (~1:5) close-in (~1 $\lambda/D$) companion candidate around CPD-66~562 and a new high-contrast (~1:170) detection separated by ~1.5 $\lambda/D$ from 2MASS~J062802.01-663738.0. The 5-$\sigma$ companion detection limits around the other two targets reach ~6.5 mag at ~200 mas and ~7 mag at ~400 mas. Comparing these limits to those obtained from the NIRISS AMI commissioning observations, we find that KPI and AMI perform similar in the same amount of observing time. Due to its 5.6 times higher throughput if compared to AMI, KPI is beneficial for observing faint targets and superior to AMI at separations >325 mas. At very small separations (<100 mas) and between ~250-325 mas, AMI slightly outperforms KPI which suffers from increased photon noise from the core and the first Airy ring of the point-spread function., Comment: 34 pages, 17 figures, accepted for publication in PASP

Because many classical Be stars may owe their nature to mass and angular-momentum transfer in a close binary, the present masses, temperatures, and radii of their components are of high interest for comparison to stellar evolution models. Kappa Dra is a 61.5-day single-lined binary with a B6 IIIe primary. With the CHARA Array instruments MIRC/MIRC-X and MYSTIC, we detected the secondary at (approximately photospheric) flux ratios of 1.49 +- 0.10% and 1.63 +- 0.09% in the H and K band, respectively. From a large and diverse optical spectroscopic database only the radial velocity curve of the Be star could be extracted. However, employing the parallaxes from Hipparcos and Gaia, which agree within their nominal 1-sigma errors, we could derive the total mass and found component masses of 3.65 +- 0.48 Msun and 0.426 +- 0.043 Msun for the Be star and the companion, respectively. Previous cross-correlation of the observed far-UV spectrum with sdO spectral model templates had not detected a companion belonging to the hot O-type subdwarf (sdO) population known from ~20 earlier-type Be stars. Guided by our full 3D orbital solution, we found a strong cross-correlation signal for a stripped subdwarf B-type companion (far-UV flux ratio of 2.3 +- 0.5%), enabling the first firm characterization of such a star, and making kappa Dra the first mid- to late-type Be star with a directly-observed subdwarf companion., Comment: Accepted to ApJ. Figure sets available on request

Achernar, the closest and brightest classical Be star, presents rotational flattening, gravity darkening, occasional emission lines due to a gaseous disk, and an extended polar wind. It is also a member of a close binary system with an early A-type dwarf companion. We aim to determine the orbital parameters of the Achernar system and to estimate the physical properties of the components. We monitored the relative position of Achernar B using a broad range of high angular resolution instruments of the VLT/VLTI (VISIR, NACO, SPHERE, AMBER, PIONIER, GRAVITY, and MATISSE) over a period of 13 years (2006-2019). These astrometric observations are complemented with a series of more than 700 optical spectra for the period from 2003 to 2016. We determine that Achernar B orbits the Be star on a seven-year period, eccentric orbit (e = 0.7255 +/- 0.0014) which brings the two stars within 2 au at periastron. The mass of the Be star is found to be mA = 6.0 +/- 0.6 Msun for a secondary mass of mB = 2.0 +/- 0.1 Msun. We find a good agreement of the parameters of Achernar A with the evolutionary model of a critically rotating star of 6.4 Msun at an age of 63 million years. We also identify a resolved comoving low-mass star, which leads us to propose that Achernar is a member of the Tucana-Horologium moving group. Achernar A is presently in a short-lived phase of its evolution following the turn-off, during which its geometrical flattening ratio is the most extreme. Considering the orbital parameters, no significant interaction occurred between the two components, demonstrating that Be stars may form through a direct, single-star evolution path without mass transfer. Since component A will enter the instability strip in a few hundred thousand years, Achernar appears to be a promising progenitor of the Cepheid binary systems., Comment: 27 pages, Astronomy & Astrophysics, in press. The full set of continuum normalized high resolution spectra of Achernar is available at https://doi.org/10.5281/zenodo.6977303

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

Despite image reconstruction becoming more widespread when interpreting OIFITS Data, model fitting in u,v space often remains the best way to interpret data, either because of the sparsity of the data, or because a quantitative measurement needs to be done. PMOIRED, is a flexible Python library to visualize, manipulate and model OIFITS data using simple geometric models. The strength of PMOIRED resides in its capability to combine linearly various simple components to create complex scenes, while linking, constraining, and adding priors to fitted parameters. The code also enables grid search to find global minima, as well as data resampling to better evaluate uncertainties. In addition to analytical functions, arbitrary radial profiles, azimuthal variations or sparse wavelet modelling of spectra are implemented., Comment: 8 pages, 3 figures. Proceeding of SPIE Telescopes+Instrumentation 2022. Code referenced at http://ascl.net/code/v/3286

The fate of a massive star during the latest stages of its evolution is highly dependent on its mass-loss history and geometry, with the yellow hypergiants being key objects to study those phases of evolution. We present near-IR interferometric observations of the famous yellow hypergiant IRC +10420 and blue spectra taken between 1994-2019. Our 2.2 $\mu$m GRAVITY/VLTI observations attain a spatial resolution of $\sim$5 stellar radii and probe the hot emission in the K-band tracing the gas via Na i double emission and the Br$\gamma$ emission. The observed configurations spatially resolve the 2.2 $\mu$m continuum as well as the Br$\gamma$ and the Na i emission lines. Our geometric modelling demonstrates the presence of a compact neutral zone (Na i) which is slightly larger than the continuum but within an extended Br$\gamma$ emitting region. Our geometric models of the Br$\gamma$ emission confirm an hour-glass geometry of the wind. To explain this peculiar geometry we investigate the presence of a companion at 7-800 au separations and find no signature at the contrast limit of our observations (3.7 mag at 3$\sigma$). We report an evolution of the ejecta over a time span of 7 years, which allows us to constrain the opening angle of the hour-glass geometry at $<$10$^\circ$. Lastly, we present the first blue optical spectra of IRC +10420 since 1994. The multi-epoch data indicate that the spectral type, and thus temperature, of the object has essentially remained constant during the intervening years. This confirms earlier conclusions that following an increase in temperature of 2000 K in less than two decades prior to 1994, the temperature increase has halted. This suggests that this yellow hypergiant has "hit" the White Wall in the HR-diagram preventing it from evolving blue-wards, and will likely undergo a major mass-loss event in the near future., Comment: Accepted for publication in MNRAS; 13 pages, 7 figures, 2 tables

Context: Stellar evolution models are highly dependent on accurate mass estimates, especially for high-mass stars in the early stages of evolution. The most direct method for obtaining model-independent masses is derivation from the orbit of close binaries. Aims: To derive the first astrometric+RV orbit solution for the single-lined spectroscopic binary MWC 166 A, based on CHARA and VLTI near-infrared interferometry over multiple epochs and ~100 archival radial velocity measurements, and to derive fundamental stellar parameters from this orbit. We also sought to model circumstellar activity in the system from K-band spectral lines. Methods: We geometrically modelled the dust continuum to derive astrometry at 13 epochs and constrain individual stellar parameters. We used the continuum models as a base to examine differential phases, visibilities and closure phases over the Br-$\gamma$ and He-I emission lines. Results: Our orbit solution suggests a period of $367.7\pm0.1$ d, twice as long as found with previous RV orbit fits, subsequently constraining the component masses to $M_1=12.2\pm2.2 M_\odot$ and $M_2=4.9\pm0.5 M_\odot$. The line-emitting gas was found to be localised around the primary and is spatially resolved on scales of ~11 stellar radii, with the spatial displacement between the line wings consistent with a rotating disc. Conclusions: The large radius and stable orientation of the line emission are inconsistent with magnetospheric or boundary-layer accretion, but indicate an ionised inner gas disk around MWC 166 Aa. We observe line variability that could be explained either with generic line variability in a Herbig star disc or V/R variations in a decretion disc. We also constrained the age of the system to ~$(7\pm2)\times10^5$ yr, consistent with the system being comprised of a main-sequence primary and a secondary still contracting towards the main sequence., Comment: 24 pages, 19 figures, 7 tables, 1 appendix. Accepted in A&A

While the number of stellar-mass black holes detected in X-rays or as gravitational wave sources is steadily increasing, the known population remains orders of magnitude smaller than predicted by stellar evolution theory. A significant fraction of stellar-mass black holes is expected to hide in X-ray-quiet binaries where they are paired with a "normal" star. Although a handful of such quiescent black hole candidates have been proposed, the majority have been challenged by follow-up investigations. A confusion that emerged recently concerns binary systems that appear to contain a normal B-type star with an unseen companion, believed to be a black hole. On closer inspection, some of these seemingly normal B-type stars instead turn out to be stars stripped of most of their mass through an interaction with their binary companion, which in at least two cases is a rapidly rotating star rather than a compact object. These contaminants in the search for quiescent black holes are themselves extremely interesting objects as they represent a rare phase of binary evolution, and should be given special attention when searching for binaries hosting black holes in large spectroscopic studies., Comment: 7 pages, 5 figures

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

Spectroscopic multiplicity surveys of O stars in young clusters and OB associations have revealed that a large portion ($\sim$ 70%) of these massive stars (M$_{i}$ $\gt$ 15 $M_{\odot}$) belong to close and short-period binaries (physical separation d $\lt$few au). Despite the recent and significant progress, the formation mechanisms leading to such close massive multiple systems remain to be elucidated. As a result, young massive close binaries (or higher-order multiple systems) are unique laboratories to figure out the pairing mechanism of high-mass stars. We present the first VLTI/GRAVITY observations of six young O stars in the M17 star-forming region ($\lesssim$ 1 Myr) and two additional foreground stars. From the interferometric model fitting of visibility amplitudes and closure phases, we search for companions and measure their positions and flux ratios. Combining the resulting magnitude difference with atmosphere models and evolutionary tracks, we further constrain the masses of the individual components. All of the six high-mass stars are in multiple systems, leading to a multiplicity fraction (MF) of 100%, yielding a 68% confidence interval of 94-100%. We detect a total number of 9 companions with separations up to 120 au. Including previously identified spectroscopic companions, the companion fraction of the young O stars in our sample reaches 2.3$\pm$0.6. The derived masses span a wide range from 2.5 to 50 $M_{\odot}$, with a great tendency towards high-mass companions. While based on a modest sample, our results clearly indicate that the origin of the high degree of multiplicity is rooted in their star formation mechanism. No clear evidence for one of the competing concepts of massive star formation (core accretion or competitive accretion) could be found. However, our results are compatible with migration as a scenario for the formation of close massive binaries., Comment: 20 pages, 14 figures, 4 tables, 3 appendices

Two scenarios have been proposed to match the existing observational constraints of the object HR 6819. The system could consist of a close inner B-type giant plus a black hole (BH) binary with an additional Be companion in a wide orbit. Alternatively, it could be a binary composed of a stripped B star and a Be star in a close orbit. Either scenario makes HR 6819 a cornerstone object as the stellar BH closest to Earth, or as an example of an important transitional, non-equilibrium phase for Be stars with solid evidence for its nature. We aimed to distinguish between the two scenarios for HR 6819. Both models predict two luminous stars but with very different angular separations and orbital motions. Therefore, the presence of bright sources in the 1-100 milliarcsec (mas) regime is a key diagnostic for determining the nature of the HR 6819 system. We obtained new high-angular resolution data with VLT/MUSE and VLTI/GRAVITY of HR 6819. The MUSE data are sensitive to bright companions at large scales, whilst the interferometric GRAVITY data are sensitive down to separations on mas scales and large magnitude differences. The MUSE observations reveal no bright companion at large separations and the GRAVITY observations indicate the presence of a stellar companion at an angular separation of ~1.2 mas that moves on the plane of the sky over a timescale compatible with the known spectroscopic 40-day period. We conclude that HR 6819 is a binary system and that no BH is present in the system. The unique nature of HR 6819, and its proximity to Earth make it an ideal system for quantitatively characterising the immediate outcome of binary interaction and probing how Be stars form., Comment: 16 pages, 11 figures, subject of ESO press release eso2204

Classical Be stars are possible products of close binary evolution, in which the mass donor becomes a hot, stripped O or B-type subdwarf (sdO/sdB), and the mass gainer spins up and grows a disk to become a Be star. While several Be+sdO binaries have been identified, dynamical masses and other fundamental parameters are available only for a single Be+sdO system, limiting the confrontation with binary evolution models. In this work, we present direct interferometric detections of the sdO companions of three Be stars 28 Cyg, V2119 Cyg, and 60 Cyg, all of which were previously found in UV spectra. For two of the three Be+sdO systems, we present first orbits and preliminary dynamical masses of the components, revealing that one of them could be the first identified progenitor of a Be/X-ray binary with a neutron star companion. These results provide new sets of fundamental parameters that are crucially needed to establish the evolutionary status and origin of Be stars., Comment: accepted to ApJ

For several transition disks (TDs), dark regions interpreted as shadows have been observed in scattered light imaging and are hypothesized to originate from misalignments between distinct disk regions. We aim to investigate the presence of misalignments in TDs. We study the inner disk geometries of 20 well-known transition disks with VLTI/GRAVITY observations and use complementary $^{12}$CO and $^{13}$CO molecular line data from ALMA to derive the orientation of the outer disk regions. We fit simple models to the GRAVITY data to derive the inner disks inclination and position angles. The outer disk geometries were derived from Keplerian fits to the ALMA velocity maps and compared to the inner disk constraints. We also predicted the locations of shadows for significantly misaligned systems. Our analysis reveals six disks to exhibit significant misalignments between their inner and outer disks. The predicted shadow positions agree well with the scattered light images of HD100453 and HD142527, and we find supporting evidence for a shadow in the disk around CQ Tau. In the other three targets for which we infer significantly misaligned disks, V1247 Ori, V1366 Ori, and RY Lup, we do not see any evident sign of shadows in the scattered light images. The scattered light shadows observed in DoAr44, HD135344B, and HD139614 are consistent with our observations, yet the underlying morphology is likely too complex to be described by our models and the accuracy achieved by our observations. Whereas we can derive precise constraints on the potential shadow positions for well-resolved inner disks around HAeBe stars, the statistical uncertainties for the marginally resolved inner disks around the TTS of our sample make it difficult to extract conclusive constraints for the presence of shadows in these systems., Comment: Accepted for publication in Astronomy \& Astrophysics

Context: FS Canis Majoris (FS CMa, HD 45677) is an unclassified B[e] star surrounded by an inclined dust disk. The evolutionary stage of FS CMa is still debated. Perpendicular to the circumstellar disk, a bipolar outflow was detected. Infrared aperture-synthesis imaging provides us with a unique opportunity to study the disk structure. Aims: Our aim is to study the intensity distribution of the disk of FS CMa in the mid-infrared L and N bands. Methods: We performed aperture-synthesis imaging of FS CMa with the MATISSE instrument (Multi AperTure mid-Infrared SpectroScopic Experiment) in the low spectral resolution mode to obtain images in the L and N bands. We computed radiative transfer models that reproduce the L- and N-band intensity distributions of the resolved disks. Results: We present L- and N-band aperture-synthesis images of FS CMa reconstructed in the wavelength bands of 3.4-3.8 and 8.6-9.0 micrometer. In the L-band image, the inner rim region of an inclined circumstellar disk and the central object can be seen with a spatial resolution of 2.7 milliarcsec (mas). An inner disk cavity with an angular diameter of 6x12mas is resolved. The L-band disk consists of a bright northwestern (NW) disk region and a much fainter southeastern (SE) region. The images suggest that we are looking at the bright inner wall of the NW disk rim, which is on the far side of the disk. In the N band, only the bright NW disk region is seen. In addition to deriving the inclination and the inner disk radius, fitting the reconstructed brightness distributions via radiative transfer modeling allows one to constrain the innermost disk structure, in particular the shape of the inner disk rim., Comment: accepted for publication in A&A