FTPI-MINN-14/18, UMN-TH-3343/14

Real-Time Feynman Path Integral Realization of Instantons

Aleksey Cherman1, ∗ and Mithat Unsal¨ 2, † 1Fine Theoretical Physics Institute, Department of Physics, University of Minnesota, Minnesota, MN 55455, USA 2Department of Physics, North Carolina State University, Raleigh, NC, 27695 In Euclidean path integrals, quantum mechanical tunneling amplitudes are associated with instan- ton 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] time ti to end up at position xf at time tf is is one of the several equivalent formulations of quantum Z x(tf )=xf mechanics (QM) and quantum field theory (QFT). The xf ,ti i S[x] ~ Axi,ti = d[x(t)] e (1) amplitude of an event is found by summing over all path x(ti)=xi histories with the contribution of each path γ weighed by its action S[γ] via the factor eiS[γ]/~. For real paths, where S is the classical action. For a non-relativistic this is a sum over pure phases. The real-time path in- bosonic particle moving in a symmetric double-well po- tegral formulation is very useful in many settings, and tential the action can be written as yields beautiful explanations of many quantum phenom- Z tf 1  2 2 2 ena, such as the interference patterns seen in double-slit S[x] = 2 dt (∂tx) − (x − 1) (2) 2g t experiments and the Aharonov-Bohm effect. Yet surpris- i ingly, some very basic quantum phenomena do not have a Our semiclassical treatment of tunneling will be justified known description in the real-time path integral. For in- if g  1. If the particle is initially near e.g. x = −1 at stance, tunneling phenomena, which are well-understood t = ti the amplitude for it to end up near x = +1 will be in the operator (Hamiltonian) formalism, as well as in zero to any order in a perturbative expansion in g. the Euclidean-time path integral formulation, are not Let us now focus on the tunneling amplitude from yet understood in the real-time Feynman path integral xi = −1 at ti = −∞ to xf = +1 at tf = +∞. The stan- formulation[2]. dard way to compute this amplitude using path integrals We focus on a particle in a symmetric double-well po- involves moving to Euclidean time by a Wick rotation tential as a paradigmatic example. The level splitting be- t → −iτ, x(t) → x(τ). The Euclidean action tween the ground state and first excited state is given by 1 Z ∞ a non-perturbative factor, ∼ exp[−A/(g2 )] where A is a S = dτ (∂ x)2 + (x2 − 1)2 . (3) ~ E 2g2 τ pure number and g is a parameter controlling the depth −∞ of the wells. It is a consequence of tunneling in the op- is then positive semi-definite and enters the path integral erator formalism or instanton saddles[3] in the Euclidean −SE 2 as e . The Euclidean equation of motion ∂τ x(τ) = path integral formulation, but its origin is mysterious in 2x(τ) x(τ)2 − 1
has two finite-action solutions, the in- the Feynman path integral. + − stanton xI and anti-instanton xI : Understanding such non-perturbative phenomena di- ± rectly in real-time path integrals is likely to be very im- xI (τ) = ± tanh(τ − τ0), (4) portant. Euclidean path integrals give deep insights into where τ is a zero mode parameter. These solutions in- the non-perturbative properties of QM and QFT, but 0 arXiv:1408.0012v2 [hep-th] 19 Aug 2014 terpolate from x = ∓1 to x = ±1 as τ goes from −∞ to there is a huge range of problems for which they are not +∞. In what follows we will focus on x+ and will drop useful, such as calculations of scattering amplitudes or I the superscript. The Euclidean action of (4) is S = 4 , non-equilibrium observables. These questions are more I 3g2 so that (4) is a saddle-point of the Euclidean path inte- naturally handled using real-time path integrals. Yet − 4 the current understanding of real-time path integrals be- gral, with contributions weighed by e 3g2 . For g  1 yond the perturbative level is extremely rudimentary, as the desired tunneling amplitude can be computed using is highlighted by the fact that even tunneling phenomena a steepest-descent approximation, and is given by are not understood in the framework. Our results are a r 4 32 − 2 step toward developing a non-perturbative understanding 3g ~ A = 2 e (5) of real-time path integrals. πg Euclidean-time instantons—To set notation we review where the prefactor is obtained by integrating over Gaus- the standard path integral treatment of tunneling in QM. sian and zero mode fluctuations around x = xI (τ). The real-time Feynman path integral representation of Real-time instantons—The Euclidean-time derivation the amplitude A for a particle initially at position xi at is not fully satisfying, as physics happens in real time. 2

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 τ.

Yet there is a good reason that instantons were never the bottoms of the wells. formulated in the real-time Feynman path integral. A (iii) It has singularities at t = (2n + 1)π/2, n ∈ Z, and natural attempt to do this involves looking for finite- these singularities are not integrable, so that the action action solutions of the real-time equation of motion is not well-defined. Until fairly recently this is where an analysis would −∂2x(t) = 2x(t) x(t)2 − 1
. (6) t have had to stop. Fortunately, there has recently been Since we are looking for an answer of the form eiS/~ = major progress in the understanding the complexification A − 2 of functional integrals [4–6] and the structure of the semi- e g ~ the action of the solution would have to be imag- inary. That means that a ‘real-time instanton’ configu- classical expansion of path integrals in QM [7] and QFT ration has to be complex. Such a notion raises major [8] via resurgence theory[9]. Motivated by these devel- conceptual issues: opments, let us see what happens if we do not insist on (A) If the field space is complexified, how do we avoid working in either purely real or imaginary time. an undesired doubling of degrees of freedom in the com- Specifically, let us start in real time t and do a complex −iα plexified path integral? Wick rotation t → τe , φ(t) → φ(τ). In contrast to (B) Relatedly, how do we avoid including pathological purely real or purely imaginary time, now the action itself directions in the complex field space which may yield ex- becomes complex: ponentially large contributions? Over what sub-space of e−iα Z +∞ S = dτ e+2iα(∂ x)2 − (x2 − 1)2 . (8) fields (complex paths) should one integrate over? α 2 τ 2g −∞ Setting these important questions aside for the mo- ment, we note that it is tempting to guess that the de- Hence issues (A) and (B) with the complexification of sired saddle field configurations might be given by the configuration space in path integrals must be deal with, Wick rotation of the Euclidean instanton back to real and we do so below. For small α, τ ≈ t and Sα ≈ S, time. Indeed, (6) does have the formal solution while for α = π/2 τ becomes the standard Euclidean time coordinate and iSπ/2 = −SE. It turns out that x(t) = i tan(t − t0) (7) an infinitesimal α acts as a crucial regulator which is which is a Wick rotation of (4). But (7) looks non- necessary to make sense of real-time instantons. 2iα 2 sensical, for many reasons: The α-dependent equation of motion −e ∂τ x(τ) = 2
(i) It is pure-imaginary-valued, so it cannot possibly give 2x(τ) x(τ) − 1 has the solutions A iS/ − 2 h i a contribution like e ~ = e g ~ , and cannot describe a −i(α−π/2) xI (τ) = ± tanh (τ − τ0)e (9) path interpolating from x = −1 to +1. (ii) It does not appear to have anything to do with tun- To make (9) a legitimate solution one must complex- neling from x = −1 to +1, since it is never localized near ify the space of path histories from x : R → R to 3

−2 π FIG. 2. (Color Online.) Upper-left: trajectory of (9) for α = 4 × 10 × 2 in complex configuration space. Upper-right: −2 π real and imaginary parts of (9) for α = 1 × 10 × 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 −2 π α = 1 × 10 × 2 ), while the real part approaches to a distribution, a signed Dirac comb. Lower-right: the trajectory for −3 π α = 1 × 10 × 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. x : Re−iα → C. While (9) looks like a simple analytic illustrated in Fig.2. The Re x < 0 and Re x > 0 half- continuation of (4), the conceptual issues raised by al- planes form two basins of attraction for the repulsive and lowing such an continuation are quite complex. attractive fixed points x = −1 and x = +1 respectively. Naively, precisely at α = 0 (9) reduces to (7), although After starting at x = −1 at τ = −∞, as α → 0 the we are about to see that the limit is subtle. Consider instanton spends more and more time whipping around the situation for 0 < α ≤ π/2 [10]. Remarkably, the in the Re x < 0 region in a huge number of more and imaginary parts of the complex instanton solution (9) more tightly spaced increasingly large arcs for τ < τ0, then jumps over to the Re x > 0 region and winds closer vanish for large |τ|, and xI approaches ∓1 as τ → ∓∞! The behavior of the complexified instanton at generic α and closer to x = +1 in the same manner, until settling values is summarized in Fig.1 and Fig.2. For generic at x = +1 at τ = +∞. From the bottom-right plot 0 < α ≤ π/2 field starts at x = −1 at τ = −∞, whips in Fig.2, the complex instanton trajectory appears to become a space-filling curve in the real-time limit. around in the left half of complex-x plane for τ < τ0, then crosses to the right half of the plane at τ = τ0, whips around x = +1 for τ > τ0, and then settles down at x = +1 at τ = +∞. The instanton remains smooth for any α 6= 0, but the In the strict α = 0 limit, the imaginary part of the behavior in the real time limit α → 0 is singular, as instanton solution is given by (7), while the real part 4 approaches to a distribution, given by of the integral thus involves finding a basis of convergent ∞ middle-dimensional integration cycles in the complexified X π Re x = πsign(t − t ) δ(t − t − (2k − 1)) (10) integral, and then working out how the original real cy- 0 0 2 k=−∞ cle decomposes into the complex cycles. Ref. [5] showed that in favorable circumstances, this can be done by us- which is a signed Dirac comb. The time average of ing Picard-Lefshetz theory[12]. Every finite-action saddle 1 R T hRe xiT and hIm xiT (which we define as h·iT = · dτ T τ0 point of the path integral is associated to a convergent with T = πk + τ0, k ∈ Z) for any α, including α = 0, es- integration cycle called a ‘Lefshetz thimble’ — a multidi- sentially follows the pattern of a Euclidean instanton, mensional generalization of a steepest-descent path. The  thimbles are solutions of the downward flow equations for +1, τ > τ0 hRe xiT = the Morse function Re [iS]. For us, the flow equation is −1, τ < τ0 hIm xi = 0 (11) +i(α−π/2) T ∂x e  −2iα 2 2
 = − e ∂τ x¯ + 2¯x x¯ − 1 (13) as illustrated by Fig.2. ∂s g2 Next, we show that the action −iS of the instanton α where x is complex and s is the flow ‘time’ along the thim- is always precisely 4/(3g2) for any α, with the action ble. The fluctuations around x are described by the so- at α = 0 defined by continuation from |α| > 0. The I lutions of (13) with the initial condition x(s → −∞, τ) = Lagrangian evaluated on (9) oscillates wildly through the x (τ). The structure of the Gaussian fluctuations can be complex plane as α → 0. Nevertheless, −iS [x ] turns I α I determined by writing x(s, τ) = x (τ) + δx(s, τ) with out to be independent of α. To show this we rewrite I δx(s → −∞, τ) = 0 and expanding (13) to O(δx): −iα Z +∞ e  +iα 0 2 +iα 0 Sα(x) = dτ e x˙ − W (x) + 2e xW˙ +i(α− π ) 2 ∂δx e 2 2g −∞ = − e−2iα∂2 + 6x2(τ) − 2 δx (14) ∂s g2 τ I 0 0 2 where ≡ ∂x,˙ ≡ ∂τ , and (W ) = V (x). But (9) satisfies the first-order equation e+iαx˙ − W 0(x) = 0, so The existence of the flow equations means that it is pos- sible to choose sensible integration cycles in the com- 1 Z +∞ 4 0 plexified path integral, and in the Gaussian-fluctuation −iS(xI ) = 2 dτx˙ I W (xI ) = 2 (12) g −∞ 3g limit the determination of our desired cycle reduces to We conjecture that the solution with minimal non-zero finding the solution of (14). By moving to the variable 0 action is (9) for any α. This is obvious at α = π/2, while τ = τ exp[−i(α − π/2)] and rescaling δx one can map at α = 0, resurgence theory implies that the minimum- (14) to the Euclidean flow equation. This implies that integrating over the thimble amounts to calculating the action saddle must have action −iS(xI ) due to the large- order structure of perturbation theory. We leave a de- standard fluctuation determinant (see [6] for a nice dis- 2 tailed exploration of this conjecture to future work. cussion of this) multiplying e−4/3g , and reproduces (5). The question of whether singular but finite-action Conclusions—We have shown how to compute tun- field configurations contribute to path integrals is neling amplitudes between low-lying states, a hallmark controversial[11]. Our results support the notion that non-perturbative feature of quantum mechanics, directly such configurations do contribute, since (9) is a finite from a real-time Feynman path integral. The amplitude action configuration whichis associated with a dramatic is given by a real-time instanton saddle point, which re- physical effect of level splitting. Hence it is physically quires a complexification of the space of paths and the clear that it must be included it in the evaluation of the addition an infinitesimal imaginary part to the time coor- path integral, even though it becomes singular and dis- dinate as a regulator for its definition as a smooth config- continuous in the strict α = 0 limit. uration. In the real-time limit the instanton trajectory Lefshetz thimbles—To find the contribution of the com- takes a wild ride through complex configuration space plex instantons to observables in the real time limit, it on its way between the minima of the potential, with is necessary to take into account fluctuations around the the path resembling a space-filling curve. Yet its action instanton saddle points. As we already noted this is a remains equal to that of the Euclidean instanton. Our subtle problem due to the complexification of the field results are a step toward developing a non-perturbative space, and involves understanding the analytic continua- understanding of physics using the real-time Feynman tion of path integrals. A program to study this issue was path integral formalism. recently initiated in [5]. Acknowledgements: We are very grateful to A sensible complexification of the integration cycle of Y. Tanizaki, G. Ba¸sar, D. Dorigoni, G. Dunne, an integral which initially has a real domain of inte- A. Joseph, and M. Shifman for discussions. A. C. gration must not double the number of degrees of free- acknowledges the partial support of the US DOE under dom. A full understanding of the analytic continuation grant number de-sc0011842. 5

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. Mis- ∗ [email protected] umi, M. Nitta, and N. Sakai, JHEP 1406, 164 (2014), † [email protected] arXiv:1404.7225 [hep-th]; T. Misumi and T. Kanazawa, [1] R. P. Feynman, Rev. Mod. Phys. 20, 367 (1948); JHEP 1406, 181 (2014), arXiv:1405.3113 [hep-ph]. R. Feynman and A. Hibbs, Quantum mechanics and path [9] There are also related developments in string theory, see integrals, International series in pure and applied physics e.g. [14]. (McGraw-Hill, 1965). [10] The instanton and anti-instanton get interchanged if [2] For an earlier attempt at grappling with instantons in −π/2 ≤ α < 0. real time see e.g. [13]. [11] D. Harlow, J. Maltz, and E. Witten, JHEP 1112, 071 [3] A. Belavin, A. M. Polyakov, A. Schwartz, and Y. Tyup- (2011), arXiv:1108.4417 [hep-th]. kin, Phys.Lett. B59, 85 (1975). [12] How to carry out this program completely for [4] F. Pham, Proc. Symp. Pure Math. 2, 319 (1983). configuration-space path integrals is an open problem, [5] E. Witten, preprint , 347 (2010), arXiv:1001.2933 [hep- but the present understanding is already sufficient for th]; preprint (2010), arXiv:1009.6032 [hep-th]. our Gaussian analysis. [6] Y. Tanizaki and T. Koike, (2014), arXiv:1406.2386 [13] D. Levkov and S. Sibiryakov, JETP Lett. 81, 53 (2005), [math-ph]. arXiv:hep-th/0412253 [hep-th]. [7] E. Bogomolny, Phys.Lett. B91, 431 (1980); J. Zinn- [14] M. Marino, R. Schiappa, and M. Weiss, Com- Justin, Nucl.Phys. B192, 125 (1981); G. V. Dunne and mun.Num.Theor.Phys. 2, 349 (2008), arXiv:0711.1954 M. Unsal, preprint (2014), arXiv:1401.5202 [hep-th]. [hep-th]; M. Marino, JHEP 0812, 114 (2008), [8] P. Argyres and M. Unsal, Phys.Rev.Lett. 109, 121601 arXiv:0805.3033 [hep-th]; M. Marino, R. Schi- (2012), arXiv:1204.1661 [hep-th]; P. C. Argyres and appa, and M. Weiss, J.Math.Phys. 50, 052301 (2009), M. Unsal, JHEP 1208, 063 (2012), arXiv:1206.1890 arXiv:0809.2619 [hep-th]; I. Aniceto, R. Schiappa, [hep-th]; G. V. Dunne and M. Unsal, JHEP 1211, and M. Vonk, 6, 339 (2012), arXiv:1106.5922 [hep-th]; 170 (2012), arXiv:1210.2423 [hep-th]; Phys.Rev. D87, R. Schiappa and R. Vaz, Commun.Math.Phys. 330, 655 025015 (2013), arXiv:1210.3646 [hep-th]; A. Cherman, (2014), arXiv:1302.5138 [hep-th]. D. Dorigoni, G. V. Dunne, and M. Unsal, Phys.Rev.Lett. 112, 021601 (2014), arXiv:1308.0127 [hep-th]; G. Basar,