Non-Gravitational Acceleration Mechanisms

• Non-gravitational acceleration mechanisms are physical processes beyond gravity, including outgassing, radiation pressure, and thermal emission.
• They are quantified using high-precision astrometry and in-situ accelerometry combined with robust thermophysical and dynamical modeling.
• Applications range from understanding comet and asteroid trajectories to ensuring accurate spacecraft navigation and probing cosmic particle acceleration.

Non-gravitational acceleration (NGA) mechanisms are physical processes that impart measurable accelerations to natural and artificial bodies beyond those predicted by Newtonian gravitation. In the context of planetary systems, these accelerations manifest in the trajectories and spin states of comets, asteroids, meteoroids, and spacecraft. Key sources include the recoil of anisotropically ejected mass (e.g., cometary outgassing), photon momentum transfer (radiation pressure), anisotropic thermal emission (Yarkovsky and YORP effects), and, for charged dust grains, Lorentz forces in magnetic fields. Robust detection and quantification of NGA require both high-precision astrometry (or direct in-situ accelerometry) and physically consistent thermophysical or dynamical models that connect the observed accelerations to underlying microphysical properties and environmental conditions.

1. Physical Origins and Principal Mechanisms

The canonical NGA mechanisms in planetary contexts are:

• Cometary Outgassing: Sublimation of volatile ices (primarily H$_2$O, CO$_2$, CO) from the sunlit hemisphere drives a directed gas flow. Momentum conservation requires a recoil force on the nucleus, quantified as

$$\vec{a}_\mathrm{ng} = \frac{1}{M}\sum_i \eta_i \dot{m}_i v_{\rm th,i} A_{\rm act,i}\hat{n}_i,$$

where $\dot{m}_i$ is the mass flux from species $i$, $v_{\rm th,i}$ the thermal outflow speed, $\eta_i$ the collimation factor, and $A_{\rm act,i}$ the area of the active patch (Neukart, 5 Nov 2025).

• Photon Momentum Transfer (Radiation Pressure): Photons impinging on a body impart momentum, with the resulting acceleration

$$a_\mathrm{rad} = \frac{3 Q_\mathrm{pr} L_\odot}{16\pi \rho c a r^2},$$

where $Q_\mathrm{pr}$ is the radiation pressure efficiency, $_2$0 the body radius, and $_2$1 the heliocentric distance (Jewitt, 2024).

• Anisotropic Thermal Emission (Yarkovsky/YORP Effects): Diurnal thermal lags cause a net force from delayed emission of absorbed solar energy, leading to gradual orbital drift (Yarkovsky effect) and spin-state changes (YORP effect) (Kanamaru et al., 2023).
• Outgassing Torques: Non-central gas emission can induce measurable torques that alter spin rates and axis orientations. The outgassing-induced torque is parameterized by a dimensionless lever arm $_2$2, connecting net force and torque as $_2$3, with typical $_2$4 to $_2$5 (Rafikov, 2018).
• Lorentz Forces: For submicron dust, charge accumulation and solar wind magnetic fields cause significant Lorentz acceleration (Jewitt, 2024), but this mechanism is negligible for macroscopic solids.

For active asteroids and spacecraft, analogous processes dominate: outgassing (volatile or non-volatile), solar radiation pressure, and residual thermal or venting forces from spacecraft systems (Magnafico et al., 2024).

2. Mathematical Formalism for NGA in Comets and Asteroids

The orbit-determination community employs the Marsden et al. (1973) “Style II” model:

$_2$6

where $_2$7 are the orbit-tied radial, transverse, and normal unit vectors, $_2$8 are constants fitted from astrometry (in AU day$_2$9), and $$\vec{a}_\mathrm{ng} = \frac{1}{M}\sum_i \eta_i \dot{m}_i v_{\rm th,i} A_{\rm act,i}\hat{n}_i,$$0 is an empirical sublimation law, e.g.,

$$\vec{a}_\mathrm{ng} = \frac{1}{M}\sum_i \eta_i \dot{m}_i v_{\rm th,i} A_{\rm act,i}\hat{n}_i,$$1

with standard parameters for water-ice: $$\vec{a}_\mathrm{ng} = \frac{1}{M}\sum_i \eta_i \dot{m}_i v_{\rm th,i} A_{\rm act,i}\hat{n}_i,$$2, $$\vec{a}_\mathrm{ng} = \frac{1}{M}\sum_i \eta_i \dot{m}_i v_{\rm th,i} A_{\rm act,i}\hat{n}_i,$$3, $$\vec{a}_\mathrm{ng} = \frac{1}{M}\sum_i \eta_i \dot{m}_i v_{\rm th,i} A_{\rm act,i}\hat{n}_i,$$4, $$\vec{a}_\mathrm{ng} = \frac{1}{M}\sum_i \eta_i \dot{m}_i v_{\rm th,i} A_{\rm act,i}\hat{n}_i,$$5 AU (Sekanina, 2021). More complex thermophysical models employ direct calculation of local surface energy balance, Clausius–Clapeyron vapor-pressure laws, and Monte Carlo sampling over rotational and spatial heterogeneities (Neukart, 5 Nov 2025).

The magnitude of NGA typically scales inversely with nucleus size ($$\vec{a}_\mathrm{ng} = \frac{1}{M}\sum_i \eta_i \dot{m}_i v_{\rm th,i} A_{\rm act,i}\hat{n}_i,$$6), and the amplitude and orientation of $$\vec{a}_\mathrm{ng} = \frac{1}{M}\sum_i \eta_i \dot{m}_i v_{\rm th,i} A_{\rm act,i}\hat{n}_i,$$7 encode details of activity distribution, rotational state, and anisotropy.

Thermal emission–driven forces for asteroids rely on accurate treatment of the 3D shape and thermal properties (thermal inertia, albedo, rotation period), implementing numerically the recoil force integration across the surface:

$$\vec{a}_\mathrm{ng} = \frac{1}{M}\sum_i \eta_i \dot{m}_i v_{\rm th,i} A_{\rm act,i}\hat{n}_i,$$8

where $$\vec{a}_\mathrm{ng} = \frac{1}{M}\sum_i \eta_i \dot{m}_i v_{\rm th,i} A_{\rm act,i}\hat{n}_i,$$9 is the local emittance and $\dot{m}_i$0 the surface normal (Kanamaru et al., 2023).

3. Observational Diagnostics and Measurement Methodologies

Detection of NGA in natural bodies traditionally exploits long-arc astrometric tracking, isolating deviations from the gravitational trajectory in the orbital elements. For comets (and active asteroids), statistical significance is established when fitted $\dot{m}_i$1, with robust signals available for kilometer-class bodies exhibiting SNR $\dot{m}_i$2 in all components (Hui et al., 2016).

Space-borne accelerometry provides direct in-situ measurements. The combination of three-axis proof-mass accelerometers with bias-rejection via mechanical or electrostatic rotation (e.g., Gravity Advanced Package, Italian Spring Accelerometer) enables unbiased low-frequency NGA determinations at sub-pm s$\dot{m}_i$3 precision (Lenoir et al., 2011, Magnafico et al., 2024). Signal modulation, frequency-domain demodulation, and careful calibration of the instrument bias are essential for absolute measurements.

For spacecraft, transient non-gravitational perturbations (e.g., venting events, thermal plumes) are identified by correlating accelerometer spikes with attitude and orbit control torques. Localization of net-force application points and subsequent elimination of spurious events or model mismatches are integral to mission radio science (Magnafico et al., 2024).

4. Model Degeneracies, Systematics, and Case Studies

Interpretation of observed NGA is subject to intrinsic model degeneracies and systematic uncertainties. For cometary and interstellar objects, both outgassing and radiation pressure can mimic similar radial $\dot{m}_i$4 scaling; dynamical fits alone often cannot discriminate between them (Spada, 2023). For 1I/'Oumuamua, dynamical solutions accommodate either solar radiation pressure (implying ultra-low density or thin-sheet geometry) or recoil from outgassing (implying pure ices, collimated jets, and exceedingly low dust/gas ratios not seen in Solar System analogs) (Katz, 2019, Spada, 2023). For 3I/ATLAS, state-of-the-art thermophysical and Monte Carlo modeling demonstrates that mixed CO/CO$\dot{m}_i$5 outgassing with localized, sub-percent active fractions can fully explain the measured NGA, ruling out the need for nonphysical compositions or area-to-mass ratios (Neukart, 5 Nov 2025). Nevertheless, varying the outgassing law’s functional form and parameter selection (symmetric vs. time-lagged vs. asymmetric power law) can introduce up to $\dot{m}_i$650% uncertainty in nucleus size estimates, underscoring the importance of systematic error exploration (Spada et al., 28 Feb 2026, Forbes et al., 20 Dec 2025).

Spin evolution provides an auxiliary diagnostic. The NGA-induced change in spin rate ($\dot{m}_i$7) depends linearly on the lever arm $\dot{m}_i$8 and the orbital-averaged NGA, calibrated empirically from a small sample to $\dot{m}_i$9 (Rafikov, 2018). The lack of comets displaying extreme spin changes at the uppermost predicted rates suggests that spin-up fission is a dominant evolutionary pathway and that systems with $i$0 are rare or unobservable due to catastrophic breakup.

5. Thermophysical Constraints and Limits on Exotic Mechanisms

Quantitative constraints from thermophysical modeling set stringent limits on viable natural NGA mechanisms. The joint solution of the diurnal and obliquity-averaged energy balance,

$i$1

combined with empirically calibrated vapor-pressure laws and jet collimation factors, tightly bounds the allowed production rates and mass loss (Neukart, 5 Nov 2025). For 3I/ATLAS, active fractions for conventional volatiles (CO/CO$i$2) below 1% suffice for nuclei of $i$3–$i$4 km radius, while less volatile species (NH$i$5, CH$i$6) cannot produce the observed NGA magnitude even with unrealistically large areas or collimation (Neukart, 5 Nov 2025).

Radiation-pressure–driven NGA requires area-to-mass ratios $i$7 orders of magnitude above those supportable by plausible densities ($i$8 kg m$i$9), precluding such models for kilometer-class nuclei (Neukart, 5 Nov 2025). Accordingly, extreme scenarios (e.g., N$v_{\rm th,i}$0 ice–dominated outgassing, light sails) are excluded when standard chemistry and heat-transport constraints are enforced in cometary NGAs.

6. Astrophysical Particle Acceleration: Non-Planetary NGA Mechanisms

In astrophysical and cosmic-ray physics, “non-gravitational” acceleration refers primarily to mechanisms converting macroscopic flows or fields into random kinetic energy of particles. First-order Fermi (diffusive shock acceleration) at supernova shocks, second-order Fermi (stochastic, turbulent scattering), magnetic reconnection, and plasma wakefield effects are the leading microphysical NGA mechanisms (Adamo et al., 2022). Scaling relations (Hillas criterion, efficiency $v_{\rm th,i}$1) and observational diagnostics (multi-messenger detection of cosmic rays, neutrinos, $v_{\rm th,i}$2-rays) allow direct probing of acceleration site parameters.

Mechanism	Max Energy (proton, eV)	Dominant Environment
Fermi-I (shock)	$v_{\rm th,i}$3	SNR, pulsar wind nebulae
Magnetic reconnection	$v_{\rm th,i}$4–$v_{\rm th,i}$5	AGN jets, Crab flares
Wakefield	$v_{\rm th,i}$6	PWN, GRB jets

This domain remains physically distinct from orbital and thermal NGAs, but shares the broad principle: momentum and energy transfer from a macroscopic driver (field, flow, emission) to a body or particle population in a non-gravitational manner.

7. Practical Implications and Future Directions

Accurately modeling and measuring NGA is essential for precise orbit determination, risk assessment (e.g., asteroid impact prediction), spacecraft navigation, and deciphering the evolutionary pathways of small bodies. For interstellar objects and Solar System comets, comprehensive multi-wavelength (including IR and radio for volatile identification), in-situ, and photometric measurements are critical for resolving the microphysical nature of NGA and breaking degeneracies inherent in dynamical-only solutions (Spada, 2023, Neukart, 5 Nov 2025).

Spacecraft missions increasingly employ high-sensitivity accelerometers with advanced signal bias rejection, enabling direct detection and characterization of unmodeled forces during critical maneuvers (Lenoir et al., 2011, Magnafico et al., 2024). Future mitigation strategies must include thorough thermal, mechanical, and outgassing characterization pre-launch to minimize systematic uncertainties in radiometric and accelerometric datasets.

For asteroidal Yarkovsky/YORP modeling, integration of high-resolution shape, roughness, and thermal inertia data—potentially supplied by thermal-IR imagers and in situ landers—will advance quantitative predictions of secular orbital and spin evolution (Kanamaru et al., 2023).

A major challenge persists in propagating systematic uncertainties from thermophysical parameter space into dynamical quantities, motivating continued development of simultaneous orbital-dynamical–thermophysical inversion frameworks, robust error-propagation methodologies, and the acquisition of baseline observations over extended heliocentric and rotational phase coverage (Spada et al., 28 Feb 2026, Forbes et al., 20 Dec 2025).