arXiv is now an independent nonprofit! Learn more
License: arXiv.org perpetual non-exclusive license
arXiv:2608.20094v1 [math.OC] 20 Aug 2026
Nonlinear Controllability and the Propagation of Local Information
From the Kalman Family to Lie Brackets, Rotation Groups, and Reachable Subgroups
Philippe Mullhaupt Affiliation: Laboratoire d’Automatique
Abstract

This article develops a self-contained Lie-theoretic route from linear controllability to nonlinear controllability on matrix Lie groups. The organizing question is whether local algebraic information can be propagated into statements about reachable sets. The linear case supplies the model: the matrix exponential, Cayley–Hamilton theorem, trajectory formula and Kalman family B,AB,,An1BB,AB,\ldots,A^{n-1}B reduce controllability to finite-dimensional linear algebra. In the nonlinear setting, Lie brackets replace matrix powers, but the associated Lie algebra may be infinite dimensional and the finite-dimensional Lie correspondence can fail. Particular emphasis is placed on rotation groups. Skew-symmetric matrices, their exponentials, commutators, one-parameter subgroups, and the Baker–Campbell–Hausdorff formula provide a concrete laboratory for seeing how infinitesimal directions propagate on SO(n)\operatorname{SO}(n). A numerical SO(3)\operatorname{SO}(3) simulation illustrates the geometry of controlled rotations, while explicit SO(4)\operatorname{SO}(4) and SO(7)\operatorname{SO}(7) commutator sequences make the propagation of a localized control direction visible. The article then develops the attainable-subgroup argument for right-invariant systems, gives a detailed proof sequence behind the Yamabe step, and works through solvable, nilpotent and ideal examples for upper-triangular matrix algebras. The finite-dimensional mechanism is finally contrasted with Sussmann’s local controllability theory for general nonlinear systems.

1 Introduction: propagating local information

A useful way to compare linear and nonlinear controllability is to ask how much information near a point can be propagated into a statement about finite-time reachability. In a linear time-invariant system, a finite family of directions — B,AB,,An1BB,AB,\ldots,A^{n-1}B — completely determines controllability. The deeper question is whether an analogous finite collection of local directions can determine reachable sets in a nonlinear system.

The natural nonlinear replacement of the Kalman family is the family of iterated Lie brackets of the vector fields defining the system. This replacement is exact in the linear case, but nonlinear vector fields introduce two essential difficulties: there is no nonlinear analogue of Cayley–Hamilton that forces the bracket family to terminate, and exact trajectories involve increasingly complicated iterated integral structures. The algebra generated by the vector fields may therefore be infinite dimensional.

Rotation matrices isolate the finite-dimensional mechanism in a setting where every object can be calculated explicitly. The Lie algebra 𝔰𝔬(n)\mathfrak{so}(n) consists of skew-symmetric matrices, the corresponding Lie group is SO(n)\operatorname{SO}(n), the exponential map is the matrix exponential, and commutators of skew-symmetric matrices are again skew symmetric. This makes the rotation group an ideal laboratory for the slogan “local information propagates.” The examples below are therefore treated as part of the main argument rather than as secondary illustrations.

2 The linear prototype

Consider the single-input LTI system

x˙=Ax+Bu,xn,u.\dot{x}=Ax+Bu,\qquad x\in\mathbb{R}^{n},\quad u\in\mathbb{R}. (1)

Let Φ(x0,t,u)\Phi(x_{0},t,u) denote the solution generated by an input u()u(\cdot). The unrestricted-time reachable set from the origin is

(0)={xn:t>0,u(),Φ(0,t,u)=x},\mathcal{R}(0)=\{x\in\mathbb{R}^{n}:\exists t>0,\exists u(\cdot),\ \Phi(0,t,u)=x\},

and the reachable set at a specified time is

(x0,t)={xn:u(),Φ(x0,t,u)=x}.\mathcal{R}(x_{0},t)=\{x\in\mathbb{R}^{n}:\exists u(\cdot),\ \Phi(x_{0},t,u)=x\}.

Small-time local controllability at the origin means that (0,t)\mathcal{R}(0,t) contains a neighborhood of the origin for every t>0t>0. Global controllability means that every target state can be reached in some finite time.

2.1 The matrix exponential and Cayley–Hamilton

The first ingredient is the exponential

eA=I+A+A22!+A33!+.e^{A}=I+A+\frac{A^{2}}{2!}+\frac{A^{3}}{3!}+\cdots.

At first sight this introduces infinitely many matrix powers. Cayley–Hamilton removes that apparent infinity: every power ApA^{p} is a linear combination of I,A,,An1I,A,\ldots,A^{n-1}. Consequently eAe^{A} itself is a linear combination of these first nn powers.

The second ingredient is the trajectory formula

Φ(x0,t,u)=eAtx0+0teA(tτ)Bu(τ)𝑑τ.\Phi(x_{0},t,u)=e^{At}x_{0}+\int_{0}^{t}e^{A(t-\tau)}Bu(\tau)\,d\tau. (2)

Together these facts give the familiar condition

rank[BABA2BAn1B]=n.\operatorname{rank}\,[B\ AB\ A^{2}B\ \cdots\ A^{n-1}B]=n. (3)

For multiple inputs, all families Bi,ABi,,An1BiB_{i},AB_{i},\ldots,A^{n-1}B_{i} are included.

2.2 The oscillator example

For

A=[0110],B=[01],A=\begin{bmatrix}0&1\\ -1&0\end{bmatrix},\qquad B=\begin{bmatrix}0\\ 1\end{bmatrix},

one has

[BAB]=[0110],[B\ AB]=\begin{bmatrix}0&1\\ 1&0\end{bmatrix},

which has rank two. Thus the position and velocity of the oscillator can be assigned arbitrarily by a suitable control over a finite interval.

3 Lie brackets as the nonlinear Kalman family

Let MM be a smooth manifold. A vector field ff assigns a tangent vector f(p)TpMf(p)\in T_{p}M at each point. In coordinates, we use the Lie bracket

[f,g]=DgfDfg.[f,g]=Dg\,f-Df\,g. (4)

It is antisymmetric and satisfies Jacobi’s identity, so the vector fields form a Lie algebra.

The linear case is recovered immediately. For f(x)=Axf(x)=Ax and the constant vector field g(x)=Bg(x)=B,

[f,g]=AB,[f,[f,g]]=A2B,[f,[f,[f,g]]]=A3B,[f,g]=-AB,\qquad[f,[f,g]]=A^{2}B,\qquad[f,[f,[f,g]]]=-A^{3}B,\ldots

Hence the Kalman family is exactly an iterated-bracket family, up to alternating signs.

For a scalar-input system

x˙=f0(x)+uf1(x),\dot{x}=f_{0}(x)+uf_{1}(x), (5)

we introduce spaces Sk(f,g)S_{k}(f,g) consisting of linear combinations of Lie monomials in ff and gg containing gg at most kk times. These spaces organize the bracket information needed in the local controllability condition.

4 Nonlinear controllability and the HLCC conditions

The following conditions are useful in the scalar-input setting, called the Hermes Local Controllability Conditions (HLCC):

  1. (HLCC 1)

    x0x_{0} is a regular equilibrium for some admissible constant control u¯\bar{u}, so f0(x0)+u¯f1(x0)=0f_{0}(x_{0})+\bar{u}f_{1}(x_{0})=0;

  2. (HLCC 2)

    dimLie(f0,f1)(x0)=dimM\dim\operatorname{Lie}(f_{0},f_{1})(x_{0})=\dim M;

  3. (HLCC 3)

    for the increasing sequence Sk(f0+u¯f1,f1)(x0)S_{k}(f_{0}+\bar{u}f_{1},f_{1})(x_{0}), equality Sk=Sk+1S_{k}=S_{k+1} is required whenever kk is odd.

Sussmann’s 1983 theorem states that these conditions imply small-time local controllability [10]. The theorem is deliberately used as motivation rather than proved: its proof is much more elaborate than the linear rank test and exposes precisely the limitations of ordinary finite-dimensional Lie theory. Sussmann’s later general theorem develops the nilpotent-approximation viewpoint further [11].

5 Rotation matrices: the central finite-dimensional example

The rotation-matrix sequence is not merely an illustration. It is the concrete finite-dimensional model through which the correspondence between a Lie algebra, its exponential map, and a Lie group becomes visible.

5.1 Skew-symmetric matrices

Consider the generic 3×33\times 3 skew-symmetric matrix

T(a,b,c)=[0aba0cbc0],T=T.T(a,b,c)=\begin{bmatrix}0&a&b\\ -a&0&c\\ -b&-c&0\end{bmatrix},\qquad T^{\top}=-T. (6)

Three numerical instances are

T1=[034305450],T2=[0613607813780],T3=[018271801427140].T_{1}=\begin{bmatrix}0&3&4\\ -3&0&5\\ -4&-5&0\end{bmatrix},\quad T_{2}=\begin{bmatrix}0&6&13\\ -6&0&78\\ -13&-78&0\end{bmatrix},\quad T_{3}=\begin{bmatrix}0&18&27\\ -18&0&14\\ -27&-14&0\end{bmatrix}.

The Lie product is the matrix commutator

[A,B]=ABBA.[A,B]=AB-BA.

For the first two matrices one obtains

[T1,T2]=[0247204247015204150],[T_{1},T_{2}]=\begin{bmatrix}0&-247&204\\ 247&0&-15\\ -204&15&0\end{bmatrix}, (7)

which is again skew symmetric. Thus the class is closed under the bracket. Antisymmetry and Jacobi’s identity follow directly from the commutator algebra.

5.2 Exponentials and orthogonal matrices

If T=TT^{\top}=-T, then

(eT)eT=eTeT=eTeT=I.(e^{T})^{\top}e^{T}=e^{T^{\top}}e^{T}=e^{-T}e^{T}=I.

Moreover deteT=etrT=1\det e^{T}=e^{\operatorname{tr}T}=1. Hence eTSO(3)e^{T}\in\operatorname{SO}(3). Direct numerical exponentiation confirms the orthogonality of rows and columns. The same construction is then repeated in dimension five, using an arbitrary 5×55\times 5 skew-symmetric matrix, to emphasize that nothing is special to dimension three.

This gives the first important propagation mechanism: the tangent directions represented by skew-symmetric matrices exponentiate into finite rotations. Studying the linear space 𝔰𝔬(n)\mathfrak{so}(n) therefore generates information about the curved manifold SO(n)\operatorname{SO}(n).

5.3 A one-parameter rotation subgroup

Consider the explicit curve

G(t)=[costsint0sintcost0001].G(t)=\begin{bmatrix}\cos t&\sin t&0\\ -\sin t&\cos t&0\\ 0&0&1\end{bmatrix}. (8)

Every G(t)G(t) is orthogonal. Differentiating at t=0t=0 gives

G˙(0)=[010100000]=:Tg.\dot{G}(0)=\begin{bmatrix}0&1&0\\ -1&0&0\\ 0&0&0\end{bmatrix}=:T_{g}.

Direct evaluation verifies that G(t)=etTgG(t)=e^{tT_{g}}. This example is the simplest manifestation of a one-parameter subgroup: an infinitesimal generator TgT_{g} determines an entire curve in the group.

5.4 The group commutator and the BCH mechanism

For two infinitesimal generators, the group commutator

eεT1eεT2eεT1eεT2e^{\varepsilon T_{1}}e^{\varepsilon T_{2}}e^{-\varepsilon T_{1}}e^{-\varepsilon T_{2}}

produces, to lowest nontrivial order, the Lie bracket direction [T1,T2][T_{1},T_{2}]. Successive truncations of the Baker–Campbell–Hausdorff expansion give

log(eAeB)=A+B+12[A,B]+112[A,[A,B]]+112[B,[B,A]]+.\log(e^{A}e^{B})=A+B+\frac{1}{2}[A,B]+\frac{1}{12}[A,[A,B]]+\frac{1}{12}[B,[B,A]]+\cdots. (9)

Numerically, adding the bracket terms progressively reduces the discrepancy between the exponential of the truncated Lie series and eAeBe^{A}e^{B}. The conceptual conclusion is explicit: local algebraic information propagates into the Lie group. The qualification is equally important: convergence limits how far a local BCH expansion can be used directly.

6 Why Lie’s third theorem matters — and why it can fail

Finite-dimensional Lie theory supplies a connected simply connected Lie group for every finite-dimensional Lie algebra. In the matrix setting the correspondence is concrete: matrix commutators live in the Lie algebra and the matrix exponential maps toward the group.

This mechanism cannot simply be transferred to the full Lie algebra of smooth vector fields. The algebra of vector fields is typically infinite dimensional. The cited example of Sergeraert [9] shows that even diffeomorphisms infinitely tangent to the identity need not embed into a one-parameter group. Thus the heuristic “vector field \leftrightarrow one-parameter subgroup” has genuine limitations outside the finite-dimensional setting.

A direct elementary warning is provided by the vector fields With

f1=(1,x12),f2=(x12,1),f_{1}=(1,x_{1}^{2})^{\top},\qquad f_{2}=(x_{1}^{2},1)^{\top},

repeated brackets produce terms of successively higher degree, including x15,x16,x17x_{1}^{5},x_{1}^{6},x_{1}^{7}, showing concretely how an infinite family of independent vector fields can arise.

7 The SO(3)\operatorname{SO}(3) controlled rotation simulation

A central example collects a nine-state nonlinear system by collecting the states into a 3×33\times 3 matrix

X=[x1x2x3x4x5x6x7x8x9].X=\begin{bmatrix}x_{1}&x_{2}&x_{3}\\ x_{4}&x_{5}&x_{6}\\ x_{7}&x_{8}&x_{9}\end{bmatrix}.

With the skew-symmetric matrices

T1=T(3,4,5),T2=T(6,13,78),T_{1}=T(3,4,5),\qquad T_{2}=T(6,13,78), (10)

the system is written

X˙=T1X+u(t)T2X,u(t)=0.2cost,X(0)=I3.\dot{X}=T_{1}X+u(t)T_{2}X,\qquad u(t)=0.2\cos t,\qquad X(0)=I_{3}. (11)

This matrix representation exposes the invariant rotation-group structure.

Since T1+u(t)T2T_{1}+u(t)T_{2} is skew symmetric for every time, the flow preserves orthogonality. Indeed,

ddt(XX)=X(T1+uT2)X+X(T1+uT2)X=0.\frac{d}{dt}(X^{\top}X)=X^{\top}(T_{1}+uT_{2})^{\top}X+X^{\top}(T_{1}+uT_{2})X=0.

With X(0)=IX(0)=I, one therefore has X(t)SO(3)X(t)\in\operatorname{SO}(3) throughout the simulation. This identity explains the geometry seen in the three Mathematica drawings: each row of X(t)X(t) remains a unit vector and hence traces a curve on the unit sphere, while the three rows remain mutually orthogonal.

Refer to caption
(a) First row of X(t)X(t).
Refer to caption
(b) Second row of X(t)X(t).
Refer to caption
(c) Third row of X(t)X(t).
Figure 1: The three trajectory drawings from the Mathematica simulation of equation 11. Each row evolves on the unit sphere because the state matrix remains in SO(3)\operatorname{SO}(3).

The corresponding Mathematica construction is included because it makes the simulation reproducible:

Listing 1: Core Mathematica simulation.
T = {{0,a,b},{-a,0,c},{-b,-c,0}};
T1 = T /. {a->3,b->4,c->5};
T2 = T /. {a->6,b->13,c->78};
XXt = {{x1[t],x2[t],x3[t]},
{x4[t],x5[t],x6[t]},
{x7[t],x8[t],x9[t]}};
xx0 = Thread[Flatten[XXt /. t->0] == Flatten[IdentityMatrix[3]]];
equDiff = Thread[Flatten[D[XXt,t]] ==
Flatten[(T1 + u T2).XXt]] /. u->0.2 Cos[t];
sols = NDSolve[Join[equDiff,xx0],Flatten[XXt],{t,0,30}][[1]];
ParametricPlot3D[Part[XXt,1]/.sols,{t,0,30},PlotRange->{{-1,1},{-1,1},{-1,1}}]

Repeating the final command for rows two and three produces figure 1.

8 Right-invariant systems on Lie groups

Consider the right-invariant control system

x˙(t)=X0(x(t))+i=1mui(t)Xi(x(t))\dot{x}(t)=X_{0}(x(t))+\sum_{i=1}^{m}u_{i}(t)X_{i}(x(t)) (12)

on a finite-dimensional Lie group GG. Let ee be the identity, (e)\mathcal{R}(e) the attainable set from ee, 𝔰\mathfrak{s} the Lie subalgebra generated by X0,,XmX_{0},\ldots,X_{m}, and SS the connected subgroup associated with 𝔰\mathfrak{s}.

Controllability reduces to two structural questions: when is (e)\mathcal{R}(e) a subgroup, and, once it is a subgroup, when is it the whole of GG?

8.1 Semigroup property

Concatenation of controls implies that (e)\mathcal{R}(e) is a semigroup. If gg is reached under one control and gg^{\prime} under a second, the concatenated control reaches the corresponding product. Right invariance is the mechanism that turns concatenation of trajectory segments into multiplication in the group.

8.2 The homogeneous case gives inverses

For the homogeneous system

x˙(t)=i=1mui(t)Xi(x(t)),\dot{x}(t)=\sum_{i=1}^{m}u_{i}(t)X_{i}(x(t)), (13)

time reversal together with sign reversal of the controls generates inverse motions. Therefore (e)\mathcal{R}(e) is a subgroup rather than merely a semigroup.

8.3 Subalgebras, ideals, solvability and nilpotency

A subalgebra 𝔥\mathfrak{h} satisfies [𝔥,𝔥]𝔥[\mathfrak{h},\mathfrak{h}]\subseteq\mathfrak{h}, while an ideal 𝔦\mathfrak{i} satisfies [𝔦,𝔤]𝔦[\mathfrak{i},\mathfrak{g}]\subseteq\mathfrak{i}. The derived series

𝔤(0)=𝔤,𝔤(k+1)=[𝔤(k),𝔤(k)]\mathfrak{g}^{(0)}=\mathfrak{g},\qquad\mathfrak{g}^{(k+1)}=[\mathfrak{g}^{(k)},\mathfrak{g}^{(k)}]

defines solvability when it reaches zero, and the lower central series

𝔤0=𝔤,𝔤k+1=[𝔤,𝔤k]\mathfrak{g}_{0}=\mathfrak{g},\qquad\mathfrak{g}_{k+1}=[\mathfrak{g},\mathfrak{g}_{k}]

defines nilpotency. The distinction is important for control because both constructions measure how rapidly new bracket directions disappear, but they do so with different nesting rules.

8.3.1 Upper-triangular matrices are solvable

Let 𝔟3\mathfrak{b}_{3} be the vector space of all real upper-triangular 3×33\times 3 matrices,

U=[abc0de00f].U=\begin{bmatrix}a&b&c\\ 0&d&e\\ 0&0&f\end{bmatrix}.

For two such matrices U1,U2U_{1},U_{2}, the commutator has the form

[U1,U2]=[0αβ00γ000],[U_{1},U_{2}]=\begin{bmatrix}0&\alpha&\beta\\ 0&0&\gamma\\ 0&0&0\end{bmatrix}, (14)

where, writing the entries of UiU_{i} as ai,bi,ci,di,ei,fia_{i},b_{i},c_{i},d_{i},e_{i},f_{i},

α=a2b1+a1b2b2d1+b1d2,\alpha=-a_{2}b_{1}+a_{1}b_{2}-b_{2}d_{1}+b_{1}d_{2},
γ=d2e1+d1e2e2f1+e1f2,\gamma=-d_{2}e_{1}+d_{1}e_{2}-e_{2}f_{1}+e_{1}f_{2},

and

β=a2c1+a1c2b2e1+b1e2c2f1+c1f2.\beta=-a_{2}c_{1}+a_{1}c_{2}-b_{2}e_{1}+b_{1}e_{2}-c_{2}f_{1}+c_{1}f_{2}.

Thus the first derived algebra lies in the strictly upper-triangular algebra 𝔫3\mathfrak{n}_{3}. If V,W𝔫3V,W\in\mathfrak{n}_{3}, then

[V,W]=[00δ000000]span{E13}.[V,W]=\begin{bmatrix}0&0&\delta\\ 0&0&0\\ 0&0&0\end{bmatrix}\in\operatorname{span}\{E_{13}\}.

Since [E13,E13]=0[E_{13},E_{13}]=0, the next derived algebra vanishes. Hence

𝔟3(0)=𝔟3,𝔟3(1)𝔫3,𝔟3(2)span{E13},𝔟3(3)=0,\mathfrak{b}_{3}^{(0)}=\mathfrak{b}_{3},\qquad\mathfrak{b}_{3}^{(1)}\subseteq\mathfrak{n}_{3},\qquad\mathfrak{b}_{3}^{(2)}\subseteq\operatorname{span}\{E_{13}\},\qquad\mathfrak{b}_{3}^{(3)}=0,

so 𝔟3\mathfrak{b}_{3} is solvable.

8.3.2 Strictly upper-triangular matrices are nilpotent

For

N=[0ab00c000]𝔫3,N=\begin{bmatrix}0&a&b\\ 0&0&c\\ 0&0&0\end{bmatrix}\in\mathfrak{n}_{3},

one computes

[N1,N2]=[00a1c2a2c1000000].[N_{1},N_{2}]=\begin{bmatrix}0&0&a_{1}c_{2}-a_{2}c_{1}\\ 0&0&0\\ 0&0&0\end{bmatrix}.

Therefore

[𝔫3,𝔫3]=span{E13},[𝔫3,E13]=0.[\mathfrak{n}_{3},\mathfrak{n}_{3}]=\operatorname{span}\{E_{13}\},\qquad[\mathfrak{n}_{3},E_{13}]=0.

The lower central series terminates after two nontrivial steps, so 𝔫3\mathfrak{n}_{3} is nilpotent of class two. The calculation also shows why nilpotency is stronger than solvability: the lower central series brackets every new layer again with the whole algebra.

8.3.3 A solvable upper-triangular algebra need not be nilpotent

Take

H=diag(1,0,0),E=E12.H=\operatorname{diag}(1,0,0),\qquad E=E_{12}.

Then H,E𝔟3H,E\in\mathfrak{b}_{3} and

[H,E]=E,adHk(E)=Efor all k1.[H,E]=E,\qquad\operatorname{ad}_{H}^{k}(E)=E\quad\text{for all }k\geq 1.

Hence the lower central series of 𝔟3\mathfrak{b}_{3} cannot terminate at zero. Thus 𝔟3\mathfrak{b}_{3} is solvable but not nilpotent. This compact computation expresses the same phenomenon as the repeated symbolic commutators of general upper-triangular matrices: diagonal directions can keep rescaling an off-diagonal direction indefinitely.

8.3.4 Two concrete ideals

Two nested subspaces make the ideal condition visible. Define

𝔦1={[abc0de000]},𝔦2={[abc00e000]}.\mathfrak{i}_{1}=\left\{\begin{bmatrix}a&b&c\\ 0&d&e\\ 0&0&0\end{bmatrix}\right\},\qquad\mathfrak{i}_{2}=\left\{\begin{bmatrix}a&b&c\\ 0&0&e\\ 0&0&0\end{bmatrix}\right\}.

For every I𝔦1I\in\mathfrak{i}_{1} and U𝔟3U\in\mathfrak{b}_{3}, the last diagonal entry of [I,U][I,U] is zero, so [𝔦1,𝔟3]𝔦1[\mathfrak{i}_{1},\mathfrak{b}_{3}]\subseteq\mathfrak{i}_{1}. Likewise the bracket of an element of 𝔦2\mathfrak{i}_{2} with an arbitrary upper-triangular matrix again has the structural form of 𝔦2\mathfrak{i}_{2}. Hence

𝔦2𝔦1𝔟3\mathfrak{i}_{2}\subset\mathfrak{i}_{1}\subset\mathfrak{b}_{3}

are ideals. These examples are useful later because reachable-set arguments naturally generate subalgebras and ideals associated with selected control directions.

9 From an attainable subgroup to the generated Lie subgroup

The central finite-dimensional statement is the following: if (e)\mathcal{R}(e) is a subgroup, then it coincides with the subgroup SS associated with the Lie algebra generated by the system fields. This is the crucial step in translating local algebraic generation into global reachability.

The subtle point is that an abstract subgroup of a Lie group need not automatically inherit the required embedded Lie-subgroup structure. Yamabe’s theorem gives the needed bridge: an arcwise connected subgroup of a Lie group is a Lie subgroup [13]. Since controlled trajectories provide arcs from the identity, (e)\mathcal{R}(e) is arcwise connected.

Let A=(e)A=\mathcal{R}(e) and let 𝔞\mathfrak{a} denote its Lie algebra. The inclusion ASA\subseteq S gives 𝔞𝔰\mathfrak{a}\subseteq\mathfrak{s}. Conversely, take constant controls with components ai{±1}a_{i}\in\{\pm 1\} and define

X(a)=X0+i=1maiXi.X(a)=X_{0}+\sum_{i=1}^{m}a_{i}X_{i}.

The positive-time curve exp(tX(a))\exp(tX(a)) lies in AA. If AA is a group, negative times lie there as well, so X(a)𝔞X(a)\in\mathfrak{a}. The resulting family of X(a)X(a) generates the original system directions, hence 𝔰𝔞\mathfrak{s}\subseteq\mathfrak{a}. Therefore

(e)=S.\mathcal{R}(e)=S. (15)

This is the classical Lie-group control mechanism developed by Jurdjevic and Sussmann [12].

10 The Yamabe mechanism: from an arcwise connected subgroup to a Lie subgroup

Yamabe’s theorem states that an arcwise connected subgroup AA of a finite-dimensional Lie group GG is a Lie subgroup [13]. For controllability this result is decisive: an attainable set may first be known only as a subgroup and as a union of controlled arcs, while the Lie-algebra argument requires a genuine Lie subgroup. The following proof sequence makes the finite-dimensional mechanism explicit.

10.1 Shrinking neighborhoods and infinitesimal directions

Let ee denote the identity of GG, and choose a nested basis of identity neighborhoods

U1U2,k1Uk={e}.U_{1}\supset U_{2}\supset\cdots,\qquad\bigcap_{k\geq 1}U_{k}=\{e\}.

Let CkC_{k} be the arcwise connected component of UkAU_{k}\cap A containing ee. Because AA is arcwise connected, every sufficiently small element of AA can be connected to ee through a curve in AA, and by shrinking the neighborhood one can keep such curves arbitrarily close to ee.

Consider sequences akCka_{k}\in C_{k}, akea_{k}\neq e, with akea_{k}\to e. In a coordinate chart around ee, the directions from ee to aka_{k} have convergent subsequences on the unit sphere. Each limiting direction determines an infinitesimal vector XX in TeGT_{e}G. Let 𝔞TeG\mathfrak{a}\subset T_{e}G be the aggregate of all infinitesimal vectors obtained in this way. The central task is to prove that 𝔞\mathfrak{a} is not merely a cone of limit directions but a Lie algebra.

For X𝔞X\in\mathfrak{a}, let

HX={exp(tX):1t1}H_{X}=\{\exp(tX):-1\leq t\leq 1\}

denote the corresponding local one-parameter subgroup. The construction of 𝔞\mathfrak{a} implies that, for every sufficiently small neighborhood VV of ee, one can find elements of AA and continuous joining curves in AA that remain arbitrarily close to HXVH_{X}V. In this sense the one-parameter subgroup is locally shadowed by arcs contained in AA.

10.2 Closure under addition

Let X,Y𝔞X,Y\in\mathfrak{a}. For large nn, choose elements of AA lying close to exp(X/n)\exp(X/n) and exp(Y/n)\exp(Y/n), with joining arcs contained in a sufficiently small neighborhood of ee. Products of these small elements remain in AA because AA is a group. The classical product limit

(exp(X/n)exp(Y/n))nexp(X+Y)\left(\exp(X/n)\exp(Y/n)\right)^{n}\longrightarrow\exp(X+Y) (16)

shows that the corresponding product arcs accumulate on the one-parameter subgroup generated by X+YX+Y. More generally, by choosing integer exponents rnr_{n} with rn/nrr_{n}/n\to r, one obtains

(exp(X/n)exp(Y/n))rnexp(r(X+Y)).\left(\exp(X/n)\exp(Y/n)\right)^{r_{n}}\longrightarrow\exp(r(X+Y)).

Hence X+YX+Y is again an infinitesimal direction generated by elements of AA, and therefore X+Y𝔞X+Y\in\mathfrak{a}. Scalar closure follows similarly by reparametrizing one-parameter subgroups. Thus 𝔞\mathfrak{a} is a vector subspace of TeGT_{e}G.

10.3 Closure under the Lie bracket

The bracket is recovered from the group commutator. For small tt,

exp(tX)exp(tY)exp(tX)exp(tY)=exp(t2[X,Y]+O(t3)).\exp(-tX)\exp(-tY)\exp(tX)\exp(tY)=\exp\!\left(t^{2}[X,Y]+O(t^{3})\right). (17)

Choose t=1/nt=1/n. Since AA is a group, the product of approximating elements and their inverses remains in AA. If integers sns_{n} are chosen with sn/n2ss_{n}/n^{2}\to s, then

(exp(X/n)exp(Y/n)exp(X/n)exp(Y/n))snexp(s[X,Y]).\left(\exp(-X/n)\exp(-Y/n)\exp(X/n)\exp(Y/n)\right)^{s_{n}}\longrightarrow\exp(s[X,Y]). (18)

The same neighborhood-control argument used for addition then shows that [X,Y]𝔞[X,Y]\in\mathfrak{a}. Therefore

X,Y𝔞X+Y𝔞,[X,Y]𝔞,X,Y\in\mathfrak{a}\quad\Longrightarrow\quad X+Y\in\mathfrak{a},\qquad[X,Y]\in\mathfrak{a},

so 𝔞\mathfrak{a} is a finite-dimensional Lie subalgebra of 𝔤=TeG\mathfrak{g}=T_{e}G.

10.4 The connected subgroup generated by 𝔞\mathfrak{a}

By finite-dimensional Lie theory there exists a unique connected immersed Lie subgroup A0GA_{0}\subset G with Lie algebra 𝔞\mathfrak{a} [1, 6, 8]. Choose a basis

X1,,XsX_{1},\ldots,X_{s}

of 𝔞\mathfrak{a} and extend it to a basis

X1,,Xs,Xs+1,,XrX_{1},\ldots,X_{s},X_{s+1},\ldots,X_{r}

of 𝔤\mathfrak{g}. Near the identity, canonical coordinates of the second kind give a local representation

g=exp(t1X1)exp(tsXs)exp(ts+1Xs+1)exp(trXr).g=\exp(t_{1}X_{1})\cdots\exp(t_{s}X_{s})\exp(t_{s+1}X_{s+1})\cdots\exp(t_{r}X_{r}). (19)

The first factor lies in A0A_{0}; the second factor is transverse to A0A_{0}.

10.5 Why AA cannot have a transverse local component

Take akCkAa_{k}\in C_{k}\subset A with akea_{k}\to e and write, using equation 19,

ak=gkhk,gkA0,a_{k}=g_{k}h_{k},\qquad g_{k}\in A_{0},

where hkh_{k} is composed only of the transverse exponentials generated by Xs+1,,XrX_{s+1},\ldots,X_{r}. Suppose infinitely many hkh_{k} were different from ee. By selecting elements fkAf_{k}\in A arbitrarily close to gkg_{k} and examining fk1akf_{k}^{-1}a_{k}, one would obtain a limiting direction belonging to 𝔞\mathfrak{a} because the elements lie in shrinking connected pieces of AA. On the other hand, the leading direction of hkh_{k} is transverse to 𝔞\mathfrak{a} by construction. This contradiction shows that sufficiently small elements of AA have no transverse component. Consequently a neighborhood of ee in AA is contained in A0A_{0}, and since both are groups this local inclusion propagates:

AA0.A\subseteq A_{0}. (20)

10.6 Recovering a neighborhood of A0A_{0} inside AA

The reverse inclusion uses the arcs that shadow the one-parameter subgroups. For each basis vector Xi𝔞X_{i}\in\mathfrak{a}, choose a continuous curve bi(t)Ab_{i}(t)\in A, 1t1-1\leq t\leq 1, which stays as close as desired to exp(tXi)\exp(tX_{i}). Define

Q={exp(t1X1)exp(tsXs):|ti|1}A0Q=\left\{\exp(t_{1}X_{1})\cdots\exp(t_{s}X_{s}):|t_{i}|\leq 1\right\}\subset A_{0}

and the map

F:QA,F(exp(t1X1)exp(tsXs))=b1(t1)bs(ts).F:Q\longrightarrow A,\qquad F\!\left(\exp(t_{1}X_{1})\cdots\exp(t_{s}X_{s})\right)=b_{1}(t_{1})\cdots b_{s}(t_{s}). (21)

When the approximating curves are chosen sufficiently close to the exact exponentials, FF is a small perturbation of the canonical-coordinate parametrization of QQ. A standard local degree/invariance-of-domain argument then implies that F(Q)F(Q) contains a neighborhood of ee in A0A_{0}. Since F(Q)AF(Q)\subset A, one obtains a neighborhood WW of ee in A0A_{0} such that WAW\subset A.

Every connected Lie group is generated by any identity neighborhood. Because A0A_{0} is connected and AA is a subgroup containing WW, it follows that

A0A.A_{0}\subseteq A. (22)

Combining equations 20 and 22 yields A=A0A=A_{0}. Thus the arcwise connected subgroup AA carries the Lie-subgroup structure associated with the infinitesimal algebra 𝔞\mathfrak{a}.

10.7 Why this theorem matters for controllability

The proof clarifies the control-theoretic role of Yamabe’s theorem. Reachability first produces paths and a semigroup or subgroup structure. The small products equation 16 and small commutators equation 18 extract addition and Lie brackets from these finite motions. Finite-dimensional Lie theory then reconstructs a connected subgroup from the resulting tangent algebra, while the local inclusion argument identifies that subgroup with the attainable subgroup itself. The conclusion is exactly the propagation principle needed in the controllability argument: infinitesimal bracket generation determines the connected subgroup that can be reached.

11 Homogeneous controllability criterion

For a homogeneous right-invariant system, (e)\mathcal{R}(e) is already a subgroup. Combining this with equation 15 gives the criterion:

A homogeneous right-invariant system is controllable on GG if and only if GG is connected and the Lie algebra generated by the controlled vector fields is all of 𝔤\mathfrak{g}.

This is the closest finite-dimensional nonlinear analogue of the Kalman rank test in this setting.

12 The homogeneous SO(4)\operatorname{SO}(4) example

Consider

X˙=AXu1+BXu2,XSO(4),\dot{X}=AXu_{1}+BXu_{2},\qquad X\in\operatorname{SO}(4), (23)

with

A=[0100101001010010],B=[0000000000010010].A=\begin{bmatrix}0&1&0&0\\ -1&0&1&0\\ 0&-1&0&1\\ 0&0&-1&0\end{bmatrix},\qquad B=\begin{bmatrix}0&0&0&0\\ 0&0&0&0\\ 0&0&0&1\\ 0&0&-1&0\end{bmatrix}.

Both are in 𝔰𝔬(4)\mathfrak{so}(4). Set

C0=B,Ck+1=[A,Ck]=adACk.C_{0}=B,\qquad C_{k+1}=[A,C_{k}]=\operatorname{ad}_{A}C_{k}.

The first seven members are

C0=[0000000000010010],C1=[0000000100000100],C_{0}=\begin{bmatrix}0&0&0&0\\ 0&0&0&0\\ 0&0&0&1\\ 0&0&-1&0\end{bmatrix},\quad C_{1}=\begin{bmatrix}0&0&0&0\\ 0&0&0&1\\ 0&0&0&0\\ 0&-1&0&0\end{bmatrix},
C2=[0001001001011010],C3=[0020000320000300],C_{2}=\begin{bmatrix}0&0&0&1\\ 0&0&1&0\\ 0&-1&0&-1\\ -1&0&1&0\end{bmatrix},\quad C_{3}=\begin{bmatrix}0&0&2&0\\ 0&0&0&-3\\ -2&0&0&0\\ 0&3&0&0\end{bmatrix},
C4=[0205205005035030],C5=[00120000131200001300],C_{4}=\begin{bmatrix}0&2&0&-5\\ -2&0&-5&0\\ 0&5&0&3\\ 5&0&-3&0\end{bmatrix},\quad C_{5}=\begin{bmatrix}0&0&-12&0\\ 0&0&0&13\\ 12&0&0&0\\ 0&-13&0&0\end{bmatrix},
C6=[012025120250025013250130].C_{6}=\begin{bmatrix}0&-12&0&25\\ 12&0&25&0\\ 0&-25&0&-13\\ -25&0&13&0\end{bmatrix}.

The algebraic pattern is clearer in graphical form. In figure 2, each cell shows a matrix entry; white cells are zero, while the signed nonzero entries are colored according to magnitude. The initial control direction B=C0B=C_{0} is supported only on the (3,4)(3,4) coordinate plane. The first commutator moves this support to the (2,4)(2,4) plane, and subsequent commutators populate further coordinate pairs until all six independent rotational planes of 𝔰𝔬(4)\mathfrak{so}(4) appear across the generated family.

Refer to caption
Figure 2: Graphical representation of Ck=adAk(B)C_{k}=\operatorname{ad}_{A}^{k}(B) for k=0,,6k=0,\ldots,6 in the SO(4)\operatorname{SO}(4) example. The expanding support makes the propagation of the localized rotational direction visible.

This is the rotation-group analogue of the linear sequence B,AB,A2B,B,AB,A^{2}B,\ldots. The relevant linear operator on the Lie algebra is now

adA:𝔰𝔬(4)𝔰𝔬(4),C[A,C].\operatorname{ad}_{A}:\mathfrak{so}(4)\to\mathfrak{so}(4),\qquad C\mapsto[A,C].

The family generated from BB by repeated application of adA\operatorname{ad}_{A} is therefore the natural finite-dimensional bracket counterpart of the Kalman family. In particular, the graphical sequence shows not merely growth of numerical coefficients but propagation of support through distinct coordinate planes.

13 The SO(7)\operatorname{SO}(7) propagation example

The same construction extends naturally to SO(7)\operatorname{SO}(7). Let A7A_{7} be the skew-symmetric nearest-neighbor chain

A7=[0100000101000001010000010100000101000001010000010],A_{7}=\begin{bmatrix}0&1&0&0&0&0&0\\ -1&0&1&0&0&0&0\\ 0&-1&0&1&0&0&0\\ 0&0&-1&0&1&0&0\\ 0&0&0&-1&0&1&0\\ 0&0&0&0&-1&0&1\\ 0&0&0&0&0&-1&0\end{bmatrix},

and localize the control generator in the final (6,7)(6,7) plane,

B7=E67E76.B_{7}=E_{67}-E_{76}.

Define D0=B7D_{0}=B_{7} and Dk+1=[A7,Dk]D_{k+1}=[A_{7},D_{k}]. The first four members are

D0=E67E76,D1=E57E75,D_{0}=E_{67}-E_{76},\qquad D_{1}=E_{57}-E_{75},
D2=(E47E74)+(E56E65)(E67E76),D_{2}=(E_{47}-E_{74})+(E_{56}-E_{65})-(E_{67}-E_{76}),
D3=(E37E73)+2(E46E64)3(E57E75).D_{3}=(E_{37}-E_{73})+2(E_{46}-E_{64})-3(E_{57}-E_{75}).

Thus the bracket depth has a direct spatial interpretation along the chain: D0D_{0} acts only at the end, D1D_{1} reaches one coordinate farther, D2D_{2} reaches the fourth coordinate, and D3D_{3} reaches the third. The full matrix values are shown graphically in figure 3.

Refer to caption
Figure 3: Graphical representation of Dk=adA7k(B7)D_{k}=\operatorname{ad}_{A_{7}}^{k}(B_{7}) for k=0,,3k=0,\ldots,3. A control generator initially confined to the (6,7)(6,7) plane propagates leftward through the nearest-neighbor chain under repeated commutation with A7A_{7}.

This example makes “propagating local information” literal. Locality refers to the support of B7B_{7} in a single rotational plane; propagation refers to the successive appearance of nonzero skew-symmetric entries in coordinate planes farther from that end. The same mechanism can be continued to higher bracket depth, and in a controllability test one asks whether the Lie algebra generated by the available directions spans all of 𝔰𝔬(7)\mathfrak{so}(7).

14 Compact groups and the nonhomogeneous case

For compact Lie groups, compactness can restore group-like reachability properties even when a drift X0X_{0} is present. In particular, in the compact Lie-group case one obtains a result of the same flavor as in the homogeneous case, whereas for a general nonhomogeneous system additional hypotheses are needed.

It is useful to distinguish three levels: finite-dimensional homogeneous systems, finite-dimensional nonhomogeneous/decomposable systems, and general infinite-dimensional vector-field algebras. The first level is completely controlled by the generated Lie algebra; the later levels require progressively more refined tools.

15 Discussion: where the analogy succeeds and where it breaks

The linear and finite-dimensional Lie-group theories share a common architecture. In the linear case, Cayley–Hamilton ensures that the directions generated by repeated action of AA form a finite family. On SO(n)\operatorname{SO}(n) and other finite-dimensional Lie groups, closure under Lie brackets and the Lie correspondence play the analogous role. BCH and group commutators explain how directions in the Lie algebra are converted into finite motions in the group.

The SO(3)\operatorname{SO}(3) simulation makes this correspondence geometric: skew-symmetric instantaneous generators preserve the orthogonal constraint, so the trajectory remains on the group and each matrix row remains on the sphere. The SO(4)\operatorname{SO}(4) and SO(7)\operatorname{SO}(7) examples make the algebraic propagation visible: a localized rotational input direction spreads through iterated commutators.

For a general nonlinear system, however, the Lie algebra of vector fields may be infinite dimensional. There is then no Cayley–Hamilton theorem to terminate the bracket generation and no unrestricted finite-dimensional Lie III correspondence to turn the algebra into a globally manageable group. This is precisely the point at which Sussmann’s formal Lie series, filtrations, nilpotent approximations, and symmetry arguments become necessary [10, 11].

16 Conclusion

The results above answer a common question at several levels: how can information known infinitesimally be propagated into a statement about finite-time reachability?

For LTI systems, the answer is the Kalman family and Cayley–Hamilton. For finite-dimensional matrix Lie groups, the answer is the Lie algebra, the exponential map, group commutators, BCH, and the subgroup generated by the controlled directions. The rotation examples are the most concrete expression of this mechanism: skew-symmetric matrices generate rotations; their commutators generate new infinitesimal rotations; periodic controlled dynamics remain on SO(3)\operatorname{SO}(3); and repeated adjoint operations on SO(4)\operatorname{SO}(4) and SO(7)\operatorname{SO}(7) visibly propagate a localized actuation direction through the algebra.

The analogy breaks when the generated vector-field algebra is infinite dimensional. The resulting failure is not incidental but structural, and it explains the technical depth of general nonlinear local controllability theorems. The finite-dimensional rotation examples therefore serve a double role: they provide exact controllability results in their own right and expose, in the cleanest possible setting, the mechanism that more general nonlinear theory attempts to recover by approximation.

Appendix A Mathematica snippets from the rotation examples

The following snippets record key calculations underlying the rotation-matrix examples.

Listing 2: Skew-symmetric matrices and commutator.
T = {{0,a,b},{-a,0,c},{-b,-c,0}};
T1 = T /. {a->3,b->4,c->5};
T2 = T /. {a->6,b->13,c->78};
T3 = T /. {a->18,b->27,c->14};
L[A_,B_] := A.B-B.A;
L[T1,T2] // MatrixForm
Listing 3: Explicit one-parameter subgroup.
Gg = {{Cos[t],Sin[t],0},{-Sin[t],Cos[t],0},{0,0,1}};
Tg = T /. {a->1,b->0,c->0};
D[Gg,t] /. t->0
MatrixExp[Tg 27456] - (Gg /. t->27456)
Listing 4: Adjoint propagation in the SO(4)\operatorname{SO}(4) example.
A = {{0,1,0,0},{-1,0,1,0},{0,-1,0,1},{0,0,-1,0}};
B = {{0,0,0,0},{0,0,0,0},{0,0,0,1},{0,0,-1,0}};
MatrixForm[#]& /@ NestList[(A.#-#.A)&,B,6]

References

  • [1] C. Chevalley, Theory of Lie Groups, Princeton University Press, 1946.
  • [2] L. S. Pontryagin, Topological Groups, Princeton University Press, 1946.
  • [3] J.-P. Serre, Lie Algebras and Lie Groups, W. A. Benjamin, 1965.
  • [4] V. S. Varadarajan, Lie Groups, Lie Algebras, and Their Representations, Springer-Verlag, 1974.
  • [5] A. A. Kirillov, Eléments de la théorie des représentations, Mir, Moscow, 1974.
  • [6] S. Helgason, Differential Geometry, Lie Groups, and Symmetric Spaces, Academic Press, 1978.
  • [7] R. Godement, Introduction à la théorie des groupes de Lie, Springer, 2004.
  • [8] A. W. Knapp, Lie Groups Beyond an Introduction, 2nd ed., Birkhäuser, 2002.
  • [9] F. Sergeraert, “Feuilletages et difféomorphismes infiniment tangents à l’identité,” Inventiones Mathematicae, vol. 39, pp. 253–275, 1977.
  • [10] H. J. Sussmann, “Lie Brackets and Local Controllability: A Sufficient Condition for Scalar-Input Systems,” SIAM Journal on Control and Optimization, vol. 21, no. 5, pp. 686–713, 1983.
  • [11] H. J. Sussmann, “A General Theorem on Local Controllability,” SIAM Journal on Control and Optimization, vol. 25, no. 1, pp. 158–194, 1987.
  • [12] V. Jurdjevic and H. J. Sussmann, “Control Systems on Lie Groups,” Journal of Differential Equations, vol. 12, pp. 313–329, 1972.
  • [13] H. Yamabe, “On an Arcwise Connected Subgroup of a Lie Group,” Osaka Mathematical Journal, vol. 2, no. 1, pp. 13–14, 1950.