Lie-transform Perturbation Theory (LPT)

• Lie-transform Perturbation Theory is a method that uses near-identity transformations to systematically separate fast oscillatory motion from slow secular dynamics in Hamiltonian systems.
• It is formulated both canonically—with scalar generating functions and Poisson brackets—and noncanonically using Lie derivatives on phase-space one-forms, applicable in fields such as plasma physics and celestial mechanics.
• Recent developments include oscillation-center reductions in unmagnetized plasmas, guiding-center formulations addressing gyrogauge invariance, and multi-parameter as well as Lie–Poisson extensions.

Searching arXiv for the cited papers and closely related Lie-transform perturbation theory work. Lie-transform perturbation theory (LPT) is a perturbative framework for constructing near-identity transformations that remove fast dependence from Hamiltonian or Lagrangian dynamical systems while preserving the underlying Poisson or symplectic structure. In plasma physics, celestial mechanics, and Hamiltonian normal-form theory, its central role is to separate fast oscillatory motion from slow secular dynamics order by order in a small parameter. In the sources considered here, LPT appears in canonical form through scalar generating functions and Poisson brackets, in noncanonical form through Lie derivatives acting on phase-space one-forms, and in Lie–Poisson form through derivations on a Lie algebra. Recent work rederives oscillation-center quasilinear theory canonically in an unmagnetized plasma (Brizard, 22 Sep 2025), defends the consistency of standard guiding-center LPT against proposed modifications (Brizard, 2023), analyzes intrinsic singularities in guiding-center Hamiltonian Lie-transform perturbation theory (Lin et al., 2024), generalizes the method to multi-parameter perturbations (Zhang, 2016), and extends the formalism to Hamiltonian systems on Lie algebras (Valvo et al., 2021).

1. Foundational construction

LPT is built from a near-identity transformation of phase-space variables generated by Lie operators. In canonical coordinates $z^\alpha=(x,p,w,t)$ on extended phase space, the canonical Poisson bracket is

$$\{f,g\}=\frac{\partial f}{\partial w}\frac{\partial g}{\partial t}-\frac{\partial f}{\partial t}\frac{\partial g}{\partial w} +\left(\nabla f\cdot\frac{\partial g}{\partial p}-\frac{\partial f}{\partial p}\cdot\nabla g\right),$$

with the six-dimensional reduction

$$\{f,g\}=\frac{\partial f}{\partial x}\cdot\frac{\partial g}{\partial p} -\frac{\partial f}{\partial p}\cdot\frac{\partial g}{\partial x}.$$

A scalar generator $S$ defines the Lie operator $L_S f\equiv \{f,S\}$, and the near-identity transformation is written as

$$T=\exp(\epsilon L_{S_1}+\epsilon^2 L_{S_2}+\cdots).$$

Applied to coordinates, this produces transformed variables such as

$$Z^\alpha=z^\alpha+\epsilon\{S_1,z^\alpha\}+\epsilon^2\{S_2,z^\alpha\} +\frac{\epsilon^2}{2}\{S_1,\{S_1,z^\alpha\}\}+\cdots$$

(Brizard, 22 Sep 2025).

In noncanonical formulations, the generators are vector fields $G_n$, and the transformation is expressed through Lie derivatives: $$z'=\exp\big(\epsilon \mathcal{L}_{G_1}+\epsilon^2\mathcal{L}_{G_2}+\cdots\big)\,z.$$ For a differential one-form $\gamma$, the Lie derivative is

$$\{f,g\}=\frac{\partial f}{\partial w}\frac{\partial g}{\partial t}-\frac{\partial f}{\partial t}\frac{\partial g}{\partial w} +\left(\nabla f\cdot\frac{\partial g}{\partial p}-\frac{\partial f}{\partial p}\cdot\nabla g\right),$$0

so the transformed one-form is obtained by pullback. In one common notation,

$$\{f,g\}=\frac{\partial f}{\partial w}\frac{\partial g}{\partial t}-\frac{\partial f}{\partial t}\frac{\partial g}{\partial w} +\left(\nabla f\cdot\frac{\partial g}{\partial p}-\frac{\partial f}{\partial p}\cdot\nabla g\right),$$1

with perturbative expansions involving $$\{f,g\}=\frac{\partial f}{\partial w}\frac{\partial g}{\partial t}-\frac{\partial f}{\partial t}\frac{\partial g}{\partial w} +\left(\nabla f\cdot\frac{\partial g}{\partial p}-\frac{\partial f}{\partial p}\cdot\nabla g\right),$$2 order by order (Brizard, 2023, Zhang, 2016).

This dual canonical/noncanonical description is not a contradiction but a change of geometric language. In canonical phase spaces one may identify generators with Hamiltonian vector fields $$\{f,g\}=\frac{\partial f}{\partial w}\frac{\partial g}{\partial t}-\frac{\partial f}{\partial t}\frac{\partial g}{\partial w} +\left(\nabla f\cdot\frac{\partial g}{\partial p}-\frac{\partial f}{\partial p}\cdot\nabla g\right),$$3, whereas in noncanonical phase spaces the Lie derivative along $$\{f,g\}=\frac{\partial f}{\partial w}\frac{\partial g}{\partial t}-\frac{\partial f}{\partial t}\frac{\partial g}{\partial w} +\left(\nabla f\cdot\frac{\partial g}{\partial p}-\frac{\partial f}{\partial p}\cdot\nabla g\right),$$4 is primary (Brizard, 2023). A plausible implication is that the practical form of LPT depends less on the specific coordinate choice than on whether the perturbative target is the Hamiltonian, the symplectic one-form, or the full Lie–Poisson derivation.

2. Canonical LPT and oscillation-center reduction

A recent canonical application is the Lie-transform derivation of oscillation-center quasilinear theory for an unmagnetized plasma with weak electrostatic waves (Brizard, 22 Sep 2025). The particle Hamiltonian is

$$\{f,g\}=\frac{\partial f}{\partial w}\frac{\partial g}{\partial t}-\frac{\partial f}{\partial t}\frac{\partial g}{\partial w} +\left(\nabla f\cdot\frac{\partial g}{\partial p}-\frac{\partial f}{\partial p}\cdot\nabla g\right),$$5

with spectral representation

$$\{f,g\}=\frac{\partial f}{\partial w}\frac{\partial g}{\partial t}-\frac{\partial f}{\partial t}\frac{\partial g}{\partial w} +\left(\nabla f\cdot\frac{\partial g}{\partial p}-\frac{\partial f}{\partial p}\cdot\nabla g\right),$$6

or, in the monochromatic eikonal form used for explicit formulae,

$$\{f,g\}=\frac{\partial f}{\partial w}\frac{\partial g}{\partial t}-\frac{\partial f}{\partial t}\frac{\partial g}{\partial w} +\left(\nabla f\cdot\frac{\partial g}{\partial p}-\frac{\partial f}{\partial p}\cdot\nabla g\right),$$7

The ordering assumes weak turbulence $$\{f,g\}=\frac{\partial f}{\partial w}\frac{\partial g}{\partial t}-\frac{\partial f}{\partial t}\frac{\partial g}{\partial w} +\left(\nabla f\cdot\frac{\partial g}{\partial p}-\frac{\partial f}{\partial p}\cdot\nabla g\right),$$8, separation between the fast phase and the slow time $$\{f,g\}=\frac{\partial f}{\partial w}\frac{\partial g}{\partial t}-\frac{\partial f}{\partial t}\frac{\partial g}{\partial w} +\left(\nabla f\cdot\frac{\partial g}{\partial p}-\frac{\partial f}{\partial p}\cdot\nabla g\right),$$9, neglect of particle trapping and nonlinear mode coupling, and a homogeneous unmagnetized plasma. The extended Hamiltonian is

$$\{f,g\}=\frac{\partial f}{\partial x}\cdot\frac{\partial g}{\partial p} -\frac{\partial f}{\partial p}\cdot\frac{\partial g}{\partial x}.$$0

and its push-forward defines the oscillation-center Hamiltonian

$$\{f,g\}=\frac{\partial f}{\partial x}\cdot\frac{\partial g}{\partial p} -\frac{\partial f}{\partial p}\cdot\frac{\partial g}{\partial x}.$$1

with

$$\{f,g\}=\frac{\partial f}{\partial x}\cdot\frac{\partial g}{\partial p} -\frac{\partial f}{\partial p}\cdot\frac{\partial g}{\partial x}.$$2

$$\{f,g\}=\frac{\partial f}{\partial x}\cdot\frac{\partial g}{\partial p} -\frac{\partial f}{\partial p}\cdot\frac{\partial g}{\partial x}.$$3

The first generator $$\{f,g\}=\frac{\partial f}{\partial x}\cdot\frac{\partial g}{\partial p} -\frac{\partial f}{\partial p}\cdot\frac{\partial g}{\partial x}.$$4 is chosen to eliminate nonresonant $$\{f,g\}=\frac{\partial f}{\partial x}\cdot\frac{\partial g}{\partial p} -\frac{\partial f}{\partial p}\cdot\frac{\partial g}{\partial x}.$$5 oscillations. For a time-independent monochromatic wave, imposing $$\{f,g\}=\frac{\partial f}{\partial x}\cdot\frac{\partial g}{\partial p} -\frac{\partial f}{\partial p}\cdot\frac{\partial g}{\partial x}.$$6 yields

$$\{f,g\}=\frac{\partial f}{\partial x}\cdot\frac{\partial g}{\partial p} -\frac{\partial f}{\partial p}\cdot\frac{\partial g}{\partial x}.$$7

At second order the phase-averaged oscillation-center Hamiltonian is the ponderomotive Hamiltonian

$$\{f,g\}=\frac{\partial f}{\partial x}\cdot\frac{\partial g}{\partial p} -\frac{\partial f}{\partial p}\cdot\frac{\partial g}{\partial x}.$$8

while $$\{f,g\}=\frac{\partial f}{\partial x}\cdot\frac{\partial g}{\partial p} -\frac{\partial f}{\partial p}\cdot\frac{\partial g}{\partial x}.$$9 removes second harmonics through

$S$0

(Brizard, 22 Sep 2025).

In the time-dependent eikonal case, the separation between resonant and nonresonant dynamics is made explicit through a window function $S$1. The resonant first-order oscillation-center Hamiltonian is

$S$2

and the nonresonant generator satisfies

$S$3

Using the Bateman–Kruskal integrating-factor method, the solution through $S$4 time derivatives is

$S$5

which vanishes at exact resonance since $S$6 faster than the denominator (Brizard, 22 Sep 2025).

This construction shows the canonical logic of LPT with particular clarity: the transformed Hamiltonian contains the nonresonant physics in an averaged potential and leaves the resonant singular structure outside that average. This suggests that, in oscillation-center theory, the separation between coherent forcing and irreversible transport is not imposed phenomenologically but encoded in the choice of generators.

3. Quasilinear transport, dispersion, and conservation laws

The oscillation-center Vlasov hierarchy follows from the pull-back $S$7 and the transformed time-evolution operator. Up to second order,

$S$8

leading to

$S$9

Without resonances, $L_S f\equiv \{f,S\}$0, so the reduced dynamics is purely ponderomotive and no diffusion occurs. In the quasilinear regime, the spatially averaged second-order equation reduces to

$L_S f\equiv \{f,S\}$1

(Brizard, 22 Sep 2025).

The quasilinear diffusion tensor in momentum space is

$L_S f\equiv \{f,S\}$2

and in velocity space

$L_S f\equiv \{f,S\}$3

For a general electrostatic spectrum,

$L_S f\equiv \{f,S\}$4

with the one-dimensional reduction

$L_S f\equiv \{f,S\}$5

The associated Fokker–Planck equation is

$L_S f\equiv \{f,S\}$6

(Brizard, 22 Sep 2025).

The distinction between nonresonant principal-value effects and resonant diffusion is organized by the Plemelj formula,

$L_S f\equiv \{f,S\}$7

which splits the ponderomotive potential from the resonant contribution. For a general spectrum,

$L_S f\equiv \{f,S\}$8

understood as a principal-value expression near resonance (Brizard, 22 Sep 2025).

The same formalism produces the linear response. The first-order particle distribution is

$L_S f\equiv \{f,S\}$9

and substitution into Poisson’s equation gives

$$T=\exp(\epsilon L_{S_1}+\epsilon^2 L_{S_2}+\cdots).$$0

This separates into $$T=\exp(\epsilon L_{S_1}+\epsilon^2 L_{S_2}+\cdots).$$1 and the amplitude equation

$$T=\exp(\epsilon L_{S_1}+\epsilon^2 L_{S_2}+\cdots).$$2

(Brizard, 22 Sep 2025).

Energy–momentum conservation becomes canonical and explicit. Defining

$$T=\exp(\epsilon L_{S_1}+\epsilon^2 L_{S_2}+\cdots).$$3

the second-order wave energy and momentum are

$$T=\exp(\epsilon L_{S_1}+\epsilon^2 L_{S_2}+\cdots).$$4

with exact balance

$$T=\exp(\epsilon L_{S_1}+\epsilon^2 L_{S_2}+\cdots).$$5

(Brizard, 22 Sep 2025).

4. Guiding-center LPT, gyrogauge structure, and the ordering controversy

For magnetized charged-particle motion, LPT is commonly formulated on the phase-space one-form rather than directly on a canonical Hamiltonian. Standard guiding-center variables are

$$T=\exp(\epsilon L_{S_1}+\epsilon^2 L_{S_2}+\cdots).$$6

with magnetic field $$T=\exp(\epsilon L_{S_1}+\epsilon^2 L_{S_2}+\cdots).$$7 and scale separation

$$T=\exp(\epsilon L_{S_1}+\epsilon^2 L_{S_2}+\cdots).$$8

A key point emphasized in the 2023 Comment is that the physical ratio $$T=\exp(\epsilon L_{S_1}+\epsilon^2 L_{S_2}+\cdots).$$9 is not itself the perturbation parameter; instead, a dimensionless ordering parameter $$Z^\alpha=z^\alpha+\epsilon\{S_1,z^\alpha\}+\epsilon^2\{S_2,z^\alpha\} +\frac{\epsilon^2}{2}\{S_1,\{S_1,z^\alpha\}\}+\cdots$$0 must be introduced explicitly, for example via mass renormalization $$Z^\alpha=z^\alpha+\epsilon\{S_1,z^\alpha\}+\epsilon^2\{S_2,z^\alpha\} +\frac{\epsilon^2}{2}\{S_1,\{S_1,z^\alpha\}\}+\cdots$$1 or charge renormalization $$Z^\alpha=z^\alpha+\epsilon\{S_1,z^\alpha\}+\epsilon^2\{S_2,z^\alpha\} +\frac{\epsilon^2}{2}\{S_1,\{S_1,z^\alpha\}\}+\cdots$$2 (Brizard, 2023).

Under charge renormalization, the guiding-center one-form is

$$Z^\alpha=z^\alpha+\epsilon\{S_1,z^\alpha\}+\epsilon^2\{S_2,z^\alpha\} +\frac{\epsilon^2}{2}\{S_1,\{S_1,z^\alpha\}\}+\cdots$$3

while under mass renormalization

$$Z^\alpha=z^\alpha+\epsilon\{S_1,z^\alpha\}+\epsilon^2\{S_2,z^\alpha\} +\frac{\epsilon^2}{2}\{S_1,\{S_1,z^\alpha\}\}+\cdots$$4

with $$Z^\alpha=z^\alpha+\epsilon\{S_1,z^\alpha\}+\epsilon^2\{S_2,z^\alpha\} +\frac{\epsilon^2}{2}\{S_1,\{S_1,z^\alpha\}\}+\cdots$$5. In both orderings the one-form is gyroangle-independent to the chosen truncation order, and the canonical gyroaction $$Z^\alpha=z^\alpha+\epsilon\{S_1,z^\alpha\}+\epsilon^2\{S_2,z^\alpha\} +\frac{\epsilon^2}{2}\{S_1,\{S_1,z^\alpha\}\}+\cdots$$6 is invariant to that order (Brizard, 2023).

The gyrogauge vector field is

$$Z^\alpha=z^\alpha+\epsilon\{S_1,z^\alpha\}+\epsilon^2\{S_2,z^\alpha\} +\frac{\epsilon^2}{2}\{S_1,\{S_1,z^\alpha\}\}+\cdots$$7

where $$Z^\alpha=z^\alpha+\epsilon\{S_1,z^\alpha\}+\epsilon^2\{S_2,z^\alpha\} +\frac{\epsilon^2}{2}\{S_1,\{S_1,z^\alpha\}\}+\cdots$$8, and the combination

$$Z^\alpha=z^\alpha+\epsilon\{S_1,z^\alpha\}+\epsilon^2\{S_2,z^\alpha\} +\frac{\epsilon^2}{2}\{S_1,\{S_1,z^\alpha\}\}+\cdots$$9

ensures gyrogauge invariance. The noncanonical bracket is

$G_n$0

with

$G_n$1

and

$G_n$2

(Brizard, 2023).

The controversy addressed by the Comment concerns a proposal that standard LPT requires modification because the gyromotion term and the parallel symplectic term should appear at the same order (Zheng, 2023). The Comment rejects this claim and argues that the apparent inconsistency arises from mixing different renormalization orderings and from not fixing the perturbation parameter $G_n$3 consistently (Brizard, 2023). It further states that the gyrogauge-invariant combination requires the relative ordering

$G_n$4

if one writes the $G_n$5-sector with a separate parameter $G_n$6; setting $G_n$7 breaks gyrogauge invariance. On this reading, the gyromotion term is one order higher than the parallel symplectic term, not at the same order (Brizard, 2023).

By contrast, the modification paper argues that conventional derivations illegitimately commute the parametric limit $G_n$8 with differentiation when the gyrophase is fast. It proposes the altered Lie action on a one-form

$G_n$9

and asserts that keeping the last term restores the $$z'=\exp\big(\epsilon \mathcal{L}_{G_1}+\epsilon^2\mathcal{L}_{G_2}+\cdots\big)\,z.$$0 contribution at first order, leading directly to

$$z'=\exp\big(\epsilon \mathcal{L}_{G_1}+\epsilon^2\mathcal{L}_{G_2}+\cdots\big)\,z.$$1

(Zheng, 2023).

The sources therefore document an explicit methodological dispute rather than a settled consensus. One source claims the standard framework is internally consistent provided the ordering is fixed and gyrogauge invariance is respected (Brizard, 2023); another claims the transformation rule itself must be altered because of heterogeneous differential ordering (Zheng, 2023). A cautious synthesis is that the disagreement centers not on the usefulness of Lie transforms for guiding-center theory, but on how the ordering is encoded in the one-form calculus.

5. Singular Hamiltonian LPT and multi-parameter generalizations

A separate development concerns singularity in the determining equations of Hamiltonian Lie-transform perturbation theory for guiding-center motion. In the differential-form formulation, the Poincaré–Cartan one-form is expanded as

$$z'=\exp\big(\epsilon \mathcal{L}_{G_1}+\epsilon^2\mathcal{L}_{G_2}+\cdots\big)\,z.$$2

and successive pullbacks generate transformed forms. At order $$z'=\exp\big(\epsilon \mathcal{L}_{G_1}+\epsilon^2\mathcal{L}_{G_2}+\cdots\big)\,z.$$3, the linear relation

$$z'=\exp\big(\epsilon \mathcal{L}_{G_1}+\epsilon^2\mathcal{L}_{G_2}+\cdots\big)\,z.$$4

shows that the determining map for $$z'=\exp\big(\epsilon \mathcal{L}_{G_1}+\epsilon^2\mathcal{L}_{G_2}+\cdots\big)\,z.$$5 is controlled by the unperturbed two-form $$z'=\exp\big(\epsilon \mathcal{L}_{G_1}+\epsilon^2\mathcal{L}_{G_2}+\cdots\big)\,z.$$6. Because extended phase space is odd-dimensional, $$z'=\exp\big(\epsilon \mathcal{L}_{G_1}+\epsilon^2\mathcal{L}_{G_2}+\cdots\big)\,z.$$7 has null vectors, so the linear problem is intrinsically singular (Lin et al., 2024).

The 2024 scheme resolves this by splitting the phase space into null and transverse subspaces and using a staggered determination of generators. The oscillatory gauge functions are constructed exactly by explicit gyro-angle integration: $$z'=\exp\big(\epsilon \mathcal{L}_{G_1}+\epsilon^2\mathcal{L}_{G_2}+\cdots\big)\,z.$$8 Because $$z'=\exp\big(\epsilon \mathcal{L}_{G_1}+\epsilon^2\mathcal{L}_{G_2}+\cdots\big)\,z.$$9 is invertible on the oscillatory subspace, the intrinsic singularity is removed by projecting onto $\gamma$0 and applying $\gamma$1. The paper argues that this avoids the uncontrolled higher-order errors that arise when gauge functions are only approximately obtained from PDEs along nontrivial null flows (Lin et al., 2024).

This singular analysis leads to a unified guiding-center formalism retaining strong $\gamma$2 shear and electromagnetic fluctuations. The transformed one-form through order $\gamma$3 has the structure

$\gamma$4

which exposes the invariant $\gamma$5 and the strong-shear corrections (Lin et al., 2024).

Another generalization addresses multiple small parameters. The 2016 multi-parameter method extends the classical single-parameter Cary–Littlejohn construction to a one-form pullback with

$\gamma$6

provided the parameters are independent (Zhang, 2016). The paper then introduces a relaxed formalism in which the parameter vector may include combinations such as $\gamma$7, with generators selected during the cancellation procedure rather than derived from autonomous PDEs.

Applied to electrostatic gyrokinetics, this method chooses generators to eliminate gyroangle dependence directly in the Lagrangian one-form. The lowest-order spatial generator is

$\gamma$8

where

$\gamma$9

and a magnetic-moment generator satisfies

$$\{f,g\}=\frac{\partial f}{\partial w}\frac{\partial g}{\partial t}-\frac{\partial f}{\partial t}\frac{\partial g}{\partial w} +\left(\nabla f\cdot\frac{\partial g}{\partial p}-\frac{\partial f}{\partial p}\cdot\nabla g\right),$$00

The resulting reduced one-form is

$$\{f,g\}=\frac{\partial f}{\partial w}\frac{\partial g}{\partial t}-\frac{\partial f}{\partial t}\frac{\partial g}{\partial w} +\left(\nabla f\cdot\frac{\partial g}{\partial p}-\frac{\partial f}{\partial p}\cdot\nabla g\right),$$01

and the paper states that finite-Larmor-radius terms are cancelled algebraically without using any gauge function (Zhang, 2016).

6. Lie-algebraic extensions, scope, and limitations

LPT also extends beyond canonical or noncanonical particle phase spaces to Hamiltonian systems posed on Lie algebras. In the 2021 formulation, the basic objects are observables in a Lie algebra $$\{f,g\}=\frac{\partial f}{\partial w}\frac{\partial g}{\partial t}-\frac{\partial f}{\partial t}\frac{\partial g}{\partial w} +\left(\nabla f\cdot\frac{\partial g}{\partial p}-\frac{\partial f}{\partial p}\cdot\nabla g\right),$$02 with bracket $$\{f,g\}=\frac{\partial f}{\partial w}\frac{\partial g}{\partial t}-\frac{\partial f}{\partial t}\frac{\partial g}{\partial w} +\left(\nabla f\cdot\frac{\partial g}{\partial p}-\frac{\partial f}{\partial p}\cdot\nabla g\right),$$03, or equivalently a Lie–Poisson algebra. A Hamiltonian flow is generated by the inner derivation $$\{f,g\}=\frac{\partial f}{\partial w}\frac{\partial g}{\partial t}-\frac{\partial f}{\partial t}\frac{\partial g}{\partial w} +\left(\nabla f\cdot\frac{\partial g}{\partial p}-\frac{\partial f}{\partial p}\cdot\nabla g\right),$$04, and for non-autonomous systems one works with

$$\{f,g\}=\frac{\partial f}{\partial w}\frac{\partial g}{\partial t}-\frac{\partial f}{\partial t}\frac{\partial g}{\partial w} +\left(\nabla f\cdot\frac{\partial g}{\partial p}-\frac{\partial f}{\partial p}\cdot\nabla g\right),$$05

A near-identity Lie transform generated by $$\{f,g\}=\frac{\partial f}{\partial w}\frac{\partial g}{\partial t}-\frac{\partial f}{\partial t}\frac{\partial g}{\partial w} +\left(\nabla f\cdot\frac{\partial g}{\partial p}-\frac{\partial f}{\partial p}\cdot\nabla g\right),$$06 acts on observables by

$$\{f,g\}=\frac{\partial f}{\partial w}\frac{\partial g}{\partial t}-\frac{\partial f}{\partial t}\frac{\partial g}{\partial w} +\left(\nabla f\cdot\frac{\partial g}{\partial p}-\frac{\partial f}{\partial p}\cdot\nabla g\right),$$07

and on derivations by conjugation

$$\{f,g\}=\frac{\partial f}{\partial w}\frac{\partial g}{\partial t}-\frac{\partial f}{\partial t}\frac{\partial g}{\partial w} +\left(\nabla f\cdot\frac{\partial g}{\partial p}-\frac{\partial f}{\partial p}\cdot\nabla g\right),$$08

(Valvo et al., 2021).

The purpose of the transformation is to preserve a chosen Lie subalgebra $$\{f,g\}=\frac{\partial f}{\partial w}\frac{\partial g}{\partial t}-\frac{\partial f}{\partial t}\frac{\partial g}{\partial w} +\left(\nabla f\cdot\frac{\partial g}{\partial p}-\frac{\partial f}{\partial p}\cdot\nabla g\right),$$09 up to a given order. Using a projector $$\{f,g\}=\frac{\partial f}{\partial w}\frac{\partial g}{\partial t}-\frac{\partial f}{\partial t}\frac{\partial g}{\partial w} +\left(\nabla f\cdot\frac{\partial g}{\partial p}-\frac{\partial f}{\partial p}\cdot\nabla g\right),$$10 onto $$\{f,g\}=\frac{\partial f}{\partial w}\frac{\partial g}{\partial t}-\frac{\partial f}{\partial t}\frac{\partial g}{\partial w} +\left(\nabla f\cdot\frac{\partial g}{\partial p}-\frac{\partial f}{\partial p}\cdot\nabla g\right),$$11, its complement $$\{f,g\}=\frac{\partial f}{\partial w}\frac{\partial g}{\partial t}-\frac{\partial f}{\partial t}\frac{\partial g}{\partial w} +\left(\nabla f\cdot\frac{\partial g}{\partial p}-\frac{\partial f}{\partial p}\cdot\nabla g\right),$$12, and a pseudo-inverse $$\{f,g\}=\frac{\partial f}{\partial w}\frac{\partial g}{\partial t}-\frac{\partial f}{\partial t}\frac{\partial g}{\partial w} +\left(\nabla f\cdot\frac{\partial g}{\partial p}-\frac{\partial f}{\partial p}\cdot\nabla g\right),$$13 of the homological operator, the first-order homological equation is

$$\{f,g\}=\frac{\partial f}{\partial w}\frac{\partial g}{\partial t}-\frac{\partial f}{\partial t}\frac{\partial g}{\partial w} +\left(\nabla f\cdot\frac{\partial g}{\partial p}-\frac{\partial f}{\partial p}\cdot\nabla g\right),$$14

with $$\{f,g\}=\frac{\partial f}{\partial w}\frac{\partial g}{\partial t}-\frac{\partial f}{\partial t}\frac{\partial g}{\partial w} +\left(\nabla f\cdot\frac{\partial g}{\partial p}-\frac{\partial f}{\partial p}\cdot\nabla g\right),$$15. At second order,

$$\{f,g\}=\frac{\partial f}{\partial w}\frac{\partial g}{\partial t}-\frac{\partial f}{\partial t}\frac{\partial g}{\partial w} +\left(\nabla f\cdot\frac{\partial g}{\partial p}-\frac{\partial f}{\partial p}\cdot\nabla g\right),$$16

The normalized Hamiltonian then takes the form

$$\{f,g\}=\frac{\partial f}{\partial w}\frac{\partial g}{\partial t}-\frac{\partial f}{\partial t}\frac{\partial g}{\partial w} +\left(\nabla f\cdot\frac{\partial g}{\partial p}-\frac{\partial f}{\partial p}\cdot\nabla g\right),$$17

(Valvo et al., 2021).

The application to a non-autonomous symmetric top on $$\{f,g\}=\frac{\partial f}{\partial w}\frac{\partial g}{\partial t}-\frac{\partial f}{\partial t}\frac{\partial g}{\partial w} +\left(\nabla f\cdot\frac{\partial g}{\partial p}-\frac{\partial f}{\partial p}\cdot\nabla g\right),$$18 shows that the standard LPT logic of homological equations, projection onto normal form, and pseudo-inverse solution survives in a noncanonical Lie–Poisson setting. This suggests that LPT is best viewed not as a technique tied to specific coordinates, but as a normal-form procedure on geometric structures that admit derivations, brackets, and a perturbative hierarchy.

The limitations across the cited works are also consistent. The oscillation-center quasilinear derivation assumes weak amplitude, timescale separation, random phases when spectral sums are used, and neglect of trapping and strong nonlinearities (Brizard, 22 Sep 2025). The standard guiding-center Comment assumes $$\{f,g\}=\frac{\partial f}{\partial w}\frac{\partial g}{\partial t}-\frac{\partial f}{\partial t}\frac{\partial g}{\partial w} +\left(\nabla f\cdot\frac{\partial g}{\partial p}-\frac{\partial f}{\partial p}\cdot\nabla g\right),$$19, consistent renormalization ordering, and explicit retention of the gyrogauge term $$\{f,g\}=\frac{\partial f}{\partial w}\frac{\partial g}{\partial t}-\frac{\partial f}{\partial t}\frac{\partial g}{\partial w} +\left(\nabla f\cdot\frac{\partial g}{\partial p}-\frac{\partial f}{\partial p}\cdot\nabla g\right),$$20 (Brizard, 2023). The singular guiding-center scheme assumes smooth fields, $$\{f,g\}=\frac{\partial f}{\partial w}\frac{\partial g}{\partial t}-\frac{\partial f}{\partial t}\frac{\partial g}{\partial w} +\left(\nabla f\cdot\frac{\partial g}{\partial p}-\frac{\partial f}{\partial p}\cdot\nabla g\right),$$21, low-frequency small-amplitude fluctuations, and local rather than global geometric treatment (Lin et al., 2024). The multi-parameter gyrokinetic construction retains only terms up to $$\{f,g\}=\frac{\partial f}{\partial w}\frac{\partial g}{\partial t}-\frac{\partial f}{\partial t}\frac{\partial g}{\partial w} +\left(\nabla f\cdot\frac{\partial g}{\partial p}-\frac{\partial f}{\partial p}\cdot\nabla g\right),$$22 with $$\{f,g\}=\frac{\partial f}{\partial w}\frac{\partial g}{\partial t}-\frac{\partial f}{\partial t}\frac{\partial g}{\partial w} +\left(\nabla f\cdot\frac{\partial g}{\partial p}-\frac{\partial f}{\partial p}\cdot\nabla g\right),$$23 and is restricted to electrostatic perturbations in the form presented (Zhang, 2016). The Lie-algebraic KAM-like construction requires nonresonance, analyticity, and smallness conditions (Valvo et al., 2021).

Taken together, these works define LPT as a family of near-identity, order-by-order normal-form constructions rather than a single formula. Its common invariant features are preservation of bracket structure, systematic removal of fast dependence, and separation of secular from oscillatory dynamics. Its technical diversity arises from the geometry of the problem: canonical Hamiltonians in oscillation-center theory, noncanonical one-forms in guiding-center and gyrokinetic reductions, singular null structures in differential-form Hamiltonian perturbation theory, and Lie–Poisson derivations in algebraic normal-form problems.