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.

The star S2 orbiting the compact radio source Sgr A* is a precision probe of the gravitational field around the closest massive black hole (candidate). Over the last 2.7 decades we have monitored the star's radial velocity and motion on the sky, mainly with the SINFONI and NACO adaptive optics (AO) instruments on the ESO VLT, and since 2017, with the four-telescope interferometric beam combiner instrument GRAVITY. In this paper we report the first detection of the General Relativity (GR) Schwarzschild Precession (SP) in S2's orbit. Owing to its highly elliptical orbit (e = 0.88), S2's SP is mainly a kink between the pre-and post-pericentre directions of motion ~ +- 1 year around pericentre passage, relative to the corresponding Kepler orbit. The superb 2017-2019 astrometry of GRAVITY defines the pericentre passage and outgoing direction. The incoming direction is anchored by 118 NACO-AO measurements of S2's position in the infrared reference frame, with an additional 75 direct measurements of the S2-Sgr A* separation during bright states ('flares') of Sgr A*. Our 14-parameter model fits for the distance, central mass, the position and motion of the reference frame of the AO astrometry relative to the mass, the six parameters of the orbit, as well as a dimensionless parameter f_SP for the SP (f_SP = 0 for Newton and 1 for GR). From data up to the end of 2019 we robustly detect the SP of S2, del phi = 12' per orbital period. From posterior fitting and MCMC Bayesian analysis with different weighting schemes and bootstrapping we find f_SP = 1.10 +- 0.19. The S2 data are fully consistent with GR. Any extended mass inside S2's orbit cannot exceed ~ 0.1% of the central mass. Any compact third mass inside the central arcsecond must be less than about 1000 M_sun.

The supermassive black hole Sagittarius A* (Sgr A*) is located at the dynamical center of the Milky Way. In a recent study of the X-ray flaring activity from Sgr A* using Chandra, XMM-Newton and Swift data from 1999 to 2015, it has been argued that the bright flaring rate has increased from 2014 Aug. 31 while the faint flaring rate decreased from around 2013 Aug. We tested the persistence of these changes in the flaring rates with new X-ray data of Sgr A* obtained in 2016-2018 (total exposure of 1.4Ms). We detected 9 flares in the Chandra data and 5 flares in the Swift data that we added to the set of 107 previously detected flares. We computed the intrinsic distribution of flare fluxes and durations corrected for the sensitivity bias using a new method that allowed us to take the error on the flare fluxes and durations into account. From this intrinsic distribution, we determined the average flare detection efficiency for each Chandra, XMM-Newton, and Swift observation. After correcting each observational exposure for this efficiency, we applied the Bayesian blocks algorithm on the concatenated flare arrival times. As in the above-mentioned study, we also searched for a flux and fluence threshold that might lead to a change in flaring rate. We improved the previous method by computing the average flare detection efficiencies for each flux and fluence range. The Bayesian block algorithm did not detect any significant change in flaring rate of the 121 flares. However, we detected an increase by a factor of about three in the flaring rate of the most luminous and most energetic flares that have occurred since 2014 Aug. 30. The X-ray activity of Sgr A* has increased for more than four years. Additional studies about the overall near-infrared and radio behavior of Sgr A* are required to draw strong results on the multiwavelength activity of the black hole.

The inner regions of the discs of high-mass young stellar objects (HMYSOs) are still poorly known due to the small angular scales and the high visual extinction involved. We deploy near-infrared (NIR) spectro-interferometry to probe the inner gaseous disc in HMYSOs and investigate the origin and physical characteristics of the CO bandhead emission (2.3-2.4 $\mu$m). We present the first GRAVITY/VLTI observations at high spectral (R=4000) and spatial (mas) resolution of the CO overtone transitions in NGC 2024 IRS2. The continuum emission is resolved in all baselines and is slightly asymmetric, displaying small closure phases ($\leq$8$^{\circ}$). Our best ellipsoid model provides a disc inclination of 34$^{\circ}$$\pm$1$^{\circ}$, a disc major axis position angle of 166$^{\circ}$$\pm$1$^{\circ}$, and a disc diameter of 3.99$\pm$0.09 mas (or 1.69$\pm$0.04 au, at a distance of 423 pc). The small closure phase signals in the continuum are modelled with a skewed rim, originating from a pure inclination effect. For the first time, our observations spatially and spectrally resolve the first four CO bandheads. Changes in visibility, as well as differential and closure phases across the bandheads are detected. Both the size and geometry of the CO-emitting region are determined by fitting a bidimensional Gaussian to the continuum-compensated CO bandhead visibilities. The CO-emitting region has a diameter of 2.74$\pm^{0.08}_{0.07}$ mas (1.16$\pm$0.03 au), and is located in the inner gaseous disc, well within the dusty rim, with inclination and $PA$ matching the dusty disc geometry, which indicates that both dusty and gaseous discs are coplanar. Physical and dynamical gas conditions are inferred by modelling the CO spectrum. Finally, we derive a direct measurement of the stellar mass of $M_*\sim$14.7$^{+2}_{-3.6}$ M$_{\odot}$ by combining our interferometric and spectral modelling results.

Context. The Event Horizon Telescope (EHT) collaboration recently obtained first images of the surroundings of the supermassive compact object M87* at the center of the galaxy M87. Aims. We want to develop a simple analytic disk model for the accretion flow of M87*. Compared to general-relativistic magnetohydrodynamic (GRMHD) models, it has the advantage of being independent of the turbulent character of the flow, and controlled by only few easy-to-interpret, physically meaningful parameters. We want to use this model to predict the image of M87* assuming that it is either a Kerr black hole, or an alternative compact object. Methods. We compute the synchrotron emission from the disk model and propagate the resulting light rays to the far-away observer by means of relativistic ray tracing. Such computations are performed assuming different spacetimes (Kerr, Minkowski, non-rotating ultracompact star, rotating boson star or Lamy spinning wormhole). We perform numerical fits of these models to the EHT data. Results. We discuss the highly-lensed features of Kerr images and show that they are intrinsically linked to the accretion-flow properties, and not only to gravitation. This fact is illustrated by the notion of secondary ring that we introduce. Our model of spinning Kerr black hole predicts mass and orientation consistent with the EHT interpretation. The non-Kerr images result in similar quality of the numerical fits and may appear very similar to Kerr images, once blurred to the EHT resolution. This implies that a strong test of the Kerr spacetime may be out of reach with the current data. We notice that future developments of the EHT could alter this situation. Conclusions. Our results show the importance of studying alternatives to the Kerr spacetime in order to be able to test the Kerr paradigm unambiguously.

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.

Our objective is to estimate the C/O ratio in the atmosphere of beta Pictoris b and obtain an estimate of the dynamical mass of the planet, as well as to refine its orbital parameters using high-precision astrometry. We used the GRAVITY instrument with the four 8.2 m telescopes of the Very Large Telescope Interferometer to obtain K-band spectro-interferometric data on $\beta$ Pic b. We extracted a medium resolution (R=500) K-band spectrum of the planet and a high-precision astrometric position. We estimated the planetary C/O ratio using two different approaches (forward modeling and free retrieval) from two different codes (ExoREM and petitRADTRANS, respectively). Finally, we used a simplified model of two formation scenarios (gravitational collapse and core-accretion) to determine which can best explain the measured C/O ratio. Our new astrometry disfavors a circular orbit for $\beta$ Pic b ($e=0.15^{+0.05}_{-0.04}$). Combined with previous results and with Hipparcos/GAIA measurements, this astrometry points to a planet mass of $M = 12.7\pm{}2.2\,M_\mathrm{Jup}$. This value is compatible with the mass derived with the free-retrieval code petitRADTRANS using spectral data only. The forward modeling and free-retrieval approches yield very similar results regarding the atmosphere of beta Pic b. In particular, the C/O ratios derived with the two codes are identical ($0.43\pm{}0.05$ vs $0.43^{+0.04}_{-0.03}$). We argue that if the stellar C/O in $\beta$ Pic is Solar, then this combination of a very high mass and a low C/O ratio for the planet suggests a formation through core-accretion, with strong planetesimal enrichment.

The formation and the evolution of protoplanetary disks are important stages in the lifetime of stars. The processes of disk evolution and planet formation are intrinsically linked. We spatially resolve with GRAVITY/VLTI in the K-band the sub au-scale region of 27 stars to gain statistical understanding of their properties. We look for correlations with stellar parameters, such as luminosity, mass, temperature and age. Our sample also cover a range of various properties in terms of reprocessed flux, flared or flat morphology, and gaps. We developed semi-physical geometrical models to fit our interferometric data. Our best models correspond to smooth and wide rings, implying that wedge-shaped rims at the dust sublimation edge are favored, as found in the H-band. The closure phases are generally non-null with a median value of ~10 deg, indicating spatial asymmetries of the intensity distributions. Multi-size grain populations could explain the closure phase ranges below 20-25 deg but other scenarios should be invoked to explain the largest ones. Our measurements extend the Radius-Luminosity relation to ~1e4 Lsun and confirm the significant spread around the mean relation observed in the H-band. Gapped sources exhibit a large N-to-K band size ratio and large values of this ratio are only observed for the members of our sample that would be older than 1 Ma, less massive, and with lower luminosity. In the 2 Ms mass range, we observe a correlation in the increase of the relative age with the transition from group II to group I, and an increase of the N-to-K size ratio. However, the size of the current sample does not yet permit us to invoke a clear universal evolution mechanism across the HAeBe mass range. The measured locations of the K-band emission suggest that these disks might be structured by forming young planets, rather than by depletion due to EUV, FUV, and X-ray photo-evaporation.

Precise measurements of the S-stars orbiting SgrA* have set strong constraints on the nature of the compact object at the centre of the Milky Way. The presence of a black hole in that region is well established, but its neighboring environment is still an open debate. In that respect, the existence of dark matter in that central region may be detectable due to its strong signatures on the orbits of stars: the main effect is a Newtonian precession which will affect the overall pericentre shift of S2, the latter being a target measurement of the GRAVITY instrument. The exact nature of this dark matter (e.g., stellar dark remnants or diffuse dark matter) is unknown. This article assumes it to be an scalar field of toroidal distribution, associated with ultra-light dark matter particles, surrounding the Kerr black hole. Such a field is a form of "hair" expected in the context of superradiance, a mechanism that extracts rotational energy from the black hole. Orbital signatures for the S2 star are computed and shown to be detectable by GRAVITY. The scalar field can be constrained because the variation of orbital elements depends both on the relative mass of the scalar field to the black hole and on the field mass coupling parameter.

We present the first fully relativistic study of gravitational radiation from bodies in circular equatorial orbits around the massive black hole at the Galactic Center, Sgr A* and we assess the detectability of various kinds of objects by the gravitational wave detector LISA. Our computations are based on the theory of perturbations of the Kerr spacetime and take into account the Roche limit induced by tidal forces in the Kerr metric. The signal-to-noise ratio in the LISA detector, as well as the time spent in LISA band, are evaluated. We have implemented all the computational tools in an open-source SageMath package, within the Black Hole Perturbation Toolkit framework. We find that white dwarfs, neutrons stars, stellar black holes, primordial black holes of mass larger than $10^{-4} M_\odot$, main-sequence stars of mass lower than $\sim 2.5\, M_\odot$ and brown dwarfs orbiting Sgr A* are all detectable in one year of LISA data with a signal-to-noise ratio above 10 during at least $10^5$ years in the slow inspiral towards either the innermost stable circular orbit (compact objects) or the Roche limit (main-sequence stars and brown dwarfs). The longest times in-band, of the order of $10^6$ years, are achieved for primordial black holes of mass $\sim 10^{-3} M_\odot$ down to $10^{-5} M_\odot$, depending on the spin of Sgr A*, as well as for brown dwarfs, just followed by white dwarfs and low mass main-sequence stars. The long time in-band of these objects makes Sgr A* a valuable target for LISA. We also consider bodies on close circular orbits around the massive black hole in the nucleus of the nearby galaxy M32 and find that, among them, compact objects and brown dwarfs stay for $10^3$ to $10^4$ years in LISA band with a 1-year signal-to-noise ratio above 10.

The rather elusive High-Frequency Quasi-Periodic Oscillations(HFQPO) observed in the X-ray lightcurve of black holes have been seen in a wide range of frequencies, even within one source. It is also notable to have been detected in "pairs" of HFQPOs with a close to integer ratio between the frequencies. The aim of this paper is to investigate some of the possible observable that we could obtain from having the Rossby Wave Instability (RWI) active in the accretion disc surrounding the compact object. Using the newly developed GR-AMRVAC code able to follow the evolution of the RWI in a full general relativistic framework, we explore how RWI can reproduce observed HFQPO frequencies ratios and if it is compatible with the observations. In order to model the emission coming from the disc we have linked our general relativistic simulations to the general relativistic ray-tracing GYOTO code and delivered synthetic observables that can be confronted to actual data from binary systems hosting HFQPOs.}{We have demonstrated in our study that some changes in the physical conditions prevailing in the part of the disc where RWI can be triggered leads to various dominant RWI modes whose ratio recovers frequency ratios observed in various X-ray binary systems. In addition, we have also highlighted that when RWI is triggered near the last stable orbit of a spinning black hole, the amplitude of the X-ray modulation increases with the spin of the black hole. Revisiting published data on X-ray binary systems, we show that this type of relationship actually exists in the five systems where an indirect measure of the spin of the black hole is available.

We present a 0.16% precise and 0.27% accurate determination of R0, the distance to the Galactic Center. Our measurement uses the star S2 on its 16-year orbit around the massive black hole Sgr A* that we followed astrometrically and spectroscopically for 27 years. Since 2017, we added near-infrared interferometry with the VLTI beam combiner GRAVITY, yielding a direct measurement of the separation vector between S2 and Sgr A* with an accuracy as good as 20 micro-arcsec in the best cases. S2 passed the pericenter of its highly eccentric orbit in May 2018, and we followed the passage with dense sampling throughout the year. Together with our spectroscopy, in the best cases with an error of 7 km/s, this yields a geometric distance estimate: R0 = 8178 +- 13(stat.) +- 22(sys.) pc. This work updates our previous publication in which we reported the first detection of the gravitational redshift in the S2 data. The redshift term is now detected with a significance level of 20 sigma with f_redshift = 1.04 +- 0.05.

To date, infrared interferometry at best achieved contrast ratios of a few times $10^{-4}$ on bright targets. GRAVITY, with its dual-field mode, is now capable of high contrast observations, enabling the direct observation of exoplanets. We demonstrate the technique on HR8799, a young planetary system composed of four known giant exoplanets. We used the GRAVITY fringe tracker to lock the fringes on the central star, and integrated off-axis on the HR8799e planet situated at 390 mas from the star. Data reduction included post-processing to remove the flux leaking from the central star and to extract the coherent flux of the planet. The inferred K band spectrum of the planet has a spectral resolution of 500. We also derive the astrometric position of the planet relative to the star with a precision on the order of 100$\,\mu$as. The GRAVITY astrometric measurement disfavors perfectly coplanar stable orbital solutions. A small adjustment of a few degrees to the orbital inclination of HR 8799 e can resolve the tension, implying that the orbits are close to, but not strictly coplanar. The spectrum, with a signal-to-noise ratio of $\approx 5$ per spectral channel, is compatible with a late-type L brown dwarf. Using Exo-REM synthetic spectra, we derive a temperature of $1150\pm50$\,K and a surface gravity of $10^{4.3\pm0.3}\,$cm/s$^{2}$. This corresponds to a radius of $1.17^{+0.13}_{-0.11}\,R_{\rm Jup}$ and a mass of $10^{+7}_{-4}\,M_{\rm Jup}$, which is an independent confirmation of mass estimates from evolutionary models. Our results demonstrate the power of interferometry for the direct detection and spectroscopic study of exoplanets at close angular separations from their stars.

During its orbit around the four million solar mass black hole Sagittarius A* the star S2 experiences significant changes in gravitational potential. We use this change of potential to test one part of the Einstein equivalence principle: the local position invariance (LPI). We study the dependency of different atomic transitions on the gravitational potential to give an upper limit on violations of the LPI. This is done by separately measuring the redshift from hydrogen and helium absorption lines in the stellar spectrum during its closest approach to the black hole. For this measurement we use radial velocity data from 2015 to 2018 and combine it with the gravitational potential at the position of S2, which is calculated from the precisely known orbit of S2 around the black hole. This results in a limit on a violation of the LPI of $|\beta_{He}-\beta_{H}| = (2.4 \pm 5.1) \cdot 10^{-2}$. The variation in potential that we probe with this measurement is six magnitudes larger than possible for measurements on Earth, and a factor ten larger than in experiments using white dwarfs. We are therefore testing the LPI in a regime where it has not been tested before.

Context. The properties of the accretion flow surrounding the supermassive central black hole of the Galaxy, Sgr A*, will be scrutinized by the new-generation instrument GRAVITY and the Event Horizon Telescope (EHT). Developing fast, robust, and simple models of such flows is thus important and very timely. Aims. We want to model the quiescent emission of Sgr A* from radio to mid-infrared by considering a magnetized compact torus and an extended jet. Results. We find perfect spectral fit both for face-on and edge-on views. These best fits give parameters values very close to that found by the most recent numerical simulations, which are much more complex than our model. The intrinsic radio size of Sgr A* is found to be in reasonable agreement with the centimetric observed constraints. Our best-fit infrared spectral index is in perfect agreement with the latest constraints. Our emission region at 1.3 mm, although larger than the Doeleman et al. (2008) Gaussian best-fit, does contain bright features at the 40 microarcsec scale. EHT-reconstructed images show that torus/jet-specific features persist after the reconstruction procedure, and that these features are sensitive to inclination. Conclusions. The main interest of our model is to give a simple and fast model of the quiescent state of Sgr A*, which gives extremely similar results as compared to state-of-the-art numerical simulations. Our model is easy to use and we publish all the material necessary to reproduce our spectra and images, so that anyone interested can use our results rather straightforwardly. We hope that such a public tool can be useful in the context of the recent and near-future GRAVITY and EHT results. Our model can in particular be easily used to test alternative compact objects models, or alternative gravity theories. The main limitation of our model is that we do not yet treat the X-ray emission.

The grand challenges of contemporary fundamental physics---dark matter, dark energy, vacuum energy, inflation and early universe cosmology, singularities and the hierarchy problem---all involve gravity as a key component. And of all gravitational phenomena, black holes stand out in their elegant simplicity, while harbouring some of the most remarkable predictions of General Relativity: event horizons, singularities and ergoregions. The hitherto invisible landscape of the gravitational Universe is being unveiled before our eyes: the historical direct detection of gravitational waves by the LIGO-Virgo collaboration marks the dawn of a new era of scientific exploration. Gravitational-wave astronomy will allow us to test models of black hole formation, growth and evolution, as well as models of gravitational-wave generation and propagation. It will provide evidence for event horizons and ergoregions, test the theory of General Relativity itself, and may reveal the existence of new fundamental fields. The synthesis of these results has the potential to radically reshape our understanding of the cosmos and of the laws of Nature. The purpose of this work is to present a concise, yet comprehensive overview of the state of the art in the relevant fields of research, summarize important open problems, and lay out a roadmap for future progress.

In this paper we describe the potential of the enhanced X-ray Timing and Polarimetry (eXTP) mission for studies related to accretion flows in the strong field gravity regime around both stellar-mass and supermassive black-holes. eXTP has the unique capability of using advanced 'spectral-timing-polarimetry' techniques to analyze the rapid variations with three orthogonal diagnostics of the flow and its geometry, yielding unprecedented insight into the inner accreting regions, the effects of strong field gravity on the material within them and the powerful outflows which are driven by the accretion process.

We report the detection of continuous positional and polarization changes of the compact source SgrA* in high states ('flares') of its variable near- infrared emission with the near-infrared GRAVITY-Very Large Telescope Interferometer (VLTI) beam-combining instrument. In three prominent bright flares, the position centroids exhibit clockwise looped motion on the sky, on scales of typically 150 micro-arcseconds over a few tens of minutes, corresponding to about 30% the speed of light. At the same time, the flares exhibit continuous rotation of the polarization angle, with about the same 45(+/-15)-minute period as that of the centroid motions. Modelling with relativistic ray tracing shows that these findings are all consistent with a near face-on, circular orbit of a compact polarized 'hot spot' of infrared synchrotron emission at approximately six to ten times the gravitational radius of a black hole of 4 million solar masses. This corresponds to the region just outside the innermost, stable, prograde circular orbit (ISCO) of a Schwarzschild-Kerr black hole, or near the retrograde ISCO of a highly spun-up Kerr hole. The polarization signature is consistent with orbital motion in a strong poloidal magnetic field.

This work presents an interferometric study of the massive-binary fraction in the Orion Trapezium Cluster with the recently comissioned GRAVITY instrument. We observe a total of 16 stars of mainly OB spectral type. We find three previously unknown companions for $\theta ^1$ Ori B, $\theta ^2$ Ori B, and $\theta ^2$ Ori C. We determine a separation for the previously suspected companion of NU Ori. We confirm four companions for $\theta ^1$ Ori A, $\theta ^1$ Ori C, $\theta ^1$ Ori D, and $\theta ^2$ Ori A, all with substantially improved astrometry and photometric mass estimates. We refine the orbit of the eccentric high-mass binary $\theta ^1$ Ori C and we are able to derive a new orbit for $\theta ^1$ Ori D. We find a system mass of 21.7 $M_{\odot}$ and a period of $53$ days. Together with other previously detected companions seen in spectroscopy or direct imaging, eleven of the 16 high-mass stars are multiple systems. We obtain a total number of 22 companions with separations up to 600 AU. The companion fraction of the early B and O stars in our sample is about 2, significantly higher than in earlier studies of mostly OB associations. The separation distribution hints towards a bimodality. Such a bimodality has been previously found in A stars, but rarely in OB binaries, which up to this point have been assumed to be mostly compact with a tail of wider companions. We also do not find a substantial population of equal-mass binaries. The observed distribution of mass ratios declines steeply with mass, and like the direct star counts, indicates that our companions follow a standard power law initial mass function. Again, this is in contrast to earlier findings of flat mass ratio distributions in OB associations. We exclude collision as a dominant formation mechanism but find no clear preference for core accretion or competitive accretion.

Eta Car is one of the most intriguing luminous blue variables in the Galaxy. Observations and models at different wavelengths suggest a central binary with a 5.54 yr period residing in its core. 2D and 3D radiative transfer and hydrodynamic simulations predict a primary with a dense and slow stellar wind that interacts with the faster and lower density wind of the secondary. The wind-wind collision scenario suggests that the secondary's wind penetrates the primary's wind creating a low-density cavity in it, with dense walls where the two winds interact. We aim to trace the inner ~5-50 au structure of Eta Car's wind-wind interaction, as seen through BrG and, for the first time, through the He I 2s-2p line. We have used spectro-interferometric observations with GRAVITY at the VLTI. Our modeling of the continuum allows us to estimate its FWHM angular size close to 2 mas and an elongation ratio of 1.06 +/- 0.05 over a PA = 130 +/- 20 deg. Our CMFGEN modeling helped us to confirm that the role of the secondary should be taken into account to properly reproduce the observed BrG and He I lines. Chromatic images across BrG reveal a southeast arc-like feature, possibly associated to the hot post-shocked winds flowing along the cavity wall. The images of He I 2s-2p served to constrain the 20 mas structure of the line-emitting region. The observed morphology of He I suggests that the secondary is responsible for the ionized material that produces the line profile. Both the BrG and the He I 2s-2p maps are consistent with previous hydrodynamical models of the colliding wind scenario. Future dedicated simulations together with an extensive interferometric campaign are necessary to refine our constraints on the wind and stellar parameters of the binary, which finally will help us predict the evolutionary path of Eta Car.

The highly elliptical, 16-year-period orbit of the star S2 around the massive black hole candidate Sgr A* is a sensitive probe of the gravitational field in the Galactic centre. Near pericentre at 120 AU, ~1400 Schwarzschild radii, the star has an orbital speed of ~7650 km/s, such that the first-order effects of Special and General Relativity have now become detectable with current capabilities. Over the past 26 years, we have monitored the radial velocity and motion on the sky of S2, mainly with the SINFONI and NACO adaptive optics instruments on the ESO Very Large Telescope, and since 2016 and leading up to the pericentre approach in May 2018, with the four-telescope interferometric beam-combiner instrument GRAVITY. From data up to and including pericentre, we robustly detect the combined gravitational redshift and relativistic transverse Doppler effect for S2 of z ~ 200 km/s / c with different statistical analysis methods. When parameterising the post-Newtonian contribution from these effects by a factor f, with f = 0 and f = 1 corresponding to the Newtonian and general relativistic limits, respectively, we find from posterior fitting with different weighting schemes f = 0.90 +/- 0.09 (stat) +\- 0.15 (sys). The S2 data are inconsistent with pure Newtonian dynamics.

SAGE (SagnAc interferometer for Gravitational wavE) is a fast track project for a space observatory based on multiple 12-U CubeSats in geostationary orbit. The objective of this project is to create a Sagnac interferometer with 73000 km circular arms. The geometry of the interferometer makes it especially sensitive to circularly polarized gravitational waves at frequency close to 1 Hz. The nature of the Sagnac measurement makes it almost insensitive to position error, allowing spacecrafts in ballistic trajectory. The light source and recombination units of the interferometer are based on compact fibered technologies, without the need of an optical bench. The main limitation would come from non-gravitational acceleration of the spacecraft. However, conditionally upon our ability to post-process the effect of solar wind, solar pressure and thermal expansion, we would detect gravitational waves with strains down to 10^-21 over a few days of observation.

We show that the rotating generalization of Hayward's non-singular black hole previously studied in the literature is geodesically incomplete, and that its straightforward extension leads to a singular spacetime. We present another extension, which is devoid of any curvature singularity. The obtained metric depends on three parameters and, depending on their values, yields an event horizon or not. These two regimes, named respectively regular rotating Hayward black hole and naked rotating wormhole, are studied both numerically and analytically. In preparation for the upcoming results of the Event Horizon Telescope, the images of an accretion torus around Sgr A*, the supermassive object at the center of the Galaxy, are computed. These images contain, even in the absence of a horizon, a central faint region which bears a resemblance to the shadow of Kerr black holes and emphasizes the difficulty of claiming the existence of an event horizon from the analysis of strong-field images. The frequencies of the co- and contra-rotating orbits at the innermost stable circular orbit (ISCO) in this geometry are also computed, in the hope that quasi-periodic oscillations may permit to compare this model with Kerr's black hole on observational grounds.

Thermal dominated X-ray spectra of neutron stars in quiescent transient X-ray binaries and neutron stars that undergo thermonuclear bursts are sensitive to mass and radius. The mass-radius relation of neutron stars depends on the equation of state that governs their interior. Constraining this relation accurately is thus of fundamental importance to understand the nature of dense matter. In this context we introduce a pipeline to calculate realistic model spectra of rotating neutron stars with hydrogen and helium atmospheres. An arbitrarily fast rotating neutron star with a given equation of state generates the spacetime in which the atmosphere emits radiation. We use the Lorene/nrotstar code to compute the spacetime numerically and the ATM24 code to solve the radiative transfer equations self-consistently. Emerging specific intensity spectra are then ray-traced through the neutron star's spacetime from the atmosphere to a distant observer with the Gyoto code. Here, we present and test our fully relativistic numerical pipeline. To discuss and illustrate the importance of realistic atmosphere models we compare our model spectra to simpler models like the commonly used isotropic color-corrected blackbody emission. We highlight the importance of considering realistic model-atmosphere spectra together with relativistic ray tracing to obtain accurate predictions. We also insist on the crucial impact of the star's rotation on the observables. Finally, we close a controversy that has been appearing in the literature in the recent years regarding the validity of the ATM24 code.

The goal of this paper is to investigate the detection by GRAVITY of different relativistic effects affecting the astrometric and/or spectroscopic observations of S2 such as the transverse Doppler shift, the gravitational redshift, the pericenter advance and higher-order general relativistic (GR) effects, in particular the Lense-Thirring effect due to the angular momentum of the black hole. We showed that the combination of S2 observations obtained with the GRAVITY instrument and the spectrograph SINFONI (Spectrograph for INtegral Field Observations in the Near Infrared) also installed at the VLT (Very Large Telescope) will lead to the detection of various relativistic effects. Such detections will be possible with S2 monitorings obtained within a few months or years, depending on the effect. Strong constraints on the angular momentum of Sgr~A* (e.g., at $1\sigma~=~0.1$) with the S2 star will be possible with a simple stellar-orbit model without using a ray-tracing code but with approximating the gravitational lensing effect. However, long monitorings are necessary, and we thus must rely on the discovery of closer-in stars near Sgr~A* if we want to efficiently constrain the black hole parameters with stellar orbits in a short time, or monitor the flares if they orbit around the black hole.

The second-generation beam combiner at the Very Large Telescope (VLT), GRAVITY, observes the stars orbiting the compact object located at the center of our galaxy, with an unprecedented astrometric accuracy of 10 $\mu$as. The nature of this compact source is still unknown since black holes are not the only candidates explaining the four million solar masses at the Galactic center. Boson stars are such an alternative model to black holes. This paper focuses on the study of trajectories of stars orbiting a boson star and a Kerr black hole. We put in light strong differences between orbits obtained in both metrics when considering stars with sufficiently close pericenters to the compact object, typically $\lesssim 30~M$. Discovery of closer stars to the Galactic center than the S2 star by the GRAVITY instrument would thus be a powerful tool to possibly constrain the nature of the central source.

To investigate the inner regions of protoplanetary disks, we performed near-infrared interferometric observations of the classical TTauri binary system S CrA. We present the first VLTI-GRAVITY high spectral resolution ($R\sim$4000) observations of a classical TTauri binary, S CrA (composed of S CrA N and S CrA S and separated by $\sim$1.4"), combining the four 8-m telescopes in dual-field mode. Our observations in the near-infrared K-band continuum reveal a disk around each binary component, with similar half-flux radii of about 0.1 au at d$\sim$130 pc, inclinations ($i=$28$\pm$3$^o$\ and $i=$22$\pm$6$^o$), and position angles (PA=0$^o\pm$6$^o$ and PA=-2$^o\pm$12$^o$), suggesting that they formed from the fragmentation of a common disk. The S CrA N spectrum shows bright HeI and Br$\gamma$ line emission exhibiting inverse P-Cygni profiles, typically associated with infalling gas. The continuum-compensated Br$\gamma$ line visibilities of S CrA N show the presence of a compact Br$\gamma$ emitting region the radius of which is about $\sim$0.06 au, which is twice as big as the truncation radius. This component is mostly tracing a wind. Moreover, a slight radius change between the blue- and red-shifted Br$\gamma$ line components is marginally detected. The presence of an inverse P-Cygni profile in the HeI and Br$\gamma$ lines, along with the tentative detection of a slightly larger size of the blue-shifted Br$\gamma$ line component, hint at the simultaneous presence of a wind and magnetospheric accretion in S CrA N.

GRAVITY is a new instrument to coherently combine the light of the European Southern Observatory Very Large Telescope Interferometer to form a telescope with an equivalent 130 m diameter angular resolution and a collecting area of 200 m$^2$. The instrument comprises fiber fed integrated optics beam combination, high resolution spectroscopy, built-in beam analysis and control, near-infrared wavefront sensing, phase-tracking, dual beam operation and laser metrology [...]. This article gives an overview of GRAVITY and reports on the performance and the first astronomical observations during commissioning in 2015/16. We demonstrate phase tracking on stars as faint as m$_K$ ~ 10 mag, phase-referenced interferometry of objects fainter than m$_K$ ~ 15 mag with a limiting magnitude of m$_K$ ~ 17 mag, minute long coherent integrations, a visibility accuracy of better than 0.25 %, and spectro-differential phase and closure phase accuracy better than 0.5{\deg}, corresponding to a differential astrometric precision of better than 10 microarcseconds ({\mu}as). The dual-beam astrometry, measuring the phase difference of two objects with laser metrology, is still under commissioning. First observations show residuals as low as 50 {\mu}as when following objects over several months. We illustrate the instrument performance with the observations of archetypical objects for the different instrument modes. Examples include the Galactic Center supermassive black hole and its fast orbiting star S2 for phase referenced dual beam observations and infrared wavefront sensing, the High Mass X-Ray Binary BP Cru and the Active Galactic Nucleus of PDS 456 for few {\mu}as spectro-differential astrometry, the T Tauri star S CrA for a spectro-differential visibility analysis, {\xi} Tel and 24 Cap for high accuracy visibility observations, and {\eta} Car for interferometric imaging with GRAVITY.

We observe the HMXB BP Cru using interferometry in the near-infrared K band with VLTI/GRAVITY. Continuum visibilities are at most partially resolved, consistent with the predicted size of the hypergiant. Differential visibility amplitude ($\Delta |V| \sim 5\%$) and phase ($\Delta \phi \sim 2 \deg$) signatures are observed across the HeI $2.059 \mu$m and Br$\gamma$ lines, the latter seen strongly in emission, unusual for the donor star's spectral type. For a baseline $B \sim 100$m, the differential phase RMS $\sim 0.2 \deg$ corresponds to an astrometric precision of $\sim 2 \mu$as. A model-independent analysis in the marginally resolved limit of interferometry reveals asymmetric and extended emission with a strong wavelength dependence. We propose geometric models based on an extended and distorted wind and/or a high density gas stream, which has long been predicted to be present in this system. The observations show that optical interferometry is now able to resolve HMXBs at the spatial scale at which accretion takes place, and therefore probe the effects of the gravitational and radiation fields of the compact object on its environment.

We present the first optical observation at sub-milliarcsecond (mas) scale of the microquasar SS 433 obtained with the GRAVITY instrument on the VLT interferometer. The 3.5 hour exposure reveals a rich K-band spectrum dominated by hydrogen Br$\gamma\ $ and \ion{He}{i} lines, as well as (red-shifted) emission lines coming from the jets. The K-band continuum emitting region is dominated by a marginally resolved point source ($<$ 1 mas) embedded inside a diffuse background accounting for 10\% of the total flux. The jet line positions agree well with the ones expected from the jet kinematic model, an interpretation also supported by the consistent sign (i.e. negative/positive for the receding/approaching jet component) of the phase shifts observed in the lines. The significant visibility drop across the jet lines, together with the small and nearly identical phases for all baselines, point toward a jet that is offset by less than 0.5 mas from the continuum source and resolved in the direction of propagation, with a typical size of 2 mas. The jet position angle of $\sim$80$^{\circ}$ is consistent with the expected one at the observation date. Jet emission so close to the central binary system would suggest that line locking, if relevant to explain the amplitude and stability of the 0.26c jet velocity, operates on elements heavier than hydrogen. The Br$\gamma\ $ profile is broad and double peaked. It is better resolved than the continuum and the change of the phase signal sign across the line on all baselines suggests an East-West oriented geometry alike the jet direction and supporting a (polar) disk wind origin.

We address the astrophysical imaging of a family of deformed Kerr black holes (BHs). These are stationary, asymptotically flat black hole (BH) spacetimes, that are solutions of General Relativity minimally coupled to a massive, complex scalar field: Kerr BHs with scalar hair (KBHsSH). Such BHs bifurcate from the vacuum Kerr solution and can be regarded as a horizon within a rotating boson star. In a recent letter, it was shown that KBHsSH can exhibit very distinct shadows from the ones of their vacuum counterparts. The setup therein, however, considered the light source to be a celestial sphere sufficiently far away from the BH. Here, we analyse KBHsSH surrounded by an emitting torus of matter, simulating a more realistic astrophysical environment, and study the corresponding lensing of light as seen by a very far away observer, to appropriately model ground-based observations of Sgr A*. We find that the differences in imaging between KBHsSH and comparable vacuum Kerr BHs remain, albeit less dramatic than those observed for the corresponding shadows in the previous setup. In particular, we highlight two observables that might allow differentiating KBHsSH and Kerr BHs. The first is the angular size of the photon ring (in a Kerr spacetime) or lensing ring (in a KBHSH spacetime), the latter being significantly smaller for sufficiently non-Kerr-like spacetimes. The second is the existence of an edge in the intensity distribution (the photon ring in Kerr spacetime). This edge can disappear for very non-Kerr-like KBHsSH. It is plausible, therefore, that sufficiently precise Very Long Baseline Interferometric observations of BH candidates can constrain this model.

Context. Quasi-periodic oscillations (QPOs) are a common feature of the power spectrums of X-ray binaries. Currently it is not possible to unambiguously differentiate the large number of proposed models to explain these phenomena through existing observations. Aims. We investigate the observable predictions of a simple model that generates flux modulation: a spiral instability rotating in a thin accretion disk. This model is motivated by the accretion ejection instability (AEI) model for low- frequency QPOs (LFQPOs). We are particularly interested in the inclination dependence of the observables that are associated with this model. Methods. We develop a simple analytical model of an accretion disk, which features a spiral instability. The disk is assumed to emit blackbody radiation, which is ray-traced to a distant observer. We compute pulse profiles and power spectra as observed from infinity. Results. We show that the amplitude of the modulation associated with the spiral rotation is a strong function of inclination and frequency. The pulse profile is quasi-sinusoidal only at low inclination (face-on source). As a consequence, a higher-inclination geometry leads to a stronger and more diverse harmonic signature in the power spectrum. Conclusions. We present how the amplitude depends on the inclination when the flux modulation comes from a spiral in the disk. We also include new observables that could potentially differentiate between models, such as the pulse profile and the harmonic content of the power spectra of high-inclination sources that exhibit LFQPOs. These might be important observables to explore with existing and new instruments.

[Abridged] Context. The high-energy radiation from black-hole binaries may be due to the reprocessing of a lamp located on the black hole axis, emitting X-rays. The observed spectrum is made of 3 components: the direct spectrum; the thermal bump; and the reflected spectrum made of the Compton hump and the iron-line complex. Aims. We aim at computing accurately the complete reprocessed spectrum (thermal bump + reflected) of black-hole binaries over the entire X-ray band. We also determine the strength of the direct component. Our choice of parameters is adapted to a source showing an important thermal component. Methods. We compute in full GR the illumination of a thin disk by a lamp along the rotation axis. We use the ATM21 radiative transfer code to compute the spectrum emitted along the disk. We ray trace this local spectrum to determine the reprocessed spectrum as observed at infinity. We discuss the dependence of the local and ray-traced spectra on the emission angle and spin. Results. We show the importance of the angle dependence of the total disk specific intensity spectrum emitted by the illuminated atmosphere when the thermal disk emission if fully taken into account. High spin implies high temperature in the inner regions, so the emitted thermal disk spectrum covers the iron-line complex. As a result we locally observe absorption lines produced in the hot disk atmosphere. Absorption lines are narrow and disappear after ray tracing the local spectrum. Conclusions. Our results mainly highlight the importance of considering the angle dependence of the local spectrum when computing reprocessed spectra, as was already found in a recent study. The main new result of our work is to show the importance of computing the thermal bump of the spectrum, as this feature can change considerably the observed iron-line complex.

Millimeter very long baseline interferometry will soon produce accurate images of the closest surroundings of the supermassive compact object at the center of the Galaxy, Sgr A*. These images may reveal the existence of a central faint region, the so-called shadow, which is often interpreted as the observable consequence of the event horizon of a black hole. In this paper, we compute images of an accretion torus around Sgr A* assuming this compact object is a boson star, i.e. an alternative to black holes within general relativity, with no event horizon and no hard surface. We show that very relativistic rotating boson stars produce images extremely similar to Kerr black holes, showing in particular shadow-like and photon-ring-like structures. This result highlights the extreme difficulty of unambiguously telling the existence of an event horizon from strong-field images.

We simulate an oscillating purely hydrodynamical torus with constant specific angular mo- mentum around a Schwarzschild black hole. The goal is to search for quasi-periodic oscil- lations (QPOs) in the light curve of the torus. The initial torus setup is subjected to radial, vertical and diagonal (combination of radial and vertical) velocity perturbations. The hydro- dynamical simulations are performed using the general relativistic magnetohydrodynamics code Cosmos++ and ray-traced using the GYOTO code. We found that a horizontal velocity perturbation triggers the radial and plus modes, while a vertical velocity perturbation trig- gers the vertical and X modes. The diagonal perturbation gives a combination of the modes triggered in the radial and vertical perturbations.

Accretion disks play an important role in the evolution of their relativistic inner compact objects. The emergence of a new generation of interferometers will allow to resolve these accretion disks and provide more information about the properties of the central gravitating object. Due to this instrumental leap forward it is crucial to investigate the accretion disk physics near various types of inner compact objects now to deduce later constraints on the central objects from observations. A possible candidate for the inner object is the boson star. Here, we will try to analyze the differences between accretion structures surrounding boson stars and black holes. We aim at analysing the physics of circular geodesics around boson stars and study simple thick accretion tori (so-called Polish doughnuts) in the vicinity of these stars. We realize a detailed study of the properties of circular geodesics around boson stars. We then perform a parameter study of thick tori with constant angular momentum surrounding boson stars. This is done using the boson star models computed by a code constructed with the spectral solver library KADATH. We demonstrate that all the circular stable orbits are bound. In the case of a constant angular momentum torus, a cusp in the torus surface exists only for boson stars with a strong gravitational scalar field. Moreover, for each inner radius of the disk, the allowed specific angular momentum values lie within a constrained range which depends on the boson star considered. We show that the accretion tori around boson stars have different characteristics than in the vicinity of a black hole. With future instruments it could be possible to use these differences to constrain the nature of compact objects.

In Spring 2011 we observed Sgr A*, the supermassive black hole at the center of our Galaxy, with XMM-Newton with a total exposure of ~226 ks in coordination with the 1.3 mm VLBI. We have performed timing analysis of the X-ray emission from Sgr A* using Bayesian blocks algorithm to detect X-ray flares observed with XMM-Newton. Furthermore, we computed X-ray smoothed light curves observed in this campaign in order to have better accuracy on the position and the amplitude of the flares. We detected 2 X-ray flares on the 2011 March 30 and April 3 which have for comparison a peak detection level of 6.8 and 5.9 sigma in the XMM-Newton/EPIC light curve in the 2-10 keV energy range with a 300 s bin. The former is characterized by 2 sub-flares: the first one is very short (~458 s) with a peak luminosity of ~9.4E34 erg/s whereas the second one is longer (~1542 s) with a lower peak luminosity of ~6.8E34 erg/s. The comparison with the sample of X-ray flares detected during the 2012 Chandra XVP campaign favors the hypothesis that the 2011 March 30 flare is a single flare rather than 2 distinct sub-flares. We model the light curve of this flare with the gravitational lensing of a simple hotspot-like structure but we can not satisfactorily reproduce the large decay of the light curve between the 2 sub-flares with this model. From magnetic energy heating during the rise phase of the first sub-flare and assuming an X-ray photons production efficiency of 1 and a magnetic field of 100 G at 2 r_g, we derive an upper limit to the radial distance of the first sub-flare of 100 r_g. We estimate using the decay phase of the first sub-flare a lower limit to the radial distance of 4 r_g from synchrotron cooling in the infrared. The X-ray emitting region of the first sub-flare is located at a radial position of 4-100 and has a corresponding radius of 1.8-2.87 in r_g unit for a magnetic field of 100 G at 2 r_g.

High frequency (millisecond) quasi-periodic oscillations (HF QPOs) are observed in the X-ray power-density spectra of several microquasars and low mass X-ray binaries. Two distinct QPO peaks, so-called twin peak QPOs, are often detected simultaneously exhibiting their frequency ratio close or equal to 3/2. Following the analytic theory and previous studies of observable spectral signatures, we aim to model the twin peak QPOs as a spectral imprint of specific dual oscillation regime defined by a combination of the lowest radial and vertical oscillation mode of optically thick slender tori with constant specific angular momentum. We examined power spectra and fluorescent K$\alpha$ iron line profiles for two different simulation setups with the mode frequency relations corresponding to the epicyclic resonance HF QPOs model and modified relativistic precession QPOs model. We use relativistic ray-tracing implemented in parallel simulation code LSDplus. In the background of the Kerr spacetime geometry, we analyze the influence of the distant observer inclination and the spin of the central compact object. Relativistic optical projection of the oscillating slender torus is illustrated by images in false colours related to the frequency shift. We show that performed simulations yield power spectra with the pair of dominant peaks corresponding to the frequencies of radial and vertical oscillation modes with the proper ratio equal to 3/2 on a wide range of inclinations and spin values. We also discuss exceptional cases of a very small and very high inclination as well as unstable high spin relativistic precession-like configuration predicting constant frequency ratio equal to 1/2. We demonstrate signifiant dependency of broadened K$\alpha$ iron line profiles on the inclination of the distant observer.

This is a White Paper in support of the mission concept of the Large Observatory for X-ray Timing (LOFT), proposed as a medium-sized ESA mission. We discuss the potential of LOFT for the study of active galactic nuclei. For a summary, we refer to the paper.

High frequency quasi-periodic oscillations (HF QPOs) are observed in the X-ray power-density spectra (PDS) of several microquasars and low mass X-ray binaries. Many proposed QPO models are based on oscillations of accretion toroidal fluid structures orbiting in the vicinity of a compact object. We study oscillating accretion tori orbiting in the vicinity of a Kerr black hole. We demonstrate that significant variation of the observed flux can be caused by the combination of radial and vertical oscillation modes of a slender, polytropic, perfect fluid, non-self-graviting torus with constant specific angular momentum. We investigate two combinations of the oscillating modes corresponding to the direct resonance QPO model and the modified relativistic precession QPO model.

Context. The supermassive black hole, Sagittarius (Sgr) A*, in the centre of our Galaxy has the largest angular size in the sky among all astrophysical black holes. Its shadow, assuming no rotation, spans ~ 50 microarcsec. Resolving such dimensions has long been out of reach for astronomical instruments until a new generation of interferometers being operational during this decade. Of particular interest is the Event Horizon Telescope (EHT) with resolution ~ 20 microarcsec in the millimeter-wavelength range 0.87 mm - 1.3 mm. Aims. We investigate the ability of the fully general relativistic Komissarov (2006) analytical magnetized torus model to account for observable constraints at Sgr A* in the centimeter and millimeter domains. The impact of the magnetic field geometry on the observables is also studied. Methods. We calculate ray-traced centimeter- and millimeter-wavelength synchrotron spectra and images of a magnetized accretion torus surrounding the central black hole in Sgr A*. We assume stationarity, axial symmetry, constant specific angular momentum and polytropic equation of state. A hybrid population of thermal and non-thermal electrons is considered. Results. We show that the torus model is capable of reproducing spectral constraints in the millimeter domain, and in particular in the observable domain of the EHT. However, the torus model is not yet able to fit the centimeter spectrum. 1.3 mm images at high inclinations are in agreement with observable constraints. Conclusions. The ability of the torus model to account for observations of Sgr A* in the millimeter domain is interesting in the perspective of the future EHT. Such an analytical model allows very fast computations. It will thus be a suitable test bed for investigating large domains of physical parameters, as well as non-black-hole compact object candidates and alternative theories of gravity.

Very long baseline interferometry at millimetre/submillimetre wavelengths (mmVLBI) offers the highest achievable spatial resolution at any wavelength in astronomy. The anticipated inclusion of ALMA as a phased array into a global VLBI network will bring unprecedented sensitivity and a transformational leap in capabilities for mmVLBI. Building on years of pioneering efforts in the US and Europe the ongoing ALMA Phasing Project (APP), a US-led international collaboration with MPIfR-led European contributions, is expected to deliver a beamformer and VLBI capability to ALMA by the end of 2014 (APP: Fish et al. 2013, arXiv:1309.3519). This report focuses on the future use of mmVLBI by the international users community from a European viewpoint. Firstly, it highlights the intense science interest in Europe in future mmVLBI observations as compiled from the responses to a general call to the European community for future research projects. A wide range of research is presented that includes, amongst others: - Imaging the event horizon of the black hole at the centre of the Galaxy - Testing the theory of General Relativity an/or searching for alternative theories - Studying the origin of AGN jets and jet formation - Cosmological evolution of galaxies and BHs, AGN feedback - Masers in the Milky Way (in stars and star-forming regions) - Extragalactic emission lines and astro-chemistry - Redshifted absorption lines in distant galaxies and study of the ISM and circumnuclear gas - Pulsars, neutron stars, X-ray binaries - Testing cosmology - Testing fundamental physical constants

The black hole at the Galactic centre exhibits regularly flares of radiation, the origin of which is still not understood. In this article, we study the ability of the near-future GRAVITY infrared instrument to constrain the nature of these events. We develop realistic simulations of GRAVITY astrometric data sets for various flare models. We show that the instrument will be able to distinguish an ejected blob from alternative flare models, provided the blob inclination is >= 45deg, the flare brightest magnitude is 14 <= mK <= 15 and the flare duration is >= 1h30.

In the context of high-frequency quasi-periodic oscillation (HF QPOs) we further explore the appearance of an observable signal generated by hot spots moving along quasi-elliptic trajectories close to the innermost stable circular orbit in the Schwarzschild spacetime. The aim of our investigation is to reveal whether observable characteristics of the Fourier power-spectral density can help us to distinguish between the two competing models, namely, the idea of bright spots orbiting on the surface of an accretion torus versus the scenario of intrinsic oscillations of the torus itself. We take the capabilities of the present observatories (represented by the Rossi X-ray Timing Explorer, RXTE) into account, and we also consider the proposed future instruments (represented here by the Large Observatory for X-ray Timing, LOFT).

Context. Some microquasars exhibit millisecond quasi-periodic oscillations (QPO) that are likely related to phenomena occuring in the immediate vicinity of the central black hole. Oscillations of accretion tori have been proposed to model these QPOs. Aims. Here, we aim at determining the observable spectral signature of slender accretion tori surrounding Kerr black holes. We analyze the impact of the inclination and spin parameters on the power spectra. Methods. Ray-traced power spectra of slender tori oscillation modes are computed in the Kerr metric. Results. We show that the power spectral densities of oscillating tori are very sensitive to the inclination and spin parameters. This strong dependency of the temporal spectra on inclination and spin may lead to observable constraints of these parameters. Conclusions. This work goes a step further in the analysis of the oscillating torus QPO model. It is part of a long-term study that will ultimately lead to comparison with observed data.

High frequency quasi-periodic oscillations (HF QPOs) appear in the X-ray variability of several accreting low-mass binaries. In a series of works it was suggested that these QPOs may have connection to inhomogeneities orbiting close to an inner edge of the accretion disc. In this paper we explore the appearance of an observable signal generated by small radiating circular hot spots moving along quasi-elliptic trajectories close to the innermost stable circular orbit in the Schwarzschild spacetime. Our consideration takes into account the capabilities of observatories that have been operating in the past two decades represented by the Rossi X-ray Timing Explorer (RXTE) and the proposed future instruments represented by the Large Observatory for X-ray Timing (LOFT). For these purposes we choose such model parameters that lead to lightcurves comparable to those observed in Galactic black hole sources, in particular the microquasar GRS 1915+105. We find that when a weak signal corresponding to the hot-spot Keplerian frequency is around the limits of the RXTE detectability, the LOFT observations can clearly reveal its first and second harmonics. Moreover, in some specific situations the radial epicyclic frequency of the spot can be detected as well. Finally, we also compare the signal produced by the spots to the signal produced by axisymmetric epicyclic disc-oscillation modes and discuss the key differences that could be identified via the proposed future technology. We conclude that the ability to recognize the harmonic content of the signal can help to distinguish between the different proposed physical models.

We investigate the possibility to distinguish the small-coupling, slow-rotation black hole solution of Chern-Simons (CS) gravity from the Kerr solution. We develop simulations of electromagnetic observables in the vicinity of CS and Kerr black holes. We show that the typical relative observable difference between CS and Kerr spacetimes is of the order of 0.1% thus beyond reach of current or near-future instruments.

We study the Rossby wave instability model of high-frequency quasi-periodic oscillations (QPO) of microquasars. We show ray-traced light curves of QPO within this model and discuss perspectives of distinguishing alternative QPO models with the future Large Observatory For X-ray Timing (LOFT) observations.

Context. One of the often discussed models for X-ray binaries high-frequency quasi-periodic oscillations is the oscillating torus model that considers oscillation modes of slender accretion tori. Aims. Here, we aim at developing this model by considering the observable signature of an optically thick slender accretion torus subject to simple periodic deformations. Methods. We compute light curves and power spectra of a slender accretion torus subject to simple periodic deformations: vertical or radial translation, rotation, expansion and shear. Results. We show that different types of deformations lead to very different Fourier power spectra and therefore could be observationally distinguished. Conclusions. This work is a first step in a longer-term study of the observable characteristics of the oscillating torus model. It gives promising perspectives on the possibility to constrain this model by studying the observed power spectra of quasi-periodic oscillations.

Context. There have been a long string of efforts to understand the source of the variability observed in microquasars, especially concerning the elusive High-Frequency Quasi-Periodic Oscillation. These oscillations are among the fastest phenomena that affect matter in the vicinity of stellar black holes and therefore could be used as probes of strong-field general relativity. Nevertheless, no model has yet gained wide acceptance. Aims. The aim of this article is to investigate the model derived from the occurrence of the Rossby wave instability at the inner edge of the accretion disk. In particular, our goal here is to demonstrate the capacity of this instability to modulate the observed flux in agreement with the observed results. Methods. We use the AMRVAC hydrodynamical code to model the instability in a 3D optically thin disk. The GYOTO ray-tracing code is then used to compute the associated light curve. Results. We show that the 3D Rossby wave instability is able to modulate the flux well within the observed limits.We highlight that 2D simulations allow us to obtain the same general characteristics of the light curve as 3D calculations. With the time resolution we adopted in this work, three dimensional simulations do not give rise to any new observable features that could be detected by current instrumentation or archive data.

The equations governing null and timelike geodesics are derived within the 3+1 formalism of general relativity. In addition to the particle's position, they encompass an evolution equation for the particle's energy leading to a 3+1 expression of the redshift factor for photons. An important application is the computation of images and spectra in spacetimes arising from numerical relativity, via the ray-tracing technique. This is illustrated here by images of numerically computed stationary neutron stars and dynamical neutron stars collapsing to a black hole.

There have been a long string of efforts to understand the source of the variability observed in microquasars but no model has yet gained wide acceptance, especially concerning the elusive High-Frequency Quasi-Periodic Oscillation (HFQPO). We first list the constraints arising from observations and how that translates for an HFQPO model. Then we present how a model based on having the Rossby Wave Instability (RWI) active in the disk could answer those constraints.

The Galactic center is an ideal laboratory to study strong-field general relativistic phenomena, as the supermassive black hole Sgr A* has the biggest angular Schwarzschild radius among all black holes. This article presents three different ways of using the immediate surroundings of Sgr A* as a probe in order to either constrain its spin, or even test the very nature of this compact object.

We calculate the "observed at infinity" image and spectrum of the accretion structure in Sgr A*, by modelling it as an optically thin, constant angular momentum ion torus in hydrodynamic equilibrium. The physics we consider includes a two-temperature plasma, a toroidal magnetic field, as well as radiative cooling by bremsstrahlung, synchrotron and inverse Compton processes. Our relativistic model has the virtue of being fully analytic and very simple, depending only on eight tunable parameters: the black hole spin and the inclination of the spin axis to our line of sight, the torus angular momentum, the polytropic index, the magnetic to total pressure ratio, the central values of density and electron temperature and the ratio of electron to ion temperatures. The observed image and spectrum are calculated numerically using the ray-tracing code GYOTO. Our results demonstrate that the ion torus model is able to account for the main features of the accretion structure surrounding Sgr A*.

The ability of the near future second generation VLTI instrument GRAVITY to constrain the properties of the Galactic center black hole is investigated. The Galactic center infrared flares are used as probes of strong-field gravity, within the framework of the hot spot model according to which the flares are the signature of a blob of gas orbiting close to the black hole's innermost stable circular orbit. Full general relativistic computations are performed, together with realistic observed data simulations, that lead to conclude that GRAVITY could be able to constrain the black hole's inclination parameter.

GYOTO, a general relativistic ray-tracing code, is presented. It aims at computing images of astronomical bodies in the vicinity of compact objects, as well as trajectories of massive bodies in relativistic environments. This code is capable of integrating the null and timelike geodesic equations not only in the Kerr metric, but also in any metric computed numerically within the 3+1 formalism of general relativity. Simulated images and spectra have been computed for a variety of astronomical targets, such as a moving star or a toroidal accretion structure. The underlying code is open source and freely available. It is user-friendly, quickly handled and very modular so that extensions are easy to integrate. Custom analytical metrics and astronomical targets can be implemented in C++ plug-in extensions independent from the main code.

The galactic central black hole Sgr A* exhibits outbursts of radiation in the near infrared (so-called IR flares). One model of these events consists in a hotspot orbiting on the innermost stable circular orbit (ISCO) of the hole. These outbursts can be used as a probe of the central gravitational potential. One main scientific goal of the second generation VLTI instrument GRAVITY is to observe these flares astrometrically. Here, the astrometric precision of GRAVITY is investigated in imaging mode, which consists in analysing the image computed from the interferometric data. The capability of the instrument to put in light the motion of a hotspot orbiting on the ISCO of our central black hole is then discussed. We find that GRAVITY's astrometric precision for a single star in imaging mode is smaller than the Schwarzschild radius of Sgr A*. The instrument can also demonstrate that a body orbiting on the last stable orbit of the black hole is indeed moving. It yields a typical size of the orbit, if the source is as bright as m_K=14. These results show that GRAVITY allows one to study the close environment of Sgr A*. Having access to the ISCO of the central massive black hole probably allows constraining general relativity in its strong regime. Moreover, if the hotspot model is appropriate, the black hole spin can be constrained.

New observations from the Hubble ultra deep field suggest that the star formation rate at z>7 drops off faster than previously thought. Using a newly determined star formation rate for the normal mode of Population II/I stars (PopII/I), including this new constraint, we compute the Thomson scattering optical depth and find a result that is marginally consistent with WMAP5 results. We also reconsider the role of Population III stars (PopIII) in light of cosmological and stellar evolution constraints. While this input may be needed for reionization, we show that it is essential in order to account for cosmic chemical evolution in the early Universe. We investigate the consequences of PopIII stars on the local metallicity distribution function of the Galactic halo (from the recent Hamburg/ESO survey of metal-poor stars) and on the evolution of abundances with metallicity (based on the ESO large program on very metal-poor stars), with special emphasis on carbon-enhanced metal-poor stars. Our most important results show that the nucleosynthetic yields of PopIII stars lead to abundance patterns in agreement with those observed in extremely metal-poor stars. In this chemical approach to cosmic evolution, PopIII stars prove to be a compulsory ingredient, and extremely metal-poor stars are inevitably born at high redshift. (Abridged)

We use measurements of cosmic shear from CFHTLS, combined with WMAP-5 cosmic microwave background anisotropy data, baryonic acoustic oscillations from SDSS and 2dFGRS and supernovae data from SNLS and Gold-set, to constrain the neutrino mass. We obtain a 95% confidence level upper limit of 0.54 eV for the sum of the neutrino masses, and a lower limit of 0.03 eV. The preference for massive neutrinos vanishes when shear-measurement systematics are included in the analysis.

The Solar Tower Atmospheric Cherenkov Effect Experiment (STACEE) is located at the National Solar Thermal Test Facility of Sandia National Laboratories in Albuquerque, New Mexico, USA. The field of solar tracking mirrors (heliostats) around a central receiver tower is used to direct Cherenkov light from atmospheric showers onto secondary mirrors on the tower, which in turn image the light onto cameras of photomultiplier tubes. The STACEE Collaboration has previously reported a detection of the Crab Nebula with approximately 7 standard deviation significance, using 32 heliostats (STACEE-32). This result demonstrates both the viability of the technique and the suitability of the site. We are in the process of completing an upgrade to 48 heliostats (STACEE-48) en route to an eventual configuration using 64 heliostats (STACEE-64) in early 2001. In this paper, we summarize the results obtained on the sensitivity of STACEE-32 and our expectations for STACEE-48 and STACEE-64.