Machine Learning Invariants of Tensors

The Big Picture

Imagine describing the shape of a crystal. No matter how you rotate or flip it, certain properties stay the same: the number of faces, the angles between them, the overall symmetry. These unchanging quantities are called invariants, and physicists have been hunting them for over a century.

The problem scales badly. A single electromagnetic field has two independent invariants. A gravitational field in ten dimensions? Dozens, tangled in hidden relationships. Sorting which invariants are truly independent and which are disguised combinations of others has traditionally required years of mathematical heroics.

A team of researchers has automated much of that heroism. Their algorithm uses random sampling, numerical linear algebra, and graph enumeration to discover the independent invariants of any tensor (a structured mathematical object, like a matrix but with potentially more dimensions). It turns a centuries-old grind into an overnight computation.

Their test case: an antisymmetric 3-form field in six-dimensional space, a structure prominent in string theory. The algorithm confirmed it has exactly five independent invariants, a result that constrains the entire space of physical theories built from that field.

Key Insight: By generating random instances of a tensor and hunting for linear dependencies between candidate scalars, this algorithm can experimentally identify which invariants are truly independent, turning a hard algebraic problem into a tractable numerical one.

How It Works

The core idea is elegant: compute two quantities for thousands of random inputs; if some weighted sum always equals zero, the quantities are secretly related. The algorithm runs in three stages.

Stage 1: Graph enumeration. Every number computable from a tensor, called a scalar, can be drawn as a tensor network: a diagram where nodes represent the tensor’s components and lines represent pairs of indices being summed. The algorithm systematically enumerates all inequivalent diagrams at each order, building a complete candidate list. For a rank-3 antisymmetric tensor (one where swapping any two indices flips the sign), the count at low orders is manageable. Higher orders grow combinatorially, but the graphs remain systematic.

Stage 2: Numerical evaluation. For each candidate scalar, the algorithm draws many random instances of the tensor, with components sampled from Gaussian distributions, and computes every invariant for every draw. This produces a large matrix: rows are random samples, columns are candidate invariants.

Stage 3: Linear algebra. The team searches for syzygies, linear relations of the form $c_1 I_1 + c_2 I_2 + \ldots = 0$ holding across all draws. Numerically, this means finding the null space of the data matrix. A nonzero null vector reveals a genuine algebraic identity: one invariant is a combination of the others. When no new independent invariants appear across several consecutive orders, the list is considered complete.

The 3-form test case. The team applied the algorithm to an antisymmetric 3-form $H_{\mu\nu\rho}$ in six dimensions. This is not an exotic curiosity; it appears centrally in string theory and supergravity as the B-field of the NS-NS (Neveu-Schwarz) sector. The algorithm builds scalars from $H$ three ways:

• Trace variables: Multiply copies of $H$ and sum over shared indices using the metric, forming combinations like $H_{\mu\nu\rho}H^{\mu\nu\rho}$ and higher-order versions.
• Hodge dual variables: Use the six-dimensional Levi-Civita symbol $\epsilon_{\mu\nu\rho\sigma\lambda\kappa}$ to construct a “mirror image” of $H$, then build invariants from that.
• Spinor variables: Rewrite the tensor in spinor indices via six-dimensional Weyl matrices, revealing invariants invisible in the original picture.

After running the algorithm across all three families, the result was unambiguous: only five invariants are truly independent. All others can be expressed as explicit algebraic functions of these five, with exact formulas provided in the paper.

This has immediate consequences. Any Lagrangian for a 3-form in six dimensions (the master function from which all equations of motion and physical predictions derive) is a function $L(x_1, x_2, x_3, x_4, x_5)$ of exactly five variables. No more, no less. Theories that look different on the surface may secretly be the same theory in disguise.

Why It Matters

In string theory, understanding the independent invariants of the Riemann tensor in ten dimensions is essential for classifying higher-derivative corrections to supergravity, the so-called $\alpha’$ expansion. Traditional approaches require painstaking manual algebra that experts have labored over for decades. Automating invariant-finding could speed that work up considerably.

In quantum information, the independent invariants of a multi-partite quantum state under local unitary transformations are exactly the entanglement measures. Identifying them is the mathematical backbone of classifying what kinds of entanglement are genuinely distinct. The same algorithm, pointed at a different tensor, could push that classification further.

More broadly, the work shows that numerical linear algebra and random sampling can serve as experimental tools in pure mathematics. Not approximations, but hypothesis generators later confirmed analytically. That’s a mode of doing mathematics that deserves more attention.

Bottom Line: By treating invariant theory as a data problem, this team confirmed that a 6D antisymmetric 3-form has exactly five independent invariants and built a reusable algorithm that can answer the same question for any tensor, with applications in string theory, quantum information, and algebraic geometry.