# The HST Large Programme on ω Centauri. III. Absolute Proper Motion

Mattia Libralato, Andrea Bellini, Luigi Bedin, Edmundo Moreno D., José Fernández-Trincado, Barbara Pichardo, Roeland Marel, Jay Anderson, Daniel Apai, Adam J. Burgasser, et al.

► To cite this version:

Mattia Libralato, Andrea Bellini, Luigi Bedin, Edmundo Moreno D., José Fernández-Trincado, et al. The HST Large Programme on ω Centauri. III. Absolute Proper Motion. The Astrophysical Journal, 2018, 854 (1), pp.45. 10.3847/1538-4357/aaa59e. hal-02066646

HAL Id: hal-02066646

https://hal.science/hal-02066646v1

Submitted on 1 Sep 2025

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

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

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

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

# The HST Large Programme on ω Centauri. III. Absolute Proper Motion

Mattia Libralato¹,*, Andrea Bellini¹,*, Luigi R. Bedin², Edmundo Moreno D.³,*, José G. Fernández-Trincado⁴,⁵,*
Barbara Pichardo³,*, Roeland P. van der Marel¹,⁶,*, Jay Anderson¹,*, Dániel Apai⁷,⁸,*, Adam J. Burgasser⁹,*,

Anna Fabiola Marino¹⁰, Antonino P. Milone¹¹, Jon M. Rees⁷, and Laura L. Watkins¹ ¹

¹Space Telescope Science Institute, 3700 San Martin Drive, Baltimore, MD 21218, USA; libra@stsci.edu

²INAF-Osservatorio Astronomico di Padova, Vicolo dell'Osservatorio 5, Padova, I-35122, Italy

³Instituto de Astronomía, Universidad Nacional Autónoma de México, Apartado Postal 70-264, 04510-México DF, Mexico
⁴Departamento de Astronomía, Universidad de Concepción, Casilla 160-C, Concepción, Chile

⁵Institut Utinam, CNRS UMR6213, Univ. Bourgogne Franche-Comté, OSU THETA, Observatoire de Besançon, BP 1615, F-25010 Besançon Cedex, France

⁶Center for Astrophysical Sciences, Department of Physics and Astronomy, Johns Hopkins University, Baltimore, MD 21218, USA

⁷Department of Astronomy and Steward Observatory, The University of Arizona, 933 N. Cherry Avenue, Tucson, AZ 85721, USA

⁸Lunar and Planetary Laboratory, The University of Arizona, 1640 E. University Boulevard, Tucson, AZ 85721, USA

⁹Center for Astrophysics and Space Science, University of California San Diego, La Jolla, CA 92093, USA

¹⁰Research School of Astronomy and Astrophysics, Australian National University, Canberra, ACT 2611, Australia

¹¹Dipartimento di Fisica e Astronomia “Galileo Galilei”—Università di Padova, Vicolo dell’Osservatorio 3, Padova, IT-35122, Italy

*Received 2017 December 11; revised 2018 January 2; accepted 2018 January 3; published 2018 February 9*

## Abstract

In this paper, we report a new estimate of the absolute proper motion (PM) of the globular cluster NGC 5139 (ω Cen) as part of the HST large program GO-14118+14662. We analyzed a field 17 arcmin southwest of the center of ω Cen and computed PMs with epoch spans of ~15.1 years. We employed 45 background galaxies to link our relative PMs to an absolute reference-frame system. The absolute PM of the cluster in our field is $(\mu_{\alpha}\cos\delta, \mu_{\delta}) = (-3.341 \pm 0.028, -6.557 \pm 0.043)$ mas yr$^{-1}$. Upon correction for the effects of viewing perspective and the known cluster rotation, this implies that for the cluster center of mass $(\mu_{\alpha}\cos\delta, \mu_{\delta}) = (-3.238 \pm 0.028, -6.716 \pm 0.043)$ mas yr$^{-1}$. This measurement is direct and independent, has the highest random and systematic accuracy to date, and will provide an external verification for the upcoming Gaia Data Release 2. It also differs from most reported PMs for ω Cen in the literature by more than 5σ, but consistency checks compared to other recent catalogs yield excellent agreement. We computed the corresponding Galactocentric velocity, calculated the implied orbit of ω Cen in two different Galactic potentials, and compared these orbits to the orbits implied by one of the PM measurements available in the literature. We find a larger (by about 500 pc) perigalactic distance for ω Cen with our new PM measurement, suggesting a larger survival expectancy for the cluster in the Galaxy.

**Key words:** astrometry – Galaxy: kinematics and dynamics – globular clusters: individual (NGC 5139 (ω Cen)) – proper motions

## 1. Introduction

NGC 5139 (ω Cen) is one of the most complex and intriguing globular clusters (GCs) in our Galaxy. It was the first GC known to contain multiple stellar populations (Anderson 1997; Bedin et al. 2004), and its populations show increasing complexity the more we study them (e.g., Bellini et al. 2017a, 2017b, 2017d; Milone et al. 2017a, and references therein). In this context, the HST large program GO-14118+14662 (PI: Bedin), aimed at studying the white-dwarf cooling sequences of ω Cen (see Bellini et al. 2013), began observations in 2015.

In Paper I of this series (Milone et al. 2017b), we traced and chemically tagged for the first time the faint main sequences of ω Cen down to the hydrogen-burning limit. In Paper II (Bellini et al. 2018), we began to investigate the internal kinematics of the multiple stellar populations of the cluster.

In this third study of the series, we present an analysis of the absolute proper motion (PM) of ω Cen. Despite the fact that it has long been considered a GC, ω Cen has also been thought to be the remnant of a tidally disrupted satellite dwarf galaxy (e.g., Bekki & Freeman 2003). As such, understanding the orbit of this object is an important step toward understanding its true nature.

One of the key pieces of information needed to trace the orbit of an object is its absolute motion in the plane of the sky, i.e., its absolute PM. Most of the current estimates of the PM of

ω Cen were obtained with ground-based photographic plates (e.g., Dinescu et al. 1999; van Leeuwen et al. 2000, hereafter D99 and vL00, respectively), which have very different systematic issues than space-based measurements.

Thanks to the large number of well-measured background galaxies in our deep HST images (Section 2), we are able to determine a new, independent measurement of the absolute PM of ω Cen (Section 3). The PM value we computed differs by more than 5σ from the most recent estimates (Section 4). Our new measurement also implies a different orbit for the cluster. As such, we also computed orbits for ω Cen with different Galactic potentials and compared them with those obtained by employing a previous PM estimate available in the literature (Section 5).

## 2. Data Sets and Reduction

This study is focused on a 2.73 × 2.73 arcmin² field (named “F1” in Paper II) at about 17 arcmin southwest of the center of ω Cen.<fnref n="12" /> The center of the field F1 is located at $(\alpha, \delta)_{ICRS} \sim (13^h25^m37.5^s32, -47^\circ40'00''21)$. The field was observed with HST’s Ultraviolet-VISible (UVIS) channel of the Wide-Field

<footnote n="12">¹² $(\alpha, \delta)_{2000} = (13^h26^m47.5^s24, -47^\circ28'46''45)$, Anderson & van der Marel (2010).</footnote>
Camera 3 (WFC3) between 2015 and 2017 during GO-14118+14662 (PI: Bedin). Additional archival images were taken with the Wide-Field Channel (WFC) of the Advanced Camera for Surveys (ACS) in 2002 (GO-9444, PI: King) and in 2005 (GO-10101, PI: King). As such, the maximum temporal baseline covered by our observations is about 15.1 yr. The list of the data sets<fnref n="13" /> employed here is summarized in Table 1 of Paper II. Note that, as part of the large program GO-14118+14662, there are also WFC3 exposures obtained with the near-infrared (IR) channel. However, we did not use near-IR data because of the worse astrometric precision of WFC3/IR channel as compared to that of ACS/WFC and WFC3/UVIS detectors.

A detailed description of the data reduction and relative-PM computation can be found in Paper II. In the following, we provide a brief overview.

We performed our analysis on the _flt exposures corrected by the official pipeline for charge-transfer-efficiency (CTE) defects (see Anderson & Bedin 2010 for the description of this empirical correction). For each camera/filter combination, we extracted positions and fluxes of all detectable sources by point-spread function (PSF) fitting in a single finding wave without removing the neighbor stars prior to the fit. We used empirical, spatial, and time-varying PSFs obtained by perturbing the publicly available<fnref n="14" /> HST “library” PSFs (see Bellini et al. 2017a). Finally, we corrected the stellar positions for geometric distortion using the solutions of Anderson & King (2006), Bellini & Bedin (2009), and Bellini et al. (2011).

We set up a common reference-frame system (“master frame”) for each epoch/camera/filter, which is used to cross-correlate the available single-exposure catalogs. We made use of the Gaia Data Release 1 (DR1, Gaia Collaboration et al. 2016a, 2016b) catalog to (i) set a master-frame pixel scale of 40 mas pixel⁻¹ (that of WFC3/UVIS detector) and (ii) orient our system with X/Y axes parallel to east/north. Then, we iteratively cross-identified the stars in each catalog with those on the master frame by means of six-parameter linear transformations.

Following Bellini et al. (2017a), we performed a “second-pass” photometry stage by using the software KS2, an evolution of the code used to reduce the ACS Globular Cluster Treasury Survey data in Anderson et al. (2008), and we measured stellar positions and fluxes simultaneously by using individual epochs/filters/exposures at once. We set KS2 to examine F606W- and F814W-filter images to detect every possible object, and then to measure them also in all other filters. Stellar positions and fluxes were measured by refitting the PSF in the neighbor-subtracted images. Differently from Bellini et al. (2017a), we are interested not only in the stellar members of ω Cen, but also field objects (stars and galaxies) that might have moved with respect to ω Cen by more than 2 pixels (the KS2 method #1 searching radius) between two epochs. Therefore, we executed KS2 independently for the GO-9444+10101 and GO-14118+14662 data to avoid mismatching objects. Finally, our instrumental magnitudes were calibrated onto the Vega magnitude system.

The KS2 software also provides positions and fluxes for detected objects in each raw exposure. Using neighbor-source subtraction, crowded-field issues that might affect the first-pass approach were reduced. Therefore, we employed these raw,

KS2-based catalogs to measure the PMs instead of those produced by the first-pass photometry.

PMs were computed using the same methodology described in Bellini et al. (2014). The cornerstones of their iterative procedure are (i) transforming the position of each object as measured in the single catalog/epoch onto the master-frame system and (ii) fitting a straight line to these transformed positions as a function of the epoch. The slope of this line, computed after several outlier-rejection stages, is a direct measurement of the PM. We refer to Bellini et al. (2014) and Paper II for a detailed description of the PM extraction and systematic corrections.

Because the transformations that relate positions measured in one exposure to those in another were based on cluster stars, by construction, the PM of each star is measured with respect to the average cluster motion. Therefore, the distribution of ω Cen stars is centered on the origin of the vector-point diagram (VPD), while foreground/background objects are located in different parts of the VPD.

The relative PMs<fnref n="15" /> of the objects in our field are shown in Figure 1. Three groups of objects are clearly distinguishable in the VPD:

(1) cluster stars, centered in the origin of the VPD, with observed dispersion of about 0.34 mas yr⁻¹ along both Δμα cos δ and Δμδ (black points);

(2) background galaxies, clustered in a well-defined location ~7.5 mas yr⁻¹ from the mean motion of ω Cen with a dispersion of 0.40 and 0.32 mas yr⁻¹ along Δμα cos δ and Δμδ, respectively (blue circles); and

(3) field stars, characterized by a broad distribution at positive Δμδ and large negative Δμα cos δ, and a narrow tail ending at the galaxy location (gray crosses).

In Paper II, we focused on the internal kinematics of the multiple stellar populations hosted in ω Cen, so that only cluster stars with high-precision PMs were analyzed. In this paper, we focus on the global, absolute motion of the cluster. Therefore, we require high accuracy but not necessarily high precision, as was the case in Paper II. As a consequence, we make use of all objects with a PM measurement in our catalog.

### 3. Absolute PMs

The first step is to link our relative PMs to an absolute reference-frame system. This can be done in two ways: indirectly (by relying on an external catalog of absolute positions), and directly (by computing PMs with respect to fixed sources in the sky, e.g., galaxies).

For the former approach, we can cross-identify cluster members with catalogs providing absolute PMs like the Hipparcos-Tycho (Høg et al. 2000) or the UCAC5 (Zacharias et al. 2017), and compute the zero-point offset to be applied to our PMs. Unfortunately, the only stars in common with the UCAC5 catalog are saturated even in the short exposures and are not suitable for astrometry.

Alternatively, we could use the Gaia catalog as a reference. Although intriguing, the major drawback is that only 130 stars are found to be in common between our catalog and Gaia, and none of them are present in the Tycho-Gaia Astrometric

<footnote n="13">¹³ DOI reference:10.17909/T9FD49.</footnote>

<footnote n="14">¹⁴ http://www.stsci.edu/~jayander/STDPSFs/</footnote>

<footnote n="15">¹⁵ To avoid confusion between relative and absolute PMs, we use the notation (Δμα cos δ, Δμδ) for the former and (μα cos δ, μδ) for the latter.</footnote>
<figure><img src="image_3.png" /><figcaption>**Figure 1.** Left: VPD of the relative PMs of the analyzed field. Black dots (within a circle of 1.5 mas yr⁻¹ from the origin of the VPD) are members of ω Cen. Field stars are plotted with gray crosses. Objects classified as galaxies in our study are shown as blue circles. The black, dashed circle enclosing these galaxies is centered at (3.4, 6.6) mas yr⁻¹ and has a radius of 1 mas yr⁻¹. Solid (open) symbols are galaxies employed (discarded) in our estimate of the absolute PM of ω Cen (see Section 3 for detail). Right: color–magnitude diagrams of all objects in the analyzed field.</figcaption></figure>

<figure><img src="image_4.png" /><figcaption>**Figure 2.** Left: VPD of the absolute PMs in equatorial coordinates of ω Cen. Black dots and gray crosses represent ω Cen stars and field objects, respectively. The blue points (with error bars) are the galaxies actually used to compute the absolute PM. Theazure arrow indicates the motion of ω Cen stars with respect to the (fixed) galaxies, the error bars indicate our σ uncertainty in the absolute-PM value. In the inset, the green and red points (with error bars) are the absolute PMs measured using only F606W- and F814W-filter exposures, respectively (see the text for details). Right: µ<sub>α</sub>cosδ and µ<sub>δ</sub> as a function of X (bottom panels) and Y (top panels) master-frame positions. We omitted field stars for clarity. The dashed, azure stripes represent the absolute-PM ±σ error in each component. No systematic trend arises from these plots, proof of the robustness of our PMs.</figcaption></figure>

Solution (TGAS), thus preventing us from indirectly registering
our PMs in an absolute reference frame. Note that TGAS
sources would be saturated in our images, as is the case of the
objects in common with UCAC5 catalog, hence they are not
useful for the registration.

The direct way consists in estimating the absolute PM from
our data set using background galaxies in the field as a
reference (see Bedin et al. 2003, 2006; Milone et al. 2006;
Bellini et al. 2010; Massari et al. 2013; Sohn et al. 2015),
which is the method we followed.
We started by identifying by eye the center of the distribution of the many galaxy-like objects in the VPD (see Figure 1). We defined a sample as all objects within a circle of 1 mas yr⁻¹ from the preliminary center of that group (black, dashed circle in Figure 1). We then visually inspected the stacked images and rejected all fuzzy galaxies (i.e., objects without a well-defined, point-source-like core), objects close to saturated stars or their bleeding columns, and sources that could not be excluded as stars. At the end of our purging process, we are left with 45 galaxies, whose trichromatic rasters are shown in the Appendix. Given the tightness of the galaxy distribution in the VPD, it is clear that the positions of these galaxies are measured quite well, even compared to the more plentiful cluster members.

We measured the 5σ-clipped median value of the PMs in each coordinate for the surviving galaxies. This median value (with the opposite sign) represents the absolute PM of ω Cen stars at the position of our field F1 (Figure 2). We find

$$ (\mu_{\alpha} \cos \delta, \mu_{\delta})@F1 = (-3.341 \pm 0.028, -6.557 \pm 0.043) \text{ mas yr}^{-1}. \quad (1) $$

The quoted errors are the standard error to the mean.

In the right-hand panels of Figure 2, we also show the absolute PMs along the $\mu_{\alpha} \cos \delta$ and $\mu_{\delta}$ directions as a function of the X and Y positions on the master frame. We find no systematic trends in the PM distribution of either galaxies or ω Cen stars. Magnitude- and color-dependent systematic errors are also not present in our PMs as shown in Figure 3 of Paper II.

The cores of our selected galaxies are resolved by HST, and the galaxies themselves may have different flux distributions and shapes if observed through different filters. As such, the centroid of a galaxy could differ if measured with different filters, resulting in biased PM measurements. Furthermore, Bellini et al. (2011) noted that blue and red photons are refracted differently when they pass through the fused-silica window of the CCD of the WFC3/UVIS detector, even if such effect is marginal for filters redder than 3000 Å.

To assess the contributions of these effects to our absolute-PM estimate, we recomputed the relative PMs as described in Paper II in two ways: (1) using only F606W images, and (2) using only F814W images. We then measured the absolute PM of the cluster in both cases, and compared the results with our original estimate. In the inset of the left-hand panel of Figure 2, we plot the original (in azure), F606W- (in green) and F814W-based (in red) absolute PMs for the cluster. The three absolute PMs are in agreement (well within the error bars), meaning that the galaxy morphology and color effects have negligible impacts on our measurements.

Our observations were obtained in different periods of the year and, as such, our absolute-PM estimate might include an annual-parallax term. However, the contribution of the annual parallax would be of about 0.025 mas yr⁻¹ at most (by assuming a distance of ω Cen of 5.2 kpc, Harris 1996, 2010 edition), a value of the order of our absolute-PM uncertainties.

Furthermore, it is worth noting that our observations were achieved in 2002 July, 2005 December, and 2015 and 2017 August. The highest contribution of the annual parallax would take place if the PM of the galaxies (the parallax effect on ω Cen stars reflects on background objects because we computed the PMs in a relative fashion) was obtained by employing 2002–2005 or 2005–2015 (or 2017) epoch pairs.

<figure><img src="image_5.png" /><figcaption>Figure 3. Comparison between observed (top panels) and simulated (bottom panels) VPDs in Galactic coordinates (left-hand panels) and color-magnitude diagrams (right-hand panels). In panels (a) and (b), we show galaxies/field stars with gray dots and highlight the two investigated groups in gold and dark green. In panels (c) and (d), we color-code the simulated halo, and thick and thin disk stars in black, light green, and pink, respectively.</figcaption></figure>

However, the PM of 43 out of the 45 galaxies we used in our absolute-PM analysis was computed by adopting all observations at our disposal, for a total temporal baseline of 15.1 years. Therefore, the fit of the PM is well constrained by the first- and last-epoch observations (at about the same phase in the Earth orbit), and the slope of the straight line (the PM) is expected to be affected by a term much lower than 0.025 mas yr⁻¹. For all of these reasons, we expect the annual parallax to be negligible and we chose not to correct it.

### 3.1. Cluster Rotation

It is well-established that ω Cen rotates in the plane of the sky (e.g., vL00; Freeman 2001; van Leeuwen & Le Poole 2002; van de Ven et al. 2006, hereafter vdV06). We can also infer the presence of rotation directly from our VPDs. Heyl et al. (2017) showed that the skewness in the PMs of 47 Tuc is a signature of the differential rotation of the cluster, confirming previous results of the multiple-field PM analysis made by Bellini et al. (2017c). As discussed in Paper II, the PMs of ω Cen are also skewed in the tangential direction, consistent with cluster rotation.

Therefore, the PM of ω Cen we obtained in our analysis is a combination of the center of mass (COM) motion and the internal cluster rotation at the position of the field F1. The four fields observed in our HST large program of ω Cen are located at different position angles and distances from the cluster center (see Figure 1 of Paper II), and observations for the other three fields are still ongoing. When all the scheduled multiple-epoch observations are completed for all the fields, we will be able to
<figure><img src="image_6.png" /><figcaption>Figure 4. VPD of the absolute PMs in equatorial coordinates highlighting the galaxies (in blue) and the two groups of field stars analyzed in Section 3.2 (in gold and dark green colors). ω Cen stars and the remaining field objects are shown with black dots and gray crosses, respectively. The red line marks the expected reflex motion of stars at different Galactocentric distances (see the text for details).</figcaption></figure>

make a precise measurement of the rotation of the cluster. Here, we use the kinematic analyses of vdV06, in which they adopted the PM catalog of vL00 to derive the rotation of ω Cen.

Our field is 17 arcmin away from the center of ω Cen, and projection effects (because of the different line of sight between the COM and our field F1) can produce an additional perspective PM component. In order to measure the PM of the COM of ω Cen, we must also remove the contribution of the projection from our estimate.

The size of our field is small as compared to its distance from the center of ω Cen. Therefore, we assumed that, at the first order, the effects of perspective and internal rotation are the same for all stars in field F1. Because of this, we corrected the absolute-PM measurement we inferred in Section 3 and not the PM of each individual star.

First, we computed the contribution of perspective using the prescription given in van der Marel et al. (2002). By adopting as the reference position of our field the average position of all its stars, we find

$$ (\mu_{\alpha} \cos \delta, \mu_{\delta})_{\text{Perspective}} = (0.008, 0.043) \text{mas yr}^{-1}. \quad (2) $$

We then assessed the contribution of the cluster rotation. Table 3 of vdV06 reports the mean PM of ω Cen in different polar sectors in the plane of the sky. These velocities size the amount of the rotation in different locations of the cluster. After we identified in which polar sector our field lies, we transformed the corresponding PM values in Table 3 of vdV06 from the axisymmetric reference frame adopted in their study to the equatorial reference frame, and decomposed these PMs in the corresponding radial and tangential components. Since we are only interested in the amount of rotation, we just need to correct our PMs by the tangential component. Finally, we transformed these PM values back to the equatorial

coordinate system. The expected contribution of cluster rotation to the absolute PM is

$$ (\mu_{\alpha} \cos \delta, \mu_{\delta})_{\text{Rot}} = (-0.111, 0.116) \text{mas yr}^{-1}. \quad (3) $$

By subtracting the perspective and the rotation components to our absolute PM at the position of our field F1, we obtain our best estimate of the absolute PM of the COM of ω Cen:

$$ \begin{aligned} (\mu_{\alpha} \cos \delta, \mu_{\delta})_{\text{COM}} \\ = (-3.238 \pm 0.028, -6.716 \pm 0.043) \text{mas yr}^{-1}. \end{aligned} \quad (4) $$

The difference between the absolute PM at the position of our field F1 and that of the COM is of about 0.19 mas yr⁻¹, dominated by cluster rotation. This correction is qualitatively in agreement with that expected by Figure 6 of vdV06, both in magnitude and direction.

It is important to note that vdV06 employed the PM catalog of vL00, which is known to have severe color- and magnitude-dependent systematic effects as described by Platais et al. (2003). The uncertainty in the COM PM listed above includes only random errors, and no systematic uncertainties in the rotation correction. If the ~4 km s⁻¹ PM rotation measured by vdV06 is inaccurate by more than ~30%, then this systematic error dominates the error budget.

To determine the corresponding Galactocentric velocity of the COM, we adopt a Cartesian Galactocentric coordinate system (X, Y, Z), with the origin at the Galactic center, the X-axis pointing in the direction from the Sun to the Galactic center, the Y-axis pointing in the direction of the Sun's Galactic rotation, and the Z-axis pointing toward the Galactic North Pole. We use the same quantities and corresponding uncertainties for the Sun and for ω Cen as in Section 5 below (see Tables 2 and 3). This implies a Galactocentric (X, Y, Z) position

$$ \mathbf{r} = (-5.1, -3.9, 1.3) \text{kpc}, \quad (5) $$

and a Galactocentric velocity vector

$$ \mathbf{v} = (97.9 \pm 3.0, -18.0 \pm 8.6, -80.2 \pm 7.2) \text{km s}^{-1}. \quad (6) $$

The corresponding Galactocentric radial and tangential velocities are

$$ (V_{\text{rad}}, V_{\text{tan}}) = (-82.0 \pm 7.8, 98.1 \pm 10.4) \text{km s}^{-1}, \quad (7) $$

and the observed total velocity with respect to the Milky Way is

$$ \mathbf{v} \equiv |\mathbf{v}| = (127.9 \pm 3.5) \text{km s}^{-1}. \quad (8) $$

The listed uncertainties were obtained from a Monte-Carlo scheme that propagates all observational distance and velocity uncertainties and their correlations, including those for the Sun.

### 3.2. Field Stars: Who is Who

The VPDs in Figures 1 and 2 reveal the presence of two main field-star components. In panel (a) of Figure 3, we highlighted them with gold and dark green crosses. These two groups of stars also show different loci in the color-magnitude diagrams (panel b). To identify the stellar populations they belong to, we compared the observed absolute PMs of field stars in our catalog to those of the Besançon models (Robin et al. 2003). We simulated stars out to 50 kpc in a 6 × 6 arcmin² region centered on our field. We assumed the photometric and PM errors to be zero since we are only interested in tagging each component of the field. The final
<table><caption>**Table 1**
List of the Absolute PM of ω Cen in the Literature and Computed in This Paper</caption><thead><tr><th>ID</th><th>References</th><th>μ<sub>α</sub> cos δ<br>(mas yr<sup>-1</sup>)</th><th>μ<sub>δ</sub><br>(mas yr<sup>-1</sup>)</th><th>Extragalactic?</th><th>Notes</th></tr></thead><tbody><tr><td>M65</td><td>Murray et al. (1965)</td><td>+0.5 ± 0.6</td><td>-7.7 ± 0.5</td><td>No</td><td></td></tr><tr><td>W66</td><td>Woolley (1966)</td><td>-3.6 ± 0.5</td><td>-6.0 ± 0.5</td><td>No</td><td>(a)</td></tr><tr><td rowspan="2">D99</td><td rowspan="2">Dinescu et al. (1999)</td><td>-4.88 ± 0.35</td><td>-3.47 ± 0.34</td><td>Yes</td><td rowspan="2">(b)</td></tr><tr><td>-5.08 ± 0.35</td><td>-3.57 ± 0.34</td><td>Yes</td></tr><tr><td>vL00</td><td>van Leeuwen et al. (2000)</td><td>-3.97 ± 0.41</td><td>-4.38 ± 0.41</td><td>No</td><td>(c)</td></tr><tr><td rowspan="2">G02</td><td rowspan="2">Geffert et al. (2002)</td><td>-4.2 ± 0.5</td><td>-5.1 ± 0.5</td><td>No</td><td rowspan="2">(d)</td></tr><tr><td>-4.1 ± 0.5</td><td>-4.8 ± 0.5</td><td>No</td></tr><tr><td>K13</td><td>Kharchenko et al. (2013)</td><td>-6.01 ± 0.25</td><td>-5.02 ± 0.25</td><td>No</td><td>(e)</td></tr><tr><td>vL00 (@F1)</td><td>This paper</td><td>-4.17 ± 0.41</td><td>-4.31 ± 0.41</td><td>No</td><td>(f)</td></tr><tr><td>vL00 (COM)</td><td>This paper</td><td>-4.06 ± 0.41</td><td>-4.43 ± 0.41</td><td>No</td><td>(g)</td></tr><tr><td>HSOY (@F1)</td><td>This paper</td><td>-4.08 ± 0.44</td><td>-6.68 ± 0.26</td><td>No</td><td></td></tr><tr><td>HSOY (All)</td><td>This paper</td><td>-3.21 ± 0.15</td><td>-6.51 ± 0.14</td><td>No</td><td>(h)</td></tr><tr><td>UCAC5 (@F1)</td><td>This paper</td><td>-4.10 ± 0.27</td><td>-7.00 ± 0.26</td><td>No</td><td></td></tr><tr><td>UCAC5 (All)</td><td>This paper</td><td>-3.30 ± 0.03</td><td>-6.80 ± 0.03</td><td>No</td><td>(h)</td></tr><tr><td>L18 (@F1)</td><td>This paper</td><td>-3.341 ± 0.028</td><td>-6.557 ± 0.043</td><td>Yes</td><td></td></tr><tr><td>L18 (COM)</td><td>This paper</td><td>-3.238 ± 0.028</td><td>-6.716 ± 0.043</td><td>Yes</td><td>(i)</td></tr></tbody></table>

**Note.** (a): Obtained by Geffert et al. (2002) by linking Murray et al. (1965) and Woolley (1966) catalog on Tycho 2 catalog. (b): The two values corresponds to the absolute PM with and without the correction for the field gradient of the mean motion of the reference stars, respectively. (c): Absolute-PM registration on Tycho 2 catalog. (d): Derived using UCAC1 (Zacharias et al. 2000) PMs. (e): Derived using PPMXL catalog (Roeser et al. 2010). (f): Obtained from van Leeuwen et al. (2000) catalog after correcting the PMs for the perspective and solid-body rotations. The errors are the original uncertainties of van Leeuwen et al. (2000) because representation of the sources of systematics in their PMs. (g): Obtained from van Leeuwen et al. (2000) catalog after removing perspective-, solid-body-, and cluster-rotation contributions. The errors are the same of vL00 (@F1). (h): PM estimate obtained considering all likely cluster members in the catalog, not only those at the position of our field F1. (i): Obtained after the correction of the perspective and cluster rotations.

simulated sample is composed of halo (shown with black crosses in the bottom panels of Figure 3), thin- (pink crosses), and thick-disk (light green crosses) stars.

Both color-magnitude diagrams and VPDs (in Galactic coordinates) in Figure 3 suggest that the broad group of field stars in our observed VPD is mainly constituted by disk stars, while the narrow stream connecting the distributions of disk stars to the galaxies is associated with halo stars.

En passant, we find broad agreement between observed and simulated PMs, an independent validation of the theoretical predictions of Besançon models.

We can reach similar conclusions by considering that, when we look at the stars along the line of sight toward ω Cen, we see reflex motion due to the motion of the Sun. By assuming a Sun motion in Galactocentric reference frame given in Section 5 below (see Table 2), we expect a reflex proper motion of

$$ \left\{ \begin{array}{l} \mu_{\alpha} \cos \delta = -6.391 \times (5.2/d) \text{ mas yr}^{-1} \\ \mu_{\delta} = -3.190 \times (5.2/d) \text{ mas yr}^{-1}, \end{array} \right. \qquad (9) $$

where *d* is the distance of a given star from the Galactic center in kiloparsecs. For hypothetical stars that stand still in the Galactocentric reference frame, this motion is equal to the observed PM. In the VPD in Figure 4, we used the equations above to mark the expected location of a given object at different Galactocentric distances.

Galaxies and very-distant halo stars have no motions in the Galactocentric reference frame. Therefore, they fall at the large distance end of the red line with a small scatter (because the PMs induced by their peculiar velocities at large distances are small). Intermediate-distance stars are in the inner disk and bulge, and these stars have an orbital motion in the Milky Way.

As such, they do not fall exactly along the line. Finally, nearby stars have a large scatter about the line because of their large PMs induced by peculiar velocities, even if the latter are relatively small.

Looking at Figure 4, it seems that gold points are much closer to the straight line and at large distances, as expected by halo and intermediate-disk stars. Green objects have instead a broad distribution far from the red line, suggesting that they are nearby disk stars.

## 4. Comparison with the Literature

The current estimates of the absolute PMs of ω Cen available in the literature are collected in Table 1 and are shown in the left-hand panel of Figure 5. With the interesting exception of the value of Geffert et al. (2002) obtained by linking the PM catalog of Murray et al. (1965) and Woolley (1966) to the Tycho 2 catalog (dark blue point with label W66), the disagreement between our result and the literature values is clear. For example, there is a difference of more than 2 mas yr⁻¹ in modulus and about 30° and 15° in orientation, resulting in 9σ and 5σ discrepancies, between our value and those of D99 and vL00, respectively.

To understand whether our methodology or the analyzed field can be at the basis of these discrepancies, as a test, we downloaded the publicly available vL00 catalog and computed the absolute PM of ω Cen using PMs therein. We limited our analysis to 45 cluster members ([stars] with a membership probability greater than 85% in the vL00 catalog) at the position of our field F1.

vdV06 pointed out that the PMs of vL00 are affected by perspective and solid-body (due to a misalignment between the
<figure><img src="image_7.png" /><figcaption>**Figure 5.** Left: comparison between the absolute PM of ω Cen obtained with HST data and the literature. The blue and red points represent the absolute PM at the distance of our field (L18 @F1) and of the COM (L18 COM, obtained with the corrections described in Section 3.1), respectively. Using the IDs listed in Table 1, we plot the values of W66 in dark blue, D99 in turquoise, vL00 in magenta, G02 in gray, and K13 in cyan. The value of M65 is shown as a light green point in the inset. The value of ω Cen COM absolute PM obtained from the corrected vL00 PMs (see the text for details) is shown in purple. The value obtained using stars in the vL00 catalog (with PMs corrected only for perspective and solid-body rotations) at the same position of our F1 field (in pink) is plotted as reference. Right: absolute-PM of ω Cen at the position of our field F1 on the same scale obtained by using UCAC5 (in orange) and HSOY (in brown) catalogs. The error bars are the standard errors to the mean. The orange and brown triangles represent the absolute PM of ω Cen computed using the PMs of all likely cluster members (not only those in the field F1) in UCAC5 and HSOY catalogs, respectively. In the inset, we zoom-in around the location of our absolute-PM estimates.</figcaption></figure>

<table><caption>**Table 2**
Parameters Employed for the Galactic Potential</caption><thead><tr><th>Parameter</th><th>Value</th><th>References</th></tr></thead><tbody><tr><td>R<sub>0</sub></td><td>8.30 ± 0.23 kpc</td><td>(a)</td></tr><tr><td>Θ<sub>0</sub></td><td>239 ± 7 km s<sup>-1</sup></td><td>(a)</td></tr><tr><td>(U, V, W)<sub>☉</sub></td><td>(-11.10 ± 1.20, 12.24 ± 2.10, 7.25 ± 0.60) km s<sup>-1</sup></td><td>(a, b)</td></tr><tr><td colspan="3"><em>Galactic Bar</em></td></tr><tr><td>Mass</td><td>1.1 × 10<sup>10</sup> M<sub>☉</sub></td><td>(c, d)</td></tr><tr><td>Present position of major axis</td><td>20°</td><td>(e)</td></tr><tr><td>Ω<sub>bar</sub></td><td>40, 45, 50 km s<sup>-1</sup> kpc<sup>-1</sup></td><td>(f, g)</td></tr><tr><td>Cut-off radius</td><td>3.28 kpc</td><td>(d, h, i)</td></tr><tr><td colspan="3"><em>Spiral Arms</em></td></tr><tr><td>M<sub>arms</sub>/M<sub>disk</sub></td><td>0.05</td><td>(l)</td></tr><tr><td>Scale length</td><td>5.0 kpc</td><td>(m)</td></tr><tr><td>Pitch angle</td><td>15.5°</td><td>(n)</td></tr><tr><td>Angular velocity</td><td>25 km s<sup>-1</sup> kpc<sup>-1</sup></td><td>(m, o)</td></tr></tbody></table>

**References.** (a): Brunthaler et al. (2011); (b): Schönrich et al. (2010); (c): Robin et al. (2012), (d): Fernández-Trincado et al. (2017b); (e): Gerhard (2002); (f): Bissantz et al. (2003); (g): Sormani et al. (2015); (h): Fernández-Trincado et al. (2017a); (i): Portail et al. (2015); (l): Pichardo et al. (2012); (m): Martinez-Medina et al. (2017); (n): Drimmel (2000); (o): Gerhard (2011).

photographic plates of different epochs) rotations. Therefore, we corrected the perspective component as described in Section 3.1; and the solid-body part as described in vdV06 (Equation (7)), using the average position of the 45 stars as reference, and transforming the corresponding PM correction from the ω Cen axisymmetric reference frame to the equatorial

reference frame. By applying these two corrections, we computed the absolute PM at the position of our field F1 (pink point in the left-hand panel of Figure 5). Finally, we removed the effect of the cluster rotation as described in Section 3.1 and obtained the absolute PM of ω Cen COM. The result is shown in the left-hand panel of Figure 5 as a purple point. Again, the difference with our absolute-PM estimates is about 5σ.

We also measured the absolute PM of ω Cen using the absolute PMs published in the UCAC5 and “Hot Stuff for One Year” (HSOY, Altmann et al. 2017) catalogs. We chose these two catalogs because both of them employ the Gaia DR1 catalog to link their astrometry to an absolute reference-frame system. The absolute-PM measurements were obtained by using 45 cluster members, selected on the basis of their location in the CMD and in the VPD, at the same location of the stars analyzed in our paper. In the right-hand panel of Figure 5 we show the absolute PM of ω Cen at the distance of our field F1 (i.e., without considering the effect of the perspective and cluster rotation) for UCAC5 (orange point) and HSOY (brown point) catalogs. Both UCAC5 and HSOY values are closer to our measurements than to those in the literature. For completeness, we also computed the PM of ω Cen using all likely cluster members in the UCAC5 and HSOY catalogs, not only those at the position of our field F1. By considering a more populated sample in the computation (2394 for UCAC5 and 339 for HSOY), the agreement with our estimate becomes tighter (orange and brown triangles for UCAC5 and HSOY, respectively, in the right-hand panel of Figure 5). Note that the errors for UCAC5 value are smaller than those for our estimate because of the larger number of considered objects. This external check further strengthens our result.
<figure><img src="image_8.png" /><figcaption>**Figure 6.** Meridional orbits of ω Cen for the last 2 Gyr for our first Galactic potential, with an angular velocity gradient of the Galactic bar of 40 (top row), 45 (middle row), 50 (bottom row) km s⁻¹ kpc⁻¹. The orbits in the first column (red lines) employed the PM value of D99, while the orbits in the second and third columns (black lines) adopted those computed in this paper. In the third column, the orbits were computed by adding the spiral arms with an angular velocity of 25 km s⁻¹ kpc⁻¹.</figcaption></figure>

## 5. The Orbit of ω Cen

In order to analyze the repercussions of the new PM determinations on the orbit of ω Cen, we computed the orbit in two detailed mass models based on the Milky Way Galaxy.

The first model is based on an escalation of the axisymmetric Galactic model of Allen & Santillan (1991), for which we converted its bulge component into a prolate bar with the mass distribution given in Pichardo et al. (2004). We also included the 3D spiral-arm model considered in Pichardo et al. (2003) to explore any possible effect of the spiral arms on ω Cen. For this first potential, we computed the minimum and maximum distances from the Galactic center, the maximum vertical distance from the Galactic plane, the orbital eccentricity, and the angular momentum. The eccentricity *e* is defined as

$$e = (R_{\text{max}} - R_{\text{min}})/(R_{\text{max}} + R_{\text{min}}), \quad (10)$$

where $R_{\text{min}}$ and $R_{\text{max}}$ are successive minimum and maximum distances from the Galactic rotation axis (i.e., *R* is the distance in cylindrical coordinates), respectively.

From a gravitational perspective, the most influential part of the Galaxy for dynamically hot (eccentric) stellar clusters or dwarf galaxies is the inner region, where the bar and the maximum mass concentration of the Galaxy are located. An accurate model of the central part of the Galaxy is paramount for an accurate determination of the orbit. For these reasons, we performed a second experiment using the new galaxy-modeling algorithm called *GravPot16*<fnref n="16" /> (Fernández-Trincado et al. 2017b; Fernández-Trincado et al. 2018, in preparation) that specifically attempts to model the inner Galactic region. *GravPot16* is a semi-analytic, steady-state, 3D gravitational potential of the Milky Way observationally and dynamically constrained. *GravPot16* is a

<figure><img src="image_9.png" /><figcaption>**Figure 7.** Same as Figure 6 but with orbits calculated assuming the second Galactic potential.</figcaption></figure>

<table><caption>**Table 3**

ω Cen Parameters Adopted in the Orbit Computation</caption><thead><tr><th>Parameter</th><th>Value</th><th>Reference</th></tr></thead><tbody><tr><td>(α, δ)<sub>J2000</sub></td><td>(13h26m47s 24, -47°28′46″45)</td><td>(a)</td></tr><tr><td>Radial Velocity</td><td>232.2 ± 0.7 km s<sup>-1</sup></td><td>(b, c)</td></tr><tr><td>Distance</td><td>5.20 ± 0.25 kpc</td><td>(d)</td></tr><tr><td>X</td><td>-5.1 kpc</td><td></td></tr><tr><td>Y</td><td>-3.9 kpc</td><td></td></tr><tr><td>Z</td><td>1.3 kpc</td><td></td></tr><tr><td colspan="3" style="text-align: center;">Dinescu et al. (1999)</td></tr><tr><td>μ<sub>α</sub> cos δ</td><td>(-5.08 ± 0.35) mas yr<sup>-1</sup></td><td></td></tr><tr><td>μ<sub>δ</sub></td><td>(-3.57 ± 0.34) mas yr<sup>-1</sup></td><td></td></tr><tr><td>ν<sub>X</sub></td><td>58.0 km s<sup>-1</sup></td><td></td></tr><tr><td>ν<sub>Y</sub></td><td>-22.8 km s<sup>-1</sup></td><td></td></tr><tr><td>ν<sub>Z</sub></td><td>0.2 km s<sup>-1</sup></td><td></td></tr><tr><td colspan="3" style="text-align: center;">Our paper (L18 @F1)<sup>a</sup></td></tr><tr><td>μ<sub>α</sub> cos δ</td><td>(-3.341 ± 0.028) mas yr<sup>-1</sup></td><td></td></tr><tr><td>μ<sub>δ</sub></td><td>(-6.557 ± 0.043) mas yr<sup>-1</sup></td><td></td></tr><tr><td>ν<sub>X</sub></td><td>95.7 km s<sup>-1</sup></td><td></td></tr><tr><td>ν<sub>Y</sub></td><td>-18.4 km s<sup>-1</sup></td><td></td></tr><tr><td>ν<sub>Z</sub></td><td>-76.1 km s<sup>-1</sup></td><td></td></tr></tbody></table>

**Note.** We adopted a distance error of ~5%, which is a reasonable guess in agreement with the most recent estimates for ω Cen distance by using RR Lyrae stars (e.g., Fernández-Trincado et al. 2015; Navarrete et al. 2017).

a We adopted the F1 PM value instead of that of the COM to remove any dependence on the correction for the PM rotation of ω Cen.

**References.** (a): Anderson & Bedin (2010); (b): Mayor et al. (1997); (c): Reijns et al. (2006); (d): Harris (1996), 2010 edition.

non-axisymmetric potential including only the bar in its non-axisymmetric components. The potential model is primarily the superposition of several composite stellar components where the density profiles in cylindrical coordinates, $ρ_i(R, z)$, are the same as those proposed in Robin et al. (2003, 2012, 2014), i.e., a

<footnote n="16">¹⁶ https://fernandez-trincado.github.io/GravPot16/</footnote>
<table><caption>**Table 4**

Monte-Carlo, Average Orbital Parameters (and Corresponding $1\sigma$ Uncertainties) of $\omega$ Cen for Our First Galactic Potential, Using D99 and Our Absolute-PM (@F1) Values</caption><thead><tr><th>Ω<sub>bar</sub><br>(km s<sup>-1</sup> kpc<sup>-1</sup>)</th><th>⟨r<sub>min</sub>⟩<br>(kpc)</th><th>⟨r<sub>max</sub>⟩<br>(kpc)</th><th>⟨|z|<sub>max</sub>⟩<br>(kpc)</th><th>⟨e⟩</th><th>⟨h⟩<br>(km s<sup>-1</sup> kpc)</th></tr></thead><tbody><tr><td colspan="6" style="text-align: center;">Dinescu et al. (1999)</td></tr><tr><td>40</td><td>0.78 ± 0.16</td><td>7.39 ± 0.35</td><td>3.86 ± 0.28</td><td>0.76 ± 0.05</td><td>-319 ± 54</td></tr><tr><td></td><td>0.83 ± 0.23</td><td>7.17 ± 0.45</td><td>3.78 ± 0.40</td><td>0.71 ± 0.07</td><td>-351 ± 64</td></tr><tr><td>45</td><td>0.83 ± 0.12</td><td>7.12 ± 0.32</td><td>3.81 ± 0.69</td><td>0.72 ± 0.05</td><td>-341 ± 41</td></tr><tr><td></td><td>0.84 ± 0.18</td><td>7.08 ± 0.26</td><td>3.70 ± 0.63</td><td>0.72 ± 0.05</td><td>-346 ± 49</td></tr><tr><td>50</td><td>0.91 ± 0.13</td><td>6.86 ± 0.29</td><td>3.67 ± 0.52</td><td>0.69 ± 0.04</td><td>-373 ± 39</td></tr><tr><td></td><td>0.92 ± 0.15</td><td>6.83 ± 0.24</td><td>3.63 ± 0.60</td><td>0.69 ± 0.04</td><td>-368 ± 41</td></tr><tr><td colspan="6" style="text-align: center;">Our paper (L18 @F1)</td></tr><tr><td>40</td><td>1.39 ± 0.14</td><td>7.21 ± 0.20</td><td>3.52 ± 0.59</td><td>0.61 ± 0.03</td><td>-513 ± 42</td></tr><tr><td></td><td>1.45 ± 0.19</td><td>7.13 ± 0.16</td><td>3.39 ± 0.66</td><td>0.57 ± 0.04</td><td>-538 ± 56</td></tr><tr><td>45</td><td>1.23 ± 0.16</td><td>7.53 ± 0.58</td><td>3.65 ± 0.78</td><td>0.64 ± 0.05</td><td>-495 ± 36</td></tr><tr><td></td><td>1.48 ± 0.14</td><td>7.06 ± 0.13</td><td>2.94 ± 0.82</td><td>0.57 ± 0.04</td><td>-565 ± 50</td></tr><tr><td>50</td><td>1.32 ± 0.16</td><td>7.40 ± 0.31</td><td>4.26 ± 0.55</td><td>0.62 ± 0.04</td><td>-475 ± 49</td></tr><tr><td></td><td>1.36 ± 0.20</td><td>7.30 ± 0.36</td><td>4.03 ± 0.76</td><td>0.59 ± 0.04</td><td>-494 ± 58</td></tr></tbody></table>

**Note.** For each angular velocity gradient of the galactic bar, the first and the second rows represent the potential with only the bar and with bar + spiral arms, respectively. From left to right, we provide the angular velocity of the galactic bar ($\Omega_{bar}$), the minimum ($\langle r_{min} \rangle$) and maximum ($\langle r_{max} \rangle$) distances to the galactic center, the maximum distance from the galactic plane ($\langle |z|_{max} \rangle$, the orbital eccentricity ($\langle e \rangle$), and the angular momentum ($\langle h \rangle$).

boxy/peanut bulge, a Hernquist stellar halo, seven stellar Einasto thin disks with spherical symmetry in the inner regions, two stellar sech² thick disks, a gaseous exponential disk, and a spherical structure associated with the dark matter halo. A new formulation for the global potential of this Milky Way density model will be described in detail in a forthcoming paper (J. G. Fernández-Trincado et al. 2018, in preparation). For the second potential, we computed the orbits in a Galactic non-inertial reference frame in which the bar is at rest, and we specifically examined the orbital Jacobi constant per unit mass, $E_j$, which is conserved in this reference frame. The orbital energy $E_j$ per unit of mass is defined as

$$E_j = \frac{1}{2}v^2 + \Phi_{\text{axi}} + \Phi_{\text{non-axi}} - \frac{1}{2}|\Omega_{\text{bar}} \times \mathbf{R}|^2, \quad (11)$$

with $\Phi_{\text{axi}}$ and $\Phi_{\text{non-axi}}$ the gravitational potential of the axisymmetric and non-axisymmetric (in this case the bar) Galactic potential, respectively; $\Omega_{\text{bar}}$ the angular velocity gradient, and $\mathbf{R}$ the vector $\mathbf{r}$ of the particles/stars.

For the potential of the bar, we assumed (i) a total mass of $1.1 \times 10^{10} M_{\odot}$, (ii) an angle of 20° for the present-day orientation of the major axis of the bar, (iii) an angular velocity gradient of the bar $\Omega_{\text{bar}}$ of 40, 45, and 50 km s⁻¹ kpc⁻¹, and (iv) a cut-off radius for the bar of $R_{\text{cut}} = 3.28$ kpc. The bar potential model was computed using a new mathematical technique that considers ellipsoidal shells with similar linear density (J. G. Fernández-Trincado et al. 2018, in preparation). Table 2 summarizes the parameters employed for the Galactic potential, the position and velocity of the Sun in the Local Standard of Rest, and the properties of the bar and of the spiral arms.

Hereafter, we discuss the results achieved by employing the absolute-PM value measured for the stars in field F1. We also ran

the corresponding calculations for the estimated COM PM and found a negligible difference (e.g., comparable perigalacticons, apogalacticons, region covered by the orbit) as expected given the small PM difference between the two estimates. The COM is physically the correct quantity to use, but by using the value for F1, we remove any dependence on the assumed, and possibly uncertain, correction for the PM rotation of the cluster (and note that the perspective and rotation corrections for $\omega$ Cen are smaller than the uncertainty in the azimuthal velocity of the Sun). The other $\omega$ Cen parameters are listed in Table 3. Finally, we also calculated, for comparison, the $\omega$ Cen orbit by assuming the PM value of D99.

Figures 6 and 7 present the orbits of $\omega$ Cen obtained by considering the first and the second Galactic potential, respectively. For each value of the angular velocity gradient of the Galactic bar (different rows) we investigated, we show the meridional orbit of $\omega$ Cen obtained using the absolute-PM value of D99 (red lines in the left columns) and that computed with HST data in our paper (black lines in the other columns). Furthermore, in the rightmost column of Figure 6, we display the orbit calculated by including the spiral arms in the computation.

Looking at Figure 6, some differences can be discerned between the different orbits. The radial excursions of the orbit of $\omega$ Cen seem to increase slightly with the new PM value, while the vertical excursions remain similar regardless of the PM adopted. The main difference between the new and the literature PM values is the larger (about half a kpc) perigalactic distance obtained with the PM computed in our paper. This more distant approach to the innermost part of the Galaxy likely means a larger survival expectancy for the cluster.

Examining the impact of the spiral arms, contrary to intuition, we find that the spiral arms can influence the orbit, although slightly. However, the differences with this and the bar-only potential model is statistical and does not change the
<table><caption>**Table 5**
Same as for Table 4 but Using the Second Galactic Potential</caption><thead><tr><th>Ω<sub>bar</sub><br>(km s<sup>-1</sup> kpc<sup>-1</sup>)</th><th>⟨r<sub>min</sub>⟩<br>(kpc)</th><th>⟨r<sub>max</sub>⟩<br>(kpc)</th><th>⟨|z|<sub>max</sub>⟩<br>(kpc)</th><th>⟨e⟩</th><th>⟨E<sub>J</sub>⟩<br>(10<sup>2</sup> km<sup>2</sup> s<sup>-2</sup>)</th></tr></thead><tbody><tr><td colspan="6" style="text-align: center;">Dinescu et al. (1999)</td></tr><tr><td>40</td><td>0.10<sup>0.20</sup><sub>0.02</sub></td><td>7.80<sup>8.13</sup><sub>7.48</sub></td><td>2.61<sup>2.77</sup><sub>2.44</sub></td><td>0.98<sup>0.99</sup><sub>0.95</sub></td><td>−1949.62<sup>−1939.96</sup><sub>−1959.09</sub></td></tr><tr><td>45</td><td>0.48<sup>0.59</sup><sub>0.38</sub></td><td>6.80<sup>6.91</sup><sub>6.68</sub></td><td>2.43<sup>2.61</sup><sub>2.29</sub></td><td>0.87<sup>0.90</sup><sub>0.84</sub></td><td>−1935.27<sup>−1924.63</sup><sub>−1945.71</sub></td></tr><tr><td>50</td><td>0.55<sup>0.63</sup><sub>0.47</sub></td><td>6.64<sup>6.73</sup><sub>6.55</sub></td><td>2.38<sup>2.48</sup><sub>2.27</sub></td><td>0.85<sup>0.87</sup><sub>0.83</sub></td><td>−1920.95<sup>−1909.32</sup><sub>−1932.38</sub></td></tr><tr><td colspan="6" style="text-align: center;">Our paper (L18 @F1)</td></tr><tr><td>40</td><td>0.90<sup>1.07</sup><sub>0.79</sub></td><td>7.12<sup>7.30</sup><sub>6.87</sub></td><td>2.81<sup>3.21</sup><sub>2.36</sub></td><td>0.77<sup>0.80</sup><sub>0.73</sub></td><td>−1842.60<sup>−1825.52</sup><sub>−1859.50</sub></td></tr><tr><td>45</td><td>0.95<sup>1.07</sup><sub>0.83</sub></td><td>7.04<sup>7.20</sup><sub>6.91</sub></td><td>3.05<sup>5.59</sup><sub>2.55</sub></td><td>0.76<sup>0.79</sup><sub>0.73</sub></td><td>−1822.04<sup>−1803.51</sup><sub>−1840.40</sub></td></tr><tr><td>50</td><td>1.04<sup>1.17</sup><sub>0.88</sub></td><td>6.83<sup>7.08</sup><sub>6.60</sub></td><td>2.84<sup>3.35</sup><sub>2.35</sub></td><td>0.74<sup>0.78</sup><sub>0.70</sub></td><td>−1801.47<sup>−1781.43</sup><sub>−1821.23</sub></td></tr></tbody></table>

Note. The orbital parameters here provided refer to the non-inertial reference frame where the bar is at rest. For this reason, we provide the orbital Jacobi energy $E_J$ instead of the angular momentum. The superscript and subscript values represent the upper (84th percentile) and lower (16th percentile) limits, respectively.

fundamental properties of the orbit (i.e., perigalacticons and apogalacticons).

Finally, by comparing Figures 6 and 7, we found that both models seem to reproduce a similar behavior in the radial and vertical excursions, and no important differences due to the details of the modeling are detectable. This means that probably the most important features of models that influence the orbital behavior of this particular cluster are the density law and the angular velocity of the bar that are basically the same in both models.

In Tables 4 and 5, we present average orbital parameters for $\omega$ Cen for the various Galactic potentials obtained by using a simple Monte Carlo procedure to vary the initial conditions. 500 orbits for the first potential and 50,000 orbits for the second potential are computed under $1\sigma$ variation of the absolute PMs, distance, radial velocity, and solar motion. The errors provided for these parameters were assumed to follow a Gaussian distribution. The parameters listed in the tables are the minimum and maximum distances to the Galactic center, the maximum distance from the Galactic plane, and the orbital eccentricity. In Table 4, we provide the angular momentum, while in Table 5 we give the orbital Jacobi energy.

## 6. Conclusions

As part of the HST large program GO-14118+14662 on $\omega$ Cen, in this third paper of the series we presented a revised estimate of the absolute PM of this cluster. Thanks to the exquisite HST-based astrometry and high-precision PMs computed in Paper II, we are able to use background galaxies to link our relative astrometry to an absolute frame. Our measured absolute PM differs significantly from the values available in the literature in magnitude, direction, and precision. A comparison with recent catalogs yields instead an excellent agreement.

We also calculated the orbit of $\omega$ Cen with two different models of the potential for the Milky Way Galaxy. The primary difference between the orbits with our new PM and those based on prior estimates is a larger (by about half a kiloparsec) perigalactic distance obtained with our new value. This may

imply a larger survival expectancy for the cluster in the inner Galactic environment.

We emphasize that the absolute PM presented here is both an independent (with the exception of the registration of the master-frame absolute scale and orientation with the Gaia DR1 catalog) and direct measurement based on the use of a reference sample of background galaxies. As such, the results presented here provide an external check for future astrometric all-sky catalogs, in particular, the upcoming Gaia DR2.

We thank the anonymous referee for useful comments and suggestions. M.L., A.B., D.A., A.J.B., and J.M.R. acknowledge support from STScI grants GO-14118 and 14662. J.G.F.-T. gratefully acknowledges the Chilean BASAL Centro de Excelencia en Astrofísica y Tecnologías Afines (CATA) grant PFB-06/2007. Funding for the GravPot16 software has been provided by the Centre national d'études spatiale (CNES) through grant 0101973 and UTINAM Institute of the Université de Franche-Comte, supported by the Region de Franche-Comte and Institut des Sciences de l'Univers (INSU). A.P.M. acknowledges support by the European Research Council through the ERC-StG 2016 project 716082 "GAL-FOR". Based on observations with the NASA/ESA Hubble Space Telescope, obtained at the Space Telescope Science Institute, which is operated by AURA, Inc., under NASA contract NAS 5-26555. This work has made use of data from the European Space Agency (ESA) mission Gaia (https://www.cosmos.esa.int/gaia), processed by the Gaia Data Processing and Analysis Consortium (DPAC, https://www.cosmos.esa.int/web/gaia/dpac/consortium). Funding for the DPAC has been provided by national institutions, in particular the institutions participating in the Gaia Multilateral Agreement.

## Appendix
Galaxy Finding Charts

In Figure 8, we show the trichromatic stacked images of the 45 galaxies used to compute the absolute PM of $\omega$ Cen. Each source is centered in a $51 \times 51$ WFC3/UVIS pixel$^2$ (~2 × 2 arcsec$^2$) stamp. We used F814W, F606W, and F336W filters for red, green, and blue channels, respectively.
<figure><img src="image_10.png" /><figcaption>**Figure 8.** Trichromatic stacked images of the galaxies adopted in our analysis.</figcaption></figure>

## ORCID iDs

Mattia Libralato https://orcid.org/0000-0001-9673-7397

Andrea Bellini https://orcid.org/0000-0003-3858-637X

Edmundo Moreno D. https://orcid.org/0000-0001-6450-2315

Barbara Pichardo https://orcid.org/0000-0002-8195-8493

Roeland P. van der Marel https://orcid.org/0000-0001-7827-7825

Jay Anderson https://orcid.org/0000-0003-2861-3995

Dániel Apai https://orcid.org/0000-0003-3714-5855

Adam J. Burgasser https://orcid.org/0000-0002-6523-9536

Laura L. Watkins https://orcid.org/0000-0002-1343-134X

## References

* Allen, C., & Santillan, A. 1991, RMXAA, 22, 255

* Altmann, M., Roeser, S., Demleitner, M., Bastian, U., & Schilbach, E. 2017, A&A, 600, L4

* Anderson, J. 1997, PhD thesis, Univ. California

* Anderson, J., & Bedin, L. R. 2010, PASP, 122, 1035

* Anderson, J., & King, I. R. 2006, Instrument Science Report ACS 2006-01, PSFs, Photometry, and Astrometry for the ACS/WFC, http://www.stsci.edu/hst/acs/documents/isrs

* Anderson, J., Sarajedini, A., Bedin, L. R., et al. 2008, AJ, 135, 1035

* Anderson, J., & van der Marel, R. P. 2010, ApJ, 710, 1032

* Bedin, L. R., Piotto, G., Anderson, J., et al. 2004, ApJL, 605, L125

* Bedin, L. R., Piotto, G., Carraro, G., King, I. R., & Anderson, J. 2006, A&A, 460, L27

* Bedin, L. R., Piotto, G., King, I. R., & Anderson, J. 2003, AJ, 126, 247

* Bekki, K., & Freeman, K. C. 2003, MNRAS, 346, L11

* Bellini, A., Anderson, J., & Bedin, L. R. 2011, PASP, 123, 622

* Bellini, A., Anderson, J., Bedin, L. R., et al. 2017a, ApJ, 842, 6

* Bellini, A., Anderson, J., Salaris, M., et al. 2013, ApJL, 769, L32

* Bellini, A., Anderson, J., van der Marel, R. P., et al. 2014, ApJ, 797, 115

* Bellini, A., Anderson, J., van der Marel, R. P., et al. 2017b, ApJ, 842, 6

* Bellini, A., & Bedin, L. R. 2009, PASP, 121, 1419

* Bellini, A., Bedin, L. R., Pichardo, B., et al. 2010, A&A, 513, A51

* Bellini, A., Bianchini, P., Varri, A. L., et al. 2017c, ApJ, 844, 167

* Bellini, A., Libralato, M., Bedin, L. R., et al. 2018, ApJ, 853, 86

* Bellini, A., Milone, A. P., Anderson, J., et al. 2017d, ApJ, 844, 164

* Bissantz, N., Englmaier, P., & Gerhard, O. 2003, MNRAS, 340, 949

* Brunthaler, A., Reid, M. J., Menten, K. M., et al. 2011, AN, 332, 461

* Dinescu, D. I., van Altena, W. F., Girard, T. M., & López, C. E. 1999, AJ, 117, 277

* Drimmel, R. 2000, A&A, 358, L13

* Fernández-Trincado, J. G., Geisler, D., Moreno, E., et al. 2017a, arXiv:1710.07433

* Fernández-Trincado, J. G., Robin, A. C., Moreno, E., Pérez-Villegas, A., & Pichardo, B. 2017b, arXiv:1708.05742

* Fernández-Trincado, J. G., Vivas, A. K., Mateu, C. E., et al. 2015, A&A, 574, A15

* Freeman, K. C. 2001, in ASP Conf. Ser. 228, Dynamics of Star Clusters and the Milky Way, ed. S. Deiters et al. (San Francisco, CA: ASP), 43

* Gaia Collaboration, Prusti, T., de Bruijne, J. H. J., et al. 2016a, A&A, 595, A1

* Gaia Collaboration, Brown, A. G. A., Vallenari, A., et al. 2016b, A&A, 595, A2

* Geffert, M., Hilker, M., Geyer, E. H., & Krämer, G.-H. 2002, in ASP Conf. Ser. 265, Omega Centauri, A Unique Window into Astrophysics, ed. F. van Leeuwen, J. D. Hughes, & G. Piotto (San Francisco, CA: ASP), 399

* Gerhard, O. 2002, in ASP Conf. Ser. 273, The Dynamics, Structure & History of Galaxies: A Workshop in Honour of Professor Ken Freeman, ed. G. S. Da Costa & H. Jerjen (San Francisco, CA: ASP), 73

* Gerhard, O. 2011, MSAIS, 18, 185

* Harris, W. E. 1996, AJ, 112, 1487

* Heyl, J., Caiazzo, I., Richer, H., et al. 2017, ApJ, 850, 186

* Høg, E., Fabricius, C., Makarov, V. V., et al. 2000, A&A, 355, L27

* Kharchenko, N. V., Piskunov, A. E., Schilbach, E., Röser, S., & Scholz, R.-D. 2013, A&A, 558, A53

* Martinez-Medina, L. A., Pichardo, B., Peimbert, A., & Carigi, L. 2017, MNRAS, 468, 3615

* Massari, D., Bellini, A., Ferraro, F. R., et al. 2013, ApJ, 779, 81

* Mayor, M., Meylan, G., Udry, S., et al. 1997, AJ, 114, 1087

* Milone, A. P., Marino, A. F., Bedin, L. R., et al. 2017b, MNRAS, 469, 800

* Milone, A. P., Piotto, G., Renzini, A., et al. 2017a, MNRAS, 464, 3636

* Milone, A. P., Villanova, S., Bedin, L. R., et al. 2006, A&A, 456, 517

* Murray, C. A., Candy, M. P., & Jones, D. H. P. 1965, RGOB, 100, 81

* Navarrete, C., Catelan, M., Contreras Ramos, R., et al. 2017, A&A, 606, C1

* Pichardo, B., Martos, M., & Moreno, E. 2004, ApJ, 609, 144

* Pichardo, B., Martos, M., Moreno, E., & Espresate, J. 2003, ApJ, 582, 230

* Pichardo, B., Moreno, E., Allen, C., et al. 2012, AJ, 143, 73

* Platais, I., Wyse, R. F. G., Hebb, L., Lee, Y.-W., & Rey, S.-C. 2003, ApJL, 591, L127

* Portail, M., Wegg, C., Gerhard, O., & Martinez-Valpuesta, I. 2015, MNRAS, 448, 713

* Reijns, R. A., Seitzer, P., Arnold, R., et al. 2006, A&A, 445, 503

* Robin, A. C., Marshall, D. J., Schultheis, M., & Reylé, C. 2012, A&A, 538, A106

* Robin, A. C., Reylé, C., & Derrière, S. 2003, A&A, 409, 523

* Robin, A. C., Reylé, C., Fliri, J., et al. 2014, A&A, 569, A13

* Roeser, S., Demleitner, M., & Schilbach, E. 2010, AJ, 139, 2440

* Schönrich, R., Binney, J., & Dehnen, W. 2010, MNRAS, 403, 1829

* Sohn, S. T., van der Marel, R. P., Carlin, J. L., et al. 2015, ApJ, 803, 56

* Sormani, M. C., Binney, J., & Magorrian, J. 2015, MNRAS, 454, 1818

* van de Ven, G., van den Bosch, R. C. E., Verolme, E. K., & de Zeeuw, P. T. 2006, A&A, 445, 513
van der Marel, R. P., Alves, D. R., Hardy, E., & Suntzeff, N. B. 2002, AJ, 124, 2639

van Leeuwen, F., & Le Poole, R. S. 2002, in ASP Conf. Ser. 265, Omega Centauri, A Unique Window into Astrophysics, ed. F. van Leeuwen, J. D. Hughes, & G. Piotto (San Francisco, CA: ASP), 41

van Leeuwen, F., Le Poole, R. S., Reijns, R. A., Freeman, K. C., & de Zeeuw, P. T. 2000, A&A, 360, 472

Woolley, R. V. D. R. 1966, ROAn, 2

Zacharias, N., Finch, C., & Frouard, J. 2017, AJ, 153, 166

Zacharias, N., Urban, S. E., Zacharias, M. I., et al. 2000, AJ, 120, 2131