Vast Polar Structure (VPOS) Dynamics

• VPOS is a flattened, nearly polar arrangement of Milky Way satellites characterized by strong spatial flattening and partial orbital-pole coherence.
• Kinematic diagnostics reveal significant orbital-pole clustering, supporting a dynamically meaningful and possibly rotationally stabilized structure.
• Recent studies test diverse origin scenarios—from group infall and tidal debris to modified dark matter halos—to explain the VPOS formation and stability.

The Vast Polar Structure (VPOS) is the thin, approximately polar arrangement of Milky Way halo substructures whose most robustly established component is the satellite-galaxy system, nearly perpendicular to the Galactic disc. In the broader historical usage of the term, the VPOS extended the older “Disc of Satellites” to include young-halo globular clusters and stellar or gaseous streams; in narrower Gaia-era usage, the strongest empirical claim concerns the spatial flattening and partial orbital-pole coherence of the satellites themselves (Pawlowski et al., 2013, Pawlowski et al., 2012, Riley et al., 2020). The VPOS is central to current work on satellite anisotropy, correlated accretion, Local Group structure, and the small-scale interpretation of galaxy formation.

1. Definition and geometric characterization

The VPOS is usually defined geometrically as a flattened, nearly polar configuration of Milky Way satellites. A conservative characterization based on the 11 classical satellites gives a best-fitting plane with $r_{\mathrm{per}}=19.6\,\mathrm{kpc}$, $r_{\mathrm{par}}=129.5\,\mathrm{kpc}$, $c/a=0.182$, $b/a=0.508$, orbital-pole concentration $\Delta_{\mathrm{std}}=29.3^\circ$, and an offset angle $\theta_{\mathrm{VPOS}}=18.9^\circ$ between the mean pole direction and the plane normal; including Canes Venatici I yields closely similar values (Pawlowski et al., 2014). In a principal-component fit to the 11 classical dwarfs, the present-day plane has $D_{\rm rms}=21.3$ kpc, $\delta=0.23$, and an inclination of about $77.3^\circ$ to the Milky Way disc, again emphasizing its strongly polar orientation (Lipnicky et al., 2016).

Later work formalized two common satellite-only realizations. For all 27 then-known Milky Way satellites, “VPOSall” has normal $(l,b)=(155.6^\circ,-3.3^\circ)$, $r_{\mathrm{par}}=129.5\,\mathrm{kpc}$0 kpc, $r_{\mathrm{par}}=129.5\,\mathrm{kpc}$1 kpc, $r_{\mathrm{par}}=129.5\,\mathrm{kpc}$2, and $r_{\mathrm{par}}=129.5\,\mathrm{kpc}$3. Excluding the three main outliers Leo I, Hercules, and Ursa Major I yields the thinner “VPOS-3,” with normal $r_{\mathrm{par}}=129.5\,\mathrm{kpc}$4, $r_{\mathrm{par}}=129.5\,\mathrm{kpc}$5 kpc, $r_{\mathrm{par}}=129.5\,\mathrm{kpc}$6 kpc, $r_{\mathrm{par}}=129.5\,\mathrm{kpc}$7, and $r_{\mathrm{par}}=129.5\,\mathrm{kpc}$8 (Pawlowski et al., 2013). These dual definitions are methodologically important because much of the later dynamical literature hinges on whether the VPOS is treated as a plane fitted to all classical satellites, to a trimmed sample, or to a kinematically selected subset.

The term “polar” refers to the orientation of the plane rather than the location of the normal. A plane normal close to Galactic latitude $r_{\mathrm{par}}=129.5\,\mathrm{kpc}$9 corresponds to a plane nearly perpendicular to the Galactic disc. This convention is used throughout the VPOS literature, including comparisons to the Great Plane of Andromeda and to Local Group-wide planar structures (Pawlowski et al., 2013, Pawlowski et al., 2014).

2. Kinematic diagnostics and orbital-pole coherence

The VPOS is not only a spatial fit. Its dynamical interpretation relies on orbital-pole clustering. In standard notation, the orbital angular momentum is $c/a=0.182$0, the orbital pole is the direction of $c/a=0.182$1, and for stream analyses the orbital-plane normal is commonly estimated from two Galactocentric anchor vectors as $c/a=0.182$2. The angle to the adopted VPOS normal is then used as the alignment diagnostic, with $c/a=0.182$3 corresponding to a cap covering 10% of the sphere; for orbital poles, the co- and counter-rotating caps together cover 20% of the sky (Riley et al., 2020).

Using proper motions for all 11 classical satellites, one influential analysis argued that the coherent orbital alignment of 7 to 9 of 11 satellites demonstrates that the VPOS is a rotationally stabilized structure rather than a merely pressure-supported flattened ellipsoid (Pawlowski et al., 2013). In that work, the 8 most concentrated orbital poles had $c/a=0.182$4 with average direction $c/a=0.182$5, while the 6 most concentrated poles gave $c/a=0.182$6 and average direction $c/a=0.182$7 (Pawlowski et al., 2013). A Gaia-focused synthesis described 8 of the 11 classical satellites as consistent with orbiting in the VPOS and identified Sculptor as the one clear counter-orbiting case, again treating orbital-pole clustering as the decisive indicator of a dynamically meaningful plane (Pawlowski, 2017).

This kinematic interpretation also led to an empirical proper-motion prediction program. If a satellite is assumed to orbit in the VPOS, its measured position and line-of-sight velocity constrain its proper motion to a narrow locus in $c/a=0.182$8 space. The method is explicitly geometric: the best possible orbital alignment is fixed by the current position vector, and the acceptable segment of the proper-motion locus is then truncated by requiring bound, non-radial motion (Pawlowski et al., 2013, Pawlowski, 2017). This framework became one of the main proposed Gaia tests of whether the VPOS is rotationally supported across the broader satellite population (Pawlowski, 2017).

3. Survey expansion, new satellites, and updated plane fits

The post-SDSS and post-DES satellite census was a critical test of whether the VPOS was an accident of northern-hemisphere survey geometry. Incorporating more than a dozen newly discovered southern systems, one update found that the revised best-fitting structure preserved the earlier orientation to within $c/a=0.182$9, retained an rms height of $b/a=0.508$0 kpc, and reduced the shortest distance of the plane from the Milky Way center to only $b/a=0.508$1 kpc (Pawlowski et al., 2015). In that analysis, “VPOS+new” had normal $b/a=0.508$2, $b/a=0.508$3 kpc, $b/a=0.508$4 kpc, $b/a=0.508$5, and $b/a=0.508$6, while the trimmed “VPOS+new-4” had $b/a=0.508$7, $b/a=0.508$8 kpc, $b/a=0.508$9 kpc, and $\Delta_{\mathrm{std}}=29.3^\circ$0 (Pawlowski et al., 2015). This was interpreted as the southern discoveries balancing the fit toward the Galactic center rather than destroying it.

The same paper tested DES-footprint bias directly by randomizing isotropic mock systems within the observed DES footprint. For the ten DES objects, the observed mean and median offsets from VPOSall were 28 and 18 kpc, whereas isotropic realizations restricted to the DES area gave expected mean and median offsets of 41 and 31 kpc; 91% of random realizations had a larger mean offset and 96% had a larger median offset than observed (Pawlowski et al., 2015). A complementary SDSS-footprint study reached a stronger statistical conclusion: the SDSS satellites genuinely added to the significance of the VPOS, and the actual survey footprint biased away from, rather than toward, the observed close alignment between the SDSS and classical planes (Pawlowski, 2015). In that analysis, the combined configuration of classical flattening and SDSS flattening-plus-alignment occurred with frequencies $\Delta_{\mathrm{std}}=29.3^\circ$1 or $\Delta_{\mathrm{std}}=29.3^\circ$2, corresponding to about $\Delta_{\mathrm{std}}=29.3^\circ$3 and $\Delta_{\mathrm{std}}=29.3^\circ$4, respectively, relative to isotropy (Pawlowski, 2015).

Specific new objects were also tested individually. Crater was found to lie close to the VPOS, and adding it slightly improved the alignment of the fitted plane with the LMC orbital pole and with the average stream-normal direction; the ATLAS stream was found to align with the VPOS strongly enough to count as an additional member, whereas the Pisces/Triangulum and PAndAS Milky Way streams did not (Pawlowski et al., 2014). These updates reinforced the view that the VPOS was not erased by subsequent discoveries, although they did not settle the separate question of whether non-satellite tracers as a population define the same structure.

4. Globular clusters, streams, and the scope of the term

Early broad formulations of the VPOS explicitly included satellite galaxies, young-halo globular clusters, and streams. In one foundational treatment, the 30 young-halo globular clusters defined a best-fit plane with normal $\Delta_{\mathrm{std}}=29.3^\circ$5, rms height 11.8 kpc, and an angular separation of less than $\Delta_{\mathrm{std}}=29.3^\circ$6 from the satellite plane; 7 of 14 analyzed streams had normals within $\Delta_{\mathrm{std}}=29.3^\circ$7 of the Disc-of-Satellites normal, with a quoted isotropic probability of $\Delta_{\mathrm{std}}=29.3^\circ$8 for obtaining at least 7 such alignments in a 14-stream sample (Pawlowski et al., 2012). In that broad sense, the VPOS extended from Galactocentric distances of about 10 kpc to 250 kpc and was presented as a correlated halo-wide structure rather than a satellite-only phenomenon (Pawlowski et al., 2012).

Gaia-era reassessments substantially narrowed this picture. Using systemic proper motions from Gaia DR2 and a 64-stream census, one later study found that globular cluster orbital poles are not clustered in the VPOS direction and that stellar stream normals are likewise not clustered around the VPOS normal (Riley et al., 2020). Quantitatively, of 150 globular clusters, 19 were classified as likely VPOS members and 113 as non-members, while for the 30 young-halo clusters the most favorable subgroup contained 9 likely members and an isotropic tail probability of $\Delta_{\mathrm{std}}=29.3^\circ$9 for at least 9 of 30 poles in the VPOS caps (Riley et al., 2020). For the 64 streams, 12 were likely members and the isotropic probability of at least 12 of 64 falling within the VPOS tolerance was $\theta_{\mathrm{VPOS}}=18.9^\circ$0, fully compatible with isotropy (Riley et al., 2020).

The resulting controversy is therefore partly terminological and partly empirical. In early work, “VPOS” referred to a multi-component polar structure containing satellites, young-halo clusters, and streams (Pawlowski et al., 2012). In later work with larger samples and better proper motions, the strongest, least controversial claim remained the satellite-galaxy plane itself, while the broader extension to globular clusters and streams was argued not to survive Gaia-era tests (Riley et al., 2020). The modern literature uses both senses, and any precise statement about the VPOS requires an explicit sample definition.

5. Stability, longevity, and dynamical status

Whether the VPOS is a long-lived dynamical structure or a present-day alignment remains unsettled, and much of the disagreement can be traced to differences in sample definition and gravitational modeling. A backward integration study based on the 11 classical satellites, HST proper motions, and a static spherical Milky Way potential concluded that the VPOS disperses well before a dynamical time. In that treatment, the classical-dwarf plane had $\theta_{\mathrm{VPOS}}=18.9^\circ$1 kpc and $\theta_{\mathrm{VPOS}}=18.9^\circ$2 at the present epoch, but its significance declined within about 0.5 Gyr, leading to the conclusion that the VPOS is not a dynamically stable structure (Lipnicky et al., 2016).

A different line of argument portrays the VPOS as a young phenomenon rather than a stable ancient plane. Using Gaia-era proper motions to classify 16 on-plane and 17 off-plane satellites, one study found an excess of bright satellites on the plane, with 8 bright systems on-plane and 2 off-plane, and reported a distinctive orbital-phase pattern: co-orbiting on-plane satellites were almost all approaching pericenter, while the two counter-orbiting on-plane satellites were leaving their last pericenters (Taibi et al., 2023). In that analysis, the on-plane systems also tended to have the lowest orbital energies for a given angular momentum, and the phase coherence was interpreted as evidence for a late-accreted, dynamically young structure, with a group-infall scenario favored over a long-lived equilibrium configuration (Taibi et al., 2023).

Recent work has also pushed in the opposite direction. In time-evolving Milky Way–LMC potentials, with membership defined kinematically from orbital-pole alignment, one study identified 15 present-day VPOS members, including 9 Milky Way and 6 likely LMC satellites, and found that the structure maintained a roughly constant thickness of $\theta_{\mathrm{VPOS}}=18.9^\circ$3 kpc, flattening $\theta_{\mathrm{VPOS}}=18.9^\circ$4, and a highly coherent normal-vector distribution with $\theta_{\mathrm{VPOS}}=18.9^\circ$5 over the last 5 Gyr (Martínez-García et al., 24 Jul 2025). That paper explicitly argued that earlier transient-VPOS results are reproduced when the full 11 classical satellites are used, but disappear when off-plane and uncertain members are removed, implying that apparent instability is highly sensitive to sample construction (Martínez-García et al., 24 Jul 2025). The present state of the question is therefore not a simple consensus but a methodological divide over which objects count as VPOS members and how the Milky Way–LMC system is modeled.

6. Origin scenarios and cosmological significance

The VPOS is one of the main “plane of satellites” problems in Local Group dynamics. In Millennium II plus semi-analytic modeling, fully VPOS-like systems satisfying both positional flattening and orbital-pole coherence criteria were found at only $\theta_{\mathrm{VPOS}}=18.9^\circ$6–$\theta_{\mathrm{VPOS}}=18.9^\circ$7 frequency, and the combined probability of obtaining both the Milky Way VPOS and the Andromeda GPoA in the same Local Group was estimated to be below $\theta_{\mathrm{VPOS}}=18.9^\circ$8 (Pawlowski et al., 2014). At the same time, a separate ELVIS-based analysis found no strong correlation between the existence or strength of a VPOS-like satellite plane and host-halo virial mass, virial radius, concentration, or paired-versus-isolated environment, implying that the mere existence of the VPOS does not robustly constrain basic Milky Way halo parameters (Pawlowski, 2017).

Several specific mechanisms have been tested and found insufficient. The reflex motion and halo distortion induced by the infall of a massive LMC do create a mild excess of orbital poles toward the LMC/VPOS direction, but the effect is far too small: the observed Milky Way dwarfs show an enhancement of about $\theta_{\mathrm{VPOS}}=18.9^\circ$9–$D_{\rm rms}=21.3$0 in the VPOS region, whereas the simulation-based enhancement is about $D_{\rm rms}=21.3$1 overall in the 50–250 kpc range and only $D_{\rm rms}=21.3$2 for high-angular-momentum particles most comparable to the observed dwarfs (Pawlowski et al., 2021). Thus the LMC may slightly bias orbital-pole distributions, but it does not explain the origin of the VPOS (Pawlowski et al., 2021).

The main origin scenarios therefore remain open. One family invokes correlated accretion or stripping associated with the LMC. In the long-term-stability study, the LMC’s orbital path remained within about $D_{\rm rms}=21.3$3–$D_{\rm rms}=21.3$4 of the VPOS plane, and the authors argued that if the LMC is on its second pericentre, the VPOS may largely originate from satellites stripped during the first passage, with later LMC-induced orbital reorientations adding further members (Martínez-García et al., 24 Jul 2025). A related “young VPOS” study instead favored the late accretion of a bound dwarf group, because the observed phase coherence should phase-mix on $D_{\rm rms}=21.3$5 Gyr timescales and is difficult to maintain in an old structure (Taibi et al., 2023).

A second family invokes tidal debris from a major interaction involving M31. Earlier work proposed that one of the tidal tails from a gas-rich major merger in M31 could have reached the Milky Way and produced a fair reproduction of the VPOS while simultaneously explaining the proximity of the Magellanic Clouds, albeit at the cost of predicting dark-matter-free dwarf galaxies (Fouquet et al., 2012). A more recent test compared Milky Way dwarfs to modeled M31 tidal-tail particles and reported that, for low-mass Milky Way models with total mass less than $D_{\rm rms}=21.3$6, the 6D separation can be less than $D_{\rm rms}=21.3$7 for most M31 modellings, describing the result as an intriguing coincidence rather than a confirmed identification (Akib et al., 1 Jan 2025).

A third family modifies the halo rather than the accretion history. A multistate scalar-field dark-matter model proposed that a spherical ground state plus a polar two-lobed first excited state can reproduce host-galaxy rotation curves and generate a large-scale polar overdensity extending to roughly 40% of the virial radius; with fitted Milky Way parameters, that lobe structure was reported to enclose about 53 of 61 known satellites (Bernal et al., 2024). In that framework, the VPOS is interpreted primarily as a manifestation of anisotropic halo structure rather than of an exceptional merger history, although the same paper explicitly did not derive the full orbital-pole coherence or plane thickness self-consistently (Bernal et al., 2024).

The VPOS therefore remains a live problem precisely because no single explanation has achieved broad closure. The observational core is a thin, nearly polar, and at least partly co-rotating satellite structure. What remains unresolved is whether that structure is best understood as a long-lived dynamical plane, a young phase-coherent accretion event, a relic of tidal debris, a consequence of LMC-linked group infall, or a signature of nonstandard halo physics.