Testing the Black Hole No-hair Theorem with Galactic Center Stellar Orbits Hong Qi,1, 2, ∗ Richard O'Shaughnessy,3 and Patrick Brady2 1School of Physics and Astronomy, Cardiff University, Cardiff CF24 3AA, United Kingdom 2Center for Gravitation Cosmology and Astrophysics, University of Wisconsin-Milwaukee, Milwaukee, WI 53201, USA 3Center for Computational Relativity and Gravitation, Rochester Institute of Technology, Rochester, NY 14623, USA (Dated: June 22, 2021) Theoretical investigations have provided proof-of-principle calculations suggesting measurements of stellar or pulsar orbits near the Galactic Center could strongly constrain the properties of the Galactic Center black hole, local matter, and even the theory of gravity itself. As in previous studies, we use a Markov chain Monte Carlo to quantify what properties of the Galactic Center environment measurements can constrain. In this work, however, we also develop an analytic model (Fisher matrix) to understand what parameters are well-constrained and why. Using both tools, we conclude that existing astrometric measurements cannot constrain the spin of the Galactic Center black hole. Extrapolating to the precision and cadence of future experiments, we anticipate that the black hole spin can be measured with the known star S2. Our calculations show that we can measure the dimensionless black hole spin to a precision of ∼0.1 with weekly measurements of the orbit of S2 for 40 years using the GRAVITY telescope's best resolution at the Galactic Center, i.e., an angular resolution of 10 µarcsecond and a radial velocity resolution of 500 m/s. An analytic expression is derived for the measurement uncertainty of the black hole spin using Fisher matrix in terms of observation strategy, star's orbital parameters, and instrument resolution. From it we conclude that highly eccentric orbits can provide better constraints on the spin, and that an orbit with a higher eccentricity is more favorable even when the orbital period is longer. We also apply it to S62, S4711, and S4714 to show whether they can constrain the black hole spin sooner than S2. If in addition future measurements include discovery of a new, tighter stellar orbit, then future data could conceivably enable tests of strong field gravity, by directly measuring the black hole quadrupole moment. Our simulations show that with a stellar orbit similar to that of S2 but at one fifth the distance to the Galactic Center and GRAVITY's resolution limits on the Galactic Center, we can start to test the no-hair theorem with 20 years of weekly orbital measurements. I. INTRODUCTION [2] can constrain the black hole properties: its mass and particularly its spin. Specifically, we wrote a Markov The supermassive black hole at the center of our galaxy chain Monte Carlo (MCMC) code and use it to compare provides unique opportunities to investigate dynamics real and synthetic astrometric and radial velocity data near a strongly-gravitating source [1,2]. Radio tele- with models for the stellar orbits and black hole mass, scopes have imaged the immediate vicinity of the black accounting for differences in reference frame between dif- hole [3,4], allowing direct constraints on the strong grav- ferent observational campaigns. Unlike previous investi- itational field regime near the black hole via imaging ac- gations, our model includes leading order post Newtonian cretion flows [5{10]. Stellar motions also constrain the corrections to the orbit from the black hole's mass, spin, number and orbits of nearby perturbers [11]. At present, and quadrupole moment, as well as the impact of un- however, the best opportunities to constrain the Galactic known non-quadrupole internal and exterior potentials. Center come from long-term monitoring of known stars Our goal is to determine whether, despite the extremely [12{16]. These measurements can also identify effects low orbital velocity v=c ' 0:02, future measurements can from the strong gravitational field [17{19] and the prop- significantly constrain strong-field features of the Galac- erties of the supermassive black hole [2, 12{14, 20{25]. tic Center black hole. We compare our MCMC results Even stronger constraints would be possible with a well- against a detailed Fisher matrix analysis, both to validate timed pulsar orbiting the Galactic Center [26{30], at sep- our results and allow the reader to easily extrapolate to arXiv:2011.02267v2 [astro-ph.GA] 21 Jun 2021 arations comparable to a recently-discovered object [31]. future measurement scenarios. High precision inference from stellar orbits ideally should This paper is organized as follows. In Section II we account for many nearby perturbers, including the local review the observations of stellar orbits near the Galac- stellar density of visible stars [32] and compact objects tic Center; review a simplified model for stellar dynamics [22]. near supermassive black holes (justified at length in Ap- Motivated by recent discoveries of new stars in close pendixA); introduce simplified and realistic models for orbits around the Galactic Center [33], we assess how the process of measuring stellar orbits, including errors. well existing and future measurements of stellar orbits In Section III we describe two techniques to assess how well measurements can constrain properties of stellar or- bits and the supermassive black hole. The first is a sim- plified, approximate Fisher matrix. The second method ∗ [email protected] uses detailed Markov chain Monte Carlo simulations of 2 synthetic data to determine how well different param- B. Simplified models of stellar orbits eters can be measured and why. After validating our procedure using analytically tractable toy models with a The approximations involved in deriving and justify- handful of parameters, we perform full-scale simulations ing our equations of motion are provided in Appendix in SectionIV to test several hypotheses including no-hair A. Neglecting the black hole's recoil or the effect of am- theorem. Using plausible choices of parameters and fu- bient material, each star's position x evolves accord- ture achievable measurement accuracy, we discuss how ing to leading-order post-Newtonian equations of motion the black hole spin and quadrupole moment can be con- [21, 39, 40] strained with the known star S2 of an orbital period of Mx Mx M Mr_ about 16 years at about 5mpc distance from the black a = − + (4 − v2) + 4 v hole and future discoverable closer stars with orbital pe- r3 r3 r r2 2J riods as small as 1-2 years at about 1 mpc separation − [2v × J^ − 3r _^n × J^ − 3^n(L · J^)=r] [34]. We also discuss how these constraints can be af- r3 fected by an intermediate-mass black hole (IMBH) and 3 Q2 + [5^n(^n · J^)2 − 2(^n · J^)J^ − ^n]; (1) a cluster of other stars in the Galactic Center. In Sec- 2 r4 tionV we summarize the conclusions we draw from the 2 where x; v = @tx; a = @t x are the harmonic coor- studies. Throughout the paper we adopt the units where dinate position, velocity, and acceleration of the star, G = c = 1. r = jxj is the coordinate distance of the star from the black hole, ^n = x=r is a unit vector pointing towards the star, L = x × v is the orbital angular momentum, M; J;Q = −J 2=M are the mass, spin angular momen- II. STATEMENT OF THE PROBLEM 2 tum, and quadrupole moment of the black hole, and the hat over a quantity denotes its unit vector, such A. Existing Observations as J^ = J=J. Each star evolves according to a post- Newtonian Hamiltonian in [41]. There are observations of stellar orbits within 1 arcsec- For the proof-of-concept analytic calculations, we sepa- ond of the Galactic Center in infrared [13, 14, 16, 33]. In rate timescales by orbit-averaging rather than work with this paper, we are analyzing two sets of long-duration ob- the full Hamiltonian, following standard practice in ce- servations reported in Ghez et al. [13] and Gillessen et al. lestial mechanics. For analytic simplicity, we will fur- [14]. The motions of stars in the immediate vicinity of Sgr thermore treat all perturbations at leading order, there- A* have been observed in infrared bands by NTT/VLT fore performing an orbit average using a Newtonian orbit; since 1992 and by Keck telescope since 1995. The two for example, at leading order an equatorial orbit has the 2 data sets we use are the Keck data from 1995 to 2007 and form r(t) = p=(1 + e cos Φ(t)), where p = a(1 − e ) is a the VLT data from 1992 to 2009. Massive young stars semilatus rectum, a is the semimajor axis, e is the ec- are found closely orbiting the black hole at the center of centricity of the orbit, and Φ(t) is the orbital phase in our Milky Way. The locations of the stars, i.e., the astro- terms of time t. Using standard methods of celestial me- metric positions, right ascensions (RA) and declinations chanics [21, 42], we find the secular equations of motion (DEC) are recorded at different epochs. Therefore, the for the orbit average (hXi) of each star's Newtonian or- relative positions of stars to the radio source Sgr A* , i.e., bital angular momentum LN ≡ µx × v and Newtonian 2 the offsets of RA and DEC are also measured. The radial Runge-Lenz vector AN ≡ µ [v × (x × v) − GM ^n]: velocities, i.e., the line of sight components of the veloc- @ hL i = Ω~ × hL i (2) ities relative the the observers, of each star at different t N N epochs are also measured. In this work, we use the stellar @t hAN i = Ω~ × hAN i (3) orbit of star S2, because it is monitored for the longest Ω~ = Ω~ + Ω~ + Ω~ (4) time, its orbit is only 16 years, and more importantly its S J Q eccentricity is high among the few closest orbits that have ~ ^ AS ^ 3 ΩS = LN = LN 3 (5) been monitored frequently for over a decade.