REVIEW
Gauge Freedom in Orbital Mechanics
Astrophysics · Submitted 3 Mar 2006 · 29 references
Pith's verdict: A solid, self-contained review of gauge freedom in orbital mechanics; the central incompatibility proof is sound, but the abstract overstates it and the numerics are deferred.
Why Pith says this
Highlight any text to comment on it.
Referee report
The paper
The paper claims
The variation-of-parameters method in orbital mechanics has a gauge freedom, and under velocity-dependent perturbations osculation and canonicity cannot coexist.
The reading
What carries the argument
The load-bearing object is the gauge function $\Phi(t,C_1,\ldots,C_6)$, the convective term $\Phi = \sum_j(\partial q/\partial C_j)\,\dot C_j$ that separates the physical velocity $\dot q$ from the unperturbed velocity function $g = \partial q/\partial t$. Choosing a functional form for $\Phi$ fixes the gauge; the Lagrange constraint $\Phi=0$ is the osculating choice. The decisive step converts this kinematic identity into a dynamical statement: inserting the canonical equations for the elements turns $\Phi$ into $\partial \Delta H/\partial p$, showing that canonicity fixes the gauge rather than leaving it free. The paper also constructs the gauge-invariant Lagrange-bracket system, Eq. (30), which shows that the maximally simplifying gauge for Hamiltonian-type equations is $\Phi = -\partial \Delta L/\partial \dot{\mathbf r}$, the contact gauge.
Extended reading notes
Core claim
On the paper's own terms, the central discovery is an incompatibility theorem for gauge choice. Starting from an unperturbed problem written in canonical variables, one transforms to elements $Q,P$ that are constants of the undisturbed motion and performs canonical perturbation theory for a disturbance $\Delta H(q,p,t)$. Because the same change of variables is used after perturbation, the convective term $\Phi \equiv \sum_j (\partial q/\partial C_j)\,\dot C_j$ is no longer free: substituting the Hamilton equations for $\dot Q,\dot P$ yields $\Phi = \partial \Delta H/\partial p$ (Eq. 44). Hence, whenever $\Delta H$ depends on momentum, the elements are nonosculating: the instantaneous conics are not tangent to the physical trajectory. The same argument applies to spin dynamics with Andoyer elements, for example when the perturbation is the passage to a noninertial frame. The paper concludes that one must choose between preserving canonicity, accepting contact and nonosculating elements, and preserving osculation, accepting non-Hamiltonian element equations.
Load-bearing premise
The proof's load-bearing premise is that, after the perturbation is switched on, the Delaunay or Andoyer 'constants' remain canonical variables satisfying Hamilton's equations under the same transformation $q=f(Q,P,t)$, $p=\chi(Q,P,t)$; if canonicity is abandoned, the paper itself shows that osculating elements can be restored by adding non-Hamiltonian terms.
Editorial extensions
If this is right
- Whenever the disturbing force depends on velocity or momentum, the standard Delaunay-type equations produce nonosculating elements; keeping the elements osculating forces extra terms that are not part of the disturbing function and that spoil Hamilton's form.
- Contact elements, defined by $\Phi = -\partial \Delta L/\partial \dot{\mathbf r}$, keep the equations simple and canonical and return the correct momentum ($p=g$), but they do not return the physical velocity and are not osculating.
- Nonosculating elements parameterize instantaneous conics that are not tangent (and need not even be coplanar) to the orbit, so features like the nonosculating inclination are not the physical inclination.
- For a satellite around a precessing oblate planet, a Goldreich-type calculation in contact elements approximates the secular parts of the osculating elements to first order in the precession rate, with differences accumulating only at very long timescales.
- Exact trajectories are gauge-invariant, but numerical solutions are not: accumulated numerical error produces gauge drift in the elements, making the choice of gauge a practical numerical issue.
Reading between the lines
- A testable extension: in any Hamiltonian system parameterized by constants of the unperturbed motion, not just orbits, a momentum-dependent perturbation should show the same split, either canonical elements whose conics or surfaces are nonosculating, or osculating elements governed by non-Hamiltonian equations.
- In long numerical integrations of dissipative or relativistic two-body problems, monitoring $\Phi = \sum_j(\partial q/\partial C_j)\dot C_j$ could reveal gauge drift masquerading as secular physical evolution; choosing the contact gauge may suppress that artifact.
- If one wants both osculation and a Hamiltonian structure for drag-like problems, one could enlarge the phase space with auxiliary variables; the paper's incompatibility theorem says this cannot happen in the original six-dimensional element space.
Original abstract
In orbital and attitude dynamics the coordinates and the Euler angles are expressed as functions of the time and six constants called elements. Under disturbance, the constants are endowed with time dependence. The Lagrange constraint is then imposed to guarantee that the functional dependence of the perturbed velocity on the time and constants stays the same as in the undisturbed case. Constants obeying this condition are called osculating elements. The constants chosen to be canonical are called Delaunay elements, in the orbital case, or Andoyer elements, in the spin case. (As some Andoyer elements are time dependent even in the free-spin case, the role of constants is played by their initial values.) The Andoyer and Delaunay sets of elements share a feature not readily apparent: in certain cases the standard equations render them non-osculating. In orbital mechanics, elements furnished by the standard planetary equations are non-osculating when perturbations depend on velocities. To preserve osculation, the equations must be amended with extra terms that are not parts of the disturbing function. In the case of Delaunay parameterisation, these terms destroy canonicity. So under velocity-dependent disturbances, osculation and canonicity are incompatible. (Efroimsky and Goldreich 2003, 2004) Similarly, the Andoyer elements turn out to be non-osculating under angular-velocity-dependent perturbation. Amendment of only the Hamiltonian makes the equations render nonosculating elements. To make them osculating, more terms must enter the equations (and the equations will no longer be canonical). In practical calculations, is often convenient to deliberately deviate from osculation by substituting the Lagrange constraint with a condition that gives birth to a family of nonosculating elements.
The compiled record
The paper's own sentences, arranged, kept beside the review and separate from it. Each label is drawn from a closed vocabulary; the number beside it is the probability the classifier gave that label. Numbers in the tables are the authors', checked by code. No sentence in this record was written by a model.
The abstract, sentence by sentence
In orbital and attitude dynamics the coordinates and the Euler angles are expressed as functions of the time and six constants called elements. Under disturbance, the constants are endowed with time dependence. The Lagrange constraint is then imposed to guarantee that the functional dependence of the perturbed velocity on the time and constants stays the same as in the undisturbed case. Constants obeying this condition are called osculating elements. The constants chosen to be canonical are called Delaunay elements, in the orbital case, or Andoyer elements, in the spin case. (As some Andoyer elements are time dependent even in the free-spin case, the role of constants is played by their initial values.) The Andoyer and Delaunay sets of elements share a feature not readily apparent: in certain cases the standard equations render them non-osculating. In orbital mechanics, elements furnished by the standard planetary equations are non-osculating when perturbations depend on velocities. To preserve osculation, the equations must be amended with extra terms that are not parts of the disturbing function. In the case of Delaunay parameterisation, these terms destroy canonicity. So under velocity-dependent disturbances, osculation and canonicity are incompatible. (Efroimsky and Goldreich 2003, 2004) Similarly, the Andoyer elements turn out to be non-osculating under angular-velocity-dependent perturbation. Amendment of only the Hamiltonian makes the equations render nonosculating elements. To make them osculating, more terms must enter the equations (and the equations will no longer be canonical). In practical calculations, is often convenient to deliberately deviate from osculation by substituting the Lagrange constraint with a condition that gives birth to a family of nonosculating elements.
background contribution method result interpretation main claim
| # | sentence | role | the authors' own work | states the main conclusion | strength |
|---|---|---|---|---|---|
| s0 | In orbital and attitude dynamics the coordinates and the Euler angles are expressed as functions of the time and six constants called elements. | background 1.00 | 0.05 | 0.02 | asserted 1.00 |
| s1 | Under disturbance, the constants are endowed with time dependence. | background 1.00 | 0.07 | 0.02 | asserted 1.00 |
| s2 | The Lagrange constraint is then imposed to guarantee that the functional dependence of the perturbed velocity on the time and constants stays the same as in the undisturbed case. | background 0.98 | 0.08 | 0.03 | asserted 1.00 |
| s3 | Constants obeying this condition are called osculating elements. | background 1.00 | 0.08 | 0.03 | asserted 1.00 |
| s4 | The constants chosen to be canonical are called Delaunay elements, in the orbital case, or Andoyer elements, in the spin case. | background 1.00 | 0.09 | 0.03 | asserted 1.00 |
| s5 | (As some Andoyer elements are time dependent even in the free-spin case, the role of constants is played by their initial values.) The Andoyer and Delaunay sets of elements share a feature not readily apparent: in certain cases the standard equations render them non-osculating. | result 0.44 | 0.43 | 0.25 | asserted 0.88 |
| s6 | In orbital mechanics, elements furnished by the standard planetary equations are non-osculating when perturbations depend on velocities. | result 0.79 | 0.32 | 0.19 | asserted 1.00 |
| s7 | To preserve osculation, the equations must be amended with extra terms that are not parts of the disturbing function. | result 0.70 | 0.40 | 0.13 | asserted 1.00 |
| s8 | In the case of Delaunay parameterisation, these terms destroy canonicity. | result 0.86 | 0.47 | 0.15 | asserted 1.00 |
| s9 | So under velocity-dependent disturbances, osculation and canonicity are incompatible. | result 0.73 | 0.53 | 0.58 | asserted 1.00 |
| s10 | (Efroimsky and Goldreich 2003, 2004) Similarly, the Andoyer elements turn out to be non-osculating under angular-velocity-dependent perturbation. | result 0.84 | 0.36 | 0.12 | asserted 0.99 |
| s11 | Amendment of only the Hamiltonian makes the equations render nonosculating elements. | result 0.89 | 0.44 | 0.09 | asserted 1.00 |
| s12 | To make them osculating, more terms must enter the equations (and the equations will no longer be canonical). | result 0.75 | 0.47 | 0.14 | asserted 0.99 |
| s13 | In practical calculations, is often convenient to deliberately deviate from osculation by substituting the Lagrange constraint with a condition that gives birth to a family of nonosculating elements. | interpretation 0.70 | 0.14 | 0.13 | hedged possibility 0.53 |
What was examined, and what was not
- arXiv listing abstract v1
- full text
- figures and tables
- references and citing sentences
- relations to other papers
- P abstract
- S sections
- T body sentences
- C citations
- R cross-paper relations
References
How the literature uses it
2 arXiv papers cite this in 3 sentences
Papers citing it · 2 sentences without a confident label
Cites it as support 1 paper
-
In fact this is a general property: under velocity-dependent disturbances, canonicity and osculation are not compatible [ref] .
Elementary derivation of the Lense-Thirring precession · Concluding remarks, Skovoroda's principle and all that · 0.97
Discussion
Disagree with the review, or know something it missed? Leave a comment, highlight any passage above to comment on it, or challenge this review formally.