Haantjes Chains in Integrable Systems

• Haantjes Chains are collections of local potentials generated via Haantjes operators within an Abelian Haantjes algebra, providing a clear criterion for Liouville-Arnold integrability.
• The vanishing Haantjes torsion in these chains guarantees the integrability of eigen-distributions and establishes commuting Hamiltonian flows and separation variables.
• An algorithmic procedure underlies the construction of Haantjes chains, facilitating the explicit determination of Darboux-Haantjes coordinates in symplectic and Jacobi-Haantjes manifolds.

Searching arXiv for the cited Haantjes-geometry papers to ground the article in the current literature. arXiv search query: "Haantjes chains symplectic Haantjes geometry integrability" Haantjes chains are collections of local potential functions $\{H_1,\dots,H_m\}$ generated by a function $H=H_1$ through the transpose action of a basis $\{K_1,\dots,K_m\}$ of an Abelian Haantjes algebra, with

$$dH_\alpha=K_\alpha^T\,dH,\qquad d\bigl(K_\alpha^T\,dH\bigr)=0.$$

In symplectic-Haantjes geometry they provide a tensorial formulation of integrals in involution, commuting Hamiltonian flows, and separation variables; in later developments the same mechanism was extended to Jacobi-Haantjes manifolds, including contact and locally conformal symplectic settings. The existence of a suitable Haantjes structure was proved to be a necessary and sufficient condition for Liouville-Arnold integrability, and under certain hypotheses it allows the determination of separation variables in an algorithmic way (Tempesta et al., 2015, Azuaje et al., 15 Jul 2025).

1. Haantjes torsion, Haantjes operators, and Haantjes algebras

The basic object is a $(1,1)$-tensor field $L:TM\to TM$. Its Nijenhuis torsion is

$$T_L(X,Y)=[LX,LY]+L^2[X,Y]-L\big([LX,Y]+[X,LY]\big),$$

and its Haantjes torsion is

$$H_L(X,Y)=L^2\,T_L(X,Y)+T_L(LX,LY)-L\big(T_L(X,LY)+T_L(LX,Y)\big).$$

An operator $L$ is called a Haantjes operator when $H_L\equiv 0$. If $H=H_1$0 is diagonal in some chart, then $H=H_1$1; for a semisimple operator with pointwise distinct eigenvalues, J. Haantjes proved that the vanishing of $H=H_1$2 is also sufficient for the integrability of each eigen-distribution (Azuaje et al., 15 Jul 2025).

A Haantjes algebra $H=H_1$3 on $H=H_1$4 is a set of Haantjes operators that is a $H=H_1$5-module and is closed under composition. In the Abelian case one also requires pairwise commutativity: $H=H_1$6 This commutative condition is central in the theory because it yields simultaneously diagonalizable families of operators in suitable coordinates and stabilizes recursive constructions of chains (Azuaje et al., 15 Jul 2025).

A recurrent source of examples is the cyclic module generated by a single Haantjes operator $H=H_1$7. If the minimal polynomial of $H=H_1$8 has degree $H=H_1$9, then

$\{K_1,\dots,K_m\}$0

is automatically an Abelian Haantjes algebra. This places polynomial recursions and multi-operator algebras within the same framework (Tempesta et al., 2014).

2. Definition of a Haantjes chain

Let $\{K_1,\dots,K_m\}$1 be a Haantjes algebra of rank $\{K_1,\dots,K_m\}$2. A function $\{K_1,\dots,K_m\}$3 generates a Haantjes chain of length $\{K_1,\dots,K_m\}$4 if there exists a basis $\{K_1,\dots,K_m\}$5 of $\{K_1,\dots,K_m\}$6 such that

$\{K_1,\dots,K_m\}$7

Equivalently, one defines

$\{K_1,\dots,K_m\}$8

and requires the $\{K_1,\dots,K_m\}$9 to be functionally independent. The functions $$dH_\alpha=K_\alpha^T\,dH,\qquad d\bigl(K_\alpha^T\,dH\bigr)=0.$$0 are the potentials of the chain (Kubů et al., 2024, Azuaje et al., 15 Jul 2025).

The geometric criterion behind this definition is the Frobenius integrability of the co-distribution generated by the transformed differentials. If $$dH_\alpha=K_\alpha^T\,dH,\qquad d\bigl(K_\alpha^T\,dH\bigr)=0.$$1 has rank $$dH_\alpha=K_\alpha^T\,dH,\qquad d\bigl(K_\alpha^T\,dH\bigr)=0.$$2 and the co-distribution

$$dH_\alpha=K_\alpha^T\,dH,\qquad d\bigl(K_\alpha^T\,dH\bigr)=0.$$3

has constant rank $$dH_\alpha=K_\alpha^T\,dH,\qquad d\bigl(K_\alpha^T\,dH\bigr)=0.$$4, and if $$dH_\alpha=K_\alpha^T\,dH,\qquad d\bigl(K_\alpha^T\,dH\bigr)=0.$$5 denotes its orthogonal annihilator distribution of rank $$dH_\alpha=K_\alpha^T\,dH,\qquad d\bigl(K_\alpha^T\,dH\bigr)=0.$$6, then $$dH_\alpha=K_\alpha^T\,dH,\qquad d\bigl(K_\alpha^T\,dH\bigr)=0.$$7 generates a Haantjes chain of length $$dH_\alpha=K_\alpha^T\,dH,\qquad d\bigl(K_\alpha^T\,dH\bigr)=0.$$8 if and only if $$dH_\alpha=K_\alpha^T\,dH,\qquad d\bigl(K_\alpha^T\,dH\bigr)=0.$$9 is Frobenius-integrable (Azuaje et al., 15 Jul 2025).

On a Poisson-Haantjes manifold, the same closure relations produce a Magri-Haantjes chain. If $(1,1)$0 is the Poisson bivector and $(1,1)$1, the Hamiltonian vector fields satisfy

$(1,1)$2

so the chain appears simultaneously at the level of 1-forms, Hamiltonians, and vector fields (Tondo, 2018).

3. Symplectic-Haantjes geometry, involution, and Liouville integrability

An $(1,1)$3 manifold is a triple $(1,1)$4 where $(1,1)$5 is symplectic of dimension $(1,1)$6, $(1,1)$7 is an Abelian Haantjes algebra of rank $(1,1)$8, and each $(1,1)$9 satisfies the compatibility condition

$L:TM\to TM$0

Within this setting, every Haantjes chain generator $L:TM\to TM$1 produces $L:TM\to TM$2 independent functions in involution: $L:TM\to TM$3 If $L:TM\to TM$4, one obtains a completely integrable Hamiltonian system (Azuaje et al., 15 Jul 2025).

For symplectic-Haantjes manifolds there is also a separation theorem. On a $L:TM\to TM$5-dimensional phase space, if $L:TM\to TM$6 is a full-length chain, then in any Darboux-Haantjes coordinate system—coordinates in which all $L:TM\to TM$7 are simultaneously diagonal and

$L:TM\to TM$8

—the Hamilton-Jacobi equation for each $L:TM\to TM$9 separates. Conversely, any family of $$T_L(X,Y)=[LX,LY]+L^2[X,Y]-L\big([LX,Y]+[X,LY]\big),$$0 functions in total separable involution in some Darboux chart arises from a Haantjes chain; in that case the operators may be reconstructed by

$$T_L(X,Y)=[LX,LY]+L^2[X,Y]-L\big([LX,Y]+[X,LY]\big),$$1

This result is presented in the magnetic-field setting as the Jacobi-Haantjes theorem (Kubů et al., 2024).

The broader structural statement is that the existence of a Haantjes structure is a necessary and sufficient condition for a Hamiltonian system to be integrable in the Liouville-Arnold sense. In that formulation, Haantjes geometry does not merely accompany integrability; it characterizes it (Tempesta et al., 2015).

4. Darboux-Haantjes coordinates, Stäckel geometry, and constructive procedures

A full Haantjes chain determines coordinates that simultaneously diagonalize the operators and separate the Hamilton-Jacobi equations. In the symplectic setting these are Darboux-Haantjes coordinates, while the associated separation data can be written in Stäckel form. Each Darboux-Haantjes chart adapted to a chain comes with a Sklyanin separation equation

$$T_L(X,Y)=[LX,LY]+L^2[X,Y]-L\big([LX,Y]+[X,LY]\big),$$2

that is, a Stäckel matrix $$T_L(X,Y)=[LX,LY]+L^2[X,Y]-L\big([LX,Y]+[X,LY]\big),$$3 and Stäckel functions $$T_L(X,Y)=[LX,LY]+L^2[X,Y]-L\big([LX,Y]+[X,LY]\big),$$4. In magnetic examples, distinct Darboux-Haantjes realizations of the same model give rise to distinct Stäckel matrices; in the inverse problem, these matrices can be upgraded to arbitrary functions $$T_L(X,Y)=[LX,LY]+L^2[X,Y]-L\big([LX,Y]+[X,LY]\big),$$5, yielding new families of integrable Hamiltonians with magnetic field on suitably curved metrics while preserving the same Haantjes web and Stäckel structure (Kubů et al., 2024).

The constructive side of the theory is explicit. A standard procedure consists of the following steps (Tempesta et al., 2014):

1. Choose the symplectic form $$T_L(X,Y)=[LX,LY]+L^2[X,Y]-L\big([LX,Y]+[X,LY]\big),$$6 and a first integral $$T_L(X,Y)=[LX,LY]+L^2[X,Y]-L\big([LX,Y]+[X,LY]\big),$$7.
2. Make an ansatz for a tensor $$T_L(X,Y)=[LX,LY]+L^2[X,Y]-L\big([LX,Y]+[X,LY]\big),$$8 with undetermined components.
3. Impose algebraic compatibility $$T_L(X,Y)=[LX,LY]+L^2[X,Y]-L\big([LX,Y]+[X,LY]\big),$$9, the chain condition $$H_L(X,Y)=L^2\,T_L(X,Y)+T_L(LX,LY)-L\big(T_L(X,LY)+T_L(LX,Y)\big).$$0, and the vanishing of the Haantjes torsion.
4. Solve the resulting linear/algebraic constraints and first-order PDEs.
5. If $$H_L(X,Y)=L^2\,T_L(X,Y)+T_L(LX,LY)-L\big(T_L(X,LY)+T_L(LX,Y)\big).$$1 is cyclic of maximal degree, construct further operators $$H_L(X,Y)=L^2\,T_L(X,Y)+T_L(LX,LY)-L\big(T_L(X,LY)+T_L(LX,Y)\big).$$2 as polynomials in $$H_L(X,Y)=L^2\,T_L(X,Y)+T_L(LX,LY)-L\big(T_L(X,LY)+T_L(LX,Y)\big).$$3.
6. Verify pairwise commutativity and recover the 1-forms $$H_L(X,Y)=L^2\,T_L(X,Y)+T_L(LX,LY)-L\big(T_L(X,LY)+T_L(LX,Y)\big).$$4, which form the chain.

This explains the recurrent description of Haantjes geometry as algorithmic: once the compatible operators are found, the corresponding potentials, commuting flows, and separation coordinates follow by tensorial recursion (Tempesta et al., 2015).

5. Representative integrable models

Several concrete systems have been exhibited as Haantjes chains in symplectic-Haantjes or Poisson-Haantjes geometry (Kubů et al., 2024, Tondo, 2018, Tempesta et al., 2014).

System	Geometric setting	Chain data
Cylindrical separation	$$H_L(X,Y)=L^2\,T_L(X,Y)+T_L(LX,LY)-L\big(T_L(X,LY)+T_L(LX,Y)\big).$$5 in $$H_L(X,Y)=L^2\,T_L(X,Y)+T_L(LX,LY)-L\big(T_L(X,LY)+T_L(LX,Y)\big).$$6	Abelian algebra $$H_L(X,Y)=L^2\,T_L(X,Y)+T_L(LX,LY)-L\big(T_L(X,LY)+T_L(LX,Y)\big).$$7, $$H_L(X,Y)=L^2\,T_L(X,Y)+T_L(LX,LY)-L\big(T_L(X,LY)+T_L(LX,Y)\big).$$8
Constant magnetic field	$$H_L(X,Y)=L^2\,T_L(X,Y)+T_L(LX,LY)-L\big(T_L(X,LY)+T_L(LX,Y)\big).$$9 with $L$0	$L$1, $L$2
Helical undulator	Rotating magnetic field	Three integrals, operators $L$3, full or partial separation
Lagrange top	Symplectic leaf $L$4	Two-generator algebra $L$5, Darboux-Haantjes coordinates $L$6
Post-Winternitz system	$L$7 after canonical change	$L$8

In the cylindrical magnetic model, the Abelian Haantjes algebra is

$L$9

with

$H_L\equiv 0$0

$H_L\equiv 0$1

For suitable magnetic vector potential $H_L\equiv 0$2 and scalar potential $H_L\equiv 0$3, the natural Hamiltonian

$H_L\equiv 0$4

generates a chain of length three, and the Darboux-Haantjes chart is precisely $H_L\equiv 0$5. In the constant-field model, the same formalism recovers separation in Cartesian, shifted Cartesian, and cylindrical coordinates through three different Darboux-Haantjes charts. In the helical undulator, the field

$H_L\equiv 0$6

admits three integrals closing a Lie algebra and commuting Haantjes operators that yield either full separation or partial separation in an adapted chart (Kubů et al., 2024).

The Lagrange top furnishes an explicit low-dimensional example. On the four-dimensional symplectic leaf $H_L\equiv 0$7, one has a Nijenhuis operator

$H_L\equiv 0$8

with minimal polynomial

$H_L\equiv 0$9

and a two-generator Abelian Haantjes algebra

$H=H_1$00

The Darboux-Haantjes coordinates are

$H=H_1$01

in which $H=H_1$02 is diagonal and the Hamilton-Jacobi equation separates (Tondo, 2018).

For the Post-Winternitz system, two nontrivial Haantjes operators satisfy

$H=H_1$03

They commute, have vanishing Haantjes torsion, and generate three independent integrals in involution; the construction is presented as a worked example of the general Haantjes-chain procedure (Tempesta et al., 2014).

6. Jacobi-Haantjes extensions, dissipation, and generalized Lenard-Magri chains

The original symplectic formulation has been generalized to Jacobi-Haantjes manifolds, where the algebra acts on Jacobi data and includes contact and locally conformal symplectic reductions. In the contact case, a contact-Haantjes manifold $H=H_1$04 satisfies

$H=H_1$05

where $H=H_1$06 is the Reeb field and $H=H_1$07 the contact Hamiltonian field of $H=H_1$08. In the locally conformal symplectic case, an LCS-Haantjes manifold $H=H_1$09 satisfies

$H=H_1$10

In both cases a Hamiltonian $H=H_1$11 generates a dissipative analogue of a Haantjes chain, with

$H=H_1$12

and generalized involution relations with respect to the Jacobi bracket (Azuaje et al., 15 Jul 2025).

For contact chains, the involution relations take the form

$H=H_1$13

The paper states that each $H=H_1$14 is constant along the reduced invariant foliation $H=H_1$15, which is described as partial integrability. A three-dimensional example on coordinates $H=H_1$16 with

$H=H_1$17

admits two Haantjes operators

$H=H_1$18

yielding the 2-chain $H=H_1$19 (Azuaje et al., 15 Jul 2025).

A related extension appears in three-dimensional Poisson quasi-Nijenhuis geometry. For an oriented $H=H_1$20-manifold with nowhere-vanishing Poisson tensor $H=H_1$21, the operator

$H=H_1$22

has vanishing Haantjes torsion, so every such involutive Poisson quasi-Nijenhuis manifold is a Haantjes manifold. With

$H=H_1$23

the powers $H=H_1$24 satisfy the generalized Lenard-Magri conditions, and each 1-form

$H=H_1$25

is closed. Locally there exist functions $H=H_1$26 such that

$H=H_1$27

and the Hamiltonians $H=H_1$28 pairwise Poisson-commute (Vizarreta et al., 23 Feb 2025).

These generalizations clarify two recurrent points. First, vanishing Haantjes torsion is weaker than vanishing Nijenhuis torsion, but it still guarantees integrability of generalized eigen-distributions. Second, Haantjes chains are not confined to conservative symplectic systems: the same tensorial mechanism governs complete integrability on $H=H_1$29 manifolds, partial integrability on contact or LCS leaves, and generalized Lenard-Magri recursions in Poisson quasi-Nijenhuis geometry (Azuaje et al., 15 Jul 2025, Vizarreta et al., 23 Feb 2025).