# Detection of orbital motions near the last stable circular orbit of the massive black hole SgrA*

R. Abuter, A. Amorim, M. Bauböck, J. Berger, H. Bonnet, W. Brandner, Y.

Clenet, V. Coudé Du Foresto, P. de Zeeuw, C. Deen, et al.

► To cite this version:

R. Abuter, A. Amorim, M. Bauböck, J. Berger, H. Bonnet, et al., Detection of orbital motions near the last stable circular orbit of the massive black hole SgrA*. Astronomy and Astrophysics - A&A, 2018, 618, pp.L10. 10.1051/0004-6361/201834294. hal-02301352

HAL Id: hal-02301352

https://hal.science/hal-02301352

Submitted on 18 Jan 2022

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
LETTER TO THE EDITOR

# Detection of orbital motions near the last stable circular orbit of the massive black hole SgrA**

GRAVITY Collaboration\*\*: R. Abuter\*\*, A. Amorim\*\*\*\*, M. Bauböck\*\*, J. P. Berger\*\*, H. Bonnet\*\*, W. Brandner\*\*, Y. Clénet\*\*, V. Coudé du Foresto\*\*, P. T. de Zeeuw\*\*\*, C. Deen\*\*, J. Dexter\*\*\*, G. Duvert\*\*, A. Eckart\*\*\*, F. Eisenhauer\*\*, N. M. Förster Schreiber\*\*, P. Garcia\*\*\*, F. Gao\*\*, E. Gendron\*\*, R. Genzel\*\*, S. Gillessen\*\*, P. Guajardo\*\*, M. Habibi\*, X. Haubois\*\*, Th. Henning\*, S. Hippler\*, M. Horrobin\*, A. Huber\*, A. Jiménez-Rosales\*, L. Jocou\*\*, P. Kervella\*\*, S. Lacour\*\*, V. Lapeyrère\*, B. Lazareff\*, J.-B. Le Bouquin\*, P. Léna\*, M. Lippa\*, T. Ott\*, J. Panduro\*\*, T. Paumard\*\*, K. Perraut\*, G. Perrin\*, O. Pfuhl\*\*, P. M. Plewa\*, S. Rabien\*, G. Rodríguez-Coira\*, G. Rousset\*, A. Sternberg\*\*, O. Straub\*, C. Straubmeier\*, E. Sturm\*, L. J. Tacconi\*, F. Vincent\*, S. von Fellenberg\*, I. Waisberg\*, F. Widmann\*, E. Wieprecht\*, E. Wiezorrek\*, J. Woillez\*, and S. Yazici\*\*\*

*(Affiliations can be found after the references)*

Received 21 September 2018 / Accepted 5 October 2018

## ABSTRACT

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 µas 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) min 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.

**Key words.** Galaxy: center – black hole physics – gravitation – relativistic processes

## 1. Introduction

The compact source SgrA* at the centre of the Milky Way harbours a concentration of 4.14 million solar masses, plausibly a massive black hole (Genzel et al. 2010; Ghez et al. 2008). SgrA* exhibits steady and continuously variable, non-thermal emission across the electromagnetic spectrum (Genzel et al. 2010; Baganoff et al. 2001; Trippe et al. 2007; Eckart et al. 2008; Yusef-Zadeh et al. 2008; Do et al. 2009; Dodds-Eden et al. 2009; Shahzamanian & Eckart 2015; Ponti et al. 2017; Witzel et al. 2018). Intercontinental microwave interferometry and polarized infrared(IR)/X-ray variability on 10–30 min timescales suggest

that this emission comes from highly relativistic electrons in a hot, magnetized accretion disk/torus of ~10 light minutes in diameter, plus perhaps a jet, just outside the innermost stable circular orbit (ISCO) of the putative massive black hole (Witzel et al. 2018; Johnson et al. 2018; Doeleman et al. 2008; Yuan et al. 2004; Markoff et al. 2001). The exploration of this innermost, relativistic accretion region with high-resolution imaging techniques promises important and fundamental information for physics and astronomy, including new stringent tests of the massive black hole paradigm.

We have been observing the Galactic centre and SgrA* with the GRAVITY instrument (Gravity Collaboration et al. 2017, 2018; Eisenhauer et al. 2011, 2008; Paumard et al. 2008) during multiple campaigns in 2017/2018\*, with the aim of testing general relativity (GR) and the massive black hole paradigm in the closest massive black hole candidate. Gravity Collaboration et al. (2018) have already reported in this journal on a high-quality measurement of the gravitational redshift in the orbit of the star S2 going through its peri-approach at 2800 $R_g$ from SgrA* ($R_g = GM_\bullet/c^2 = 6.1 \times 10^{11}$ cm, or 5 µas) in May 2018. Another

<footnote>\* The data are only available at the CDS via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via http://cdsarc.u-strasbg.fr/viz-bin/qcat?J/A+A/618/L10</footnote>

<footnote>\*\* GRAVITY is developed in a collaboration by the Max Planck Institute for Extraterrestrial Physics, LESIA of Paris Observatory/Université PSL/CNRS/Sorbonne Université/Univ. Paris Diderot/Sorbonne Paris Cité, IPAG of Université Grenoble Alpes/CNRS, the Max Planck Institute for Astronomy, the University of Cologne, the CENTRA – Centro de Astrofísica e Gravitação, and the European Southern Observatory.</footnote>

<footnote>\*\*\* Corresponding authors: O. Pfuhl (e-mail: pfuhl@mpe.mpg.de), J. Dexter (e-mail: jdexter@mpe.mpg.de) and T. Paumard (e-mail: thibaut.paumard@obspm.fr).</footnote>

<footnote n="1">¹ ESO Telescopes at the La Silla Paranal Observatory programme IDs 099.B-0162, 0100.B-0731, 0101.B-0195, and 0101.B-0576.</footnote>
<figure><img src="image_2.png" /><figcaption>**Fig. 1.** Bottom left (panel a): time evolution of the east-west (east positive, blue) and north-south (red) position offset of the July 22 flare (MJD = 58321.9954) centroids from their medians, as well as the flux density evolution (right y-axis, black) in units of the flux of S2 (14.0 mag). Error bars are $1\sigma$. For this purpose the total intensity was computed from the sum of the two polarization directions. The points represent the average of the “Pfuhl” and “Waisberg” analyses (see Appendix B). Bottom right (panel b): projected orbit of the flare centroid on the sky (colour ranging from brown to dark blue as a qualitative marker of time through the 30 min observation, relative to their medians (small black cross) and after removal of the S2 motion and differential refraction between S2 and SgrA*). The orange square and $1\sigma$ uncertainty is the long-term astrometric position of the mass centre of the S2 orbit (approximately the orbital centroid, although shifts between apparent and true centroids can be introduced by lensing, relativistic beaming, and azimuthal shearing of an initially compact “hot spot”). Top left (panel c) and top right (panel d): comparison of the data of the *bottom two panels* with a realization of a simple hot spot model in the Schwarzschild metric, including light bending, lensing, time dilation and other effects of GR and/or special relativity (SR), computed from the NERO relativistic ray tracing code (Bauböck et al., in prep.). Similar results were obtained with the GYOTO code (Vincent et al. 2011b; Grould et al. 2016). The purple and cyan continuous curves in (c) show the same orbit in $x(t)$ and $y(t)$, compared to the data in blue and red. The continuous blue curve in (d) denotes a hot spot on a circular orbit with $R = 1.17 \times R(\text{ISCO}, a = 0, M = 4.14 \times 10^6 M_{\odot})$, seen at inclination $160^\circ$ (clockwise on the sky, as for the data in (d)) and with the line of nodes at $\Omega = 160^\circ$ ($\chi_r^2 = 1.2$). Open blue circles and grey bars connect the data points to their locations on the best fit orbit.</figcaption></figure>

main goal of our observations is to search for orbital motions of “hot spots” of relativistic gas in the innermost accretion zone around the black hole’s ISCO (Broderick & Loeb 2005, 2006; Hamaus et al. 2009).

Such hot spots have been proposed to originate from magnetic shocks or re-connection events in the innermost accretion zone (Eckart et al. 2008; Zamaninasab et al. 2010; Dexter & Fragile 2013; Chan et al. 2009; Dodds-Eden et al. 2010; Ponti et al. 2017) leading to local acceleration of electrons to relativistic $\gamma$-factors of $10^{3 \dots 6}$, sufficient to generate the variable IR (and X-ray) emission, in analogy to solar flares (Dodds-Eden et al. 2010; Lin et al. 2001; Ponti et al. 2017). This is the subject of the current paper.

## 2. Observations

The GRAVITY instrument combines the four 8 m telescopes of the European Southern Observatory (ESO)-Very Large Telescope (VLT) interferometrically for 3 milli-arcsec (mas) resolution imaging and ~20–70 microarcsec ($\mu$as) astrometry in the K-band (2.2 $\mu$m) continuum. For details of the instrument and the data analysis and positional extraction we refer to Gravity Collaboration et al. (2017) and Appendix A. Briefly, the light of the four telescopes is extracted into mono-mode fibres for two positions on the sky and then interfered in the beam combiner for all six baselines of the interferometer. One fibre is placed on the bright ($K_s = 10$) star IRS16C about 1" N-E
of SgrA*, and the other is on SgrA*, plus the orbiting star S2. In 2018 these two sources were conveniently separated by only 14–20 mas, which is less than the fibre diameter (50 mas) and the diffraction beam of an 8 m UT at 2.2 µm (56 mas). This allowed precise, continuous measurements of the positional separation vector between the K-band continuum emission of SgrA* and S2, $\delta r(t)_{SgrA*-S2}$, while the fringes were detected and stabilised with the second fibre on IRS16C. Since IR flares are polarized, we recorded the Stokes component $Q$, or $Q$ and $U$.

The astrometric position precision of GRAVITY measurements in good atmospheric conditions depends on the magnitude of SgrA* (relative to S2, which has $K_s = 14$) and the integration time. While SgrA* is detected in ~90% of our science frames in 2018, its median magnitude is ~17, which requires a co-addition of more than an hour to reach an rms astrometric precision of $\sigma_{astromet} \sim 30$ µarcsec (µas). However, during bright “states” (henceforth referred to as “flares”, with $K_s \le 15$, whose probability of occurrence is $<10^{-2}$, Dodds-Eden et al. 2011; Witzel et al. 2018) this positional precision can be reached in integrations of 2–10 min. Three additional corrections to the raw astrometry are required. The orbital motion of the reference (S2) over an hour observation is about 10 µas, which is subtracted from $\delta r(t)_{SgrA*-S2}$. Since the near-IR spectral energy distribution (SED) of SgrA* is significantly redder than that of S2, its flux weighted wavelength across the instrumental bandpass is larger than that of S2, and $\delta r(t)_{SgrA*-S2}$ has to be corrected by 1.001. Likewise, differential refraction results in corrections of 10–20 µas. With these aspects taken into account, it is possible to detect the expected motions of a compact, hot spot orbiting a black hole of 4 million solar masses near the ISCO, which for a non-rotating black hole corresponds to an orbital diameter of 60 µas and an orbital period of 31 min.

# 3. Results

## 3.1. Emission centroid motions

We observed two bright flares with a peak approaching the flux of S2 on July 22 and July 28, 2018, as well as a fainter flare (0.3 – 0.5× S(S2)) on May 27, 2018. These flares lasted for 30–90 min. Figure 1a shows for the July 22 flare the RA (blue) and Dec (red) positional offsets (from their medians) of the SgrA* emission centroid as a function of time, along with the flux evolution in units of S2’s flux ($K_s = 14$, or 15 mJy, black). At all times, K-band emission from SgrA* is spatially unresolved. We detect significant and continuous positional changes of the emission centroid in both coordinates, of ~120 µas over ~30 min. This corresponds to ~0.3 times the speed of light. The motions appear to trace out 50–70% of a closed, clockwise loop (Fig. 1b,d). In RA/Dec the data can be fitted well by sine-curves of the same amplitude, characteristic for a closed circular orbit observed near face-on (Fig. 1c). Four independent analyses of these data with different software by different GRAVITY team members give comparable and consistent results on these motions. Figure 1 uses the average of the “Pfuhl (P)” and “Waisberg (W)” analyses, which are also shown separately in Fig. B.1 (Appendix B). The position of SgrA*’s mass centroid derived from the S2 orbit (Gravity Collaboration et al. 2018) is consistent with the centroid of the flare orbit, to within the ±50 µas uncertainties (orange square and cross in Figs. 1b and d).

We discuss in Appendix B measurements of a second, similarly bright flare on July 28, 2018, (MJD = 58328.0841) and a fainter flare on May 27, 2018 (MJD = 58266.3420). The data in Fig. B.4 broadly exhibit the same properties as those in

<table><caption>**Table 1.** GYOTO fit parameters for hot spot orbits of the three flares.</caption><thead><tr><th>Flare</th><th>Analysis</th><th>a</th><th>R (R<sub>g</sub>)</th><th>i (°)</th><th>Ω (°)</th><th>χ<sub>r</sub><sup>2</sup> (N<sub>d.o.f.</sub>)</th></tr></thead><tbody><tr><td rowspan="6">Jul 22nd</td><td>P</td><td>0</td><td>7.3</td><td>164</td><td>118</td><td>1.1 (18)</td></tr><tr><td>W</td><td>0</td><td>7.2</td><td>151</td><td>122</td><td>1.9 (14)</td></tr><tr><td>P+W</td><td>0</td><td>8.0</td><td>149</td><td>115</td><td>1.3 (14)</td></tr><tr><td>P</td><td>-1</td><td>7.7</td><td>151</td><td>109</td><td>1.07 (18)</td></tr><tr><td>L</td><td>0</td><td>7.1</td><td>167</td><td>185</td><td>1.80 (16)</td></tr><tr><td>L</td><td>-1</td><td>7.4</td><td>157</td><td>164</td><td>1.65 (16)</td></tr><tr><td rowspan="2">Jul 28</td><td>P</td><td>0</td><td>9.1</td><td>164</td><td>103</td><td>4.5 (30)</td></tr><tr><td>P</td><td>-1</td><td>9.1</td><td>152</td><td>110</td><td>4.3 (30)</td></tr><tr><td rowspan="2">May 27</td><td>P</td><td>0</td><td>8.3</td><td>179</td><td>131</td><td>2.7 (16)</td></tr><tr><td>P</td><td>-1</td><td>7.1</td><td>163</td><td>127</td><td>2.1 (16)</td></tr></tbody></table>

**Notes.** Analysis “Pfuhl” (P), “Waisberg” (W). “LESIA” (L) and average of P and W (P+W).

<table><caption>**Table 2.** NERO fit parameters for hot spot orbits of the three flares.</caption><thead><tr><th>Flare</th><th>Analysis</th><th>a</th><th>R (R<sub>g</sub>)</th><th>i (°)</th><th>Ω (°)</th><th>χ<sub>r</sub><sup>2</sup> (N<sub>d.o.f.</sub>)</th></tr></thead><tbody><tr><td rowspan="4">Jul 22nd</td><td>P</td><td>0</td><td>8</td><td>145</td><td>126</td><td>1.6 (16)</td></tr><tr><td>W</td><td>0</td><td>7</td><td>140</td><td>160</td><td>1.6 (16)</td></tr><tr><td>P+W</td><td>0</td><td>7</td><td>160</td><td>160</td><td>1.2 (16)</td></tr><tr><td>L</td><td>0</td><td>7</td><td>160</td><td>160</td><td>1.6 (16)</td></tr><tr><td>Jul 28th</td><td>P</td><td>0</td><td>9</td><td>135</td><td>137</td><td>4.5 (30)</td></tr><tr><td>May 27th</td><td>P</td><td>0</td><td>7</td><td>170</td><td>34</td><td>2.1 (16)</td></tr></tbody></table>

**Notes.** For description see caption of Table 1.

Fig. 1. Again the emissions for these two flares appear to follow incomplete clockwise loops of total maximum to minimum amplitude ~110–140 µas over a time period of ~40–70 min, and again the loop centroids are consistent within the uncertainties with SgrA*’s mass position. Overall the data quality is somewhat poorer than in Fig. 1, owing to the less favourable atmospheric conditions and lower fluxes. The short lifetimes observed for all three best flares are expected in models of hot spots (Broderick & Loeb 2005, 2006; Hamaus et al. 2009) since the differential rotation in an accretion disk would in a single orbital period shear out an initially compact gas cloud into an elongated arc.

The apparent orbits of the flares, in addition to their common clockwise motion, appear to be relatively face on (Tables 1 and 2). Our best fits of the data in Fig. 1 with the NERO and GYOTO relativistic ray tracing codes (Bauböck et al., in prep., Vincent et al. 2011b) including light bending, lensing, time dilation, etc.) yield $R \approx 7 \pm 0.5$ R<sub>g</sub>, a line of nodes at $\Omega \approx 115^{\circ}-160^{\circ}$ and an inclination $i \approx 160^{\circ} \pm 10^{\circ}$, with an orbital period of $P = 40 \pm 8$ min, for orbits in an angular momentum parameter $a = 0$ (Schwarzschild) space time (Tables 1 and 2). See also Appendix B, Fig. B.2, and Table 2. Tables 1 and 2 summarise the fit results for all three flares obtained with two independent analyses of the data and two independent relativistic fit codes. All three flares can in principle be accounted for by a single orbit model. Interestingly this common orbit shows a similar orientation and angular momentum direction as the clockwise stellar disk and the G2 object ($\Omega \approx 99^{\circ}$, $i \approx 129^{\circ}$, Bartko et al. 2009; Genzel et al. 2010; Gillessen et al. 2012; Pfuhl et al. 2015; Plewa et al. 2017). Early estimates of the orientation of the hot gas in the innermost accretion zone based on radio data suggest a similar orientation as well ($\Omega \sim 128...157^{\circ}$ E. of N. (note 180^{\circ}$
<figure><img src="image_3.png" /><figcaption>**Fig. 2.** Bottom left (panel a): total flux I/S2 (blue, relative to S2) and Q/I flux ratio as a function of time for the May 27 flare. Top left (panel b): as in (a) but for the July 22 flare. Top right (panel c): I/S2 (blue), $Q/(Q^2 + U^2)^{1/2}$ (red) and $U/(Q^2 + U^2)^{1/2}$ (green) evolution during the July 28 flare. Bottom right (panel d): evolution of the July 28 flare in the plane of normalized Stokes parameters $Q/(Q^2 + U^2)^{1/2}$ (horizontal) and $U/(Q^2 + U^2)^{1/2}$ (vertical). The red arrows denote the polarization directions on sky. The blue cross in the lower right denotes a typical error bar.</figcaption></figure>

degeneracy), yet higher inclination $i \sim 50...68^\circ$, Dexter et al. 2010; Broderick et al. 2011; Wang et al. 2013).

Figures 1 and B.4 show that flux distribution in all three flares is not strongly peaked on any one side of the line of nodes, as would be expected from relativistic beaming for an inclined orbit. We show in Appendix C that inclinations of >155° (or equivalently <25°) would reduce the expected impact of relativistic beaming of radiation to below a factor of 1.5 at one to two times ISCO, which is consistent with our data and the relativistic orbit modelling (Tables 1 and 2), since we see little evidence for rapid or even superluminal centroid motions due to strong light bending or multiple images (Fig. B.2, Broderick & Loeb 2005, 2006; Hamaus et al. 2009). In cases of such low inclinations, the observed flux variations do not primarily reflect variations in beaming and lensing during an orbit but are more determined by the balance of local heating and cooling of the relativistic electrons. The synchrotron cooling time is $\tau_{synchr} = 15 \times (B/20\text{G})^{-1.5} (d/2.2\mu\text{m})^{0.5}$ min (Gillessen et al. 2006). Theoretical estimates of the magnetic field in the innermost accretion zone around SgrA* vary between 5 and 100 G (Yuan et al. 2003; Mościbrodzka et al. 2009; Dodds-Eden et al. 2010; Dexter et al. 2010; Ponti et al. 2017). If the hot spot is created by magnetic reconnection, the local magnetic energy density is tapped for heating of the electrons, and as a result $B^2/8\pi$ drops considerably, lengthening the cooling time

(Dodds-Eden et al. 2010; Ponti et al. 2017). It is thus uncertain whether the observed flare durations are determined by the cooling time and/or the shearing of the hot spot due to differential rotation.

## 3.2. Rotation of polarization angles

In addition to the astrometric motions we observe systematic and continuous variations of the position angle of the polarization of the IR synchrotron emission. Figure 2a shows a full-turn change of Q/I projection over $P_{pol} \approx 73 \pm 15$ min during the May 27 flare, and Fig. 2b indicates half a turn of polarization change in $26 \pm 8$ min during the July 22 flare (in the latter two cases we only recorded one polarization). Figures 2c,d shows the evolution of the $U/(Q^2 + U^2)^{1/2}$ ($\pm45^\circ$) versus $Q/(Q^2 + U^2)^{1/2}$ (0/90°) normalized Stokes components during the July 28 flare. The data appear to trace out a loop structure in time, with a period of $P_{pol} \approx 48 \pm 6$ min, again suggestive of orbital motion of an ordered field/polarization structure (see Appendix D). We note that significant swings of polarization angle have previously been observed in NACO polarimetry of a few SgrA* flares (Trippe et al. 2007; Zamaninasab et al. 2010; Shahzamanian & Eckart 2015), albeit not as impressive and complete as in Fig. 2. We show in Appendix D that such smooth polarization swings (with $P_{pol} \sim P_{orbit}$) can be accounted for if the magnetic field axis is orthogonal to the orbital axis,
<figure><img src="image_4.png" /></figure>

**Fig. 3.** Radius (R) – orbital period (P<sub>orb</sub>) relations for flare loops on July 22 (black), July 28 (blue), and May 27 (pink) compared to circular orbits around a black hole of 4.14 (solid lines), 3, and 5.5 (dashed grey lines) million solar masses. The solid red curve is for orbits around a spin 0 (Schwarzschild) black hole (Schwarzschild 1916; Bardeen et al. 1972). The solid green and black lines are for prograde and retrograde orbits around a maximum (spin parameter |a| = 1) Kerr black hole (Kerr 1963). Small green, red, and black circles mark the innermost stable circular orbits (ISCO) in these cases (Bardeen et al. 1972). For the three flare loops, the filled black circle, blue square, and open pink triangles mark the (R, P) estimates from the best Schwarzschild (a = 0) orbit fits with the relativistic ray tracing codes (Figs. 1, B.3, B.4, Tables 1 and 2). In principle all three flares can be accounted for by the same Schwarzschild R ≈ 6–10 R<sub>g</sub> circular orbit with a period of ~45 min. Interestingly the data can also be fit by a retrograde orbit on ISCO for an extreme Kerr hole (a = −1). The periods P<sub>pol</sub> for the full swings of the polarization orientation of the July 22, July 28 and May 27 flares are given as right pointing black, blue, and pink arrows (with their rms uncertainties marked as dotted vertical lines; see Fig. 2). For a polarized hot spot in a poloidal field the polarization rotation periods are equal to the spot orbital period, in agreement with our data.

as in a poloidal field configuration, and if the orbiting hot spot
is observed at low inclination. The *Q* − *U* loops are caused by
light-bending effects. If the field instead were dominated by a
toroidal configuration in the plane of the motions, one would
expect *P*<sub>pol</sub> = 0.5 × *P*<sub>orbit</sub>.

4. Discussion – Evidence for orbital motion in the deep relativistic zone

Figure 3 summarises our constraints on orbit radii and periods of the flares, interpreted as orbiting compact hot spots, and on the expected orbit kinematics around a massive black hole. As discussed above, in this simple hot spot scenario, the centroid motions as well as the polarization rotation in all three observed flares (July 22, July 28, May 27, 2018) are broadly consistent with similar circular orbits of a polarized hot spot. To obtain more quantitative estimates, the impact of orbit broadening due to astrometric errors, due to the incomplete orbital coverage and the effects of GR and SR need to be included. When these effects are included in circular orbit models (Figs. 1, B.3, B.4), our analysis shows that the flare centroid motions and polarization swings of all three flares are plausibly consistent with the same circular orbit of R = 6–10 R<sub>g</sub> and P = 33–65 min (Tables 1 and 2). This common orbit is at 1.17(±0.3) ISCO for a low

spin (a ≈ 0, “Schwarzschild”) black hole (for 4 million solar
masses, Schwarzschild 1916; Bardeen et al. 1972). These data
constrain the mass of the Schwarzschild hole to be the same
to within ±30% as that obtained from the precision S2 orbit
(Gravity Collaboration et al. 2018). Alternatively, the two flares
could also be on the retrograde ISCO of a highly spun-up Kerr
hole (a ≈ −1, Fig. B.3, Kerr 1963; Bardeen et al. 1972).

The Bardeen & Petterson (1975) effect is not expected to operate in a hot, geometrically thick accretion flow. So the face-on geometry we infer may be unrelated to the black hole spin direction, unless the accretion flow is brought into alignment for example by magnetic torques (McKinney et al. 2012; Sorathia et al. 2013).

Our analysis and conclusions are fairly simple and empiri-
cally driven. The errors are substantial and the modelling by a
single compact spot might be naïve. If continuous energy injec-
tion dominates, the location of the hot spot may be affected by
the propagation of that injection. It is certainly possible that in
addition to rotation in a disk or torus, the hot spot may also
have a line-of sight motion, perhaps due to a combined disk-jet
geometry (Markoff et al. 2001; Mościbrodzka et al. 2014). How-
ever, the detection of orbital motion in bright flares is poten-
tially evidence against radial motion in a jet being the dominant
component (upper limit on motion perpendicular to orbit <0.1c;
<30 µas); it also argues against Rossby wave instabilities with
m ≥ 2 in the accretion flow (Falanga et al. 2007) as being the
origin of the observed flux modulation. In summary, precision
astrometry with GRAVITY at the VLTI provides strong support
from spatially resolved measurements of the motions and polar-
ization rotation of three strong near-IR flares in the relativistic
accretion zone that SgrA* is indeed a massive black hole in the
Schwarzschild–Kerr metric.

**Acknowledgements.** We are very grateful to our funding agencies (MPG, ERC, CNRS, DFG, BMBF, Paris Observatory, Observatoire des Sciences de l'Univers de Grenoble, and the Fundação para a Ciência e Tecnologia), to ESO and the ESO/Paranal staff, and to the many scientific and technical staff members in our institutions who helped to make NACO, SINFONI, and GRAVITY a reality. S.G., P.P. and C.D. acknowledge support from ERC starting grant No. 306311. F.E. and O.P. acknowledge support from ERC synergy grant No. 610058 (Black-HoleCam). J.D., M.B., and A.J.-R. were supported by a Sofja Kovalevskaja award from the Alexander von Humboldt foundation. A.A. & P.G. acknowledge support from FCT-Portugal with reference UID/FIS/00099/2013.

**References**

Baganoff, F. K., Bautz, M. W., Brandt, W. N., et al. 2001, Nature, 413, 45

Balbus, S. A., & Hawley, J. F. 1991, ApJ, 376, 214

Bardeen, J. M., & Petterson, J. A. 1975, ApJ, 195, L65

Bardeen, J. M., Press, W. H., & Teukolsky, S. A. 1972, ApJ, 178, 347

Bartko, H., Martins, F., Fritz, T. K., et al. 2009, ApJ, 697, 1741

Broderick, A. E., & Loeb, A. 2005, MNRAS, 363, 353

Broderick, A. E., & Loeb, A. 2006, MNRAS, 367, 905

Broderick, A. E., Fish, V. L., Doeleman, S. S., & Loeb, A. 2011, ApJ, 735, 110

Bromley, B. C., Melia, F., & Liu, S. 2001, ApJ, 555, L83

Buchholz, R. M., Witzel, G., Schödel, R., & Eckart, A. 2013, A&A, 557, A82

Chan, C.-K., Liu, S., Fryer, C. L., et al. 2009, ApJ, 701, 521

Connors, P. A., Piran, T., & Stark, R. F. 1980, ApJ, 235, 224

Dexter, J. 2016, MNRAS, 462, 115

Dexter, J., & Agol, E. 2009, ApJ, 696, 1616

Dexter, J., & Fragile, P. C. 2013, MNRAS, 432, 2252

Dexter, J., Agol, E., Fragile, P. C., & McKinney, J. C. 2010, ApJ, 717, 1092

Do, T., Ghez, A. M., Morris, M. R., et al. 2009, ApJ, 691, 1021

Dodds-Eden, K., Porquet, D., Trap, G., et al. 2009, ApJ, 698, 676

Dodds-Eden, K., Sharma, P., Quataert, E., et al. 2010, ApJ, 725, 450

Dodds-Eden, K., Gillessen, S., Fritz, T. K., et al. 2011, ApJ, 728, 37

Doeleman, S. S., Weintroub, J., Rogers, A. E. E., et al. 2008, Nature, 455, 78

Eatough, R. P., Falcke, H., Karuppusamy, R., et al. 2013, Nature, 501, 391

Eckart, A., Baganoff, F. K., Zamaninasab, M., et al. 2008, A&A, 479, 625
Eisenhauer, F., Perrin, G., & Rabien, S. 2008, in *The Power of Optical/IR Interferometry: Recent Scientific Results and 2nd Generation*, eds. A. Richichi, F. Delplancke, F. Paresce, & A. Chelli, 431

Eisenhauer, F., Perrin, G., Brandner, W., et al. 2011, *The Messenger*, 143, 16

Falanga, M., Melia, F., Tagger, M., Goldwurm, A., & Bélanger, G. 2007, *ApJ*, 662, L15

Fritz, T. K., Gillessen, S., Dodds-Eden, K., et al. 2011, *ApJ*, 737, 73

Genzel, R., Eisenhauer, F., & Gillessen, S. 2010, *Rev. Mod. Phys.*, 82, 3121

Ghez, A. M., Salim, S., Weinberg, N. N., et al. 2008, *ApJ*, 689, 1044

Gillessen, S., Eisenhauer, F., Quataert, E., et al. 2006, *ApJ*, 640, L163

Gillessen, S., Genzel, R., Fritz, T. K., et al. 2012, *Nature*, 481, 51

Gravity Collaboration (Abuter, R., et al.) 2017, *A&A*, 602, A94

Gravity Collaboration (Abuter, R., et al.) 2018, *A&A*, 615, L15

Grould, M., Paumard, T., & Perrin, G. 2016, *A&A*, 591, A116

Hamaus, N., Paumard, T., Müller, T., et al. 2009, *ApJ*, 692, 902

Johnson, M. D., Narayan, R., Psaltis, D., et al. 2018, *ApJ*, 865, 104

Kerr, R. P. 1963, *Phys. Rev. Lett.*, 11, 237

Lachaume, R., & Berger, J.-P. 2013, *MNRAS*, 435, 2501

Lapeyrere, V., Kervella, P., Lacour, S., et al. 2014, in *Optical and Infrared Interferometry IV*, Proc. SPIE, 9146, 91462D

Lin, R. P., Feffer, P. T., & Schwartz, R. A. 2001, *ApJ*, 557, L125

Markoff, S., Falcke, H., Yuan, F., & Biermann, P. L. 2001, *A&A*, 379, L13

Mathar, R. J. 2005, *Baltic Astron.*, 14, 277

McCrea, W. H. 1972, *MNRAS*, 157, 359

McKinney, J. C., Tchekhovskoy, A., & Blandford, R. D. 2012, *MNRAS*, 423, 3083

Mościbrodzka, M., Gammie, C. F., Dolence, J. C., Shiokawa, H., & Leung, P. K. 2009, *ApJ*, 706, 497

Mościbrodzka, M., Falcke, H., Shiokawa, H., & Gammie, C. F. 2014, *A&A*, 570, A7

Narayan, R., Igumenshchev, I. V., & Abramowicz, M. A. 2003, *PASJ*, 55, L69

Ott, T., Eckart, A., & Genzel, R. 1999, *ApJ*, 523, 248

Paumard, T., Perrin, G., & Eckart, A. 2008, in *The Power of Optical/IR Interferometry: Recent Scientific Results and 2nd Generation*, eds. A. Richichi, F. Delplancke, F. Paresce, & A. Chelli, 313

Pfuhl, O., Gillessen, S., Eisenhauer, F., et al. 2015, *ApJ*, 798, 111

Plewa, P. M., Gillessen, S., Eisenhauer, F., et al. 2015, *MNRAS*, 453, 3234

Plewa, P. M., Gillessen, S., Pfuhl, O., et al. 2017, *ApJ*, 840, 50

Ponti, G., George, E., Scaringi, S., et al. 2017, *MNRAS*, 468, 2447

Reid, M. J., Menten, K. M., Brunthaler, A., et al. 2014, *ApJ*, 783, 130

Schwarzschild, K. 1916, *Seite*, 1916, 189

Shahzamanian, B., Eckart, A., Valencia-S., M., et al. 2015, *A&A*, 576, A20

Shao, M., & Colavita, M. M. 1992, *A&A*, 262, 353

Sorathia, K. A., Krolik, J. H., & Hawley, J. F. 2013, *ApJ*, 777, 21

Tatulli, E., Millour, F., Chelli, A., et al. 2007, *A&A*, 464, 29

Thompson, A. R., Moran, J. M., Swenson, Jr., G. W., et al. 2017, *Interferometry and Synthesis in Radio Astronomy*, 3rd Edition (Springer)

Trippe, S., Paumard, T., Ott, T., et al. 2007, *MNRAS*, 375, 764

Vincent, F. H., Paumard, T., Gourgoulhon, E., & Perrin, G. 2011a, *Class. Quant. Grav.*, 28, 225011

Vincent, F. H., Paumard, T., Perrin, G., et al. 2011b, *MNRAS*, 412, 2653

Vincent, F. H., Paumard, T., Perrin, G., et al. 2014, *MNRAS*, 441, 3477

Wang, Q. D., Nowak, M. A., Markoff, S. B., et al. 2013, *Science*, 341, 981

Witzel, G., Martinez, G., Hora, J., et al. 2018, *ApJ*, 863, 15

Wyant, J. C. 1975, *Appl. Opt.*, 14, 2622

Yang, X.-L., & Wang, J.-C. 2014, *A&A*, 561, A127

Yuan, F., Quataert, E., & Narayan, R. 2003, *ApJ*, 598, 301

Yuan, F., Quataert, E., & Narayan, R. 2004, *ApJ*, 606, 894

Yusef-Zadeh, F., Wardle, M., Heinke, C., et al. 2008, *ApJ*, 682, 361

Zamaninasab, M., Eckart, A., Witzel, G., et al. 2010, *A&A*, 510, A3

¹ Max Planck Institute for Extraterrestrial Physics (MPE), Giessen-bachstr. 1, 85748 Garching, Germany

² LESIA, Observatoire de Paris, Université PSL, CNRS, Sorbonne Université, Univ. Paris Diderot, Sorbonne Paris Cité, 5 place Jules Janssen, 92195 Meudon, France

³ Max Planck Institute for Astronomy, Königstuhl 17, 69117 Heidelberg, Germany

⁴ 1st Institute of Physics, University of Cologne, Zülpicher Straße 77, 50937 Cologne, Germany

⁵ Univ. Grenoble Alpes, CNRS, IPAG, 38000 Grenoble, France

⁶ Universidade de Lisboa – Faculdade de Ciências, Campo Grande, 1749-016, Lisboa, Portugal

⁷ Faculdade de Engenharia, Universidade do Porto, rua Dr. Roberto Frias, 4200-465 Porto, Portugal

⁸ European Southern Observatory, Karl-Schwarzschild-Str. 2, 85748 Garching, Germany

⁹ European Southern Observatory, Casilla 19001, Santiago 19, Chile

¹⁰ Sterrewacht Leiden, Leiden University, Postbus 9513, 2300 RA Leiden, The Netherlands

¹¹ Departments of Physics and Astronomy, Le Conte Hall, University of California, Berkeley, CA 94720, USA

¹² School of Physics and Astronomy, Tel Aviv University, Tel Aviv 69978, Israel

¹³ Max Planck Institute for Radio Astronomy, Auf dem Hügel 69, 53121 Bonn, Germany

¹⁴ CENTRA – Centro de Astrofísica e Gravitação, IST, Universidade de Lisboa, 1049-001 Lisboa, Portugal

¹⁵ Center for Computational Astrophysics, Flatiron Institute, 162 5th Ave., New York, NY 10010, USA
# Appendix A

## A.1. Observations and analysis of GRAVITY data

In this section, we summarise the May and July 2018 interferometric observations in the near-IR (2.2 µm) of SgrA* relevant to this paper.

The observations were taken at the VLTI in Chile using the recently deployed instrument GRAVITY (Gravity Collaboration et al. 2017). The instrument coherently combines the light of the four 8 m Unit Telescopes (UT) of the ESO Paranal site. We chose the most sensitive low resolution mode of GRAVITY. In this mode the science spectrum is dispersed across 14 pixels, with five independent spectral elements (R ~ 20). All data sets were obtained in polarisation split mode, that is, a Wollaston prism is inserted in the optical train and the two linear polarisations are recorded independently.

Each observation followed the same sequence; all four UTs locked their Coudé IR adaptive optics (CIAO) module on the brightest source in the field, the red supergiant IRS7 ($m_K \approx 6.5$, distance from SgrA* ~5.5"). The interferometric observations started with IRS16NW feeding the fringe-tracker operating at a frame rate of 1 kHz and IRS16C feeding the science channel with an integration time (DIT) of 1 s. Those two bright stars ($m_K \approx 10.0-10.5$, separation from SgrA* ~1") were used to find fringes and to zero the optical delay of the science channel. After this initial bright pair, we kept IRS16C as fringe-tracking star and only changed the science target (e.g. S2 and SgrA*) by moving the internal fibre actuators and rotating the field suitably. The science observations were interleaved exposures of NDIT = 30 and DIT = 10 s each. We repeatedly pointed to the object R2, a moderately bright ($m_K \approx 12.1$, separation ~1.5") nearby unresolved giant star, which served as a local calibrator and the S2/SgrA* binary. On average we took 5 × 5 min SgrA* exposures before moving to the calibrator R2 and the science target S2/SgrA*. Every four or five exposures we interleaved a sky exposure. To do so, we offset the star separator field actuator located close to the Coudé focus by several arcseconds, pointing to a location devoid of stars in the Galactic centre.

## A.2. Data reduction

We used the standard GRAVITY pipeline to process the data (Lapeyrere et al. 2014; Gravity Collaboration et al. 2017). Each exposure consists of NDIT science frames, which are averaged after processing. Each individual frame is flat-fielded and corrected for a background bias by subtracting the closest sky exposure, detector noise, and wavelength calibrated on the internal calibration source. The data are then reduced based on a pixel-to-visibility matrix (P2VM; Tatulli et al. 2007), which represents the matrix encoded instrument transfer function, including throughput, coherence, phase-shift, and cross-talk information for each individual pixel. In a second step the science complex visibilities are phase-referenced to the fringe-tracker complex visibilities using the laser metrology and accounting for the fibre dispersion to get phase-referenced visibilities. The observatory transfer function (i.e. coherence loss due to vibrations, uncorrected atmosphere, birefringence, etc.) was calibrated by observing a local calibrator (in our case the nearby giant star R2 at a distance of distance ~1.5").

We investigated the optimum NDIT number of science frames to average for the subsequent model fitting. We found that for this study, NDIT = 3 frames (i.e. 30 s) is a good compromise, which provides sufficient signal-to-noise ratio (S/N) and a reliable noise estimate, while at the same time allows the fast

flux evolution of SgrA* to be resolved over timescales of minutes. Longer averages tend to smear out the signal due to rapid changes in brightness and/or separation.

## A.3. Model fitting

The reported astrometric positions are based on a two-component binary fitting ($S_{SgrA*}$ and $S_{S2}$). We have developed several independent fitting codes, employing least-square minimisation with start parameter variation, Markov-chain-Monte-Carlo (MCMC) optimisation, and a combination of both techniques. The goal was to have independent consistency checks, ensure robustness of the results based on various optimisation techniques, and to quickly explore new model parameters. All codes are based on a binary model, which can be expressed in its most simple form by a complex visibility V as

$$V(u_k, v_k) = \frac{1 + f_k e^{-2\pi i(u_k \cdot \Delta\alpha + v_k \cdot \Delta\delta)}}{1 + f_k}, \qquad (A.1)$$

where $f_k$ is the flux ratio of the two sources in the spectral channel, $\lambda_k$, $u_k$, and $v_k$ are the spatial frequencies ($u_k = u/\lambda_k$ and $v_k = v/\lambda_k$, with u, v being the physical separation of the telescopes in east, north direction) and ($\Delta\alpha$, $\Delta\delta$) denotes the source separation vector (in right ascension and declination). The red colour of SgrA* is usually expressed as a spectral slope $vF_v \sim v^\beta$ or alternatively $F_\lambda \sim \lambda^{-(1+\beta)}$ with typically $\beta \approx 0^{+0.6}_{-2}$ (Genzel et al. 2010; Witzel et al. 2018; Dodds-Eden et al. 2011). The slope of the early-type star S2 in the Rayleigh-Jeans limit is assumed to be $vF_v \sim v^3$ or $F_\lambda \sim \lambda^{-4}$.

Therefore, we account for the wavelength dependent flux ratio by

$$f_k = F_{SgrA*}/F_{S2} \sim \left( \frac{\lambda_k}{2.2 \mu m} \right)^\gamma, \qquad (A.2)$$

where $\gamma$ relates to the intrinsic spectral slope $\beta$ of SgrA* by $\gamma = 3 - \beta$.

The short coherence length of the low-spectral-resolution mode and the comparably large separation of SgrA* and S2 leads to a coherence loss, which is approximated with an optical delay d dependent factor $\Gamma$ (Lachaume & Berger 2013)

$$\Gamma(d_k, R_k) = \mathrm{sinc}\left(\frac{d_k}{2R_k}\right), \qquad (A.3)$$

with $d_k = 2\pi \times (u\Delta\alpha + v\Delta\delta)/\lambda_k$ and the instrument resolution of a spectral channel $k, R_k = \lambda_k/\Delta\lambda_k$.

The angular separation of S2 and SgrA* is comparable to the beam-diameter of the single-mode fibres (~50 mas). Therefore, we have to consider the relative injection of S2 and SgrA* per telescope (alignment errors can lead to a different injection ratio). Based on auxiliary data we can infer the telescope t dependent injection $I_{t,SgrA*}$ and $I_{t,S2}$ based on the fibre- and object separation. We write the injected flux ratio of telescope t and spectral channel k as

$$f_{t,k} = \frac{I_{t,SgrA*}}{I_{t,S2}} f_k, \qquad (A.4)$$

with the telescope t and the intrinsic (wavelength dependent) flux ratio $f_k = f_{2.2}(\lambda_k/2.2 \mu m)^\gamma$ and obtain, in the most basic model, the complex visibility for telescope 1,2 for one spectral channel k,

$$V(u_k, v_k) = \frac{1 + \Gamma(d_k, R_k) \sqrt{f_{1,k}f_{2,k}f_k} e^{-2\pi i/\lambda_k(u\Delta\alpha + v\Delta\delta)}}{\sqrt{1 + f_{1,k}} \sqrt{1 + f_{2,k}}}. \qquad (A.5)$$
<figure><img src="image_5.png" /><figcaption>**Fig. A.1.** Example corner plot from the “Waisberg” MCMC analysis showing the marginalised likelihood distributions over the five parameters of the fit for one single 30 s frame during the July 22 flare: ΔRA and ΔDec positions, SgrA*/S2 flux ratio *f* in the photometric fibre, the SgrA* spectral index *α*, and the unresolved background level *V*₀ in the data.</figcaption></figure>

Some of the four analysis codes integrate over the bandpass and allow for a colour-dependant background. Overall the results are in excellent agreement and are independent of the optimization technique and the detailed implementation. Figure A.1 shows one example MCMC fit of the July 22 “Waisberg” analysis.

## A.4. Linear polarization analysis

During all observations, a Wollaston prism was introduced into the optical beam between the integrated optics chip and the dispersive grism and the detector (Gravity Collaboration et al. 2017). This setup splits the light into two orthogonal polarization components, which are imaged onto the detector. The Wollaston prism is aligned with the axes of the integrated optics chip to avoid polarization crosstalk from the intrinsically birefringent chip. At the entrance of the GRAVITY instrument, a half-wave plate is located in the beam. During observations the half-wave plate co-rotates with the field de-rotator device. This ensures that the incident angle of the internal laser metrology is fixed with respect to the derotator mirrors and thus avoids birefringence-related phase-shifts and corresponding astrometric biases. This however leads to a rotation of the detected polarization components with respect to the sky N-E coordinate system. Therefore the measured polarization angle needs to be corrected for the field rotation at the instrument.

In order to analyse the linear polarization properties of a source on sky with an arbitrary position angle, it is common practice to split the polarization and to probe the polarization in at least two rotated states. In our case we offset the co-rotating half-wave plate by 0° and 45° during science and calibrator exposures. Each 0° and 45° science exposure is calibrated with the corresponding calibrator exposures.

A single exposure provides two orthogonal polarization states. Consequently, the half-wave plate at 0° yields the Stokes parameter $Q = P_0 - P_{90}$ and the 45° state provides the parameter $U = P_{45} - P_{135}$. The total intensity is $I = P_0 + P_{90} = P_{45} + P_{135}$. The degree of polarization can be inferred from $\Pi = \sqrt{U^2 + Q^2/I}$ and the polarization angle is defined as $\theta = 1/2 \arctan(U/Q)$.

We tested and calibrated the linear polarization measurement with GRAVITY on sources with known polarization properties (taken from Ott et al. 1999; Buchholz et al. 2013), namely GCIRS21 ($P \approx 14\%$; $\theta \approx 15°$), IRS16SW ($P \approx 3.1\%$; $\theta \approx 20°$), GCIRS33E ($P \approx 5.7\%$; $\theta \approx 35°$) and IRS1W ($P \approx 1.8\%$; $\theta \approx -37°$). Here the polarization angle is defined in the range [-90°, 90°] with the angle increasing east of north. Based on comparison with the aforementioned reference stars, we estimate the systematic uncertainty of the polarization degree measurement to be $\Delta P \approx 1.5\%$ and the angular uncertainty to be $\Delta \theta \approx 10°$. It should be noted that the polarization measurement of SgrA* is actually a differential measurement, since it is based on the flux ratio between S2 and SgrA*. Instrumental polarization effects, which are common to both sources should cancel to first order.

One caveat of measuring the polarization angle and degree with the technique of rotating a half-wave plate is that the two measurements are done sequentially. This requires interpolating $Q$ and $U$ over the full observing sequence to fill the respective gaps. A typical exposure takes 5 min per state and the overheads to rotate the half-wave plate take another ~1 min. This means that the gaps are roughly 6 min. The polarization measurement is effectively smoothed over short timescales. On timescales of a flare (≥30 min), the time resolution is sufficient to capture the evolution of the polarization.

## A.5. Field rotation

One potential source of systematic error is the field rotation during an observation. In particular polarization can be affected due to changing angles of incidence as the mirrors track the source on sky. In principle this is calibrated by repeated exposures on a calibrator star, yet some signature might still be present in the data. Figure A.2 shows the evolution of the field angle at the instrument (subtracted by the mean angle) during the three flares. The field angle at the instrument changes only slowly during all three flares. In particular, during the flare with polarization information (July 28) the angle changed by less than 14°, and on July 22 by only 3°. On the short timescales of the flares the field rotation is unlikely to impact the polarization and astrometry at a significant level.

## A.6. Astrometric precision and accuracy

We derived the astrometric precision based on the scatter of the individual fit results. We fitted individual 30s exposures (3 frames with 10s DIT). The results were binned over 3 min (i.e. ~6 positions) resulting in the final astrometric positions presented in Figs. 1 and 2. Consequently we estimated the astrometric precision from the scatter within each bin. Figure A.3 shows the precision as a function of SgrA* brightness. For the brightest states, we achieve a 1D precision of 20 µas and a 2D precision of ~30 µas as rms.

The long-term stability of the wavelength calibration of GRAVITY in low spectral resolution has been measured to be better than 0.45 nm rms (over 6 months). This corresponds to $2.0 \times 10^{-4}$ in relative terms; a negligible contributor to the uncertainties.
<figure><img src="image_6.png" /><figcaption>Fig. A.2. Field rotation (mean angle subtracted) at the instrument during the three flares.</figcaption></figure>

<figure><img src="image_7.png" /><figcaption>Fig. A.3. Astrometric precision as a function of the flux density of SgrA*. Each point represents the 2D standard deviation of six individual 30 s position fits within 3 min, as in Figs. 1 and 2, divided by $\sqrt{6}$. Here both polarizations were averaged. The contiguous line indicates a photon-noise-limited phase noise (e.g. Wyant 1975) $\sigma_{\phi} \propto 1/\sqrt{N_{ph}}$, where $N_{ph}$ is the number of photons.</figcaption></figure>

## A.7. Systematic astrometric error sources

### A.7.1. Baseline uncertainty

From Eq. (A.1) it is obvious that the measured separation $\Delta\alpha$, $\Delta\delta$ depends on the precise knowledge of the baseline vector $(u, v)$. Any uncertainty in the baseline ($\delta B$) directly translates into a separation error $\delta S$ as

$$ \frac{\delta B}{B} = \frac{\delta S}{\Delta S} \qquad (A.6) $$

The knowledge of the baseline vector is limited by several terms such as the global uncertainty of the telescope array geometry and the knowledge on the actual pointing vector to the object (i.e. an uncertainty in RA and Dec). The global array geometry only refers to the separation of the pivot points of the telescopes. Unlike in radio interferometry, where the separation of the telescopes is much larger than the size of the telescopes and the telescope can be considered point-like in the baseline $u-v$ space, in optical interferometry the telescopes size can play a significant role. This means that optical errors such as pupil mis-registration or pupil vignetting can contribute to baseline errors.

The typical optical path accuracy of the VLTI is <1 cm (i.e. the fringes are found within <1 cm from the predicted position) during a preset to objects, where accurate positions are

available. This means that the global array geometry has to be accurate to the same level. The absolute positions of SgrA* is extremely well determined based on a radio reference frame, which is tied to nearby Quasars (accuracy <1 mas, Reid et al. 2014; Plewa et al. 2015). This leads to a negligible baseline error of 1 mas × 100 m ≈ 10⁻⁴ cm. In order to minimize the pupil error, the GRAVITY instrument is equipped with a dedicated laser guiding system, which controls pupil runout. The guiding system centres the pupil with a residual 1D scatter over 5 min of ~0.5 cm rms. Constant monitoring of the pupil illumination argues against baseline errors due to vignetting. Overall the baseline error is at maximum on the order of ~1 cm. For a 100 m baseline the corresponding systematic astrometric error for two sources separated by ~20 mas is $\delta S \approx 20 \text{ mas} \cdot 10^{-2}/100 \approx 2 \mu\text{as}$. The small separation between SgrA* and S2 in 2018 leads to a negligible contribution of baseline errors in the relative astrometry.

### A.7.2. Effective wavelength accuracy

Another systematic uncertainty is related to the effective wavelength (see Eq. (A.1)). The wavelength determines the image scale for an interferometer. This means that any wavelength uncertainty leads to a proportional astrometric uncertainty. In the low-resolution mode used for the observations the instrument features 14 spectral channels across the K-band (2.0–2.45 μm). The effective bandwidth of each channel varies from FWHM ~70–140 nm, with significant overlap between the channels. The large bandwidth of the spectral channels can lead to a shift of the effective wavelength for objects, which are significantly redder than the relatively blue calibration lamp. The average $H-K$ colour in the Galactic centre is 1.8 mag (Fritz et al. 2011). This steep colour slope from extinction can be approximated by a power-law dependence of $E(\lambda)_{2,2} \propto \lambda^{-5}$ in the K-band. The strong extinction leads to an effective colour temperature of a B0V type star such as S2 of only ~1000 K. In comparison the calibration lamp with a colour temperature of ~2800 K appears blue. Considering the exact bandpass shape, we calculated an overall effective wavelength shift of ~2 nm, that is, an image scale change of ≈2 nm/2200 nm = 0.1%. This overall shift is taken into account for the final astrometry. However, the spectral slope of SgrA* is intrinsically redder than S2 and potentially varies with the state of SgrA* between $\beta = -2... + 0.6$ (Genzel et al. 2010). This leads to an additional effective wavelength uncertainty, which is of the order ~1 nm, that is, 0.05%. The corresponding astrometric uncertainty is ~10 μarcsec, assuming 20 mas separation between S2 and SgrA*.

### A.7.3. Atmospheric limitations

Due to a small non-common atmospheric path, the two objects are subject to differential tip-tilt. The order of magnitude can be approximated from (Shao & Colavita 1992),

$$ \sigma_{tt} = 540 B^{-2/3} \Theta \cdot t^{-1/2} \text{ [arcsec]} \qquad (A.7) $$

with $B = 100$ m, object separation $\Theta = 1 \times 10^{-7}$ rad (20 mas) and $t = 100$ s, the residual tip-tilt is insignificant at the level $\sigma_{tt} \approx 0.2 \mu\text{arcsec}$.

### A.7.4. Atmospheric refraction

Atmospheric refraction leads to a wavelength-dependent shift of the observed zenith angle relative to the true zenith angle. For a single telescope this shift needs to be considered for an
<figure><img src="image_8.png" /><figcaption>Fig. B.1. Comparison of the Pfuhl (bottom) and Waisberg (top) analyses of the July 22 flare data. Nomenclature as in Fig. 1.</figcaption></figure>

accurate pointing and astrometry. The atmospheric refraction is defined as $R := z_t - z_a$ with the true and apparent zenith distance $z_t$ and $z_a$. The wavelength-dependent refraction causes astronomical objects to appear dispersed into a spectrum along the parallactic angle. The refractive index of air at 2.2 µm for typical Paranal conditions (pressure = 744 hPa, temperature = 10°C, humidity=10%) is $(n_{2.2} - 1) \times 10^6 = 203.95$ and the dispersion is $dn/d\lambda|_{2.2} = -2.21 \times 10^{-10} [\text{nm}^{-1}]$. The corresponding differential refraction across the K-band (2.0–2.45 µm) leads to an angular dispersion of $\Delta R \approx 20 [\text{mas}] \tan(z_t)$. This can lead to a loss of injection into the single-mode fibres for large zenith angles.

For interferometers such as the VLTI the atmosphere can be considered uniform and plane parallel (earth curvature can be neglected, e.g. Mathar 2005). The refraction as the rays travel through the atmosphere is governed by Snell's law, that is, $n_0 \sin(z_0) = n \sin(z)$ is constant. The corresponding delay measured in an interferometer $D = b n_0 \sin(z_0) = b n \sin(z)$ is constant and is equal to the delay in free space (Thompson et al. 2017). In the case of an interferometer with evacuated delay lines, atmospheric refraction only needs to be considered for the pointing of the telescopes but has no effect on the interferometric measurement.

The delay lines of the VLTI are however in air and therefore dispersion is introduced due to the optical path compensation. The effect of dispersion in the delay lines is exactly the same as differential refraction on sky.

The colour difference between S2 and SgrA* leads to a differential refraction between the two objects for observations at non-zero zenith angles. Across the K-band the differential refraction in the direction of the parallactic angle (with stronger refraction in the blue) is

$$ \Delta R = 45 (\mu\text{as nm}^{-1}) \cdot \Delta\lambda \tan(z_t), \qquad (A.8) $$

for a wavelength shift of $\Delta\lambda$. The three flares discussed here were all observed at zenith angles $z_t < 35^\circ$ and changed during the observation by less than 10°. We assume an average colour for SgrA* of $\beta = 0$, which corresponds to an effective wavelength shift of 1 nm relative to S2 for 100 nm the bandwidth of the instrument resolution. We take the corresponding astrometric shift into account. The relative shift during the flares due to changing zenith angle is in all cases smaller than 15 µarcsec. Relevant for the flare astrometry is only a change in differential refraction. The effective wavelength uncertainty of SgrA* of up to ~1 nm (due to its varying spectral slope), translates into a maximum differential refraction variation of ~26 µas.
<figure><img src="image_9.png" /><figcaption>Fig. B.2. Light curves (top) and centroid tracks (bottom) of a compact hot spot orbiting a Schwarzschild black hole for four different orbital radii and inclinations ($R/R_g=6$ (solid blue), 8 (red dashed), 10 (short-dashed green) and 12 (dash-dotted red)). From left to right the inclinations are 20, 40, 70 and 90°. Since the ordinary centroid for the 90° case simply yields a horizontal line, we instead plot the orbital phase $\Phi(t)$ against $x(t)$ (based on the NERO code (Bauböck et al., in prep.); see also Hamaus et al. 2009).</figcaption></figure>

# Appendix B

## B.1. Data modelling

### B.1.1. Model limitations

We developed four independent codes to fit the coherent signal of SgrA* and S2 and at the same time account for instrumental properties. The codes differ in the fitting approach (MCMC, least-square with parameter grids, etc.), the number of free parameters (e.g. intrinsic source spectral slope is fixed or fitted), and the weighting of closure phase to visibility data. The different choices in the codes are all well motivated and can thus be considered as an exploration of the parameter space. Figure B.1 shows the results of two codes named “Waisberg” and “Pfuhl” for the July 22 flare (Fig. 1). The “LESIA” results are nearly identical to the Waisberg results. The “Horrobin” code shows differences on a similar scale to Pfuhl and Waisberg. While the four codes agree on the main results and features, they show also differences, which reflect another source of systematic uncertainty on the order 20–30 µarcsec.

## B.2. Astrometric modelling of orbits

Figure B.2 shows models of hot spots on circular orbits in the Schwarzschild metric (angular momentum parameter $a = 0$) obtained with ray tracing methods of the geodesics, including secondary, tertiary, and quaternary images (Hamaus et al. 2009, see also Broderick & Loeb 2005, 2006; Vincent et al. 2011b). The same reference also shows models for the Kerr metric. Similar models with or without polarization have been published elsewhere (Broderick & Loeb 2005, 2006; Hamaus et al. 2009; Vincent et al. 2011b, 2014). The $a = 0$ orbits for inclinations $<20°$ near ISCO and at somewhat higher inclinations for

$R>\text{ISCO}$ result in spot centroid motions that are little influenced by light bending and multiple images. The reverse is true for high spin, highly inclined orbits. Given the smooth near-circular spot motions at least for the July 22 data (Fig. 1), the data are consistent with the low inclination orbits in Fig. B.2.

To simulate the effects of astrometry noise on these simple near-circular orbits, we added a randomly drawn $\delta x$ and $\delta y$ noise component to each of the model data points, with a magnitude comparable to the average empirically determined noise in each of the flares 20–25 µas 1D for July 22, and 45 µas for July 28 and May 27. We assumed that the orbits obey the radius R-period P relations of circular orbits in the Schwarzschild–Kerr metric,

$$P = 2.137(\text{min}) \times \left( \frac{M_{\bullet}}{4.14 \times 10^6 M_{\odot}} \right) \times \left( a + \left( \frac{R}{R_{\text{g}}} \right)^{3/2} \right), \quad (\text{B.1})$$

where $a$ is the spin parameter $\in [-1, 1]$, and $R_g = 6.11 \times 10^{11} [M_{\bullet}/4.14 \times 10^6 M_{\odot}]$ cm = 5 µas. We then varied $R$ to obtain a reasonable match with the data (and their uncertainties) in the three flares. As we discuss in the main text, we find that an orbit on a common radius, $R \approx 1.17 \pm 0.25 R_{\text{ISCO}}(a = 0)$ can match all three flares (Fig. 3). However, due to the fact that only a partial orbit is observed, we cannot rule out non-circular orbits.

We model the astrometry data with fully relativistic models of a compact distribution of gas orbiting at a common velocity in the Kerr spacetime (a “hotspot”, Hamaus et al. 2009; Vincent et al. 2011b, 2014; Connors et al. 1980). When the orbit is in the equatorial plane of the black hole, the motion can be calculated analytically with angular frequency $d\phi/dt = \Omega = (r^{3/2} + a)^{-1}$ in Boyer–Lindquist (BL) coordinates. Ray tracing is used to account for relativistic effects. From an observer’s era photon geodesics are traced backwards towards the black
hole. The hotspot density is assumed to fall off as a Gaussian with distance from its centre, with a characteristic size $R_{spot} \ll R$. The intensity along the ray is then calculated assuming the emission is optically thin and taking into account the frequency shifts due to Doppler beaming and gravitational redshift.

We have used three different codes for this technique. All codes allowed for optimisation of the orbit radius, inclination, position angle, and phase angle of the spot at $t = 0$. The spin parameter was fixed to either 0 or -1, corresponding to counter-rotating around a maximally spinning black hole. In grtrans<fnref n="2" /> (Dexter 2016; Dexter & Agol 2009), the emissivity is polarized synchrotron radiation from a power-law distribution of non-thermal electrons (Broderick & Loeb 2006). The ray tracing technique accounts both for the emitted polarization and its parallel transport to the observer.

In GYOTO<fnref n="3" /> (Vincent et al. 2011b,a; Grould et al. 2016), the source is an optically thick coordinate-sphere of radius 0.5 (in BL coordinates) with emissivity $I_v(v) = \text{constant}$. In this set-up, the power-law index would only matter for weighing the contribution of the secondary image, which is minor at such low inclination. The source is set on a circular orbit in the equatorial plane of a Kerr black-hole and we ray-trace the appearance of the sphere at the actual observing dates over a field of view of 200 µas with a pixel size of 1 µas. The model astrometry is then the centroid of these ray-traced images. We minimise the free parameters with the standard Python procedure scipy.optimize.curve_fit. The distance of the GC and mass of the central object are fixed at $M_{BH} = 4.14 \times 10^6 M_\odot$ and $R_0 = 8.127$ kpc (Gravity Collaboration et al. 2018).

In the NERO code (Bauböck et al., in prep.), we combine the YNOGKM<fnref n="4" /> (Yang & Wang 2014) code to calculate the timelike geodesics of particle orbits with the geokerr<fnref n="5" /> (Dexter & Agol 2009) code to find the null geodesics of photon trajectories. Here we consider only models of a small hotspot with a Gaussian density profile on a circular orbit around a non-spinning ($R \ge 6 R_g$) or maximally spinning ($R < 6 R_g$) black hole. We image the resulting orbit by ray-tracing photon trajectories in a field of view corresponding to 1.5 times the size of the orbital radius. We calculate the flux at each point by integrating the density along the photon path, that is, the emissivity is proportional to the density. By accounting for the time of travel of the photons as well as the hotspot motion, we can calculate a time-dependent image of the region near the black hole. The position-weighted average of the flux of this image gives the centroid of the emission.

We fit the resulting centroid tracks to the data over a parameter space that spans the radius and inclination of the orbits, the position angle of the angular momentum vector, and the phase of the point at $t = 0$. We employ a grid fit of models between $R = 3 R_g$ and $R = 12 R_g$, where those orbits that fall outside theISCO are modelled with a zero-spin black hole, while at $R < 6 R_g$ we use prograde orbits around a maximally spinning black hole.

<figure><img src="image_10.png" /><figcaption>**Fig. B.3.** Combined constraints of the location (astrometry) of the July 22 flare in the $R/R_g$ – inclination plane, from fits with the NERO code (Bauböck et al., in prep. but also Yang & Wang 2014; Dexter & Agol 2009). For $R \ge 6 R_g$ we use $a = 0$, and for $R < 6 R_g$, $a = 1$. The three astrometry contours represent the 1-, 2- and 3-$\sigma$ constraints from the model fitting, with the dark blue being the most favoured. Since we do not have any prior knowledge of the inclination of the orbit, we have weighted the probabilities by $\sin(i)$ to account for the geometric bias. The polarization measurement of July 28 with a period of $P_{pol} = 48 \pm 6$ min yields the right-hatched vertical constraint. The requirement of a Doppler boost (“contrast”) <2.5 is given by the green area, with a lower contrast pushing further down. The presence of a single and full loop in the polarization data yields the left-hatched, mostly left-right constraint. The combined constraints favour $R = 7.6 \pm 0.5 R_g$ and inclination $i \le 30$ degrees (black encircled area).</figcaption></figure>

All codes find the same relativistic effects. Lensing causes the orbit to appear somewhat larger on the sky. At higher inclination, strong Doppler beaming causes large flux modulations between the approaching and receding sides of the orbit with an additional peak from lensing behind the black hole. The spot appears to move faster when approaching than when receding.

## B.3. Hot spot models of the May 27 and July 28 flares

We have analysed the centroid spot motions of the May 27 and July 28 flares with the same techniques as for the July 22 flare. The results are given in Fig. B.4.

## Appendix C: Limit on inclination of orbits based on the lack of Doppler boosting signal

We have discussed in the main text the lack of any obvious pronounced brightening of the flux evolution in one particular section of the three orbits. If a relativistic hot spot is moving with $\beta = v/c$ and $\Gamma = (1 - \beta^2)^{-1/2}$ at an angle $90 - i$ relative to the line of sight, emitting a spectrum with spectral index $\alpha$, $S_v \propto v^{-\alpha}$, then the combination of relativistic aberration, time dilation, and Doppler frequency shifting boosts or de-boosts the observed flux density in the stationary observer’s frame by a factor (McCrea 1972):

$$S_{v,\text{obs}} = v^{-\alpha} \times [\Gamma \times (1 - \beta \cos(90^{\circ} - i))]^{-(3+\alpha)}. \quad (C.1)$$

<footnote n="2">² grtrans is available from: https://github.com/jadexter/grtrans</footnote>

<footnote n="3">³ GYOTO from: https://gyoto.obspm.fr</footnote>

<footnote n="4">⁴ YNOGKM from: http://www1.ynao.ac.cn/{\sim}yangxl/yxl.html</footnote>

<footnote n="5">⁵ geokerr from: https://faculty.washington.edu/agol/geokerr/index.html</footnote>
<figure><img src="image_11.png" /><figcaption>**Fig. B.4.** Bottom left (panel a): time evolution of the east-west (east positive, blue) and north-south (red) position offset of the July 28 (MJD = 58328.0841) centroids from their medians, as well the flux density evolution (right y-axis, black), in units of the flux of S2 (14.0 mag). Error bars are $1\sigma$. For this purpose the total intensity was computed from the sum of the two polarization directions. The data points use the “Pfuhl” analysis (see Appendix B). Bottom right (panel b): projected orbit of the flare centroid on the sky with colours from red to blue marking the time evolution. The orange square indicates the black hole position including long-term astrometric uncertainty. The thin continuous blue curve denotes a simple model of a hot spot on a circular orbit with $R = 9R_g$ ($a = 0$, $M = 4.14 \times 10^6 M_\odot$), seen at inclination $135^\circ$ (clockwise on the sky, as for the data (black curved arrow)) and with the line of nodes at $\Omega = 120^\circ$ (same as in Fig. 1), fitted with NERO code (Bauböck et al., in prep. and Yang & Wang 2014; Dexter & Agol 2009). The cyan and violet curves show the orbit in $x(t)$ and $y(t)$. Top (panels c and d): as in (panel a) and (panel b) but for the May 27 flare (MJD 58266.3420). Here the best fitting orbit is face on with $R = 7R_g$.</figcaption></figure>

Figure C.1 shows these “boosting” factors as a function of angle $90^\circ - i$ for an ISCO ($a = 0$) orbit at $R = 6R_g$ (blue) and an orbit at $R = 10R_g$ (red), each for two spectral indices, $\alpha = 0.6$ (solid) and $\alpha = 1.5$ (dotted). A boosting factor <1.5 plausibly inferred from the light curves in Figs. 1, B.2, and B.3 suggests $i < 27^\circ$. At such low inclinations,

almost face on, one would also expect little lensing and multiple images for $R \ge 6R_g$ (see Broderick & Loeb 2006; Hamaus et al. 2009), again consistent with our observations. We have found a consistent inclination constraint using the flux modulation in light curves from the relativistic hot spot models described above.
<figure><img src="image_12.png" /></figure>

**Fig. C.1.** Doppler boosting factor in the observer's frame *o* relative to the hot spot's moving frame *hs*, $S_o \sim S_{hs} \times D$, with $D = (\Gamma(1-\beta))^{-(3+\alpha)}$, for a moving hot spot at $\beta = v/c$, $\Gamma = (1-\beta^2)^{-1/2}$, and at an angle $90^\circ - i$ relative to the line of sight (McCrea 1972). We show *D* for two radii, ISCO $a=0$ (blue) and $R=10R_g$ (red). Solid curves assume spectral index $\alpha = 0.6$ ($S_v \propto v^{-\alpha}$, appropriate for very bright flares) and dotted curves assume $\alpha = 1.5$ ($S_v \propto v^{-\alpha}$, appropriate for average bright states, Witzel et al. 2018). If one assumes that $D < 2.5$, the inclination has to be less than about $45^\circ$.

## Appendix D: Polarization loops

During the July 28 flare we observed both polarization components with the half wave plate, as discussed in the instrumental description above. Figure 2(d) shows the location of the July 28 flare in the *Q* (horizontal) and *U* (vertical) plane, normalized to the total intensity and with a typical rms error shown in the lower left. The fractional polarization is 20–40%, consistent with past single telescope results (Trippe et al. 2007; Eckart et al. 2008; Shahzamanian & Eckart 2015). During the 100 min duration of the flare the *Q*-*U* Stokes components trace out about 1.3 times a full “loop” or 180° rotation of the polarization angle. We estimate the period for the polarization loop to be $P_{pol} \sim 48 \pm 6$ min, comparable to the orbital period of the centroid motion in that

flare (Fig. 3). For the May 27 and July 22 flares we only recorded the *Q*-polarization. Figure 2a shows that *Q* varied over a full period during the ~70 min duration of the May 27 flare, again comparable to the orbital period, while the polarization on July 22 rises from zero to a peak and goes back to zero, consistent with half of a period.

The high Sgr A* near-IR polarization arises from optically thin synchrotron emission from energetic electrons (energy $\gamma = E/(mc^2) \approx 10^3$). In the hotspot model described above, the polarization angle as a function of time traces out the magnetic field geometry in the accreting gas around the black hole. If the field is initially weak, it will be sheared into a toroidal configuration (Balbus & Hawley 1991) while any accreted magnetic flux will form a magnetically dominated atmosphere near the pole. If sufficient flux can be advected inwards, the fields can become dynamically important with predominantly vertical magnetic field at all latitudes near the black hole (Narayan et al. 2003). Figure D.1 shows models of the motion of the *Q*-*U* Stokes parameters at three different radii on a circular orbit around a Schwarzschild black hole, where the dominant magnetic field direction was poloidal, or perpendicular to the orbital plane. In this case the polarization angle varies in loops with a period $P_{pol} = P_{orbit}$, as in the data of Fig. 3. The loops are more pronounced and centered on $(Q,U) = (0,0)$ for smaller orbital radii. In contrast, for a toroidal field configuration, one would expect loops with half the orbital period, $P_{pol} = 0.5 \times P_{orbit}$, or two loops per orbital period (Bromley et al. 2001; Dexter 2016). For the range of inferred radii from the astrometry data, producing a loop passing through (0,0) as observed requires an inclination of $i \sim 15-30^\circ$. This low inclination is consistent with both the apparent orbit on the sky and the lack of strong flux modulation from Doppler beaming or gravitational lensing.

The loops in the poloidal field case are caused by light bending, which adds an inwards radial component to the wave vector and an azimuthal component to the polarization map. At low inclination and small radius this effect becomes strong enough to produce one large loop per orbit. The observed polarization signature is therefore further evidence that the material is orbiting close to the black hole. Our finding of a likely poloidal field geometry further suggests that the SgrA* accretion zone may be magnetically dominated, consistent with the high field strength inferred at ~0.1 pc from the rotation measure of the magnetar SGR J1745-2900 (Eatough et al. 2013).

Alternatively, such a field configuration could arise in a jet. If part of a collimated outflow, the emission radius would need to be substantially larger than what we observe and be only mildly relativistic to obey our observed limits on Doppler beaming and apparent out-of-plane motion.
<figure><img src="image_13.png" /><figcaption>**Fig. D.1.** Left (panel a): simulation of the motion of a polarized, orbiting hot spot for a poloidal magnetic field configuration, for 30° inclination and three different orbital radii $R = 4, 8, 12 R_g$ (panel b) and $R = 8 R_g$ at varying inclination. The observed July 28 flare polarization data are consistent with the $Q/U$ loop signature and period for the $8 R_g, i = 30°$ curve (but at lower polarization fraction).</figcaption></figure>