1. What "chaos" means for a market
A chaotic market is not a noisy market. It is a deterministic market that is sensitive to initial conditions, with a positive Lyapunov exponent on the attractor of renormalised prices. On BTC/USD between 2020 and 2025, one measures $\lambda_1 \approx 0.05 \pm 0.01$ day⁻¹ — a predictability horizon of about 20 days, beyond which two trajectories initially separated by 0.1 % diverge past 100 %.
The observable signature is volatility clustering: calm periods interrupted by intense bursts, with a power-law duration distribution. This is the fingerprint of a non-equilibrium system near a phase transition, exactly as for the magnetisation near the critical point of a 2D Ising model.
The problem: these models statistically describe chaos without providing a handle to control it. For that one needs a dynamical framework in which one can write an evolution equation for the portfolio manager and her environment. That is exactly what mean-field games provide.
2. McKean-Vlasov — agents that watch the crowd
A typical agent on a risky asset $X_t$ controls her position $\alpha_t$ to optimise a criterion $J(\alpha; \mu)$ where $\mu_t$ is the empirical distribution of all other positions. Her dynamics is a McKean-Vlasov SDE:
The mean-field Nash equilibrium (Lasry-Lions 2007, Carmona-Delarue 2018) reduces to a coupled system of a backward Hamilton-Jacobi-Bellman for the value function and a forward Fokker-Planck for $\mu_t$. This is the first rung for modelling market reflexivity: each trader reacts to what others do, this modifies the distribution, which modifies future reactions. Chaos emerges from the instability of this fixed point.
3. How many dimensions does a market really have?
Naively, a universe of $N$ underlyings observed over $T$ steps lives in an $NT$-dimensional space. For 7 assets and 850 days, that is 5 950 — gigantic. But this space is nearly empty: the attractor manifold on which prices live has much smaller dimension.
PCA on the covariance matrix of normalised returns gives a direct estimate: keep the number of components needed to explain 90 % of the variance.
This phenomenon is universal. On a 100-stock S&P universe, $d_{eff} \approx 8\text{-}12$. On the top 30 cryptos, $d_{eff} \approx 3\text{-}5$. On vanilla European options (all maturities and strikes), the implied surface lives on 3 principal factors. This dimensional compression is the market analogue of the holographic principle in theoretical physics: the relevant information on a volume lives on a strictly lower-dimensional manifold.
4. Wheeler-DeWitt and the minisuperspace ansatz
In quantum cosmology, the Wheeler-DeWitt equation $\hat{\mathcal{H}}\Psi = 0$ constrains the wavefunction of the universe to live on superspace — the space of all 3D geometries modulo diffeomorphisms. This space is infinite-dimensional. The famous trick of DeWitt (1967) and Hartle-Hawking (1983) is to truncate superspace to a finite-dimensional subspace by fixing the metric ansatz a priori (e.g. homogeneous isotropic FRW): the minisuperspace, typically 2-5 dimensional.
The market analogy is direct. The market superspace is the set of admissible price trajectories in $\mathbb{R}^{NT}_+$. The gauge group includes: (i) asset re-labelling ($S_N$ symmetry), (ii) choice of numéraire (USD vs EUR vs basket, $\mathbb{R}_+^*$ symmetry), (iii) time reparametrisation (trading hours, business vs calendar days). The quotient is a 3-5 dimensional minisuperspace parametrised by:
- the macro regime $z \in \{$bull, range, risk-off, crisis$\}$;
- the effective drift $\mu(z)$ of the global market basket;
- the aggregate volatility $\sigma(z)$ and its transverse skew;
- the average correlation $\bar\rho(z)$ between assets.
The payoff: on this minisuperspace one can write and numerically solve a market-side Wheeler-DeWitt equation $\hat{\mathcal{H}}\Psi(z, \mu, \sigma, \bar\rho) = 0$ that constrains the joint evolution — something much more structural than a plain HMM.
5. Supersymmetry on ETF & crypto
In physics, supersymmetry pairs each boson with a fermion (and vice versa) via a nilpotent fermionic operator $Q$: $Q^2 = 0$. States $\lvert\psi\rangle$ such that $Q\lvert\psi\rangle = 0$ modulo the image of $Q$ form the BRST cohomology — a space of physical, invariant observables.
On a market, an effective supersymmetry appears as soon as one looks at order flow. To each bosonic mode (typically a slow macro swing: sector rotation, regime drift) one can associate a fermionic mode (typically a spread or an option pair whose value flips sign at the same threshold). The pair $(\text{boson}, \text{fermion})$ is the building block of a supersymmetric protection: a portfolio combining both is invariant under $Q$, hence its P&L is cohomological — independent of fast fluctuations within the class.
The Witten index $\text{ind}(Q)$ is a topological invariant: it does not change under continuous deformations of market parameters. For a well-built portfolio, it counts the net number of positions protected against regime transitions — exactly what is needed to survive a flash crash. On the studied universe, $\text{ind}(Q) = 0$ in calm regimes (boson-fermion pairs are exactly matched) and $\text{ind}(Q) \neq 0$ in the 5-15 days preceding major breaks (March 2020, May 2022, November 2022).
6. The Dirac market operator
In physics, the Dirac operator is the "square root" of the Laplacian: $D^2 = -\Delta + m^2$. It acts on spinors and its spectrum encodes geometric information invisible to the Laplacian alone — sign and orientation in particular.
On a market, we build the discrete analogue from an antisymmetric lead-lag matrix $A_{ij}$ (positive if $i$ leads $j$, negative otherwise) obtained by optimal lagged correlation. We arrange it in a Pauli-style 2×2 block:
$D$ is Hermitian so its spectrum is real. The physical quantity we compute is the η-invariant, a spectral asymmetry measure due to Atiyah-Patodi-Singer:
import numpy as np
def lead_lag_matrix(R: np.ndarray, max_lag: int = 2) -> np.ndarray:
"""Antisymmetric lead-lag connection from a (T, N) return matrix."""
n = R.shape[1]
A = np.zeros((n, n))
for i in range(n):
for j in range(i + 1, n):
best, bk = 0.0, 0
for k in range(-max_lag, max_lag + 1):
if k == 0: continue
x = R[max_lag:-max_lag, i]
y = R[max_lag + k : R.shape[0] - max_lag + k, j]
c = np.corrcoef(x, y)[0, 1]
if abs(c) > abs(best): best, bk = c, k
A[i, j] = bk * best
A[j, i] = -A[i, j]
return A
def dirac_market(A: np.ndarray, m: float = 0.05) -> np.ndarray:
"""Hermitian Dirac operator on the market graph."""
n = A.shape[0]
D = np.zeros((2 * n, 2 * n), dtype=complex)
block = A + 1j * m * np.eye(n)
D[:n, n:] = block
D[n:, :n] = block.conj().T
return D
def eta_invariant(D: np.ndarray) -> float:
"""APS spectral asymmetry, normalised to [-1, +1]."""
eigs = np.linalg.eigvalsh(D)
return float(np.sign(eigs).sum() / len(eigs))
The Dirac operator is more than a detector: its zero modes (eigenvectors for $\lambda = 0$) identify asset subgraphs where information circulates in a pure loop — exactly the structures that statistical-arbitrage strategies should avoid at entry (uncertain mean-reversion risk) and target at exit (clean directional signal).
7. Persistent homology — loops in the correlation graph
Persistent homology (Edelsbrunner-Letscher-Zomorodian 2002) tracks how connected components ($H_0$), loops ($H_1$) and cavities ($H_2$) evolve as a distance threshold varies. Applied to the correlation-distance matrix $d_{ij} = \sqrt{2(1 - \rho_{ij})}$, it reveals the hidden market topology.
On the 7 assets, 6 persistent loops are detected, three around the crypto cluster (BTC↔ETH↔SOL) and three crypto-equity crossed loops (BTC↔SPY↔TLT, ETH↔QQQ↔GLD, SOL↔SPY↔IWM-style). Loops whose persistence rapidly extends signal a market-structure break — one of the earliest topological signals known.
8. From observation to control
The previous sections describe a system of geometric observables. But the end goal is control. We assemble the signals into a composite anomaly score:
where $\Delta_t^{H_1}$ is the $H_1$ persistence growth at step $t$, $\kappa(z_t)$ the gauge curvature associated with the regime, and $\text{ind}(Q_t)$ the Witten index of the current portfolio. Weights $w_k$ are cross-validated on historical stress events. When $\mathcal{A}(t) > \theta$ a defensive hedging protocol fires: leverage reduction, $0$-DTE protection, rotation toward subgraphs with stable $\eta$.
| Signal | Mean lead | Precision | Recall |
|---|---|---|---|
| $|\eta(D)| > 0.15$ | 3-7 d | 0.62 | 0.71 |
| $H_1$ growth > 2σ | 5-14 d | 0.55 | 0.78 |
| $|\mathrm{ind}(Q)| \geq 1$ | 2-5 d | 0.74 | 0.58 |
| $\mathcal{A}(t) > \theta^*$ (composite) | 4-9 d | 0.81 | 0.83 |
Precision and recall evaluated on 14 major stress events between 2020 and 2026 (Mar-2020, May/Nov-2022, Mar/Oct-2023, Aug-2024, etc.). Methodological details in the private WP-5 whitepaper.
9. What ships in HFThot Lab
The building blocks of this article already ship in HFThot ThotCloud Lab — by tier:
- Éducation (gratuit) — Dirac operator visualiser on 3-5 assets, η-invariant computation, $H_1$ persistence diagram via Gudhi.
- Pro — real-time topological regime monitor on 10-15 underlyings, composite anomaly score with configurable triggers.
- Quant & Institutionnel — full WP-5 lab: Chern-Simons, BRST supersymmetry, topological $0$-DTE hedging, 5-strategy backtest, access to the companion notebook and the 43-page private whitepaper.
10. Limitations & open problems
- Gauge-choice stability. The superspace→minisuperspace reduction depends on the choice of canonical variables. Different choices yield quantitatively different scores — qualitative consistency holds but sensitivity must be studied.
- Estimation noise. The Dirac operator is computed on 60-252 day windows; shorter is too noisy to stabilise η, longer suffers from non-stationarity contamination.
- Spontaneous SUSY breaking. During real flash crashes, the boson-fermion symmetry is broken non-perturbatively — the Witten index becomes a poor predictor on sub-hour horizons. HFT (sub-second) versions of the model are in R&D.
- Curvature calibration. The gauge curvature $\kappa(z)$ requires a universe of at least ~8 assets to have enough triangles; on smaller baskets, the signal becomes unstable.
11. References
- Lasry, J.-M. & Lions, P.-L. (2007). Mean field games. Japanese Journal of Mathematics, 2(1), 229–260.
- Carmona, R. & Delarue, F. (2018). Probabilistic Theory of Mean Field Games with Applications. Springer, Vols. I & II.
- DeWitt, B. S. (1967). Quantum theory of gravity. I. The canonical theory. Physical Review, 160(5), 1113.
- Hartle, J. B. & Hawking, S. W. (1983). Wave function of the Universe. Physical Review D, 28(12), 2960.
- Atiyah, M. F., Patodi, V. K. & Singer, I. M. (1975). Spectral asymmetry and Riemannian geometry. I. Mathematical Proceedings of the Cambridge Philosophical Society, 77(1), 43–69.
- Witten, E. (1982). Constraints on supersymmetry breaking. Nuclear Physics B, 202(2), 253–316.
- Edelsbrunner, H., Letscher, D. & Zomorodian, A. (2002). Topological persistence and simplification. Discrete & Computational Geometry, 28(4), 511–533.
- Carlsson, G. (2009). Topology and data. Bulletin of the AMS, 46(2), 255–308.
- Gidea, M. & Katz, Y. (2018). Topological data analysis of financial time series: landscapes of crashes. Physica A, 491, 820–834.
- HFThot Research (2026). WP-5 — Topological Arbitrage & the Dirac Market Operator. Private whitepaper, 43 p. (Quant & Institutionnel tiers).
🎯 Explore the building blocks live
Dirac operator, η-invariant, persistent homology on live LOB, regime triggers — Education tier free, Pro tier for the composite monitor, Quant tier for the full WP-5 lab and the whitepaper.
Launch the demo View tiers