Effective Potential for Complex Langevin Equations

Effective Potential for Complex Langevin Equations

Gerald Guralnik gerry@het.brown.edu    Cengiz Pehlevan cengiz@het.brown.edu Department of Physics, Brown University, Box 1843, Providence, RI 02912
July 26, 2019
Abstract

We construct an effective potential for the complex Langevin equation on a lattice. We show that the minimum of this effective potential gives the space-time and Langevin time average of the complex Langevin field. The loop expansion of the effective potential is matched with the derivative expansion of the associated Schwinger-Dyson equation to predict the stationary distribution to which the complex Langevin equation converges.

Complex Langevin Equation, Schwinger-Dyson Equations
pacs:
11.15.Ha,02.50.Ey,05.10.Gg
preprint: Brown-HET-1571

I Introduction

Sampling complex integral weights of the form is important for many applications in high energy physics. Complex Langevin equations have been proposed for this purpose Parisi83 (); Klauder83 (). Many theoretical questions are still to be answered for complex Langevin equations. Most importantly, it is not known a priori if a complex Langevin simulation will converge to a stationary distribution. Even if it does, there are multiple stationary distributions that satisfy stationary limit equations (Fokker-Planck equation Salcedo93 () or Schwinger-Dyson equation Pehlevan07 ()) and it is not known how to identify to which one of these that the complex Langevin equation converges. See Pehlevan07 () for a discussion of these issues and references to literature.

Identification of the particular stationary distribution sampled by the complex Langevin equation needs input from the dynamics of the complex Langevin simulation. Analysis depending only on the quantities in the stationary limit has not been successful in differentiating between stationary distributions. To account for this fact, we introduce an effective potential approach to complex Langevin equations. Even though this effective potential does not have all the properties of the conventional quantum field theory effective action, we show that it governs the probability that the spacetime and Langevin time average of the fields takes on specific values. In the infinite Langevin time limit, it is shown that the spacetime and Langevin time average of the field must be equal to the minimum of the effective potential. As in quantum field theory, a loop expansion for this effective action can be defined. The parameter that defines loops is a factor that multiplies the noise in the complex Langevin equation. The same parameter is shown to count derivatives in the Schwinger-Dyson equation derived from the action . This point is important, because it is known that the stationary distributions of the complex Langevin equation are complexified path integral solutions of the Schwinger-Dyson equation, i.e. see Pehlevan07 (). The main idea of this paper is to use the minimum of the effective action to identify the stationary distribution to which the complex Langevin equation converges by matching the loop expansion with the derivative expansion. This provides an approach to addressing the multiple stationary distributions problem.

In two recent studies Aarts08 (); Berges0708 (), the classical fixed point structure of complex Langevin equations were shown to be important for the convergence behavior. Our results will also parallel this result since to leading order the minima of the effective potential will be given by classical fixed points. The next order will be used to remove the degeneracy and this will lead us to the identification of the stationary distribution that is sampled by the complex Langevin simulation.

We start by reviewing relevant results in Section II. We then introduce the effective potential in a lattice setting, discuss its properties and calculate it to one loop order (Section III). Before we conclude, we give a simple application of these ideas (Section IV).

One point we omit in this paper is the discussion of the continuum limit. The question of renormalizability is a nontrivial problem and deserves a separate study of its own.

Ii Stationary Distributions of Complex Langevin Equations

In this section we review the relevant results of Pehlevan07 (). Even though, our ultimate goal in this paper is to do calculations in a lattice setting, the arguments in this section are given in the continuum for ease of notation. The lattice versions of these results are either in Pehlevan07 () or Berges07 ().

Given a complex action of a scalar field in a Euclidean space,

(1)

we want to study the stochastic dynamics given by

(2)

where and are independent Gaussian noises with correlators given by

(3)

and

(4)

and are real constants and . corresponds to the complex Langevin equation. We will not set to any particular value for the moment. will define orders in the approximations we do below. Equation (2) can be written as a stochastic system of two real fields, . In fact, this is what we will do in the rest of this paper.

Note that most of the time in the literature is set to 1. There is only a real noise in the complex Langevin equation. As shown in Pehlevan07 (), this may lead to initial condition dependence of the stationary distribution sampled by the complex Langevin equation. This happens because for some actions the complex Langevin dynamics may be confined to sampling in lower or upper half planes depending on the initial condition, not having access to all of the complex plane. In this paper, we want to avoid such initial condition dependence and assume that the stationary distribution sampled by the complex Langevin equation is unique. Furthermore, setting to some value larger then 1 will be important for deriving a path integral representation of the generating functional for the complex Langevin process in the next section.

Given some initial conditions for the field, Langevin equation (2) has a time dependent probability distribution, see e.g Justin02 (),

(5)

which can be shown to satisfy a Fokker-Planck equation,

(6)

The existence and uniqueness of stationary distributions for this Fokker-Planck equation is still an unanswered question.

We continue by defining a generating functional for equal time correlation functions

(7)

where is a real source field. Here the brackets denote an average over the noise. Field variables are considered as functionals of the noise through equation (2), hence we include a superscript . The effect of noise appears in the stochastic average through the probability distribution (see e.g. Justin02 ()),

(8)

This is now an integral over all field configurations weighted by . This point will be treated in more detail in the next section. Using the Fokker-Planck equation (II), one can show that the following equation holds

(9)

If there is a stationary distribution, the left hand side is zero. When , functional derivatives of the right hand side with respect to produces Schwinger-Dyson identities for the action , since the integrand must vanish due to the arbitrariness of . plays the role of in quantum field theory. It is easy to see this by defining a making a change of variables .111In terms of this new variable, the correlation functions are given by functional derivatives of with respect to . Then, the stationary distribution condition becomes

(10)

counts the number of functional derivatives, as would. As is well known, see Bender77 (), a derivative expansion of Schwinger-Dyson equations is equivalent to a loop expansion of the associated path integral around a local minimum.

If one considers zero dimensional theories, an explicit construction of the stationary equal time correlation functional can be given. Generalization to a lattice is straightforward, see Pehlevan07 (). Now, the action is a function of . Equation (9) becomes

(11)

In particular, the generating function for the stationary distribution satisfies

(12)

We can solve this equation following Garcia96 (). Our ansatz is

(13)

where is a contour over the complex plane. Inserting this into equation (12) one gets:

(14)

This equation can be solved for

(15)

and is contour that connects the zeros of on the complex plane. As advertised, plays the role of . Since equation (12) is linear, any linear combination of the (independent) solutions will also be a solution,

(16)

where define an independent subset of paths .222Essentially the same result is derived by another method in Salcedo93 (). The only constraint on the coefficients at this point is the normalization condition for correct distribution interpretation, namely that . We restricted the form of the generating function, but did not identify it completely. Examples are given in Pehlevan07 ().

In this paper, we use an effective action approach to identify the stationary distribution. The expansion gives a small noise expansion for the complex Langevin equation and a loop expansion of the stationary distribution generating function. Matching of both series order by order will tell us the exact form of the stationary distribution.

Iii Effective Action for Complex Langevin Equation

In this section, we will define the effective action for complex Langevin equation and do an expansion to gather information about the stationary distribution. We choose to use a lattice regularization. We work with an Ito type discretization for equation (2).

First we define our notation. We introduce a Euclidean spacetime lattice with isotropic lattice spacing . Lattice points are denoted by a dimensional vector , where . The coordinate of a lattice point is given by . Our convention is to denote Euclidean time with . Periodic boundary conditions are assumed for the space-time lattice. We discretize Langevin time with intervals. Langevin simulation is assumed to start at and run for time steps untill . In terms of dimensionless lattice variables , , , , , where the noise correlators are now

(17)

the discretized equation (2) in Ito calculus reads

(18)

In the rest of this paper we assume . To avoid confusions, we note that is the lattice variable at the spacetime point and Langevin time .Here, the dimensionless lattice Laplacian is defined by

(19)

where is a unit vector pointing along the -direction. We assume an initial condition is set at by specifying the field on the spacetime lattice.

iii.1 Generating Functional for Dynamic Correlation Functions

First step in our anaylsis is to define a generating functional for dynamic correlation functions,

(20)

is a functional of the noise fields through the Langevin equation (18). The sum in parameter starts from , since we assume an initial condition is given. Unless explicitly stated, sums and products run from to . We use subscript to differentiate from the generating functional of equal time correlation functions defined by equation (7). is more general than since it generates correlation functions of the field variable at different times as well as equal times. We will define an effective action using this generating functional. Before doing that, we derive a path integral representation for . Although we will not follow their work directly, see Kawara90 () for a discussion of path integral representations of generating functionals based on Ito discretization and Gozzi83 () for continuum formulation. By definition, is an average over the probability density functions of the Gaussian random variables,

(21)

is a proper normalization constant for each equation it appears in. The next step is to rewrite this expression as an integral over field configurations. Using equation (18), we can change integration variables from the Gaussian random variables and to and . Note that due to our discretization, the Jacobian of this transformation is a constant. Defining

(22)

the generating functional becomes

(23)

where the measure is

(24)

There is no integration over , since we assume it to be set by initial conditions. The exponent in the second line, with the normalization factor , is the probability density of a particular realization of the Langevin process given some initial . We will denote it by . The previously defined probability distribution (5) and are very closely related, see e.g. Namiki92 (). For , is an overall factor for the exponent, which will allow us to define a loop expansion for this generating functional.

For later convenience, we define the quantity in the exponent of the second line of equation (III.1) as

(25)

Forgetting the complex Langevin equation roots for one moment, one can think of equation (III.1) as a the generating functional of a field theory on a dimensional lattice with an action given by equation (25). This is not a usual field theory, because it does not have translational invariance due to terms appearing in equation (25). One end of the lattice () is fixed, whereas the other end () is free. This picture will be useful later when we discuss the effective potential.

Before continuing, we note that having was a great convenience for this calculation. Otherwise, we would start with averages over only the real noise variables and this would lead to a difficulty when changing variables in the path integral from noise to field variables. The number of field variables would be half of . Therefore, one would need to choose half of the field variables as independent ones and express the remaining in terms of the independent ones. Having saved us from these difficulties.

iii.2 Effective Action, Effective Potential and Their Interpretation

Now we do the standart quantum field theoretic procedure to define an effective action for this generating functional,

(26)

where

(27)

By definition, it follows that

(28)

For field configurations that are homogenous and static in Langevin time, the effective action reduces to “effective potential”,

(29)

In the next section, we will calculate the effective potential to one loop order (first order in ).

We note that the effective action defined above does not have all the properties one remembers from quantum field theory. In particular, it does not have the energy interpretation Coleman88 () because the complex Langevin equation (2) as a system of two real Langevin equations is not derivable from an action.333A Langevin system that is derived from an action will have component equations . See Hochberg99 () for a detailed discussion. However, one useful property is still valid. The effective potential governs the probability that the spacetime and Langevin time average of the field takes specific values. To show this, we follow closely the argument of Hochberg99 ().

We ask the probability of spacetime and Langevin time average of the field variables taking a specific value,

(30)

At this point we will make an approximation, though we will argue that our approximation becomes exact in the infinite Langevin time limit. We are looking for a simplification of the effective action evaluated in the presence of a constant source field. For this purpose, we consider a derivative expansion of the effective action. The first term will be the effective potential, as defined above. Other terms will have factors of (lattice) derivatives of the field variable . Note that we are now regarding the field variable as the expactation value of a “quantum” operator of the lattice field theory given by (III.1). Consider complex Langevin initial conditions that have no space-time dependence. Then, for a constant source field, one expects the derivatives along the space-time directions to vanish due to periodic boundary conditions and translation invariance. On the Langevin time direction translation invariance breaks, because of the existence of an initial condition. Even though this is the case, we will argue that in the long Langevin time limit, one can ignore higher order terms and safely approximate the effective action by the effective potential. Consider for a moment a complex Langevin simulation starting with space-time independent initial conditions. Assuming convergence to a stationary distribution, in the long Langevin time limit the field expectation values will “forget” the initials conditions and approach a stationary value. This means that the Langevin time derivatives of the field expectation values will be small, in fact will vanish in the long time limit. By definition, one would be able to calculate this effect using the generating functional given by (III.1) by setting the source fields to zero. We are interested in constant source fields that are not vanishing. Now let’s go back to the Euclidean lattice field theory picture defined by (III.1). For zero source fields, in the very large region the lattice looks pretty much translationally invariant as argued above. Adding a constant source field will not break this invariance. Therefore, vanishing of the Langevin time derivatives must still hold. We argue that if this is the case, the effective potential will dominate over all higher order terms in the derivative expansion. To motivate this last comment a little more rigorously, consider a constant (like the effective potential) and a function (like the higher order terms in the derivative expansion) that dies as . Now consider the quantity

(31)

We will show that this limit is zero, which will prove that our comment is valid, that the effective potential dominates over higher terms in the derivative expansion. To prove that it is zero, we will show that it is smaller than any positive . We will need to recall the definition of a convergent sequence. A sequence converges to the limit if, for any , there exists a such that for . Now we apply this definition to sequence , such that . Now there exists a such that for . Then

(32)

Therefore, we can safely assume that a constant source field implies constant . As discussed above, this constant is the expectation value of the field variable in the long Langevin time limit in the presence of a constant source field. Then

(33)

Here the superscpript is used to denote that fields are functions of the constant source field. We will look at infinite Langevin time limit. Then using the method of stationary phase

(34)

Thus, the effective potential governs the probability distribution of the spacetime and Langevin time average of the field. As the simulation time goes to infinity, the space-time and Langevin time average of the fields must equal to the minimum of the effective potential444We assume no degeneracy among the minima.. Our argument for effective potential dominance assumed space-time indenpendent initial conditions. However, the stationary distribution that the complex Langevin equation converges should be independent of the initial condition. Therefore this result can be generalized to any initial conditions. Note that if our argument on the dominance of the effective potential is valid, this is an exact result in the infinite limit.

This interpretation of the effective potential will enable us to determine the stationary distribution to which the complex Langevin equation converges.

iii.3 Effective Potential in One Loop

In this section, we will calculate effective potential to first order in in the infinite Langevin time limit. As we showed in the previous section, in this limit the space-time and Langevin time average of the field must equal to the minimum of the effective potential.

Calculation of the effective potential for a complex Langevin simulation that runs for a finite time is not easy. The calculation requires, as we will see, the evaluation of a -by- matrix determinant and to the authors, there does not seem to be a way of simplifying this calculation. Instead we will do our calculation in the infinite Langevin time limit, but we will take this limit in a different way than what we have been considering until now. Until now, we let the complex Langevin simulation to start at with a fixed initial condition and run until a finite time . Then we took the limit . Instead, we now start our simulation at time with a fixed initial condition and run it until . Then we let go to infinity. In the large limit, due to ergodicity and Langevin time independence of the complex Langevin equation, the Langevin time averages of these two processes will converge to the same stationary distribution. Another way to see the equivalence is by shifting the latter process by steps in Langevin time and thus reaching a process of the former type. Note that this changing of limiting process does not effect the Schwinger-Dyson equation related arguments. Since stationary distributions sampled by the complex Langevin equation for both processes are the same, we can calculate the effective potential for the latter process, and its minimum will match the minimum of the former one. Matching of minima is sufficient for our purposes, but it is very likely true that the effective actions themselves match also. From now on, we work in the limit where the limit is taken in the latter way.

We start by expanding the integrand of the generating functional for dynamical correlation functions around constant field configurations and . First, we choose the constant source fields as

(35)

Using these values of the source fields, we get

(36)

Here the sums go from to . The matrix is given by

(37)

where

(38)

When deriving these equations, we made use of the Cauchy-Riemann equations. They are evaluated at the constant field configurations and . Now, using the definition of the effective potential, it is easy to see that

(39)

Before continuing our calculation, we stop to discuss this result. As expected, the leading order term sets the minima of the effective potential at the classical fixed points, which are saddle points of the action . This is just saying that if there were no noise (and the simulation did not diverge), the fields will evolve to a fixed point and stay there for ever. Averages over space-time and Langevin time in the infinite Langevin time limit will be set by the classical fixed points. At the classical level, all saddle points of are degenerate, meaning they are all minima of the effective potential. Adding noise will remove the degeneracy through higher terms in the expansion. We will calculate the determinant term below and show in a specific example how it removes the degeneracy.

On a different line of argument, we showed above that if the complex Langevin equation converges to a stationary distribution, this distribution will be a solution of the Schwinger-Dyson equation for the theory with complex action . Moreover, we showed that counts the number of derivatives (loops) for this theory. Now we merge this result with the discussion of the previous paragraph as follows:

Assuming ergodicity, averages over the stationary distribution is equivalent to averages over Langevin time in the infinite Langevin time limit. Since we set periodic boundary conditions on the space-time lattice, the resulting stationary distribution will have translation invariance in space-time. In particular, the Langevin time averages for the field variable will be constant in space-time. Therefore, we can say that Langevin time average of the field variable is equal to the Langevin time and space-time average of . Based on our results above, we know that the latter is equal to the minimum of the effective action for the complex Langevin equation which is given in equation (39) to order in . Since simultaneously counts loop order for the stationary distribution, which is a solution to the Schwinger-Dyson equations, we conclude that the complex Langevin equation converges to a stationary distribution that is quantized around a saddle point of the complex action . This saddle point is given by the minimum of the effective action of the complex Langevin equation.

The last statement needs more explanation. First of all, by quantization of a complex action we mean not the Feynman path integral quantization but a Schwinger action principle quantization. From Schwinger action principle, one derives Schwinger-Dyson equations, which have multiple solutions Guralnik07 (); Garcia96 ().555Other differential equations also follow from Schwinger action principle, but it was shown in Guralnik07 () that a study of Schwinger-Dyson equations is sufficient since in the continuum limit solutions of Schwinger-Dyson equations satisfy the remaining differential equations. A solution of Schwinger-Dyson equations defines a generating functional for a quantization of . This is consistent with equation (16), where we showed that in zero-dimensions solutions of Schwinger-Dyson equations can be represented by linear combinations of complexified path integrals. This result trivially generalizes to a lattice Guralnik07 (); Pehlevan07 ().666To our knowledge generalizations of generating functionals with complex integration paths are not yet constructed on the continuum. Complex Langevin equations can be of help in defining these continuum generating functionals. Complex Langevin simulations have only one stationary distribution, the question is to find which quantization of this distribution gives. Our argument shows that the stationary distribution will contain contributions only from one of the saddle points of .777Here we assumed that the order term of the equation (39) breaks the degeneracy among the minima completely. If not, one may need to go to higher orders. If the degeneracy is not broken at all, it may be that the stationary distribution contains contributions from all degenerate minima. The generating functional for this distribution may be defined by a linear combination of the generating functionals that contain only contributions from one of the degenerate minima. See Guralnik07 (). One way of making this statement precise is through the results of Guralnik07 (). There it was shown that for zero-dimensional theories, loop expansions of Schwinger-Dyson equations can be Borel resummed to give path integrals of type (13) with the path being equivalent to a steepest descent path passing through a dominant saddle point of . It is these generating functionals for zero-dimensions, and their generalizations to higher dimensions that we argue to be the stationary distributions of the complex Langevin equations. Examples of this will be given in the next section.

Now we go back to the determinant calculation. Calculation of the determinant of matrix is most easily done in the Fourier basis. The discrete time Fourier transform and the inverse transform are given by

(40)

The nonzero elements of the matrix in this new basis are: