Introduction to Supersymmetric Theory of Stochastics

Many natural and engineered dynamical systems, including all living objects, exhibit signatures of what can be called spontaneous dynamical long-range order (DLRO). This order's omnipresence has long been recognized by the scientific community, as evidenced by a myriad of related concepts, theoretical and phenomenological frameworks, and experimental phenomena such as turbulence, $1/f$ noise, dynamical complexity, chaos and the butterfly effect, the Richter scale for earthquakes and the scale-free statistics of other sudden processes, self-organization and pattern formation, self-organized criticality, etc. Although several successful approaches to various realizations of DLRO have been established, the universal theoretical understanding of this phenomenon remained elusive. The possibility of constructing a unified theory of DLRO has emerged recently within the approximation-free supersymmetric theory of stochastics (STS). There, DLRO is the spontaneous breakdown of the topological or de Rahm supersymmetry that all stochastic differential equations (SDEs) possess. This theory may be interesting to researchers with very different backgrounds because the ubiquitous DLRO is a truly interdisciplinary entity. The STS is also an interdisciplinary construction. This theory is based on dynamical systems theory, cohomological field theories, the theory of pseudo-Hermitian operators, and the conventional theory of SDEs. Reviewing the literature on all these mathematical disciplines can be time consuming. As such, a concise and self-contained introduction to the STS, the goal of this paper, may be useful.


A. Dynamical Long-Range Order
It is well established experimentally and numerically that many seemingly unrelated sudden processes in astrophysics [1], geophysics [2], neurodynamics [3,4], econodynamics [5], and other branches of modern science exhibit power-law statistics, the very reason why the Richter scale is logarithmic. This is simply one example of the spontaneous long-range dynamical behavior (LRDB) that emerges in many nonlinear dynamical systems (DSs) with no underlying long-range interactions that could potentially explain such behavior. Two other well-known examples of LRDB are the infinitely long memory of perturbations known as the butterfly effect [6], and the algebraic power-spectra commonly known as 1/f noise or the long-term memory effect [7] found in many existing DSs, including apparently all living objects [8,9].
It was understood that the LRDB must be a signature of some type of spontaneous dynamical long-range order (DLRO). The existence and omnipresence of this DLRO has long been recognized by the scientific community, as evidenced by a myriad of related concepts, including chaos [10][11][12], turbulence [13,14], dynamical complexity [15], self-organization [16], pattern formation [17], and self-organized criticality [18].
Several successful approaches to various realizations of DLRO have been established. For example, the concept of deterministic chaos is a centerpiece of the well-developed dynamical system (DS) theory. Nevertheless, there existed no universal theoretical understanding of DLRO. In particular, no rigorous stochastic generalization of the concept of deterministic chaos existed previously, whereas all natural DSs are never completely isolated from their environments and are thus always stochastic.
A class of models with the potential to reveal the mathematical essence of the ubiquitous DLRO is the stochastic (partial) differential equations (SDEs). Indeed, SDEs most likely have the widest applicability in modern science. In physics, for example, SDEs are the effective equations of motion (EoM) for all physical systems above the scale of quantum degeneracy/coherence. In quantum models, SDEs are used in a variety of ways. For example, SDEs play a central role in quantum optics (see, e.g., [19] and the references therein). In many-body quantum models, SDEs are used in the investigation of non-equilibrium quantum dynamical phenomena in the form of the effective EoM of the collective quantum modes [20] and order parameters [21]. They also represent a useful tool for quantum statistics [22]. In other scientific disciplines, SDEs are even more fundamental, as they appear at the level of the very formulation of dynamics, unlike the EoM in physics, which descend from least action principles.
Nevertheless, the mathematical essence of DLRO remained elusive.

B. Topological Supersymmetry of Continuous Time Dynamics
One way of deducing the potential theoretical origin of DLRO is provided by the following qualitative yet solid argument. From the field-theoretic point of view, LRDB is indicative of the presence of a gapless excitation with an infinite correlation length/time. There are only two possible scenarios for such a situation to occur: the accidental or critical scenario and the Goldstone scenario. In the accidental scenario, the parameters of the model can be fine-tuned to ensure that a certain excitation has zero gap. This is exactly the situation with (structural) phase transitions, where an excitation called the soft mode becomes gapless exactly at the transition temperature (or other parameter). This allows the system to move effortlessly from a previously stable vacuum to a new vacuum. Immediately following the transition, the soft mode "hardens" again, i.e., it acquires a finite gap that signifies the dynamical stability of the new vacuum. In other words, only at exactly the transition point the soft mode is gapless and thus has an infinite correlation length/time.
The accidental scenario for DLRO contradicts the fact that DLRO is robust against moderate variations in the parameters of the model. For example, a slight variation in the magnitude of the electric current flowing through a dirty conductor will not destroy the 1/f noise. In other words, in phase diagrams, DLRO occupies full-dimensional phases and not the lower dimensional transitions/boundaries between different full-dimensional phases.
This observation unambiguously suggests that the Goldstone scenario is the only possibility for the field-theoretic explanation of DLRO. More specifically, the Goldstone theorem states that, under the conditions of the spontaneous breakdown of a global continuous symmetry, the ground state is degenerate and that, in spatially extended models, this degeneracy tailors the existence of a gapless excitation called the Goldstone-Nambu particle. As such, DLRO may be the result of the spontaneous breakdown of some global continuous symmetry.
It is understood that the symmetry responsible for DLRO cannot be a conventional bosonic symmetry because DSs with no bosonic symmetries can also exhibit DLRO, e.g., be chaotic. In other words, DLRO must be a result of the spontaneous breakdown of some fermionic symmetry or supersymmetry. It has long been known that supersymmetries are indeed present in some classes of SDEs. The work on supersymmetric theories of SDEs began with the Parisi-Sourlas stochastic quantization procedure [31][32][33][34][35][36][37][38][39][40][41], which leads from a Langevin SDE, i.e., an SDE with a gradient flow vector field, to a model with N = 2 supersymmetry. The Parisi-Sourlas quantization procedure was later identified as a realization of the concept of Nicolai maps [42,43] and "half" of the N = 2 supersymmetry as a corresponding Becchi-Rouet-Stora-Tyutin (BRST) or topological supersymmetry, which is a definitive feature of Witten-type topological or cohomological field theories [44][45][46][47][48][49][50].
From the perspective of the theory of ubiquitous DLRO, the consideration of specific models is clearly insufficient. In reality, EoM are never exactly Langevin or classical mechanical, and the generalization of the discussion to all or at least most of SDEs is necessary.
In other words, the supersymmetry responsible for DLRO must be an attribute of all SDEs to be able to account for omnipresence of DLRO in nature. The remaining question in the Goldstone scenario of the theory of DLRO now is whether such supersymmetry exists.
Traces of this supersymmetry can be found in the literature. In [70], for example, the authors considered a non-potential generalization of the Langevin stochastic dynamics and noted that half of the N = 2 supersymmetry survives a non-potential perturbation. Nevertheless, to the best of the knowledge of the present author, this N = 1 supersymmetry in the context of all SDEs has not been addressed previously. One possible reason for this is the pseudo-Hermitianity of the stochastic evolution operator of a general SDE. Specifically, the theory of pseudo-Hermitian evolution operators appeared only relatively recently [60][61][62][63][64] as a generalization of the theory of P T -symmetric evolution operators [65][66][67][68][69]. It was only after the theory of pseudo-Hermitian operators became available that studies on topological supersymmetry in the context of the general SDEs could be resumed. The idea that the spontaneous breakdown of this supersymmetry pertinent to all SDEs may be the mathematical essence of DLRO, or rather of one of its realizations known previously as self-organized criticality, was reported in [71]. Further work in this direction [72][73][74][75] resulted in the formulation of what can be called the supersymmetric theory of stochastics (STS). The goal of this paper is to present the current state of the STS in a self-consistent manner. This paper can be viewed as a compilation of a few previous works and as a compilation that corrects several mistakes made during the early stages of the development of the STS and that clarifies a couple of points that were previously swept under the carpet. This paper also presents a few new results, including a discussion of the pseudo-time reversal symmetry.
Given the multidisciplinary character of STS, it would take an enormous amount of work to review all the relevant results from DS theory, cohomological field theory, the classical theory of SDEs, and physics. This goal is not pursued in this paper, and references are provided on only the most relevant results that are known to the author and that the material presented here is directly based on. The author would like to apologize in advance if some important related works have escaped his attention.

C. Relation to Existing Theories
The topological supersymmetry breaking picture of DLRO aligns well with the previous understanding of the concept of dynamical chaos. For example, the nontrivial connection between chaos and topology is at the heart of the topological theory of chaos [76]. Furthermore, it was also known that, in some cases, the transition into chaos must be a phase transition of some sort, as evident from certain universal features of the onset of chaotic behavior [77]. The only unexpected insight from the supersymmetry breaking picture of DLRO is the fact that its mathematical essence is in a sense opposite to the semantics of the word chaos. Indeed, chaos literally means "absence of order", whereas the phase with the spontaneously broken supersymmetry is the low-symmetry or "ordered" phase. This is why DLRO may be a more accurate identifier for this phenomenon than, say, stochastic chaos.
In this paper, both terms will be used interchangeably.
STS in a nutshell is the following. An SDE defines the noise-configuration-dependent trajectories in the phase space. The collection of all these trajectories can be viewed as a family of noise-configuration-dependent phase space diffeomorphisms. Instead of studying the trajectories, one can equivalently study the actions, called the pullbacks, that these diffeomorphisms induce on the exterior algebra of the phase space. The original trajectories can be reconstructed from these pullbacks so that the later contain all the information on the SDE-defined dynamics.
The pullbacks have one very important advantage over the trajectories. Unlike trajectories in the general case of a nonlinear phase space, the pullbacks are linear objects and can thus be averaged over the noise configurations. Such a stochastically averaged pullback is the finite-time stochastic evolution operator (SEO). Thus, it becomes clear where the supersymmetry originates from: all the diffeomorphism-induced pullbacks and consequently the finite-time SEO are commutative with the exterior derivative, which is thus a (super-)symmetry of any SDE. In other words, the existence of this supersymmetry in all SDEs is merely the algebraic version of the most fundamental and indisputable statement that continuous dynamics preserves the continuity of the phase space.
The idea to study pullbacks induced by random maps averaged over noise configurations appeared first, to the best of this author's knowledge, in DS theory, where the analogue of the finite-time SEO is known as the generalized transfer operator [78]. From this perspective, the STS can be viewed as a continuation of DS theory. On the side of the Parisi-Sourlas quantization procedure, the path integral representation of the Witten index of the STS is a member of the cohomological field theories. Furthermore, the SEO of the general SDE is pseudo-Hermitian; thus, the STS is within the domain of applicability of the theory of pseudo-Hermitian operators. In other words, the STS is a multidisciplinary mathematical construction. It combines a few major mathematical disciplines that are naturally synergetic within the STS. This synergy can ensure fruitful cross-fertilization during future work on the STS. To date, the STS has already provided a few novel findings, therein making it interesting from several points of view, as discussed below.
For DS theory, an interesting result from the STS is the established equivalence between the so-called sharp-trace of the generalized transfer operator, the stochastic Lefschetz index of the corresponding SDE-defined diffeomorphisms, and the Witten index of the STS.
From the perspective of the conventional theory of SDEs, a valuable result from the STS is the demonstrated equivalence between the Stratonovich interpretation of SDEs and the (bi-graded) Weyl symmetrization procedure. For a field theorist, there are two potentially interesting results from the STS. First, the cohomological field theories, or rather the methodology developed within them (e.g., the localization principle and topological invariants as expectation values on instantons), together with the theory of pseudo-Hermitian evolution operators, may find multiple applications in almost all branches of modern science. Second, there are very few known analytical mechanisms that can result in the spontaneous breakdown of supersymmetry [80], which is basically one of the main reasons behind the introduction of the concept of explicit (or soft) supersymmetry breaking [81]. The STS provides yet another such mechanism: the topological supersymmetry in (deterministic) chaotic DSs is spontaneously broken by the non-integrability of the flow vector field.
From a wider perspective, SDEs find applications in almost all modern scientific disciplines, ranging from social sciences and econodynamics to astrophysics and high-energy physics. Therefore, the STS in general and this paper in particular may be interesting to specialists working in any of these areas of science.

D. Models of Interest and the Structure of This Paper
The following class of SDEs that covers most of the models in the literature will be of primary interest:ẋ (t) = F (x(t)) + (2Θ) 1/2 e a (x(t))ξ a (t) ≡ F(t). (1) Here and in the following, summation over repeated indexes is assumed; x(t) : R → X is a trajectory of the DS in a D-dimensional topological manifold called the phase space, X; F (x) ∈ T X x is the flow vector field from the tangent space of X at the point x; ξ = {ξ a ∈ R, a = 1, 2...} are noise variables; and e a (x) ∈ T X x is a set of vector fields. The positiondependent/independent e are often called multiplicative/additive noise. The notation F is introduced to separate the flow perturbed by the noise from the deterministic flow, F . As will be discussed in Section III B, the SDE in Equation (1) is the Stratonovich SDE along the lines of stochastic calculus on manifolds (see, e.g., [82] and the references therein). It will also be argued that the STS appears to point to the possibility that the Stratonovich approach is the only correct choice for continuous time models.
The parameter Θ represents the temperature or rather the intensity of the noise. As will be made clear below in Section III, the vector fields e a define the noise-induced metric on X: g ij (x) = e i a (x)e j a (x). Therefore, in situations wherein the number of vector fields e a equals the dimensionality of the phase space, these vector fields can be identified as veilbeins (see, e.g., Chapter 7 of [83]). In the general case, however, the number of es must not necessarily be equal to the dimensionality of the phase space.
Most of the discussion will be directed toward models with Gaussian white noise. The probability distribution of its configurations is with C being a normalization constant such that Here, the functional or infinitely dimensional integration is over all the configurations of the noise. The stochastic expectation value of some functional f (ξ) is defined as The fundamental correlator of the Gaussian white noise is The infinite-dimensional integrations in Equations (4) and (5) can be given a more concrete meaning by splitting the time domain into a large number of intervals with infinitesimal duration ∆t and then taking the continuous time limit, namely, ∆t → 0. Before taking this limit, each noise configuration can be viewed as a piece-wise constant function (see Figure   1) on each interval, i.e., the value of the noise variable ξ a (t) = ξ a n for t n > t > t n−1 . The discrete-time version of the probability distribution of the Gaussian white noise in Equation (2) is P Ns (ξ) ∝ e −∆t n ξ a n ξ a n /2 , and that of the correlator in Equation (5) is ξ a n ξ b n Ns = ∆t −1 δ ab δ nn , whereas all the other (even) order correlators are The theory of stochastic dynamics defined by Equation (1) can be constructed in two steps. The first step is to understand the deterministic temporal evolution defined by the ordinary differential equation (ODE) obtained from the SDE in Equation (1) by fixing the noise configuration. This problem will be addressed in Section II, where a few concepts closely related to the continuous-time dynamics will also be introduced. The second step is the stochastic generalization of this deterministic evolution, which will be addressed in Section III.
The realistic noises are more complicated than Gaussian white noise, which is, of course, a mathematical idealization. In Section IV, the path integral representation of the theory will enable the generalization to noise of any form. Further generalization to the spatially extended models with infinite-dimensional phase spaces will also be discussed briefly in Section IV C. Having established general technical aspects of the STS, the discussion will concentrate on the analysis of the structure of the ground states in Section V. The classification of ergodic stochastic models on the most general level related to topological supersymmetry breaking will be proposed. This in particular will help reveal the theoretical picture of the stochastic dynamics on the border of "ordinary chaos", known previously under such names as intermittency, complexity, and self-organized criticality. Finally, in Section VI, the paper will be concluded with a brief discussion of a few potentially fruitful directions for future work.

A. Dynamics as Maps
For a fixed noise configuration, Equation (1) is an ODE with a time-dependent flow vector field in its Right-Hand Side (R.H.S.). This ODE defines a two-parameter family of maps of the phase space onto itself, namely, M tt : X → X: These maps have straightforward interpretations: x(t) = M tt (x ) is the solution of the ODE with the condition x(t ) = x . Clearly, The only difference here with the stationary flows described, e.g., in Chapter 5 of [83] is that the maps depend on both the initial and final moments of evolution, i.e., t and t, and not only on the duration of the evolution, i.e., t − t . This is the result of the dependence of the noise configuration on time, which breaks the time-translation symmetry. Following stochastic averaging over the Gaussian white noise, which does possess time-translation symmetry, this symmetry of the model will be restored (see Section III).
Only physical models in which the maps (for finite time evolution) are invertible and differentiable will be considered. On the mathematical level, this means that F and e's are sufficiently smooth in X such that the Picard-Lindelöf theorem (see, e.g., [84]) on the existence and uniqueness of the solution of an ODE for any initial condition is applicable.
In other words, all maps are diffeomorphisms.
To avoid the necessity of addressing various subtle mathematical aspects not directly related to the subject of interest, the fixed noise configuration will be assumed as a continuous function of time. However, this continuity is not necessary. The noise configuration only needs to be integrable in the sense that there must exist a W a (t) such that dW a (t)/dt = ξ a (t).
For Gaussian white noise, W a (t) is called the Wiener process.
A physicist's proof of the invertibility of maps defined by Equation (10) is as follows. A physical ODE provides only one outcome x at t for each initial condition x at t . The same must be true for the time-reversed physical ODE, which provides only one x at t for any x at t. In other words, the map M tt is a one-to-one map, i.e., it is invertible.
If at time t the DS is described by a total probability function P (x), the expectation According to Equation (9), this expectation value at a later time moment t > t is This view on dynamics can be clarified through the following example. Consider X = R D and an ODE of the simple formẋ = v, where v ∈ R D is a constant vector field. The For f (x) being one of the coordinates, i.e., f (x) = x i , Equation (12) states that One can now make the transformation of the variable of integration in Equation (12), i.e., x → M t t (x), Here, M * t t is the operation of the variable transformation applied to the coordinate-free object consisting of P (x) and the collection of all the differentials dx 1 ...dx D , where J is the Jacobian of the tangent map, with being the coordinate representation of the tangent map.
Equation (13) suggests that the forward temporal evolution of the variables of the DS is equivalent to the backward temporal evolution of the coordinate-free object representing the total probability distribution (TPD) In algebraic topology, this object is known as a top differential form (D-form), the infinitedimensional linear space of all D-forms is denoted as Ω D (X), and the operation M * t t in Equation (14) is called the action or the pullback induced by M t t on ψ (D) .
Note that the diffeomorphism in Equation (13) is for the inverse temporal evolution as compared to the time flow in the SDE. This seeming confusion of the time direction can be clarified as follows. The pullbacks act in the opposite direction compared to the diffeomorphisms inducing them. This is the reason for the term pullback. The graphical representation of this situation is given in Figure 2. There, one introduces an infinite number of copies of the phase space for each time moment, X(t), and dynamics is defined as a twoparameter family of diffeomorphisms between these copies: M tt : Continuous-time deterministic dynamics with a fixed noise configuration can be viewed as a two-parameter family of diffeomorphisms of the phase space onto itself or between the copies of the phase space: M tt : X(t ) → X(t). The temporal evolution of a differential form is a pullback In this path-integral-like picture of dynamics, the pullback in Equation (14) can be given as or as opposed to The relation between the direction of the flow of time for maps and the corresponding pullbacks can be expressed via the following diagram: This diagram particularly suggests that the composition law for pullbacks is

B. Differential Forms as Wavefunctions
The description of a stochastic model in terms of only TPDs as in the previous subsection is insufficient in the general case. This can be observed from the following qualitative example. Consider the simplest Langevin SDE with X = R, F = ∂U (x)/∂x, e = 1.
Consider also the case of the stable Langevin potential U , as shown in Figure 3a. It is clear that, after a sufficiently long temporal evolution, this DS will forget its initial condition, and its (only) variable will be distributed according to some steady-state TPD, which is the ground state of this DS (see Section V B for details). In contrast, when the Langevin potential is unstable, as in Figure 3b, the DS will never forget its initial condition because a small difference in the initial conditions will grow exponentially. No meaningful steady-state TPD can be prescribed to its unstable variable. This example signifies that the steady-state probability distributions make sense only for stable variables.
b ) i.e., U (x) (blue parabola oriented up), exhibits dynamics (broken dashed arrow) that can be characterized as the gradual settling to a steady-state probability distribution P (x) (bell-shaped curve). These dynamics exhibit a loss of the dynamical memory of the initial condition; (b) In the case of an unstable Langevin potential (blue parabola oriented down), the dynamics escape to infinity (dashed arrow pointing left). The dynamics is sensitive to the initial condition. No meaningful steady-state probability distribution can be associated with the ground state in this case.
The previous example may not look physical because, for any initial condition, the DS escapes to infinity and never returns. Perhaps a better example for the same purpose is a (deterministic) chaotic DS, in which the unstable variables exist even after the infinitely long temporal evolution, i.e., even in the ground state of the DS. In DS theory, the existence of these unstable variables is revealed by positive (global) Lyapunov exponents. Such a chaotic ground state must not be a probability distribution in its unstable variables. That this is indeed so will be observed in Section V C 2 below. The DS theory predecessors of such ground states are the Sinai-Ruelle-Bowen conditional probability functions on the global unstable manifolds [85].
Section III F will demonstrate on a more rigorous level that it is a mathematical necessity that the Hilbert space of a stochastic DS be not only the space of the TPDs but rather the entire exterior algebra of X: with the elements being the differential forms of all degrees (see, e.g., Chapter 5 of [83]) Here, 0 ≤ k ≤ D, ψ is an antisymmetric tensor, ∧ is the wedge or antisymmetrized product of differentials, e.g., dx 1 ∧dx 2 = −dx 2 ∧dx 1 = dx 1 ⊗dx 2 −dx 2 ⊗dx 1 , and Ω k (X) is the space of all differential forms of degree k (k-forms).
This by no means contradicts the intuitive understanding that it must be possible to associate a TPD with any wavefunction. As will be clear later, the TPD associated with a wavefunction is not the wavefunction itself but rather, as in quantum theory, is the bra-ket combination, which is a D-form and/or a TPD (see, e.g., Section III C and the discussion following Equation (96)).
One possible interpretation of the differential forms is the generalized (total, conditional, marginal) probability distributions in the coordinate-free setting. The following example demonstrates how the conditional probability distribution can be represented as a differential form (the dimensionality of X is D = 3): 13 = P (x 1 x 3 |x 2 ), and ψ (2) 32 = P (x 2 x 3 |x 1 ). Similarly, the TPD introduced previously is being the Levi-Civita antisymmetric tensor and with p(i 1 ...i D ) being the parity of the permutation of indexes.
The geometrical meaning of a k-form is a differential of a k-dimensional oriented volume.
Therefore, a k-form can be integrated over a k-dimensional submanifold or a k-chain, c k , This quantity can be interpreted as follows. If one introduces local coordinates such that the k-chain belongs to the k-dimensional manifold cut out by (x k+1 , ..., dx D ) = (Const (k+1) , ..., Const (D) ), then Equation (25) is the probability of finding variables (x 1 , ..., x k ) within this k-chain given that all the other variables are known with certainty to be equal (Const (k+1) , ..., Const (D) ).
It is worth stressing that the interpretation of the differential forms as the generalized probability distributions is valid only locally in the general case. Only in a neighborhood of a given point and with a properly chosen local coordinates can a differential form be thought of as a conditional probability distribution. Globally, however, this may not be possible because there may not exist global coordinates such that a given differential form is positive everywhere on X and normalizable. Moreover, if it were possible to interpret all differential forms as conditional probability distributions in the global sense, then there would be no reason to consider the extended Hilbert space in the first place. Indeed, a conditional probability distribution can be constructed from the TPD so that it does not contain any additional information, and it would suffice to describe the DS in terms of the TPD only.
The exact physical meaning of the wavefunctions in the STS is an open question. At this moment, as a working interpretation of the wavefunction, one can adopt the point of view on the wavefunction from quantum theory. Namely, the ket of the wavefunction at a given moment of time is an abstract object that contains information about the system's past, whereas the bra-ket combination of a wavefunction has the meaning of the TPD.
To finalize the above justification for the use of the extended Hilbert space, it must be stressed that the idea of using the entire exterior algebra as a Hilbert space of a DS is by no means new. This is a well-known method in the supersymmetric theory of Hamilton models in references [52][53][54], where it was even demonstrated to a certain degree that the information of chaoticity of a Hamilton model is better represented by differential forms.
Moreover, the mathematical object known as the generalized transfer operator that will play a central role in Section III was designed in the DS theory to probe chaos, and this object was defined on the entire exterior algebra [78].
To establish the law of the temporal evolution of k-forms, one assumes that the DS is described by ψ (k) (x) at time moment t . By analogy with Equation (13), the quantity in Equation (25) at a later time moment t > t is Here, M * t t : Ω k (X) → Ω k (X) is the generalization of the pullback in Equation (14) to pullbacks acting on Ω (k) (X). Explicitly, where the k-form is from Equation (24) and d(M t t (x)) i is from Equation (15).
C. Operator Algebra

Lie Derivative
The infinitesimal pullback is known as the physical or Lie derivativê The infinitesimal map defined by Equation (1) can be given as with F being the R.H.S. of Equation (1). Accordingly, the infinitesimal tangent map defined in Equation (16) is with Using Equations (28)- (30) and the definition of the pullback in Equation (27), one arrives at the following expression for the Lie derivative: with ψ (k) being from Equation (24).
The finite-time pullback satisfies the following equation: where , which follows from Equation (22), has been used together with the definition of the Lie derivative in Equation (28). The integration of this equation where T denotes the operator of chronological ordering. This operator is necessary becausê L F (τ ) at different τ s do not commute. Equation (34) can be represented in the form of a Taylor series as

Exterior Derivative
One of the fundamental operators of the exterior algebra is the exterior multiplication dx i ∧ : Ω k (X) → Ω k+1 (X). This operator can be defined via its action on a k-form from Equation (24): Viewing the differentials in the definition of a k-form in Equation (24) as the operators of exterior multiplication, one can also define the operation of the exterior product of differential forms: The other fundamental operator of the exterior algebra is the interior multiplication where dx i j denotes a missing element. As can be readily verified, the (anti)commutation relations for these operators are Here and in the following, the square brackets denote the bi-graded commutator: with deg(X) = #(dx∧) − #(î) being the degree of operatorX, i.e., the difference be- The centerpiece of the theory under consideration is the exterior derivative or de Rahm The exterior derivative is a bi-graded differentiation, i.e., for any operatorsX andŶ , In the new notations, the Lie derivative can be given via the Cartan formulâ whereî F = F iî i is the interior multiplication by F, and [d,

Hodge Dual
Yet another operation that will be used later is the Hodge star defined as where is the Levi-Civita antisymmetric tensor and g = det g ij is the determinant of the metric on X. As previously mentioned, the natural choice of metric on X is the noiseinduced metric g ij = e i a e j a . In Section III C 3, it will be noted that, for certain purposes, other metrics on X can be used. The Hodge star has the following property: wherek is the operator of the degree of the differential form so that This can easily be verified using where the indexes are lowered by the Euclidean metric and In other words, the square of the Hodge star is a unity operator up to a sign. Up to the same sign, the Hodge star is its own inverse: and The Hodge star naturally defines an internal product on Ω(X) for φ, ψ ∈ Ω. The internal product is Hermitian positive definite, i.e., Thus, it may serve as a Hermitian metric on Ω. As will be discussed in Section III C 3, the eigensystem of the pseudo-HermitianĤ provides its own non-trivial metric on the Hilbert space. It is this metric that must be viewed as the fundamental metric of the Hilbert space of the model and for which the standard notation ·|· must be reserved, whereas the round brackets can be used for Equation (54).
Equation (54) can also be used for the definition of the concept of the Hermitian conjugate of an operator for any φ, ψ ∈ Ω and any operatorÂ : Ω(X) → Ω(X). Using this definition, it is straight- and the explicit expression for the so-called codifferential, which is the Hermitian conjugate of the exterior derivative, Finally, in the forthcoming discussion, the concept of the Hodge Laplacian will also be recalled:∆

D. Fermionic Variables
The exterior algebra has an alternative field-theoretic representation in terms of the fermionic variables that will be used in the path integral representation of the theory in Section IV A as well as at the end of the next section.
Following reference [49], one notes that the (anti-)commutation relations in Equation (39) are equivalent to those of Grassmann or anticommuting variables, i.e., χ i , and derivatives over them, i.e., ∂/∂χ j : Therefore, one can make the formal substitution and a wavefunction in the new notations becomes whereas the expression for the exterior derivative, reveals whyd is a "super" operator: it destroys a bosonic or commuting variable x i and creates a fermionic or anticommuting variable χ i .
Equation (61) can be viewed as a k-th term of the Taylor expansion of a wavefunction which is now a function of a pair of variables that are the supersymmetric partners with respect to the Operator (62).
Some properties of fermionic variables are similar to those of bosonic variables. For example, one can introduce the fermionic δ-function Here, f (χ) is an arbitrary function of a fermionic variable, χ is yet another fermionic variable, the differential is d D χ = dχ D ...dχ 1 , and Note that the definition of the fermionic δ-function depends on the relative position of has been used.
The above property of fermionic variables and their δ-function can be established using Berezin rules of integration over Grassmann numbers. The latter include identities such as Another property of fermionic variables that has a straightforward bosonic analogue is the exponential representation of a fermionic delta function whereχ is yet another additional fermionic variable.
Other properties of fermionic variables may be in a sense opposite to their bosonic counterparts, e.g., whereas for bosonic variables, one would have d D xδ D (Âx) = |detÂ| −1 for x ∈ R D . There are many other interesting properties and relations associated with fermionic variables (see, e.g., [86]). In the forthcoming discussion, however, only those introduced so far will be used.

A. Stochastic Generalization of Dynamics
In the previous section, the noise configuration was assumed to be fixed, and the dynamics was essentially deterministic. The next step is to account for all possible realizations of the noise. This goal can be achieved as follows.
The stochastic generalization of Equation (13) is and that of Equation (26) is where the notation for the stochastic average is from Equation (4) and the new operator The operation of the stochastic averaging here is legitimate because the pullbacks are linear operators on Ω(X). Now, the possibly highly nonlinear (stochastic) dynamics is described by linear operators acting on the linear Hilbert space. The price one pays for this "linearization" is the infinitely larger dimensionality of Ω(X) compared to the dimensionality of X. One may wonder now if this dramatic increase in the dimensionality of the objects of interest could be an unnecessary complication. It is not. As long as one is interested in the stochastic dynamics, he has to consider the infinite-dimensional space of the probability distributions regardless.
The finite-time stochastic evolution operator (SEO) in Equation (70) is known in DS theory as the generalized transfer operator [78]. The only new element in Equation (70) is that the pullbacks in Equation (70) are those of the inverse maps.
In the case of white (not necessarily Gaussian) noise, the noise variables at different times do not correlate, and thus, where the composition law for pullbacks in Equation (22) has been used. Unlike in Equation The quantum mechanical analogue of the finite-time SEO is the finite-time quantum evolution operator, denoted typically asÛ = e −itĤq , withĤ q being some Hermitian Hamiltonian. The capitalized U signifies here that the quantum evolution is unitary. The stochastic evolution is not unitary; thus, the current notation for the finite-time SEO borrowed from reference [78] is more suitable.
One can now use the picture of the time intervals as discussed in Section I D. Explicitly, the time domain of the temporal evolution, i.e., (t, t ), is a union of a large number, i.e., , where t n = t + n∆t and ∆t = (t − t )/N . The law of the infinitesimal stochastic evolution can now be found as where the following equality has been used: If one now introduces the "infinitesimal" SEO aŝ one obtains The integration of this differential equation with the conditionM tt =1 Ω(X) leads tô One can now introduce the time-dependent wavefunction ψ(t) =M tt ψ(t ) by analogy with the Schrödinger representation of quantum theory. Equation (75) To establish the explicit expression for the SEO in Equation (74) and thatM as follows from the definition of the Lie derivative in Equation (28). Using the linearity of the Lie derivative in its vector field, Equations (7), (8), (74) and (79), one arrives at The physical meaning of the two terms in the SEO (81) is clear. The first term is the deterministic flow along F , and the second term represents the noise-induced diffusion. The with g 0,1 being some functions of its arguments, can be called the diffusion Laplacian. This operator is a member of the family of Laplace operators. In the general case, however, this operator is neither Hodge (or de Rahm) Laplacian (58) Using Equation (35), the above equation results in Now, using the continuous-time analogue of Equation (80) together with ξ a (t) Ns = 0 and Equation (5), one readily arrives at where The subtle point here is that the upper limit of the integration over τ 2 is exactly at the peak of the δ-function, i.e., τ 2 = τ 1 . The δ-function, no matter how narrow, is a symmetric function of its argument, i.e., δ(τ 1 − τ 2 ) = δ(τ 2 − τ 1 ), and its integral over its entire domain is unity. Therefore, τ 1 t δ(τ 1 − τ 2 )dτ 2 must be interpreted as 1/2. Consequently, C = 1, and the SEO in Equation ( 86) is the same operator obtained earlier in Equation (81) within the piece-wise constant picture of the noise.
The conventional Fokker-Planck (FP) equation for the TPD is simply Equation (77) for the D-forms. Its explicit expression can be readily found. Using the Cartan Formula (43) and noting thatdψ (D) = 0 for any This is the well-known FP equation in the so-called Stratonovich interpretation of SDEs.
This brings the discussion to the Ito-Stratonovich dilemma addressed next.

B. Ito-Stratonovich Dilemma
The Ito-Stratonovich dilemma [88][89][90][91][92] is a well-known ambiguity in the exact form of the FP operator that appears when an SDE is looked upon as a continuous time limit of a related stochastic difference equation. A discussion of this issue can begin by rewriting the original SDE (1) in an equivalent form where W a (t) = dξ a (t)/dt is introduced to emphasize that ξ is "integrable" in a certain mathematical sense so that the SDE is well defined.
The classical view on SDEs is through the continuous-time limit of the related stochastic difference equations (SdE): Here, ∆x n = x n − x n−1 , with x n and x n−1 being the DS variables at two consecutive time moments; t n = t n−1 + ∆t; and ∆W a n = ξ a n ∆t, with ξ a n being the noise variable acting between t n and t n−1 ; moreover, the parameter α ∈ i.e., α = 1/2, of the mid-point (x n + x n−1 )/2; and the Kolmogorov or "isothermal" choice, i.e., α = 1 [87], of the final point x n .

As shown in Appendix A 1, the FP equation of the continuous-time limit of the SdE is
Equation (88) with the shifted flow vector field: It is also shown in Appendix A 1 that if, instead of observing the SdE as a formal equation This equation is just as well defined as the original SDE (89). Indeed, the transition from Equation (89) to (92) uses only the two following fundamental properties: diffeomorphisms defined asẋ = F induce the (time-reversed) action on Ω defined by ∂ t ψ = −L F ψ, and the Lie derivative is linear in its argument,L C 1 F 1 +C 2 F 2 = C 1LF 1 + C 2LF 2 , where F 1,2 are some (differentiable) vector fields on X and C 1,2 are constants. Therefore, if the SDE together with dW a is well defined, then the SFE is also well defined.
One way to understand why the FP operator for Equation (92) does not depend on the interpretation of the noise is as follows. As is clear from Equation (91), the FP operator for Equation (89) is independent of the interpretation of the SDE when the noise is additive, i.e., when the e are independent of the position of X so that the factor (2Θ) 1/2 e a , to which the noise is coupled in Equation (89), is independent of time. Similarly, the factor (2Θ) 1/2L ea , to which the noise is coupled in Equation (92), is always independent of time even for the position-dependent e, which renders the FP operator derived from Equation (84) noiseinterpretation independent.
Within the classical view on SDEs as a continuous-time limit of SdEs, no interpretation can in principle have a qualitative mathematical advantage over the others. This follows from the mere fact that different interpretations can be transformed among themselves by a mere shift of the flow vector field in accordance with Equation (91). Indeed, consider, for example, an Ito SDE and a Stratonovich SDE with the flow appropriately shifted such that the FP operators of the two models are the same. These two models define the same stochastic model, and there is simply no room to accommodate any mathematical advantage of one model over the other.
In other words, the classical theory of SDEs has an intrinsic redundancy in the sense that each stochastic model has infinitely many representatives corresponding to different interpretations of the SDE. This redundancy can only be removed if there existed a reference point outside the classical view of SDEs as a continuous-time limit of α-parameterized SdEs.
The SFEs introduced above may serve as this external reference point. If this point of view is adopted, it must be said that the Ito-Stratonovich dilemma is resolved in favor of the Stratonovich interpretation of SDEs because Stratonovich SDEs provide the same FP operator as the corresponding SFEs. Within this framework, one would state that the continuous-time limit of two SdEs with the same F and e but with different α is two different SDEs with shifted flow vector fields. This point of view is summarized in Figure 4.
Supporting this proposition, it can be mentioned that the Stratonovich approach is more natural from a mathematical point of view (see, e.g., [93] and the references therein). On the side of physics, the Stratonovich interpretation is superior because all noise sources are never white; rather, they have finite correlation times, and in the white-noise limit, the colorednoise SDEs become Stratonovich SDEs [91]. Furthermore, it is experimentally established that the Stratonovich approach is more accurate in regard to the numerical simulations of physical models [94]. Simultaneously, the only known advantage of the Ito interpretation of being "respective" of the Markovian property is a misinterpretation of the fact that the SdEs with α = 0 define the increment as a function of the initial point implicitly, as discussed at the end of Appendix A 1.
The resolution of the Ito-Stratonovich dilemma, however, is a purely mathematical problem, and the above picture should only be viewed as one of the scenarios of its possible solution. Fortunately, it is not important for further discussion whether the Ito-Stratonovich dilemma is resolved at this point. To liberate the STS from the burden of the Ito-Stratonovich dilemma, one can always think that the STS is a theory of SFEs and not of SDEs. If one is interested in the STS of the Ito or any other interpretation of the SDE,  all he has to do is to shift the flow vector field accordingly.
The Ito-Stratonovich dilemma will be readdressed in Section IV D, where it will be shown that the Stratonovich interpretation of SDEs is equivalent to the Weyl symmetrization rule, whereas the so-called Martingale property of the Ito interpretation is equivalent to the unphysical convention of placing all the momentum operators after all the position operators.
It will also be reestablished that different interpretations of SDEs can be transformed among themselves by a shift of the flow vector field and that this result is correct for the entire SEO and not only for the FP operator as discussed here.

C. Properties of the Stochastic Evolution Operator
In this subsection, some of the most important properties of the eigensystem of the SEO are discussed. For simplicity, it is assumed that the phase space is closed, the noise-induced metric g ij is positive definite everywhere on X, Θ > 0 so that the SEO is elliptic, andĤ is diagonalizable, with a discrete spectrum bounded from below. It is natural to believe that most of the claims here hold true or are at least transformative to more general classes of models.

Fermion Number Conservation
The SEO is of "zeroth" degree, i.e., degĤ = 0 (see the definition of the degree of an operator in Equation (40)). This means that it does not mix wavefunctions of different degrees. As a result, the operator of the degree of the differential form in Equation (47) is commutative with the SEO: As noted in Section II D, the differentials of wavefunctions can be viewed as fermions.
Therefore,k can be interpreted as the number of fermions, and Equation (93) reflects the conservation of this quantity.
Yet another method of expressing the idea of the conservation of the number of fermions is to note the block-diagonal structure ofĤ: is the projection ofĤ onto Ω k (X).

Completeness
The SEO is a real operator. Therefore, its spectrum consists of real eigenvalues and pairs of complex conjugate eigenvalues that are called Ruelle-Pollicott resonances in DS theory.
This form of the spectrum is a sufficient condition forĤ to be pseudo-Hermitian [60]. As a pseudo-Hermitian operator,Ĥ has a complete bi-orthogonal eigensystem where ψ andψ are the right and left eigenfunctions (or rather eigenforms), and the bra-ket notation is |ψ n ≡ ψ n ∈ Ω kn (X) and ψ n | ≡ψ n ∈ Ω D−kn (X).
Here, k n is the eigenvalue of the operator of the degree of a wavefunction: As noted in the previous subsection, the degree of a wavefunction is a "good quantum number" becausek is commutative with H and consequently the both operators can be simultaneously diagonalized.
The eigenstates of a given degree, say, k, provide a complete bi-orthogonal eigensystem on Ω k (X) such that the resolution of unity on Ω k (X) is The bra-ket combination of an eigenstate of any degree is a D-form, which has the meaning of the TPD associated with this eigenstate.

Pseudo-Time-Reversal Symmetry
The pseudo-Hermitianity of the SEO is closely related to the pseudo-time-reversal symmetry of the model. To begin the discussion of this symmetry, one needs the explicit expression forĤ † that can be established as follows. The defining property ofĤ † is given by Equation (55): One notes now that the Lie derivative is a differentiation, i.e., for any ψ 1,2 ∈ Ω(X) and any vector field G ∈ T X. In addition, becausedψ (D) = 0 for any ψ (D) ∈ Ω D (X). The last two formulas lead to the conclusion that, for any ψ 1 ∧ ψ 2 ∈ Ω D , This equality can now be used to rewrite the Left Hand Side (L.H.S.) of Equation (101) as whereĤ Here, −1 is the inverse of the Hodge star defined in Section II C 3, and the following has been used: because bothĤ T and are real. Equations (101) and (105) givê The next goal is to examine the properties of the model with respect to the time-reversal operation. One first notes that Equation (106) is the SEO of the SDE obtained from Equation (1) by reversing the flow of time:ẋ In other words, the "naive"reversal of time has the following effect on the SEO: in the original stochastic evolution in Equation (77).
Let us recall now that the time reversal in quantum mechanics is also accompanied by swapping bras and kets. This is needed because the information of the system's past and future are stored in the kets and bras, respectively, whereas the past and future are interchanged by the time reversal. This bra-ket swapping in the STS is accomplished using the Hodge star where T = T denotes the composition of the operations.
In Section V, it will be discussed that the structure of the wavefunction of the supersymmetric ground states is such that the coordinate directions with/without differentials correspond to the stable/unstable (local) variables. On the other hand, the time reversal makes the stable variables unstable and vice versa. Therefore, this operation must act on a wavefunction in such a way that, in the directions with/without differentials, the wavefunction loses/acquires differentials. This understanding strengthens the above relation between the time-reversal operation and the Hodge star operation, which acts on a wavefunction in exactly this manner.
Clearly, an operation T does not seem to be a symmetry of the model becauseĤ =Ĥ † for pseudo-HermitianĤ. Nevertheless, the stochastic evolution defined byĤ † may turn out to be physically equivalent to that defined byĤ. On a mathematical level, this "equivalency" means that there exist such that The existence of such η is the definitive property of all pseudo-Hermitian operators. Therefore, the model does possess the so-called ηT -symmetry: The operator η can be called the Hilbert space metric. It relates the bras and the Hodge duals of the eigenstates as where the notation (ψ k | ≡ ψ * k was previously introduced in Equation (53). From the orthogonality property in Equation (96b), one finds In other words, η is the inverse of the "overlap" matrix, i.e., of (the transpose of) the matrix of the inner products of the kets of the eigenstates.
The ηT operation acts on a wavefunction as If ψ T,n is an eigenstate of the time-reversed SDE (109), then ηT (ψ T,n ) is an eigenstate of the original SDE (1) although with a complex conju- It can be said that the eigenstates with complex eigenvalues (the Ruelle-Pollicott resonances) break the ηT symmetry. If one such eigenstate is a ground state of the model, the ηT symmetry can be said to be spontaneously broken because the ground state of the time-reversed SDE has a different eigenvalue.
Equation (115) can also be rewritten aŝ That is,Ĥ (k) andĤ are related by the similarity transformation η. This immediately suggests that the two operators are isospectral: This result will be used later in the discussion of the possible forms of the SEO spectra in Section III C 8.
Up to this moment, it was not specified what phase-space metric g is being used. This phase space metric, or X-metric for short, enters the above formulas through the definition of the Hodge dual in Equation (45). What has been said so far in this subsection is correct for any "good enough" X-metric. In other words, one has a freedom in choosing g. This freedom can in principle be used for the simplification of the explicit expressions for g-dependent objects such asĤ † and η. For g-independent objects, such asĤ,Ĥ T , the Eigensystem (96), and the composition η, the choice of g is unimportant.

Topological Supersymmetry
The SEO isd-exact, i.e., it has the form of a bi-graded commutator The exterior derivative is commutative with the SEO: This can be observed from the nilpotency property of the exterior derivative, i.e.,d 2 = 0, leading to the conclusion thatd commutes with anyd-exact operator: The commutativity of an operator with the SEO indicates that this operator is a symmetry of the model. The reason whyd is a symmetry can be explained as follows. The finite-time SEO is a stochastically averaged pullback induced by the SDE-defined diffeomorphisms. Therefore, the finite-time SEO commutes withd because any pullback induced by a diffeomorphism is commutative withd. In other words, this symmetry is a consequence of the fact that continuous(-time) dynamics preserves the continuity of the phase space, as previously mentioned in Section I C.
Note also that not all possible evolution operators that commute withd are necessarilŷ d-exact as in Equation (123). Ad-exact evolution operator implies more than simply the commutativity withd. As will be discussed below, the additional implication of ad-exact evolution operator is that all thed-symmetric eigenstates have a zero eigenvalue.
In terms of the fermionic variables of Section II D, the exterior derivative substitutes commuting or bosonic variables with anticommuting or fermionic ones. Therefore, it can be identified as a supersymmetry.

Topological Supersymmetry vs. N = 2 Supersymmetry
The topological supersymmetry operator is the same for all SDEs. The operator does not contain any information on the specifics of dynamics, which in turn are solely encoded in the other fermionic operator, i.e.,d, in Equation (124). One way to look at this operator is as the operator of the current of the probability density. This point of view is at least partially correct, as is evident from the stochastic evolution equation for top differential forms can the model be said to be N = 2 supersymmetric. In particular, the evolution operator is a square of N = 2 mixed-degree fermionic operators, In the general case, however, Equation (129) is incorrect, and the model can only be identified as a member of the family of the cohomological field theories for which the evolution operator is a bi-graded commutator with the topological supersymmetry operator, as in Equation (123). In other words, the difference between general form SDEs and Langevin SDEs and/or Hamilton models is the same as the difference between topological quantum mechanics (see, e.g., [96]) and N = 2 supersymmetric quantum mechanics (see, e.g., [97]).
In addition, note that, to identify a model as a cohomological field theory, one must also require that the ground states of the model be supersymmetric. In the STS, however, the ground states are not supersymmetric in the most interesting situations with the spontaneously broken topological supersymmetry. Therefore, the identification of the STS as a cohomological field theory is technically inaccurate.
It is well known that N = 2 supersymmetry leads to the pairing of the non-supersymmetric eigenstates into boson-fermion doublets, whereas all the supersymmetric eigenstates are singlets with exactly zero eigenvalues. In the next Subsection III C 6, it will be shown that the topological supersymmetry tailors the same structure of the eigensystem.

Boson-Fermion Pairing of Eigenstates
From the group-theoretic point of view, the topological supersymmetry is a continuous one-parameter group of transformationŝ of which the SEO is invariant:Ĝ sĤĜ−s =Ĥ.
As in the case of any other symmetry, the eigenstates must be irreducible representations of this group. There are only two types of irreducible representations of this symmetry: most of the eigenstates are non-d-symmetric "bosonic-fermionic" doublets or pairs of eigenstates, and some of the eigenstates ared-symmetric singlets.
Each pair of non-d-symmetric eigenstates, which will be denoted as |ϑ and |ϑ , can be defined via a single bra-ket pair, i.e., |ϑ n and ϑ n |, such that so the bra-ket pairs of the non-d-symmetric pairs of eigenstates can be given as and |ϑ n =d|ϑ n , ϑ n | = ϑ n |.
Here, the expression ϑ n |d must be understood as a differential formθ n such that Xθ n ∧ dγ = Xθ n ∧ γ for any γ, where the barred notation for the bras in Equation (96b) has been used. Using the standard relation d (γ 1 ∧ γ 2 ) = 0, valid for all γ 1,2 , it can be easily established that, up to a sign, the differential formθ n =dθ n .
The orthogonality relations for the ϑ are These relations can be derived from the structure of the non-d-symmetric eigenstates in Equations (133) and (134) as follows. Consider an eigen-bra ϑ n | ∈ Ω D−k . There are two types of eigen-kets of degree k that can potentially overlap with ϑ n | non-trivially: |ϑ n 1 ∈ Ω k and thed-exact eigen-ketsd|ϑ n 2 ∈ Ω k with |ϑ k ∈ Ω k−1 . The eigensystem is bi-orthogonal; thus, only one eigen-ket can provide a non-zero overlap with ϑ n |. As observed from Equation (132), this eigen-ket isd-exact. All the other kets must have zero overlap with ϑ n | because they correspond to different eigenvalues of the SEO, which has no spectral degeneracy in the most general situation. As long as this argument is valid in the most general case of no spectral degeneracy, it is always valid even in the case where there is spectral degeneracy due to some additional symmetry of the model or accidentally.
The pairing of the non-supersymmetric eigenstates into the boson-fermion doublets can be demonstrated as follows. Consider the complete set of the eigenstates ofĤ (0) , i.e., |ϑ n 0 ∈ Ω 0 , where n 0 is the label running over all these eigenstates. Consider also the set ofd-exact states, i.e., |ϑ n 0 =d|ϑ n 0 ∈ Ω 1 . These are the eigenstates ofĤ (1) with the same eigenvalues as |ϑ n 0 becaused commutes withĤ. This set is incomplete in Ω 1 , and there exists another set of eigenstates ofĤ (1) , i.e., |ϑ n 1 ∈ Ω 1 , with n 1 being yet another label running over this new set of eigenstates. Unlike |ϑ n 0 , |ϑ n 1 are notd-closed, i.e., the operatord does not annihilate them in general. Once again, one considers the set of thê d-exact eigenstates ofĤ (2) , |ϑ n 1 =d|ϑ n 1 ∈ Ω 2 . This recurrent procedure terminates at Ω D , the eigenstates of which are alld-closed. At the end of this procedure, almost all (see below) the eigenstates are non-supersymmetric pairs related byd.
A similar procedure exists for bras. The only difference is that the procedure runs in the opposite direction, i.e., from ϑ n D | to ϑ n 0 |, becauseθ n k ∈ Ω D−k , as observed from Equation (96b), and consequently ϑ n k |d is an eigen-bra ofĤ (k−1) and not ofĤ (k+1) , as could be incorrectly expected.
The boson-fermion pairing procedure does not count eigenstates that are non-trivial in the de Rahm cohomology, in other words, the eigenstates that ared-closed but that are not These eigenstates are the supersymmetric ord-symmetric singlets. Their bras are also non-trivial in the de Rahm cohomology: An important property of the supersymmetric eigenstates is that the expectation value of anyd-exact operator vanishes on these eigenstates: Note now that the SEO is ad-exact operator. Therefore, all thed-symmetric eigenstates have zero eigenvalues: Each de Rahm cohomology class must provide oned-symmetric eigenstate of the form where |h k is the harmonic differential form from this particular class of the de Rahm cohomology. That this is true can be demonstrated using perturbation theory, as done in Appendix A 2. The statement that each de Rahm cohomology class provides oned-symmetric eigenstate must be correct even outside the domain of the applicability of the perturbation theory. This follows from the "completeness argument": if it is not so, then the eigensystem of the pseudo-Hermitian SEO is incomplete, which contradicts the theory of pseudo-Hermitian operators.

Topological Supersymmetry and N = 2 Pseudo-Supersymmetry
In Section III C 5, it was discussed that, in the general case,d is not a supersymmetry of the model. Nevertheless, the second supercharge does exist. Although its explicit form may not be easy to establish, the second supercharge can be easily constructed out of the eigensystem of the model. Using the notations in Equations (133) and (134), the SEO can be given asĤ where the second supercharge of the model iŝ The two operatorsd ‡ andd must differ by ad-closed piece, i.e., because Equations (123) and (140) where the notations from Equations (133) and (134) have been used. This operator can be visualized in Figure 5 as reverse arrows representingd. In other words, the second supercharge does not contain any additional information on the supersymmetric structure of the SEO, and for this reason, it will not be considered in this paper any further. It must be noted, however, that in situations in which there exists an operator η such thatd ‡ = η −1d † η In Section IV F, it is discussed that it is possible to view only one of these eigenstates as the ground state. Thin dotted arrows represent the action ofd that couples all the non-d-symmetric eigenstates into boson-fermion pairs.

Thermodynamic Equilibrium and Stochastic Poincaré-Bendixson Theorem
In Section III C 6 above, it was argued that each de Rahm cohomology class must provide one supersymmetric eigenstate. Whether this is true is not important for further discussion. What is important is the existence of the supersymmetric state of thermodynamic equilibrium (TE), i.e., the steady-state (zero-eigenvalue) TPD: ψ T E ∈ Ω (D) (in DS theory the TE state is known as the invariant measure). Its presence can be established through the physical version of the completeness argument. Indeed, all the non-d-symmetric eigen-states from Ω D (X) ared-exact, i.e., they are of Type (134). This means that the integral of all such eigenstates over X is zero: Xd ϑ = 0. On the other hand, a wavefunction from Ω D (X) has the meaning of the TPD. The integral of a meaningful TPD over X must not vanish. This suggests that at least oned-symmetric eigenstate from Ω D (X) must exist for physical models. It can also be shown that the bra of the TE state is a constant function: It is also an important piece of understanding that the TE state is always the "ground state" for Ω D . In other words, among all the eigenstates from Ω D , it has the smallest real part of its eigenvalue, which is zero, of course. Indeed, imagine that this is not true and that there exists an eigenstate from Ω D (X) such that its eigenvalue is real and negative.
As mentioned in the previous paragraph, the integral of the ket of such a non-d-symmetric eigenstate over X vanishes. Thus, the ket of this eigenstate must be negative somewhere on X. One can now take a random TPD and evolve it in time sufficiently long. This will eventually lead to the situation whereby the TPD will become negative somewhere on X due to the dominant contribution from this presumably existing non-d-symmetric eigenstate.
The negative TPD is not physical. Thus, it can be concluded that, for physical models, this situation is not realizable, and there are no eigenstates in Ω D with real negative eigenvalues.
In a similar manner, one can rule out the possibility that there is a pair of Ruelle-Pollicott resonances in Ω D with a negative real part of their eigenvalues. Thus, for physical models, the TE state is the "ground state" in Ω D .
The same reasoning applies to the SEO (106) of the time-reversed SDE (109). Combined with the fact that specĤ (0) = specĤ (D) T , which follows from Equation (122), this observation suggests thatĤ (0) also never breaks the supersymmetry, i.e., the supersymmetric zeroeigenvalue eigenstate is the ground state in Ω 0 . The ket of this ground state of Ω 0 is a constant function on X.
This brings the discussion to the Poincaré-Bendixson theorem stating that smooth deterministic flows can be chaotic only in three-plus dimensions (see, e.g., [95]). The above analysis of the SEO spectra seemingly leads to the stochastic version of this theorem. Indeed, as long asĤ (D) andĤ (0) do not break the topological supersymmetry as discussed in the two previous paragraphs, the overall supersymmetry cannot be spontaneously broken unless the dimensionality of the phase space is three or higher. This can be straightforwardly deduced from Figure 5. Clearly, if the dimensionality of the phase space is less than three, at least one eigenstate of a pair of the non-supersymmetric eigenstates (the degrees of which differ by one) with a negative real part of their eigenvalue must be either in Ω D or Ω 0 , which contradicts the above properties of the spectra ofĤ (D) andĤ (0) .

Realizable Spectra
The properties of the SEO discussed previously limit its possible spectra to only the three types given in Figure 5. Furthermore, in DS theory, there are theorems stating that, for a certain class of models that mimic chaotic behavior, namely, the so-called expanding models, the eigenvalues of the ground states must be real [78].
There are no other limitations. To convince oneself, one needs at least one example for each of the two types of SEO spectra with spontaneously broken supersymmetry (Figure 5b,c). These examples must not necessarily be analytical. Well-established numerical examples are sufficient. Such examples exist in the theory of the magnetohydrodynamical phenomenon of kinetic dynamo (KD), as was very recently found in [98] and discussed briefly in Appendix A 3. Thus, both types of the supersymmetry breaking spectra in Figure 5 are realizable.

D. Witten Index
One of the fundamental partition-function-like objects is the "sharp trace" of the finitetime SEO known in supersymmetric quantum theory as the Witten index: It was previously established that all the eigenstates with non-zero eigenvalues are nond-symmetric. They come in pairs of even and odd degrees. Thus, their contributions cancel out from the Witten index. Onlyd-symmetric eigenstates with zero eigenvalue contribute to the Witten index, which is thus independent of the duration of the time evolution: where b k is the number ofd-symmetric states of degree k. If one believes that each de Rahm cohomology class provides oned-symmetric state, then b k are the Betti numbers, and W equals the Euler characteristic of X, Eu(X).
The next goal is to discuss the physical meaning of Equations (144) and (155)  This can be performed using the fermionic variables in Section II D. The pullback in Equation (27) can be given as where the bosonic and fermionic δ-functions substitute the arguments x and χ by M t t (x) and T M t t (x)χ, respectively, with T M tt (x) being the tangent map in Equation (16).
Following stochastic averaging, one arrives at the finite-time SEO in the following repre- The Witten index in Equation (144) takes the form of the trace of the finite-time SEO with periodic boundary conditions for both the commuting and anticommuting variables: where the characteristic polynomial formula det(1 + λT M t t (x)) = D k=0 λ k m k (x) has been utilized, with The denominator in Equation (149) originates from the integration over the bosonic variables, whereas m k can be viewed as a fermionic trace over Ω k (X). Indeed, the basis of the differentials in Ω k (X), which the fermionic variables represent, is given by the C k D ordered combinations of the differentials: dx i 1 ∧ ... ∧ dx i k , i 1 < ... < i k . Thus, the trace of the fermionic variables over Ω k (X) is where d(M tt (x)) is defined in Equation (15). One arrives now at which is yet another version of Equation (144). (149) can be given as

The Witten index in Equation
where is known as the Lefschetz index of the map M t t . The Lefschetz-Hopf theorem states that, under some general conditions, where the trace is over the de Rahm cohomology H k (X). In the limit t → t, when M tt → Id X , the Lefschetz index reduces to the signed sum of the Betti numbers, i.e., to the Euler characteristics of X. On the other hand, it was previously established that W is independent of t. This leads to the conclusion that the Witten index equals the Euler characteristic of X for any duration of temporal evolution.
The topological character of W can be qualitatively understood in the following manner (see Figure 6).  and such a situation is clearly somewhat suspicious because, as previously mentioned, the partition function of the noise is one of the fundamental objects of the model, i.e., the object that appears at the level of the very formulation of stochastic dynamics.

E. Dynamical Partition Function
Yet another fundamental object is the dynamical partition function (DPF) and/or the "counting trace' of the finite-time SEO: By analogy with Equation (149), the DPF is the trace of the finite-time SEO with periodic/anti-periodic boundary conditions for the bosonic/fermionic variables: where the m are defined in Equation (150).
The physical meaning of the DPF is observed in the limit of the infinitely long temporal evolution. Consider models in which the absolute values of the eigenvalues of the tangent map Spec(T M t t ) = (µ 1 (t, t ), ..., µ D (t, t )), in the long-time limit t − t → ∞, are such that , with λ being the stochastic versions of the (global) Lyapunov exponents. The class of models that satisfy this condition must certainly exist in the deterministic limit as follows from the classical DS theory.
In the assumption that none of the λ vanish, one has, in the limit of t − t → ∞, In other words, in this class of models, the DPF grows slower than the stochastically averaged number of fixed points of the SDE-induced diffeomorphisms or, equivalently, of the number of periodic solutions of the SDE.
At this point, it must be stressed that, in this and previous subsections, the summation over the fixed points of the SDE-defined diffeomorphisms (see, e.g., Equation (149)) has been used as if these fixed points were isolated in X. This is not true in general. The fixed points of the diffeomorphisms may appear in submanifolds of X, with the Morse-Bott flow vector fields being one example of this situation. How to count the fixed points in this general situation using the standard methodology of DS theory used in this subsection is not clear. This problem, however, does not exist in the operator representation of the theory considered previously as well as in its path integral representation in the next section.
It is also important to discuss the fundamental difference between the dynamical partition function of the STS and the thermodynamic partition function in statistical (quantum) physics. The latter is defined as T re −βĤq , whereĤ q is a Hermitian Hamiltonian of a quantum model and β is the inverse temperature. Equation (155) has a very similar appearance.
Furthermore, in the literature on, e.g., N = 2 supersymmetric quantum mechanics, it is often said Equation (155) is the result of the Wick rotation of the "real" time of evolution, i.e., t → i × t ∼ β. This is not so from the point of view of the STS. The time t in Equation (155) is the original time of the stochastic evolution and not that of the Schrödinger evolution.
This explains the absence of the imaginary unity in the exponent. The direct quantum analogue of Equation (155) is the generating functional T re −itĤq .

F. Topological Supersymmetry Breaking, Chaos, and Dynamical Entropy
In models with the type of SEO spectra given in Figure 5b, the DPF grows exponentially in the long-time limit: where the factor of 2 comes from thed-degeneracy of the non-d-symmetric ground state and E g is the ground state's eigenvalue, i.e., the eigenvalue with the smallest and negative real part. For the type of spectra in Figure 5c, one has where the pair of the Ruelle-Pollicott resonances with the least real part of their eigenvalues provides the dominant contribution in the long time limit.
It can be recalled that, in deterministic chaotic DSs, the number of periodic solutions grows exponentially with time in the long time limit. The rate of this exponential growth is related to the concept of dynamical entropy (see, e.g., [95]). This exponential growth is provided by the infinite number of unstable periodic orbits with arbitrary large periods, which constitute strange or fractal attractors [76]. This exponential growth is basically the reason why chaotic dynamics is sometimes identified as complex dynamics. This term is borrowed from information theory. There, a problem is identified as complex if the number of elementary operations needed to obtain its solution grows exponentially with the "size" of the problem.
As demonstrated in Section III E, for a wide class of models, the number of periodic solutions grows faster than the DPF. For spectra given in Figure 5b,c, the stochastically averaged number of periodic solution grows at least exponentially and the supersymmetry is spontaneously broken because the ground states are non-d-symmetric as they have nonzero eigenvalues. Thus, one concludes that the stochastic generalization of the concept of deterministic chaos is the spontaneous breakdown of the topological supersymmetry. Even more convincing evidence that this is indeed so will be provided in Section IV G, where it will be shown that the spontaneous breakdown of topological supersymmetry must always be accompanied by the emergence of the long-term memory of perturbations that must be associated with the famous butterfly effect.
In DS theory, there exists the so-called Shub conjecture (see, e.g., [109] and the references therein), stating that, for a sufficiently smooth map M : X → X, the spectral radius of M * : H * (X) → H * (X), where H * denotes the homology group, provides a lower bound for the topological entropy, i.e., the central measure of chaos (see, e.g., [110,111] and the references therein). The spectral radius of the finite-time stochastic evolution operator M * 0t Ns = e −tĤ is (up to a sign) the real part of its ground-state eigenvalue, which can therefore be recognized as the stochastic generalization of the lower bound for the topological entropy in the Shub conjecture. Thus, if the real part of the ground-state eigenvalue is negative, the topological entropy is positive, and the model must be identified as chaotic.
As a conclusion of this section, it can be stressed that the stochastic chaos is the opposite or rather complementary concept of that of the thermodynamic equilibrium, i.e., to the situation with unbroken supersymmetry when the supersymmetric state of the thermodynamic equilibrium discussed in Section III C 8 is among the ground states of the model, as in Figure 5a. Indeed, the supersymmetry can be spontaneously broken or unbroken but not both at the same time.

IV. PATH INTEGRAL REPRESENTATION
Path integrals are a powerful analytical tool that can greatly simplify various tasks that otherwise would be tedious. They represent a component of the mathematical foundation of quantum theory and have also been used in the studies of stochastic dynamics (see, e.g., [112] and the references therein). In the case of the STS, the path integral representation of the theory addressed in this section (see Figure 7) allows particularly for the generalization of the theory to models with noise of any form, not simply Gaussian white noise.

FIG. 7. Path integral representation of the finite-time stochastic evolution operator (SEO). Each
time slice t n hosts a boson-fermion pair of variables x(t n ) ∈ X and χ(t n ) ∈ T X x (t n ). In between the time slices, there are pairs of Lagrange multipliers and fermionic momenta from the cotangent space B(t n ),χ(t n ) ∈ T X * x(tn) or T X * x(t n−1 ) (in the continuous time limit, this choice makes no difference), as well as the noise variables ξ(t n ). The time flows from right to left, as explained at the end of Section II C 1. The finite-time SEO is obtained by integrating out all the variables except x(t), χ(t) and x(t ), χ(t ). Further integration over x(t), χ(t) with the periodic boundary conditions x(t) = x(t ), χ(t) = χ(t ) results in the Witten index W , whereas the integration with the anti-periodic boundary conditions for the fermionic variables x(t) = x(t ), χ(t) = −χ(t ) results in the dynamical partition function Z.

A. Finite-Time Stochastic Evolution Operator
In the discrete-time picture introduced in Section III A, the domain of the temporal evolution is split into N 1 segments with boundaries at t n = t + n∆t, ∆t = (t − t )/N , t N ≡ t and t 0 ≡ t . The finite-time SEO can be given as the stochastically averaged composition of pullbacks at each time segment: where the composition law for the pullbacks from Equation (22) has been used multiple times.
At each time moment, one can now introduce a copy of the phase space and a pair of bosonic and fermionic variables x(t n ), χ(t n ) so that Equation (159) can be given aŝ where the pullbacks on the R.H.S. are combinations of the bosonic and fermionic δ-functions, as in Equation (147).
To exponentiate the bosonic δ-function, one introduces an additional bosonic variable called the Lagrange multiplier or bosonic momentum from the cotangent space of X B(t n ) ∈ T X * x(tn) : In a similar manner, one can exponentiate the fermionic δ-function using Identity (66) and introducing the fermionic momentum from the cotangent spaceχ ∈ T X * x : δ D (χ(t n−1 ) −M t n−1 tn χ(t n )) = d D (iχ(t n ))e iχ i (tn)(χ(t n−1 )−Mt n−1 tn (x(tn))χ(tn)) i .
Here, the imaginary unity is needed to bring the model to the form conventional in the literature on cohomological field theories.
In the continuous-time limit, i.e., N → ∞, ∆t → 0, one has where T F(x) is introduced in Equation (30). Combining the above representation of the δ-functions and Equation (162), one arrives at Here,S where The action can be expressed in the so-called Q-exact form whereΨ is known as the gauge fermion and the curly brackets denote the operator of the topological supersymmetry: This operator is the path integral version of the exterior derivative. In particular, it has similar properties: it is nilpotent, i.e., {Q, {Q, X(Φ)}} = 0 for any X(Φ), and it is a bi-graded differentiation, i.e., where X and Y are some functionals of Φ, and deg(X) is the degree of X defined as the difference between the numbers of χs andχs in X. Equation (170) is the path integral version of Equation (42).
To perform the stochastic averaging, one can first separate the noise term in the action with and As a next step, one can integrate out the noise field, which is no longer assumed to be Gaussian white. The noise, however, remains assumed physical so that the path integration over all the noise configurations remains well defined.
Integrating out the noise transforms Equation (163) into the following form: where the arguments of the finite-time SEO are dropped for brevity and the new action with c (k) being the irreducible correlators of the noise. The zeroth-order term in the Taylor series vanishes because the partition function of the noise is assumed normalized, i.e., for y = 0, one has log e 0 Ns = log1 = 0. Now, using the property of the nilpotency of Q, the differentiation rule in Equation (170) and the fact that all y in Equation (173) are Q-exact, one arrives at where the new gauge fermion with Ψ 0 defined in Equation (172) and for k ≥ 1. In other words, even after stochastic averaging, the action is Q-exact.

B. Interpretations of Stochastic Quantization
A Q-exact action as in Equation (174) is a definitive feature of cohomological field theories [44][45][46][47][48]. Their standard path integral representations include periodic boundary conditions (PBC) for the fermionic variables. This is related to the fact that the PBC are consistent with the Q-operator: This equality and Equations (149) and (174) lead to where the path integration is over As a result, the topological character is lost for the DPF Note that this does not mean that the model described by the DPF no longer has the topological supersymmetry. The topological supersymmetry is not a property of the DPF but rather of the most fundamental object in the theory, the finite-time SEO. Equation 179 simply means that this particular object, i.e., the DPF, that we construct from the finite-time SEO is not of topological character.
The fundamental difference between the DPF and the Witten index can be best revealed in the limit of the infinitely long temporal evolution. There (see Section IV E below), only the ground states contribute to the DPF, whereas only thed-symmetric states contribute to W . Thus, the difference between W and Z is particularly pronounced under the conditions of the spontaneously broken topological supersymmetry when the ground states of the model are non-d-symmetric.

C. Generalization to Spatially Extended Models
The class of models under consideration can be generalized further to spatially extended models. These models are defined by the following stochastic (partial) (integro-)differential equations:ẋ where r is the spatial coordinate of the "base-space". In the general case, the flow vector The stochastic quantization procedure of Equation (180) is along the same lines, the action is Q-exact with the topological supersymmetry operator and the gauge fermion before integrating away the noise variables is After integrating out the noise, one arrives at a model with a Q-exact action.

D. Weyl-Stratonovich Symmetrization and Martingale
The story of stochastic quantization would not be complete without a discussion on how the path integral representation can be turned back into the operator representation.
The exercise to be conducted in this subsection will reveal a close relation between the Ito-Stratonovich dilemma (see Appendix A 1 and the end of Section III A) and the Weyl symmetrization rule of quantum theory.
Models with Gaussian white noise are of interest herein. The gauge fermion can be acquired from Equation (174d), and recalling that the only irreducible correlator of the Gaussian white noise is the one given in Equation (5), whered is the path integral version of Operator (124).
Accordingly, the action is where This function is the stochastic analogue of the Hamilton function in the path integral representation of quantum mechanics and can thus be called the stochastic Hamilton function (SHF). The SHF (186) is the path integral version of the SEO (81). In particular, Equations (43), (81) and (186) reveal that Q is the path integral version of the commutator with the exterior derivative, whereas iχ i is that of the interior multiplication and/or the fermionic momentum operator in Equation (60).
Consider now infinitesimal temporal evolution between t n−1 and t n = t n−1 + ∆t. In the continuous-time limit ∆t → 0, the path integral representation of the infinitesimal stochastic evolution is Here, for the sake of brevity, the notations are different from those in Figure 7. The relation with the notations in Figure 7 is as follows: and y, ϕ ≡ x(t n−1 ), χ(t n−1 ). Other notations introduced in Equation (187)  representation that different choices of α lead to different SEOs, as will be observed below.
Using Equation (187), the infinitesimal evolution of the wavefunction can be given as The Taylor expansion of the exponent in the R.H.S. of Equation (187) in ∆t results in The sign in front of the bosonic momentum is unambiguous because this operator acts on the ordinary commuting variables. In contrary, the sign in front of the fermionic momentum operator is ambiguous or rather can be established unambiguously only by a tedious exercise of carefully tracking all the signs associated with the relative positions of the anticommuting variables and their differentials. There is a simpler way, however, of accomplishing this task.
One can demand that the bi-graded commutator of the exterior derivative in the operator representation act in the same way as the operator of the Q-differentiation in the path integral representation. In other words, because which is just another version of Equation (169), It can be readily verified that the choice of sign in Equation (190) satisfies this requirement.
Another important observation from Equation (189) is related to the order of operators.
Both the variables x and y (and χ and ϕ) in the path integral representation correspond to the operatorx (andχ) in the operator representation of the theory. The difference between, say, x and y is that y acts on the wavefunction before the momentum operators B, whereas x acts after B. In particular, if H included a term B i x j α , this term in the operator representation would become The same can be said about fermionic operators, with the only correction being that the symmetrization must be bi-graded, i.e., This reveals that if, in the Ito case (α = 0), all the momentum operators must act after all the position operators, in the Stratonovich case (α = 1/2), the operators must be sym- This is indeed true, as will be demonstrated next. To facilitate the procedure of establishing the expression for the SEO, one can utilize the "commutativity" of the operation of the bi-graded symmetrization and the substitution of the momenta fields by their corresponding operators: where the bi-graded symmetrization of H(Φ) follows the same rules described above, with the only difference being that the fields Φ are not operators; rather, they are c-numbers.
Two other useful observations are that where X are some arbitrary functions of Φ, and as follows from Equations (191) and (192). Using these properties and Equation (186), one withd(Φ) given in Equation (184).
One can now proceed straightforwardly starting with the expression, The symmetrization in the first term related to the flow vector field is trivial because all the operators commute so thatd being the Lie derivative along e a (see Equation (43)). As can be straightforwardly verified, Using this martingale operator ordering rule, one can now find the SEO for the Ito interpretation of SDEs. Again, to facilitate the derivation, one notes that Equation (194) are also correct for the martingale ordering. Thus, the Ito SEO is alsod-exact, i.e., and the Ito version of Operator (198) iŝ with the subscript "mart" denoting the martingale ordering. One can now proceed straightforwardly as follows: where F 0 is the shifted flow vector field defined in Equation (91).
The previous equation and Equations (199) and ( The importance of this result is that it was obtained for the entire SEO and not only for the FP operator acting on only top differential forms, as done in Appendix A 1.

E. Generating Functional and Correlators
Various correlators and expectation values in the theory can be established through the introduction of the generating functional Here, the periodic/anti-periodic boundary conditions for the bosonic/fermionic fields are used, and J α is a set of external "probing" fields coupled to the system via a set of operators O α , sometimes called observables. We now present a few examples of observables with the corresponding probing fields: where f (x) is a function on X. In models with linear phase spaces, an observable of the form can also be used. The generating Functional (206) can be thought of as the DPF perturbed by probing fields. In particular, Z tt (0) = Z tt .
In practice, what one is interested in is the limit of the infinitely long evolution The limit here does not exist for interesting cases of the spontaneously broken supersymmetry because the DPF grows exponentially in this limit (see Equation (157)). Taking this limit must always be the very last operation after other manipulations are already performed. For example, the following notation for the family of correlators must be understood as The next goal now is to pass to the operator representation of the theory. Again, models with the Gaussian white noise will be considered.
Previously, the following relations were established for the DPF: Their analogue for the generating functional is The chronologically ordered exponent in Equation (212) can be given in the form of the formal Taylor series similar to Equation (35), To establish the operator representation of Correlators (210), one first performs the chronological ordering in the denominator of Equation (210): where t k > t k−1 ... > t 1 is the chronologically ordered permutation of t 1 ...t k . Accordingly, α 1 ...α k is the same permutation of α 1 ...α k , and (−1) P is the sign that may appear if this permutation is odd for fermionic operators. Further, where ψ n |M t + t k = ψ n |e −(t + −t 1 )En andM t 1 t − |ψ n = e −(t 1 −t − )En |ψ n have been used.
The contribution from each eigenstate includes a factor e −(t + −t − )En . In the limit of the infinitely long temporal evolution t ± → ±∞, only the contribution from the "ground" states with the least real part of their eigenvalue, i.e., ReE n = Γ g = min n ReE n , survive.
All other eigenstates provide exponentially vanishing contributions that can be neglected. Therefore, For the spectra presented in Figures 5a and 5b, the eigenvalue of the ground states is unique, and the situation is relatively simple. For example, the partition function takes the following form: where N g is the number of the groundstates. When E g = 0, N g = 2 because the ground states are the bosonic-fermionic pair. In the case E g = 0, i.e., the situation of unbroken supersymmetry, the ground states ared-symmetric, and N g must be the sum of the Betti numbers because each de Rahm cohomology class must provide oned-symmetric eigenstate.
Using Equations (214) and (215), the Correlators (210) take the following form: where the time arguments of the operatorsÔ(t) do not suggest that these operators have explicit dependence on time but rather indicate the moments of time that these operators act on the wavefunction.
In the Heisenberg representation, the wavefunctions are viewed independent of time, and the temporal evolution is passed onto the operators that now have explicit time dependences: with t * being some reference time moment. In this representation, Equation (216) is even simpler: gives The time independence of these expectation values can be interpreted as the ergodicity of the model. Thus, the ergodicity in STS is the property that the stochastic expectation values in the limit of infinitely long temporal evolution are those over the ground state(s).
Models with the spectra in Figure 5a,b are automatically ergodic.
For operators that are functions on X, i.e., O(t) = f (x(t)) ∈ Ω 0 (X), the expectation value has the following form: where P g (x) =ψ g (x) ∧ ψ g (x)/N g ∈ Ω D (X) is the TPD averaged over the ground states.
The time independence of P g (x) can be misinterpreted here as though the model is in the state of thermodynamic equiliubrium. In other words, the ergodicity can be mistaken for thermodynamic equilibrium. It is actually very common in the literature that the ergodicity and be confused with thermodynamic equilibrium. If the ergodicity and thermodynamic equilibrium were indeed equivalent, the concept of the "ergodic theory of chaos" [85] would not make sense (see also the last paragraph of Section III F). In other words, ergodicity is not equivalent to thermodynamic equilibrium.
The point here is that the TPD in Equation (220) is not the wavefunctions themselves but rather the bra-ket combination. This situation is similar to that in quantum theory, where it is the bra-ket combinations of the eigenstates that are the TPDs. The fact that this combination (the diagonal element of the density matrix) is stationary in time by no means implies that the eigenstate itself has zero eigenvalue.
For models with spectra given in Figure 5c, the situation is more subtle because there is a pair of Ruelle-Pollicott resonances with two different eigenvalues, i.e., E g and E * g , with the same "attenuation rate" ReE g = ReE * g = min n ReE n . These states are equally good candidates for the title of the ground state of the model. In the limit of the infinitely long temporal evolution, the DPF is Z t + −t − ≈ 4cos(t + − t − )ImE g e |ReEg|(t + −t − ) . This invalidates Equations (215) and (218) unless some additional arguments can circumvent this problem.
The same problem concerning the identification of the ground state exists in quantum theory. There, the finite-time quantum evolution operator isÛ t + t − = e −(iĤq)(t + −t − ) , whereĤ q is some Hermitian Hamiltonian with real spectra. To ensure that the generating functional in the long time limit receives contribution only from the ground state(s) with the lowest possible eigenvalue ofĤ q , one can Wick-rotate time "a little", i.e., t → t + i0 + , with 0 + being a vanishingly small positive constant. This approach can be borrowed for the STS, as illustrated in Figure 8. After Wick rotating time a little, only the ground states with the least "energy", i.e., ImE n , survive the limit of the infinitely long temporal evolution. Once this is done, the passage from Equations (210) to (216) and all the later formulas become valid.
As already mentioned in Section III C 3, the SEO possesses the ηT-symmetry, and each eigenstate with a complex eigenvalue must break this symmetry. By declaring one of the   8. (a) In quantum theory, the ground state is the one with the lowest energy. This can be justified by Wick rotating time a "little", i.e., t → t + i0, so that in the limit of infinitely long temporal evolution, the generating functional Z t = T re −iĤt receives a contribution only from the ground state with the least eigenvalue: Z t | t→∞ → e −iHgt ; (b) A similar approach can be used in the supersymmetric theory of stochastics (STS). If a Ruelle-Pollicott resonance has the smallest attenuation rate ReE, the ground state can be thought to be the one with the smallest ImE, i.e., the parameter analogous to the quantum energy.

F. One Way to a Unique Ground State
At this point of the discussion of the STS, the ground states are not unique for all three types of spectra given in Figure 5. Indeed, for models with unbroken topological supersymmetry, there may be many supersymmetric states, each of which may be viewed as a ground state of the model. For models with spontaneously broken topological supersymmetry (Figure 5b,c), the non-d-symmetric ground state is doubly degenerate because it is a boson-fermion pair of eigenstates.
The ground state can be made unique using yet another additional reasoning. This reasoning follows from the analysis of the supersymmetric states of the integrable models in the deterministic limit in Section V C 1. There, it will be discussed that the supersymmetric states of integrable deterministic models are the so-called Poincaré duals of the global unstable manifolds of the flow. One example of this situation is given in Figure 9 for where the fact that the bra-ket combination of each of these supersymmetric states is a δ-functional TPD on the corresponding critical points has been used.
On the other hand, it is intuitively clear that this expectation value must equal f (A).
To bypass this controversy, one can propose to view the ground state with the maximal number of fermions as the true ground state of the model. This rule can be called the principle of "minimal knowledge" for the following reason. The presence of a fermion in a wavefunction means that the wavefunction is a distribution in the corresponding bosonic variable, whereas the absence of a fermion suggests that the corresponding bosonic variable is not "thermalized" so that something else (e.g., the bra of the wavefunction or an external observer) must know with certainty the value of this bosonic variable. In other words, the more fermions a wavefunction has, the less external knowledge one needs to view the wavefunction as a "complete" probability distribution.
With this principle at hand, the ground state of the model is unique. When the topological supersymmetry is unbroken, the ground state is the TE state. For the broken supersymmetry case, the ground state isd-exact, i.e., |ϑ G =d|ϑ G , where the notations of Equation (134) have been used. Now that the ground state is unique, the correlators take the familiar field-theoretic form of the "vacuum" correlators: where J c (t), c = 1, 2..., is a set of probing fields and f is a set of predetermined vector fields.
The action of the model transforms accordingly: The methodology of Section IV E applies, with the perturbation operators or observables being Here, the notation is switched from O to L to note that these perturbation operators are Lie derivatives in the operator representation: withî f c = f c,i (x(t))∂/∂χ i . In the Heisenberg representation, becaused is commutative withĤ so thatd H (t) = e (t−t * )Ĥd e (t−t * )Ĥ =d. The expression for the response correlators follows now from Equation (222): Here, the fact that a product of thed-exact Lie derivatives from Equation (227)  On the contrary, if the topological supersymmetry is spontaneously broken, some of the perturbation correlators do not vanish. This can be interpreted as though the model "remembers" perturbations even in the limit of the infinitely long temporal evolution. This is how the STS reveals the famous butterfly effect. It is worth noting that the butterfly effect was previously often viewed as an intrinsic part of the definition of (deterministic) chaos, whereas within the STS, it is a derivable consequence.
The butterfly effect derived above is a part of a more general statement known as the Goldstone theorem. This theorem states that a model must exhibit a long-range order under the conditions of the spontaneous breakdown of a continuous global symmetry. In spatially extended models, this tailors the existence of a gapless excitation called the Goldstone-Nambu boson for bosonic symmetries and the goldstino for supersymmetries. This long-range order associated with thed-symmetry breaking is the DLRO discussed in the Introduction. Transient dynamics is often referred to as out-of-equilibrium dynamics. The same term is often used for the characterization of chaotic behavior. In these two situations, the term "out-of-equilibrium" has two different meanings. In one case, it means non-ergodic dynamics out of the global ground state of the DS, whereas in the second case, it denotes dynamics out of thed-symmetric state of the thermodynamic equilibrium but within the global nond-symmetric ground state. This second type of the "out-of-equilibrium" dynamics is often called "self-sustained" dynamics, i.e., happening forever.
In this section, only ergodic or self-sustained dynamics is addressed. Transient dynamics is beyond the scope of this paper. It is worth mentioning, however, that in some cases (glasses, for example) it must be possible to map transient dynamics in a model onto an ergodic dynamics in another model with the spontaneously broken supersymmetry.

B. Unstable Manifolds and Ground States: Langevin SDEs
In this subsection, the relation betweend-symmetric ground states in the weak noise limit and unstable manifolds of flow vector fields will be discussed. It is convenient to start the discussion with Langevin SDEs -the most studied class of SDEs closely related to N = 2 supersymmetric quantum mechanics (see, e.g., [99]). For simplicity, the noise-induced metric is assumed to be Euclidean: e i a = δ i a , g ij = δ ij . The flow vector field F i (x) = −δ ij U j (x) is defined via the Langevin potential U (x), and U j = ∂U/∂x j . The SEO of this model is given by Equation (123) The similarity transformationÂ →Â U = e U/(2Θ)Â e −U/(2Θ) acts on the SEO aŝ with (χ i ) † = δ ij ∂/∂χ j and (∂/∂x i ) † = −∂/∂x i .
BecauseĤ andĤ U are related via a similarity transformation, their spectra are identical.
As to the eigenstates, they are related as Up to the factor Θ, the operatorĤ U is the Hermitian Hamiltonian of N = 2 supersymmetric quantum mechanics. Its spectrum is real and non-negative. This implies that the topological supersymmetry is never broken in this class of models as long as there exists at least oned-symmetric ground state of the thermodynamic equilibrium (see Section III C 8).
In the single-variable case with the harmonic potential U = ωx 2 , the zero-eigenvalue ground state ofĤ U is (see, e.g., Section 10.2.4 in [99]) Here, the relation between bras and kets is trivial becauseĤ U is Hermitian. In terms of the eigensystem of the original non-HermitianĤ, the bra and ket are different. Using Equation (234), one has These are the ground state wavefunctions of the two models in Figure 3. For the stable variable case (ω > 0), the ket of the ground state is the narrow distribution around the stationary position x = 0, and the bra is not a distribution; rather, it is a constant function.
In the unstable case (ω < 0), the bra and ket are switched. As a result, the wavefunction is a narrow distribution in stable variables and is a constant function in unstable variables of the unstable manifold of this critical point, as illustrated in Figure 10. The small gray arrows in the transverse directions represent differentials/fermions. For integrable (non-chaotic) flow vector fields, local unstable/stable manifolds can be glued into the global unstable/stable manifolds indicated as closed dashed curves from the two first homology classes of the phase space, which is assumed here to be a 2D torus. The exterior derivative annihilates the Poincaré duals of the closed global unstable manifolds (see Figure 11), which in this case is the wavefunction of (one of) the globald-symmetric ground state(s). The so-emerged wavefunctions are known as Poincaré duals. They appear in one version of Poincaré duality stating that, for each k-dimensional submanifold c k , there exists a differential form ψ c k ∈ Ω (D−k) such that c k ϕ (k) = X ϕ (k) ∧ ψ c k for all ϕ (k) ∈ Ω (k) . Using this terminology, the bra and ket of the locald-symmetric ground state on a non-degenerate critical point of a Langevin SDE and in the weak noise limit are the Poincaré duals of the local stable and unstable manifolds, respectively.
The local unstable manifolds have boundaries on the lower dimensional local unstable manifolds of more stable critical points. For example, in Figure  Each de Rahm cohomology class may contain more than one globald-symmetric ground state. This can be the case only in the strict deterministic limit for the following reason.
Each of suchd-symmetric ground states is a superposition of oned-symmetric ground state and ad-exact piece. This means that pairs of non-d-symmetric states accidentally have zero eigenvalues. Any noise will introduce exponentially weak tunneling effects that must lift this accidental degeneracy, leaving only oned-symmetric ground state in each de Rahm cohomology class.
The example of the global supersymmetric eigenstates for the Langevin SDE in the deterministic limit on a 2D torus is given in Figure 9. In Section IV E, this example was used to argue that, among all the supersymmetric states of a model with unbroken supersymmetry, one should choose the state of the thermodynamic equilibrium as the true ground state, within which various correlators and observables should be calculated.

Chaotic Models
The next goal is to analyze qualitatively the structure of the ground states in chaotic or non-integrable deterministic models. These ground state(s) must represent the dynamics on fractal or strange attractors. Just like in the integrable models above, strange attractors are formed by the intersection of the stable and unstable manifolds. The bra/ket of the ground state must represent (or rather be) the Poincaré duals of these manifolds. The (un)stable manifolds in chaotic deterministic models are not well-defined topological manifolds however.
They can fold on themselves in a recursive manner, as illustrated for the class of models known as "homoclinic tangle" in Figure 12a. The straightforward attempt to construct a Poincaré dual for such an unstable manifold leads to the ambiguity in the orientation of the manifold at the point of the accumulation of self-folding. The would-be Poincaré dual must have such a coordinate dependence that circling around the orbit produces a sign change. This functional dependence on the position of the orbit also suggests that the topological supersymmetry is broken.
This ambiguity has its analogues in quantum theory. For example, a non-rotationally symmetric electron wavefunction on a rotationally symmetric atom (p, d, f, ... orbitals) would be ambiguous at the origin if it did not vanish there, which is always the case. For the same reason, in the theory of superfluids, the superfluidic order parameter of the Bose condensate at the core of a vortex must vanish.
In the case of the homoclinic tangle in Figure 12a, the ambiguity of the Poincaré dual of the unstable manifold can be remedied by modifying it with a continuous function that vanishes at the origin. This will introduce the coordinate dependence along the unstable manifold, and this coordinate dependence automatically suggests that the wavefunction is not annihilated byd; thus, the ground state wavefunction representing the unstable manifold is non-d-symmetric.
Another way to see that the ground state in a chaotic model is non-d-symmetric can be borrowed from the topological theory of chaos [76]. There, the global unstable manifold is qualitatively represented by a branched manifold that has self-intersections (see Figure   12b). The action byd on the Poincaré dual of this branched manifold is the Poinaré dual of its self-intersection. Thus, such a wavefunction is non-d-symmetric.
Yet another way to convince oneself that chaotic deterministic models have non-dsymmetric ground states is to recall that strange attractors contain an infinite number of unstable periodic orbits with arbitrary large periods. Some of these orbits have nonorientable local unstable manifolds, as illustrated in Figure 12c. The Poincaré duals of these local unstable manifolds must be modified by such coordinate dependence that going along the periodic orbit changes the sign of the wavefunction. Again, such a wavefunction is non-d-symmetric.
The above qualitative analysis of the ground states of the deterministic chaotic models is only an indication that the topological supersymmetry breaking must be the field-theoretic essence of deterministic chaos. The rigorous proof of this statement is given by Equation (157), which establishes the exponential growth of periodic solutions, being definitive for chaos, as well as by the emergence of the butterfly effect discussed in Section IV G.
D. Stochastic Models: Two Types of "Border of Chaos" In deterministic models, thed-symmetry is spontaneously broken or not depending on whether its flow vector field is non-integrable (chaotic) or integrable in the sense of DS theory.
The stochastic generalization of this picture is the subject of interest in this subsection.
One important thing to note is that in the high-temperature limit, the SEO (81) is dominated by the diffusion Laplacian. In a wide class of models (e.g., torsion-free vielbeins [74]) the diffusion Laplacian equals the Hodge Laplacian (58). The latter has real and nonnegative spectra, which correspond to the unbrokend-symmetry, so that thed-symmetry must always be unbroken at sufficiently large temperatures. Only models of this type are of interest here. It can be said that the noise destroys the DLRO at sufficiently high temperatures in this class of models.
Two qualitatively different types of the "border of chaos" exist for this class of models (see Figure (13)). For the first type, thed-broken phase gradually narrows with increasing temperature, which corresponds to the situation discussed, e.g., in [29]. The second type (Figure 13b) is more involved. There, thed-broken phase first widens with increasing temperature before shrinking. In other words, there exists a phase with an integrable flow vector field on one hand and withd-symmetry spontaneously broken on the other. This peculiar phase can be called noise-induced chaos (N-phase) because the supersymmetry can be restored by decreasing the temperature. In the deterministic limit, the N-phase collapses into the boundary of the deterministic chaos. 13. Two types of the "border of chaos". (a) In type I, there are only two phases: the chaotic phase (C) and the phase of the thermodynamics equilibrium (T). In the low-temperature limit, the topological supersymmetry of the C-phase is broken by the non-integrable flow vector field. As the temperature increases, the border moves to the "right" because the noise has the tendency to destroy the DLRO; (b) In type-II phase diagrams, there is an additional phase of the noise-induced chaos (N) where the flow vector field is integrable but where the topological supersymmetry is spontaneously broken by the condensation of (anti-)instantons, i.e., the noiseinduced tunneling processes between, e.g., different attractors. One type of dynamical behavior in the N-phase is such that an external observer sees a sequence of unpredictable jumps between patterns of "regular" behavior and/or attractors. This type of dynamics can be recognized as that of self-organized criticality. At higher temperatures, the sharp boundary between the N-and C-phases must smear out into a crossover because the perturbative supersymmetric ground states overlap significantly, and it is not possible for an external observer to tell one (anti-)instantonic process from another. The N-and C-phases must merge into a complicated phase (X) with the spontaneously broken topological supersymmetry.

Low-Temperature Regime and Self-Organized Criticality
In the low-temperature regime of the type-II phase diagram, thed-broken phase consists of two major subphases: the ordinary chaotic phase (C-phase), where thed-symmetry is broken by the non-integrable flow vector field, and the N-phase, where thed-symmetry is broken by some other mechanism. There are two other known mechanisms for the spontaneous breakdown of a symmetry. The first one is an anomaly, i.e., the possibility that a symmetry is broken by perturbative or fluctuational corrections. Supersymmetries, however, are difficult to break via anomaly. This fact is related to the so-called supersymmetry non-renormalization theorems [100,101]. This suggests that thed-symmetry breaking in the N-phase must be due to the other mechanism. This other mechanism of supersymmetry breaking is known as the condensation of (anti-)instantonic configurations [50]. Due to the renormalization theorems, the "dynamical" supersymmetry breaking by (anti-)instantons is considered as one of the most reliable mechanisms of supersymmetry breaking in high-energy physics models.
In case of stochastic dynamics, these (anti-)instantonic configurations are the tunneling processes that appear due to the noise-induced and exponentially weak overlap between perturbative ground states on unstable manifolds. One of the effects that the noise-induced tunneling processes will provide is the removal of the degeneracy of the deterministic zeroeigenvalue eigenstates representing Poincaré duals of "parallel" global unstable manifolds within the same de Rahm cohomology class discussed in the last paragraph of Section V C 1.
As a result, each de Rahm cohomology class will have only oned-symmetric eigenstate, whereas other eigenstates will acquire (exponentially small) non-zero eigenvalues. This removal of degeneracy does not necessary suggest thatd-symmetry is spontaneously broken.
Clearly, the very existence of the noise-induced tunneling processes is insufficient. Indeed, for a Langevin SDE (see Section V B) with a Langevin potential with multiple local minima, the tunneling processes between these local minima certainly exist at non-zero temperatures.
Nevertheless, thed-symmetry is never broken for this class of models. In other words, the weak-noise tunneling processes can only help the spontaneousd-symmetry breaking in models with flow vector fields that are close to being chaotic on their own. This is why the N-phase resides on the "border of chaos".
The physical picture of one type of dynamics in the N-phase is as follows. The fluctuating dynamics is mostly around unstable manifolds such as point attractors or limit cycles. The dynamics is sporadically interrupted by noise-induced tunneling processes or jumps between different attractors. Because it is the noise-induced tunneling processes that break thedsymmetry, the jumps must exhibit signatures of long-range dynamical behavior such as the power-law statistics. This power-law statistics of jumps, or avalanches as they also called in the literature, is a well-established phenomenon with the Richter scale for earthquakes being perhaps the best known example.
The ubiquitous power-law statistics of avalanches in nature was previously proposed to explain via the concept of self-organized criticality [18]. There, the power-law statistics is believed to be the signature of a gapless soft mode (see discussion in Section I B) associated with the ongoing phase transition into chaos, whereas the conspicuous contradiction with the fact that, unlike phase transitions, the N-phase has a finite width is circumvented by postulating of the existence of a mysterious force that fine-tunes the parameters of the stochastic model into the phase transition into chaos. This understanding of the essence of stochastic dynamics in the N-phase is all but scientific. The Goldstone mode explanation by the STS discussed above resolves this issue.

High-Temperature Regime
In the previous discussion of the weak-noise regime, the concept of noise-induced tunneling processes is well defined because the overlap between the perturbative ground states is exponentially weak. As a result, an external observer will be able to differentiate between tunneling process. At higher temperatures, the overlap is no longer weak, and it may become difficult for an external observer to differentiate between tunneling events. This suggests that the sharp boundary between the C-and N-phases must smear out into a crossover.
Note that the boundary between the N-and C-phases is not ad-symmetry-breaking phase transition; thus, its disappearance does not contradict any symmetry-based argument.
It can be said that, above a certain temperature, the C-and N-phases must merge into a complicated phase with spontaneously brokend-symmetry. In Figure 13b, this phase is indicated as an X-phase. Borrowing from the terminology of high-energy physics, one way to identify this phase is as stochastic chaos in the "strongly coupled" regime, where strong coupling would mean the strong overlap between the perturbative ground states.

VI. CONCLUSIONS AND OUTLOOK
This paper offers a brief introduction to the current state of the recently proposed approximation-free supersymmetric theory of stochastic differential equations (STS). This theory provides several novel theoretical insights into stochastic dynamics. It establishes a rigorous stochastic generalization of the concept of dynamical chaos, which is found to be the phenomenon of the spontaneous breakdown of topological or de Rahm supersymmetry that all SDEs possess. This paper also reveals that stochastic chaos is complementary to thermodynamic equilibrium, corresponding in turn to the unbroken topological supersymmetry.
Being the low-symmetry or the ordered phase, the chaotic phase has what can be called a dynamical long-range order, whereas the phase of the thermodynamic equilibrium does not. The presence of this order is the reason why many natural, engineered, and social DSs exhibit emergent long-range dynamical behavior such as 1/f noise, i.e., long-term memory effects; the butterfly effect, i.e., sensitivity to the initial conditions; and the algebraic, i.e., scale-free, statistics of sudden or instantonic processes. These and a few other qualitative findings, such as the clarification of the concept of ergodicity, are the main outcomes of the STS so far. Further work on the STS may lead to more specific and valuable results.
As is discussed next, one of the most fruitful directions of further investigation is the work on the methodology of the identification of the dynamical long-range order parameter and construction of the low-energy effective theory (LEET) for it in spatially extended models such as hydrodynamical models.
The most important qualitative aspect of dynamics under the conditions of the spontaneous breakdown of a global symmetry occupies a reduced phase space. In other words, this aspect occurs in a reduced number of "low-energy" variables called the order parameter, and the LEETs are the theories describing this dynamics. In ferromagnets, for instance, the order parameter is the (local) magnetization of the electron liquid, and the LEET (or rather the equations of motion of the LEET) is the Landau-Lifshitz-Gilbert equation. In superconductors, the order parameter is the wavefunction of the Bose-Einstein condensate of the Cooper pairs, and the LEET is the corresponding Ginzburg-Landau theory. In solids, in which the global translational symmetry is broken by the lattice structure, the order parameter is the local displacement of atoms from their average positions in the lattice, and the LEET is the low-energy theory describing, say, the propagation of transverse sound, which is the Goldstone-Nambu particle in this case.
Concerning chaotic DSs, the most important variables are the unstable and/or unthermalized variables of the wavefunction of the non-d-symmetric ground state. It is in these variables that a chaotic DS exhibits the infinite memory of perturbations. The local order parameter must be the gapless fermions or goldstinos that are the supersymmetric partners of the unstable bosonic variables. In spatially extended nonlinear models, the unstable variables must be the moduli of the solitonic configurations consisting of fundamental solitons such as kinks, domain walls, and vortices. The processes of the creation/annihilation of (pairs of) the fundamental solitons are the (anti-)instatonic processes, the condensation of which is the essence of the noise-induced chaotic phase discussed in Section V D 1. One of the candidates for such models is a two-dimensional vortex-mediated turbulence, wherein the goldstinos must be the supersymmetric partners of the spatial positions of the (anti-)vortices.
It is intuitively appealing to believe that the dynamical long-range order of the spontaneously broken topological supersymmetry must at least partially possess some global or topological features. These features may appear on the level of the inter-goldstino interactions. In the above example of the vortex-mediated 2D turbulence, the interaction that remembers the braining between the vortices may as well be approximated as a (chiral) gauge field. Such an LEET would appear somewhat reminiscent of the Schwartz-type topological field theories used in models related to the concept of topological quantum computing.
Somewhere down this line of thinking it may turn out that some complex DSs may be useful distribution. Therefore, one may as well expect that the wavefunction may change suddenly upon observation in a Bayesian update manner. If this is indeed true, yet another question arises of how this change is related to wavefunction collapse in quantum theory. Hopefully, future work will reveal answers to these and other open questions in the STS. It would also be interesting to see how the STS relates to other modern approaches to stochastics such as Stochastic Thermodynamics [113][114][115][116].

ACKNOWLEDGMENTS
The author's sincerest gratitude is offered to Kang L. Wang for valuable discussions, support, and encouragement.
with dots denoting terms of higher order in ∆t andĤ (D) being the sought after FP operator.
In the continuous-time limit, the above expression can be given the familiar form of the FP equation: The task now it to establish the explicit expression forĤ (D) .
It is understood that, for small ∆t, ∆x is also small. In other words, the Taylor expansion of ∆x in ∆t begins with the first-order term. Equation (A6) also seemingly implies that it suffices to retain only terms of first order in ∆t. This would indeed be true if it was not for the stochastic averaging over ξ n . This averaging will transform the terms that are second order in ∆t and contain two ξ into the first-order terms in ∆t, as is clear from Equation (7).
Therefore, the Taylor expansion up to second order in ∆x must suffice for the derivation of the FP operator.
The Taylor  (A8) Here, f j ≡ ∂f /∂x i and similar for f ij = ∂ 2 f /∂x i ∂x j , and the partial integration has been used.
The next step is to Taylor expand ∆x up to second order in ∆t and substitute this expansion into the above expression. Using Equation (A1), one has ∆x i = F i n ∆t + α(F i n ) j F j n ∆t 2 + ...
Substituting this expression into Equation (A8), using Equation (A2), and performing the stochastic averaging over ξ n with the help of Equation (7), one arrives at with the FP operator beinĝ and with the α-dependent flow vector field from Equation (91).
In the above derivation of the FP operator, SdE (A9) was used as a formal equation defining ∆x. One can take an alternative view on stochastic dynamics in which the dynamics is continuous in time and the noise is piece-wise constant, as given in Figure 1. For a fixed noise configuration, one has a continuous trajectory x(t), defined byẋ = F n (x(t)) with the initial condition x(t n−1 ) = x n 1 . Now, there is no freedom in choosing α because ∆x is uniquely defined by the evolution according to the Picard-Lindelöf theorem. In particular, ∆x has a unique Taylor expansion in ∆t: The first coefficient here is determined from the SDE itself, whereas the second coefficient is obtained via one differentiation of the SDE over time thus, the quantity in Equation (A12) becomes ∆x i = F i n (x)∆t + 1 2 F i n j (x)F j n (x)∆t 2 + ...
Comparing this equation with Equation (A9), one concludes that the Stratonovich choice of α = 1/2 must always be used for the continuous-time picture of temporal evolution.
Concerning the Ito interpretation of SDEs, it is often said that, unlike all other interpretations, the Ito approach respects the Markovian property in the sense that the increment ∆x n or, equivalently, the final point x n = x n−1 + ∆x n is a function of only x n−1 and not of x n . This advantage of Ito SDEs, however, is a misinterpretation. Indeed, the very statement that x n is a function of itself for α > 0 does not make sense from the point of view of functional dependence. This sentence only tells us that x n as a function of x n−1 is given only implicitly by Equation (90). For a fixed noise variable ξ n , the final point x n together with the increment ∆x n is always a function of x n−1 only. Its explicit expression is given by Equation (A9) up to second order in ∆t, the only accuracy relevant in the continuous-time limit.
Furthermore, the Markovian property of stochastic processes is concerned not with the trajectories (the variables x n and x n−1 ) but rather with the temporal evolution of TPDs. In application to the SdEs (90), the Markovian property means that the TPD at time moment t n depends on the TPD at the previous time moment t n−1 only and not on the TPD at earlier time moments. As clearly observed from Equation (A6), which is correct for all α, all the interpretations of SDEs satisfy this requirement of Markovianity. In other words, Ito SDEs are just as Markovian as SDEs in all the other interpretations.
In other words, the only advantage of the Ito interpretation is the relative ease of its numerical implementation because the increment as a function of x n−1 is given explicitly by the Ito SdE. This convenience for numerical implementations, however, does not have any significance from the mathematical point of view.

Perturbative Supersymmetric Eigenstates
The correspondence between supersymmetric states and de Rahm cohomology classes can be established using standard perturbation theory. The first step is to recall that the Hodge Laplacian from Equation (58) and |ζ n =d|ζ n , ζ n | = ζ n |, and their eigenvalues ∆ n > 0.
One can now split the SEO into two parts aŝ and viewV as a perturbation. The "zeroth-order" SEO, i.e.,Ĥ 0 , is elliptic, whereas the perturbation operator is only linear in spatial derivatives:V =f i ∂ ∂x i +ĝ. This implies that the perturbation series must be well defined, e.g., convergent, at least for some class of models and for sufficiently large Θ.
Similarly, the bra of thisd-symmetric state is θ n | = h n | + θ n |d.
Thus, within the domain of the applicability of the perturbation theory, each de Rahm cohomology class provides one supersymmetric eigenstate.

Spectra
In this Appendix, it is discussed how the theory of KD is related to STS and how this relation provides examples of the two supersymmetry-breaking spectra in Figure 5b,c.
The KD is a part of the more general hydromagnetodynamical phenomenon of the magnetic dynamo. The latter is the ability of a moving conducting medium to generate and/or sustain a magnetic field [102]. Many astrophysical objects exhibit magnetic dynamos, including galaxies [103,104], galaxy clusters [105], stars [106], and planets, including the Earth [107]. In turn, the KD is the linear regime of a magnetic dynamo when a relatively weak magnetic field is generated by a stationary flow of the conducting medium. The KD is realized, e.g., in the early stages of the formation of galaxies.
The temporal evolution of the magnetic field within the KD effect is governed by the induction equation Here,∂ is the gradient operator of the Euclidean space X = R 3 ;ˆ =∂ i∂i is the standard Laplace operator; × denotes the vector product;∂× is the curl of a vector; B is the magnetic field vector; v is the vector field of the underlying flow of the conducting medium; and η = 1/σµ is the magnetic diffusivity, with σ and µ being the electrical conductivity and permeability, respectively. The first term represents the well-known magnetohydrodynamical phenomenon of the "freezing" of the magnetic field into the conducting medium, whereas the second term is the magnetic field diffusion.
In the theory of KD, the "phase space" is non-compact R 3 . On the other hand, in this paper, compact (and closed) phase spaces are under consideration. This problem can be circumvented. The point is that the spatial structures of the KD magnetic fields have local support, and one can always compactify the phase space into the 3D sphere at spatial infinity without affecting the structure of the KD magnetic fields.
Equation (A22) can be presented in a coordinate-free form. Instead of the vector B, one can equivalently use the 2-form representing the magnetic field where A = A i dx i is the 1-form of the vector potential. In components, B ij , called the magnetic field tensor, is where ijk is the antisymmetric Levi-Civita tensor. Equation (A22) can now be expressed as ∂ t B i = e ipq∂ p e qkl v k B l + ηˆ B i .
Lowering and raising the indexes in the Euclidean space has no effect on the values of the components of the antisymmetric tensor, e.g., e ijk = e ijk . Using the identity and ∂ i B i = 0, Equation (A22) can be rewritten as where v i j = ∂ j v i . Now, using the induction equation can be further transformed as Multiplying both sides of this equation by e iab and summing over index i, one arrives at Using the identity e iab e jkl = det The first term here is the Lie derivative applied to the 2-form (A23). Therefore, Equation (A22) can also be given as This result is rather natural. As previously mentioned, the first term in the R.H.S. of Equation (A22) is the infinitesimal temporal evolution of the magnetic field "frozen" into the conducting medium. This freezing is the evolution solely due to the flow along v, and such an evolution is given by the Lie derivative. This also explains why the Lie derivative is also known as the physical derivative.
The Laplacian in the Euclidean space is given aŝ whered † = −ı i δ ij ∂ j is the codifferential operator defined in Equation (57). Thus, the KD evolution operator isĤ whered =î v − ηî i δ ij∂ j and where the identity δ ij = δ i a δ j a has been used. It is now clear that Equation (A27) is the SEO of the following SDE:ṙ with ξ(t) ∈ R 3 being Gaussian white noise.
This is the result needed to establish that the supersymmetry-breaking spectra of both types in Figure 5b,c are realizable. Indeed, it is well established that the eigenvalues of the KD operator with the lowest real part can be not only negative but also complex (see, e.g., [108] and the references therein). The complex eigenvalues indicate that the spatial structure of the growing magnetic field is also rotating.