Λ-symmetries for first-order ODE systems

See the works of Cicogna: @cicognareduction.

Overview

A Λ-symmetry is a generalized notion of symmetry for an autonomous system (system of first order ODEs)
u˙a=fa(u),A=fa(u)∂ua,
generated by a vector field

X=φa(u),∂ua.

1. Definition (bracket form)

Compute the Lie bracket
[X,A]a=φb∂bfa−fb∂bφa.
We say that X is a Λ-symmetry of the system if there exists a matrix field Λab(u) such that
[X,A]a=(Λφ)a=Λabφb.
Equivalently, in vector notation,
[X,A]=ΛX.

Remarks:


2. Geometric meaning: why reduction may fail or be weaker

Classical symmetry reduction uses Frobenius: if [X,A] lies in the distribution generated by X and A (in particular if [X,A] is a linear combination of X and A) then the distribution is involutive and projects to a well-defined reduced vector field on the quotient by the flow of X.

For a general Λ-symmetry, [X,A]=ΛX is not necessarily a linear combination of X and A. Consequently:

Nonetheless, Λ controls precisely how the invariants of X evolve under A, so an algebraic reduction is possible: write the system in coordinates adapted to X (invariants + orbit coordinate) and use the information encoded in Λ to understand the dependence on the orbit coordinate.


3. Adapted coordinates and the reduction recipe (practical)

Step 0 — choose adapted coordinates. Locally choose coordinates (wi,z) such that
X(wi)=0 (i=1,…,n−1),X(z)=1.
Thus the wi are invariants of X and z parametrizes the orbits.

Write the dynamical vector field in these coords:
A=Fi(w,z)∂wi+G(w,z)∂z.

Step 1 — bracket in adapted coordinates. Because X=∂z, the bracket is
[X,A]=[∂z,Fi∂wi+G∂z]=(∂zFi)∂wi+(∂zG)∂z.

Step 2 — impose the Λ-condition. The Λ-equation becomes
∂zFi=(Λφ)wi,∂zG=(Λφ)z.
But in the adapted frame φ=(0,…,0,1)T (last component =1), so the right-hand sides are simply the last column of Λ in these coordinates:
∂zFi=Λwiz,∂zG=Λzz.

Interpretation / reduction criterion.

Practical recipe to attempt reduction:

  1. Pick a candidate generator X.
  2. Compute adapted coordinates (w,z) (solve X(w)=0, X(z)=1 locally).
  3. Compute F(w,z),G(w,z) from f.
  4. Compute ∂zF and read off the last column of Λ. If these vanish (or are simple), the w-subsystem reduces.

4. Formula for Λ in general coordinates

Given X=φa∂ua and A=fa∂ua, compute Ba:=[X,A]a. Then the defining linear system for Λab is
Ba=Λabφb.
This is a system of n linear equations (index a) for the n2 unknown entries of Λ. It is underdetermined in general.

Two remarks:


5. Small 2D example (illustrative)

System: harmonic oscillator
x˙=y,y˙=−x.
Take X=x∂x+y∂y (radial scaling). Then
[X,A]=−2(x∂y+y∂x).
Expressed as [X,A]=ΛX with φ=(x,y)T, one choice is
Λ=(0−2−20).

Adapted invariants: w=y/x. Compute
w˙=−(1+w2),
so w evolves autonomously in this example. That autonomy is a special feature of the example (the last column of Λ produces a closed w-equation).


6. When is a Λ-symmetry useful?


7. Short summary