2019ApJ...872...24V

The 3D velocities of M31 and M33 are important for understanding the evolution and cosmological context of the Local Group. Their most massive stars are detected by Gaia, and we use Data Release 2 (DR2) to determine the galaxy proper motions (PMs). We select galaxy members based on, e.g., parallax, PM, color-magnitude diagram location, and local stellar density. The PM rotation of both galaxies is confidently detected, consistent with the known line-of-sight rotation curves: V_rot=-206±86kms^–1 (counterclockwise) for M31, and V_rot=80±52kms^–1 (clockwise) for M33. We measure the center-of-mass PM of each galaxy relative to surrounding background quasars in DR2. This yields that (µ_α*_,µδ_) equals (65±18,-57±15)µasyr^–1 for M31 and (31±19,-29±16)µasyr^–1 for M33. In addition to the listed random errors, each component has an additional residual systematic error of 16µasyr^–1. These results are consistent at 0.8σ and 1.0σ with the (2 and 3 times higher accuracy) measurements already available from Hubble Space Telescope (HST) optical imaging and Very Long Baseline Array water maser observations, respectively. This lends confidence that all these measurements are robust. The new results imply that the M31 orbit toward the Milky Way (MW) is somewhat less radial than previously inferred, V_tan,DR2+HST₌₅₇-31_⁺³⁵kms^–1, and strengthen arguments that M33 may be on its first infall into M31. The results highlight the future potential of Gaia for PM studies beyond the MW satellite system.