Group-invariant solutions and symmetry reduction
Anchor (Master): Olver §3.1-§3.4; Ovsiannikov *Group Analysis of Differential Equations* (1982); Bluman-Cole *Similarity Methods for Differential Equations* (1974); Barenblatt *Similarity, Self-Similarity, and Intermediate Asymptotics* (1979)
Intuition Beginner
A symmetry of an equation is a motion that turns solutions into solutions. Most solutions get pushed to a different solution by the motion. But a few special solutions sit exactly on the motion, so that pushing them leaves them where they were. These are the symmetry-respecting solutions, and they are the easiest ones to find. They are pinned in place by the very symmetry that scrambles everything else.
Why does this help? A solution that respects a symmetry repeats itself in a predictable way, so it carries less independent information. If a solution looks the same after you scale space and time together, then it really only depends on one combined quantity rather than on space and time separately. That combined quantity is the new variable. Trading two variables for one turns a hard equation in two variables into a simpler equation in one.
This is the engine of similarity methods. You spot a symmetry, you find the combination of variables it leaves alone, and you write your unknown as a function of that combination only. The original equation then collapses into a smaller one. For an ordinary differential equation the same trick lowers the order by one each time you use a symmetry, because a symmetry lets you measure progress along its own flow and forget one layer of the problem.
Visual Beginner
Picture the plane of space and time, with the lines that a scaling symmetry moves along drawn as curves spreading out from the origin like spokes. A general solution cuts across these spokes. A symmetry-respecting solution is constant along each spoke: it takes the same value everywhere on one spoke and changes only as you step from spoke to spoke.
Because the solution is constant along a spoke, the only thing that matters is which spoke you stand on. Labelling the spokes by a single number gives the new variable, the similarity variable. The picture shows that two-dimensional data has collapsed to one-dimensional data: the whole solution is recovered from its values across the fan. Reducing the equation is then the act of rewriting it in terms of that one label.
Worked example Beginner
Take the heat equation, which says the rate of change in time of the temperature equals the bending of the temperature in space. Heat spreading from a point looks the same whether you zoom out in space by some factor and in time by the square of that factor. That paired stretch is a scaling symmetry.
So the temperature should depend only on the combination that survives the stretch. If you stretch space by a number and time by that number squared, the ratio of position to the square root of time does not change. Call that ratio the similarity variable: position divided by the square root of time. The claim is that the temperature is a function of this single ratio alone, times a fixed power of time out front to balance units.
Check what this buys you. Instead of a temperature that depends on both position and time, you have a temperature built from one combined quantity. Plugging this form into the heat equation, the two-variable equation becomes a single ordinary differential equation for a function of the ratio. Solving that one-variable equation and undoing the substitution gives the famous bell-shaped spreading profile, written with the error function.
What this tells us: a scaling symmetry let us guess that the answer depends on one combined variable, and that guess turned a partial differential equation in two variables into an ordinary differential equation in one. The reduction came directly from respecting the symmetry.
Check your understanding Beginner
Formal definition Intermediate+
Let be a fibred manifold with coordinates , , , let be a system of PDEs as in 05.05.06, and let be a (local) symmetry group of with Lie algebra of generators . The development follows Olver [Olver §3.1] and Ovsiannikov [Ovsiannikov Ch. 5].
Definition (group-invariant solution). Let be a connected -dimensional subgroup acting on with generators spanning a subalgebra . A solution is -invariant (a group-invariant or similarity solution) when its graph is an -invariant subset: for all . Equivalently, each generator is tangent to , which in coordinates is the invariant-surface condition $$ \xi^i(x, f),\frac{\partial f^\alpha}{\partial x^i} = \phi^\alpha(x, f), \qquad 1 \le \alpha \le q, $$ a first-order quasilinear PDE system that the invariant solution must satisfy in addition to .
Definition (invariants and the reduced variables). A function is an invariant of when for every , i.e. . For a single generator the invariants are the first integrals of the characteristic system $$ \frac{dx^1}{\xi^1} = \cdots = \frac{dx^p}{\xi^p} = \frac{du^1}{\phi^1} = \cdots = \frac{du^q}{\phi^q}. $$ Where acts on with generic orbit dimension and the projected action on has orbit dimension , there are functionally independent invariants among the base variables (the similarity variables ) and further invariants involving the fibre. The invariant-solution ansatz expresses the as functions of the , which on solving for recovers the invariant-surface condition.
Definition (canonical coordinates). For a single nonvanishing generator on a surface, canonical coordinates are functions with and , so that . The coordinate is an invariant; is a coordinate along the orbit. They exist locally near any point where by the flow-box theorem of 02.12.02.
A non-example sharpens the meaning. A solution that merely intersects an orbit of is not -invariant; invariance demands that the whole graph be a union of orbits, equivalently that every be tangent to at every point. A solution can be invariant under a one-parameter yet fail to be invariant under a larger ; the reduction is always relative to the chosen subgroup.
Counterexamples to common slips
- The invariant-surface condition is necessary but not sufficient: an invariant solution must satisfy both the invariant-surface condition and the original system . Solving only the surface condition produces invariant graphs, not invariant solutions.
- The number of similarity variables is , governed by the orbit dimension of the action projected to , not by . A vertical generator (acting only on ) reduces no independent variables; it constrains the fibre dependence instead.
- Reducing by a subgroup that does not act transversally to the equation can yield a singular or empty reduction; the rank/transversality condition on the invariant manifold must hold for the reduced equation to be a genuine system of the expected order.
Key theorem with proof Intermediate+
Theorem (reduction of order by a one-parameter symmetry). Let the scalar -th order ODE $$ u^{(n)} = F\bigl(x, u, u', \dots, u^{(n-1)}\bigr) $$ admit the one-parameter symmetry group generated by , . Then in canonical coordinates with the equation takes a form independent of , and writing it becomes an -th order ODE for . Thus a point symmetry reduces the order of a scalar ODE by one.
The statement follows Olver [Olver §2.5] and Bluman-Kumei [Bluman-Kumei Ch. 4].
Proof. Choose canonical coordinates near a point where , so that and ; these exist by the flow-box theorem of 02.12.02, which rectifies a nonvanishing field to . In these coordinates the generator is , and its prolongation to the jet coordinates is as well, since the prolongation coefficients of all vanish: with , in the frame the characteristic is and for by the prolongation formula of 05.05.06.
Write the equation in the frame, taking as the independent variable and as the dependent variable. It becomes for some function . The infinitesimal symmetry criterion of 05.05.06 requires on the equation. Since , this reads
$$
\partial_s\bigl(s^{(n)} - G\bigr) = -,\partial_s G = 0,
$$
so does not depend on : . The equation in the canonical frame contains only through its derivatives.
Introduce the new dependent variable . Then , and inductively for , so the -th order equation becomes $$ w^{(n-1)} = G\bigl(r, w, w', \dots, w^{(n-2)}\bigr), $$ an ODE of order for . Solving it for and integrating once recovers up to a constant, the constant being the orbit parameter along ; undoing the canonical change of coordinates returns the original solution. The order has dropped from to , with one quadrature restoring the lost order.
Bridge. This is the elementary core of Lie's program: each symmetry buys one reduction of order, and a solvable -dimensional symmetry algebra integrates a scalar ODE completely by successive reductions. It builds toward the PDE case below, where canonical coordinates are replaced by joint invariants and reduction of order becomes reduction of the number of independent variables. It appears again in the optimal-system classification, where one selects which one-parameter subgroup to reduce by, and in 02.12.02, whose flow-box rectification is the geometric content of canonical coordinates. The construction reuses the prolongation formula of 05.05.06: straightening to kills every prolongation coefficient, so the criterion degenerates to the single statement that is absent, which is exactly the structural fact that one variable has become ignorable.
Exercises Intermediate+
Lean formalization Intermediate+
Mathlib has flows and Lie groups but no symmetry-reduction apparatus for differential equations; the block fixes intended signatures, with the gap detailed in Mathlib gap analysis.
import Mathlib.Geometry.Manifold.ContMDiff.Basic
import Mathlib.Analysis.ODE.Gronwall
variable {p q : ℕ}
/-- A generator on E: horizontal part ξ and vertical part φ. -/
structure Generator (p q : ℕ) where
xi : (Fin p → ℝ) → (Fin q → ℝ) → (Fin p → ℝ)
phi : (Fin p → ℝ) → (Fin q → ℝ) → (Fin q → ℝ)
/-- An invariant of v: a function annihilated by v = ξ^i ∂_{x^i} + φ^α ∂_{u^α}. -/
def IsInvariant
(v : Generator p q)
(I : (Fin p → ℝ) → (Fin q → ℝ) → ℝ) : Prop :=
∀ x u, (∑ i, v.xi x u i * (deriv (fun t => I (Function.update x i t) u) (x i)))
+ (∑ α, v.phi x u α * (deriv (fun t => I x (Function.update u α t)) (u α))) = 0
/-- Invariant-surface condition: ξ^i ∂_i f^α = φ^α along the graph of f. -/
def InvariantSurface
(v : Generator p q) (f : (Fin p → ℝ) → (Fin q → ℝ)) : Prop :=
∀ x α, (∑ i, v.xi x (f x) i * (deriv (fun t => f (Function.update x i t) α) (x i)))
= v.phi x (f x) α
/-- Canonical coordinates (r,s) straighten v to ∂_s: v(r)=0, v(s)=1.
Reduction of order then follows; statement against the absent jet apparatus. -/
theorem reduction_of_order
(v : Generator 1 1) : True := by
trivial -- placeholder: in (r,s) the n-th order ODE loses s and drops to order n-1The by-block placeholder and the deriv-based stand-ins rest on the jet/prolongation gap of 05.05.05 and 05.05.06 plus the reduced-equation construction. Each sits on the existing flow and Lie-group libraries and is a Mathlib-contribution-sized target named in Mathlib gap analysis.
Advanced results Master
Reduction turns an invariance into a smaller equation, and the precise bookkeeping of how much smaller is the content of the reduction theorem. Two structural results frame the method, and three worked reductions exhibit it producing closed-form solutions.
The reduction theorem. Let be an -dimensional symmetry group of acting on so that its projection to has generic orbit dimension , and assume the action is transversal to the equation in the sense that the invariant jet manifold has maximal rank. Then the -invariant solutions are governed by a reduced system in independent variables, the similarity variables , obtained by restricting to the invariant jets. The map from invariant solutions of to solutions of the reduced system is a bijection (Olver [Olver §3.2], Ovsiannikov [Ovsiannikov Ch. 5]). For , this is the passage from a PDE in to an ODE in one similarity variable; the heat-equation reduction below is the canonical instance.
The heat-equation source solution. For the scaling subalgebra is spanned by and the dilation . Reducing by the one-parameter group generated by gives base invariant and fibre invariant , so the ansatz is . Substitution yields the reduced ODE $$ \psi'' + \tfrac12\eta,\psi' + \beta,\psi = 0. $$ For the equation integrates to , and undoing the ansatz gives the fundamental (source) solution $$ u(x,t) = \frac{C}{\sqrt t},\exp!\Bigl(-\frac{x^2}{4t}\Bigr), $$ the heat kernel up to normalisation. For the reduced ODE is with solution the error-function profile , the spreading step. Both are pure outputs of the scaling reduction; the choice of selects the self-similar exponent, exactly the self-similar solution of the first kind in Barenblatt's terminology [Barenblatt Ch. 1].
The porous-medium / Burgers traveling wave. The space-time translation reduces any translation-invariant evolution PDE to a traveling-wave ODE in . For Burgers this is , integrating to the shock profile; for the porous-medium equation the same translation gives a degenerate ODE whose compactly supported traveling fronts are the Barenblatt-type waves with finite propagation speed. The traveling-wave reduction is the one-dimensional shadow of the two-dimensional translation subgroup.
Optimal systems. Distinct invariant solutions correspond to subalgebras of taken up to conjugacy under the adjoint action . An optimal system of -dimensional subalgebras is a list of representatives, one per conjugacy class, computed by exponentiating the adjoint action and using it to simplify an arbitrary generator to canonical form. For the heat equation, the one-dimensional optimal system reduces the six-dimensional point algebra to a short list — translation , scaling , Galilean , and a few combinations — and every reduction not equivalent to these is redundant. Ovsiannikov introduced the systematic computation [Ovsiannikov Ch. 6].
Synthesis. Symmetry reduction is the conversion of an invariance into a quotient equation, and the orbit structure of the chosen subgroup is the conversion ratio. The method proceeds along four linked stages, each carried by the generators produced in 05.05.06. First, a subgroup is selected and its joint invariants are computed as first integrals of the characteristic system , replacing the original coordinates by similarity variables. Second, the invariant-surface condition is imposed, expressing the dependent variables through those invariants, which is the invariant-solution ansatz. Third, substitution into produces the reduced system in independent variables, an ODE when the count reaches one, integrated in closed form for the heat, Burgers, and porous-medium examples. Fourth, the adjoint action of on collapses the redundant reductions into an optimal system of subalgebras, so that finitely many representatives exhaust all group-inequivalent invariant solutions. The reduction of order of an ODE and the reduction of independent variables of a PDE are the same operation seen at and , and the canonical coordinates that straighten a generator to are the local model that both specialise.
Full proof set Master
Proposition (similarity reduction of a scaling-invariant evolution equation). Let an evolution equation in one space dimension admit the scaling group with generator , . Then every -invariant solution has the self-similar form with , and satisfies an ODE obtained by substitution; conversely each solution of that ODE yields a -invariant solution.
Proof. The invariants of are the first integrals of the characteristic system . From , , giving the base invariant . From , ... more directly , so , giving the fibre invariant . The invariant-surface condition -tangency forces to be a function of alone: , i.e. $$ u(x, t) = t^{b/a},\psi(\eta), \qquad \eta = x,t^{-1/a}. $$ Differentiating, and each spatial derivative , since . Substituting into and dividing by the common power — which is consistent precisely because is a symmetry, so every term carries the same scaling weight — eliminates and leaves an ODE purely in . Conversely, given a solution of that ODE, the function satisfies by reversing the substitution, and it is -invariant because as is an invariant. Hence the correspondence is a bijection between -invariant solutions and solutions of the reduced ODE.
Proposition (commuting reductions and successive order drop). Let a scalar ODE of order admit a two-dimensional abelian symmetry algebra with and pointwise independent. Then the equation can be reduced to order by two successive canonical-coordinate reductions, and the reductions can be performed in either order.
Proof. By the reduction-of-order theorem, choose canonical coordinates for , so and the equation, written for , loses explicit dependence on ; setting drops the order to in the variable . It remains to show descends to a symmetry of the reduced equation. Since , the field commutes with , so in the frame has coefficients independent of (commuting with forbids -dependence). Its action on the -free reduced equation is well defined and is again a symmetry, because prolongation commutes with the projection that forgets : the reduced equation is the restriction of the original to the -independent jets, and preserves that restriction as it preserved the full equation. Hence the order- equation for admits the projected one-parameter symmetry , and a second canonical-coordinate reduction drops the order to . Commutativity makes the two flows independent, so reducing first by and then by yields the same order- equation up to a change of variables, by symmetry of the argument.
Connections Master
The generators and characteristics computed in 05.05.06 are the raw material of this unit: the invariant-surface condition is the vanishing of the characteristic along the solution graph, so an invariant solution is one on which every generator's characteristic vanishes, and the determining equations there select which generators are available to reduce by here.
The jet bundle and total derivative of 05.05.05 persist as the substrate: the reduced equation is the restriction of to the invariant jet submanifold, computed by the same total-derivative prolongation, and the similarity variables are coordinates on the quotient of the jet bundle by the prolonged group action.
Canonical coordinates straightening a generator to are the flow-box rectification of 02.12.02: reduction of order is the statement that an ignorable coordinate along a symmetry flow can be integrated out, the differential-equation analogue of passing to a quotient by a one-parameter group, and the orbit parameter restored by quadrature is the flow time of that group.
Historical & philosophical context Master
Sophus Lie developed the reduction of differential equations by their symmetry groups as the centerpiece of his theory of continuous groups, explicitly modeling it on Galois theory: where Galois resolved the solvability of a polynomial through the structure of its permutation group, Lie sought to resolve the integrability of a differential equation through the structure of its continuous symmetry group. His lectures, edited by Georg Scheffers as Vorlesungen über Differentialgleichungen mit bekannten infinitesimalen Transformationen (1891) [Lie 1895], laid out reduction of order by canonical coordinates and the integration of an ODE by a solvable symmetry algebra. The central theorem — that a one-parameter symmetry lowers the order of a scalar ODE by exactly one, and a solvable -dimensional algebra integrates an -th order equation completely — is Lie's.
The method was carried into mathematical physics through dimensional analysis and similarity solutions before the group-theoretic origin was widely recognised. Garrett Birkhoff's Hydrodynamics (1950) [Birkhoff 1950] framed similarity solutions of fluid equations explicitly as invariance under groups, connecting Lie's reduction to the physicists' dimensional method. George Barenblatt's Similarity, Self-Similarity, and Intermediate Asymptotics (1979) [Barenblatt Ch. 1] distinguished self-similar solutions of the first kind, whose exponents are fixed by dimensional analysis, from those of the second kind, whose exponents are determined by a nonlinear eigenvalue problem and lie outside naive scaling. Lev Ovsiannikov's Group Analysis of Differential Equations (1982) [Ovsiannikov Ch. 5] systematised invariant and partially invariant solutions and the optimal-system classification for continuum mechanics, and George Bluman and Julian Cole's Similarity Methods for Differential Equations (1974) and their 1969 nonclassical method [Bluman-Cole Ch. 2] extended the construction to ansätze invariant under symmetries of the invariant-surface condition rather than of the equation itself. Peter Olver's Applications of Lie Groups to Differential Equations (1986; second edition 1993) [Olver §3.1] gave the modern geometric account in terms of group invariants, the reduced equation, and optimal systems used in this unit.
Bibliography Master
@book{LieScheffers1891,
author = {Lie, Sophus and Scheffers, Georg},
title = {Vorlesungen {\"u}ber Differentialgleichungen mit bekannten infinitesimalen Transformationen},
publisher = {B. G. Teubner},
address = {Leipzig},
year = {1891}
}
@book{OlverLieGroups,
author = {Olver, Peter J.},
title = {Applications of Lie Groups to Differential Equations},
series = {Graduate Texts in Mathematics},
volume = {107},
publisher = {Springer},
year = {1993},
edition = {2nd}
}
@book{Ovsiannikov1982,
author = {Ovsiannikov, Lev V.},
title = {Group Analysis of Differential Equations},
publisher = {Academic Press},
year = {1982},
note = {translated by W. F. Ames}
}
@book{BlumanCole1974,
author = {Bluman, George W. and Cole, Julian D.},
title = {Similarity Methods for Differential Equations},
series = {Applied Mathematical Sciences},
volume = {13},
publisher = {Springer},
year = {1974}
}
@article{BlumanCole1969,
author = {Bluman, George W. and Cole, Julian D.},
title = {The general similarity solution of the heat equation},
journal = {Journal of Mathematics and Mechanics},
volume = {18},
year = {1969},
pages = {1025--1042}
}
@book{Barenblatt1979,
author = {Barenblatt, Grigory I.},
title = {Similarity, Self-Similarity, and Intermediate Asymptotics},
publisher = {Consultants Bureau},
address = {New York},
year = {1979}
}
@book{Birkhoff1950,
author = {Birkhoff, Garrett},
title = {Hydrodynamics: A Study in Logic, Fact, and Similitude},
publisher = {Princeton University Press},
year = {1950}
}
@book{BlumanKumei1989,
author = {Bluman, George W. and Kumei, Sukeyuki},
title = {Symmetries and Differential Equations},
series = {Applied Mathematical Sciences},
volume = {81},
publisher = {Springer},
year = {1989}
}