On Stochastic Quantisation of Supersymmetric Theories

TL;DR

This paper explores how to reconcile stochastic quantization with space supersymmetry in supersymmetric theories, demonstrating that the stochastic process converges to the standard path integral regardless of certain kernel choices.

Contribution

It introduces supersymmetry-compatible kernels in stochastic quantization for N=2 supersymmetric quantum mechanics, ensuring convergence to the path integral.

Findings

Kernels depend on a parameter M that influences stochastic evolution

The infinite stochastic time limit is independent of M

Compatibility of stochastic TQFT supersymmetry with space supersymmetry is established

Abstract

We explain how stochastic TQFT supersymmetry can be made compatible with space supersymmetry. Taking the case of N=2 supersymmetric quantum mechanics, (the proof would be the same for the Wess-Zumino model), we determine the kernels that ensure the convergence of the stochastic process toward the standard path integral, under the condition that they are covariant under supersymmetry. They depend on a massive parameter M that can be chosen at will and modifies the course of the stochastic evolution, but the infinite stochastic time limit of the correlation functions is in fact independent on the choice of M.

Equations

$$<f(\phi )>\equiv \int [\mathrm{d}\phi {]}_{x}f(\phi )\exp -S[\phi ]=\underset{\tau \to \infty }{\lim }<<f({\phi }^{\eta }(x,\tau )>{>}^{K,S}$$

$$\frac{\mathrm{d}{\phi }_a}{\mathrm{d}\tau }$$

$$\displaystyle<<f({\varphi}^{\eta})>>^{K}\$$

$$<<\phi (x_1,\tau ),..,\phi (x_n,\tau )>>=\int [\mathrm{d}\phi {]}_x\phi (x_1),..,\phi (x_n)P^{S,K}(\phi ,\tau )$$

$$\displaystyle\frac{\partial P^{S,K}({\varphi},\tau)}{\partial\tau}=\int{\rm d}x\frac{\delta}{\delta{\varphi}_{a}x}K_{ab}\Big{(}\frac{\delta}{\delta{\varphi}_{b}(x)}+\frac{\delta S}{\delta{\varphi}_{b}(x)}\Big{)}P^{S,K}({\varphi},\tau)$$

$$P^{S,K}(\phi ,\tau =\infty )=\exp -\int S[\phi ]$$

$$\displaystyle P^{S,K}({\varphi}(x),\tau)=\Big{(}\exp\tau\int{\rm d}x\frac{\delta}{\partial{\varphi}_{a}(x)}K_{ab}(\frac{\delta}{\delta{\varphi}_{b}(x)}+\frac{\delta S}{\delta{\varphi}_{b}(x)})\Big{)}P^{S,K}(V,\tau=\tau_{0})$$

$$\displaystyle{\varphi}_{a}(x)\to\Big{(}{\varphi}_{a}(x,\tau),\Psi_{{}_{a}}(x,\tau),\overline{\Psi}_{a}(x,\tau),b_{a}(x,\tau\Big{)}$$

$$s_{stoc}\phi =\Psi ,s_{stoc}\Psi =0,s_{stoc}\overline{\Psi }=b,s_{stoc}b=0$$

$$\displaystyle\int[{\rm d}m]_{\tau}[{\rm d}n]_{\tau}[{\rm d}{p}]_{\tau}[{\rm d}q]_{\tau}\exp\int d\tau\Big{[}n\ ({\cal M}(m)+\alpha)\mp{p}\frac{\delta{\cal M}(m)}{\delta m}q\ \Big{]}=1$$

$$\displaystyle{s_{stoc}}\Big{(}p({\cal M}(m)+\alpha)\Big{)}=n({\cal{M}}(m)+\alpha)\mp p\frac{\delta{\cal M}}{\delta m}{q}$$

$$<<f(\phi )>{>}^K=\int [\mathrm{d}{\eta }_a{]}_{\tau }f(\phi )\exp -\frac{1}{2}\int \mathrm{d}t{\eta }_a{K}_{ab}^{-1}{\eta }_b$$

$$=$$

$$X(t,\theta ,\hat{\theta })=q(t)+\theta \lambda (t)+\hat{\theta }\hat{\lambda }(t)+\theta \hat{\theta }A(t)$$

$$D_{\theta }\equiv {\partial }_{\theta }+\theta {\partial }_t{D}_{\hat{\theta }}\equiv {\partial }_{\hat{\theta }}+\hat{\theta }{\partial }_t$$

$$\{D_{\theta },D_{\theta }\}=\{D_{\hat{\theta }},D_{\hat{\theta }}\}=0\frac{1}{2}\{D_{\theta },D_{\hat{\theta }}\}={\partial }_t$$

$$\delta q$$

$$\hat{\delta }q$$

$${\Sigma }_0\equiv \int \mathrm{d}t\mathrm{d}\theta \mathrm{d}\hat{\theta }\mathrm{d}\tau (\frac{1}{2}X{D}_{\hat{\theta }}{D}_{\theta }X)=\int \mathrm{d}t(\frac{1}{2}{A}^2-\frac{1}{2}q{\partial }_t^{2}q-\frac{1}{2}\lambda {\partial }_t\lambda -\frac{1}{2}\hat{\lambda }{\partial }_t\hat{\lambda })$$

$${\Sigma }_W\equiv \int \mathrm{d}t\mathrm{d}\theta \mathrm{d}\hat{\theta }\mathrm{d}\tau (W(X))=\int \mathrm{d}t(A{W}_q(q)-\lambda W_{qq}(q)\hat{\lambda })$$

$$I_W\equiv {\Sigma }_0+{\Sigma }_W=\int \mathrm{d}t(\frac{1}{2}(A+W_q{)}^2-\frac{1}{2}{W}_q^2+\frac{1}{2}({\partial }_{t}q{)}^2-\frac{1}{2}\lambda {\partial }_t\lambda -\frac{1}{2}\hat{\lambda }{\partial }_t\hat{\lambda }-\lambda W_{qq}\hat{\lambda })$$

$$I_W^{t\to it}$$

$$=$$

$$\eta ={\eta }_q+\theta {\eta }_l+\hat{\theta }{\eta }_{\hat{l}}+\theta \hat{\theta }{\eta }_A$$

$$\dot{q}$$

$$\dot{A}$$

$$\dot{\lambda }$$

$$\dot{\hat{\lambda }}$$

$$X_W^{\ast }=\frac{\delta I_w}{\delta A}-\theta \frac{\delta I_w}{\delta \hat{\lambda }}+\hat{\theta }\frac{\delta I_w}{\delta \lambda }+\theta \hat{\theta }\frac{\delta I_w}{\delta q},$$

$$X_W^{\ast }=-A-W_q+\theta (-i{\partial }_t\hat{\lambda }+W_{qq}\lambda )-\hat{\theta }(-i{\partial }_t\lambda -W_{qq}\hat{\lambda })-\theta \hat{\theta }({\partial }_t^{2}q+A{W}_{qq}-\lambda W_{qq}\hat{\lambda })$$

$$\dot{X}(t,\theta ,\hat{\theta },\tau )$$

$$B\equiv b_A+\theta b_{\hat{\lambda }}+\hat{\theta }{b}_{\lambda }+\theta \hat{\theta }{b}_q$$

Full text

1 introduction

Giving the existence of stochastic quantisation, an interesting problem is to understand the way stochastic evolution can preserve space supersymmetry, or, more mathematically, how physical space supersymmetry and stochastic TQFT supersymmetry can be made compatible.

Stochastic quantisation computes the Euclidean quantum correlation functions of a given field theory in a $d$-space $\{x^{\mu}\}$ by interpreting this space as the $\tau=\infty$ boundary of a $d+1$-space, $\{x^{\mu}\}\to\{x^{\mu},\tau\}$. The stochastic time $\tau$ is in fact a bulk coordinate and stochastic quantisation gives a microscopic understanding of the Euclidean path integral measure $\exp-\frac{S}{\hbar}$ as an equilibrium distribution at $\tau=\infty$, analogously as Langevin equations determine the thermal equilibrium Boltzmann formula [1][2]. It is also a fact that a BRST-TQFT localisation procedure can enforce Langevin equations in a supersymmetric way so that the whole apparatus of stochastic quantisation can be expressed as a TQFT in $d+1$ dimensions [3]. This gives a good handle on the arguments in [4].

The drift force of the stochastic process is basically proportionate to the Euclidean equations of motion of the theory on the boundary, modulo some other forces that can possibly improve the convergence of the process. It is often the case that kernels that realise factors of multiplication of the equations of motion must be carefully chosen to improve, and, sometimes, to define, the convergency of the stochastic process [1]. For maintaining space supersymmetry, if any, appropriate kernels are needed. The goal of this note is to explain how they can be determined.

It is natural to ask how symmetries, which are realised on the boundary at $\tau=\infty$, are represented in the bulk, that is, at finite values of the stochastic time. For gravity and gauge symmetries, and even string theory, one solves the question rather elegantly in the framework of equivariant cohomology  [3][5]. It leads one to define a stochastic equivariant topological BRST invariance. One then obtains a universal frame for the stochastic quantisation of any given theory with a gauge invariance, including a method to fix the gauge invariance in the bulk with a BRST symmetry construction.

The case of theories with space supersymmetry is puzzling. The equations of motion have a different supersymmetry covariance than the fields. Thus, they cannot be used as drift forces without introducing appropriate kernels to modify consistently their covariance. In this note, we study in details the case of theories with global supersymmetry. We take the case of the simple $N=2$ supersymmetric mechanics and explain how the supersymmetry of the theory on the boundary can be enforced in the bulk. The case of the Wess and Zumino model can be handled in exactly the same way, modulo some notational complications, and we let the reader make the correspondence. The superfield formalism is handy to get Langevin equations that are compatible with supersymmetry for all values of the stochastic time. The method we follow to compute the appropriate kernels and enforce supersymmetry in the stochastic bulk can be actually applied to more delicate cases, including when supersymmetry mixes with gauge symmetries, as will be shown in a separate publication.

2 A reminder of kernels in stochastic quantisation

Consider a quantum field theory of fields ${\varphi}_{a}(x)$ with a local action $S[{\varphi}]=\int{\rm d}xL({\varphi}_{a}(x))$, where $a$ is some index that labels the fields, and the $x$’s denote Euclidean coordinates in $d$ dimensions. Stochastic quantisation introduces the stochastic time coordinate $\tau$, with ${\varphi}_{a}(x)\to{\varphi}_{a}(x,\tau)$. The $\tau$ evolution of the field ${\varphi}_{a}(x,\tau)$ is defined by a Langevin process, such that the path integral weight $\exp-S[{\varphi}]$ is its equilibrium distribution, when it exists. Given a correlation function $<f({\varphi}(x))>$ in the space $\{x^{\mu}\}=\{x^{\mu},\tau=\infty\}$, which can be computed by the usual Euclidean path integral on the space $\{x^{\mu}\}$, the fundamental property of stochastic quantisation is that

[TABLE]

The $(x,\tau)$ dependent correlators $<<f({\varphi}^{\eta}(x,\tau)>>^{K,S}$ are computed as an average over the noises $\eta_{a}(x,\tau)$. Here, ${\varphi}^{\eta}$ is computed in function of the noises $\eta_{a}$ by solving the following Langevin differential equation with some any given initial condition ${\varphi}^{\eta}(x,\tau=\tau_{0})={\varphi}_{0}(x)$,

[TABLE]

The dependance in $\hbar$ is through a rescaling of the noise $\eta$. The correlations of the noise are in fact defined by the following formula, true for any given sufficiently regular functional $f$

[TABLE]

It is often the case that $K_{ab}$ must be different from the trivial identity to ensure a well-behaved drift force in the r.h.s of the Langevin equation (3), and a proper behaviour of the stochastic time evolution at $\tau\to\infty$. One must indeed have a uniformly attractive force in the Langevin process for all components of fields. Thus $\eta$ is not necessarily a white noise, and this happens for instance for the quantisation of a spinor *** For a free Dirac action $S=\int{\rm d}x\overline{\lambda}(\gamma\cdot\partial-m){\lambda}$ one can use the following Langevin equation $\dot{\lambda}=(\gamma\cdot\partial+m)((\gamma\cdot\partial-m){\lambda}+\eta_{\lambda}=(\partial^{2}-m){\lambda}+\eta_{\lambda},$ with fermionic noise correlation functions $<<\eta_{\lambda}(x,\tau)\eta_{\lambda}(x,\tau)>>=\delta(x-x^{\prime})\delta(\tau-\tau^{\prime})\frac{1}{\gamma\cdot\partial+m}$. It gives the correct result $P^{K=\gamma\cdot\partial+m}({\lambda},\tau=\infty)=\exp-\int{\rm d}x\overline{\lambda}(\gamma\cdot\partial-m){\lambda}$. If, instead, we take $\dot{\lambda}=(\gamma\cdot\partial-m){\lambda}+\eta_{\lambda}$ the Langevin process is ill-defined. In the former case the eigenvalues of the drift force $(\partial^{2}-m){\lambda}$ are positive, and the Langevin process converges. In the later case, the drift force is $(\gamma\cdot\partial-m){\lambda}$, with positive and negative eigenvalues, so that the limit $\tau=\infty$ cannot be reached. .

The dependence on $K_{ab}$ evaporates when one reaches the limit $\tau=\infty$. Indeed, given different kernels $K$ that ensure a convergence of the stochastic process, the value of $\lim_{\tau\to\infty}<<f({\varphi}(x,\tau)>>^{K,S}$ is independent on the choice of $K$, and is given by the standard part integral formula in the r.h.s. of Eq. (1).

This property can be demonstrated from the Fokker–Planck equation that is implied by the $K$ dependent Langevin equation. Indeed, the Langevin equation (2) implies the existence of a (Fokker–Planck) kernel $P^{S,K}({\varphi}(x),\tau)$, which permits one to computes the time evolution of equal time stochastic correlators as follows

[TABLE]

Here, $P^{S,K}$ is the solution of the Fokker–Planck equation

[TABLE]

This equation needs an initial condition, eg: $P^{S,K}({\varphi}(x),\tau=\tau_{0})=\delta({\varphi}(x)-{\varphi}_{0}(x))$. For sufficiently regular theories, one can check an exponential damping $\sim O(\exp-\tau)$ in the dependence on the initial condition ${\varphi}_{0}(x)$ when $\tau\to\infty$. The details of the evolution depend on the choice of $K$, but not the limit. If a stationary distribution $\lim_{\tau\to\infty}P^{S,K}({\varphi}(x,\tau)$ exists, Eq. (5) implies

[TABLE]

independently of the choice of $K$. At finite $\tau$, the $K$-dependance of $P^{S,K}$ is, explicitly :

[TABLE]

This standard result of statistical mechanics can be simplified by the use of stochastic TQFT supersymmetry [6][4]. There is indeed a supersymmetric path integral that expresses the stochastic process in the bulk $\{x,\tau\}$, where ${\varphi}_{a}(x)\to{\varphi}_{a}(x,\tau)$ [1][4]. It enlarges the phase space of the fields as follows:

[TABLE]

The stochastic supersymmetry acts on this enlarged space of fields and it is defined by a nilpotent graded differential operator ${s_{stoc}}$ [3]. ${s_{stoc}}$ is alike a BRST topological symmetry operation, which justifies the name stochastic BRST supersymmetry, where ${\varphi}_{a},\Psi_{a},\overline{\Psi}_{a},b_{a}$ are identified as two trivial BRST multiplets under ${s_{stoc}}$ and

[TABLE]

If there a gauge symmetry acting on ${\varphi}_{a}(x)$, some refinements are necessary and one must build an equivariant supersymmetry [5]. The physical conclusions remain the same for the relaxation to a stable equilibrium and they justify the gauge theory path integral formula on the boundary.

We now indicate a universal formula that will be useful in the following. It leaves aside the unessential indices $x^{\mu}$ that can be discretised. Given a pair of boson fields $(m(\tau),n(\tau))$ and a pair of fermion fields $({{p(\tau)}},{q(\tau)})$, or the reversed case where on interchanges bosons and fermions, determinant formula imply that one has indeed the following identity

[TABLE]

where ${\cal M}(m)$ is any given local functional of $m$. $\alpha$ is independent on $m,n,p,q$ and the formula holds true whatever $\alpha$ is. Defining ${s_{stoc}}m=q,{s_{stoc}}q=0,{s_{stoc}}p=n,{s_{stoc}}n=0$ and ${s_{stoc}}$, one has

[TABLE]

${s_{stoc}}$ acts as a graded differential operator, so ${s_{stoc}}^{2}=0$ on any functional of $m,n,p,q$.

The identity (10) is important, as it explains the chain of identities in the next equation (12). Indeed, replacing $\alpha\to\eta$, $m\to{\varphi},n\to b,p\to=\overline{\Psi},q\to\Psi$ and ${\cal M}({\varphi})\to\frac{{\rm d}{\varphi}}{{\rm d}t}+K\frac{\delta S}{\delta{\varphi}}$ in the functional identity (10), one can write the following succession of equations :

[TABLE]

These equalities hold modulo normalising factors for each line. The first line enforces the definition of $<<f>>^{K}$ as an average over noises, using the Langevin equation $\frac{{\rm d}{\varphi}_{a}}{{\rm d}t}=-K_{ab}\frac{\delta{S}}{\delta{\varphi}_{b}}+\eta_{a}$ and an arbitrarily given initial condition ${\varphi}(\tau_{0})$ at some reference time $\tau_{0}$ that one must introduce to solve  ${\varphi}$ in function of $\eta$. The last line teaches us that the correlator $<<f({\varphi})>>^{K}$ is computed by the path integral of a topological quantum field theory, with a localising gauge function related to the Langevin equation. In other words, the path integral of the ${s_{stoc}}$-exact Lagrangian in the last line defines the same correlation functions as those computed either from the Langevin equations, or from their corresponding Fokker-Planck equations (if one restricts oneself to equal stochastic time correlators). This supposes that the possible zero modes of the operators that stand between the ghosts $\Psi$ and $\overline{\Psi}$ are handled properly with periodic boundary conditions for $\Psi$.

The Hamiltonian that defines the stochastic time evolution is the $\tau$ Legendre transform of the supersymmetric Lagrangian expressed in the third or fourth lines of Eq. (12), as generically explained in [4].

Eq. (12) is valid whether ${\varphi}_{a}$ is bosonic or fermionic, so it can applied when ${\varphi}_{a}$ stands for any one of the field components $x,{\lambda},{\hat{\lambda}},A$ of the superfield $X(t)$ that we will shortly introduce.

In the stochastic BRST-TQFT supersymmetric formulation, the proof of the K independence of the $\tau=\infty$ limit of correlators, if they exist, follows from the properties of determinants, for the bosonic case as well as for the fermionic case for ${\varphi}$.

3 Supersymmetric kernels

Let us now consider a supersymmetric theory with a set of auxiliary fields that gives a closed system of transformations. We will prove that supersymmetry is compatible with the stochastic time evolution.

For the notational simplicity of the proof we consider the $N=2$ supersymmetric quantum mechanics. One has $\{x^{\mu}\}=\{t\}$, two supersymmetry generators $\delta_{i}=(\delta,\hat{\delta})$ and one multiplet (1,2,1), that is, one propagating boson ${q}({t})$, two fermions ${\lambda}({t})$ and $\hat{\lambda}({t})$ and one auxiliary field $A({t})$. The construction we will detail can be generalised to all other multiplets. The proof can be generalised to the Wess and Zumino model with no difficulty but notational complications. We first consider that $t$ is the real time and we will shortly shift to the Euclidean time, $t\to it$. Call $X(t,\theta,{\hat{\theta}})$ the superfield

[TABLE]

where $\theta$ and ${\hat{\theta}}$ are Grassman coordinates, with $\int{\rm d}{\hat{\theta}}=\int{\rm d}\theta=0$ and $\int{\rm d}{\hat{\theta}}{\hat{\theta}}=\int{\rm d}\theta\theta=1$. Define also

[TABLE]

The action of supersymmetry transformations on the superfield $X$ is given by $Q=\partial_{\theta}-\theta\partial_{t}$ and $\hat{Q}=\partial_{\hat{\theta}}-{\hat{\theta}}\partial_{t}$. By expansion in components, one finds the following action of both supersymmetries.

[TABLE]

They satisfy $\frac{1}{2}\{\delta_{i},\delta_{j}\}=\delta_{ij}\partial_{t}$ and $\delta_{i}{\rm d}+{\rm d}\delta_{i}=0$. The massless supersymmetric free action is

[TABLE]

A general supersymmetric interaction, including a mass term, is determined by the prepotential $W(X)$

[TABLE]

($W_{q},\ W_{qq}$ mean $\partial_{q}W,\ \partial_{q}^{2}W$). So, the generic standard interacting supersymmetric action that we consider is

[TABLE]

With $t\to it$, one gets the following equivalent (modulo boundary terms) expressions of the Euclidian action †††As all supersymmetric theories, this model has a twisted formulation, where one can rewrite the supersymmetry algebra as a nilpotent one. In this case, one redefines $\psi=\frac{{\lambda}-i{\hat{\lambda}}}{\sqrt{2}}$ and $=\frac{{\lambda}+i{\hat{\lambda}}}{\sqrt{2}i}$ By doing this change of variables, the Euclidean action becomes ${I_{{\rm{twisted}}}}=\int{\rm d}t(\frac{1}{2}A^{2}+W_{q}A-\frac{1}{2}(\partial_{t}{q})^{2}+\overline{\psi}\partial_{t}\Psi-\overline{\psi}W_{{q}{q}}\Psi)\sim\int{\rm d}t(-\frac{1}{2}(\partial_{t}{q})^{2}-\frac{1}{2}W_{q}^{2}-\overline{\psi}\partial_{t}\Psi-\overline{\Psi}W_{{q}{q}}\psi)$. One can reintroduce another auxiliary field $b$, and get the following twisted formulation for the $N=2$ supersymmetric action, ${I_{{\rm{twisted}}}}=\int{\rm d}t(\frac{1}{2}b^{2}+b(\partial_{t}{q}+W_{q})-\overline{\psi}(\partial_{t}\psi+W_{{q}{q}}\Psi))$, and the twisted nilpotent supersymmetry is $\delta_{\rm twisted}{q}=\psi,\delta_{\rm twisted}\psi=0,\delta_{\rm twisted}\overline{\psi}=b,\delta_{\rm twisted}b=0$, with ${I_{{\rm{twisted}}}}=\delta_{\rm twisted}\int{\rm d}t\overline{\psi}(\frac{1}{2}b+\partial_{t}{q}+W_{q})$. In this expression, the $N=2$ supersymmetry appears twisted, with some interested properties. As for its generalisation in the stochastic formulation, one can proceed as in the untwisted case, and, once one has computed the stochastic TQFT supersymmetry, one gets an interesting 4-simplex, because both supersymmetries of the stochastically quantised action can be described as nilpotent ones.

[TABLE]

From now on we will denote $I_{\small w}\equiv-I_{W}^{t\to it}$, keeping in mind that we compute Euclidean correlation functions in the $\tau=\infty$ limit. In order to achieve the stochastic quantisation, one uses the Euclidean action. It is unclear if the Wick rotation on $t$ can be done at finite values of the stochastic time. The possibility of a Wick rotation can be however checked in the limit $\tau=\infty$.

For polynomial interactions, it is convenient to separate the mass term $\frac{1}{2}m^{2}{q}^{2}$ from the rest of the interactions in $W(q)$, which are of degree higher than 2, with $W\equiv\frac{1}{2}m{q}^{2}+V({q})$ and $m=W_{{q}{q}}(0)$. This decomposition will be used shortly, to check perturbatively the possible stochastic equilibrium of the supersymmetric model.

We must define the supersymmetry covariance of all fields of the stochastic process. Consider the noise fields $\eta_{a}$. Given a component field $a$ of the multiplet $X$, its noise $\eta_{a}$ is a random fluctuations of $a$, modulo terms proportional to some equations of motion. Thus, if stochastic quantisation is compatible with supersymmetry, the noises of ${q},A,{\lambda},{\hat{\lambda}}$ build a superfield

[TABLE]

and transform accordingly. As we will demonstrate, non-trivial kernels $K_{ab}$ are necessary in order to obtain supersymmetry covariant Langevin equations as well as a proper convergence of the Langevin process. The Langevin equation that respects supersymmetry must be written as

[TABLE]

We will use superfield arguments to compute the appropriate kernels $K_{ab}$.

4 Determination of the kernels to enforce supersymmetry at all values of stochastic time

We will use the following remark : $X$ being a superfield, the supersymmetric invariance of the action ${I_{\small w}}(X)$ implies that its equations of motion build the following ”dual” superfield

[TABLE]

independently of the details of the supersymmetric action and thus of $W$. The proof is simple: if one varies $X\to X+\delta X$, the variation of ${I_{\small w}}$ is the supersymmetric term $\int dx(\frac{\delta{I_{\small w}}}{\delta{q}}\delta{q}+\frac{\delta{I_{\small w}}}{\delta{\lambda}}\delta{\lambda}+\frac{\delta{I_{\small w}}}{\delta{\hat{\lambda}}}\delta{\hat{\lambda}}+\frac{\delta{I_{\small w}}}{\delta A}\delta A)$. This variation is nothing but the integral over ${\rm d}\theta{\rm d}{\hat{\theta}}$ of the superfield $\delta X$ times $X^{*}_{W}$, which picks out the coefficients term of $\theta{\hat{\theta}}$ from the product $X^{*}_{W}\delta X$. Thus $X^{*}_{W}$ is a superfield, and one has:

[TABLE]

The stochastic time derivative of $X$ is also a superfield,

[TABLE]

We will need more independent superfields to define the stochastic BRST supersymmetry framework.

Let $a,b,...$ denote the components $x,{\lambda},{\hat{\lambda}},A$ of the superfield $X$. The goal is to write a supersymmetric action of the form ${s_{stoc}}{{\overline{\Psi}}_{a}}[K_{ab}b_{b}+\frac{{\rm d}Xa}{{\rm d}\tau}+K_{ab}\frac{\delta{S}}{\delta X_{b}}]$, where ${s_{stoc}}$ is the stochastic BRST symmetry operation we defined for Eq. (12). In the case we discuss, one has ${s_{stoc}}X_{a}=\Psi_{a},\ {s_{stoc}}\Psi_{a}=0,\ {s_{stoc}}\overline{\Psi}_{a}=b_{a},\ {s_{stoc}}b_{a}=0$. We will express the ${s_{stoc}}$-exact action depending on $X$ and its BRST topological partners as the integral of a super-Lagrangian.

Then, we will integrate this super-Lagrangian over superspace $(\theta,{\hat{\theta}})$ and get the supersymmetric ${s_{stoc}}$-exact action in components. Using Eq. (12), we will identify in a reversed way the Langevin equation for each component in $X$. Indeed, they are encoded in Eq. (12) as the coefficient of the terms linear in the $b$’s. Such Langevin equations will be covariant under supersymmetry, by construction.

To impose the supersymmetry transformation of the $b_{a}$’s, we simply declare that

[TABLE]

is a superfield. This determines the supersymmetry transformations of all components $b_{a}$.

Thus, the program is to first write a super-action depending on $B$, which has the form of the $b$-dependant terms in Eq.(12) after integration over superspace.

To get the stochastic ghost depending part of  Eq.(12), and enforce the compatibility of the supersymmetry of our model with the stochastic quantisation, we also introduce the stochastic antighost superfield $\overline{\Psi}_{X}$, with ${s_{stoc}}\overline{\Psi}_{X}=B$, ${s_{stoc}}B=0$,

[TABLE]

and the stochastic ghost superfield $\Psi_{X}$ of $X$,with ${s_{stoc}}X=\Psi_{X}$, ${s_{stoc}}\Psi_{X}=0$ and the decomposition

[TABLE]

For the sake of the definition of the kernels, we now introduce a mass parameter $M$. This parameter may have to be fine-tuned, to ensure the eventual convergence of the stochastic process, using the kernels that we shall shortly determine. The values of the correlators at $\tau=\infty$, if the limit exists, will not depend on the choice of $M$ according to the general proof of the $K$ independence of the limit.

One then defines the following supersymmetric quadratic functional of B :

[TABLE]

Whatever the value of $M$ is in this supersymmetric auxiliary free action, the equations of motion of $B$ determine the ”dual” superfield :

[TABLE]

Since $\dot{X},X^{*}_{W},B,B_{I_{M}}^{*}$ are all superfields, the following action is supersymmetric :

[TABLE]

It can be expressed in components as

[TABLE]

where $b_{x}^{\prime}=b_{q}-Mb_{A}$. This supersymmetric action has a quadratic dependance in the fields $b_{a}=(b_{q},b_{\lambda},b_{\hat{\lambda}},b_{A})$, which are the auxiliary fields of the topological BRST symmetry of stochastic quantisation, as in Eq. (12). The coefficients of the linear terms in the $b$’s provide the supersymmetry covariant Langevin equations

[TABLE]

The quadratic terms in $b$ define the Gaussian noise distribution by

[TABLE]

The general theorem applies and the limit of the field correlators in the limit $\tau\to\infty$ is independent on the chosen value of M.

To see perturbatively the $\tau\to\infty$ convergence of the solutions of Langevin equations, one expands $W_{qq}(q)={m}{}+V_{qq}(q)$ to verify that the drift forces are negative for all fields for $W=\frac{m}{2}{q}^{2}$ ‡‡‡The stochastic equation for $A$ is consistent if one eliminates $A$ by its equation of motion, $A=-m{q}$, so there is no evolution for $A$ provided $\eta_{A}=M\eta_{q}$. But the the various closure relations hold modulo equations of motion.

[TABLE]

Both operators $\partial_{t}^{2}-mM$ and $\partial_{t}^{2}-M^{2}$ have negative eigenvalues in the Fourier transform space. (This requires the Euclidian formulation for the stochastic quantisation). One must verify the the negativity of the eigenvalues of the differential operators that define the $\tau$ evolution of $A,{q}$ and ${\lambda},{\hat{\lambda}}$. They are

[TABLE]

The eigen values are non-degenerate for $M\neq m$ and degenerate for the simplest and easiest choice $M=m$. The determinant of the first matrix is $(\partial_{t}^{2}-M^{2})(\partial_{t}^{2}-m^{2})>0$ and for the second matrix it is $(\partial_{t}^{2}-mM)^{2}-(m-M)^{2}\partial_{t}^{2}>0$. Moreover the traces of both matrices are negative. Thus, whatever the value of $M$ is, we have the non positivity requirement and a normalisable vacuum for the supersymmetric Fokker–Planck process, and the stochastic process is converging at $\tau\to\infty$ for the fields ${q},A,{\lambda},{\hat{\lambda}}$. The values of correlators at $\tau\to\infty$ is independent on the choice of $M$, although the details of the evolution are $M$ dependent. The solutions have a different dependence on $\tau$ for $M\neq m$ and $M=m$, but the limit $M\to m$ is continuous.

We have thus obtained a set of Langevin equations with a well defined convergence at $\tau=\infty$ that are by construction covariant under the $N=2$ supersymmetry. For $M=m$, the stochastic process is simplest because no mixing occurs between the drift forces of all fields at the free level, and the Fokker–Planck Hamiltonian only involves the eigenvalues of the single operator $-\partial_{t}^{2}+M^{2}$ for all fields.

We can now write the complete ${s_{stoc}}$-exact action associated to these supersymmetric Langevin equations, that is, a stochastic BRST-TQFT supersymmetric action in the 2d-space $\{t,\tau\}$, which is the integral of a superfield, and thus a $N=2$ supersymmetric 2d action.

We use the already anticipated stochastic topological superghosts $\Psi_{X}$ and $\overline{\Psi}_{X}$, upon which the graded differential operator ${s_{stoc}}$ acts as ${s_{stoc}}X=\Psi_{X},\ {s_{stoc}}\Psi_{X}=0,\ {s_{stoc}}\overline{\Psi}_{X}=B,\ {s_{stoc}}B=0$. Defining,

[TABLE]

one has ${s_{stoc}}{\overline{\Psi}}_{I_{M}}^{*}=B_{I_{M}}^{*}$ and ${s_{stoc}}B_{I_{M}}^{*}=0$. The complete N=2 supersymmetric action that is invariant under the stochastic BRST supersymmetry and expresses the Langevin equations (53) is thus

[TABLE]

It can be simply expanded in components. Then, all the details of the stochastic evolution of the correlators of all components of the superfield can be studied while maintaining the supersymmetry.

5 Conclusion

Using a superfield construction, we explained the way the stochastic time evolution of the stochastically quantised $N=2$ supersymmetric quantum mechanics correlation functions is compatible with the $N=2$ supersymmetry. The stochastic TQFT supersymmetry and the space supersymmetry satisfy

[TABLE]

for every finite value of the stochastic time, which is the required property. We introduced kernels that ensure the convergence of the stochastic process toward the standard path integral of $N=2$ supersymmetric quantum mechanics. These kernels are covariant under supersymmetry and depend on a massive parameter that can be chosen at will. Its value modifies the course of the stochastic quantisation, but the limit of the correlation functions when the stochastic time runs to $\tau=\infty$ is independent on its choice.

The same demonstration can be mutatis mutandi repeated for the Wess and Zumino model. It is also possible to repeat this construction to more sophisticated supersymmetric theories, including those with a gauge invariance, provided one uses the method of equivariant cohomology.

Acknowledgment. I wish to thank the NCTS in Hsinchu as well as the Max Planck Institute in Golm for their generous and warm hospitality.