Proper Motion Anomalies Technique

• Proper motion anomalies technique is an astrometric method that examines discrepancies between instantaneous and long-term proper motions to infer the presence of unseen companions.
• It leverages precise data from missions such as Hipparcos and Gaia to probe intermediate orbital separations (1–100 au) that are often inaccessible to transits or radial velocity methods.
• The method is instrumental in exoplanet, binary, and galactic dynamics studies, offering actionable insights into system architecture and unresolved nuclear events.

Proper motion anomalies (PMa) is an astrometric technique that detects or constrains the presence of unseen companions—stellar, substellar, or planetary—through the analysis of discrepancies between a star’s “instantaneous” (short-term) and “long-term” proper motions. The method exploits the high-precision astrometry available from space missions such as Hipparcos, Gaia, and their time baselines, leveraging the reflex acceleration that orbiting secondary bodies impress on the stellar photocenter. Modern implementations also extend to detecting apparent proper motion shifts in unresolved extragalactic sources due to phenomena such as variable or transient events. PMa provides a uniquely powerful way to probe the intermediate-separation regime (1–100 au) largely inaccessible to transits, radial velocities, or direct imaging, and is instrumental in exoplanet, binary, and multiplicity surveys, as well as the study of galactic structure and transient nuclear phenomena.

1. Mathematical and Observational Foundations

The core of the PMa technique is the comparison of astrometric measurements taken at multiple epochs separated by a well-characterized time interval. Given two epochs (e.g., Hipparcos at $t_H=1991.25$ and Gaia EDR3 at $t_G=2016.0$), the long-term proper motion $\mu_{\mathrm{long}}$ is computed by differencing the positions $\theta_G$ and $\theta_H$ and dividing by baseline $\Delta t = t_G - t_H$: $$\mu_{\mathrm{long}} = \frac{\theta_G - \theta_H}{\Delta t}$$ The “instantaneous” proper motion $\mu_{\mathrm{inst}}$ at each epoch comes from the corresponding catalog values (e.g., $\mu_{G,\alpha*}, \mu_{G,\delta}$ for Gaia; $\mu_{H,\alpha*}, \mu_{H,\delta}$ for Hipparcos). The proper motion anomaly vector is then: $t_G=2016.0$0 For stars, this anomaly is interpreted as the signature of an unmodeled acceleration imparted by a companion. The PMa vector can be converted to a tangential velocity anomaly at distance $t_G=2016.0$1 via: $t_G=2016.0$2 Error propagation is handled by quadrature addition of uncertainties from both epochs.

In the VIM (Variability-Induced Movers) context for AGNs and quasars, the technique focuses on the apparent displacement of the photocenter caused by transient or strongly variable events within the Gaia resolution element. The observed proper motion is linked directly to astrometric centroid shifts parameterized by the flare’s brightness and position with respect to the unresolved core (Khamitov et al., 2023).

2. Sensitivity Regimes, Detection Criteria, and Statistical Properties

The PMa method is fundamentally sensitive to companions whose orbital periods are comparable to, but less than, the astrometric time baseline. For Hipparcos–Gaia comparisons ($t_G=2016.0$3 yr), this corresponds to separations of a few to several tens of au. Detection sensitivity (in mass) rises linearly with distance and declines for very wide orbits (smaller accelerations) and very close-in companions (perturbations averaged out or diluted).

The detection significance $t_G=2016.0$4 is computed as: $t_G=2016.0$5 A common threshold is $t_G=2016.0$6, corresponding to a false alarm rate $t_G=2016.0$7 (Kervella et al., 2021). For the nearest stars, Gaia DR2/EDR3 achieves $t_G=2016.0$8 cm s$t_G=2016.0$9 pc$\mu_{\mathrm{long}}$0, enabling planet detection at the sub-Jovian level for the most favorable cases.

Correct error modeling—including catalog covariances and frame corrections—is critical. Systematics such as perspective acceleration, reference-frame tie, and time-window smearing for short-period orbits must be applied. The method is mostly orthogonal to radial velocity and direct imaging in sensitivity, providing unique access to intermediate-period companions (Kervella et al., 2018).

3. Derivation of Physical Parameters and Exclusion Curves

Given a measured or constrained PMa, the stellar tangential acceleration is expressed as: $\mu_{\mathrm{long}}$1 Assuming a single companion of mass $\mu_{\mathrm{long}}$2 at orbital separation $\mu_{\mathrm{long}}$3,

$\mu_{\mathrm{long}}$4

This links the observable $\mu_{\mathrm{long}}$5 to the companion’s mass and orbital radius: $\mu_{\mathrm{long}}$6 Consequently, for non-detections, exclusion curves in the $\mu_{\mathrm{long}}$7 plane can be constructed: $\mu_{\mathrm{long}}$8 Assumptions typically include a circular, coplanar orbit (relative to any resolved disk) and the dominant influence of a single companion. Monte Carlo sampling across phase and inclination yields detection probability curves (Bendahan-West et al., 5 Jan 2026).

4. Applications: Exoplanet Detection, Multiplicity, and Galactic Dynamics

Stellar Companions and Exoplanet Detection

The Hipparcos–Gaia PMa technique has yielded high-precision constraints on the companion population for $\mu_{\mathrm{long}}$9 of Hipparcos stars, with over 30% showing significant ($\theta_G$0) anomalies (Kervella et al., 2018). Notable results include:

• For Proxima Centauri: Planets $\theta_G$1 M$\theta_G$2 beyond 1 au and $\theta_G$3 M$\theta_G$4 beyond 10 au are excluded (Kervella et al., 2018).
• For $\theta_G$5 Eridani: The technique recovers the known $\theta_G$6 Eri b and excludes additional $\theta_G$7 M$\theta_G$8 objects inside $\theta_G$9 au.
• For debris disk gaps (e.g., HD 92945, HD 107146): PMa exclusion curves, combined with JWST coronagraphy and radial velocity data, tightly localize possible planets responsible for observed structures to inner $\theta_H$020–30 au, with mass ranges from $\theta_H$1–$\theta_H$2 M$\theta_H$3 depending on system and separation (Bendahan-West et al., 5 Jan 2026).

Architecture of Planet-Hosting Binaries

PMa is used to build volume-limited catalogs of planet-host stars with significant astrometric accelerations (e.g., 66 out of several thousand Kepler/TESS planet candidates with $\theta_H$4), allowing the identification of unresolved stellar, substellar, or planetary companions at separations of a few to tens of au. This application reveals increased false positive rates for eclipsing binaries among high-PMa TESS objects and informs selection for follow-up (Zhang et al., 2023).

Galactic and Solar Kinematics

A related “proper motion anomaly” formalism enables model-independent measurement of the Sun’s velocity vector relative to the Galactic rest frame by analyzing the mean perpendicular proper motions of kinematically cold streams. This approach is independent of the Galactic potential and achieves $\theta_H$53 km/s accuracy for each velocity component with current data, provided systematic distance uncertainties are controlled (Malhan et al., 2017).

Nuclear Transients and AGN Astrometry (VIM-Induced PMa)

In unresolved AGN and quasars, significant apparent proper motions can arise from luminous transient events (e.g., supernovae, TDEs, luminous blue variable outbursts) offset from the galactic nucleus. The centroid displacement is modeled as

$\theta_H$6

and the resulting $\theta_H$7 is measured over the Gaia time baseline. Candidates are identified with significance criterion $\theta_H$8 (Khamitov et al., 2023). Flares with $\theta_H$9 at $\Delta t = t_G - t_H$0–$\Delta t = t_G - t_H$1 mas are required to match observed PMa. This method allows detection and classification of nuclear SNe and TDEs on spatial scales of $\Delta t = t_G - t_H$2100–200 pc out to $\Delta t = t_G - t_H$3.

5. Methodological Variants and Data Products

Multiple implementations of PMa exist, reflecting both catalog availability and specific scientific goals:

Reference	Baseline (yr)	Typical Target Sample	Key PMa Signal	Companion Sensitivity
Hipparcos–Gaia	$\Delta t = t_G - t_H$425	Bright, nearby stars ($\Delta t = t_G - t_H$5300 pc)	$\Delta t = t_G - t_H$6	$\Delta t = t_G - t_H$7 M$\Delta t = t_G - t_H$8 (nearest M dwarfs), $\Delta t = t_G - t_H$9few M$$\mu_{\mathrm{long}} = \frac{\theta_G - \theta_H}{\Delta t}$$0 at $$\mu_{\mathrm{long}} = \frac{\theta_G - \theta_H}{\Delta t}$$120–100 au (Kervella et al., 2021)
Hipparcos–URAT1	$$\mu_{\mathrm{long}} = \frac{\theta_G - \theta_H}{\Delta t}$$221–23	Northern hemisphere stars	UrHip–Hipparcos difference	Improved sensitivity to binaries (Frouard et al., 2015)
Streams/Solar	variable	Tidal streams & groups	Perpendicular PM anomaly	Solar reflex velocity (model-independent) (Malhan et al., 2017)
Extragalactic	$$\mu_{\mathrm{long}} = \frac{\theta_G - \theta_H}{\Delta t}$$32.8	AGN, quasars	Photocenter shift under flare	Detection of SNe/TDEs, not orbital companions (Khamitov et al., 2023)

Catalog products include tabulated component-wise $$\mu_{\mathrm{long}} = \frac{\theta_G - \theta_H}{\Delta t}$$4, significance $$\mu_{\mathrm{long}} = \frac{\theta_G - \theta_H}{\Delta t}$$5, tangential velocities, and, where appropriate, companion mass constraints or exclusion curves. RUWE (Renormalized Unit Weight Error) augment PMa detections by flagging non-single behaviour (Kervella et al., 2021).

6. Limitations, Systematic Uncertainties, and Extensions

The PMa technique is subject to several limitations:

• Degeneracy: PMa alone constrains $$\mu_{\mathrm{long}} = \frac{\theta_G - \theta_H}{\Delta t}$$6, and cannot distinguish mass from separation without additional data.
• Sensitivity Range: Primarily sensitive to companions with semimajor axes in the $$\mu_{\mathrm{long}} = \frac{\theta_G - \theta_H}{\Delta t}$$71–30 au regime; much wider orbits induce accelerations too small to detect over the baseline, while much shorter orbits are diluted or averaged out (Bendahan-West et al., 5 Jan 2026).
• Unresolved multiplicity: For multiple companions, secular perturbations complicate the interpretation; the formalism assumes a single dominant perturber (Kervella et al., 2018).
• Systematics: Distance uncertainties dominate for galactic applications; for streams, $$\mu_{\mathrm{long}} = \frac{\theta_G - \theta_H}{\Delta t}$$8 distance systematics produce km/s-scale biases in inferred Solar velocities (Malhan et al., 2017). Catalog systematics (frame rotations, saturation for $$\mu_{\mathrm{long}} = \frac{\theta_G - \theta_H}{\Delta t}$$9) necessitate careful calibration.
• Time-window smearing: Gaia PMs average over mission-length windows ($\mu_{\mathrm{inst}}$01–2 yr), suppressing detectability for $\mu_{\mathrm{inst}}$1. Correction factors ($\mu_{\mathrm{inst}}$2) are applied for window smearing and mean subtraction (Kervella et al., 2018).
• Perspective acceleration: Corrections are required for high-proper-motion, nearby stars over the $\mu_{\mathrm{inst}}$325 yr baseline (Kervella et al., 2021).

Ongoing improvements in Gaia data releases—extending time baselines and including epoch astrometry—are expected to further enhance sensitivity, enable full-orbit solutions, and extend the method toward lower mass companions and shorter orbital periods. Synergy with radial velocity, direct imaging, and high-resolution astrometry will further constrain companion properties and system architectures.

7. Practical Implementation and Scientific Impact

Implementation of PMa techniques follows a modular approach—derivation of $\mu_{\mathrm{inst}}$4 with rigorous error propagation, calculation of detection metrics, conversion to physical constraints, and integration with auxiliary indicators (e.g., RUWE, common proper motion pair analysis). The method has enabled:

• A census of the multiplicity fraction in local stellar populations, with the combination of PMa, RUWE, and resolved common proper motion yielding $\mu_{\mathrm{inst}}$5 binarity indicators in the Hipparcos sample (Kervella et al., 2021).
• Exclusion and characterization of long-period giant planets in debris disk and exoKuiper belt systems when direct imaging fails to detect candidates, dramatically reducing the parameter space for hidden planets (Bendahan-West et al., 5 Jan 2026).
• False-positive analysis in the context of transit planet searches, permitting statistical correction for eclipsing binaries and informing observational strategies (Zhang et al., 2023).
• Purely astrometric discovery of luminous transients in extragalactic sources, enabling independent nuclear supernova and TDE rate studies across broad redshift ranges (Khamitov et al., 2023).

A plausible implication is that, as time baselines lengthen and future data releases improve astrometric precision and epoch coverage, PMa will become increasingly valuable for full dynamical characterization of planetary, stellar, and nuclear transient systems across the Milky Way and beyond.