Center-of-Mass Coupling in Physics

• Center-of-mass coupling is the phenomenon where collective motion and internal degrees of freedom become intertwined under non-ideal conditions such as anharmonicity or trap anisotropy.
• It arises from mechanisms like anharmonic confinement, magnetic field anisotropies, and relativistic corrections, impacting systems from ultracold atoms to nuclear density functionals.
• Engineered center-of-mass coupling is exploited as a resource in precision metrology, molecule formation control, and quantum interface designs in trapped-ion and nanophotonic systems.

Searching arXiv for recent and foundational papers on center-of-mass coupling across relevant subfields. Center-of-mass coupling denotes the failure, modification, or controlled exploitation of the usual separation between collective translational motion and internal or relative degrees of freedom. In many-body and few-body physics, the center of mass (CoM) is ordinarily introduced to factor out trivial global motion, but this decoupling is exact only under restrictive conditions: translational invariance, harmonic confinement, identical trapping tensors, or specific gauge structures. When those conditions are violated—or when relativistic, electromagnetic, anharmonic, or symmetry-restoration effects are retained—the CoM becomes dynamically entangled with internal coordinates, relative motion, spin, fields, or collective excitations. The resulting phenomena span asymptotic gravitational-wave charges, excitons in magnetic and electrostatic environments, ultracold-atom resonances, quantum metrology, polymer transport, nuclear density functionals, and relativistic composite dynamics (Khairnar et al., 25 Mar 2026, Le et al., 2 Mar 2026, Sala et al., 2015, Gietka, 2022, Kafker et al., 12 Mar 2025, Krause et al., 2016).

1. Conceptual scope and exact versus approximate separation

The canonical CoM transformation replaces particle coordinates by a collective coordinate and a set of internal or relative coordinates. In translationally invariant systems this often diagonalizes the kinetic energy and can isolate an internal Hamiltonian. A representative statement appears in internal density-functional theory for isolated self-bound systems, where the Hamiltonian is decomposed as

$$\hat H=\hat H_{\rm cm}+\hat H_{\rm int}, \qquad \hat H_{\rm cm}=-\frac{\hbar^2}{2mN}\Delta_{\mathbf R},$$

with the full wave function separating as $$\Psi(\{\mathbf r_i\})=\Gamma(\mathbf R)\,\psi_{\rm int}(\{\boldsymbol\xi_\alpha\})$$ (Messud, 2012).

That ideal separation is not generic. In ultracold atoms trapped in anharmonic single-well potentials, the total Hamiltonian takes the form

$$H(\mathbf r,\mathbf R)=H_{\rm rel}(\mathbf r)+H_{\rm cm}(\mathbf R)+W(\mathbf r,\mathbf R),$$

where the coupling term $W(\mathbf r,\mathbf R)$ is generated by quartic and sextic trap corrections and drives inelastic confinement-induced resonances (Sala et al., 2015). In merged optical tweezers, different trap frequencies or anisotropies produce an explicit bilinear coupling,

$$V_{\rm cpl}({\bf R},{\boldsymbol{\mathcal R}})= \mu\,\bigl[{\boldsymbol{\mathcal R}}-{\boldsymbol{\mathcal R}}_0\bigr]^T \Delta\boldsymbol{\omega}^2 \bigl[{\bf R}-{\bf R}_0\bigr],$$

which vanishes only when the two traps are identical and aligned (Bird et al., 2024). In spin-orbit-coupled Bose gases and Coulomb crystals, the CoM mode is intentionally coupled to spin variables and used as a metrological resource (Gietka, 2022).

A broad misconception is that CoM separation is a merely kinematic preliminary step that can always be imposed without dynamical consequence. The literature shows otherwise. Exact separation may survive in one formulation but fail in another approximation, or may require nontrivial conserved quantities such as pseudomomentum. Conversely, in several settings the coupling itself is the physically relevant object rather than a nuisance to be removed (Le et al., 2 Mar 2026, Bird et al., 2024).

2. Mechanisms that generate center-of-mass coupling

Several distinct mechanisms recur across subfields.

Anharmonic confinement generates mixed terms between CoM and internal coordinates. For two identical atoms in an optical-lattice well expanded to sixth order, the coupling term contains contributions such as

$$W(\mathbf r,\mathbf R)=\sum_j V_j\Bigl[-\,k_j^4\,r_j^2 R_j^2 +\tfrac13\,k_j^6\,r_j^2R_j^4 +\tfrac1{12}k_j^6\,r_j^4R_j^2\Bigr],$$

which vanish in the purely harmonic approximation but mediate inelastic confinement-induced resonances once anharmonicity is kept (Sala et al., 2015).

Trap inhomogeneity and anisotropy couple relative and CoM motion in merged traps. When two atoms experience different trapping tensors, the potential cannot be written as a sum of pure-relative and pure-CoM pieces, and the bilinear $V_{\rm cpl}$ above shifts crossing positions and splittings relevant for molecule formation (Bird et al., 2024).

Magnetic fields with anisotropic masses obstruct naive CoM separation for 2D magnetoexcitons. In anisotropic materials the conserved quantity is not the canonical momentum but the pseudomomentum

$$\hat{\mathbf K} =\hat{\mathbf P} -\tfrac{e}{2}\,\mathbf B\times(\mathbf r_e-\mathbf r_h),$$

and exact separation requires a pseudomomentum-dependent gauge transformation (Le et al., 2 Mar 2026). Approximate factorized wave functions are only strictly valid in the heavy-hole limit $(\rho_x,\rho_y\to0)$ (Le et al., 2 Mar 2026).

Spatially varying external potentials couple exciton internal and CoM motion in 2D semiconductors. For a weak periodic electrostatic potential, the leading coupling is

$$V_{\rm coup}(\mathbf R,\mathbf r)\simeq \mathbf r\!\cdot\!\nabla U(\mathbf R) = i\sum_{\mathbf G\neq0}U_{\mathbf G}(\mathbf G\cdot \mathbf r)e^{i\mathbf G\cdot\mathbf R},$$

so that CoM momentum is shifted by a reciprocal vector $$\Psi(\{\mathbf r_i\})=\Gamma(\mathbf R)\,\psi_{\rm int}(\{\boldsymbol\xi_\alpha\})$$0 while internal angular momentum changes by $$\Psi(\{\mathbf r_i\})=\Gamma(\mathbf R)\,\psi_{\rm int}(\{\boldsymbol\xi_\alpha\})$$1 (Lu et al., 2024).

Relativistic mass-energy equivalence couples CoM and internal dynamics even in free motion. Starting from

$$\Psi(\{\mathbf r_i\})=\Gamma(\mathbf R)\,\psi_{\rm int}(\{\boldsymbol\xi_\alpha\})$$2

the weakly relativistic expansion contains the term

$$\Psi(\{\mathbf r_i\})=\Gamma(\mathbf R)\,\psi_{\rm int}(\{\boldsymbol\xi_\alpha\})$$3

which is the leading-order CoM–internal coupling and yields time dilation of internal dynamics (Krause et al., 2016).

Field-mediated and spin-dependent couplings arise when a CoM coordinate interacts through internal degrees of freedom or collective modes. In a 1+1D open quantum system, the coupling

$$\Psi(\{\mathbf r_i\})=\Gamma(\mathbf R)\,\psi_{\rm int}(\{\boldsymbol\xi_\alpha\})$$4

produces a nonlinear tripartite interaction between CoM, internal oscillator, and field, leading after elimination of the environment to a non-Markovian quantum Langevin equation for the CoM (2101.07440). In trapped-ion crystals, an oscillating optical dipole force yields

$$\Psi(\{\mathbf r_i\})=\Gamma(\mathbf R)\,\psi_{\rm int}(\{\boldsymbol\xi_\alpha\})$$5

thereby coupling collective spin to CoM displacement (Affolter et al., 2020).

3. Exact formulations and conserved-quantity methods

A central divide in the literature is between approximate factorization and exact or symmetry-consistent treatments.

For anisotropic 2D magnetoexcitons in a perpendicular magnetic field, exact separation is achieved by using eigenfunctions of the conserved pseudomomentum,

$$\Psi(\{\mathbf r_i\})=\Gamma(\mathbf R)\,\psi_{\rm int}(\{\boldsymbol\xi_\alpha\})$$6

which yields a relative Hamiltonian depending only on $$\Psi(\{\mathbf r_i\})=\Gamma(\mathbf R)\,\psi_{\rm int}(\{\boldsymbol\xi_\alpha\})$$7, derivatives, and $$\Psi(\{\mathbf r_i\})=\Gamma(\mathbf R)\,\psi_{\rm int}(\{\boldsymbol\xi_\alpha\})$$8 after subtracting the pure CoM kinetic energy (Le et al., 2 Mar 2026). The resulting exact relative Hamiltonian contains anisotropy-dependent orbital-Zeeman and diamagnetic-renormalization terms absent from the stationary-CoM approximation (Le et al., 2 Mar 2026).

In nuclear structure, translational symmetry is restored by projection rather than by simply subtracting an average kinetic term. The Peierls–Yoccoz projector

$$\Psi(\{\mathbf r_i\})=\Gamma(\mathbf R)\,\psi_{\rm int}(\{\boldsymbol\xi_\alpha\})$$9

produces the projected state

$$H(\mathbf r,\mathbf R)=H_{\rm rel}(\mathbf r)+H_{\rm cm}(\mathbf R)+W(\mathbf r,\mathbf R),$$0

and defines a CoM correction energy

$$H(\mathbf r,\mathbf R)=H_{\rm rel}(\mathbf r)+H_{\rm cm}(\mathbf R)+W(\mathbf r,\mathbf R),$$1

with $$H(\mathbf r,\mathbf R)=H_{\rm rel}(\mathbf r)+H_{\rm cm}(\mathbf R)+W(\mathbf r,\mathbf R),$$2 obtained from overlap and Hamiltonian kernels (Kafker et al., 12 Mar 2025). This treatment incorporates both kinetic and interaction-induced CoM correlations, unlike the common estimate $$H(\mathbf r,\mathbf R)=H_{\rm rel}(\mathbf r)+H_{\rm cm}(\mathbf R)+W(\mathbf r,\mathbf R),$$3 (Kafker et al., 12 Mar 2025).

In asymptotic general relativity, the CoM is encoded in Bondi–van der Burg–Metzner–Sachs charges rather than a finite-particle coordinate. For gravitational radiation, the spatial center-of-mass charge is

$$H(\mathbf r,\mathbf R)=H_{\rm rel}(\mathbf r)+H_{\rm cm}(\mathbf R)+W(\mathbf r,\mathbf R),$$4

and under an infinitesimal boost $$H(\mathbf r,\mathbf R)=H_{\rm rel}(\mathbf r)+H_{\rm cm}(\mathbf R)+W(\mathbf r,\mathbf R),$$5 it transforms as

$$H(\mathbf r,\mathbf R)=H_{\rm rel}(\mathbf r)+H_{\rm cm}(\mathbf R)+W(\mathbf r,\mathbf R),$$6

This formulation turns CoM fixing into a problem of selecting a physically meaningful BMS frame for numerical-relativity waveforms (Khairnar et al., 25 Mar 2026).

These examples indicate that exact CoM treatment often requires enlarging the symmetry framework rather than merely changing coordinates. This suggests a general methodological point: whether CoM coupling is removable depends on the conserved generators and asymptotic structure of the theory, not only on the algebraic form of the Hamiltonian.

4. Observable consequences across physical systems

The physical manifestations of CoM coupling are diverse but structurally related.

In ultracold atoms, anharmonic CoM–relative coupling creates avoided crossings between unbound trap states and CoM-excited molecular states. Resonance positions follow

$$H(\mathbf r,\mathbf R)=H_{\rm rel}(\mathbf r)+H_{\rm cm}(\mathbf R)+W(\mathbf r,\mathbf R),$$7

and the coupling strength is set by matrix elements

$$H(\mathbf r,\mathbf R)=H_{\rm rel}(\mathbf r)+H_{\rm cm}(\mathbf R)+W(\mathbf r,\mathbf R),$$8

These inelastic confinement-induced resonances lead to coherent molecule formation, losses, and heating, and occur for both positive and negative $$H(\mathbf r,\mathbf R)=H_{\rm rel}(\mathbf r)+H_{\rm cm}(\mathbf R)+W(\mathbf r,\mathbf R),$$9-wave scattering lengths (Sala et al., 2015).

In mergoassociation, CoM coupling shifts the avoided crossings used for molecule formation in merged tweezers. For large trap separation, the splitting at the ground-state crossing can be estimated by an effective Rabi frequency $W(\mathbf r,\mathbf R)$0, while crossings involving motional excitations receive matrix elements proportional to $W(\mathbf r,\mathbf R)$1 times overlap factors (Bird et al., 2024). For RbCs with $W(\mathbf r,\mathbf R)$2 kHz, the ground-state crossing peaks at $W(\mathbf r,\mathbf R)$3 kHz when $W(\mathbf r,\mathbf R)$4 and decreases as the trap-frequency mismatch grows (Bird et al., 2024).

In 2D semiconductors, periodic electrostatic potentials hybridize $W(\mathbf r,\mathbf R)$5, $W(\mathbf r,\mathbf R)$6, and $W(\mathbf r,\mathbf R)$7 exciton Rydberg states with different CoM momenta, producing multiple absorption peaks and valley-coherent bright excitons when combined with exchange and in-plane electric fields (Lu et al., 2024). In anisotropic 2D magnetoexcitons, exact separation reveals anisotropy-dependent magnetic couplings that significantly influence magnetoexciton energies, diamagnetic coefficients, and probability densities for the ten lowest states in monolayer black phosphorus and titanium trisulfide (Le et al., 2 Mar 2026).

In polymer solutions, nanoparticles exhibit short-time subdiffusion not explained solely by coupling to polymer segmental relaxations. Simulations indicate an additional coupling between nanoparticle dynamics and polymer CoM motion, retained even without many-body hydrodynamic interactions when long-time dynamics are matched (Chen et al., 2017). The long-time diffusivity still follows $W(\mathbf r,\mathbf R)$8, but the short-time subdiffusive exponent declines smoothly rather than showing the sharp transition predicted by the original coupling theory (Chen et al., 2017).

In classical mechanics, even simple rolling systems can display CoM participation in internal normal modes. For a hoop and a cylinder joined by a spring, the Lagrangian contains the kinetic cross term $W(\mathbf r,\mathbf R)$9, and external static friction can drive the CoM so that it does not remain at rest during normal modes unless $$V_{\rm cpl}({\bf R},{\boldsymbol{\mathcal R}})= \mu\,\bigl[{\boldsymbol{\mathcal R}}-{\boldsymbol{\mathcal R}}_0\bigr]^T \Delta\boldsymbol{\omega}^2 \bigl[{\bf R}-{\bf R}_0\bigr],$$0 or the spring geometry is tuned to make the friction forces vanish (Nyitray, 2 Aug 2025).

5. CoM coupling as a resource: sensing, metrology, and actuation

Not all CoM coupling is a correction to be removed; several papers exploit it operationally.

In trapped-ion crystals, a driven CoM displacement $$V_{\rm cpl}({\bf R},{\boldsymbol{\mathcal R}})= \mu\,\bigl[{\boldsymbol{\mathcal R}}-{\boldsymbol{\mathcal R}}_0\bigr]^T \Delta\boldsymbol{\omega}^2 \bigl[{\bf R}-{\bf R}_0\bigr],$$1 enters the spin-dependent optical-dipole-force Hamiltonian and generates an induced spin precession

$$V_{\rm cpl}({\bf R},{\boldsymbol{\mathcal R}})= \mu\,\bigl[{\boldsymbol{\mathcal R}}-{\boldsymbol{\mathcal R}}_0\bigr]^T \Delta\boldsymbol{\omega}^2 \bigl[{\bf R}-{\bf R}_0\bigr],$$2

on resonance (Affolter et al., 2020). The resulting signal enables phase-coherent detection of CoM motion with a measured $$V_{\rm cpl}({\bf R},{\boldsymbol{\mathcal R}})= \mu\,\bigl[{\boldsymbol{\mathcal R}}-{\boldsymbol{\mathcal R}}_0\bigr]^T \Delta\boldsymbol{\omega}^2 \bigl[{\bf R}-{\bf R}_0\bigr],$$3 pm displacement at single-measurement signal-to-noise ratio $$V_{\rm cpl}({\bf R},{\boldsymbol{\mathcal R}})= \mu\,\bigl[{\boldsymbol{\mathcal R}}-{\boldsymbol{\mathcal R}}_0\bigr]^T \Delta\boldsymbol{\omega}^2 \bigl[{\bf R}-{\bf R}_0\bigr],$$4, corresponding to $$V_{\rm cpl}({\bf R},{\boldsymbol{\mathcal R}})= \mu\,\bigl[{\boldsymbol{\mathcal R}}-{\boldsymbol{\mathcal R}}_0\bigr]^T \Delta\boldsymbol{\omega}^2 \bigl[{\bf R}-{\bf R}_0\bigr],$$5 pm/$$V_{\rm cpl}({\bf R},{\boldsymbol{\mathcal R}})= \mu\,\bigl[{\boldsymbol{\mathcal R}}-{\boldsymbol{\mathcal R}}_0\bigr]^T \Delta\boldsymbol{\omega}^2 \bigl[{\bf R}-{\bf R}_0\bigr],$$6, $$V_{\rm cpl}({\bf R},{\boldsymbol{\mathcal R}})= \mu\,\bigl[{\boldsymbol{\mathcal R}}-{\boldsymbol{\mathcal R}}_0\bigr]^T \Delta\boldsymbol{\omega}^2 \bigl[{\bf R}-{\bf R}_0\bigr],$$7 yN/ion/$$V_{\rm cpl}({\bf R},{\boldsymbol{\mathcal R}})= \mu\,\bigl[{\boldsymbol{\mathcal R}}-{\boldsymbol{\mathcal R}}_0\bigr]^T \Delta\boldsymbol{\omega}^2 \bigl[{\bf R}-{\bf R}_0\bigr],$$8, and $$V_{\rm cpl}({\bf R},{\boldsymbol{\mathcal R}})= \mu\,\bigl[{\boldsymbol{\mathcal R}}-{\boldsymbol{\mathcal R}}_0\bigr]^T \Delta\boldsymbol{\omega}^2 \bigl[{\bf R}-{\bf R}_0\bigr],$$9V/m/$$W(\mathbf r,\mathbf R)=\sum_j V_j\Bigl[-\,k_j^4\,r_j^2 R_j^2 +\tfrac13\,k_j^6\,r_j^2R_j^4 +\tfrac1{12}k_j^6\,r_j^4R_j^2\Bigr],$$0 sensitivities (Affolter et al., 2020). Because this displacement is $$W(\mathbf r,\mathbf R)=\sum_j V_j\Bigl[-\,k_j^4\,r_j^2 R_j^2 +\tfrac13\,k_j^6\,r_j^2R_j^4 +\tfrac1{12}k_j^6\,r_j^4R_j^2\Bigr],$$1 times smaller than the zero-point fluctuations, the work isolates measurement imprecision free from COM-mode noise under the specific off-resonant protocol employed (Affolter et al., 2020).

In quantum metrology with spin-orbit-coupled Bose-Einstein condensates and Coulomb crystals, the CoM oscillator is coupled to spin through Rabi- or Dicke-type interactions. For the spin-orbit system,

$$W(\mathbf r,\mathbf R)=\sum_j V_j\Bigl[-\,k_j^4\,r_j^2 R_j^2 +\tfrac13\,k_j^6\,r_j^2R_j^4 +\tfrac1{12}k_j^6\,r_j^4R_j^2\Bigr],$$2

while for the Coulomb crystal,

$$W(\mathbf r,\mathbf R)=\sum_j V_j\Bigl[-\,k_j^4\,r_j^2 R_j^2 +\tfrac13\,k_j^6\,r_j^2R_j^4 +\tfrac1{12}k_j^6\,r_j^4R_j^2\Bigr],$$3

In both cases the unknown parameter is imprinted onto the squeezing of the CoM mode, and generalized quadrature measurements can saturate Cramér–Rao bounds with Heisenberg-like scaling in the number of CoM excitations (Gietka, 2022).

In nanophotonics, CoM acceleration of coupled waveguides is produced by injecting a frequency-degenerate superposition of modes with opposite transverse parity. The cross term

$$W(\mathbf r,\mathbf R)=\sum_j V_j\Bigl[-\,k_j^4\,r_j^2 R_j^2 +\tfrac13\,k_j^6\,r_j^2R_j^4 +\tfrac1{12}k_j^6\,r_j^4R_j^2\Bigr],$$4

generates a transverse beating force of the same sign on both waveguides, unlike ordinary eigenmode forces which are equal and opposite (Fernandes et al., 2017). The beating-force amplitude is comparable to the eigenmode forces, and correct Maxwell-stress-tensor evaluation requires including the end faces of the infinitesimal slice because the stress tensor cross terms vary along $$W(\mathbf r,\mathbf R)=\sum_j V_j\Bigl[-\,k_j^4\,r_j^2 R_j^2 +\tfrac13\,k_j^6\,r_j^2R_j^4 +\tfrac1{12}k_j^6\,r_j^4R_j^2\Bigr],$$5 (Fernandes et al., 2017).

These cases show that CoM coupling can be engineered as a readout transducer, force amplifier, or design degree of freedom. A plausible implication is that future precision platforms will increasingly treat CoM modes as controllable carriers of information rather than passive spectators.

6. Symmetry restoration, frame fixing, and theoretical controversies

A major theme in the literature is that neglecting CoM coupling can produce quantitatively significant or conceptually inconsistent results.

In nuclear density-functional theory, translational symmetry breaking occurs for all nuclei, and the magnitude of the CoM correction to binding energies varies between $$W(\mathbf r,\mathbf R)=\sum_j V_j\Bigl[-\,k_j^4\,r_j^2 R_j^2 +\tfrac13\,k_j^6\,r_j^2R_j^4 +\tfrac1{12}k_j^6\,r_j^4R_j^2\Bigr],$$6 and $$W(\mathbf r,\mathbf R)=\sum_j V_j\Bigl[-\,k_j^4\,r_j^2 R_j^2 +\tfrac13\,k_j^6\,r_j^2R_j^4 +\tfrac1{12}k_j^6\,r_j^4R_j^2\Bigr],$$7 MeV from light to heavy nuclei (Kafker et al., 12 Mar 2025). This scale exceeds both the $$W(\mathbf r,\mathbf R)=\sum_j V_j\Bigl[-\,k_j^4\,r_j^2 R_j^2 +\tfrac13\,k_j^6\,r_j^2R_j^4 +\tfrac1{12}k_j^6\,r_j^4R_j^2\Bigr],$$8 MeV RMS error of the Bethe–Weizsäcker mass formula and the typical $$W(\mathbf r,\mathbf R)=\sum_j V_j\Bigl[-\,k_j^4\,r_j^2 R_j^2 +\tfrac13\,k_j^6\,r_j^2R_j^4 +\tfrac1{12}k_j^6\,r_j^4R_j^2\Bigr],$$9–$V_{\rm cpl}$0 MeV RMS energy deviations of unprojected DFT mass fits (Kafker et al., 12 Mar 2025). The conclusion is not merely that a correction should be added, but that CoM fluctuations should be incorporated uniformly into the EDF, potentially via variation after projection or an explicit functional guided by the generator-coordinate treatment (Kafker et al., 12 Mar 2025).

In internal DFT for self-bound systems, a local CoM-correlations potential is derived as a parameter-free alternative to projection techniques. The pure CoM-correlation contribution $V_{\rm cpl}$1 is defined within the exchange-correlation sector and leads to a local potential

$V_{\rm cpl}$2

The proposed implementation uses multiconvolutions evaluated through FFTs and is intended for both fermions and bosons (Messud, 2012). This suggests a computational route intermediate between exact projection and crude kinetic subtractions.

In gravitational-wave modeling, CoM fixing is a BMS-frame problem. Past work approximated the numerical-relativity CoM motion by a linear fit,

$V_{\rm cpl}$3

where $V_{\rm cpl}$4 and $V_{\rm cpl}$5 represent boost and translation parameters. The post-Newtonian treatment augments this by the leading oscillatory CoM charge,

$V_{\rm cpl}$6

whose tangential direction $V_{\rm cpl}$7 captures the physical out-spiralling oscillations (Khairnar et al., 25 Mar 2026). The improved fitting ansatz substantially reduces sensitivity to the fitting window, with the largest reported robustness gains being by a factor of $V_{\rm cpl}$8 for the boost vector and $V_{\rm cpl}$9 for the translation vector (Khairnar et al., 25 Mar 2026).

A more far-reaching controversy appears in the COM-separation critique of standard many-body treatments in quantum chemistry and solid-state physics. That work argues that beyond the Born–Oppenheimer approximation, a COM-compatible description requires “hypervibrations” including phonons, rotons, and translons, and that pure electron-phonon models are justified only in adiabatic systems (Svrcek, 2010). It further claims that rotonic and translonic couplings are necessary for symmetry breaking and for the proper interpretation of Jahn–Teller effects and superconductivity (Svrcek, 2010). Since these claims are stronger than the more conventional treatments elsewhere in the literature, they are best regarded as a specific theoretical program rather than a consensus position.

Overall, the recurring controversy is not whether CoM variables exist, but whether approximate elimination of CoM motion is benign. The record across nuclear, condensed-matter, ultracold, and gravitational applications indicates that it often is not.

7. Unifying perspective and open directions

Across otherwise unrelated systems, CoM coupling appears in three recurrent roles.

First, it can be a spurious degree of freedom to be projected out or fixed. This is the case for translational symmetry restoration in nuclei and for BMS-frame fixing of numerical-relativity waveforms (Kafker et al., 12 Mar 2025, Khairnar et al., 25 Mar 2026).

Second, it can be an intrinsic dynamical channel that modifies spectra, transport, or resonance structure. Exact magnetoexciton separation, inelastic confinement-induced resonances, exciton miniband hybridization, relativistic time-dilation coupling, and polymer-mediated nanoparticle subdiffusion all fall into this category (Le et al., 2 Mar 2026, Sala et al., 2015, Lu et al., 2024, Krause et al., 2016, Chen et al., 2017).

Third, it can be a resource for control and measurement, as in trapped-ion displacement sensing, CoM-excitation metrology, and optical beating forces in nanowaveguides (Affolter et al., 2020, Gietka, 2022, Fernandes et al., 2017).

A useful synthesis is that CoM coupling becomes unavoidable whenever the effective separation theorem is weakened: by anharmonicity, anisotropy, gauge fields, exchange, relativistic inertia, open-system backaction, or symmetry breaking. This suggests that future work will continue along two complementary directions. One is to improve exact or symmetry-respecting elimination schemes, especially in EDFs, waveform modeling, and many-body numerical methods. The other is to deliberately engineer CoM coupling for sensing, molecule formation, and coherent control. In both directions, the central lesson is the same: the center of mass is not universally a passive coordinate, but a dynamical sector whose coupling structure can encode essential physics.