Search

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, submitted to Astronomy & Astrophysics

Precise mass constraints are vital for the characterisation of brown dwarfs and exoplanets. Here, we present how the combination of data obtained by Gaia and GRAVITY can help enlarge the sample of substellar companions with measured dynamical masses. We show how the Non-Single-Star (NSS) two-body orbit catalogue contained in Gaia DR3 can be used to inform high angular resolution follow-up observations with GRAVITY. Applying the method presented in this work for eight Gaia candidate systems, we detect all eight predicted companions, seven of which being previously unknown and five of substellar nature. Among the sample is Gaia DR3 2728129004119806464 B, which - detected at a angular separation of (34.01 $\pm$ 0.15) mas from the host - is the closest substellar companion ever imaged. In combination with the system's distance this translates to a physical separation of (1.054 $\pm$ 0.002) AU. The GRAVITY data are 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 enables us to further characterise them in terms of their age, effective temperature and radius by the application of evolutionary models. 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 hold great promise for efficiently increasing the sample of directly-imaged companions in the future, especially in the light of the Gaia's upcoming DR4 and the advent of GRAVITY+., Comment: submitted to Astronomy and Astrophysics

The ALMA [CII] Resolved Ism in STar-forming gALaxies (CRISTAL) survey is a Cycle 8 ALMA Large Programme that studies the cold gas component of high-redshift galaxies. Its sub-arcsecond resolution observations are key to disentangling physical mechanisms that shape galaxies during cosmic dawn. In this paper, we explore the morphology and kinematics of the cold gas, star-forming, and stellar components in the star-forming main-sequence galaxy CRISTAL-05/HZ3, at z = 5.54. Our analysis includes 0.3" spatial resolution (~2 kpc) ALMA observations of the [CII] line. While CRISTAL-05 was previously classified as a single source, our observations reveal that the system is a close interacting pair surrounded by an extended component of carbon-enriched gas. This is imprinted in the disturbed elongated [CII] morphology and the separation of the two components in the position-velocity diagram (~100 km/s). The central region is composed of two components, named C05-NW and C05-SE, with the former being the dominant one. A significant fraction of the [CII] arises beyond the close pair up to 10 kpc, while the regions forming new massive stars and the stellar component seem compact (r_[CII] ~ 4 r_UV), as traced by rest-frame UV and optical imaging obtained with the Hubble Space Telescope and the James Webb Space Telescope. Our kinematic model, using the DYSMALpy software, yields a minor contribution of dark matter of C05-NW within a radius of ~2x Reff. Finally, we explore the resolved [CII]/FIR ratios as a proxy for shock-heating produced by this merger. We argue that the extended [CII] emission is mainly caused by the merger, which could not be discerned with lower-resolution observations. Our work emphasizes the need for high-resolution observations to fully characterize the dynamic stages of infant galaxies and the physical mechanisms that drive the metal enrichment of the circumgalactic medium., Comment: Submitted to A&A - comments are welcome! - 19 pages, 13 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

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

The broad-line region (BLR) of active galactic nuclei (AGNs) traces gas close to the central supermassive black hole (BH). Recent reverberation mapping (RM) and interferometric spectro-astrometry data have enabled detailed investigations of the BLR structure and dynamics, as well as estimates of the BH mass. These exciting developments motivate comparative investigations of BLR structures using different broad emission lines. In this work, we have developed a method to simultaneously model multiple broad lines of the BLR from a single-epoch spectrum. We apply this method to the five strongest broad emission lines (H$\alpha$, H$\beta$, H$\gamma$, Pa$\beta$, and He $I\;\lambda$5876) in the UV-to-NIR spectrum of NGC 3783, a nearby Type I AGN which has been well studied by RM and interferometric observations. Fixing the BH mass to the published value, we fit these line profiles simultaneously to constrain the BLR structure. We find that the differences between line profiles can be explained almost entirely as being due to different radial distributions of the line emission. We find that using multiple lines in this way also enables one to measure some important physical parameters, such as the inclination angle and virial factor of the BLR. The ratios of the derived BLR time lags are consistent with the expectation of theoretical model calculations and RM measurements., Comment: 16 pages, 9 figures, 4 tables, accepted for publication in A&A

Massive star-forming galaxies in the high-redshift universe host large reservoirs of cold gas in their circumgalactic medium (CGM). Traditionally, these reservoirs have been linked to diffuse H I Lyman-$\alpha$ (Ly$\alpha)$ emission extending beyond $\approx 10$ kpc scales. In recent years, millimeter/submillimeter observations are starting to identify even colder gas in the CGM through molecular and/or atomic tracers such as the [C II] $158\,\mu$m transition. In this context, we study the well-known J1000+0234 system at $z=4.54$ that hosts a massive dusty star-forming galaxy (DSFG), a UV-bright companion, and a Ly$\alpha$ blob. We combine new ALMA [C II] line observations taken by the CRISTAL survey with data from previous programs targeting the J1000+0234 system, and achieve a deep view into a DSFG and its rich environment at a 0.2" resolution. We identify an elongated [C II]-emitting structure with a projected size of 15 kpc stemming from the bright DSFG at the center of the field, with no clear counterpart at any other wavelength. The plume is oriented $\approx 40^{\circ}$ away from the minor axis of the DSFG, and shows significant spatial variation of its spectral parameters. In particular, the [C II] emission shifts from 180 km/s to 400 km/s between the bottom and top of the plume, relative to the DSFG's systemic velocity. At the same time, the line width starts at 400-600 km/s but narrows down to 190 km/s at top end of the plume. We discuss four possible scenarios to interpret the [C II] plume: a conical outflow, a cold accretion stream, ram pressure stripping, and gravitational interactions. While we cannot strongly rule out any of these with the available data, we disfavor the ram pressure stripping scenario due to the requirement of special hydrodynamic conditions., Comment: 17 pages (14 main text, 2 for references and 1 appendix page), 7 figures and 4 tables. Submitted to A&A

By using the GRAVITY instrument with the near-infrared (NIR) Very Large Telescope Interferometer (VLTI), the structure of the broad (emission-)line region (BLR) in active galactic nuclei (AGNs) can be spatially resolved, allowing the central black hole (BH) mass to be determined. This work reports new NIR VLTI/GRAVITY interferometric spectra for four type 1 AGNs (Mrk 509, PDS 456, Mrk 1239, and IC 4329A) with resolved broad-line emission. Dynamical modelling of interferometric data constrains the BLR radius and central BH mass measurements for our targets and reveals outflow-dominated BLRs for Mrk 509 and PDS 456. We present an updated radius-luminosity (R-L) relation independent of that derived with reverberation mapping (RM) measurements using all the GRAVITY-observed AGNs. We find our R-L relation to be largely consistent with that derived from RM measurements except at high luminosity, where BLR radii seem to be smaller than predicted. This is consistent with RM-based claims that high Eddington ratio AGNs show consistently smaller BLR sizes. The BH masses of our targets are also consistent with the standard $M_\mathrm{BH}$-$\sigma_*$ relation. Model-independent photocentre fitting shows spatial offsets between the hot dust continuum and the BLR photocentres (ranging from $\sim$17 $\mu$as to 140 $\mu$as) that are generally perpendicular to the alignment of the red- and blueshifted BLR photocentres. These offsets are found to be related to the AGN luminosity and could be caused by asymmetric K-band emission of the hot dust, shifting the dust photocentre. We discuss various possible scenarios that can explain this phenomenon., Comment: Accepted for publication at A&A

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 report high-quality H${\alpha}$/CO, imaging spectroscopy of nine massive (log median stellar mass = 10.65 $M_{\odot}$), disk galaxies on the star-forming, main sequence (henceforth `SFGs'), near the peak of cosmic galaxy evolution ($z\sim$1.1-2.5), taken with the ESO-Very Large Telescope, IRAM-NOEMA and Atacama Large Millimeter/submillimeter Array. We fit the major axis position-velocity cuts with beam-convolved, forward models with a bulge, a turbulent rotating disk, and a dark matter (DM) halo. We include priors for stellar and molecular gas masses, optical light effective radii and inclinations, and DM masses from our previous rotation curve analyses of these galaxies. We then subtract the inferred 2D model-galaxy velocity and velocity dispersion maps from those of the observed galaxies. We investigate whether the residual velocity and velocity dispersion maps show indications for radial flows. We also carry out kinemetry, a model-independent tool for detecting radial flows. We find that all nine galaxies exhibit significant non-tangential flows. In six SFGs, the inflow velocities ($v_r\sim$30-90 km s$^{-1}$, 10%-30% of the rotational component) are along the minor axis of these galaxies. In two cases the inflow appears to be off the minor axis. The magnitudes of the radial motions are in broad agreement with the expectations from analytic models of gravitationally unstable, gas-rich disks. Gravitational torques due to clump and bar formation, or spiral arms, drive gas rapidly inward and result in the formation of central disks and large bulges. If this interpretation is correct, our observations imply that gas is transported into the central regions on ~10 dynamical time scales., Comment: 31 pages, 27 figures, 3 tables. The Astrophysical Journal in press

ERIS, the Enhanced Resolution Imager and Spectrograph, is an instrument that both extends and enhances the fundamental diffraction limited imaging and spectroscopy capability for the VLT. It replaces two instruments that were being maintained beyond their operational lifetimes, combines their functionality on a single focus, provides a new wavefront sensing module for natural and laser guide stars that makes use of the Adaptive Optics Facility, and considerably improves on their performance. The observational modes ERIS provides are integral field spectroscopy at 1-2.5 {\mu}m, imaging at 1-5 {\mu}m with several options for high contrast imaging, and longslit spectroscopy at 3-4 {\mu}m, The instrument is installed at the Cassegrain focus of UT4 at the VLT and, following its commissioning during 2022, has been made available to the community., Comment: 20 pages with 29 figures; accepted for A&A (minor changes)

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

The 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

We present a kinematic analysis of the main-sequence galaxy HZ4 at $z=5.5$. Our study is based on deep, spatially resolved observations of the [CII] 158 $\mu$m transition obtained with the Atacama Large Millimeter/Submillimeter Array (ALMA). From the combined analysis of the disk morphology, the two-dimensional velocity structure, and forward-modeling of the one-dimensional velocity and velocity dispersion profiles, we conclude that HZ4 has a regular rotating disk in place. The intrinsic velocity dispersion in HZ4 is high ($\sigma_{0}=65.8^{+2.9}_{-3.3}$ km s$^{-1}$), and the ratio between the rotational velocity and the intrinsic velocity dispersion is $V_{\rm rot}/\sigma_{0}=2.2$. These values are consistent with the expectations from the trends of increasing $\sigma_{0}$ and decreasing $V_{\rm rot}/\sigma_{0}$ as a function of redshift observed in main-sequence galaxies up to $z\approx4$. Galaxy evolution models suggest that the high level of turbulence observed in HZ4 can only be achieved if, in addition to stellar feedback, there is radial transport of gas within the disk. Finally, we find that HZ4 is baryon dominated on galactic scales ($\lesssim2\times R_{\rm e}$), with a dark matter fraction at one effective radius of $f_{\rm DM}(R_{\rm e})=0.41^{+0.25}_{-0.22}$. This value is comparable to the dark matter fractions found in lower redshift galaxies that could be the descendants of HZ4: massive ($M_{\star}\approx10^{11}~M_{\odot}$), star-forming galaxies at $z\sim2$, and passive, early type galaxies at $z\approx0$., Comment: Accepted for publication in A&A Letters, 10 pages, 6 figures

We present a high-resolution kinematic study of the massive main-sequence star-forming galaxy (SFG) SDSS J090122.37+181432.3 (J0901) at z=2.259, using 0.36 arcsec ALMA CO(3-2) and 0.1-0.5 arcsec SINFONI/VLT H-alpha observations. J0901 is a rare, strongly-lensed but otherwise normal massive (log(M_star/M_sun)~11) main sequence SFG, offering a unique opportunity to study a typical massive SFG under the microscope of lensing. Through forward dynamical modeling incorporating lensing deflection, we fit the CO and H-alpha kinematics in the image plane out to about one disk effective radius (R_e ~ 4 kpc) at a ~600pc delensed physical resolution along the kinematic major axis. Our results show high intrinsic dispersions of the cold molecular and warm ionized gas (sig0_mol ~ 40 km/s and sig0_ion ~ 66 km/s) that remain constant out to R_e; a moderately low dark matter fraction (f_DM(R_e) ~ 0.3-0.4) within R_e; and a centrally-peaked Toomre Q-parameter -- agreeing well with the previously established sig0 vs. z, f_DM vs. Sig_baryon, and Q's radial trends using large-sample non-lensed main sequence SFGs. Our data further reveal a high stellar mass concentration within ~1-2 kpc with little molecular gas, and a clumpy molecular gas ring-like structure at R ~ 2-4 kpc, in line with the inside-out quenching scenario. Our further analysis indicates that J0901 had assembled half of its stellar mass only ~400 Myrs before its observed cosmic time, and cold gas ring and dense central stellar component are consistent with signposts of a recent wet compaction event of a highly turbulent disk found in recent simulations., Comment: 29 pages and 21 figures in total (14 pages and 10 figures in main text and the rest in appendix). Accepted for publication in ApJ

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

We analyze Ha or CO rotation curves (RCs) extending out to several galaxy effective radii for 100 massive, large, star-forming disk galaxies (SFGs) across the peak of cosmic galaxy star formation (z~0.6-2.5), more than doubling the previous sample presented by Genzel et al. (2020) and Price et al. (2021). The observations were taken with SINFONI and KMOS integral-field spectrographs at ESO-VLT, LUCI at LBT, NOEMA at IRAM, and ALMA. We fit the major axis kinematics with beam-convolved, forward models of turbulent rotating disks with bulges embedded in dark matter (DM) halos, including the effects of pressure support. The fraction of dark to total matter within the disk effective radius ($R_e ~ 5 kpc$), $f_DM (R_e)=V_{DM}^2 (R_e)/V_{circ}^2 (R_e)$, decreases with redshift: At z~1 (z~2) the median DM fraction is $0.38\pm 0.23$ ($0.27\pm 0.18$), and a third (half) of all galaxies are "maximal" disks with $f_{DM} (R_e)<0.28$. Dark matter fractions correlate inversely with the baryonic surface density, and the low DM fractions require a flattened, or cored, inner DM density distribution. At z~2 there is ~40% less dark matter mass on average within $R_e$ compared to expected values based on cosmological stellar-mass halo-mass relations. The DM deficit is more evident at high star formation rate (SFR) surface densities ($\Sigma_{SFR}>2.5 M_{\odot} yr^{-1} kpc^{-2}$) and galaxies with massive bulges ($M_{bulge}>10^{10} M_{\odot}$). A combination of stellar or active galactic nucleus (AGN) feedback, and/or heating due to dynamical friction, either from satellite accretion or clump migration, may drive the DM from cuspy into cored mass distributions. The observations plausibly indicate an efficient build-up of massive bulges and central black holes at z~2 SFGs., Comment: Submitted to ApJ (34 pages, 12 figures, 4 tables)

Using kinematics to decompose galaxies' mass profiles, including the dark matter contribution, often requires parameterization of the baryonic mass distribution based on ancillary information. One such model choice is a deprojected S\'ersic profile with an assumed intrinsic geometry. The case of flattened, deprojected S\'ersic models has previously been applied to flattened bulges in local star-forming galaxies (SFGs), but can also be used to describe the thick, turbulent disks in distant SFGs. Here we extend this previous work that derived density ($\rho$) and circular velocity ($v_{\rm circ}$) curves by additionally calculating the spherically-enclosed 3D mass profiles ($M_{\rm sph}$). Using these profiles, we compare the projected and 3D mass distributions, quantify the differences between the projected and 3D half-mass radii ($R_{\rm e}; r_{\rm 1/2,mass,3D}$), and present virial coefficients relating $v_{\rm circ}(R)$ and $M_{\rm sph}(

GN-108036 is a star-forming galaxy at $z=7.21$, and one of the most distant known sources in the Northern hemisphere. Based on observations from the NOrthern Extended Millimeter Array (NOEMA), here we report the tentative detection of the [CII] line at $\approx4\sigma$ significance. The integrated [CII] line emission is spatially offset about $\sim4$ kpc from the rest-frame ultraviolet (UV) emission. The total [CII] luminosity ($L_{\rm [CII]}=2.7\times10^8~L_{\odot}$) is consistent with the relation between [CII] luminosity and star formation rate (SFR) observed in nearby and high-$z$ star forming galaxies. More interestingly, the [CII] line is blueshifted with respect to the Ly$\alpha$ line by $980\pm10$ km s$^{-1}$. If confirmed, this corresponds to the largest velocity offset reported to date between the Ly$\alpha$ line and a non-resonant line at $z\gtrsim6$. According to trends observed in other high redshift galaxies, the large Ly$\alpha$ velocity offset in GN-108036 is consistent with its low Ly$\alpha$ equivalent width and high UV absolute magnitude. Based on Ly$\alpha$ radiative transfer models of expanding shells, the large Ly$\alpha$ velocity offset in GN-108036 could be interpreted as the presence of a large column density of hydrogen gas, and/or an outflow with a velocity of $v_{\rm out}\sim\Delta v_{\rm Ly \alpha}/2\sim500$ km s$^{-1}$. We also report the 3$\sigma$ detection of a potential galaxy companion located $\sim30$ kpc east of GN-108036, at a similar systemic velocity, and with no counterpart rest-frame UV emission., Comment: Accepted for publication in Astronomy & Astrophysics Letters, 8 pages, 6 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

We present a kinematic analysis of the main-sequence galaxy HZ4 at $z=5.5$. Our study is based on deep, spatially resolved observations of the [CII] 158 $\mu$m transition obtained with the Atacama Large Millimeter/Submillimeter Array (ALMA). From the combined analysis of the disk morphology, the two-dimensional velocity structure, and forward-modeling of the one-dimensional velocity and velocity dispersion profiles, we conclude that HZ4 has a regular rotating disk in place. The intrinsic velocity dispersion in HZ4 is high ($\sigma_{0}=65.8^{+2.9}_{-3.3}$ km s$^{-1}$), and the ratio between the rotational velocity and the intrinsic velocity dispersion is $V_{\rm rot}/\sigma_{0}=2.2$. These values are consistent with the expectations from the trends of increasing $\sigma_{0}$ and decreasing $V_{\rm rot}/\sigma_{0}$ as a function of redshift observed in main-sequence galaxies up to $z\approx4$. Galaxy evolution models suggest that the high level of turbulence observed in HZ4 can only be achieved if, in addition to stellar feedback, there is radial transport of gas within the disk. Finally, we find that HZ4 is baryon dominated on galactic scales ($\lesssim2\times R_{\rm e}$), with a dark matter fraction at one effective radius of $f_{\rm DM}(R_{\rm e})=0.41^{+0.25}_{-0.22}$. This value is comparable to the dark matter fractions found in lower redshift galaxies that could be the descendants of HZ4: massive ($M_{\star}\approx10^{11}~M_{\odot}$), star-forming galaxies at $z\sim2$, and passive, early type galaxies at $z\approx0$., Comment: Accepted for publication in A&A Letters, 10 pages, 6 figures

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

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

We present a follow-up analysis examining the dynamics and structures of 41 massive, large star-forming galaxies at z~0.67-2.45 using both ionized and molecular gas kinematics. We fit the galaxy dynamics with models consisting of a bulge, a thick, turbulent disk, and a NFW dark matter halo, using code that fully forward models the kinematics, including all observational and instrumental effects. We explore the parameter space using Markov Chain Monte Carlo (MCMC) sampling, including priors based on stellar and gas masses and disk sizes. We fit the full sample using extracted 1D kinematic profiles. For a subset of 14 well-resolved galaxies, we also fit the 2D kinematics. The MCMC approach robustly confirms the results from least-squares fitting presented in Paper I (Genzel et al. 2020): the sample galaxies tend to be baryon-rich on galactic scales (within one effective radius). The 1D and 2D MCMC results are also in good agreement for the subset, demonstrating that much of the galaxy dynamical information is captured along the major axis. The 2D kinematics are more affected by the presence of non-circular motions, which we illustrate by constructing a toy model with constant inflow for one galaxy that exhibits residual signatures consistent with radial motions. This analysis, together with results from Paper I and other studies, strengthens the finding that massive, star-forming galaxies at z~1-2 are baryon-dominated on galactic scales, with lower dark matter fractions towards higher baryonic surface densities. Finally, we present details of the kinematic fitting code used in this analysis., Comment: Accepted for publication in ApJ (23 pages, 8 figures, 5 tables)

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

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

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

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

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

We present a follow-up analysis examining the dynamics and structures of 41 massive, large star-forming galaxies at z~0.67-2.45 using both ionized and molecular gas kinematics. We fit the galaxy dynamics with models consisting of a bulge, a thick, turbulent disk, and a NFW dark matter halo, using code that fully forward models the kinematics, including all observational and instrumental effects. We explore the parameter space using Markov Chain Monte Carlo (MCMC) sampling, including priors based on stellar and gas masses and disk sizes. We fit the full sample using extracted 1D kinematic profiles. For a subset of 14 well-resolved galaxies, we also fit the 2D kinematics. The MCMC approach robustly confirms the results from least-squares fitting presented in Paper I (Genzel et al. 2020): the sample galaxies tend to be baryon-rich on galactic scales (within one effective radius). The 1D and 2D MCMC results are also in good agreement for the subset, demonstrating that much of the galaxy dynamical information is captured along the major axis. The 2D kinematics are more affected by the presence of non-circular motions, which we illustrate by constructing a toy model with constant inflow for one galaxy that exhibits residual signatures consistent with radial motions. This analysis, together with results from Paper I and other studies, strengthens the finding that massive, star-forming galaxies at z~1-2 are baryon-dominated on galactic scales, with lower dark matter fractions towards higher baryonic surface densities. Finally, we present details of the kinematic fitting code used in this analysis., Comment: Accepted for publication in ApJ (23 pages, 8 figures, 5 tables)

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

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

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

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

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

We report high quality, Halpha or CO rotation curves (RCs) to several Re for 41 large, massive, star-forming disk galaxies (SFGs), across the peak of cosmic galaxy evolution (z~0.67-2.45), taken with the ESO-VLT, the LBT and IRAM-NOEMA. Most RC41 SFGs have reflection symmetric RCs plausibly described by equilibrium dynamics. We fit the major axis position-velocity cuts with beam-convolved, forward modeling with a bulge, a turbulent rotating disk, and a dark matter (DM) halo. We include priors for stellar and molecular gas masses, optical light effective radii and inclinations, and DM masses from abundance matching scaling relations. Two-thirds or more of the z>1.2 SFGs are baryon dominated within a few Re of typically 5.5 kpc, and have DM fractions less than maximal disks (=0.12). At lower redshift (z<1.2) that fraction is less than one-third. DM fractions correlate inversely with the baryonic angular momentum parameter, baryonic surface density and bulge mass. Inferred low DM fractions cannot apply to the entire disk & halo but more plausibly reflect a flattened, or cored, inner DM density distribution. The typical central 'DM deficit' in these cores relative to NFW distributions is ~30% of the bulge mass. The observations are consistent with rapid radial transport of baryons in the first generation massive gas rich halos forming globally gravitationally unstable disks, and leading to efficient build-up of massive bulges and central black holes. A combination of heating due to dynamical friction and AGN feedback may drive DM out of the initial cusps., Comment: 15 figures, 4 tables, accepted for publication in ApJ

We use deep adaptive optics assisted integral field spectroscopy from SINFONI on the VLT to study the spatially resolved properties of ionized gas outflows driven by active galactic nuclei (AGN) in three galaxies at z~2.2 -- K20-ID5, COS4-11337 and J0901+1814. These systems probe AGN feedback from nuclear to circumgalactic scales, and provide unique insights into the different mechanisms by which AGN-driven outflows interact with their host galaxies. K20-ID5 and COS4-11337 are compact star forming galaxies with powerful $\sim$1500 km s$^{-1}$ AGN-driven outflows that dominate their nuclear H$\alpha$ emission. The outflows do not appear to have any impact on the instantaneous star formation activity of the host galaxies, but they carry a significant amount of kinetic energy which could heat the halo gas and potentially lead to a reduction in the rate of cold gas accretion onto the galaxies. The outflow from COS4-11337 is propagating directly towards its companion galaxy COS4-11363, at a projected separation of 5.4 kpc. COS4-11363 shows signs of shock excitation and recent truncation of star formation activity, which could plausibly have been induced by the outflow from COS4-11337. J0901+1814 is gravitationally lensed, giving us a unique view of a compact (R = 470 $\pm$ 70 pc), relatively low velocity ($\sim$650 km s$^{-1}$) AGN-driven outflow. J0901+1814 has a similar AGN luminosity to COS4-11337, suggesting that the difference in outflow properties is not related to the current AGN luminosity, and may instead reflect a difference in the evolutionary stage of the outflow and/or the coupling efficiency between the AGN ionizing radiation field and the gas in the nuclear regions., Comment: Accepted for publication in ApJ. Main text 23 pages, 15 figures and 4 tables, plus Appendix (3 pages, 3 figures, 1 table)

The Galactic Center black hole Sagittarius A* is a variable NIR source that exhibits bright flux excursions called flares. The low-flux density turnover of the flux distribution is below the sensitivity of current single-aperture telescopes. We use the unprecedented resolution of the GRAVITY instrument at the VLTI. Our light curves are unconfused, overcoming the confusion limit of previous photometric studies. We analyze the light curves using standard statistical methods and obtain the flux distribution. We find that the flux distribution of SgrA* turns over at a median flux density of (1.1\pm0.3)mJy. We measure the percentiles of the flux distribution and use them to constrain the NIR K-band SED. Furthermore, we find that the flux distribution is intrinsically right-skewed to higher flux density in log space. Flux densities below 0.1mJy are hardly ever observed. In consequence, a single powerlaw or lognormal distribution does not suffice to describe the observed flux distribution in its entirety. However, if one takes into account a power law component at high flux densities, a lognormal distribution can describe the lower end of the observed flux distribution. We confirm the RMS-flux relation for Sgr~A* and find it to be linear for all flux densities in our observation. We conclude that Sgr~A* has two states: the bulk of the emission is generated in a lognormal process with a well-defined median flux density and this quiescent emission is supplemented by sporadic flares that create the observed power law extension of the flux distribution.

Infrared observations of Sgr A* probe the region close to the event horizon of the black hole at the Galactic center. These observations can constrain the properties of low-luminosity accretion as well as that of the black hole itself. The GRAVITY instrument at the ESO VLTI has recently detected continuous circular relativistic motion during infrared flares which has been interpreted as orbital motion near the event horizon. Here we analyze the astrometric data from these flares, taking into account the effects of out-of-plane motion and orbital shear of material near the event horizon of the black hole. We have developed a new code to predict astrometric motion and flux variability from compact emission regions following particle orbits. Our code combines semi-analytic calculations of timelike geodesics that allow for out-of-plane or elliptical motions with ray tracing of photon trajectories to compute time-dependent images and light curves. We apply our code to the three flares observed with GRAVITY in 2018. We show that all flares are consistent with a hotspot orbiting at R$\sim$9 gravitational radii with an inclination of $i\sim140^\circ$. The emitting region must be compact and less than $\sim5$ gravitational radii in diameter. We place a further limit on the out-of-plane motion during the flare., Comment: Accepted in A&A

We investigate what drives the redshift evolution of the typical electron density ($n_e$) in star-forming galaxies, using a sample of 140 galaxies drawn primarily from KMOS$^{\rm 3D}$ ($0.6\lesssim{z}\lesssim{2.6}$) and 471 galaxies from SAMI ($z<0.113$). We select galaxies that do not show evidence of AGN activity or outflows, to constrain the average conditions within H II regions. Measurements of the [SII]$\lambda$6716/[SII]$\lambda$6731 ratio in four redshift bins indicate that the local $n_e$ in the line-emitting material decreases from 187$^{+140}_{-132}$ cm$^{-3}$ at $z\sim$ 2.2 to 32$^{+4}_{-9}$ cm$^{-3}$ at $z\sim$ 0; consistent with previous results. We use the H$\alpha$ luminosity to estimate the root-mean-square (rms) $n_e$ averaged over the volumes of star-forming disks at each redshift. The local and volume-averaged $n_e$ evolve at similar rates, hinting that the volume filling factor of the line-emitting gas may be approximately constant across $0\lesssim{z}\lesssim{2.6}$. The KMOS$^{\rm 3D}$ and SAMI galaxies follow a roughly monotonic trend between $n_e$ and star formation rate, but the KMOS$^{\rm 3D}$ galaxies have systematically higher $n_e$ than the SAMI galaxies at fixed offset from the star-forming main sequence, suggesting a link between the $n_e$ evolution and the evolving main sequence normalization. We quantitatively test potential drivers of the density evolution and find that $n_e$(rms) $\simeq{n_{H_2}}$, suggesting that the elevated $n_e$ in high-$z$ H II regions could plausibly be the direct result of higher densities in the parent molecular clouds. There is also tentative evidence that $n_e$ could be influenced by the balance between stellar feedback, which drives the expansion of H II regions, and the ambient pressure, which resists their expansion., Comment: Main text 24 pages, 8 figures. Accepted for publication in ApJ

The spin of the supermassive black hole that resides at the Galactic Centre can in principle be measured by accurate measurements of the orbits of stars that are much closer to SgrA* than S2, the orbit of which recently provided the measurement of the gravitational redshift and the Schwarzschild precession. The GRAVITY near-infrared interferometric instrument combining the four 8m telescopes of the VLT provides a spatial resolution of 2-4 mas, breaking the confusion barrier for adaptive-optics-assisted imaging with a single 8-10m telescope. We used GRAVITY to observe SgrA* over a period of six months in 2019 and employed interferometric reconstruction methods developed in radio astronomy to search for faint objects near SgrA*. This revealed a slowly moving star of magnitude 18.9 in K band within 30mas of SgrA*. The position and proper motion of the star are consistent with the previously known star S62, which is at a substantially larger physical distance, but in projection passes close to SgrA*. Observations in August and September 2019 easily detected S29, with K-magnitude of 16.6, at approximately 130 mas from SgrA*. The planned upgrades of GRAVITY, and further improvements in the calibration, hold the promise of finding stars fainter than magnitude 19 at K., Comment: 12 pages, 11 figures, published in A&A. Here we updated Fig. 1 with the latest motions laws of the S-stars as obtained from VLT/NACO observations until end 2019

Ever deeper and wider lookback surveys have led to a fairly robust outline of the cosmic star formation history, which culminated around z~2 -- a period often nicknamed "cosmic noon." Our knowledge about star-forming galaxies at these epochs has dramatically advanced from increasingly complete population censuses and detailed views of individual galaxies. We highlight some of the key observational insights that influenced our current understanding of galaxy evolution in the equilibrium growth picture: $\bullet$ scaling relations between galaxy properties are fairly well established among massive galaxies at least out to z~2, pointing to regulating mechanisms already acting on galaxy growth; $\bullet$ resolved views reveal that gravitational instabilities and efficient secular processes within the gas- and baryon-rich galaxies at z~2 play an important role in the early build-up of galactic structure; $\bullet$ ever more sensitive observations of kinematics at z~2 are probing the baryon and dark matter budget on galactic scales and the links between star-forming galaxies and their likely descendants; $\bullet$ towards higher masses, massive bulges, dense cores, and powerful AGN and AGN-driven outflows are more prevalent and likely play a role in quenching star formation. We outline emerging questions and exciting prospects for the next decade with upcoming instrumentation, including the James Webb Space Telescope and the next generation of Extremely Large Telescopes., Comment: 77 pages, 11 figures. Accepted by and posted with permission from Annual Review of Astronomy and Astrophysics, Vol. 58

We present new near-infrared VLTI/GRAVITY interferometric spectra that spatially resolve the broad Br$\gamma$ emission line in the nucleus of the active galaxy IRAS 09149-6206. We use these data to measure the size of the broad line region (BLR) and estimate the mass of the central black hole. Using an improved phase calibration method that reduces the differential phase uncertainty to 0.05 degree per baseline across the spectrum, we detect a differential phase signal that reaches a maximum of ~0.5 degree between the line and continuum. This represents an offset of ~120 $\mu$as (0.14 pc) between the BLR and the centroid of the hot dust distribution traced by the 2.3 $\mu$m continuum. The offset is well within the dust sublimation region, which matches the measured ~0.6 mas (0.7 pc) diameter of the continuum. A clear velocity gradient, almost perpendicular to the offset, is traced by the reconstructed photocentres of the spectral channels of the Br$\gamma$ line. We infer the radius of the BLR to be ~65 $\mu$as (0.075 pc), which is consistent with the radius-luminosity relation of nearby active galactic nuclei derived based on the time lag of the H$\beta$ line from reverberation mapping campaigns. Our dynamical modelling indicates the black hole mass is $\sim 1\times10^8\,M_\odot$, which is a little below, but consistent with, the standard $M_{\rm BH}$-$\sigma_*$ relation., Comment: Accepted for publication in A&A. 14 pages, 11 figures in main text

We study the time-variable linear polarisation of Sgr A* during a bright NIR flare observed with the GRAVITY instrument on July 28, 2018. Motivated by the time evolution of both the observed astrometric and polarimetric signatures, we interpret the data in terms of the polarised emission of a compact region ('hotspot') orbiting a black hole in a fixed, background magnetic field geometry. We calculated a grid of general relativistic ray-tracing models, created mock observations by simulating the instrumental response, and compared predicted polarimetric quantities directly to the measurements. We take into account an improved instrument calibration that now includes the instrument's response as a function of time, and we explore a variety of idealised magnetic field configurations. We find that the linear polarisation angle rotates during the flare, which is consistent with previous results. The hotspot model can explain the observed evolution of the linear polarisation. In order to match the astrometric period of this flare, the near horizon magnetic field is required to have a significant poloidal component, which is associated with strong and dynamically important fields. The observed linear polarisation fraction of $\simeq 30\%$ is smaller than the one predicted by our model ($\simeq 50\%$). The emission is likely beam depolarised, indicating that the flaring emission region resolves the magnetic field structure close to the black hole., Comment: 13 pages, 14 figures, to be published in A&A

We report high quality, Halpha or CO rotation curves (RCs) to several Re for 41 large, massive, star-forming disk galaxies (SFGs), across the peak of cosmic galaxy evolution (z~0.67-2.45), taken with the ESO-VLT, the LBT and IRAM-NOEMA. Most RC41 SFGs have reflection symmetric RCs plausibly described by equilibrium dynamics. We fit the major axis position-velocity cuts with beam-convolved, forward modeling with a bulge, a turbulent rotating disk, and a dark matter (DM) halo. We include priors for stellar and molecular gas masses, optical light effective radii and inclinations, and DM masses from abundance matching scaling relations. Two-thirds or more of the z>1.2 SFGs are baryon dominated within a few Re of typically 5.5 kpc, and have DM fractions less than maximal disks (=0.12). At lower redshift (z<1.2) that fraction is less than one-third. DM fractions correlate inversely with the baryonic angular momentum parameter, baryonic surface density and bulge mass. Inferred low DM fractions cannot apply to the entire disk & halo but more plausibly reflect a flattened, or cored, inner DM density distribution. The typical central 'DM deficit' in these cores relative to NFW distributions is ~30% of the bulge mass. The observations are consistent with rapid radial transport of baryons in the first generation massive gas rich halos forming globally gravitationally unstable disks, and leading to efficient build-up of massive bulges and central black holes. A combination of heating due to dynamical friction and AGN feedback may drive DM out of the initial cusps., Comment: 15 figures, 4 tables, accepted for publication in ApJ

The Galactic Center black hole Sagittarius A* is a variable NIR source that exhibits bright flux excursions called flares. The low-flux density turnover of the flux distribution is below the sensitivity of current single-aperture telescopes. We use the unprecedented resolution of the GRAVITY instrument at the VLTI. Our light curves are unconfused, overcoming the confusion limit of previous photometric studies. We analyze the light curves using standard statistical methods and obtain the flux distribution. We find that the flux distribution of SgrA* turns over at a median flux density of (1.1\pm0.3)mJy. We measure the percentiles of the flux distribution and use them to constrain the NIR K-band SED. Furthermore, we find that the flux distribution is intrinsically right-skewed to higher flux density in log space. Flux densities below 0.1mJy are hardly ever observed. In consequence, a single powerlaw or lognormal distribution does not suffice to describe the observed flux distribution in its entirety. However, if one takes into account a power law component at high flux densities, a lognormal distribution can describe the lower end of the observed flux distribution. We confirm the RMS-flux relation for Sgr~A* and find it to be linear for all flux densities in our observation. We conclude that Sgr~A* has two states: the bulk of the emission is generated in a lognormal process with a well-defined median flux density and this quiescent emission is supplemented by sporadic flares that create the observed power law extension of the flux distribution.

We use deep adaptive optics assisted integral field spectroscopy from SINFONI on the VLT to study the spatially resolved properties of ionized gas outflows driven by active galactic nuclei (AGN) in three galaxies at z~2.2 -- K20-ID5, COS4-11337 and J0901+1814. These systems probe AGN feedback from nuclear to circumgalactic scales, and provide unique insights into the different mechanisms by which AGN-driven outflows interact with their host galaxies. K20-ID5 and COS4-11337 are compact star forming galaxies with powerful $\sim$1500 km s$^{-1}$ AGN-driven outflows that dominate their nuclear H$\alpha$ emission. The outflows do not appear to have any impact on the instantaneous star formation activity of the host galaxies, but they carry a significant amount of kinetic energy which could heat the halo gas and potentially lead to a reduction in the rate of cold gas accretion onto the galaxies. The outflow from COS4-11337 is propagating directly towards its companion galaxy COS4-11363, at a projected separation of 5.4 kpc. COS4-11363 shows signs of shock excitation and recent truncation of star formation activity, which could plausibly have been induced by the outflow from COS4-11337. J0901+1814 is gravitationally lensed, giving us a unique view of a compact (R = 470 $\pm$ 70 pc), relatively low velocity ($\sim$650 km s$^{-1}$) AGN-driven outflow. J0901+1814 has a similar AGN luminosity to COS4-11337, suggesting that the difference in outflow properties is not related to the current AGN luminosity, and may instead reflect a difference in the evolutionary stage of the outflow and/or the coupling efficiency between the AGN ionizing radiation field and the gas in the nuclear regions., Comment: Accepted for publication in ApJ. Main text 23 pages, 15 figures and 4 tables, plus Appendix (3 pages, 3 figures, 1 table)

Infrared observations of Sgr A* probe the region close to the event horizon of the black hole at the Galactic center. These observations can constrain the properties of low-luminosity accretion as well as that of the black hole itself. The GRAVITY instrument at the ESO VLTI has recently detected continuous circular relativistic motion during infrared flares which has been interpreted as orbital motion near the event horizon. Here we analyze the astrometric data from these flares, taking into account the effects of out-of-plane motion and orbital shear of material near the event horizon of the black hole. We have developed a new code to predict astrometric motion and flux variability from compact emission regions following particle orbits. Our code combines semi-analytic calculations of timelike geodesics that allow for out-of-plane or elliptical motions with ray tracing of photon trajectories to compute time-dependent images and light curves. We apply our code to the three flares observed with GRAVITY in 2018. We show that all flares are consistent with a hotspot orbiting at R$\sim$9 gravitational radii with an inclination of $i\sim140^\circ$. The emitting region must be compact and less than $\sim5$ gravitational radii in diameter. We place a further limit on the out-of-plane motion during the flare., Comment: Accepted in A&A