---
title: Hamiltonian Cubical Tensors
url: https://www.emergentmind.com/topics/hamiltonian-cubical-tensors
type: topic
---

# Hamiltonian Cubical Tensors

Hamiltonian cubical tensors are a multilinear generalization of Hamiltonian matrices, capturing the fundamental algebraic and geometric properties of Hamiltonian systems in the context of higher-order tensors. These structures play a central role in the emerging theory of tensor-based polynomial Hamiltonian systems and in the extension of classical symplectic geometry to the algebra of third-order tensors via the T-product formalism. They enable direct analysis and synthesis of Hamiltonian properties, invariants, and stability criteria for polynomial dynamical systems with fundamentally multilinear interactions.

## 1. Definitions and Algebraic Characterizations

Let $n$ be even and $J \in \R^{n \times n}$ the canonical symplectic matrix,
\[
J = \begin{pmatrix}
0 & I_{n/2} \\
- I_{n/2} & 0
\end{pmatrix}
\]
A $k$th-order $n$-dimensional cubical tensor is a multi-array $A \in \R^{[k, n]} = \R^{n \times n \times \cdots \times n}$ (with $k$ modes). Hamiltonian cubical tensors are defined as follows [2503.21487]:

A tensor $A \in \R^{[k, n]}$ is called Hamiltonian cubical (with respect to $J$) if for every nontrivial permutation $\sigma \in S_k$,
\[
\left( J^\top A \right)^{\top_\sigma} + J A = 0 \in \R^{[k, n]}
\]
where $A^{\top_\sigma}$ denotes the mode permutation of $A$ according to $\sigma$. For $k=2$ this recovers the classical Hamiltonian matrix condition $A^\top J + J A = 0$.

The following are equivalent for $A \in \R^{[k, n]}$:
1. $A$ is Hamiltonian cubical.
2. There exists a supersymmetric tensor $R \in \R^{[k, n]}$ (i.e., $R^{\top_\sigma} = R$ for all $\sigma \in S_k$) such that $A = J R$.
3. The tensor $J A$ is supersymmetric.

This equivalence generalizes the classical result that every Hamiltonian matrix can be written as $A = J R$ with $R$ symmetric [2503.21487].

## 2. Hamiltonian Cubical Tensors in the T-product Algebra

In the T-product framework, the theory is developed for third-order (“cubical”) tensors $\mathcal{A} \in \mathbb{C}^{m \times n \times p}$, represented via their block-circulant matricization and manipulated using the circular T-product. The T-product is defined so that $\text{bcirc}(\mathcal{A} * \mathcal{B}) = \text{bcirc}(\mathcal{A}) \cdot \text{bcirc}(\mathcal{B})$, or equivalently via a diagonalization by discrete Fourier transform (DFT) along the third mode [2605.20829].

Within this formalism, the symplectic unit tensor $\mathcal{J}$ is defined so that each of its DFT frontal slices is $J$. A tensor $\mathcal{H} \in \mathbb{C}^{2n \times 2n \times p}$ is T-Hamiltonian if
\[
(\mathcal{J} * \mathcal{H})^H = \mathcal{J} * \mathcal{H}
\]
where $^H$ denotes the T-conjugate transpose.

In the Fourier domain, this equates to each slice $\hat{\mathcal{H}}^{(i)}$ satisfying the classical Hamiltonian matrix condition:
\[
(\mathbf{J} \hat{\mathcal{H}}^{(i)})^{H} = \mathbf{J} \hat{\mathcal{H}}^{(i)}
\]
or, equivalently,
\[
\hat{\mathcal{H}}^{(i)H} J + J \hat{\mathcal{H}}^{(i)} = 0
\]
Therefore, T-Hamiltonian tensors in this setting are precisely those whose DFT slices are Hamiltonian matrices [2605.20829].

## 3. Hamiltonian Polynomial Systems and Tensor-based Structure

A polynomial vector field on $\R^n$ admits the representation:
\[
\dot{x} = A_k x^{k-1} + A_{k-1} x^{k-2} + \cdots + A_2 x
\]
where $A_j \in \R^{[j,n]}$ are tensors. A polynomial Hamiltonian is similarly expressed as
\[
H(x) = B_k x^k + B_{k-1} x^{k-1} + \cdots + B_2 x^2
\]
with $B_j$ supersymmetric cubical tensors. The multivariate gradient structure of $H$ yields
\[
\nabla H = \sum_{j=2}^k j B_j x^{j-1}
\]
and the induced Hamiltonian flow is $\dot x = J \nabla H(x)$. The polynomial system is Hamiltonian with Hamiltonian $H$ if and only if all system tensors $A_j$ are Hamiltonian cubical. In this case,
\[
A_j = j J B_j,\qquad B_j = \frac{1}{j} J^\top A_j, \quad j=2,\ldots,k
\]
This criterion extends the matrix-based test for the Hamiltonian property and allows algorithmic recovery of the Hamiltonian structure from the polynomial flow [2503.21487].

## 4. Spectral Theory, T-eigenvalues, and Normal Forms

The spectral properties of Hamiltonian cubical tensors generalize those of classical Hamiltonian matrices. In the T-product algebra, the T-Jordan canonical form shows every third-order tensor is T-similar to a tensor whose DFT slices are in ordinary Jordan form, with diagonal entries termed T-eigenvalues [2605.20829].

For T-Hamiltonian tensors, the spectrum is symmetric under $\lambda \mapsto -\overline{\lambda}$: every T-eigenvalue appears paired with its negative conjugate. For classical Hamiltonian tensors (mode-$k$ cubical tensors or with $k=2$), the sum of H-eigenvalues vanishes and the product equals $\det R$, where $R$ is the associated supersymmetric tensor [2503.21487].

A key structural result is the T-Williamson normal form: For a real symmetric positive-definite third-order tensor $\mathcal{M}$ ($\hat{\mathcal{M}}^{(i)}$ symmetric positive-definite for all $i$), there exists a T-symplectic tensor $\mathcal{S}$ and T-diagonal $\mathcal{D}$ such that
\[
\mathcal{M} = \mathcal{S}^H * \mathcal{D} * \mathcal{S}
\]
In the Fourier domain, each $\hat{\mathcal{D}}^{(i)} = \operatorname{blkdiag}(\Lambda^{(i)}, \Lambda^{(i)})$ with $\Lambda^{(i)} > 0$. The construction relies on slice-wise classical Williamson factorization, DFT and inverse DFT [2605.20829].

## 5. Stability, Analytical Criteria, and Numerical Methods

For tensor-based polynomial Hamiltonian systems, Lyapunov stability of an equilibrium $x^*$ is determined by the Hessian of the Hamiltonian:
\[
\nabla^2 H(x^*) = \sum_{j=2}^k j(j-1) B_j (x^*)^{j-2}
\]
An equilibrium is stable if this matrix is sign-definite. For $x^* = 0$, this reduces to examining $B_2$ [2503.21487].

Numerical computation of the T-Williamson form is efficient: it employs FFT on the third tensor mode, slice-wise classical Williamson factorization (cost per slice $O(n^3)$), and inverse FFT. Overall complexity is $O(p n^3)$ for a tensor of size $2n \times 2n \times p$ [2605.20829]. Numerical results confirm that the residuals of the defining relations are at machine precision and the runtime scales linearly in $p$ and cubically in $n$.

## 6. Illustrative Examples and Applications

Hamiltonian cubical tensors provide robust methods to identify and construct Hamiltonian structure in polynomial dynamical systems:
- For a cubic 2D system $\dot x_1 = x_1^2 + 2x_2$, $\dot x_2 = -2 x_1 x_2$, system tensors $A_2, A_3$ can be directly checked to be Hamiltonian cubical, yielding $H(x_1, x_2) = x_1^2 x_2 + x_2^2$ [2503.21487].
- In higher-dimensional examples, recovery of the polynomial Hamiltonian and stability analysis via the tensor-form Hessian are computationally tractable and more efficient than symbolic algebra as $n$ increases.
- In the T-product framework, applications include Fourier-domain encoding of covariance-matrix families as arise in continuous-variable quantum dynamics, enabled by T-Hamiltonian and T-symplectic tensors [2605.20829].

These results link the theory of cubical tensors with polynomial first integrals and generalized Killing tensors, as appear in covariant algorithms for higher-order invariants in nonrelativistic Hamiltonian systems [1404.3422].

## 7. Connections and Extensions

The classical theory of Hamiltonian matrices is fully recovered as the $k=2$ (matrix) case, while for higher-order ($k \geq 3$) tensors, Hamiltonian structure is captured by the interplay between symplectic algebra and tensor symmetries (“supersymmetry”). The block structure, spectral symmetry, and normal form results parallel those in symplectic geometry and linear algebra but require fundamentally multilinear generalizations.

Hamiltonian cubical tensors are essential in the analysis of integrability, conservation laws, and the construction of polynomial invariants in both finite and infinite-dimensional Hamiltonian dynamics. They also provide the algebraic backbone for constructive algorithms in numerical and symbolic computation of invariants and normal forms in modern tensor-based approaches to dynamical systems [2503.21487, 2605.20829].

Source: https://www.emergentmind.com/topics/hamiltonian-cubical-tensors