A purely acceleration-based measurement of the fundamental Galactic parameters Jo Bovy∗ David A. Dunlap Department of Astronomy and Astrophysics, University of Toronto, 50 St. George Street, Toronto, Ontario, M5S 3H4, Canada (Dated: July 29, 2021) Klioner et al. have used the Gaia EDR3 data to directly measure the solar system's accelera- tion within the Milky Way using the apparent proper motions of quasars. Here I show that this single absolute acceleration measurement in combination with relative accelerations obtained from pulsar orbital decay allows one to determine all of the parameters describing the dynamics of our −1 local Galactic environment, including the circular velocity at the Sun V0 = 244 ± 8 km s and its 0 −1 −1 derivative V0 = 2 ± 9 km s kpc , the local angular frequency, the Oort constants, and the Sun's motion with respect to the LSR. This is the first determination of these parameters that only uses the general theory of relativity without the need for additional assumptions. Introduction|A good determination of the local gravi- where V0 and R0 are the local circular velocity and the tational and velocity fields near the Sun in the Milky Way distance to the Galactic center. This acceleration points is of fundamental importance to our understanding of the towards the Galactic center. local matter distribution [1]|including its dark matter Galactic gravitational fields are often expressed in content [2{4]|for correcting observations for the effects terms of the velocity of a circular orbit and the circu- of the Sun's velocity and acceleration, for predicting the lar velocity at the Sun, V0, is one of the most funda- expected annual modulation signal in dark-matter de- mental parameters of the Milky Way. It is clear that tection experiments [5], and for many other applications the acceleration measurement in Equation (1) provides in astro(-particle )physics. However, our knowledge of a stringent constraint on V0, but to actually obtain V0 the gravitational field within any galaxy, including our (and all of the other fundamental parameters that I will own, has been fundamentally limited by the fact that discuss below), we need to combine Equation (1) with a gravitational theories (i.e., that of the general theory of measurement of the Sun's distance to the Galactic cen- relativity, or its usually applicable Newtonian limit) only ter, R0. Luckily, we are fortunate enough to now have directly connect accelerations to the mass distribution, a measurement of R0 that essentially only employs fun- while we cannot measure the absolute galactic accelera- damental physics, because a highly precise measurement tion of any star. of R0 was obtained from post-Newtonian modeling of the This changed recently, because Klioner et al. [6] per- orbit of the star S2 around the Milky Way's supermassive formed a precise and, more importantly, accurate deter- black hole. Precise measurements of S2's orbit using the mination of the acceleration of the solar system barycen- GRAVITY instrument on the VLT have allowed its grav- ter with astrometric data from Gaia EDR3 [7, 8]. This itational redshift [12] and Schwarzschild precession [13] measurement uses the apparent proper-motion pattern of to be measured, allowing a precise determination of R0 distant quasars induced by the aberration effect that the [14], because R0 connects the observed motion on the sky acceleration gives rise to. This measurement agrees with with the observed line-of-sight velocity projection of the earlier determinations using Very Long Baseline Interfer- orbit. I use the value for R0 from the recent re-analysis ometry [9, 10], but is based on a set of reference sources of the GRAVITY data [15] that is orders of magnitude larger, allowing an analy- sis that is much freer of systematics. The Milky Way's R0 = 8:275 0:034 kpc ; (2) ± gravitational field is the largest contributor to the total acceleration [11] and at the current measurement pre- but note that the GRAVITY collaboration's preferred cision, plausible non-axisymmetric contributions to the value for R0 changes by more than its uncertainties be- acceleration are too small to be detected. The largest tween different analyses (a situation addressed in Ref. arXiv:2012.02169v3 [astro-ph.GA] 28 Jul 2021 non-Galactic contribution comes from the Large Magel- [15]) and disagrees with the analysis of the gravitational- lanic Cloud [6], but it is only at the few percent level and redshift of S2 by Ref. [16], who find R0 = 7:95 0:06 kpc. ± almost perpendicular to the Galactic contribution. The Together, the measurements in Equations (1) and (2) al- Klioner et al. measurement is therefore a direct mea- low the local circular velocity to be determined surement of the centripetal acceleration at the Sun, with s magnitude a0, which I express as a proper motiona ~ −1 R0 V0 = 244 8 km s : (3) 8:275 kpc a V 2 ± a~ = 0 = 0 = 5:05 0:35 µas yr−1 ; (1) c c R0 ± For the value of R0 from Ref. [16], the central value shifts −1 only to V0 = 239 km s . Another parameter that immediately follows is the an- ∗ [email protected] gular frequency of the circular orbit at the Sun, Ω0 = 2 6 ) 1 5 − 4 as yr µ ( ˜ a 3 2 R 0 )& (2013) 0 . (R V c Oort constants IAU MWPotential2014McMillan (2017) Irrgang et al Dehnen & Binney (1998) FIG. 1. Comparison of the Gaia EDR3 measurement of the local acceleration, expressed as a proper motiona ~ (see Eqn. 1), with previous indirect determinations and popular Milky Way models. V0=R0, which is constrained to be This discrepancy, if real, is an indication that the solar neighborhood as a whole may be perturbed away from a r −1 −1 8:275 kpc circular orbit by non-axisymmetric forces, at the level of Ω0 = 29:4 1:0 km s kpc : (4) 10 km s−1. ± R0 ≈ By combining the measurements of R0 and V0 in the However, we can do much more with these mea- introduction with the observed proper motion µSgr A∗ of surements if we combine them with a few other key the supermassive black hole at the center of the Milky measurements of velocities and relative accelerations in Way, Sgr A∗|assumed to be at rest at the Galactic the Milky Way and in fact determine essentially all fun- center|we can now directly determine the Sun's peculiar damental parameters describing the local gravitational velocity, because and velocity field. I explore this in this brief paper. V While the resulting measurements do not currently µSgr A∗ = Ω0 + : (5) have high precision, they stand out as being the only R0 ones that do not assume dynamical equilibrium or any Combining the latest measurement of µSgr A∗ = 6:411 other method for connecting observed velocities to the 0:008 mas yr−1 [20] with the measurements from− the in-± gravitational field. troduction, we get −1 The Sun's motion with respect to the LSR|Generally V = 8:0 8:4 km s : (6) a nuisance parameter, the Sun's peculiar velocity with re- ± spect to the circular orbit through its location (the LSR) This measurement is at the lower end of the previ- is an important parameter when studying the local ve- ously considered range and it agrees with the accepted −1 locity distribution and its gradients and is, for better or V = 12 2 km s from Ref. [18]. It is about 2σ ± −1 worse, often used as part of the correction of observed ve- off from the global value of V 25 km s from Ref. ≈ locities to the Galactocentric reference frame (this is, in [2, 19]. The uncertainty in V is largely driven by fact, unnecessary; however, it is often done). While in the the uncertainty ina ~ and it should therefore improve radial and vertical directions, the assumption of axisym- significantly with future Gaia data releases, which will metry allows the Sun's motion to be determined simply allow us to get to the root of the discrepancy between as the opposite of the mean velocity of local stars in those local and global values. directions|because on average these velocities should be zero in the Galactocentric frame|the asymmetric drift The local acceleration compared to previous constraints makes the determination of the Sun's peculiar velocity in and models|Before continuing the determination of the the direction of Galactic rotation, V , much trickier [17]. local fundamental dynamical parameters, we briefly dis- The generally accepted value derived from observations cuss how the observed acceleration compares to that −1 of local stars is V 12 km s faster than Galactic ro- derived from earlier determinations of V0 and R0 and tation [18]. However,≈ stellar velocities on larger scales also how often-used models for the Milky Way's grav- −1 are most easily understood if V 25 km s [2, 19]. itational potential predicted the observed acceleration. ≈ 3 300 ) 250 1 − 200 ) (km s Observed pulsar data R ( c Gaia EDR3 a˜ V 150 MWPotential2014 McMillan (2017) 1 1 1 100 V = 242 7 km s− + (2 9 km s− kpc− )(R R ) c ± ± − 0 7.0 7.5 8.0 8.5 9.0 R (kpc) FIG. 2. Local circular velocity curve determined solely from observed absolute accelerations for the solar system (Gaia EDR3; black point) and 6 pulsars relatively close to the Galactic mid-plane whose relative acceleration obtained from pulsar timing can be turned into an absolute determination using the solar system's acceleration. The rotation curve is approximately flat over the 2 kpc surrounding the solar radius.