Multidimensional Lorenz Attractor

• Multidimensional Lorenz attractor is a chaotic invariant set in high-dimensional systems characterized by robust pseudo-hyperbolic properties.
• It emerges through high-codimension bifurcations, such as homoclinic and heteroclinic tangencies, ensuring effective three-dimensional dynamics within complex phase spaces.
• Canonical models include discrete Hénon maps and sectional-hyperbolic flows, demonstrating clear geometric splitting with strong contraction and sectional expansion.

A multidimensional Lorenz attractor is a robust, genuinely chaotic invariant set arising in higher-dimensional dynamical systems—typically in discrete-time (diffeomorphism) or continuous-time (flow) settings—whose phase-space structure and bifurcation mechanisms generalize the classical Lorenz attractor from three to higher dimensions. These attractors arise naturally in the context of sectionally-hyperbolic and pseudo-hyperbolic systems, particularly through renormalization of local/global bifurcations (homoclinic or heteroclinic tangencies) at fixed or periodic points with appropriate degeneracies, and admit a tangent bundle splitting with uniform strong contraction in a high-codimensional subbundle, and robust sectional (two-dimensional) volume expansion in a lower-codimension complement. The organizing principle is that for a $d$-dimensional map or flow, the effective dynamics on the attractor is at least three-dimensional ($d_e=3$), ensuring the absence of global trapping by lower-dimensional invariant manifolds (Ovsyannikov, 2023, Gonchenko et al., 2014, Ovsyannikov, 2017, Ovsyannikov, 2021, Pacifico et al., 2022, Li et al., 2023, Karatetskaia et al., 2024).

1. Structure and Definitions

The fundamental feature of a multidimensional Lorenz attractor is a geometric structure supporting robust, non-uniformly hyperbolic chaos with pseudo-hyperbolic splitting:

• Strong stable subbundle $E^{ss}$ of dimension $d-2$, with $Df$ (or $Df_t$ for flows) contracting vectors exponentially;
• Center-unstable subbundle $E^{uc}$ (dimension 2), in which all two-dimensional subspaces are strictly volume-expanding;
• Absence of global invariant 2D manifolds: the attractor is not contained in any lower-dimensional foliation.

For discrete maps, the canonical model is a $d$-dimensional extension of the three-dimensional Hénon (discrete Lorenz) map: $$\begin{aligned} \bar{X}_1 &= X_2, \ \bar{X}_2 &= X_3, \ &\vdots \ \bar{X}_{d-1} &= X_d, \ \bar{X}_d &= M_1 + \sum_{j=2}^d M_j X_j + Q(X) - X_d^2, \end{aligned}$$ where $Q(X)$ is a quadratic form, and $\det D\bar{X} = 1$ ensures volume preservation at leading order (Gonchenko et al., 2014, Ovsyannikov, 2023).

For continuous flows, the prototypical model is the sectional-hyperbolic attractor—generalizing singular-hyperbolicity to higher dimensions. In this regime, $T_{\Lambda}M = E^{ss}\oplus F^{cu}$ with $E^{ss}$ uniformly contracted and $F^{cu}$ sectionally expanded (i.e., all 2-planes expand volume exponentially) (Pacifico et al., 2022, Li et al., 2023).

2. Bifurcation Mechanisms and Normal Forms

Multidimensional Lorenz attractors typically arise through high-codimension local or global bifurcations:

• Homoclinic/heteroclinic tangencies: At a $(d-1, 1)$-saddle or periodic orbit with $\text{codim}=d$, a non-simple quadratic tangency (non-transversal intersection) between invariant manifolds yields a first-return map whose rescaling and normal form converges to the multidimensional Hénon type (Gonchenko et al., 2014, Ovsyannikov, 2017, Ovsyannikov, 2021, Ovsyannikov, 2023).
• Degeneracies in multipliers: For flows or maps with multiple neutral or triple-resonant eigenvalues (e.g., $(-1, i, -i)$), triple/fourfold winged Lorenz-like attractors can emerge, as shown via codimension-3 normal forms and explicit $Z_4$-symmetric ODE normal forms (Karatetskaia et al., 2024).

The key parameter regime is structured by the requirement that the effective dynamics is genuinely three-dimensional (i.e., no lower-dimensional reduction due to global invariant manifolds), with the parameter space containing infinitely many domains (windows) converging to the organizing bifurcation (Gonchenko et al., 2014, Ovsyannikov, 2023).

3. Paradigmatic Models and Realizations

Model Type	Canonical Map/ODE	Effective Dimensionality	Reference
Discrete Hénon	$\bar{x}=y, \bar{y}=z, \bar{z}=M_1+B x+M_2 y - z^2$	3	(Ovsyannikov, 2023)
Sectional-hyperbolic flow	$T_\Lambda M=E^{ss}\oplus F^{cu}$	$n\geq3$	(Pacifico et al., 2022)
4D Lorenz-type DA-surgery	Special vector fields on $B^3\times I$	4	(Li et al., 2023)
$Z_4$-symmetric ODE	See equation (NF)	3	(Karatetskaia et al., 2024)

In the discrete case, the multidimensional Lorenz attractor is realized as a wild chaotic set of the multidimensional Hénon map, with robust pseudo-hyperbolic properties (one positive Lyapunov exponent, volume expansion in $E^{uc}$, and contraction in $E^{ss}$) (Gonchenko et al., 2014, Ovsyannikov, 2021).

For flows, e.g., four-dimensional DA-modified Lorenz flows, the attractor exhibits a chain-recurrence class with a unique Lyapunov-stable structure, robust heterodimensional cycles, and an exceptional, sectionally expanding 2D subbundle that is strictly embedded in a higher-dimensional dominated splitting (Li et al., 2023).

4. Dynamical, Statistical, and Geometric Properties

• Pseudo-hyperbolicity: Uniform contraction in $E^{ss}$ and volume (area) expansion in $E^{uc}$, verified by Lyapunov spectrum $$(\Lambda_1 > 0, \Lambda_2 \approx 0, \Lambda_3 \ll 0)$$ and by continuity of invariant subbundle splitting (angle criteria) (Karatetskaia et al., 2024, Ovsyannikov, 2021).
• Robustness: These attractors persist under $C^r$-small perturbations; the folding/stretching mechanism and volume expansion preclude the appearance of stable periodic sinks, making them wild or genuinely strange in the sense of Afraimovich–Bykov–Shilnikov (Gonchenko et al., 2014, Ovsyannikov, 2021).
• Existence and Uniqueness of Measures: In the context of flows, such as higher-dimensional sectional-hyperbolic attractors, there is a unique measure of maximal entropy for $C^1$-open/dense families (Pacifico et al., 2022). The unique non-atomic equilibrium state for Hölder continuous potentials is guaranteed when singularity measures are ruled out.
• Geometry: The phase space is foliated by strong stable leaves, and the global dynamics organizes around a "butterfly" shape, with orbits alternating between "wings" due to return map foldings (Gonchenko et al., 2014, Ovsyannikov, 2017, Karatetskaia et al., 2024).

5. Classification, Higher-Dimensional Generalizations, and Open Problems

Recent advances have resolved almost all codimension-3 bifurcation scenarios in dimension 3, incorporating not only resonance and saddle-focus cases but also Belyakov-type transitions and alternating leading stable subsets (3DL-bifurcations). For higher dimensions ($d > 3$), there is evidence and conjecture—supported by rescaling procedures, limiting renormalization, and the structure of the Jacobian—that similar wild attractors arise in $d$-dimensional Hénon-type maps and in flows with sectionally-hyperbolic splitting (Ovsyannikov, 2023, Ovsyannikov, 2021, Ovsyannikov, 2017, Gonchenko et al., 2014).

Open research questions include:

• Rigorous classification of all generic high-codimension tangency and resonance bifurcations in $d > 3$ that produce robust multidimensional Lorenz attractors.
• Systematic study of limiting maps with more than two neutral multipliers, possibly leading to $4$D or higher pseudo-hyperbolic attractors (Ovsyannikov, 2017).
• Computer-assisted proofs of wild hyperbolicity and robust mixing properties in multidimensional Lorenz maps, especially for parameter regimes near resonance points (Ovsyannikov, 2023).
• Applied realization in multi-mode physical systems (e.g., lasers, coupled oscillators) where such high-dimensional bifurcation networks may be observable in experiments (Ovsyannikov, 2017).
• Understanding exceptional behaviors such as the existence of robustly heterodimensional cycles without homoclinic tangencies and new partially hyperbolic phenomena for $n\geq4$ (Li et al., 2023).

6. Numerical, Analytical, and Experimental Characterization

Numerical studies in both normal forms and discrete quadratic models demonstrate existence of open regions in parameter space supporting robustly chaotic, pseudo-hyperbolic attractors, confirmed by computation of maximal Lyapunov exponents and angles between invariant bundles (Karatetskaia et al., 2024). Analytical construction relies on normal form reductions, renormalization/rescaling lemmas, geometric models, and local/global bifurcation theory. There is a systematic methodology for reducing neighborhood dynamics near high-degeneracy bifurcations to a universal Hénon-type prototype (Gonchenko et al., 2014, Ovsyannikov, 2023).

Experimental confirmation remains highly challenging, but the theoretical robustness under perturbation suggests that such structures could, in principle, be realized in higher-dimensional dissipative or conservative physical systems, especially those admitting high-codimension bifurcations or symmetries leading to multiple neutral directions (Ovsyannikov, 2017).

---

In summary, the theory of multidimensional Lorenz attractors unifies a large class of robust, wild, structurally stable chaotic sets arising in $n$-dimensional systems—constructed via high-codimension bifurcations and invariant decomposition mechanisms generalizing the classical Lorenz model. This area remains a focal point for research on high-dimensional chaos, singular flows, pseudo-hyperbolicity, nonlinear resonance, and the topological/dynamical classification of hyperchaotic attractors (Ovsyannikov, 2023, Gonchenko et al., 2014, Ovsyannikov, 2021, Ovsyannikov, 2017, Pacifico et al., 2022, Li et al., 2023, Karatetskaia et al., 2024).