---
title: Proper Motion Anomaly in Astrometry
url: https://www.emergentmind.com/topics/proper-motion-anomaly
type: topic
---

# Proper Motion Anomaly in Astrometry

Proper motion anomaly denotes a discrepancy between an observed proper motion and a reference proper motion defined by a kinematic model, a long-term astrometric baseline, or an assumed inertial frame. In the cited literature, the term is used in several distinct but related ways: as the residual between instantaneous and long-term stellar proper motions in binary detection, as a bias in a field-averaged center-of-mass proper motion, as a systematic aberrational term induced by Solar System or stellar acceleration, and as an apparent motion produced by photocenter variability rather than bulk translation [1811.08902] [1102.1504] [2407.19182]. The common structure is a comparison between measured angular motion and an expected linear or bulk motion; the physical interpretation depends on which reference motion is adopted and on which processes are omitted from the model.

## 1. Operational definitions across astrometric subfields

In nearby-star multiplicity work, the proper-motion anomaly is usually defined as the difference between an instantaneous proper motion and a long-term center-of-mass estimate. Kervella and collaborators write
\[
\Delta\boldsymbol\mu=\boldsymbol\mu_{\rm inst}-\boldsymbol\mu_{\rm long\text{-}term},
\]
with the long-term motion derived from Hipparcos and Gaia positions and the instantaneous term taken from Hipparcos or Gaia catalog proper motions [1811.08902]. In the Gaia DR2 Cepheid and RR Lyrae analysis, the long-term motion is
\[
\mu_{\rm HG}=(\vartheta_{\rm G2}-\vartheta_{\rm Hip})/(t_{\rm G2}-t_{\rm Hip}),
\]
and the anomaly vector is
\[
\Delta\mu=\mu_{\rm G2}-\mu_{\rm HG},
\]
with signal-to-noise ratio
\[
S=|\Delta\mu|/\sigma_{\Delta\mu}
\]
used to classify detections [1903.03632].

In the UrHip catalog, the anomaly is defined more specifically as the difference between a short-baseline Hipparcos+URAT1 motion and the Hipparcos proper motion alone,
\[
\Delta\mu\equiv \mu_{\rm UrHip}-\mu_{\rm Hipparcos},
\]
evaluated per coordinate and tested against the combined uncertainty \(\sigma_{\Delta\mu}\) [1509.05626]. There the anomaly is interpreted primarily as evidence for non-linear photocenter motion caused by an unseen companion.

In the Large Magellanic Cloud context, the quantity of interest is not a stellar reflex wobble but a bias in the recovered center-of-mass proper motion. After subtracting modeled rotation and perspective terms from the observed field motions, Bekki defines
\[
\Delta\mu\equiv \langle\mu_{\rm obs}\rangle-\mu_{\rm true}
= \frac{1}{N}\sum_{i=1}^{N}\delta\mu_i,
\]
with variance
\[
{\rm Var}(\Delta\mu)=\sigma^2/N,
\]
so that the anomaly is a field-sampling bias driven by residual local random motions [1102.1504].

In Galactic and extragalactic reference-frame work, “proper motion anomaly” can refer to a coherent systematic term. For the newly derived line-of-sight contribution to secular aberration drift, the extra stellar proper motion is
\[
\Delta\vec\mu=-\,\frac{1}{c}(\vec r\cdot\vec V_0)\,\vec\mu,
\]
where \(\vec V_0\) is the Solar System barycenter velocity and \(\vec r\) is the stellar direction vector [2407.19182]. In this usage, the anomaly is not evidence for binarity but for an unmodeled frame effect.

These definitions are mathematically similar—each is a residual—but physically heterogeneous. A common misconception is to treat all proper-motion anomalies as direct signs of companions. The cited literature shows that anomalies may instead reflect internal kinematics of a galaxy, observer acceleration, stellar acceleration in the Galactic potential, or variability-induced photocenter shifts [1102.1504] [2407.19182] [2309.11308].

## 2. Field-averaged biases in the Large Magellanic Cloud

Recent Hubble Space Telescope studies of the Large Magellanic Cloud measured two-component proper motions in typically small sets of HRC/ACS fields and corrected each field for disk rotation and perspective. Bekki emphasized that these corrections do not remove local random stellar motions in a thick, velocity-dispersed disk, so the mean of \(\sim 10\) fields can differ significantly from the true center-of-mass proper motion [1102.1504]. In his formulation,
\[
{\rm PM(field)}_i={\rm PM(CM)}_{\rm true}+{\rm PMres(field)}_i,
\]
and after subtracting the modeled correction \(VC_i\),
\[
{\rm PMest(CM)}_i={\rm PM(field)}_i-VC_i.
\]
The observational center-of-mass proper motion is then the mean over fields.

The key result is that the residual field term need not average to zero when \(N\) is small. In Bekki’s \(N\)-body models, a fiducial LMC with a thick stellar disk, scale height \(\sim 0.45\,{\rm kpc}\), Toomre \(Q\sim 1.5\), and maximum circular velocity \(V_c\approx 120\,{\rm km\,s^{-1}}\) has central one-dimensional dispersions \(\sigma_x\approx \sigma_y\approx 55\,{\rm km\,s^{-1}}\) and \(\sigma_z\approx 43\,{\rm km\,s^{-1}}\), corresponding at \(D\approx 50\,{\rm kpc}\) to \(\sigma_\mu\approx 0.21\,{\rm mas\,yr^{-1}}\) [1102.1504]. The deviation between the observed and true center-of-mass proper motions can be as large as \(\sim 50\,{\rm km\,s^{-1}}\) (\(\sim 0.21\,{\rm mas\,yr^{-1}}\)).

Because \({\rm Var}(\Delta\mu)=\sigma^2/N\), the root-mean-square bias scales as \(1/\sqrt{N}\). Bekki gives explicit values for the fiducial case: for \(N=10\), \(\Delta\mu_{\rm rms}\approx 0.066\,{\rm mas\,yr^{-1}}\) (\(\sim 15\,{\rm km\,s^{-1}}\)); for \(N=100\), \(\approx 0.021\,{\rm mas\,yr^{-1}}\) (\(\sim 5\,{\rm km\,s^{-1}}\)); and for \(N=1000\), \(\approx 0.007\,{\rm mas\,yr^{-1}}\) (\(\sim 1.6\,{\rm km\,s^{-1}}\)) [1102.1504]. The \(N\)-body experiments confirm the same \(1/\sqrt{N}\) scaling and indicate that several hundred fields are required to suppress the anomaly below \(\sim 10\,{\rm km\,s^{-1}}\).

This usage of proper motion anomaly is methodologically important because it separates deterministic bulk corrections from stochastic residuals. A common misconception is that accurate rotation and perspective modeling is sufficient. Bekki’s analysis shows that “perfect” bulk correction still leaves a residual \(\delta\mu_i\) set by the local dispersion, so reliable recovery of the LMC center-of-mass motion requires large \(N\) as well as accurate kinematic modeling [1102.1504].

## 3. Astrometric binarity, reflex motion, and multiplicity censuses

In stellar multiplicity studies, proper-motion anomaly is primarily a tracer of orbital reflex motion of the photocenter. The UrHip catalog combined Hipparcos positions at epoch \(1991.25\) with URAT1 positions at mean epoch \(2013.45\), over an effective baseline of \(\Delta t\simeq 22.6\,{\rm yr}\), to compute improved proper motions for \(67{,}340\) stars [1509.05626]. The anomaly criterion was
\[
\Delta\mu/\sigma_{\Delta\mu}\ge 3.5
\]
in at least one coordinate. After excluding \(10{,}151\) Hipparcos binaries and \(2{,}130\) Tycho-2–Hipparcos \(\Delta\mu\)-binaries, the single-star subset contained \(56{,}542\) objects, of which \(5{,}054\) new candidates were flagged as likely astrometric binaries [1509.05626].

The Hipparcos–Gaia framework generalized this approach. For nearby stars within \(50\,{\rm pc}\), Kervella et al. analyzed \(6{,}741\) stars and also presented a catalog for \(>99\%\) of the Hipparcos catalog (\(>117{,}000\) stars) [1811.08902]. They adopted a \(3\sigma\) detection criterion and reported, within \(50\,{\rm pc}\), \(\Delta\mu>2\sigma\) for \(38.5\%\) of stars, \(\Delta\mu>3\sigma\) for \(30.6\%\), and \(\Delta\mu>5\sigma\) for \(24.6\%\). The median proper-motion-anomaly uncertainty was \(\sigma(\mu)=234\,\mu{\rm as\,yr^{-1}}\), corresponding to \(\sigma(\Delta v_{\rm tan})\approx 1.0\,{\rm m\,s^{-1}\,pc^{-1}}\), where
\[
\Delta v_{\rm tan}=4.74047\,\frac{\Delta\mu\,[{\rm mas/yr}]}{\varpi\,[{\rm mas}]}\;({\rm m/s})
\]
[1811.08902]. This made it possible both to detect companions and to place upper limits on their masses and separations.

Gaia EDR3 improved the same methodology. In the EDR3 PMa catalog, \(\Delta\mu=\mu_{\rm G3}-\mu_{\rm HG}\) used a \(24.75\,{\rm yr}\) Hipparcos-to-Gaia baseline, Monte Carlo propagation of Hipparcos and Gaia covariance matrices, and a significance threshold of \(S/N>3\) [2109.10912]. The median PMa uncertainty was \(\sigma(\Delta\mu_{\rm G3})\approx 56\,\mu{\rm as\,yr^{-1}}\), corresponding to \(\sigma(\Delta v_{\rm tan})\approx 26\,{\rm cm\,s^{-1}\,pc^{-1}}\), a factor \(\sim 2.5\) improvement in PMa \(S/N\) relative to DR2 [2109.10912]. The catalog contained \(116{,}343\) Hipparcos stars with \(\Delta\mu_{\rm G3}\) measurements; \(37{,}515\) stars had significant PMa detections (\(32\%\)), \(12{,}914\) had common proper-motion bound candidate companions, \(25{,}067\) had \(RUWE>1.4\), and \(50{,}720\) stars (\(43\%\) of the Hipparcos sample) exhibited at least one binarity signal when PMa, CPM, and RUWE were combined [2109.10912].

Variable-star applications demonstrate the astrophysical reach of PMa beyond nearby dwarfs. In Gaia DR2, \(254\) classical Cepheids and \(198\) RR Lyrae stars present in both Hipparcos and Gaia were tested [1903.03632]. Using \(S>5\) as strong detections, \(3<S\le 5\) as clear detections, and \(2<S\le 3\) as suspicious candidates, the study found \(57\) significant PMa detections and \(75\) additional candidates among Cepheids, and \(13\) significant detections and \(61\) candidates among RR Lyrae stars [1903.03632]. For \(28\) binary Cepheids with known spectroscopic orbits, the PMa vectors were combined with radial velocities to infer the orientation of the orbital plane and derive companion masses. The inferred true binary fraction was \(>80\%\) for classical Cepheids, while RR Lyrae stars showed a raw detection fraction of \(\ge 13/198\approx 7\%\) [1903.03632].

This family of methods shares a single premise: a single star should move approximately linearly, whereas a companion induces non-linear photocenter motion. The main limitations are also consistent across studies: sensitivity declines for small parallaxes, for orbital periods shorter than the Hipparcos or Gaia observing windows, and when luminous companions invalidate the assumption that the photocenter coincides with the primary [1903.03632] [2109.10912].

## 4. Secular aberration, Galactic acceleration, and reference-frame systematics

Proper-motion anomalies are not confined to orbital perturbations. At microarcsecond-per-year precision, coherent observer-induced or Galactic-potential-induced terms become measurable. The VLBA Extragalactic Proper Motion Catalog showed that quasars and other distant radio sources are not fixed on the sky at this level. Using \(713\) extragalactic proper motions measured over \(\sim 30\,{\rm yr}\) with VLBI, the secular aberration drift was detected at \(6.3\sigma\) significance [1710.02099]. Modeled as a spheroidal dipole, the fit yielded a square-root dipole power of \(4.89\pm 0.77\,\mu{\rm as\,yr^{-1}}\), a Cartesian amplitude of \(1.69\pm 0.27\,\mu{\rm as\,yr^{-1}}\), and an apex at \((\ell,b)=(275.2^\circ\pm 10.0^\circ,-29.4^\circ\pm 8.8^\circ)\), consistent within \(1\sigma\) with the Galactic center direction [1710.02099]. The paper also emphasized that a tight no-net-rotation constraint can partially absorb the dipole.

For stars in the Milky Way, Liu, Xie, and Zhu derived aberrational proper motions generated by both the Solar System barycenter acceleration and stellar acceleration in the Galactic potential [1306.1287]. With the Galactic aberration constant \(A^{(B)}\simeq 5\,\mu{\rm as\,yr^{-1}}\), they found that within \(200\,{\rm pc}\) of the Galactic center the systematic proper motion can exceed \(1000\,\mu{\rm as\,yr^{-1}}\) under a flat rotation curve, but is limited to about \(150\,\mu{\rm as\,yr^{-1}}\) for a more realistic linearly rising core rotation curve [1306.1287]. They also showed that the Kovalevsky formulation is not appropriate when the orbital period is only a fraction of the light time from the star to the Solar System barycenter, and that for such short-period stars the aberration shifts orbital phase on the sky rather than the inferred orbit shape.

A further refinement was derived for stellar proper motions in Gaia DR3. In addition to the classical galactocentric-acceleration term,
\[
\mu_{\rm GA}=\frac{1}{c}\Bigl[\vec a-(\vec r\cdot \vec a)\vec r\Bigr]\sim 5\,\mu{\rm as/yr},
\]
an extra secular aberration drift arises from the change in the line-of-sight direction itself,
\[
\Delta\vec\mu=-\,\frac{1}{c}(\vec r\cdot \vec V_0)\,\vec\mu
\]
[2407.19182]. This term tends to decrease observed proper motions for stars with Galactic longitudes between \(0^\circ\) and \(180^\circ\), and increase them in the remaining region. If ignored, it induces an additional proper motion of \(>1\,{\rm mas\,yr^{-1}}\) for \(84\) stars and \(>0.02\,{\rm mas\,yr^{-1}}\) for \(5{,}944{,}879\) stars in Gaia DR3, comparable to or larger than typical formal uncertainties at \(G<13\) [2407.19182].

In this domain, a proper-motion anomaly is a correlated frame effect rather than a source-specific perturbation. The practical consequence is that both secular aberration contributions must be modeled if the stellar reference frame is to remain consistent with the extragalactic realization of the ICRS [2407.19182].

## 5. Variability-induced motion and extragalactic “proper motion imitations”

A distinct class of apparent proper-motion anomalies arises when the source photocenter moves because the brightness distribution changes within the instrumental resolution element. For unresolved or marginally resolved AGNs and quasars in Gaia, this is the variability-induced motion effect [2309.11308]. Gaia fits a single centroid to each transit; when a transient occurs within the \(\simeq 60\,{\rm mas}\) Gaia resolution element, the instantaneous photocenter is displaced toward the transient, and over the \(2.8\)-year EDR3 baseline the centroid motion can mimic a linear proper motion.

The formal model in the AGN transient analysis uses a fast-rise, exponential-decay flare profile. With transient flux \(F_T(t)\), quiescent AGN flux \(\langle F_{\rm AGN}\rangle\), and angular offset \(X_T\), the centroid shift is
\[
\Delta\theta(t)=X_T\,\frac{F_T(t)}{\langle F_{\rm AGN}\rangle+F_T(t)}
= X_T\,\frac{R(t)}{1+R(t)},
\]
where \(R(t)=F_T(t)/\langle F_{\rm AGN}\rangle\), and the apparent proper motion is
\[
\mu_{\rm app}(t)=\frac{d\Delta\theta(t)}{dt}
= X_T\,\frac{R'(t)}{(1+R(t))^2}
\]
[2309.11308]. A transient at \(X_T\simeq 50\,{\rm mas}\) with peak brightness equal to the AGN and decay time \(\alpha\simeq 100\,{\rm d}\) yields a peak photocenter shift \(\Delta\theta_{\max}\simeq 25\,{\rm mas}\), and a drop of \(\sim 20\,{\rm mas}\) over the first year of decay, corresponding to \(\mu_{\rm app}\sim 20\,{\rm mas/yr}\) [2309.11308].

The SRG/eROSITA–Gaia EDR3 cross-match produced a catalog of \(248\) spectroscopically confirmed extragalactic sources in the eastern Galactic hemisphere, including \(18\) sources with redshifts measured at the RTT-150 telescope [2309.11308]. The proper motions of these objects are formally significant in Gaia, but their extragalactic nature is confirmed. The paper argues that such anomalies can be explained by transient events near AGN nuclei or quasars, including supernovae outbursts, tidal destruction events in AGNs with double nuclei, variability of large-mass supergiants, and OB associations in the field of view of a variable-brightness AGN [2309.11308].

This usage directly counters a common misinterpretation: significant Gaia proper motions of extragalactic sources need not indicate real bulk transverse motion. They may instead be astrometric artifacts of time-variable photocenters.

## 6. Extended and context-dependent usages

The phrase also appears in specialized contexts where the anomaly is defined relative to a population expectation rather than a catalog baseline. In high proper motion X-ray binaries, Maccarone et al. used “proper motion anomaly” to denote a peculiar transverse speed significantly larger than the typical peculiar speeds of comparable binaries and, in some cases, larger than the local Galactic rotation speed [1402.3616]. For IGR J17544–2619, the measured proper motion components \(\mu_{\rm RA}\cos{\rm Dec}=-13.86\pm 1.34\,{\rm mas\,yr^{-1}}\) and \(\mu_{\rm Dec}=+7.31\pm 1.41\,{\rm mas\,yr^{-1}}\) imply \(\mu_{\rm tot}\simeq 15.6\,{\rm mas\,yr^{-1}}\) and \(v_\perp\simeq 266\,{\rm km\,s^{-1}}\) at \(3.6\,{\rm kpc}\), with the paper quoting a peculiar velocity of about \(275\,{\rm km\,s^{-1}}\) [1402.3616]. For 2A 1822–371, \(\mu_{\rm tot}\simeq 14.8\,{\rm mas\,yr^{-1}}\) gives \(v_\perp\simeq 70\,{\rm km\,s^{-1}}\) at \(1\,{\rm kpc}\) and \(\simeq 350\,{\rm km\,s^{-1}}\) at \(5\,{\rm kpc}\) [1402.3616]. Here the anomaly traces supernova-driven mass loss, natal kicks, and residual eccentricity.

In microlensing, the term can denote an unexpectedly large lens-source relative proper motion,
\[
\mu_{\rm rel}=|\mu_L-\mu_S|,
\]
rather than a residual between catalogs [2108.02499]. In OGLE-2019-BLG-1058, the finite-source model yielded \(t_E=2.53\pm 0.10\,{\rm d}\), \(\rho_*=0.0092\pm 0.0004\), \(\theta_E=0.122\pm 0.011\,{\rm mas}\), and \(\mu_{\rm rel}=17.58\pm 1.66\,{\rm mas\,yr^{-1}}\), unusually large compared with typical bulge–bulge events of \(\sim 5\)–\(10\,{\rm mas\,yr^{-1}}\) [2108.02499]. A marginal terrestrial parallax signal initially suggested a free-floating-planet candidate in the disk, but direct measurement of the source proper motion,
\[
\mu_S(\ell,b)=(-12.46\pm 2.41,\,+0.24\pm 2.41)\,{\rm mas\,yr^{-1}},
\]
showed that the large \(\mu_{\rm rel}\) was caused by an extreme source proper motion; the posterior then shifted to a low-mass \(\sim 0.08\)–\(0.09\,M_\odot\) stellar lens in the bulge [2108.02499].

These examples show that proper motion anomaly is not a uniquely standardized observable across astronomy. The phrase consistently identifies a mismatch between observed angular motion and an adopted reference expectation, but the reference may be a catalog baseline, a bulk center-of-mass model, an inertial frame, a population kinematic prior, or a fixed-photocenter assumption. The interpretation is therefore inseparable from the astrometric model against which the motion is judged [1811.08902] [1102.1504] [2309.11308].

Source: https://www.emergentmind.com/topics/proper-motion-anomaly