## Real-Time Feynman Path Integral Realization of Instantons

Aleksey Cherman 1, ∗ and Mithat Unsal ¨ 2, †

1 Fine Theoretical Physics Institute,

Department of Physics, University of Minnesota, Minnesota, MN 55455, USA

2 Department of Physics, North Carolina State University, Raleigh, NC, 27695

In Euclidean path integrals, quantum mechanical tunneling amplitudes are associated with instanton configurations. We explain how tunneling amplitudes are encoded in real-time Feynman path integrals. The essential steps are borrowed from Picard-Lefschetz theory and resurgence theory.

Introduction -The real-time Feynman path integral[1] is one of the several equivalent formulations of quantum mechanics (QM) and quantum field theory (QFT). The amplitude of an event is found by summing over all path histories with the contribution of each path γ weighed by its action S γ [ ] via the factor e iS γ / [ ] glyph[planckover2pi1] . For real paths, this is a sum over pure phases. The real-time path integral formulation is very useful in many settings, and yields beautiful explanations of many quantum phenomena, such as the interference patterns seen in double-slit experiments and the Aharonov-Bohm effect. Yet surprisingly, some very basic quantum phenomena do not have a known description in the real-time path integral. For instance, tunneling phenomena, which are well-understood in the operator (Hamiltonian) formalism, as well as in the Euclidean-time path integral formulation, are not yet understood in the real-time Feynman path integral formulation[2].

We focus on a particle in a symmetric double-well potential as a paradigmatic example. The level splitting between the ground state and first excited state is given by a non-perturbative factor, ∼ exp[ -A/ g ( 2 glyph[planckover2pi1] )] where A is a pure number and g is a parameter controlling the depth of the wells. It is a consequence of tunneling in the operator formalism or instanton saddles[3] in the Euclidean path integral formulation, but its origin is mysterious in the Feynman path integral.

Understanding such non-perturbative phenomena directly in real-time path integrals is likely to be very important. Euclidean path integrals give deep insights into the non-perturbative properties of QM and QFT, but there is a huge range of problems for which they are not useful, such as calculations of scattering amplitudes or non-equilibrium observables. These questions are more naturally handled using real-time path integrals. Yet the current understanding of real-time path integrals beyond the perturbative level is extremely rudimentary, as is highlighted by the fact that even tunneling phenomena are not understood in the framework. Our results are a step toward developing a non-perturbative understanding of real-time path integrals.

Euclidean-time instantons -To set notation we review the standard path integral treatment of tunneling in QM. The real-time Feynman path integral representation of the amplitude A for a particle initially at position x i at time t i to end up at position x f at time t f is

$$\mathcal { A } _ { x _ { i }, t _ { i } } ^ { x _ { f }, t _ { i } } = \int _ { x ( t _ { i } ) = x _ { i } } ^ { x ( t _ { f } ) = x _ { f } } d [ x ( t ) ] \, e ^ { \frac { i } { h } S [ x ] } \, \quad \quad ( 1 )$$

where S is the classical action. For a non-relativistic bosonic particle moving in a symmetric double-well potential the action can be written as

$$S [ x ] = \frac { 1 } { 2 g ^ { 2 } } \int _ { t _ { i } } ^ { t _ { f } } d t \left [ ( \partial _ { t } x ) ^ { 2 } - ( x ^ { 2 } - 1 ) ^ { 2 } \right ] \quad ( 2 )$$

Our semiclassical treatment of tunneling will be justified if g glyph[lessmuch] 1. If the particle is initially near e.g. x = -1 at t = t i the amplitude for it to end up near x = +1 will be zero to any order in a perturbative expansion in g .

Let us now focus on the tunneling amplitude from x i = -1 at t i = -∞ to x f = +1 at t f = + ∞ . The standard way to compute this amplitude using path integrals involves moving to Euclidean time by a Wick rotation t →-iτ, x ( ) t → x τ ( ). The Euclidean action

$$S _ { E } = \frac { 1 } { 2 g ^ { 2 } } \int _ { - \infty } ^ { \infty } d \tau \left [ ( \partial _ { \tau } x ) ^ { 2 } + ( x ^ { 2 } - 1 ) ^ { 2 } \right ]. \quad ( 3 )$$

is then positive semi-definite and enters the path integral as e -S E . The Euclidean equation of motion ∂ x τ 2 τ ( ) = 2 x τ ( ) ( x τ ( ) 2 -1 ) has two finite-action solutions, the instanton x + I and anti-instanton x -I :

$$x _ { I } ^ { \pm } ( \tau ) = \pm \tanh ( \tau - \tau _ { 0 } ), \text{ \quad \ \ } ( 4 )$$

where τ 0 is a zero mode parameter. These solutions interpolate from x = ∓ 1 to x = ± 1 as τ goes from -∞ to + ∞ . In what follows we will focus on x + I and will drop the superscript. The Euclidean action of (4) is S I = 4 3 g 2 , so that (4) is a saddle-point of the Euclidean path integral, with contributions weighed by e -4 3 g 2 . For g glyph[lessmuch] 1 the desired tunneling amplitude can be computed using a steepest-descent approximation, and is given by

$$\mathcal { A } = \sqrt { \frac { 3 2 } { \pi g ^ { 2 } } } e ^ { - \frac { 4 } { 3 g ^ { 2 } h } }$$

where the prefactor is obtained by integrating over Gaussian and zero mode fluctuations around x = x I ( τ ).

Real-time instantons -The Euclidean-time derivation is not fully satisfying, as physics happens in real time.

FIG. 1. (Color Online.) The complexified instanton configuration (9). Left column: α = π/ 2 (Euclidean instanton). Center column: α = π/ 4. Right column: α = π/ 8. The top row shows the trajectories of (9) for τ ∈ -∞ ( , + ∞ ), while the bottom row shows the real and imaginary parts of (9) as a function of τ .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0012v2-figure-3111529a-b2e1-4a5a-a76a-8cd687e13fcb.jpg)

Yet there is a good reason that instantons were never formulated in the real-time Feynman path integral. A natural attempt to do this involves looking for finiteaction solutions of the real-time equation of motion

$$- \partial _ { t } ^ { 2 } x ( t ) = 2 x ( t ) \left ( x ( t ) ^ { 2 } - 1 \right ). \quad \quad ( 6 ) \quad \text{$\mathfrak { a }$}$$

Since we are looking for an answer of the form e iS/ glyph[planckover2pi1] = e -A g 2 glyph[planckover2pi1] the action of the solution would have to be imaginary. That means that a 'real-time instanton' configuration has to be complex. Such a notion raises major conceptual issues:

- (A) If the field space is complexified, how do we avoid an undesired doubling of degrees of freedom in the complexified path integral?
- (B) Relatedly, how do we avoid including pathological directions in the complex field space which may yield exponentially large contributions? Over what sub-space of fields (complex paths) should one integrate over?

Setting these important questions aside for the moment, we note that it is tempting to guess that the desired saddle field configurations might be given by the Wick rotation of the Euclidean instanton back to real time. Indeed, (6) does have the formal solution

$$x ( t ) = i \tan ( t - t _ { 0 } ) \quad \quad \ ( 7 ) \quad \text{an}$$

which is a Wick rotation of (4). But (7) looks nonsensical, for many reasons:

(i) It is pure-imaginary-valued, so it cannot possibly give a contribution like e iS/ glyph[planckover2pi1] = e -A g 2 glyph[planckover2pi1] , and cannot describe a path interpolating from x = -1 to +1.

(ii) It does not appear to have anything to do with tunneling from x = -1 to +1, since it is never localized near the bottoms of the wells.

(iii) It has singularities at t = (2 n +1) π/ , n 2 ∈ Z , and these singularities are not integrable, so that the action is not well-defined.

Until fairly recently this is where an analysis would have had to stop. Fortunately, there has recently been major progress in the understanding the complexification of functional integrals [4-6] and the structure of the semiclassical expansion of path integrals in QM [7] and QFT [8] via resurgence theory[9]. Motivated by these developments, let us see what happens if we do not insist on working in either purely real or imaginary time.

Specifically, let us start in real time t and do a complex Wick rotation t → τe -iα , φ ( ) t → φ τ ( ). In contrast to purely real or purely imaginary time, now the action itself becomes complex:

$$\stackrel { \dots } { \text{of} } \\ \alpha = 0$$

Hence issues (A) and (B) with the complexification of configuration space in path integrals must be deal with, and we do so below. For small α , τ ≈ t and S α ≈ S , while for α = π/ 2 τ becomes the standard Euclidean time coordinate and iS π/ 2 = -S E . It turns out that an infinitesimal α acts as a crucial regulator which is necessary to make sense of real-time instantons.

The α -dependent equation of motion -e 2 iα ∂ x τ 2 τ ( ) = 2 x τ ( ) ( x τ ( ) 2 -1 ) has the solutions

$$x _ { I } ( \tau ) = \pm \tanh \left [ ( \tau - \tau _ { 0 } ) e ^ { - i ( \alpha - \pi / 2 ) } \right ] \quad ( 9 )$$

To make (9) a legitimate solution one must complexify the space of path histories from x : R → R to

FIG. 2. (Color Online.) Upper-left: trajectory of (9) for α = 4 × 10 -2 × π 2 in complex configuration space. Upper-right: real and imaginary parts of (9) for α = 1 × 10 -2 × π 2 , with a green α = π/ 2 trajectory shown for comparison. Lower-left: illustration of the manner in which the imaginary part of (9) approaches the singular configuration (7) for small α (here α = 1 × 10 -2 × π 2 ), while the real part approaches to a distribution, a signed Dirac comb. Lower-right: the trajectory for α = 1 × 10 -3 × π 2 , with red dots marking the centers of the wells x = ± 1. In the real-time limit the real-time instanton appears to become a space-filling curve in complexified configuration space.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0012v2-figure-9fd3d403-5a06-4362-a47f-e1cd2f97d39c.jpg)

x : R e -iα → C . While (9) looks like a simple analytic continuation of (4), the conceptual issues raised by allowing such an continuation are quite complex.

Naively, precisely at α = 0 (9) reduces to (7), although we are about to see that the limit is subtle. Consider the situation for 0 &lt; α ≤ π/ 2 [10]. Remarkably, the imaginary parts of the complex instanton solution (9) vanish for large | τ | , and x I approaches ∓ 1 as τ →∓∞ ! The behavior of the complexified instanton at generic α values is summarized in Fig. 1 and Fig. 2. For generic 0 &lt; α ≤ π/ 2 field starts at x = -1 at τ = -∞ , whips around in the left half of complexx plane for τ &lt; τ 0 , then crosses to the right half of the plane at τ = τ 0 , whips around x = +1 for τ &gt; τ 0 , and then settles down at x = +1 at τ = + ∞ .

illustrated in Fig. 2. The Re x &lt; 0 and Re x &gt; 0 halfplanes form two basins of attraction for the repulsive and attractive fixed points x = -1 and x = +1 respectively. After starting at x = -1 at τ = -∞ , as α → 0 the instanton spends more and more time whipping around in the Re x &lt; 0 region in a huge number of more and more tightly spaced increasingly large arcs for τ &lt; τ 0 , then jumps over to the Re x &gt; 0 region and winds closer and closer to x = +1 in the same manner, until settling at x = +1 at τ = + ∞ . From the bottom-right plot in Fig. 2, the complex instanton trajectory appears to become a space-filling curve in the real-time limit.

glyph[negationslash]

The instanton remains smooth for any α = 0, but the behavior in the real time limit α → 0 is singular, as

In the strict α = 0 limit, the imaginary part of the instanton solution is given by (7), while the real part

approaches to a distribution, given by

$$\text{Re} \, x = \pi \text{sign} ( t - t _ { 0 } ) \sum _ { k = - \infty } ^ { \infty } \delta ( t - t _ { 0 } - \frac { \pi } { 2 } ( 2 k - 1 ) ) \ \ ( 1 0 ) \quad \text{integr} \quad \text{cls de} \quad \text{that} \ \ ;$$

which is a signed Dirac comb. The time average of 〈 Re x 〉 T and 〈 Im x 〉 T (which we define as 〈·〉 T = 1 T ∫ T τ 0 · dτ with T = πk + τ 0 , k ∈ Z ) for any α , including α = 0, essentially follows the pattern of a Euclidean instanton,

$$\langle \text{Re} \, x \rangle _ { T } & = \begin{cases} \, + 1, \quad & \tau > \tau _ { 0 } & \quad & \text{th} \\ \, - 1, \quad & \tau < \tau _ { 0 } & \quad & \text{the} \end{cases} \\ \langle \text{Im} \, x \rangle _ { T } & = 0 & \quad & \quad & \text{(11)}$$

as illustrated by Fig. 2.

Next, we show that the action -iS α of the instanton is always precisely 4 / (3 g 2 ) for any α , with the action at α = 0 defined by continuation from | α | &gt; 0. The Lagrangian evaluated on (9) oscillates wildly through the complex plane as α → 0. Nevertheless, -iS α [ x I ] turns out to be independent of α . To show this we rewrite

$$S _ { \alpha } ( x ) = \frac { e ^ { - i \alpha } } { 2 g ^ { 2 } } \int _ { - \infty } ^ { + \infty } d \tau \left ( \left [ e ^ { + i \alpha } \dot { x } - W ^ { \prime } ( x ) \right ] ^ { 2 } + 2 e ^ { + i \alpha } \dot { x } W ^ { \prime } \right ) \int _ { \frac { \partial \ell } { A } } ^ { \nu \omega \nu }$$

where ′ ≡ ∂ , x ˙ ≡ ∂ τ , and ( W ′ ) 2 = V ( x ). But (9) satisfies the first-order equation e + iα ˙ x -W x ′ ( ) = 0, so

$$- i S ( x _ { I } ) = \frac { 1 } { g ^ { 2 } } \int _ { - \infty } ^ { + \infty } d \tau \dot { x } _ { I } W ^ { \prime } ( x _ { I } ) = \frac { 4 } { 3 g ^ { 2 } } \quad ( 1 2 ) \quad \text{ simple}$$

We conjecture that the solution with minimal non-zero action is (9) for any α . This is obvious at α = π/ 2, while at α = 0, resurgence theory implies that the minimumaction saddle must have action -iS x ( I ) due to the largeorder structure of perturbation theory. We leave a detailed exploration of this conjecture to future work.

The question of whether singular but finite-action field configurations contribute to path integrals is controversial[11]. Our results support the notion that such configurations do contribute, since (9) is a finite action configuration whichis associated with a dramatic physical effect of level splitting. Hence it is physically clear that it must be included it in the evaluation of the path integral, even though it becomes singular and discontinuous in the strict α = 0 limit.

Lefshetz thimbles -To find the contribution of the complex instantons to observables in the real time limit, it is necessary to take into account fluctuations around the instanton saddle points. As we already noted this is a subtle problem due to the complexification of the field space, and involves understanding the analytic continuation of path integrals. A program to study this issue was recently initiated in [5].

A sensible complexification of the integration cycle of an integral which initially has a real domain of integration must not double the number of degrees of freedom. A full understanding of the analytic continuation of the integral thus involves finding a basis of convergent middle-dimensional integration cycles in the complexified integral, and then working out how the original real cycle decomposes into the complex cycles. Ref. [5] showed that in favorable circumstances, this can be done by using Picard-Lefshetz theory[12]. Every finite-action saddle point of the path integral is associated to a convergent integration cycle called a 'Lefshetz thimble' - a multidimensional generalization of a steepest-descent path. The thimbles are solutions of the downward flow equations for the Morse function Re [ iS ]. For us, the flow equation is

$$\frac { \partial x } { \partial s } = - \frac { e ^ { + i \left ( \alpha - \pi / 2 \right ) } } { \overline { g ^ { 2 } } } \left [ e ^ { - 2 i \alpha } \partial _ { \tau } ^ { 2 } \bar { x } + 2 \bar { x } \left ( \bar { x } ^ { 2 } - 1 \right ) \right ] \quad ( 1 3 )$$

where x is complex and s is the flow 'time' along the thimble. The fluctuations around x I are described by the solutions of (13) with the initial condition x s ( →-∞ , τ ) = x I ( τ ). The structure of the Gaussian fluctuations can be determined by writing x s, τ ( ) = x I ( τ ) + δx s, τ ( ) with δx s ( →-∞ , τ ) = 0 and expanding (13) to O ( δx ):

$$\Big | _ { \infty } \Big | _ { \infty } \Big | _ { \infty } \Big | _ { \infty } \Big | _ { \infty } \Big | _ { \infty } \Big | _ { \infty } \Big | _ { \infty } \Big | _ { \infty } \Big | _ { \infty } \Big | _ { \infty } \Big | _ { \infty } \Big | _ { \infty } \Big | _ { \infty } \Big | _ { \infty } \Big | _ { \infty } \Big | _ { \infty } \Big | _ { \infty } \Big | _ { \infty } \Big | _ { \infty } \Big | _ { \infty } \Big | _ { \infty } \Big | _ { \infty } \Big \| _ { \infty } \Big | _ { \infty } \Big | _ { \infty } \Big | _ { \infty } \Big | _ { \infty } \Big | _ { \infty } \Big | _ { \infty } \Big | _ { \infty } \Big | _ { \infty } \Big | _ { \infty } \Big | _ { \infty } \Big | _ { \infty } \Big | _ { \infty } \Big | _ { \infty } \Big | _ { \infty } \Big | _ { \infty } \Big | _ { \infty } \Big | _ { \infty } \Big | _ { \infty } \Big | _ { \infty } \Big | _ { \infty } \Big | _ { \infty } \Big | _{ \infty } \Big | _ { \infty } \Big | _ { \infty } \Big | _ { \infty } \Big | _ { \infty } \Big | _ { \infty } \Big | _ { \infty } \Big | _ { \infty } \Big | _ { \infty } \Big | _ { \infty } \Big | _ { \infty } \Big | _ { \infty } \Big | _ { \infty } \Big | _ { \infty } \Big | _ { \infty } \Big | _ { \infty } \Big | _ { \infty } \Big | _ { \infty } \Big | _ { \infty } \Big | _ { \infty } \Big | _ { \infty } \Big | _ { \infty } \Big | _ { \outline } \Big | _ { \infty } \Big | _ { \infty } \Big | _ { \infty } \Big | _ { \infty } \Big | _ { \infty } \Big | _ { \infty } \Big | _ { \infty } \Big | _ { \infty } \Big | _ { \infty } \Big | _ { \infty } \Big | _ { \infty } \Big | _ { \infty } \Big | _ { \infty } \Big | _ { \infty } \Big | _ { \infty } \Big | _ { \infty } \Big | _ { \infty } \Big | _ { \infty } \Big | _ { \infty } \Big | _ { \infty } \Big | _ { \infty } \Big | _ { \infty }$$

The existence of the flow equations means that it is possible to choose sensible integration cycles in the complexified path integral, and in the Gaussian-fluctuation limit the determination of our desired cycle reduces to finding the solution of (14). By moving to the variable τ ′ = τ exp[ -i ( α -π/ 2)] and rescaling δx one can map (14) to the Euclidean flow equation. This implies that integrating over the thimble amounts to calculating the standard fluctuation determinant (see [6] for a nice discussion of this) multiplying e -4 / g 3 2 , and reproduces (5).

Conclusions -We have shown how to compute tunneling amplitudes between low-lying states, a hallmark non-perturbative feature of quantum mechanics, directly from a real-time Feynman path integral. The amplitude is given by a real-time instanton saddle point, which requires a complexification of the space of paths and the addition an infinitesimal imaginary part to the time coordinate as a regulator for its definition as a smooth configuration. In the real-time limit the instanton trajectory takes a wild ride through complex configuration space on its way between the minima of the potential, with the path resembling a space-filling curve. Yet its action remains equal to that of the Euclidean instanton. Our results are a step toward developing a non-perturbative understanding of physics using the real-time Feynman path integral formalism.

Acknowledgements: We are very grateful to Y. Tanizaki, G. Ba¸ sar, D. Dorigoni, G. Dunne, A. Joseph, and M. Shifman for discussions. A. C. acknowledges the partial support of the US DOE under grant number DE-SC0011842.

- ∗ acherman@umn.edu
- † unsal.mithat@gmail.com
- [1] R. P. Feynman, Rev. Mod. Phys. 20 , 367 (1948); R. Feynman and A. Hibbs, Quantum mechanics and path integrals , International series in pure and applied physics (McGraw-Hill, 1965).
- [2] For an earlier attempt at grappling with instantons in real time see e.g. [13].
- [3] A. Belavin, A. M. Polyakov, A. Schwartz, and Y. Tyupkin, Phys.Lett. B59 , 85 (1975).
- [4] F. Pham, Proc. Symp. Pure Math. 2 , 319 (1983).
- [5] E. Witten, preprint , 347 (2010), arXiv:1001.2933 [hepth]; preprint (2010), arXiv:1009.6032 [hep-th].
- [6] Y. Tanizaki and T. Koike, (2014), arXiv:1406.2386 [math-ph].
- [7] E. Bogomolny, Phys.Lett. B91 , 431 (1980); J. ZinnJustin, Nucl.Phys. B192 , 125 (1981); G. V. Dunne and M. Unsal, preprint (2014), arXiv:1401.5202 [hep-th].
- [8] P. Argyres and M. Unsal, Phys.Rev.Lett. 109 , 121601 (2012), arXiv:1204.1661 [hep-th]; P. C. Argyres and M. Unsal, JHEP 1208 , 063 (2012), arXiv:1206.1890 [hep-th]; G. V. Dunne and M. Unsal, JHEP 1211 , 170 (2012), arXiv:1210.2423 [hep-th]; Phys.Rev. D87 , 025015 (2013), arXiv:1210.3646 [hep-th]; A. Cherman, D. Dorigoni, G. V. Dunne, and M. Unsal, Phys.Rev.Lett. 112 , 021601 (2014), arXiv:1308.0127 [hep-th]; G. Basar,
- G. V. Dunne, and M. Unsal, JHEP 1310 , 041 (2013), arXiv:1308.1108 [hep-th]; A. Cherman, D. Dorigoni, and M. Unsal, (2014), arXiv:1403.1277 [hep-th]; T. Misumi, M. Nitta, and N. Sakai, JHEP 1406 , 164 (2014), arXiv:1404.7225 [hep-th]; T. Misumi and T. Kanazawa, JHEP 1406 , 181 (2014), arXiv:1405.3113 [hep-ph].
- [9] There are also related developments in string theory, see e.g. [14].
- [10] The instanton and anti-instanton get interchanged if -π/ 2 ≤ α &lt; 0.
- [11] D. Harlow, J. Maltz, and E. Witten, JHEP 1112 , 071 (2011), arXiv:1108.4417 [hep-th].
- [12] How to carry out this program completely for configuration-space path integrals is an open problem, but the present understanding is already sufficient for our Gaussian analysis.
- [13] D. Levkov and S. Sibiryakov, JETP Lett. 81 , 53 (2005), arXiv:hep-th/0412253 [hep-th].
- [14] M. Marino, R. Schiappa, and M. Weiss, Commun.Num.Theor.Phys. 2 , 349 (2008), arXiv:0711.1954 [hep-th]; M. Marino, JHEP 0812 , 114 (2008), arXiv:0805.3033 [hep-th]; M. Marino, R. Schiappa, and M. Weiss, J.Math.Phys. 50 , 052301 (2009), arXiv:0809.2619 [hep-th]; I. Aniceto, R. Schiappa, and M. Vonk, 6 , 339 (2012), arXiv:1106.5922 [hep-th]; R. Schiappa and R. Vaz, Commun.Math.Phys. 330 , 655 (2014), arXiv:1302.5138 [hep-th].
Aleksey Cherman Mithat Unsal
Real-Time Feynman Path Integral Realization of Instantons
In Euclidean path integrals, quantum mechanical tunneling amplitudes are
associated with instanton configurations. We explain how tunneling amplitudes
are encoded in real-time Feynman path integrals. The essential steps are
borrowed from Picard-Lefschetz theory and resurgence theory.