---
title: Center-of-Mass Frame Fixing
url: https://www.emergentmind.com/topics/center-of-mass-com-frame-fixing
type: topic
---

# Center-of-Mass Frame Fixing

Searching arXiv for recent and directly relevant papers on center-of-mass frame fixing across domains.
Center-of-mass (CoM) frame fixing denotes the class of procedures by which a physical description is reformulated so that bulk translational motion, frame drift, or off-center coordinate choices do not contaminate intrinsic observables. The operation takes different technical forms in different fields: projection onto zero total momentum in nuclear density functional theory, estimation of an offset vector between a sensor and a spacecraft’s true center of mass, BMS-frame transformations of numerical-relativity waveforms, laboratory-to-CoM kinematic corrections in associated particle imaging, and explicit alignment of a manipulator pose to an object CoM in grasp planning. Despite these differences, the common aim is to separate intrinsic dynamics from coordinate artifacts or mechanically induced bias [2503.09470][2105.02300][2307.01724].

## 1. Formal scope and common mathematical structure

At the most basic level, CoM frame fixing begins with the decomposition of motion into collective and relative components. For a two-body system with positions \(\mathbf x_1,\mathbf x_2\), masses \(m_1,m_2\), CoM coordinate \(\mathbf R\), and relative coordinate \(\mathbf r\),
\[
\mathbf R=\frac{m_1\mathbf x_1+m_2\mathbf x_2}{m_1+m_2},
\qquad
\mathbf r=\mathbf x_1-\mathbf x_2,
\]
and the equations of motion separate into
\[
(m_1+m_2)\ddot{\mathbf R}=0
\]
together with the relative dynamics. In an exact CoM frame, \(\mathbf R^{(\mathrm{cm})}\equiv 0\), so only intrinsic motion remains [2009.04920][1408.4128].

In many-body and field-theoretic settings the same idea is expressed through symmetry generators rather than elementary coordinates. In nuclear DFT the issue is that the trial Slater determinant or Bogoliubov vacuum \(|\Phi\rangle\) is not an eigenstate of the total momentum operator \(\hat P=\sum_{k=1}^A \hat p_k\), so translational invariance is broken. In numerical relativity at future null infinity, the relevant object is the CoM charge
\[
G^i(u)\equiv \frac{K^i(u)+uP^i(u)}{P^t(u)},
\]
whose nonzero intercept or slope signals residual translation or boost freedom in the waveform’s BMS frame [2503.09470][2208.04356].

A concise cross-domain summary is:

| Domain | Quantity being fixed | Representative operation |
|---|---|---|
| Nuclear EDF/DFT | Total momentum of the many-body state | Projection with \(\hat P_{\mathrm{CM}}=\int d^3a\,U(a)\) |
| Numerical relativity | CoM charge \(G^i(u)\) at \(\mathscr I^+\) | Boost/translation fit and BMS transform |
| Space instrumentation | Offset \(d=[d_x,d_y,d_z]^T\) between sensor and spacecraft COM | EKF–RTS smoothing with \(\chi^2\) outlier rejection |
| API reconstruction | Lab-frame effect of nonzero \(\vec v_{\rm com}\) | Kinematic correction using average CoM velocity |
| Grasp planning | Misalignment between gripper origin and object CoM proxy | Closed-form translation \(t^*=\bar p_o-\bar p_g\) |

These examples suggest that “frame fixing” is not a single algorithm but a structural requirement: one must either impose the correct CoM frame, estimate it, or transform to it before interpreting observables.

## 2. Translational symmetry restoration in nuclear many-body theory

In nuclear energy density functional theory, translational symmetry breaking is universal: unlike rotational or gauge symmetry breaking, it occurs for all nuclei. The isolated nucleus must have zero total momentum, yet the standard mean-field energy
\[
E_{MF}=\langle\Phi|H|\Phi\rangle
\]
contains spurious CoM kinetic energy because \(|\Phi\rangle\) is not an eigenstate of \(\hat P\). The intrinsic state is obtained by projection,
\[
|\Psi_{\rm intr}\rangle \equiv \hat P_{CM}|\Phi\rangle,
\qquad
\hat P_{CM}=\int d^3a\,U(a),
\]
with \(U(a)=\exp[(i/\hbar)\hat P\cdot a]\). In coordinate representation,
\[
\Psi(r_1,\ldots,r_A)=\int d^3a\,\Phi(r_1+a,\ldots,r_A+a),
\]
so the projected wave function is translationally invariant by construction, with \(\langle\Psi_{\rm intr}|\hat P|\Psi_{\rm intr}\rangle=0\) and \([\hat P_{CM},H]=0\) [2503.09470].

The associated energy correction is
\[
\Delta E_{CM}\equiv E_0-E_{MF}\le 0,
\qquad
E_0=\frac{\langle\Psi_{\rm intr}|H|\Psi_{\rm intr}\rangle}{\langle\Psi_{\rm intr}|\Psi_{\rm intr}\rangle}.
\]
Using the generator-coordinate form,
\[
E_0=
\frac{\int d^3a\,\langle\Phi(a)|H|\Phi(0)\rangle}
{\int d^3a\,\langle\Phi(a)|\Phi(0)\rangle}.
\]
A Gaussian overlap approximation, valid up to a few fm of shift, recovers the usual renormalized kinetic-energy prescriptions at lowest order, but the full GCM integral captures both kinetic and interaction contributions to \(\Delta E_{CM}\) [2503.09470].

Using the SeaLL1 functional and exact CoM projection, the correction is large across the nuclear chart:

| Nucleus | \(\Delta E_{CM}\) (MeV) |
|---|---:|
| \(^{16}\mathrm O\) | \(-10.05\) |
| \(^{40}\mathrm{Ca}\) | \(-8.95\) |
| \(^{90}\mathrm{Zr}\) | \(-8.46\) |
| \(^{208}\mathrm{Pb}\) | \(-7.36\) |

These values exceed the Bethe–Weizsäcker mass-formula RMS error of about \(3.5\) MeV and also exceed the typical \(2\)–\(3\) MeV RMS deviations of uncorrected DFT mass fits. The same work reports slight but systematic changes in rms radii, of order \(0.002\)–\(0.06\) fm, comparable to experimental uncertainties. Implementation requires overlap and Hamiltonian kernels over \(O(10\)–\(30)\) shifts and increases cost by a modest factor of a few relative to a single mean-field run, provided box size, mesh handling, and FFTs are treated carefully [2503.09470].

Historically, the method is traced to the CoM projection suggested by Peierls in 1957. In current EDF practice its importance is methodological rather than optional: the spurious translational contribution is numerically larger than the accuracy targets of modern mass models.

## 3. BMS-frame fixing of numerical-relativity waveforms

For gravitational-wave calculations, CoM frame fixing is a problem in asymptotic symmetry rather than ordinary mechanics. Numerical waveforms at future null infinity carry the freedom of the Bondi–van der Burg–Metzner–Sachs group, so comparison among post-Newtonian, numerical-relativity, and perturbative waveforms requires a common BMS frame. Early practice often used a Newtonian CoM trajectory built from horizon coordinates,
\[
\vec x_{\rm CoM}(t)=\frac{m_a\vec x_a(t)+m_b\vec x_b(t)}{M},
\]
followed by a least-squares fit for a translation \(\vec\alpha\) and boost \(\vec\beta\). That approach was shown to be less effective than previously thought, because coordinate black-hole motion need not track the true asymptotic CoM motion, especially for high mass ratio or precessing runs [2105.02300].

Charge-based methods replace bulk gauge-dependent coordinates with asymptotic quantities extracted from the waveform itself. In the Moreschi–Boyle convention one defines the mass aspect \(m\), Lorentz aspect \(N\), and energy-moment aspect \(E\), from which the Poincaré charges \(P_\Psi\), \(J_\Psi\), \(K_\Psi\), and \(E_\Psi\) are constructed by sphere integrals. The CoM charge
\[
G^i(u)=\frac{E^i(u)}{P^t(u)}=\frac{K^i(u)+uP^i(u)}{P^t(u)}
\]
is then fit over an inspiral window. In the linearized method one writes
\[
G^i(u)\approx g_0^i+v^i u,
\]
takes the boost parameter as \(\beta^i=v^i\), and the translation as \(a^i=-g_0^i\), iterating until residuals are small. The same framework is integrated with supertranslation fixing through the Moreschi supermomentum, enabling a complete BMS-frame specification [2208.04356].

On a set of 13 binary black-hole systems, the asymptotic-charge CoM fix reduced residual CoM charge from typical \(|\vec G|\sim 10^{-3}\)–\(10^{-2}\) to \(<10^{-5}\). It also reduced leakage at twice the \((2,2)\) frequency in the \((2,1)\) mode by a factor of \(\sim 10\)–\(100\), lowered mode-by-mode mismatches from \(\mathcal O(10^{-3})\) to \(\mathcal O(10^{-6})\) or below, improved NR–PN \(L^2\) alignment errors over a four-orbit window typically by factors \(5\)–\(20\), and restored the correct memory offset in the \((2,0)\) mode [2105.02300]. A separate charge-based BMS-fixing framework reported a method that is \(20\) times faster than previous optimization-based approaches when mapping to the superrest frame [2208.04356].

For quasicircular, nonprecessing binaries, a later refinement replaced the pure linear ansatz with a post-Newtonian model of the boosted CoM charge that captures physical out-spiraling oscillations. Across 20 SXS simulations, the largest improvement in robustness to fitting-window choice was by a factor of \(\sim 25\) for the boost vector and \(\sim 20\) for the translation vector, with the maximum robustness obtained when the window is centered in the inspiral. That method was incorporated into the `scri` frame-fixing workflow for waveforms produced with Cauchy-characteristic evolution [2603.24661].

A persistent misconception in this area is that Poincaré-only alignment or Newtonian CoM tracking is sufficient. The charge-based results indicate otherwise: waveform-intrinsic asymptotic charges are needed to suppress physically spurious mode mixing and to make hybridization with PN data reliable.

## 4. Estimation and correction in sensing and measurement systems

In experimental systems, CoM frame fixing often appears as an estimation problem. For TaiJi-1, the gravitational reference sensor requires the test-mass center to coincide with the satellite’s center of gravity in order to avoid disturbances from angular acceleration and gradient. The relevant frames are the body frame \(\{B\}\), the test-mass frame \(\{T\}\), and the center-of-mass frame \(\{C\}\), all with identical orientation, while the offset is
\[
\Delta r=d=[d_x,d_y,d_z]^T
\]
from spacecraft COM to test-mass COM, expressed in \(\{B\}\). Over the calibration interval the state is treated as constant,
\[
x_k=[d_{x,k},d_{y,k},d_{z,k}]^T,\qquad f(x)\equiv x,
\]
and the measured GRS output obeys
\[
z_k=\tilde A_k d+n_k.
\]
The estimation pipeline combines an Extended Kalman Filter with a Rauch–Tung–Striebel smoother, uses a \(\chi^2\) residual test with \(\gamma=10^{-3}\) and \(N_z=3\) for outlier rejection, and cross-checks the result with nonlinear least squares solved by Levenberg–Marquardt [2307.01724].

The final TaiJi-1 offsets were reported as \(d_x=-189\pm 27~\mu\mathrm m\), \(d_y=638\pm 68~\mu\mathrm m\), and \(d_z=-818\pm 20~\mu\mathrm m\), equivalently \(dx\approx -0.19\) mm, \(dy\approx 0.64\) mm, and \(dz\approx -0.82\) mm. After in-orbit CoM calibration, the modulation peak in the GRS ASD at the applied torque frequency was suppressed, and the residual acceleration noise in the \(0.001\)–\(0.1\) Hz band was reduced by a factor \(\gtrsim 2\)–\(3\) [2307.01724].

Associated Particle Imaging presents a different measurement problem. In the DT fusion reaction used by API, the \(\alpha\) particle and neutron are exactly back-to-back in the CoM frame, but in the laboratory frame the opening angle is slightly less than \(180^\circ\) because the reacting ion has nonzero \(\vec v_{\rm com}\). The CoM velocity is
\[
\vec v_{\rm com}=\frac{m_i\vec v_i}{m_i+m_t},
\]
and in a thick Ti target it varies with depth through the stopping-power equation \(dE_i/dx=-S(E_i)\). The reconstruction therefore uses a cross-section-weighted average CoM speed obtained from stopping powers and DT fusion cross sections [2204.06124].

When the CoM effect is included in API reconstruction, the mean of reconstructed locations becomes a correctable systematic shift or tilt, but the distribution retains an irreducible spread because the CoM velocity fluctuates event by event. Reported consequences include a systematic tilt or shift of up to \(10\)–\(12\) cm at \(1\) m and an irreducible broadening in the beam-axis direction with a 90% interval of order \(5\)–\(6\) cm at \(1\) m [2204.06124]. This distinction between correctable bias and non-correctable spread is central to experimental CoM correction.

## 5. Off-center actuation, control, and manipulation

In control problems, CoM frame fixing may require either moving the reference frame away from the CoM or aligning the controlled body to the CoM of another object. For a UAV with an off-center slung load, Lv et al. formulate the dynamics about the suspension point rather than the UAV CoM. The frames are the inertial frame \(\mathcal I\), quadrotor body frame \(\mathcal B_q\) at the quadrotor CoM, suspension-point frame \(\mathcal B\), and payload frame \(\mathcal B_p\). With \(L=[l_x,l_y,l_z]^T\) the vector from suspension point to quadrotor CoM and \(l_c=[0,0,l]^T\) the vector from suspension point to payload CoM, the inertial positions are
\[
\xi_q=\xi-R_b^iL,\qquad \xi_p=\xi+R_p^i l_c.
\]
The resulting equations of motion show explicit coupling terms induced by the offset \(L\), including the \(M_{\eta 2}\ddot\xi\) term in the attitude dynamics [2601.03386].

The control design is cascaded. In the middle loop, the virtual input \(M_\sigma\ddot\xi\) regulates the swing angle through a locally exponentially stable subsystem. In the inner loop, the torque
\[
\tau_\eta=M_{\eta1}\ddot\eta_{tr}+M_{\eta2}\ddot\xi+C_\eta\dot q+G_\eta-D_\eta
\]
cancels the coupling exactly, without neglecting the \(M_{\eta2}\ddot\xi\) term or treating it as a disturbance. The paper proves local exponential stability by Lyapunov analysis and argues that fixing the frame at the suspension point simplifies swing-angle control while preserving exact treatment of inertial coupling [2601.03386]. This suggests that CoM frame fixing is not always equivalent to “placing the origin at the CoM”; in some systems, the technically correct move is to choose a frame in which the physically relevant coupling becomes controllable.

In robotic grasp planning, DISF uses CoM alignment explicitly as a contact-stability term. Because the true volumetric CoM is unavailable from point clouds, both object and gripper CoM are approximated by centroids,
\[
\mathrm{centroid}(\{y_k\})=\frac1K\sum_{k=1}^K y_k,
\]
with \(\bar p_o=\mathrm{centroid}(O)\) and \(\bar p_g=\mathrm{centroid}(RF+t)\). The translation-refinement stage minimizes
\[
E_{\mathrm{CoM}}(t)=\|\bar p_g+t-\bar p_o\|^2,
\]
whose closed-form solution is
\[
t^*=\bar p_o-\bar p_g.
\]
This step is embedded between rotation optimization and finger-aperture optimization in the sequence RotOpt \(\rightarrow\) TransRefine \(\rightarrow\) FingerOpt [2512.24550].

The reported effect is substantial. In the Known-shape regime, VISF achieved median \(E_{\mathrm{CoM}}\approx 1.2\times 10^{-2}\) m and DISF reduced it to about \(4.5\times 10^{-3}\) m. In the Observed-shape regime, VISF had median misalignment \(\approx 5.4\times 10^{-2}\) m and DISF reduced it to \(\approx 2.1\times 10^{-2}\) m. Across three robots, average success rose from \(67\%\) to \(93\%\) in Known-shape and from \(33\%\) to \(70\%\) in Observed-shape; on a real UR3e with observed point clouds, DISF achieved \(8/9\) successes versus \(2/9\) for VISF [2512.24550]. Here CoM fixing functions as a stability prior rather than a symmetry restoration.

## 6. Artificial fixing, geometric interpretation, and broader implications

A recurring failure mode is to hold a body fixed in a frame that is only approximately inertial. Gómez et al. examined the Milky Way–Large Magellanic Cloud system and showed that when the LMC is massive, artificially fixing the Milky Way center of mass biases both the LMC orbit and the phase-space structure of other tracers. For \(M_{\rm LMC}\) in the range \(3\times 10^{10}\,M_\odot\) to \(2.5\times 10^{11}\,M_\odot\), the Milky Way CoM within \(50\) kpc can be displaced by \(\Delta R_1\approx 5\)–\(30\) kpc and \(\Delta V_1\approx 10\)–\(75\) km/s over the last \(0.3\)–\(0.5\) Gyr [1408.4128].

The orbital consequences are large. In a representative model with \(M_{\rm vir}=1.5\times 10^{12}\,M_\odot\) and \(M_{\rm LMC}=10^{11}\,M_\odot\), the fixed-MW treatment yields orbital period \(P\approx 13.1\) Gyr and apocenter \(R_{\rm apo}\approx 1.94\,R_{\rm vir}\), whereas the free-MW treatment gives \(P\approx 6.8\) Gyr and \(R_{\rm apo}\approx 1.2\,R_{\rm vir}\). For the Sagittarius stream, including MW recoil as well as LMC torque reduces the angular separation between apocenters by up to \(\sim 17^\circ\), tilts the orbital plane by \(\sim 9^\circ\), and changes debris predictions by tens of degrees on the sky and dozens of km/s in velocity [1408.4128]. In this context, CoM frame fixing is a requirement for dynamical fidelity, not a bookkeeping preference.

The geometric value of the CoM frame is especially transparent in repulsive Rutherford scattering. In the fixed-target frame the shadow caustic has the universal paraboloidal form
\[
z(\rho)=\frac{\rho^2}{8}-2
\]
in scaled units, with the target at the focus. In the CoM frame the projectile and target each cast their own paraboloidal shadow,
\[
z_p(\rho_p)=\frac{\rho_p^2}{8\eta_2}-2\eta_2,\qquad
z_t(\rho_t)=-\frac{\rho_t^2}{8\eta_1}+2\eta_1,
\]
and the focal points of the two shadows coincide at the CoM itself [2009.04920]. This example gives a precise geometrical meaning to “fixing the CoM frame”: the intrinsic symmetry of the two-body interaction becomes manifest only after the collective coordinate is removed.

Across these domains, a common misconception is that CoM fixing is merely a postprocessing convenience. The evidence points in the opposite direction. In nuclei it removes a \(7\)–\(10\) MeV spurious energy bias; in waveform modeling it suppresses mode mixing and improves PN–NR compatibility; in precision instrumentation it lowers acceleration noise; in Galactic dynamics it changes inferred orbital histories; and in manipulation it alters grasp stability [2503.09470][2105.02300][2307.01724][1408.4128][2512.24550]. A plausible implication is that CoM frame fixing should be regarded as part of model definition whenever translational degrees of freedom are not directly observable but still enter the computation.

Source: https://www.emergentmind.com/topics/center-of-mass-com-frame-fixing