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

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.726 +/- 0.002 mas and microlens parallax, we determine that the lens system consists of two M dwarfs with masses of M_1 = 0.261 +/- 0.009 M_sun and M_2 = 0.252 +/- 0.017 M_sun, a projected separation of r_\perp = 7.42 +/- 0.33 AU and a distance of D_L = 2.31 +/- 0.09 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: see the ancillary file for animation associated with Fig. 8

Studying the orbital motion of stars around Sagittarius A* in the Galactic Center provides a unique opportunity to probe the gravitational potential near the supermassive black hole at the heart of our Galaxy. Interferometric data obtained with the GRAVITY instrument at the Very Large Telescope Interferometer (VLTI) since 2016 has allowed us to achieve unprecedented precision in tracking the orbits of these stars. GRAVITY data have been key to detecting the in-plane, prograde Schwarzschild precession of the orbit of the star S2, as predicted by General Relativity. By combining astrometric and spectroscopic data from multiple stars, including S2, S29, S38, and S55 - for which we have data around their time of pericenter passage with GRAVITY - we can now strengthen the significance of this detection to an approximately $10 \sigma$ confidence level. The prograde precession of S2's orbit provides valuable insights into the potential presence of an extended mass distribution surrounding Sagittarius A*, which could consist of a dynamically relaxed stellar cusp comprised of old stars and stellar remnants, along with a possible dark matter spike. Our analysis, based on two plausible density profiles - a power-law and a Plummer profile - constrains the enclosed mass within the orbit of S2 to be consistent with zero, establishing an upper limit of approximately $1200 \, M_\odot$ with a $1 \sigma$ confidence level. This significantly improves our constraints on the mass distribution in the Galactic Center. Our upper limit is very close to the expected value from numerical simulations for a stellar cusp in the Galactic Center, leaving little room for a significant enhancement of dark matter density near Sagittarius A*., Comment: Submitted to A&A on September 17, 2024

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

We present the Wavefront Sensor units of the Gravity Plus Adaptive Optics (GPAO) system, which will equip all 8m class telescopes of the VLTI and is an instrumental part of the GRAVITY+ project. It includes two modules for each Wavefront Sensor unit: a Natural Guide Star sensor with high-order 40x40 Shack-Hartmann and a Laser Guide Star 30x30 sensor. The state-of-the-art AO correction will considerably improve the performance for interferometry, in particular high-contrast observations for NGS observations and all-sky coverage with LGS, which will be implemented for the first time on VLTI instruments. In the following, we give an overview of the Wavefront Sensor units system after completion of their integration and characterization., Comment: Optical and Infrared Interferometry and Imaging IX

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.

TeNPy (short for 'Tensor Network Python') is a python library for the simulation of strongly correlated quantum systems with tensor networks. The philosophy of this library is to achieve a balance of readability and usability for new-comers, while at the same time providing powerful algorithms for experts. The focus is on MPS algorithms for 1D and 2D lattices, such as DMRG ground state search, as well as dynamics using TEBD, TDVP, or MPO evolution. This article is a companion to the recent version 1.0 release of TeNPy and gives a brief overview of the package., Comment: v2: updated funding acknowledgement

Transient grating spectroscopy has become a mainstay among metal and semiconductor characterization techniques. Here we extend the technique towards the shortest achievable timescales by using tabletop high-harmonic generation of attosecond extreme ultraviolet (XUV) pulses that diffract from transient gratings generated with 500-1000 nm sub-5 fs near-infrared (NIR) pulses. We demonstrate the power of attosecond transient grating spectroscopy (ATGS) by investigating the ultrafast photoexcited dynamics in an Sb semimetal thin film. ATGS provides an element-specific, background-free signal, unfettered by spectral congestion, in contrast to transient absorption spectroscopy. With the ATGS measurements in Sb polycrystalline thin films, we observe the generation of coherent phonons and investigate the lattice and carrier dynamics. Among the latter processes, we extract carrier thermalization, hot carrier cooling, and electron-hole recombination, which are on the order of 20 fs, 50 fs, and 2 ps, timescales, respectively. Furthermore, simultaneous collection of transient absorption and transient grating data allows us to extract the total complex dielectric constant in the sample dynamics, including the real-valued refractive index, from which we are also able to investigate carrier-phonon interactions and longer-lived phonon dynamics. The outlined experimental technique expands the capabilities of transient grating spectroscopy and attosecond spectroscopies by providing a wealth of information concerning carrier and lattice dynamics with an element-selective technique, at the shortest achievable timescales.

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

A systematic study, based on the third-moment structure function, of Sgr A*'s variability finds an exponential rise time $\tau_{1,\rm{obs}}=14.8^{+0.4}_{-1.5}~\mathrm{minutes}$ and decay time $\tau_{2,\rm{obs}}=13.1^{+1.3}_{-1.4}~\mathrm{minutes}$. This symmetry of the flux-density variability is consistent with earlier work, and we interpret it as caused by the dominance of Doppler boosting, as opposed to gravitational lensing, in Sgr~A*'s light curve. A relativistic, semi-physical model of Sgr~A* confirms an inclination angle $i<45$ degrees. The model also shows that the emission of the intrinsic radiative process can have some asymmetry even though the observed emission does not. The third-moment structure function, which is a measure of the skewness of the light-curve increments, may be a useful summary statistic in other contexts of astronomy because it senses only temporal asymmetry, i.e., it averages to zero for any temporally symmetric signal., Comment: Accepted for publication in A&A letters

The peaceable queens problem asks to determine the maximum number $a(n)$ such that there is a placement of $a(n)$ white queens and $a(n)$ black queens on an $n \times n$ chessboard so that no queen can capture any queen of the opposite color. In this paper, we consider the peaceable queens problem and its variant on the toroidal board. For the regular board, we show that $a(n) \leq 0.1716n^2$, for all sufficiently large $n$. This improves on the bound $a(n) \leq 0.25n^2$ of van Bommel and MacEachern. For the toroidal board, we provide new upper and lower bounds. Somewhat surprisingly, our bounds show that there is a sharp contrast in behaviour between the odd torus and the even torus. Our lower bounds are given by explicit constructions. For the upper bounds, we formulate the problem as a non-linear optimization problem with at most $100$ variables, regardless of the size of the board. We solve our non-linear program exactly using modern optimization software. We also provide a local search algorithm and a software implementation which converges very rapidly to solutions which appear optimal. Our algorithm is sufficiently robust that it works on both the regular and toroidal boards. For example, for the regular board, the algorithm quickly finds the so-called Ainley construction. Thus, our work provides some further evidence that the Ainley construction is indeed optimal., Comment: 38 pages with 7 figures. Main body is 13 pages. Appendix A contains our best odd torus examples for all odd board sizes between 13 and 63

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

Understanding initial electron thermalization has relevance to both fundamental scientific knowledge and application to the construction of novel devices. In this study, attosecond transient absorption is used to directly measure initial electron thermalization times of 38 $\pm$ 8 fs, 15 $\pm$ 3 fs, 4.2 $\pm$ 1 fs, and 2.0 $\pm$ 0.3 fs for Mg, Pt, Fe, and Co, respectively. Through time dependent density function theory calculations, it is shown that the fast electron thermalization observed in Fe and Co is correlated with a strong local field effect. We find that a simple analytical model can be used to calculate the initial electron thermalization time measured by the transient extreme ultraviolet absorption spectroscopy method performed here. Our results suggest that the most significant contributions to the initial electron thermalization times are the basic metal properties of the density of states volume available for scattering and screened electron interaction. Many-body effects contribute less, but still significantly to the initial electron thermalization time. Ultimately the information gained through this study shows the unique view that attosecond transient absorption spectroscopy contributes to unraveling and monitoring electron dynamics and its connection to many-body effects in metals and beyond., Comment: 32 pages, 13 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

We investigate the dynamics of an ion moving through a homogeneous Bose-Einstein condensate (BEC) after an initial momentum is imparted. For this, we derive a master equation in the weak-coupling limit and Lamb-Dicke approximation for the reduced density matrix of the ion. We study the time evolution of the ion's kinetic energy and observe that its expectation value, identified as the ion temperature $T_\mathrm{ion}$, is reduced by several orders of magnitude in a time on the order of microseconds for a condensate density in the experimentally relevant range between $10^{13}\,\mathrm{cm}^{-3}$ and $10^{14}\,\mathrm{cm}^{-3}$. We characterize this behavior by defining the duration at half maximum as the time required by $T_\mathrm{ion}$ to reach half of its initial value, and study its dependence on the system parameters. Similarly, we find that the expectation value of the ion's momentum operator is reduced by nine orders of magnitude on the same timescale, making the ion's position converge to a final value. Based on these results, we conclude that the interaction with the bosonic bath allows for cooling and pinning of the ion by decreasing the expectation value of its kinetic energy and velocity, which constitutes a result of direct relevance for current atom-ion experiments., Comment: 16 pages, 11 figures

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

The functional determinant approach (FDA) is a simple method to compute exactly certain observables for ideal quantum systems and has been successfully applied to the Fermi polaron problem to obtain the dynamical overlap and spectral function. Unfortunately, its application to Bosonic ultracold gases is prohibited by the failure of the grand canonical ensemble for these systems. In this paper, we show how to circumvent this problem and develop a Bosonic FDA. This yields exact injection and ejection spectra for ideal Bose polarons at arbitrary temperatures. While coherent features visible at absolute zero quickly smear out with rising temperature as expected, the line width of the main peak is, counterintuitively, found to decrease near unitarity. Furthermore, we provide explicit formulas for the overlap operator, which allow to compute the necessary determinants for both Bose and Fermi polarons more efficiently than previously possible., Comment: 10 pages, 3 figures

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

Few-shot prompting and step-by-step reasoning have enhanced the capabilities of Large Language Models (LLMs) in tackling complex tasks including code generation. In this paper, we introduce a prompt selection and augmentation algorithm aimed at improving mathematical reasoning and robot arm operations. Our approach incorporates a multi-stage example augmentation scheme combined with an example selection scheme. This algorithm improves LLM performance by selecting a set of examples that increase diversity, minimize redundancy, and increase relevance to the question. When combined with the Program-of-Thought prompting, our algorithm demonstrates an improvement in performance on the GSM8K and SVAMP benchmarks, with increases of 0.3% and 1.1% respectively. Furthermore, in simulated tabletop environments, our algorithm surpasses the Code-as-Policies approach by achieving a 3.4% increase in successful task completions and a decrease of over 70% in the number of examples used. Its ability to discard examples that contribute little to solving the problem reduces the inferencing time of an LLM-powered robotics system. This algorithm also offers important benefits for industrial process automation by streamlining the development and deployment process, reducing manual programming effort, and enhancing code reusability., Comment: 17 pages, 4 figures

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

Aims. As part of the ongoing GRAVITY+ upgrade of the Very Large Telescope Interferometer infrastructure, we aim to improve the performance of the GRAVITY Fringe-Tracker, and to enable its use by other instruments. Methods. We modify the group delay controller to consistently maintain tracking in the white light fringe, characterised by a minimum group delay. Additionally, we introduce a novel approach in which fringe-tracking is performed in the non-observable Optical Path Length state-space, using a covariance-weighted Kalman filter and an auto-regressive model of the disturbance. We outline this new state-space representation, and the formalism we use to propagate the state-vector and generate the control signal. While our approach is presented specifically in the context of GRAVITY/GRAVITY+, it can easily be adapted to other instruments or interferometric facilities. Results. We successfully demonstrate phase delay tracking within a single fringe, with any spurious phase jumps detected and corrected in less than 100 ms. We also report a significant performance improvement, as evidenced by a reduction of about 30 to 40% in phase residuals, and a much better behaviour under sub-optimal atmospheric conditions. Compared to what was observed in 2019, the median residuals have decreased from 150 nm to 100 nm on the Auxiliary Telescopes and from 250 nm to 150 nm on the Unit Telescopes. Conclusions. The improved phase-delay tracking combined with whit light fringe tracking means that from now-on, the GRAVITY Fringe-Tracker can be used by other instruments operating in different wavebands. The only limitation remains the need for an optical path dispersion adjustment., Comment: 16 pages, 8 figures. Accepted for publication in 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)

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

Attosecond extreme-ultraviolet (XUV) transient absorption spectroscopy measurements on the Peierls-distorted phase of the semimetal antimony (Sb) are presented. After excitation by an ultrashort, broad band near-infrared (NIR) pulse, the distortion is (partly) lifted causing the well-known coherent phonon motion of the lattice. While the overall observed dynamics generally follow a displacive excitation model, a delayed onset of the pump-induced carrier dynamics due to hot-carrier thermalization is observed, as well as a large spectral phase dependence in the coherent phonon oscillation. This is attributed to significantly different carrier relaxation timescales for carrier energies above and near the Fermi level of the semimetal and corroborated by a simple theoretical model that considers the carrier relaxation timescales in the displacive phonon model to explain the observed dynamics. Our results provide direct experimental evidence about the role of carrier-relaxation in the origin of displacive coherent phonon motion., Comment: Revised manuscript version, 19 pages, 13 figures, 1 table

A complicating factor in the realization and observation of quantum spin liquids in materials is the ubiquitous presence of other degrees of freedom, in particular lattice distortion modes (phonons). These provide additional routes for relieving magnetic frustration, thereby possibly destabilizing spin-liquid ground states. In this work, we focus on triangular-lattice Heisenberg antiferromagnets, where recent numerical evidence suggests the presence of an extended U(1) Dirac spin liquid phase which is described by compact quantum electrodynamics in 2+1 dimensions (QED$_3$), featuring gapless spinons and monopoles as gauge excitations, and believed to flow to a strongly-coupled fixed point with conformal symmetry. Using complementary perturbation theory and scaling arguments, we show that a symmetry-allowed coupling between (classical) finite-wavevector lattice distortions and monopole operators of the U(1) Dirac spin liquid generally induces a spin-Peierls instability towards a (confining) 12-site valence-bond solid state. We support our theoretical analysis with state-of-the-art density matrix renormalization group simulations. Away from the limit of static distortions, we demonstrate that the phonon energy gap establishes a parameter regime where the spin liquid is expected to be stable, and show that the monopole-lattice coupling leads to softening of the phonon in analogy to the Kohn anomaly. We discuss the applicability of our results to similar systems, in particular the Dirac spin liquid on the Kagome lattice., Comment: 26 pages, 14 figures

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

We describe a simple parallel-friendly lightweight graph reordering algorithm for COO graphs (edge lists). Our ``Batched Order By Attachment'' (BOBA) algorithm is linear in the number of edges in terms of reads and linear in the number of vertices for writes through to main memory. It is highly parallelizable on GPUs\@. We show that, compared to a randomized baseline, the ordering produced gives improved locality of reference in sparse matrix-vector multiplication (SpMV) as well as other graph algorithms. Moreover, it can substantially speed up the conversion from a COO representation to the compressed format CSR, a very common workflow. Thus, it can give \emph{end-to-end} speedups even in SpMV\@. Unlike other lightweight approaches, this reordering does not rely on explicitly knowing the degrees of the vertices, and indeed its runtime is comparable to that of computing degrees. Instead, it uses the structure and edge distribution inherent in the input edge list, making it a candidate for default use in a pragmatic graph creation pipeline. This algorithm is suitable for road-type networks as well as scale-free. It improves cache locality on both CPUs and GPUs, achieving hit rates similar to the heavyweight techniques (e.g., for SpMV, 7--52\% and 11--67\% in the L1 and L2 caches, respectively). Compared to randomly labeled graphs, BOBA-reordered graphs achieve end-to-end speedups of up to 3.45. The reordering time is approximately one order of magnitude faster than existing lightweight techniques and up to 2.5 orders of magnitude faster than heavyweight techniques.

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

We consider two heavy particles immersed in a Bose-Einstein condensate in three dimensions and compute their mutual interaction induced by excitations of the medium. For an ideal Bose gas, the induced interaction is Newtonian up to a shift in distance which depends on the coupling strength between impurities and Bosons. For a real BEC, we find that on short distances, the induced potential is dominated by three-body physics of a single Boson bound to the impurities, leading to an Efimov potential. At large distances of the order of the healing length, a Yukawa potential emerges instead. In particular, we find that both regimes are realized for all impurity-boson couplings and determine the corresponding crossover scales. The transition from the real to the ideal condensate at low gas parameters is investigated., Comment: 7 pages, 5 figures