Issue 
A&A
Volume 608, December 2017



Article Number  A60  
Number of page(s)  22  
Section  Astrophysical processes  
DOI  https://doi.org/10.1051/00046361/201731148  
Published online  07 December 2017 
General relativistic effects on the orbit of the S2 star with GRAVITY
^{1} LESIA, Observatoire de Paris, PSL Research University, CNRS UMR 8109, Université Pierre et Marie Curie, Université Paris Diderot, 5 place Jules Janssen, 92190 Meudon, France
email: marion.grould@obspm.fr
^{2} LUTh, Observatoire de Paris, PSL Research University, CNRS UMR 8102, Université Paris Diderot, 5 place Jules Janssen, 92190 Meudon Cedex, France
Received: 10 May 2017
Accepted: 22 August 2017
Context. The first observations of the GRAVITY instrument obtained in 2016, have shown that it should become possible to probe the spacetime close to the supermassive black hole Sagittarius A* (Sgr A*) at the Galactic center by using accurate astrometric positions of the S2 star.
Aims. 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 higherorder general relativistic (GR) effects, in particular the LenseThirring effect due to the angular momentum of the black hole.
Methods. We implement seven stellarorbit models to simulate both astrometric and spectroscopic observations of S2 beginning near its next pericenter passage in 2018. Each model takes into account a certain number of relativistic effects. The most accurate one is a fully GR model and is used to generate the mock observations of the star. For each of the six other models, we determine the minimal observation times above which it fails to fit the observations, showing the effects that should be detected. These threshold times are obtained for different astrometric accuracies as well as for different spectroscopic errors.
Results. Transverse Doppler shift and gravitational redshift can be detected within a few months by using S2 observations obtained with pairs of accuracies (σ_{A},σ_{V}) = (10−100 μas, 1−10 km s^{1}) where σ_{A} and σ_{V} are the astrometric and spectroscopic accuracies, respectively. Gravitational lensing can be detected within a few years with (σ_{A},σ_{V}) = (10 μas, 10 km s^{1}). Pericenter advance should be detected within a few years with (σ_{A},σ_{V}) = (10 μas, 1−10 km s^{1}). Cumulative highorder photon curvature contributions, including the Shapiro time delay, affecting spectroscopic measurements can be observed within a few months with (σ_{A},σ_{V}) = (10 μas, 1 km s^{1}). By using a stellarorbit model neglecting relativistic effects on the photon path except the major contribution of gravitational lensing, S2 observations obtained with accuracies (σ_{A},σ_{V}) = (10 μas, 10 km s^{1}), and a black hole angular momentum (a,i′,Ω′) = (0.99,45°,160°), the 1σ error on the spin parameter a is of about 0.4, 0.2, and 0.1 for a total observing run of 16, 30, and 47 yr, respectively. The 1σ errors on the direction of the angular momentum reach σ_{i′} ≈ 25° and σ_{Ω′} ≈ 40° when considering the three orbital periods run. We found that the uncertainties obtained with a less spinning black hole (a = 0.7) are similar to those evaluated with a = 0.99.
Conclusions. 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σ = 0.1) with the S2 star will be possible with a simple stellarorbit model without using a raytracing code but with approximating the gravitational lensing effect. However, long monitorings are necessary, and we thus must rely on the discovery of closerin 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.
Key words: black hole physics / relativistic processes / Galaxy: center / astrometry / infrared: stars
© ESO, 2017
1. Introduction
Decades of studies have demonstrated the presence of a compact object of several million solar masses at the center of the Galaxy (Wollman et al. 1977; Genzel et al. 1996; Eckart & Genzel 1997; Ghez et al. 1998, 2008; Gillessen et al. 2009b, 2017). One of the finest pieces of evidence supporting the existence of this compact source was obtained with the monitoring of S stars in the central parsec over a dozen years by Ghez et al. (2008) and Gillessen et al. (2009b). In particular, a complete orbit of the closest star to the Galactic center, named S2, has been obtained. Such observations combined with monitorings of other S stars led to a confident constraint of the mass of the compact object of ≈(4.31 ± 0.42) × 10^{6}M_{⊙} (Gillessen et al. 2009b; see Boehle et al. 2016; or Gillessen et al. 2017 for a recent improvement on the estimation of this mass). Nowadays, the assumption is that this object is probably a supermassive black hole described by general relativity (GR; see e.g., Broderick et al. 2009, 2011, and references therein). Several methods for proving the existence of this GR black hole are investigated, such as observing the accretion disk around the object, the flares occurring near it, or the stellar orbits of stars close to Sagittarius A* (Sgr A*; Will 2008; Doeleman et al. 2009; Merritt et al. 2010; Grandclément et al. 2014; Broderick et al. 2014; Vincent et al. 2016; Johannsen 2016). The second generation instrument at the VLT (Very Large Telescope), GRAVITY, is expected to better constrain the nature of this object (Eisenhauer et al. 2003). By using its astrometric accuracy of about 10 μas, it will probe spacetime in strong gravitational fields by observing stars and gas located near Sgr A*.
Different theoretical studies have been performed in order to determine whether it will be possible to detect GR effects with stellar orbits. The main purpose is to prove that observations of stars orbiting Sgr A* are affected by GR effects induced by the presence of a Kerr black hole. Several authors have shown that loworder GR effects should be detectable using astrometric and/or spectroscopic measurements, such as the pericenter advance, the transversal Doppler shift or the gravitational redshift (Jaroszynski 1998; Fragile & Mathews 2000; Weinberg et al. 2005; Zucker et al. 2006; Parsa et al. 2017; Nishiyama et al. 2017). In particular, Jaroszynski (1998) showed that it will be easy to reject a stellarorbit model neglecting the pericenter advance when considering astrometric measurements obtained on stars whose semimajor axis can reach 2 mpc (the semimajor axis of the S2 star is 5 mpc), and whether the astrometric accuracy on data is ≲0.1 mas. If we consider an accuracy of about 20 μas, it should be possible to detect this effect with a star whose semimajor axis can reach 5 mpc. Furthermore, Zucker et al. (2006) considered both astrometric and spectroscopic measurements of several S stars and showed that loworder relativistic effects affecting spectroscopy could be detected when considering monitorings of 10 yr with instruments reaching accuracies of 1.5 mas for astrometry and 25 km s^{1} for spectroscopy. Other theoretical investigations have been done in order to constrain highorder effects such as the LenseThirring effect and the quadrupole moment of the black hole (Kraniotis 2007; Will 2008; Kannan & Saha 2009; Merritt et al. 2010; Angélil & Saha 2010; Angélil et al. 2010). More precisely, Will (2008) demonstrated the possibility of constraining the quadrupole moment of the black hole by observing astrometric positions of at least two stars whose eccentricity and orbital period satisfy and yr, respectively. Angélil et al. (2010) estimated the spectroscopic accuracies necessary to detect various relativistic effects with measurements obtained during one orbital period of S2 (≈15.8 yr). For instance, these authors showed that an accuracy of 10 m/s is needed to constrain the LenseThirring effect in one orbital period. More recently, Zhang et al. (2015) and Yu et al. (2016) have shown it possible to constrain the angular momentum of the black hole by using astrometric and spectroscopic observations of the S2 star. Both authors used a sophisticated stellarorbit model including both the computation of null and timelike geodesics. In addition, works on the detection of gravitational lensing have been performed by Bozza & Mancini (2012) and have shown that this effect is sufficiently important to be detected by GRAVITY; in particular, with the S17 star whose gravitational lensing induces an astrometric shift of about 30 μas in 2018. Besides, studies performed by Bozza & Mancini (2004, 2005, 2012) and Jørgensen et al. (2016) showed that gravitational lensing is affected by the angular momentum of the black hole and can thus lead to a constraint on the LenseThirring effect. However, we will need very accurate instruments to detect such deviations since the astrometric shift reaches only a few microarcseconds.
The aim of this paper is to further develop the investigations performed by these various authors but only focusing on S2. The choice to only work on the S2 star is motivated by the fact that we do not know whether closerin stars will be observed with GRAVITY. It is thus important to extract the maximum information from this star. Considering different astrometric and spectroscopic accuracies for the S2 observations, we estimate the different minimal observation times above which it is possible to detect different relativistic effects. We are thus capable of determining the threshold times needed for the GRAVITY instrument to detect different relativistic effects. The astrometric accuracy of this instrument will improve the detection of relativistic effects and thus better constrain the nature of Sgr A*. This paper can be considered as a first step in the development of the numerical tools necessary to interpret the forthcoming accurate GRAVITY data. The study performed in this paper is done by implementing different models allowing us to describe the future S2 data with different degrees of refinement in the implementation of the various relativistic effects. With such models we will be able to determine whether we can detect the different effects for a given pair of astrometric and spectroscopic accuracies, and thus determine which model can be used to interpret the forthcoming S2 data with minimal computing time. A part of this paper is exclusively devoted to the angular momentum of the central black hole candidate, where we discuss the constraint on its norm and direction obtained by using simulated GRAVITY observations of the S2 star. Contrary to Zhang et al. (2015) and Yu et al. (2016), we use various stellarorbit models and determine which one allows us to investigate the constraint on the black hole angular momentum parameters with minimal computing time whilst still obtaining goodquality fits.
The paper is organized as follows: Sect. 2 is devoted to the explanation on how S2 astrometric and spectroscopic observations are simulated. Section 3 defines the models fitted to the simulated observations and used to detect relativistic effects, and Sect. 4 explains the procedures used to estimate both the threshold times and the angular momentum of the black hole; the different results are also given in this section. Finally, conclusions and discussions are given in Sect. 5.
2. Mock observations of the S2 star
In our study, we consider two observables: the astrometric positions and the radial velocities (spectroscopy) of the S2 star. To generate mock observations of the S2 star, we consider a fully GR model by using the raytracing code Gyoto^{1} (Vincent et al. 2011; Grould et al. 2016). The Kerr metric is considered in this model, and thus all relativistic effects are taken into account, such as pericenter advance, transverse Doppler shift, gravitational redshift, gravitational lensing, LenseThirring effect and the Shapiro time delay (see Appendix E.1 for a brief definition of these effects). The Roemer effect is also naturally taken into account in Gyoto (see Appendix E.1). In this fullGR model, gravitational lensing is obtained by making two assumptions: we neglect multiple images of the star (e.g., the secondary image) and we consider only one photon of the primary image. The first assumption is valid since the influence of multiple images on the astrometric positions of the S2 star is negligible (≲0.5 μas if only considering the secondary image). The second assumption means that we do not compute the flux of the primary image. It is valid since the star is far enough from the black hole to neglect the amplification effect due to lensing. One photon is thus sufficient to recover the astrometric position of the star. For more technical explanations about how we compute astrometric positions and radial velocities of S2 in the fullGR model with Gyoto, see Appendix A.
In this paper, we consider three different reference frames in Fig. 1: the blackhole frame (x_{bh},y_{bh},z_{bh}) centered on the black hole and labeled in KerrSchild coordinates (Visser 2007); the z_{bh}axis is taken along the angular momentum axis of the black hole; the orbit frame (x_{orb},y_{orb},z_{orb}) centered on the black hole and such that (x_{orb},y_{orb}) spans the plane of the orbit and z_{orb} is along the angular momentum of the orbit; and the observer frame (α,δ,z_{obs}) located at the observer position with (α,δ) spanning the observer screen whose origin is located at the center of the screen and which corresponds to the apparent position of the black hole; z_{obs} being directed towards the black hole.
Fig. 1 Illustration of the different reference frames: the blackhole frame (in black) (x_{bh},y_{bh},z_{bh}) labeled in KerrSchild coordinates, where the z_{bh}axis corresponds to the angular momentumaxis; the orbit frame (in blue) (x_{orb},y_{orb},z_{orb}); and finally, the observer frame (in green) (α,δ,z_{obs}). The angles i′ and Ω′ allow to recover the direction of the angular momentum of the black hole; they give the position of the blackhole frame relative to the observer frame. The point P on the orbit of the star denotes the pericenter. 
Fig. 2 Astrometric and spectroscopic observations of the monitorings of one period (upper plots) and 6 months (lower plots) of the S2 star, simulated with the fullGR model. For one period of monitoring we consider one week of observation at pericenter passage (blue points), 6 months of observations around the pericenter passage (red points) and 14 yr of observations for the rest of the orbit (black points). The 6 months of monitoring corresponds to the first nineteen points of the oneperiod run. 
To compute the orbit of S2 with Gyoto we need the initial position and velocity of the star in KerrSchild coordinates. These coordinates are obtained using the two steps listed below:

Step 1: get the threedimensional positions(α_{s},δ_{s},z_{s,obs}) and velocities (v_{αs},v_{δs},v_{zs,obs}) of the star in the observer frame at a given observation date t_{obs}, by using the ThieleInnes formulas (Taff 1985) and the Keplerian orbital parameters: period T, semimajor axis a_{sma}, eccentricity e, time of the pericenter passage t_{p}, inclination i, angle of the line of nodes Ω and angle from ascending node to pericenter ω (see Fig. A.2 for an illustration of these three angles). See Appendix B for a demonstration of how we compute the star coordinates in the observer frame.

Step 2: get the initial coordinates of the star in the blackhole frame. To do so, we apply a rotation matrix to the coordinates of the star obtained in the observer frame. The matrix depends on two angles Ω′ and i′ giving the direction of the angular momentum of the black hole with respect to the observer (see Fig. 1). The rotation matrix is given by (1)where i′′ = 3π/ 2 + i′. We assume that the obtained position and velocity correspond to the initial KerrSchild coordinates of the star.
Knowing the initial coordinates of the star, we can integrate the timelike geodesic in Gyoto, that is, obtain its GR orbit. We note that the Newtonian orbit can be considered as an osculating orbit of the GR one.
The initial coordinates of S2 are obtained considering the orbital parameters equal to the bestfit values evaluated by Gillessen et al. (2009b): T = 15.8 yr, a_{sma} = 0.123′′, e = 0.88, t_{p} = 2002.32 yr, Ω = 225.39°, ω = 63.56°, i = 135.25°. We also consider the distance between the observer and the center of our galaxy equal to the bestfit found by these authors: R_{0} = 8.33 kpc. For the angular momentum parameters of the black hole, we choose a = 0.99, i′ = 45° and Ω′ = 160° where a is the dimensionless spin of the black hole.
The S2 noisy data are obtained by adding a Gaussian random noise to the fullGR observations whose distribution is parametrized by a standard deviation σ_{A} for the astrometry and σ_{V} for the spectroscopy. In this study, we consider different values for σ_{A}: 10 μas, 30 μas, 50 μas and 100 μas; and σ_{V}: 1 km s^{1}, 10 km s^{1} and 100 km s^{1}. We note that several of those accuracies can be reached by current and future instruments (see Table 1, whose principal accuracies are taken from Weinberg et al. 2005; Gillessen et al. 2009a; Genzel et al. 2010; Eisenhauer et al. 2011).
Astrometric and spectroscopic accuracies of various current and future instruments, capable of observing in the near infrared.
As shown in Angélil & Saha (2010, 2011), Angélil et al. (2010), and Zucker et al. (2006), spectroscopic measurements obtained during the pericenter passage are a powerful tool to detect relativistic effects. In particular, Zucker et al. (2006) showed that in the case of S2, the transverse Doppler shift and gravitational redshift represent a significant contribution to radial velocity of about ≈200 km s^{1} near pericenter. These are the reasons why we choose to better sample the mock observations at S2 pericenter passage. We consider different runs of observation ranging from one month to three periods (≈47 yr). All shorter runs are subsets of the threeperiod run. The latter run is sampled, for each period, as follows:

two points per night during one week at pericenter passage(starting in 2018.11 for the first period);

one point per month during six months at pericenter passage (between 2017.78 and 2018.29 for the first period);

one point every four months for the rest of the orbit (between 2018.63 and 2033.20 for the first period).
All runs start between 2017.78 and 2018.11. As an illustration, the monitoring of one period and six months are visible in Fig. 2. We note that astrometric and spectroscopic observations are supposed to be at the same dates, which could be more difficult in practice since radial velocity measurements should not be done with GRAVITY (because of its poorer spectroscopic accuracy, see Table 1).
Fig. 3 Illustration of each model described in Sect. 3.1, only focusing on the different effects affecting the astrometric measurements of the S2 star. The positions of both the star (blue line) and the photon (red line) at two dates t_{1} and t_{2} are represented for models B to G. At some particular observation time t_{obs}, the apparent position of the star on the observer screen, is highlighted by the symbols A to G. The corresponding emitting position of the star along its orbit is also represented. In model A, the difference Δt = t_{obs}−t_{em} is zero since the photon has an infinite speed. In models B to G, the Roemer time delay, noted Δt_{R}, is included, thus Δt depends on it. In model D, the star trajectory is described by GR but only accounts for the pericenter advance (PA) since the LenseThirring effect is neglected. Contrary to model D, model E takes into account this latter effect on the timelike geodesic (LTS for LenseThirring on the Star). Model F is analogous to model E but approximates the gravitational lensing effect (GL). Model G is the fullGR model; gravitational lensing is thus better estimated, the LenseThirring effect on the photon path (LTP) is naturally taken into account, and the quantity Δt depends on the Shapiro time delay, Δt_{S}, in addition to Δt_{R}. 
3. Models
3.1. Definitions
We want to estimate the minimal observation times needed to detect different effects affecting the S2 star observations, and considering different astrometric and spectroscopic accuracies. For doing so, we implement different models, in order of increasing complexity and computing time needed, including three Keplerian models (whose star orbit is Keplerian) and four Relativistic models (whose star orbit is relativistic) listed here:

Model A: Keplerian model without the Roemer time delay due tothe finite speed of light, and without adding relativistic effects.

Model B: Keplerian model only considering the Roemer effect.

Model C: Keplerian model considering the Roemer effect, the transverse Doppler shift, and the gravitational redshift.

Model D: GR model without ray tracing. The orbit of the star is relativistic but not the photon trajectory. The Roemer effect, the transverse Doppler shift, the gravitational redshift and the pericenter advance are taken into account.

Model E: similar to model D but takes into account the LenseThirring effect on the star trajectory.

Model F: similar to model E but considers a supplement effect corresponding to an approximation of gravitational lensing.

Model G: FullGR model described in Sect. 2. All previous effects are naturally taken into account in this model and additional effects are considered such as the Shapiro time delay and the LenseThirring effect on the photon trajectory.
An illustration of each model can be seen in Fig. 3. More detail on how the observations are generated in the different models are given below.
The astrometric and spectroscopic observations generated in model A are simply obtained by using the orbital parameters and the ThieleInnes formulas (see Appendix B).
For model B, we use the same protocol as model A but we include the Roemer time delay solving the following equation (see Appendix C): (2)where t_{em} is the emission date, z_{s,obs}(t_{em}) is the position of the star along the lineofsight at emission, and c is the speed of light. Knowing the observation date t_{obs} and using this equation, we can determine the emission date t_{em} at which we need to compute the astrometric and spectroscopic data (obtained with formulas of Appendix B) to take into account the Newtonian timetraveling of the photon.
Model C is similar to model B but considers two loworder relativistic effects only affecting the spectroscopy: the transverse Doppler shift and gravitational redshift. To implement these effects, we compute an approximated radial velocity allowing to simulate both the Newtonian Doppler shift (also called the longitudinal Doppler shift) and the two relativistic redshifts. This formula is given by (see Appendix D for a demonstration) (3)where ϵ = 2GM/ (c^{2}r_{em}) with r_{em} being the radial coordinate of the star at emission in the blackhole frame given by Eq. (B.1), being the velocity of the star in the blackhole frame and the projection of along the lineofsight in the observer frame (α,δ,z_{obs}). As model C is a Keplerian model, the velocity corresponds to the orbital velocity given by (4)We mention that the Roemer time delay affecting the spectroscopy is included since the three quantities r_{em}, and are computed considering the same protocol as model B.
To generate the star orbit in model D, we use a Schwarzschild metric and the procedure described in Sect. 2: compute a GR orbit with the raytracing code Gyoto considering an initial position of the star generated with the Keplerian orbital parameters and the ThieleInnes formulas. In this model, we do not compute null geodesics, only the star trajectory is relativistic. The pericenter advance is thus naturally taken into account. As in model C, we simulate the Roemer time delay, the transverse Doppler shift and the gravitational redshift by using the Eqs. (2)and (3). More precisely, the emission date given by Eq. (2)is used to recover the star coordinates (evaluated with Gyoto) allowing to simulate astrometric positions of this star (by projecting these coordinates in the plane of the sky) affected by the Newtonian travel time of the photon. In Eq. (3), the orbital velocity is computed through the GR expression where g_{ij} are the spatial metric coefficients and is the threevelocity of the star defined as u^{i}/u^{t} with u being the fourvelocity of the star. The projection of this threevelocity, , is obtained considering a photon not affected by the spacetime curvature (i.e., we simply project the star velocity in Euclidian space along the lineofsight). We note that the three effects mentioned above are better estimated in model D than in model C since the coordinates of the star (position and velocity) used in both Eqs. (2)and (3)are affected by the pericenter advance.
Model E is equivalent to model D but the star orbit is obtained by using a Kerr metric instead of a Schwarzschild metric. The LenseThirring effect on the star trajectory is thus naturally considered in such a model.
Model F is similar to model E but approximatively simulates the astrometric shift of the star induced by gravitational lensing. It is taken into account by using analytical approximations developed by Sereno & de Luca (2006). These formulas are obtained in weakdeflection limit, meaning that the minimal distance between a photon and the black hole is higher than the Schwarzschild radius R_{S} = 2GM/c^{2}. In addition, they are developed in the weakfield regime, the observer and the emitter are thus considered in flat spacetime.
The different effects taken into account in each model are summarized in Table 2.
Effects considered in each model.
3.2. Effects affecting the astrometry of the S2 star
In this section we discuss the different effects that change the astrometric position of the S2 star. In particular, we focus on the pericenter advance, the Roemer and Shapiro time delay, the LenseThirring effect, and gravitational lensing. All plots in this section and the following one, devoted to spectroscopy, are obtained with black hole and orbital parameters given in Sect. 2. See Appendix E.2 for a brief explanation on how we evaluate the astrometric contribution (and spectroscopic one) of each effect.
The pericenter advance is an effect which increases with the number of turns made by the star around the black hole. To evaluate this effect we can compare a Keplerian orbit with a GR orbit. That is what is presented in Fig. 4. We can see that the maximal astrometric difference is located near pericenter passages. At first pericenter passage, the difference is weak since it corresponds to the first data points: the apparent positions start to differ when the orbit has evolved. During the first S2 period, the maximal magnitude of the pericenter advance is reached near the second pericenter passage (last magenta dot in Fig. 4) and is equal to ≈3 mas. At second and third periods, the maximal impact is of about 8 mas and 16 mas, respectively.
Fig. 4 Astrometric impact of the pericenter advance on the S2 star observed during three orbital periods. Solid magenta, blue, and black circles correspond to the first, second, and third S2 period, respectively. Open red circles and triangles represent the apocenter and pericenter passages, respectively. 
Fig. 5 Astrometric impact of different effects on the S2 star observations obtained during three periods. First plot: Roemer time delay. Second plot: Shapiro time delay. Third plot: LenseThirring effect considering the angular momentum parameters a = 0.99, i′ = 45° and Ω′ = 160°. Fourth plot: gravitational lensing. Open red circles and triangles denote the position of the apocenter and the pericenter, respectively. See Appendix E.2 for the precise definition of the various quantities used in these plots. 
The two first plots in Fig. 5 present the time delay effects on the S2 astrometric positions: the left plot shows the influence of the Roemer effect and the right plot gives the impact of the Shapiro effect. We can see that the influence of this first effect is always higher than the best astrometric accuracy of GRAVITY (10 μas) and can reach ≈450 μas. It shows that this contribution cannot be neglected in stellarorbit models used to interpret the GRAVITY data. We remind that this effect is not a GR effect. We note that the Roemer time delay has already been considered to treat the data obtained on the S stars by Ghez et al. (2008) and Gillessen et al. (2009b). The shift induced by the Shapiro time delay is always lower than 10 μas. The maximal offsets are reached near pericenter at ≈8 μas.
On the third plot in Fig. 5 we can see the LenseThirring effect, where we considered a = 0.99, i′ = 45° and Ω′ = 160°. It takes into account the black hole angular momentum effect on both the S2 stellar orbit and the photon path. However, the contribution of the LenseThirring effect on the null geodesic is negligible: the astrometric impact of this effect on both Shapiro time delay and gravitational lensing is ≪1 μas. Thus, the plot only shows the LenseThirring impact on the S2 timelike geodesic. This effect is similar to pericenter advance since it increases with the number of turns made by the star around the black hole. During the first, second and third period it reaches ≈10 μas, ≈25 μas and ≈40 μas near apocenter passage, respectively. Near pericenter passages, the offset is always lower than 10 μas. It shows that in the black hole configuration considered here (a = 0.99, i′ = 45° and Ω′ = 160°), the LenseThirring effect is negligible near pericenter meaning that it is important to observe near apocenter if we want to investigate a constraint on the black hole angular momentum parameters. A deeper analysis, considering different values for the parameters a, i′ and Ω′, will be given in Sect. 4.2.2. We note that the results discussed here are consistent with those given in Zhang et al. (2015) and Yu et al. (2016).
The last plot in Fig. 5 shows the gravitational lensing effect on apparent positions of the S2 star. The maximal shifts are reached near pericenter passages and are of about 20 μas. Most of the time, gravitational lensing is as low as 2 μas; close to apocenter passages, they reach ≈1 μas.
To summarize, maximal astrometric offsets due to each relativistic effect are listed in Table 3.
Maximal astrometric offsets in μas reached at first, second and third periods of the S2 star, due to different relativistic effects and considering the black hole angular momentum parameters a = 0.99, i′ = 45°, Ω′ = 160°.
3.3. Effects affecting the spectroscopy of the S2 star
Here, we want to show the influence of different effects on the S2 radial velocity measurements. We focus on the pericenter advance, the Roemer time delay, the transverse Doppler shift, the gravitational redshift and the LenseThirring effect. We also look for the impact of using an approximated radial velocity (see Eq. (3)). More precisely, it means that we are interested in the cumulative spectroscopic impact of highorder effects affecting the photon trajectory such as the Shapiro time delay; we denote these cumulative contributions HighOrder Photon Curvature (HOPC).
As for the astrometry, the impact of the pericenter advance on measured radial velocity increases with the number of orbits made by S2. Figure 6 shows the spectroscopic shift obtained when comparing radial velocities evaluated with a Keplerian orbit with those estimated with a GR orbit. The first maximal shift appears near the second pericenter passage (last magenta dot in Fig. 6) and is of about 140 km s^{1}. At second and third periods, the maximal offset reaches ≈1520 km s^{1} and ≈2800 km s^{1}, respectively.
Fig. 6 Spectroscopic impact of the pericenter advance on the S2 star observed during three orbital periods. Solid magenta, blue, and black circles correspond to the first, second, and third S2 period, respectively. Open red circles and triangles represent the apocenter and pericenter passages, respectively. 
Fig. 7 Spectroscopic impact of different effects on the S2 star observations obtained during three periods. First plot: Roemer time delay. Second plot: transverse Doppler shift and gravitational redshift. Third plot: LenseThirring effect considering the angular momentum parameters a = 0.99, i′ = 45° and Ω′ = 160°. Fourth plot: HOPC contributions. Open red circles and triangles denote the position of the apocenter and the pericenter, respectively. See Appendix E.2 for the precise definition of the various quantities used in these plots. 
The first plot in Fig. 7 presents the influence of the Roemer time delay on radial velocity of the S2 star. The maximal shifts reach ≈50 km s^{1} near pericenter passages.
The second plot in Fig. 7 shows the cumulative influence of the transverse Doppler shift and the gravitational redshift. Here again, the maximal offsets are reached near pericenter. The highest values are ≈200 km s^{1} which is consistent with Zucker et al. (2006). Near apocenter passages the shifts are lower than 10 km s^{1}.
The LenseThirring effect on radial velocities is visible on the third plot in Fig. 7. It essentially shows the shift due to its impact on the star trajectory. We see that at each period, the absolute maximal shifts are reached near pericenter passages, but the influence of this effect is very low and always below 1 km s^{1}. These results are similar to those found by Angélil & Saha (2010), Zhang et al. (2015) and Yu et al. (2016).
The last plot in Fig. 7 corresponds to the impact of the HOPC contributions. The maximal offsets are reached near pericenter and are ≈5 km s^{1}. As for astrometry, the LenseThirring effect on the photon trajectory essentially does not modify the radial velocities of the S2 star (≈10^{2} km s^{1}). The observed shift here is thus due to other approximations such as the Shapiro time delay.
Maximal spectroscopic offsets obtained with each relativistic effect are listed in Table 4.
We can conclude that all relativistic effects on radial velocity are maximal near S2 pericenter passages (but not exactly at pericenter). However, considering the astrometry, the LenseThirring effect is maximal near apocenter passages. Moreover, the LenseThirring effect on the photon path is negligible for both astrometry and spectroscopy. Concurrently, the Schwarzschild spacetime curvature and the LenseThirring effect on the star trajectory are not negligible. In particular, this latter effect should be marginally detected at first apocenter passage by the GRAVITY instrument (≈10 μas − in 2026).
3.4. Differences between models A to F and the fullGR model
We report in Table 5 the maximal astrometric and spectroscopic differences between models A to F and model G.
The maximal astrometric differences between Keplerian models (A–C) and the fullGR model are due to the fact that the pericenter advance is not taken into account in those models. Spectroscopic differences are also mainly dominated by the absence of the pericenter advance. However, we note variations between Keplerian models, which are explained by the fact that models B and C take into account different effects: model B includes the Roemer time delay, and model C considers the Roemer time delay, the transverse Doppler shift and the gravitational redshift.
At the first S2 period, the astrometric difference between models D and G is dominated by the absence of gravitational lensing. At second and third periods, it is dominated by the absence of the LenseThirring effect on the star trajectory. For model E, the astrometric difference is only due to the fact that it neglects gravitational lensing. Finally, for model F, the astrometric shift is induced by the fact that gravitational lensing is not completely reproduced by the approximations of Sereno & de Luca (2006). Radial velocities of models D to F are shifted from those computed with model G by ≈5 km s^{1}. This offset is due to the HOPC contributions.
Maximal astrometric and spectroscopic differences of each model with respect to the fullGR model (model G) considering three orbital periods of the S2 star.
4. Fitting
In the following sections, we give an estimation of the threshold times above which we can detect the different effects discussed above with observations of the S2 star, and considering various astrometric and spectroscopic accuracies. We also focus on the constraint on the black hole angular momentum parameters with this star. In this section, we explain the methods used to estimate these different threshold times, and to constrain the angular momentum.
4.1. Procedures
Each model is described by eight parameters corresponding to the seven orbital parameters and the distance R_{0} between the observer and the black hole. The mass of the black hole is not an individual parameter to fit but varies through the third Kepler’s law: . Models E–G are described by three other parameters since they take into account the LenseThirring effect. Those parameters are the norm and direction of the angular momentum of the black hole: a, i′ and Ω′.
To estimate the different minimal observation times needed to detect the various effects, we determine the threshold times above which the models A to F fail to fit the fullGR observations generated with model G. We mention that the fitting method used for Keplerian models and relativistic models will be different. The reason for this is explained below.
The fitting procedure used to fit models A–C is the LevenbergMarquardt method (Levenberg 1944) based on the leastsquares method. The quantity to minimize is (5)where N is the number of observation dates. The quantities labeled obs correspond to the mock observations generated with the fullGR model. Those labeled m are obtained with models A–C. The different accuracies (σ_{A},σ_{V}) considered to estimate the different threshold times are listed in Sect. 2. The initial parameters (initial guess) considered for the fitting of models A−C are given by P_{init} = P_{Gillessen} + 1σ_{Gillessen} where P_{Gillessen} is a vector containing the bestfit parameters of Gillessen et al. (2009b) given in Sect. 2, and 1σ_{Gillessen} is the vector containing the 1σ error of each parameter also estimated by Gillessen et al. (2009b): σ_{T} = 0.11 yr, σ_{asma} = 0.001′′, σ_{e} = 0.003, σ_{tp} = 0.01 yr, σ_{Ω} = 0.84°, σ_{ω} = 0.84°, σ_{i} = 0.47°, σ_{R0} = 0.48 kpc.
If the model fits the data well, the distribution of χ^{2} obtained by evaluating this quantity several times must follow a χ^{2} law with k degrees of freedom^{2} where k = 3N−n, with n being the number of fitted parameters and where the factor three corresponds to the fact that we consider both astrometry (α, δ) and spectroscopy; see the red curve in Fig. 8 for an illustration of such a distribution, when considering a reduced χ^{2} defined by (6)To distinguish a bad fit from a good fit we can use the χ^{2} test. This test allows us to determine whether one evaluation of the quantity χ^{2} is consistent with the assumption that it is drawn from a χ^{2} law. In other words, it tests the null hypothesis given by (7)To test the null hypothesis H_{0}, we determine a limit χ^{2} noted by solving the equation (8)where is the cumulative distribution function (CDF) of the χ^{2} law with k degrees of freedom, and p is the probability of incorrectly rejecting H_{0}, which we fix at 5%. Thus, if the condition is verified, knowing H_{0}, we incorrectly reject the model in 5% of cases. It means that if we fit for instance a model 1000 times to data generated with this same model, there will be 5% of the 1000 estimations of which will be superior to (see the red curve in Fig. 8). We specify that the 1000 are obtained by fitting the model to 1000 runs of observation which have the same duration in time, but whose noise differs by using a random draw of a normal law. On Table 6 several values of are evaluated for the models A to F and different runs of observation.
Fig. 8 Red curve: distribution of obtained with a model describing the observations well. The distribution follows a reduced χ^{2} law. The value (see the text for a definition of this quantity) is marked by the dashed line and is obtained considering p = 5% in Eq. (8). Blue curve: example of a distribution when the model does not describe the observations well. 90% of the are higher than the . For both curves, we consider k = 100. 
The χ^{2} test requires only one estimation of χ^{2}. However, we decide to use a method allowing to meaningfully reject the null hypothesis. More precisely, we chose to fit a Keplerian model 100 times to observations generated with model G, and to reject the model if 90% of the are superior to (see the blue curve in Fig. 8). The minimal observation time above which we consider that the model fails to describe the data is thus equal to the time duration of the run where 90% of the satisfy . Strictly speaking, if we consider an infinite number of realizations and that more than 5% of the satisfy , then the model can be rejected. However, in this paper we consider a finite number of fits, and we thus choose to be conservative and claim that a model fails when it reaches 90%. In addition, we chose to report the minimal observation times obtained when 60% of the verify , since the model can already be rejected with sufficiently high confidence at such a percentage.
In order to validate the procedure defined here, we fit model C 100 times to observations obtained also with this model. The aim is to verify that we find 5% of the satisfying . For doing so, we chose two runs of observation: 32 yr (≈2 periods of S2) with (σ_{A},σ_{V}) = (10μas, 1 km s^{1}) and 4 months with (σ_{A},σ_{V}) = (100μas, 100 km s^{1}). In both cases, we find ≈6% of the satisfying which is very close to the expected value. Such results show that 100 fits are sufficient to get a first estimation of the minimal observation times.
estimated for different models and runs of observation.
We also apply this test to model E: fit 100 times model E considering the LevenbergMarquardt method, to observations generated with this model. In this case, we find that the fitting method is not appropriate since we obtain ≈35% of the which satisfy when considering for instance the run of 32 yr and (σ_{A},σ_{V}) = (10μas, 1 km s^{1}). This percentage is explained by the fact that the parameter space is more difficult to probe with relativistic models. Indeed, if we repeat the test but consider the initial guess to be equal to the parameters used to generate the observations instead of P_{init} = P_{Gillessen} + 1σ_{Gillessen}, we find ≈6%. The LevenbergMarquardt method is thus not appropriate for highly nonlinear models such as relativistic models, since the initial guess needs to be chosen close to the solution. This is the reason why we chose another fitting method for models D–F. We point out that we do not use the same method for all models because the computing time needed by the fitting method used for relativistic models is more important.
To estimate the minimal observation times above which the relativistic models fail to describe the observations obtained with model G, we use the same protocol as for Keplerian models but we consider a Monte Carlo Markov chain (MCMC) method to fit the different models (Binder 2002). More precisely, we use the emcee^{3} software allowing MCMC simulations using the Affine Invariant Ensemble Sample method proposed by ForemanMackey et al. (2013). We specify here that if we note the observations, O, and the vector containing the parameters of a model chosen during the MCMC, P, the posterior probability density π(P  O) is given by (using the Bayesian theorem) (9)where f(O  P) is the likelihood function expressed as (10)and π(P) is the prior probability density of the parameters P which we chose distributed according to a uniform law: the orbital parameters and R_{0} are uniformly chosen between P_{Gillessen} − 6σ_{Gillessen} and P_{Gillessen} + 6σ_{Gillessen}, and the angular momentum parameters are uniformly chosen in there own domain of variation: a ∈ [0,1]^{4}, i′ ∈ [0°,180°] and Ω′ ∈ [0°,360°].
In order to check the fitting method used for models D to F, we compare the percentages of verifying the condition and obtained with both the LevenbergMarquardt and MCMC methods. For doing so, we fit model C 100 times to observations generated with model G and obtained during one period with (σ_{A},σ_{V}) = (30μas, 10 km s^{1}). For both methods we find that 45% of are superior to the . These results show that the fitting obtained with emcee is in accordance with that obtained with LevenbergMarquardt.
Besides, we apply the same test as done in the Keplerian part: fit model E 100 times to observations generated with this model and obtained during ≈2 periods of S2 with (σ_{A},σ_{V}) = (10μas, 1 km s^{1}). In this case, we find ≈6% instead of the 35% obtained with the LevenbergMarquardt method. We can thus say that the MCMC method is more appropriate for relativistic models than that of LevenbergMarquardt.
In the section devoted to the black hole angular momentum constraint, we also use the emcee software. Since the majority of the parameters are already constrained by Gillessen et al. (2009b), we chose to vary the different parameters in a domain whose bounds are defined by P_{Gillessen} ± 1σ_{Gillessen}. For the black hole angular momentum parameters, they vary in their domain of variation: a ∈ [0,1], i′ ∈ [0°,180°] and Ω′ ∈ [0°,360°]. The different fittings will be done considering accuracies that can be used with current instruments:10−30μas and 10 km s^{1}. The model used to constrain the eleven parameters is model F. The advantage of this model is the computing time since it does not use ray tracing. We remind that the maximum differences between models F and G is of about 7 μas and 5 km s^{1}. In spite of those differences, the percentage of verifying when considering the runs of observation of three periods with (σ_{A},σ_{V}) = (10μas, 10 km s^{1}) is only ≈10%. This shows that model F seems sufficient to describe observations obtained with model G and thus appropriate to investigate the constraint on the norm and direction of the angular momentum of the black hole.
Estimations of the threshold times needed to detect different effects with the S2 star, considering various astrometric and spectroscopic accuracies.
4.2. Results
4.2.1. Constraint on various effects
The aim of this section is to estimate the minimal observation times required to detect different effects acting on the S2 star astrometric and spectroscopic observations, at 12 given pairs of accuracies (σ_{A},σ_{V}). We will thus be capable of determining the threshold times needed for GRAVITY to detect relativistic effects. We remind that the different results are obtained by fitting models A to F to fullGR S2 observations generated with model G. We mention that we do not compute threshold times of relativistic models for all pairs of accuracies (σ_{A},σ_{V}) because the accuracies that we consider are sufficient to make conclusions on the possibility of constraining different effects with these models.
Table 7 gives the different threshold times obtained for models A to F and various pairs of accuracies. In what follows, we consider that an effect is detectable when a model X, neglecting this effect, fails after a significantly shorter period than a model X+1 more sophisticated, taking into account this effect. The detection of an effect is thus obtained by comparing two different models. The left column of Table 7 gives some of the effects missing in a model and that cause its failure to reproduce the S2 observations. We mention that a conclusion of all results discussed in the following paragraphs is given at the end of this section.
First, Table 7 clearly shows that the time threshold for telling an effect grows with poor spectral and astrometric accuracies. This is obvious: more observation time is needed to demonstrate that data with low quality are at odds with a given model.
Second, provided the spectroscopic accuracy is ≲10 km s^{1}, the threshold times increase mostly drastically for models C to F compared to models A and B. This means that models C to F become much better at describing the observations than models A and B. Models A and B are the only ones that do not contain any relativistic effects. This shows that even for rather poor astrometric accuracies (of order 100 μas), lowestorder relativistic effects will be at hand after only a few months of monitoring, provided spectroscopic accuracy is ≲10 km s^{1}.
Finally, let us now compare the time thresholds for successive models in order to determine the minimum observation times needed to tell the various effects.
Roemer effect: this effect is tested by comparing the results of models A and B. We note that for spectroscopic accuracies ≲10 km s^{1}, models A and B have similar behaviors. This means that at such high spectroscopic accuracies, some relativistic effects dominate the Roemer time delay. Only at the lowest spectroscopic accuracy (100 km s^{1}) is the Roemer effect strong enough to be detectable; within a few months to a few years depending on the astrometric accuracy. Indeed, the threshold times differ most significantly between models A and B at such accuracies. For instance, at (σ_{A},σ_{V}) = (10 μas, 100 km s^{1}), when model A reaches 90% rejection after 4 months, model B has only reached 40% rejection. This shows that the detection of the effect is strong.
Relativistic redshifts: we now compare models B and C. For spectroscopic accuracies ≲10 km s^{1}, relativistic redshifts appear after only a few months of observations whatever the astrometric accuracy. Again, the very different time thresholds between the two models allows us to get a strong detection of the Doppler transverse shift and gravitational redshift. We add that the rejection percentages estimated with model C at runs of observation where model B fails to reproduce the observations at (σ_{A},σ_{V}) = (10−100μas, 1−10 km s^{1}) are all inferior to 20%, which also supports the strong detection of the relativistic redshifts. This discussion also holds for model A since it behaves very similarly to model B at these accuracies.
Pericenter advance: this effect is tested by comparing models C and D. At very high astrometric accuracy (10 μas), models C and D behave either very similarly (for high spectral resolution) or rather similarly (for medium spectral resolution). For this latter case, the threshold times at 90% rejection are different by a factor of three, but the thresholds at 60% rejection are the same. This shows that the effect is detected only weakly. Moreover, it highlights that higherorder effects such as the gravitational lensing and HOPC contributions dominate the pericenter advance since threshold times are mainly equal for both models. Only at a lower astrometric accuracy of 30 μas can the pericenter advance be detected clearly, within ~15 to 20 yr (i.e., around one orbital period of S2) depending on the spectroscopic accuracy. If we want to conclude on the detection of the pericenter advance at accuracies (σ_{A},σ_{V}) = (10 μas, 10 km s^{1}), we need to compare models C and D taking into account all highorder effects in both models by using the raytracing code Gyoto. In such a case, model D becomes a fullGR model and the only missing effect in model C is the pericenter advance. In order to get an estimation of the threshold times that could be obtained with such models at those accuracies, we fit the current model C to observations of the S2 star generated with the current model D. Indeed, the only missing effect needed to reproduce the observations of S2 with model C will be the pericenter advance. The different results are given in square brackets in Table 7. They show that the pericenter advance should be detectable within a few years when considering accuracies (σ_{A},σ_{V}) = (10 μas, 10 km s^{1}).
LenseThirring: Let us now compare models D and E. The threshold times are very similar for all accuracies. This demonstrates that the LenseThirring effect is not detectable for runs of observation 30 yr, and for the considered astrometric and spectroscopic accuracies. The obtained threshold times also show that both models D and E fail to reproduce the data due to highorder effects corresponding to the gravitational lensing and the HOPC contributions.
Gravitational lensing: this effect can be detected by comparing models E and F and by providing a very high astrometric accuracy (10 μas). A strong detection within a few years (p = 60%) is possible for a spectroscopic accuracy of 10 km s^{1}. We add that at (σ_{A},σ_{V}) = (10μas, 10 km s^{1}), when model E reaches 60% rejection after 6 yr, model F has only reached 20% rejection which confirms the strong detection of gravitational lensing. Its detection is weaker at 1 km s^{1}, showing that higherorder effects, corresponding to the HOPC contributions, are involved.
Fig. 9 LenseThirring effect on astrometric (upper plots) and spectroscopic (lower plots) observations of the S2 star, considering a = 0.99 and various Ω′ and i′. The angular momentum direction i′ = 45° and Ω′ = 160° considered for the mock observations are visible on each plot and are denoted with the solid red circles. The different types of curves on each plot correspond to different values of the angle Ω′. For i′ = [0°,180°], we obtain the same curves whatever Ω′ for the astrometry but they are different when considering the spectroscopy: thin and thick lines correspond to i′ = 180° and i′ = 0°, respectively. On the other plots, solid, dotted, dashed and dashdotdotted curves correspond to Ω′ equal to 0°, 45°, 90° and 135°, respectively. Open circles on all plots denote the position of the apocenter and open triangles denote the pericenter passages. 
In order to determine whether the threshold times obtained with models A–C depend on the initial guess used in the LevenbergMarquardt method, we chose these parameters far from the solution: P_{init} = P_{Gillessen} + 5σ_{Gillessen} (instead of P_{init} = P_{Gillessen} + 1σ_{Gillessen}, see Sect. 4.1). For this test, we consider model C and the pairs of accuracies (σ_{A},σ_{V}) = (10μas, 1−100 km s^{1}). We find similar results to those estimated considering the previous initial guess: 10 months for σ_{V} = 1 km s^{1} and 6 yr for σ_{V} = 10−100 km s^{1}. This test tends to show that fittings of models A–C have converged to the global minimum since the threshold times are similar for both initial guesses considered. We also reiterate that we showed in Sect. 4.1 that the percentage found with model C and the LevenbergMarquardt method, at (σ_{A},σ_{V}) = (30μas, 10 km s^{1}), and the run of observation of one orbital period, was similar to the percentage found with the MCMC method (45%), which supports the fact that the fitting method used for Keplerian models converges to the global minimum.
Moreover, we also investigate the influence of the pericenter passages sampling on the threshold time. For doing so, we again use model C and accuracies (σ_{A},σ_{V}) = (10μas, 1−100 km s^{1}). The new sampling is similar to the previous one described in Sect. 2 but we consider only one point per night during the three weeks where the pericenter passages are observed, instead of two points: there are a total of 7 data points instead of 14 at S2 pericenter passages. The minimal observation times obtained are similar to those estimated with the previous sampling: 6 months for σ_{V} = 1 km s^{1} and 8 yr for σ_{V} = 10−100 km s^{1} with the new sampling, and 10 months for σ_{V} = 1 km s^{1} and 6 yr for σ_{V} = 10−100 km s^{1} with the previous sampling. This shows that sampling at pericenter weakly impacts the results. However, it remains essential to correctly sample during pericenter since the majority of relativistic effects are maximal near S2 pericenter passages (see Tables 3 and 4).
To summarize, we can say that if we consider S2 astrometric and spectroscopic observations starting in 2018, we can detect the Roemer effect within 4 months by using models A and B at (σ_{A},σ_{V}) = (10μas, 100 km s^{1}). Relativistic redshifts can be detected within 2 months by using models B (or A) and C with (σ_{A},σ_{V}) = (10μas, 10 km s^{1}). Gravitational lensing can be detected by using models C (or D, or E) and F within ≈4 yr at (σ_{A},σ_{V}) = (10μas, 10 km s^{1}). Pericenter advance should be detected within 8 yr with (σ_{A},σ_{V}) = (10μas, 1−10 km s^{1}) and using modified models C and D taking into account the computation of null geodesics. HOPC contributions can be detected within 6−10 months with (σ_{A},σ_{V}) = (10μas, 1 km s^{1}) and with model C (or D, or E, or F). Regarding the LenseThirring effect, it is not detectable if we consider observations obtained during 2 periods of S2. However, this result does not exclude the possibility of getting a first constraint on the angular momentum parameters of the black hole, which is discussed in the following section.
4.2.2. Constraint on the black hole angular momentum
As noticed in Tables 3 and 4, the LenseThirring effect impacts weakly the photon trajectory. However, these results are obtained only considering the angular momentum direction (i,Ω′′) = (45°,160°). In order to validate the fact that null geodesics are weakly affected by this effect we compute its impact on both astrometry and spectroscopy considering various pairs (i′,Ω′). We found a negligible shift on both observables: <1 μas and <1 km s^{1}, which shows that ray tracing is not primordial to constrain the angular momentum parameters of the black hole with S2. These results are in accordance with Zhang et al. (2015) and Yu et al. (2016).
Now, we are interested in the LenseThirring effect on the star trajectory, and its effect on the astrometric observations of the S2 star. We remind that the astrometric shift is maximal near the three apocenter passages (see the third plot in Fig. 5). However, we again considered one direction for the angular momentum of the black hole. The upper plots in Fig. 9 give astrometric shifts induced by the LenseThirring effect and obtained for several directions (i′,Ω′). We note in each case that the shift is maximal near apocenter passages. It is thus important to observe near to S2 apocenter to constrain the parameters a, i′ and Ω′ with this star. However, as mentioned previously, the majority of the relativistic effects are maximal near pericenter. Astrometric observations near pericenter are thus also necessary to investigate a constraint on angular momentum parameters. More precisely, we need to get a strong constraint on orbital parameters and both distance and mass of the black hole by using relativistic effects observed near pericenter if we want to constrain the LenseThirring effect near apocenter. On the upper plots in Fig. 9, we note that the LenseThirring effect reaches between 10 μas and 40 μas during the three orbital periods for some values of (i′,Ω′). However, the astrometric shift can be less than 10 μas throughout the three S2 periods. This shows that the detection of this effect will be possible only for particular parameters (a,i′,Ω′). See for instance Yu et al. (2016) to get pairs of angles (i′,Ω′) that could be favorable for detecting the LenseThirring effect when considering a = 0.99. If we refer to the results found by Broderick et al. (2011), the angular momentum parameters of the black hole candidate located at the center of our galaxy should be: a = 0,0^{+ 0,64 + 0,86}, and where the errors are those obtained at 1σ and 2σ. We therefore determined the maximal astrometric shift induced by the LenseThirring effect and obtained considering the values of (a,i′,Ω′) allowed by Broderick et al. (2011). First, if we consider the set of values of (a,i′,Ω′) contained in the interval ± 1σ; the impact of this effect is maximal for (a,i′,Ω′) = (0.64,68° + 5°,−52°−15°). The shift reaches 5 μas, 10 μas and 15 μas at first, second and third apocenter passages, respectively. If we consider now the set of values of (a,i′,Ω′) contained in the interval ± 2σ, the impact of this effect is maximal for (a,i′,Ω′) = (0.86,68° + 9°,−52°−24°): the shift reaches 5 μas, 12 μas and 18 μas at first, second and third apocenter passages, respectively. These results show that it seems difficult to strongly constrain the black hole angular momentum parameters with the GRAVITY instrument, if we consider astrometric observations obtained on three orbital periods of S2 and values of (a,i′,Ω′) predicted by Broderick et al. (2011).
Degeneracies of the astrometric shift associated to the angles i′ and Ω′ are also observed, but not all of these degeneracies are visible on the upper plots in Fig. 9. First, if we observe the rotation axis of the black hole from the top (i′ = 0°) or from the bottom (i′ = 180°), the shift is the same whatever the angle Ω′ (see solid curve in Fig. 9 for i′ = [0°,180°]). If we observe the rotation axis edge on (i′ = 90°), there is a central symmetry with respect to the center of the plane of the sky. Indeed, the shifts at (i′ = 90°,Ω′) and (i′ = 90°,180° + Ω′) are similar. Finally, the astrometric shifts are similar when considering (i′,Ω′) = (45°,135°) and (i′,Ω′) = (135°,315°): there is a degeneracy between (i′,Ω′) = (45°,Ω′) and (i′,Ω′) = (135°,180° + Ω′). To summarize, there are three degeneracy groups present whatever the norm of the angular momentum of the black hole (expect for a = 0):

(i′ = 0°,∀Ω′) and (i′ = 180°,∀Ω′);

(i′ = 90°,Ω′) and (i′ = 90°,180° + Ω′);

(i′,Ω′) and (180°−i′,180° + Ω′) if 0°<i′< 90° or 90°<i′< 180°.
These degeneracies have also been noticed independently by Yu et al. (2016). In addition to these degeneracies, shifts obtained with different triplets (a,i′,Ω′) are also similar. For instance, both shifts obtained with (a,i′,Ω′) = (0.99,135°,45°) and (a,i′,Ω′) = (0.7,90°,135°) reach ≈6 μas, 15 μas and 20 μas at first, second and third periods, respectively.
Constraints obtained by fitting model F to astrometric and spectroscopic observations of the S2 star, generated with model G.
Fig. 10 1D and 2D probability distributions obtained by fitting model F to observations during three orbital periods of the S2 star generated with model G and (σ_{A},σ_{V}) = (10μas, 10 km s^{1}). The angular momentum parameters of the black hole are (a,i′,Ω′) = (0.99,45°,160°). Blue dots and lines on each plot correspond to the parameters used to generate the mock observations.The darker the delineated area on 2D distributions, the denser the area. 
On the bottom plots in Fig. 9 the spectroscopic shifts are visible, due to the LenseThirring effect for various directions (i′,Ω′). The maximal shift is in this case located near pericenter of the S2 star. Moreover, whatever the angles i′ and Ω′, the shift is always inferior to 1 km s^{1}. It is thus necessary to get an accuracy better than (or close to) 1 km s^{1} if we want to investigate the constraint on the angular momentum parameters with spectroscopy. Such accuracy is not available with current instruments. We mention nonetheless that the spectroscopic accuracy intended for the EELT is of about 1 km s^{1}. Degeneracies listed above are almost all broken when considering spectroscopic measurements. For instance, the shifts are different (up to ≈2 km s^{1}) for (a,i′,Ω′) = (0.99,45°,135°) and (a,i′,Ω′) = (0.99,180°−45°,135°), which is not the case for astrometry. The only degeneracy is observed for i′ = 0° or i′ = 180°. Indeed, we obtain the same shifts whatever the angle Ω′ at those i′ (see the thick and thin lines on the first bottom plot in Fig. 9).
To summarize, to constrain the black hole parameters (a,i′,Ω′) with current instruments it is better to consider both S2 astrometric and spectroscopic observations obtained at both pericenter and apocenter. These two observables obtained at pericenter should allow for optimization of the constraint on the orbital parameters, and the mass and distance of the black hole. The astrometric data obtained at apocenter should allow us to constrain the angular momentum parameters of the black hole. However, the detection of the LenseThirring effect will only be possible for privileged triplets (a,i′,Ω′) such as (0.99,45°,160°). By using results of Yu et al. (2016), we can estimate that if we consider a high norm for the angular momentum, ≈40% of pairs (i′,Ω′) seem favorable to constrain the LenseThirring effect.
The aim now is to determine whether it is possible to strongly constrain the angular momentum parameters of the black hole by using astrometric and spectroscopic observations of the S2 star, and considering accuracies available with current instruments: 10−30μas and 10 km s^{1}. Contrary to Zhang et al. (2015) and Yu et al. (2016), we want to estimate these parameters without using ray tracing, which means using a stellarorbit model that does not require an important amount of time for the computation. More precisely we chose to use model F allowing to reproduce at best the model G observations (see the end of Sect. 4.1).
Table 8 gives the results of the different fittings, where we considered various norm a, observation time, and astrometric accuracies for the run. In Figs. 10 and 11 we can see the 1D and 2D probability distributions of the different parameters for (σ_{A},σ_{V}) = (10μas, 10 km s^{1}), and a = 0.99 and a = 0.7, respectively. We can see that all parameters on both figures are well fitted. In particular, the 1σ error of the norm a for both fittings is of about 0.1. For the angles i′ and Ω′ the 1σ error is around 15° and 30°, respectively. Other black hole parameters such as the mass and the distance are also correctly constrained by model F: σ_{R0}/R_{0} ≈ 4 × 10^{4} and σ_{M}/M ≈ 10^{3} with σ_{R0} and σ_{M} the 1σ errors; the constraints are improved by a factor ≈100 with respect to Gillessen et al. (2009b). These results are similar to those estimated by Yu et al. (2016) whose constraints are obtained by fitting a stellarorbit model including ray tracing, on mock observations also generated with this model. We mention that the parameters used by these authors to simulate the S2 observations are nearly similar to those used for our study. The difference between these authors and this study is the number of data points: 120 in Yu et al. (2016) and 192 in this paper, and the fact that we neglect the photon path curvature but approximate the gravitational lensing effect by using analytical formulas.
The 1σ errors of the angular momentum parameters obtained when considering (a,N_{p},σ_{A},σ_{V}) = (0.99,3,30 μas, 10 km s^{1}) are σ_{a} ≈ 0.26, σ_{i′} ≈ 30° and σ_{Ω′} ≈ 40°. Those estimated for (a,N_{p},σ_{A},σ_{V}) = (0.99,2,10 μas, 10 km s^{1}) are σ_{a} ≈ 0.2, σ_{i′} ≈ 25° and σ_{Ω′} ≈ 40°. For both runs of observation we get σ_{R0}/R_{0} ≈ 5 × 10^{4} and σ_{M}/M ≈ 10^{3}. Such results show that improving the astrometric accuracy or increasing the duration of the run of observation mainly allows us to obtain a better constraint on the black hole angular momentum. The results obtained for the mass and the distance of the black hole are only weakly (or not) modified. In addition, as demonstrated by Yu et al. (2016) , a better spectroscopic accuracy (e.g., σ_{V} = 1 km s^{1} instead of 10 km s^{1}) allows us to improve the constraints on both the mass and the distance of the black hole, but does not allow us to obtain better estimations of the angular momentum parameters. This is due to the fact that spectroscopic impact of the LenseThirring effect is weak (see the bottom plots in Fig. 9). In addition, we mention that the 1σ error obtained on the norm a reaches ≈0.4 when considering a run of observation of one orbital period and parameters (a,σ_{A},σ_{V}) = (0.99,10 μas, 10 km s^{1}). Such constraint thus needs to be considered since it allows us to decrease the parameter space of a for future fittings made with longer monitorings of S2 or with other S stars data. However, to strongly constrain the norm a through the observation of stellar orbits located at the Galactic center in a suitable time, it is necessary to detect new stars with GRAVITY closer to Sgr A* than S2. As mentioned by Yu et al. (2016), it should be possible to reach a 1σ uncertainty of σ_{a} ≲ 0.02 if we observe stars during a period ≲10 yr with (σ_{A},σ_{V}) = (10μas, 1−10 km s^{1}), and whose semimajor axis and eccentricity satisfy a_{sma} ≲ 40 mas and e ≳ 0.95, respectively.
The study performed here shows that, even with a stellarorbit model, without using a raytracing code but with considering approximated formulas to simulate gravitational lensing, it is possible to strongly constrain the norm of the angular momentum, the mass and the distance of the black hole, and to get non negligible constraints on the direction (i′,Ω′). In particular, we found similar results as Yu et al. (2016) which uses a model with ray tracing. However, astrometric and spectroscopic differences between models F and G are not negligible when considering mock observations obtained for (σ_{A},σ_{V}) = (10 μas, 1 km s^{1}). Indeed, as mentioned in Table 7, model F fails to describe the observations within 18 yr for p = 90%. With such spectroscopic accuracy it is thus necessary to use model G.
5. Conclusions and discussions
To conclude, the studies made in this paper show that various relativistic effects can be detected with current instruments by using astrometric and spectroscopic observations of the S2 star beginning near its next pericenter passage in 2018. In particular, transverse Doppler shift and gravitational redshift can be detected within a few months by using observations obtained with accuracies (σ_{A},σ_{V}) = (10−100 μas, 1−10 km s^{1}). If we consider the pair of accuracies (σ_{A},σ_{V}) = (10 μas, 10 km s^{1}), these effects can be detected by combining observations obtained with GRAVITY and SINFONI at VLT. Spectroscopic measurements can also be obtained with NIRSPEC at Keck. Gravitational lensing can be detected within a few years for (σ_{A},σ_{V}) = (10 μas, 10 km s^{1}). The GRAVITY and SINFONI instruments can also be used to detect this effect. It should be possible to detect the pericenter advance within a few years for (σ_{A},σ_{V}) = (10 μas, 10 km s^{1}), also by using these two instruments. The HOPC contributions can be observed within several months when considering accuracies (σ_{A},σ_{V}) = (10 μas, 1 km s^{1}). These effects can be detected using GRAVITY and MICADO. The first light of this latter instrument is expected to be in 2024, detection of the HOPC contributions will thus be obtained later than the previous effects.
Contrary to the other relativistic effects, the astrometric impact of the LenseThirring effect is maximal near apocenter passages of the S2 star, whatever the direction of the angular momentum of the black hole. The influence of this effect on the photon path is negligible which is not the case for the timelike geodesic, when considering high spin rate and some specific directions of the angular momentum. Indeed, the trajectory of the star is affected by this effect that can lead to a maximal astrometric shift of about 10 μas, 25 μas and 40 μas at the first, second and third orbital period of S2, respectively. The maximal spectroscopic shift is of about 1 km s^{1} near the three pericenter passages. Strong constraint on the angular momentum parameters of the black hole with S2 observations generated by current instruments is thus only possible using astrometric measurements (near apocenter passages) since spectrographs reaching accuracies of about 1 km s^{1} or better are not yet available. However, spectroscopic (and astrometric) observations are also important to constrain the angular momentum since the other relativistic effects are maximal near pericenter passages. More precisely, such observations will allow us to bring strong constraint on the orbital parameters, and the mass and the distance of the black hole, and thus lead to constrain the LenseThirring effect.
As null geodesics are weakly affected by the LenseThirring effect, we investigated the possibility of constraining the parameters of the angular momentum of the black hole without considering ray tracing, and by using analytical approximations of the gravitational lensing from Sereno & de Luca (2006) since such an effect can reach an amount of ≈20 μas near pericenter of the S2 star. We showed that with such a stellarorbit model it is possible to constrain the parameter space with small uncertainties. In particular, if we consider observations obtained during three S2 orbital periods and (a,Ω′,i′) = (0.99,160°,45°) with accuracies (σ_{A},σ_{V}) = (10 μas, 10 km s^{1}), we find very good 1σ errors of σ_{a} = 0.1, σ_{i′} = 15° and σ_{Ω′} = 30°. We found similar constraints even with a smaller norm of a = 0.7. If we consider σ_{A} = 30 μas or an observing time of one or two orbital periods instead of three, the constraints on the norm a are still very good: σ_{a} = 0.2−0.4. The constraint on the LenseThirring effect is thus possible with S2 observations obtained with GRAVITY and SINFONI, and by using a model that does not use ray tracing of photons, and is thus much faster. However, long monitorings are necessary to highly constrain the angular momentum parameters with S2. Detection by GRAVITY of closer stars to the Galactic center than S2 would allow us to obtain similar constraints within a few years.
We showed that the GRAVITY instrument will be capable of detecting low and highorder relativistic effects with the S2 star. This instrument will thus allow to test GR in the strongfield regime, and possibly constrain the properties of the compact source Sgr A*. We note however that the detected effects could be explained by other theories of gravitation. It is thus necessary to go further ahead in developing methods used to test GR and to investigate how it is possible to highlight deviations of those alternative theories from GR. Moreover, even if the central object is described by GR, other exotic objects also described by this theory can explain the mass at the Galactic center, such as the boson stars or the gravastar (Ruffini & Bonazzola 1969; Mazur & Mottola 2001). Comparisons between Kerr and exotic metrics have already been explored but further studies need to be performed (Grandclément et al. 2014; Sakai et al. 2014). In particular to determine whether degeneracies between the two types of objects can appear. For instance, the possibility that stellar orbits obtained in strongfield regime with a Kerr black hole and a boson star are similar has not been ruled out, even when the compact object parameters and the orbital parameters used to generate such orbits are different.
It is important to state that the work done in this paper is a preliminary study on detection of relativistic effects with the S2 star. Indeed, it neglects several contributions such as the extended mass that may be present in the Galactic center and composed of stars, stellar remnants or dark matter. It is possible that this mass modifies the star trajectory and induces a Newtonian precession which is opposite to the one due to the pericenter advance. Thus, the astrometric and spectroscopic impacts of this latter effect on S2 observations will be decreased. Several works have been done on this topic in order to determine whether this mass could be constrained and to evaluate its impact on stellar orbits observed at the Galactic center.
In particular, Rubilar & Eckart (2001) showed that the Newtonian effect can either partially or entirely compensate for the relativistic precession. The authors determined that in the particular case of the S2 star, an extended mass equal to 0.1% of the black hole mass is needed to dominate the pericenter advance effect. Weinberg et al. (2005) showed that it should be possible to constrain this extended mass if we observe the motion of 100 stars located in the 0.4 central parsec for a period of 10 yr. In their study these authors considered an extended mass of about 10^{3}M_{⊙}, and uncertainties of 0.5 mas and 10 km s^{1}. Moreover, they claimed that we will be able to detect relativistic effects such as the pericenter advance in spite of the Newtonian precession induced by this mass. This will require the use of the same conditions as used to constrain the mass. The detection of the LenseThirring effect will require consideration of astrometric accuracies less than 0.05 mas. Other studies have been performed on the detection of the LenseThirring effect in the presence of gravitational perturbations generated by stars and stellar remnants (Merritt et al. 2010; Sadeghian & Will 2011). In particular, Merritt et al. (2010) demonstrated that the detection of this effect should be possible when considering astrometric observations of GRAVITY, if the semimajor axis of the star is less than 0.5 mpc (the semimajor of S2 is around 5 mpc). In a more recent work made by Zhang & Iorio (2017), the authors studied the influence of the S102 star on S2 astrometric and spectroscopic observations. They concluded that this star will very likely obscure the angular momentuminduced effect of the black hole. However, the authors found that this effect dominates the stellar perturbations if the observed stars have semi majoraxes smaller than 0.5−2 mpc and if the black hole is maximally spinning.
All of these studies highlight the importance of taking into account the hypothetical extended mass in future stellarorbit models used to interpret the S2 star observations. If GRAVITY does not discover stars closer to the Galactic center, the detection of the different relativistic effects with S2 will be possible only if we succeed in constraining the extended mass by using observations of several stars obtained at the Galactic center, and if this mass is sufficiently weak to not dominate all relativistic effects. Besides, as mentioned in various papers (see e.g., Merritt et al. 2010; Angélil et al. 2010), we expect that the Newtonian effect can be dissociated from the relativistic effects if we consider runs of observation that include several orbital periods of the star, because those types of effects have different temporal evolutions.
See Andrae et al. (2010) for more details on how to estimate the degrees of freedom of a model.
Acknowledgments
This work was supported by ASHRA (Action Spécifique Haute Résolution Angulaire). We thank Didier Pelat for helpful discussions on statistical tests used in modelfitting processes. We also thank Guillaume Schworer for both having presented us the MCMC method used in this paper, and for advice and explanations regarding this method.
References
 Andrae, R., SchulzeHartung, T., & Melchior, P. 2010, ArXiv eprints [arXiv:1012.3754] [Google Scholar]
 Angélil, R., & Saha, P. 2010, ApJ, 711, 157 [NASA ADS] [CrossRef] [Google Scholar]
 Angélil, R., & Saha, P. 2011, ApJ, 734, L19 [NASA ADS] [CrossRef] [Google Scholar]
 Angélil, R., Saha, P., & Merritt, D. 2010, ApJ, 720, 1303 [NASA ADS] [CrossRef] [Google Scholar]
 Binder, K., & Heermann, D. W. 2002, Monte Carlo simulation in statistical physics: an introduction (New York: Springer) [Google Scholar]
 Boehle, A., Ghez, A. M., Schödel, R., et al. 2016, ApJ, 830, 17 [NASA ADS] [CrossRef] [Google Scholar]
 Bozza, V., & Mancini, L. 2004, ApJ, 611, 1045 [NASA ADS] [CrossRef] [Google Scholar]
 Bozza, V., & Mancini, L. 2005, ApJ, 627, 790 [NASA ADS] [CrossRef] [Google Scholar]
 Bozza, V., & Mancini, L. 2012, ApJ, 753, 56 [NASA ADS] [CrossRef] [Google Scholar]
 Broderick, A. E., Loeb, A., & Narayan, R. 2009, ApJ, 701, 1357 [NASA ADS] [CrossRef] [Google Scholar]
 Broderick, A. E., Fish, V. L., Doeleman, S. S., & Loeb, A. 2011, ApJ, 735, 110 [NASA ADS] [CrossRef] [Google Scholar]
 Broderick, A. E., Johannsen, T., Loeb, A., & Psaltis, D. 2014, ApJ, 784, 7 [NASA ADS] [CrossRef] [Google Scholar]
 Catanzarite, J. H. 2010, ArXiv eprints [arXiv:1008.3416] [Google Scholar]
 Doeleman, S. S., Fish, V. L., Broderick, A. E., Loeb, A., & Rogers, A. E. E. 2009, ApJ, 695, 59 [NASA ADS] [CrossRef] [Google Scholar]
 Eckart, A., & Genzel, R. 1997, MNRAS, 284, 576 [NASA ADS] [Google Scholar]
 Eisenhauer, F., Abuter, R., Bickert, K., et al. 2003, in Instrument Design and Performance for Optical/Infrared Groundbased Telescopes, eds. M. Iye, & A. F. M. Moorwood, Proc. SPIE, 4841, 1548 [Google Scholar]
 Eisenhauer, F., Perrin, G., Brandner, W., et al. 2011, The Messenger, 143, 16 [NASA ADS] [Google Scholar]
 ForemanMackey, D., Hogg, D. W., Lang, D., & Goodman, J. 2013, PASP, 125, 306 [CrossRef] [Google Scholar]
 Fragile, P. C., & Mathews, G. J. 2000, ApJ, 542, 328 [NASA ADS] [CrossRef] [Google Scholar]
 Genzel, R., Thatte, N., Krabbe, A., Kroker, H., & TacconiGarman, L. E. 1996, ApJ, 472, 153 [NASA ADS] [CrossRef] [Google Scholar]
 Genzel, R., Eisenhauer, F., & Gillessen, S. 2010, Rev. Mod. Phys., 82, 3121 [NASA ADS] [CrossRef] [Google Scholar]
 Ghez, A. M., Klein, B. L., Morris, M., & Becklin, E. E. 1998, ApJ, 509, 678 [NASA ADS] [CrossRef] [Google Scholar]
 Ghez, A. M., Salim, S., Weinberg, N. N., et al. 2008, ApJ, 689, 1044 [NASA ADS] [CrossRef] [Google Scholar]
 Gillessen, S., Eisenhauer, F., Fritz, T. K., et al. 2009a, ApJ, 707, L114 [NASA ADS] [CrossRef] [Google Scholar]
 Gillessen, S., Eisenhauer, F., Trippe, S., et al. 2009b, ApJ, 692, 1075 [NASA ADS] [CrossRef] [Google Scholar]
 Gillessen, S., Plewa, P. M., Eisenhauer, F., et al. 2017, ApJ, 837, 30 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Grandclément, P., Somé, C., & Gourgoulhon, E. 2014, Phys. Rev. D, 90, 024068 [NASA ADS] [CrossRef] [Google Scholar]
 Grould, M., Paumard, T., & Perrin, G. 2016, A&A, 591, A116 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Jaroszynski, M. 1998, Acta Astron., 48, 653 [NASA ADS] [Google Scholar]
 Johannsen, T. 2016, Class. Quant. Grav., 33, 113001 [NASA ADS] [CrossRef] [Google Scholar]
 Jørgensen, J. H., Bjælde, O. E., & Hannestad, S. 2016, MNRAS, 458, 3614 [NASA ADS] [CrossRef] [Google Scholar]
 Kannan, R., & Saha, P. 2009, ApJ, 690, 1553 [NASA ADS] [CrossRef] [Google Scholar]
 Kraniotis, G. V. 2007, Class. Quant. Grav., 24, 1775 [NASA ADS] [CrossRef] [Google Scholar]
 Levenberg, K. 1944, Quart. Applied Math., 2, 164 [Google Scholar]
 Mazur, P. O., & Mottola, E. 2001, General Relativity and Quantum Cosmology [arXiv:grqc/0109035] [Google Scholar]
 Merritt, D., Alexander, T., Mikkola, S., & Will, C. M. 2010, Phys. Rev. D, 81, 062002 [NASA ADS] [CrossRef] [Google Scholar]
 Nishiyama, S., Saida, H., Takamori, Y., et al. 2017, PASJ, submitted [Google Scholar]
 Parsa, M., Eckart, A., Shahzamanian, B., et al. 2017, ApJ, 845, 22 [NASA ADS] [CrossRef] [Google Scholar]
 Rubilar, G. F., & Eckart, A. 2001, A&A, 374, 95 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Ruffini, R., & Bonazzola, S. 1969, Phys. Rev., 187, 1767 [NASA ADS] [CrossRef] [Google Scholar]
 Sadeghian, L., & Will, C. M. 2011, Class. Quant. Grav., 28, 225029 [NASA ADS] [CrossRef] [Google Scholar]
 Sakai, N., Saida, H., & Tamaki, T. 2014, Phys. Rev. D, 90, 104013 [NASA ADS] [CrossRef] [Google Scholar]
 Sereno, M., & de Luca, F. 2006, Phys. Rev. D, 74, 123009 [NASA ADS] [CrossRef] [Google Scholar]
 Taff, L. G. 1985, Celestial mechanics: A computational guide for the practitioner (New York: WileyInterscience) [Google Scholar]
 Vincent, F. H., Paumard, T., Gourgoulhon, E., & Perrin, G. 2011, Class. Quant. Grav., 28, 225011 [NASA ADS] [CrossRef] [Google Scholar]
 Vincent, F. H., Meliani, Z., Grandclément, P., Gourgoulhon, E., & Straub, O. 2016, Class. Quant. Grav., 33, 105015 [NASA ADS] [CrossRef] [Google Scholar]
 Visser, M. 2007, ArXiv eprints [arXiv:0706.0622] [Google Scholar]
 Weinberg, N. N., Milosavljević, M., & Ghez, A. M. 2005, ApJ, 622, 878 [NASA ADS] [CrossRef] [Google Scholar]
 Will, C. M. 2008, ApJ, 674, L25 [NASA ADS] [CrossRef] [Google Scholar]
 Wollman, E. R., Geballe, T. R., Lacy, J. H., Townes, C. H., & Rank, D. M. 1977, ApJ, 218, L103 [NASA ADS] [CrossRef] [Google Scholar]
 Yu, Q., Zhang, F., & Lu, Y. 2016, ApJ, 827, 114 [NASA ADS] [CrossRef] [Google Scholar]
 Zhang, F., & Iorio, L. 2017, ApJ, 834, 198 [NASA ADS] [CrossRef] [Google Scholar]
 Zhang, F., Lu, Y., & Yu, Q. 2015, ApJ, 809, 127 [NASA ADS] [CrossRef] [Google Scholar]
 Zucker, S., Alexander, T., Gillessen, S., Eisenhauer, F., & Genzel, R. 2006, ApJ, 639, L21 [NASA ADS] [CrossRef] [Google Scholar]
Appendix A: Obtaining the astrometric positions and radial velocities of a moving star with the raytracing code Gyoto
Appendix A.1: Astrometry
Gyoto is a raytracing code integrating null geodesics backwards in time. It works on the basis that each photon is initially located at the observer screen: one pixel of the screen corresponds to the final direction of one photon. When the photon reaches the star during its integration, the pixel corresponding to this photon illuminates. It is thus possible to obtain the image of the star. However, in our study we only use one photon of the primary image and consider the final direction of this photon (α_{γ},δ_{γ}) as being the astrometric position of the star (α_{s},δ_{s}). The aim of this method is to decrease the computation time without any loss of precision on the simulated apparent position of the star. This position is obtained when the Euclidian distance between the radial coordinate of the photon and the surface of the star reaches ≈10^{2}M. The corresponding maximal astrometric and spectroscopic errors of the star are about 10^{1}μas and 10^{4} km s^{1}, respectively.
As we consider a fullGR model to generate the S2 observations, we need to consider a moving star in the raytracing code Gyoto. At a given observation date, we do not know the position of the star relative to the black hole. In other words, we cannot predict where the image of the star will be in the observer screen. We thus need to implement a procedure to get the star position (α_{s},δ_{s}). To do so, we use a mathematical function named MinDistance in Gyoto. This function gives the squared minimum Euclidian distance between the photon and the surface of the star. Zeroes of this function correspond to photons that have reached the surface of the star. The procedure that we use to converge to the closest photon from the surface of the star (≈10^{2}M) with the MinDistance function is described below:

step 1: we compute a MinDistance map considering an initialscreen with a sufficiently big fieldofview to contain the full S2orbit. A resolution of 5 × 5 pixels is taken for this first image (see the topimage in Fig. A.1). Then, we locate the pixel wherethe value of the MinDistance function is the smallest: thisminimal pixel is associated with the photon which passes theclosest to the star. This corresponds to the black pixel marked by awhite square on the top image in Fig. A.1.

Step 2: we compute another MinDistance map also considering a resolution of 5 × 5 pixels but with a fieldofview both centered on the pixel previously found and equal to the size of this pixel. This new map is thus a zoom of the former map. In this image, we again locate the minimal pixel (see the bottom image in Fig. A.1).

Step 3: we repeat the step 2 until the minimal distance between the photon and the surface of the star is of about 10^{2}M. When this distance is reached, we obtain the astrometric position (α_{s},δ_{s}) of the star which is equal to the final direction of the photon (α_{γ},δ_{γ}).
The computing time needed to reach the distance ≈10^{2}M is of about 3 s.
Fig. A.1 Maps of the logarithmic MinDistance function. Top: map obtained considering a fieldofview of 200 mas. The darker the pixels, the smaller the distance between the photon and the surface of the star. The pixel where the MinDistance function is minimal is marked by a white square. Bottom: map obtained considering the previous white square as fieldofview. The colorbars are labeled in M unit. 
Appendix A.2: Spectroscopy
The energy E of a photon measured by an observer is given by (A.1)where u and p are the fourvelocity of the observer and the fourmomentum of the photon along its geodesic, respectively. The emitted E_{em} and received E_{obs} energies of the photon are related as (A.2)with (A.3)where u_{obs(em)} and p_{obs(em)} are the fourvelocity of the observer and the fourmomentum of the photon at reception(emission), respectively. To evaluate the quantities at emission we need to know the coordinates of the star and the photon at emission. To do so, we consider the photon used to estimate the astrometric position of the star. Indeed, by using this photon we have access to its emission date and we thus can obtain the position of both the photon and the star at this date. The radial velocity of the star can also be obtained since it is related to the g factor by (A.4)where is the redshift of the star.
Fig. A.2 Illustration of the orbital parameters i, Ω and ω corresponding to the inclination of the orbit, the angle of the line of nodes and the argument of pericenter, respectively. The frame (x′,y′,z′) is a second observer frame and is related to the observer frame (α,δ,z_{obs}) defined in Sect. 2. 
Appendix B: Recovering the star coordinates in the observer coordinates (α, δ, z_{obs}) from its orbital parameters
We define a new frame (x′,y′,z′) related to the observer frame defined in Sect. 2 as: x′ = δ, y′ = α, and z′ = −z_{obs} (see Fig. A.2). Knowing the orbital parameters of the star, it is possible to get its position and velocity in the frame (x′,y′,z′). To do so, we use the usual trajectory equation of a star orbiting a central mass, and originating from the equation of motion in the twobody problem (B.1)where ν is the true anomaly corresponding to the angle between the pericenter position and the star (see Fig. B.1). This angle is obtained by using the following formula (B.2)where ℰ is the eccentric anomaly (see Fig. B.1). To evaluate this angle we need to solve the Kepler equation given by (B.3)where ℳ is the averaged anomaly (see Fig. B.1). This last angle depends on the period, the time of the pericenter passage and a date t as (B.4)The Cartesian coordinates of the star are expressed in the orbital plane (x_{orb},y_{orb}), at a given t, as By using these coordinates and the ThieleInnes formulas given by (Catanzarite 2010) we can obtain the coordinates of the star in the new frame (x′,y′,z′)Finally, the coordinates of the star in the observer frame (α,δ,z_{obs}) are given by (B.7)where (α_{s},δ_{s}) is the astrometric position of the star and v_{zs,obs} = V is its radial velocity.
Fig. B.1 Illustration of the true ν and eccentric ℰ anomalies. The black circle has a radius equal to the semimajor axis a_{sma}. The point O is the origin of the circle, F is the focus of the orbit which corresponds to the black hole in our case, and P is the pericenter. The position r(ν) of the star is also present. 
Appendix C: Resolution of the Roemer equation
In order to simulate the Roemer effect we solve the Eq. (2)by using a dichotomy method. It is resolved to within ≈10^{4} s which corresponds to a negligible astrometric shift: <10^{6}μas. It also implies an error inferior to 10^{7} km s^{1} for radial velocities.
Appendix D: Radial velocity of the star computed in models C to F
In this section, we consider G = M = c = 1. As in Appendix A.2, we consider an emitter with a fourvelocity u_{em}, emitting photons with a fourmomentum p_{em}. These photons are received by a static observer possessing a fourvelocity u_{obs}. The fourmomentum of the photons at reception is p_{obs}. As the observer is static we have . Besides, it is far from the black hole thus g_{tt}  _{obs}^{5}≈−1, so we get (D.1)which leads to . The energy of the photon as measured by this observer is thus equal to (D.2)We decide to decompose the quantities u_{em} and p_{em} in the 3 + 1 formalism of GR. To do so, we consider another observer called ZAMO (for Zero Angular Momentum Observer) with a fourvelocity u_{ZAMO}. The fourvelocity of the emitter and fourmomentum of the photon, at emission, can be decomposed into a part parallel to u_{ZAMO} (the temporal part) and a part orthogonal to u_{ZAMO} (the spatial part) as (D.3)where (D.4)is the Lorentz factor of the star as measured by the ZAMO. The unit fourvector n_{ZAMO} is the direction of emission of the photon as seen by the ZAMO, and V_{ZAMO} is the fourvelocity of the emitter as observed by the ZAMO. Finally, E_{ZAMO} is the energy of the photon at emission given by (D.5)where the fourvelocity of the ZAMO is approximated by since it is far from the black hole (the norm of the angular momentum of the black hole a tends towards zero). As seen in Appendix A.2, the total redshift is expressed as (D.6)where (D.7)Using the different previous expressions and the fact that the quantity ∂_{t}·p is conserved along the null geodesic, g becomes (D.8)Now we want to obtain an expression of the velocity V_{ZAMO} of the emitter as observed by the ZAMO. For doing so, we use the first equation of (D.3)and the approximation of the fourvelocity of the ZAMO u_{ZAMO}. We get (D.9)The Lorentz factor can be approximated as (D.10)The new components of V_{ZAMO} thus reduce to (D.11)where is the fourvelocity of the emitter in the blackhole frame with a null time coordinates. For simplicity, we nominate this velocity a threevelocity in the rest of the paper. Using the fact that the emitter is far from the black hole (a → 0), we can write the quantity g_{tt}  _{em} as (D.12)where ϵ is small compared to 1 and is equal to ϵ ≈ 2 /r_{em}, with r_{em} the radial coordinate of the emitter in the blackhole frame. Finally, the GR redshift can be expressed as (D.13)where is the velocity of the emitter projected along the lineofsight.
Appendix E: Measuring the impact of the different effects on both astrometric and spectroscopic measurements
This Appendix defines the various effects that can impact the S2 observations, and present the methods used for computing these effects.
Appendix E.1: Definition of each effect
The different effects are:

The Roemer time delay affecting both the astrometry and thespectroscopy. It is the only nonrelativistic effect of this list. It isdue to the finite speed of light propagating in a Newtonianspacetime. Depending on where the star is located along its orbit,the time needed by photons to reach the observer (followingEuclidian straight lines) will defer.

The pericenter advance affecting both the astrometry and the spectroscopy. It is due to the spacetime curvature on the star trajectory. The orbit precesses because of the gravitational field caused by the central mass (the black hole in this paper). Thus, pericenter and apocenter of the star are shifted from one period to another.

The LenseThirring effect affecting both the astrometry and the spectroscopy. It is due to the angular momentum of the black hole. This effect varies with respect to the norm and the direction of the angular momentum. The LenseThirring effect affects both the star and the photon trajectories.

The gravitational lensing effect only affecting the astrometry. It is due to the curvature of the photon geodesic that changes the apparent position of the star on the plane of sky.

The Shapiro time delay affecting both the astrometry and the spectroscopy. This is due to the slowdown of the proper time of the photon with respect to the proper time of the observer when the photon crosses a gravitational field.

The relativistic redshifts only affecting the spectroscopy. The transversal Doppler shift appears in special relativity and is due to the relative motion between the emitter and the observer. The gravitational redshift appears in GR and is due to the spacetime curvature.
Appendix E.2: Methods used to evaluate the impact of each effect on observations
The astrometric contribution of the Roemer effect is obtained by using the formula where ΔA_{α} = α_{A}−α_{B} et ΔA_{δ} = δ_{A}−δ_{B}, with (α_{A},δ_{A}) and (α_{B},δ_{B}) the apparent positions simulated by models A and B, respectively. The radial velocity contribution is obtained by using the formula ΔV_{Roemer} = V_{B}−V_{A} where V_{A} and V_{B} are radial velocities estimated with models A and B, respectively.
The astrometric and spectroscopic contributions of the pericenter advance presented in Figs. 4 and 6 are obtained by comparing models C and D. For the astrometry, we use the formula where ΔA_{α} = α_{C}−α_{D} and ΔA_{δ} = δ_{C}−δ_{D}, with (α_{C},δ_{C}) and (α_{D},δ_{D}) the apparent positions simulated by models C and D, respectively. For the spectroscopy, we compute the difference ΔV_{PA} = V_{C}−V_{D} where V_{C} and V_{D} are radial velocities estimated with models C and D, respectively.
The LenseThirring effect on astrometric positions is computed with the formula where ΔA_{α} = α_{G,a = 0}−α_{G,a = 0.99} and ΔA_{δ} = δ_{G,a = 0}−δ_{G,a = 0.99}. The index G, a = 0 and G, a = 0.99 denote positions computed considering a spin of 0 and 0.99 in model G, respectively. The radial velocity contribution of this effect is obtained with the formula ΔV_{LT} = V_{G,a = 0.99}−V_{G,a = 0} where V_{G,a = 0.99} and V_{G,a = 0} are radial velocities estimated considering a = 0.99 and a = 0 in model G, respectively.
The gravitational lensing effect is computed by using the following formula: where ΔA_{α} = α_{G, GL = 0}−α_{G,GL = 1} and ΔA_{δ} = δ_{G, GL = 0}−δ_{G, GL = 1}. The index G, GL = 0 and G, GL = 1 denote positions computed without and with gravitational lensing in model G, respectively. Let us better explain how the astrometric impact of this effect is reestimated. First we consider the star position (α_{G, GL = 1},δ_{G, GL = 1}) on the observer screen as computed by the fullGR model (i.e., our model G, see Appendix A.1). A photon is thus ray traced backward in time from the observer screen until it reaches the star at some spacetime position . Let us now consider the modified position on the observer screen (α_{G, GL = 0},δ_{G, GL = 0}) the star would have if light was propagating along Euclidian straight lines from to the observer screen. For doing so, we project the coordinates of the star at in the plane of sky by using the ThieleInnes formulas. The leading astrometric shift (α_{G, GL = 1}−α_{G, GL = 0},δ_{G, GL = 1}−δ_{G, GL = 0}) is due to both the gravitational lensing and the LenseThirring effect on the photon path. However, given that this latter effect is negligible for S2, this procedure gives access to the pure astrometric impact of the gravitational lensing effect.
To explain how we recover the astrometric impact of the Shapiro time delay let us still consider the position (α_{G, GL = 0},δ_{G, GL = 0}) as obtained following the procedure above. In this procedure, the photon reaches the observer at a time t_{obs} that takes into account the Shapiro time delay. Let us now consider the position of the star on the observer screen as computed by the modified GR model without considering GR effects on photon path (i.e. model E). In this model, the photon reaches the observer at a time slightly different from the previous t_{obs}, because at this time the photon path is an Euclidian straight line not affected by the gravitational field generated by the black hole. As a consequence, the star position at t_{obs} in model E, (α_{E},δ_{E}), will differ from (α_{G, GL = 0},δ_{G, GL = 0}). The astrometric shift (α_{G, GL = 0}−α_{E},δ_{G, GL = 0}−δ_{E}) is due to the Shapiro effect. The formula used to recover this effect is thus given by where ΔA_{α} = α_{G, GL = 0}−α_{E} and ΔA_{δ} = δ_{G,GL = 0}−δ_{E}.
The spectroscopic HOPC contributions, including the Shapiro time delay, are obtained computing the difference between radial velocities estimated with model E and those found in model G.
Finally, the transverse Doppler shift and gravitational redshift are obtained by using the formula ΔV_{TD,Grav} = V_{E,(TD,Grav) = 1}−V_{E,(TD,Grav) = 0} where V_{E,(TD,Grav) = 1} and V_{E,(TD,Grav) = 0} are radial velocities estimated with and without implementing the transverse Doppler shift and gravitational redshift in model E, respectively.
All Tables
Astrometric and spectroscopic accuracies of various current and future instruments, capable of observing in the near infrared.
Maximal astrometric offsets in μas reached at first, second and third periods of the S2 star, due to different relativistic effects and considering the black hole angular momentum parameters a = 0.99, i′ = 45°, Ω′ = 160°.
Maximal astrometric and spectroscopic differences of each model with respect to the fullGR model (model G) considering three orbital periods of the S2 star.
Estimations of the threshold times needed to detect different effects with the S2 star, considering various astrometric and spectroscopic accuracies.
Constraints obtained by fitting model F to astrometric and spectroscopic observations of the S2 star, generated with model G.
All Figures
Fig. 1 Illustration of the different reference frames: the blackhole frame (in black) (x_{bh},y_{bh},z_{bh}) labeled in KerrSchild coordinates, where the z_{bh}axis corresponds to the angular momentumaxis; the orbit frame (in blue) (x_{orb},y_{orb},z_{orb}); and finally, the observer frame (in green) (α,δ,z_{obs}). The angles i′ and Ω′ allow to recover the direction of the angular momentum of the black hole; they give the position of the blackhole frame relative to the observer frame. The point P on the orbit of the star denotes the pericenter. 

In the text 
Fig. 2 Astrometric and spectroscopic observations of the monitorings of one period (upper plots) and 6 months (lower plots) of the S2 star, simulated with the fullGR model. For one period of monitoring we consider one week of observation at pericenter passage (blue points), 6 months of observations around the pericenter passage (red points) and 14 yr of observations for the rest of the orbit (black points). The 6 months of monitoring corresponds to the first nineteen points of the oneperiod run. 

In the text 
Fig. 3 Illustration of each model described in Sect. 3.1, only focusing on the different effects affecting the astrometric measurements of the S2 star. The positions of both the star (blue line) and the photon (red line) at two dates t_{1} and t_{2} are represented for models B to G. At some particular observation time t_{obs}, the apparent position of the star on the observer screen, is highlighted by the symbols A to G. The corresponding emitting position of the star along its orbit is also represented. In model A, the difference Δt = t_{obs}−t_{em} is zero since the photon has an infinite speed. In models B to G, the Roemer time delay, noted Δt_{R}, is included, thus Δt depends on it. In model D, the star trajectory is described by GR but only accounts for the pericenter advance (PA) since the LenseThirring effect is neglected. Contrary to model D, model E takes into account this latter effect on the timelike geodesic (LTS for LenseThirring on the Star). Model F is analogous to model E but approximates the gravitational lensing effect (GL). Model G is the fullGR model; gravitational lensing is thus better estimated, the LenseThirring effect on the photon path (LTP) is naturally taken into account, and the quantity Δt depends on the Shapiro time delay, Δt_{S}, in addition to Δt_{R}. 

In the text 
Fig. 4 Astrometric impact of the pericenter advance on the S2 star observed during three orbital periods. Solid magenta, blue, and black circles correspond to the first, second, and third S2 period, respectively. Open red circles and triangles represent the apocenter and pericenter passages, respectively. 

In the text 
Fig. 5 Astrometric impact of different effects on the S2 star observations obtained during three periods. First plot: Roemer time delay. Second plot: Shapiro time delay. Third plot: LenseThirring effect considering the angular momentum parameters a = 0.99, i′ = 45° and Ω′ = 160°. Fourth plot: gravitational lensing. Open red circles and triangles denote the position of the apocenter and the pericenter, respectively. See Appendix E.2 for the precise definition of the various quantities used in these plots. 

In the text 
Fig. 6 Spectroscopic impact of the pericenter advance on the S2 star observed during three orbital periods. Solid magenta, blue, and black circles correspond to the first, second, and third S2 period, respectively. Open red circles and triangles represent the apocenter and pericenter passages, respectively. 

In the text 
Fig. 7 Spectroscopic impact of different effects on the S2 star observations obtained during three periods. First plot: Roemer time delay. Second plot: transverse Doppler shift and gravitational redshift. Third plot: LenseThirring effect considering the angular momentum parameters a = 0.99, i′ = 45° and Ω′ = 160°. Fourth plot: HOPC contributions. Open red circles and triangles denote the position of the apocenter and the pericenter, respectively. See Appendix E.2 for the precise definition of the various quantities used in these plots. 

In the text 
Fig. 8 Red curve: distribution of obtained with a model describing the observations well. The distribution follows a reduced χ^{2} law. The value (see the text for a definition of this quantity) is marked by the dashed line and is obtained considering p = 5% in Eq. (8). Blue curve: example of a distribution when the model does not describe the observations well. 90% of the are higher than the . For both curves, we consider k = 100. 

In the text 
Fig. 9 LenseThirring effect on astrometric (upper plots) and spectroscopic (lower plots) observations of the S2 star, considering a = 0.99 and various Ω′ and i′. The angular momentum direction i′ = 45° and Ω′ = 160° considered for the mock observations are visible on each plot and are denoted with the solid red circles. The different types of curves on each plot correspond to different values of the angle Ω′. For i′ = [0°,180°], we obtain the same curves whatever Ω′ for the astrometry but they are different when considering the spectroscopy: thin and thick lines correspond to i′ = 180° and i′ = 0°, respectively. On the other plots, solid, dotted, dashed and dashdotdotted curves correspond to Ω′ equal to 0°, 45°, 90° and 135°, respectively. Open circles on all plots denote the position of the apocenter and open triangles denote the pericenter passages. 

In the text 
Fig. 10 1D and 2D probability distributions obtained by fitting model F to observations during three orbital periods of the S2 star generated with model G and (σ_{A},σ_{V}) = (10μas, 10 km s^{1}). The angular momentum parameters of the black hole are (a,i′,Ω′) = (0.99,45°,160°). Blue dots and lines on each plot correspond to the parameters used to generate the mock observations.The darker the delineated area on 2D distributions, the denser the area. 

In the text 
Fig. 11 As in Fig. 10 but for a = 0.7. 

In the text 
Fig. A.1 Maps of the logarithmic MinDistance function. Top: map obtained considering a fieldofview of 200 mas. The darker the pixels, the smaller the distance between the photon and the surface of the star. The pixel where the MinDistance function is minimal is marked by a white square. Bottom: map obtained considering the previous white square as fieldofview. The colorbars are labeled in M unit. 

In the text 
Fig. A.2 Illustration of the orbital parameters i, Ω and ω corresponding to the inclination of the orbit, the angle of the line of nodes and the argument of pericenter, respectively. The frame (x′,y′,z′) is a second observer frame and is related to the observer frame (α,δ,z_{obs}) defined in Sect. 2. 

In the text 
Fig. B.1 Illustration of the true ν and eccentric ℰ anomalies. The black circle has a radius equal to the semimajor axis a_{sma}. The point O is the origin of the circle, F is the focus of the orbit which corresponds to the black hole in our case, and P is the pericenter. The position r(ν) of the star is also present. 

In the text 
Current usage metrics show cumulative count of Article Views (fulltext article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.