Download The Power of Perturbation Theory

Survey
yes no Was this document useful for you?
   Thank you for your participation!

* Your assessment is very important for improving the workof artificial intelligence, which forms the content of this project

Document related concepts

Interpretations of quantum mechanics wikipedia , lookup

EPR paradox wikipedia , lookup

Coupled cluster wikipedia , lookup

Quantum state wikipedia , lookup

Propagator wikipedia , lookup

Feynman diagram wikipedia , lookup

Orchestrated objective reduction wikipedia , lookup

Quantum field theory wikipedia , lookup

Asymptotic safety in quantum gravity wikipedia , lookup

Canonical quantization wikipedia , lookup

Quantum chromodynamics wikipedia , lookup

Quantum electrodynamics wikipedia , lookup

History of quantum field theory wikipedia , lookup

Hidden variable theory wikipedia , lookup

Topological quantum field theory wikipedia , lookup

Renormalization wikipedia , lookup

Renormalization group wikipedia , lookup

Scalar field theory wikipedia , lookup

Perturbation theory (quantum mechanics) wikipedia , lookup

Path integral formulation wikipedia , lookup

Yang–Mills theory wikipedia , lookup

Perturbation theory wikipedia , lookup

Instanton wikipedia , lookup

Transcript
The Power of Perturbation Theory
Marco Serone, SISSA, Trieste
based on 1612.04376 and to appear, with G. Spada and G. Villadoro
February 2017
1
Main result:
All observables in a one-dimensional quantum
mechanical system with a bounded potential can be
entirely computed from a single perturbative series
2
Introduction and Motivation
It is known that perturbative expansions in QM and QFT
are only asymptotic with zero radius of convergence
[Dyson, 1952]
In special cases, like the anharmonic oscillator in QM and 42,3 QFT
it has been proved that perturbation theory is Borel resummable
and leads to the exact results (no non-perturbative contributions)
[Loeffel et al, 1969; Simon and Dicke, 1970; Eckmann, Magnen and
Seneor, 1975; Magnen and Seneor, 1977]
Most perturbative expansions are not Borel resummable
Technically, this is due to the appearance of singularities of the Borel
function along an integral needed to perform to get the answer
The singularity can be avoided by deforming the contour
but the result present an ambiguitiy of order exp( a/ n )
3
If one is able to find a semi-classical instanton like configuration (and its whole
series) leading to the same factors, one might hope to remove the ambiguity
A systematic way to proceed along these lines uses the theory of resurgence
[Ecalle, 1981]
Unfortunately, one encounters various difficulties using this approach
• Even if Borel summability is somehow proven, convergence to
the exact result requires the knowledge of some analytic
properties of the observable, in general hard to verify
•
Beyond the weak coupling regime, one needs to consider an
infinite number of instantons (impractical)
4
Alternative geometric picture starts from path integral
Perturbation theory is infinite dimensional generalization of
steepest-descent method to evaluate ordinary integrals
Understanding which instantons contribute to a given observable amounts to
understand which saddles of the action should be considered in the path integral
We can hope to address this point using generalization to
infinite dimensions of known mathematical methods.
[Witten, 2011]
They provide a geometric alternative to the resurgence program
Key question:
Under what conditions only one saddle point (trivial vacuum)
contributes so that perturbation theory gives the full answer?
We answered this question in QM.
5
Plan
One dimensional integrals
Lefschetz thimbles
Borel sums
Exact perturbation theory
Path integral
Examples
Conclusions
6
One dimensional integrals
1
Z( ) ⌘ p
Convergent for any
Z
1
dx g(x) e
f (x)/
1
0
We can compute the integral exactly using steepest-descent methods:
1.
2.
Determine the saddle points contributing to Z
Resum the “perturbative” expansion around each saddle
Consider first point 1.
Analytically continue in complex plane
1
Z( ) = p
Z
dz g(z) e
Cx
7
f (z)/
Saddle points z : f 0 (z ) = 0
For each saddle we call downward and upward flows J and K the
lines where Re F decrease and increase, respectively, and Im F is constant
F (z) ⌘
f (z)/
J is the path of steepest descent. If it flows to Re F =
1
it is called a Lefschetz thimble (or just thimble).
By construction, the integral over a thimble is always well defined and convergent.
If J hits another saddle point, the flow splits int two branches
and ambiguity arises (Stokes line).
Picard-Lefschetz theory: deform the contour Cx in a contour C :
C=
h. . .i
X
n = hCx , K i
J n
denotes intersection pairings
hJ , K⌧ i =
8
⌧
[Witten, 1001.2933]
If J = K⌧ intersection not well defined.
Ambiguity resolved by assigning a small imaginary part to
Z( ) =
1
Z ( )⌘ p
X
Z
n Z ( )
dz g(z) e
f (z)/
J
Not all saddles contribute to the integral, only those with n 6= 0
Example: i)
1 2 1 4
f (x) = x + x ,
2
4
g(x) = 1
Three saddles: z0 = 0, z± = ±i
Degenerate situation. Give a small imaginary part to
9
1
Im
1
>0
Im
4
4
J+
K+
K+
2
J+
2
z+
J0
z0
0
z≠
K0
≠2
≠2
J≠
0
J0
z≠
≠2
K≠
K0
z+
z0
0
≠4
≠4
<0
2
≠4
≠4
4
n0 = 1, n± = 0
J≠
≠2
0
K≠
2
Only saddle at the origin contributes to the integral.
No need to deform the original contour of integration (thimble itself)
No “non-perturbative” contributions
10
4
1 2 1 4
Example: ii)
f (x) =
x + x ,
2
4
Three saddles: z0 = 0, z± =1 ±1
Im
4
2
>0
Im
0
J≠
z0
K0
≠2
C+ = J
J+
z+
≠2
K+
J0
0
0
2
≠4
≠4
4
J0 + J+
<0
K0
2
z≠
≠4
≠4
Again degenerate situation. 1
4
K≠
≠2
g(x) = 1
J≠
z≠
J0
z0
K+
z+
J+
K≠
≠2
0
2
C = J + J0 + J+
The integral is on a Stokes line and the intersection numbers depend on the
deformation. Ending result of the integral will eventually be unambiguous
There are “non-perturbative” contributions
11
4
Borel Sums
1
Z ( )⌘ p
Z
f (z)/
dz g(z) e
J
can be computed using a saddle-point expansion. Resulting series is asymptotic
Z( )
N
X
Zn
n
n=0
= O(
N +1
),
!0
as
Borel transform to reconstruct function: assume Zn ⇠ n!an nc
Bb Z(t) =
1
X
n=0
Zn
tn
(n + b + 1)
ZB ( ) =
(Borel Le Roy function)
Z
1
0
dt e t tb Bb Z(t )
If no singularities for t>0 and under certain assumptions
ZB ( ) = Z( )
12
Analytic structure of Borel function close to the origin can be
determined by the large-order behaviour of the series coefficients
n
for Zn ⇠ n!a , B0 Z(t) =
P
n
(at)
⇠
n
1
1 at
Singularity dangerous or not depending on the sign of a:
a>0 (same sign series) singularity for t>0
☹
a<0 (alternating series) singularity for t < 0
😀
Lateral Borel summation: move the contour off the real positive
axis to avoid singularity for t>0.
Result is ambiguous and ambiguity signals presence of
non-perturbative contributions to Z not captured by ZB
Deformation to define the lateral Borel sum corresponds to
the one we did to avoid Stokes line
13
Thimbles are Borel Resummable
Z
1
Z ( )⌘ p
dz g(z) e
f (z)/
J
is Borel resummable to the exact result.
Proof: Change variable
t=
f (z)
f (z )
[Berry, Howls 1991]
For any z 2 J , t is real and non-negative
Z ( )=e
f (z )
Z
1
dt t
1/2
t
e B ( t)
0
X g(zk ( t)) p
B ( t) ⌘
t .
0
|f (zk ( t)|
k=1,2
By definition, on a thimble f 0 (z) 6= 0 for z 6= z
=) B (t) regular for t > 0
But this is exactly the expression one would get by
considering the Borel transform with b=-1/2!
q.e.d.
14
Summarizing
1
Z( ) ⌘ p
Z
1
dx g(x) e
f (x)/
=
1
X
n Z ( )
Each Z ( ) is Borel resummable to the exact result
1 saddle contributes: Z( ) reconstructable from perturbation theory
More saddles contribute: Z( ) given by multi-series
(non-perturbative corrections)
1 2 1 4
f (x) = x + x
2
4
Example first case:
Zn =
p
2( )
(2n +
n!
n
Z
1
0
dt t
1
2)
1
2
B
1/2 Z(t)
=
s
p
1 + 1 + 4t
1 + 4t
⇣ 1 ⌘
1
1
e t B 1/2 Z( t) = p e 8 K 14
8
2
15
Exact result
Example second case:
Coefficients at minima:
f (x) =
Z±,n =
1/2 Z± (
For real
2n +
n!
1
2
iZ0,n
Coefficients at maximum:
B
1 2 1 4
x + x
2
4
t) =
B
s
p
1+ 1 4 t
2(1 4 t)
1/2 Z
Recall we need Im
not Borel resummable
6= 0 to avoid Stokes line
The thimble decomposition instructs us how to consistently perform lateral
Borel sums and consistently sum up the different contributions:
Im
Im
>0
<0
Z ( )
Z0 ( ) + Z+ ( )
1
⇡
= p e8
4
Z ( ) + Z0 ( ) + Z+ ( )

I
1
4
✓
1
8
◆
+ I 14
✓
1
8
◆
Exact result
16
Exact Perturbation Theory
When the decomposition in thimbles is non-trivial, the
intersection numbers n have to be determined
This is easy for one-dimensional integrals, but very complicated in
the path integral case, where in general they will be infinite
Key idea: the thimble decomposition can be modified by a simple trick
The decomposition of the integral Z in terms of thimbles
is governed by the function f(z) and not by g(z). Define
1
p
Ẑ( , 0 ) ⌘
fˆ(x) ⌘ f (x) + f (x) ,
Z
ĝ(x,
1
dx e
fˆ(x)/
0)
0
0)
1
⌘ g(x)e
f (x)/
By construction
Ẑ( ,
ĝ(x,
= ) = Z( )
17
0
lim
|x|!1
f (x)
=0
f (x)
Saddle-point expansion in
at fixed
0:
f in fˆ: “classical” deformation
f in ĝ: “quantum” deformation”
Thimble decomposition is determined by fˆ
By a proper choice of f , we can generally trivialize the thimble decomposition
=) one real saddle and no need to deform the contour of integration
Series expansion of Ẑ( ,
0)
in
at fixed
0:
Exact Perturbation Theory (EPT)
Non-perturbative contributions of Z are all contained
in the perturbative expansion of Ẑ
18
Consider again
and choose
1 2 1 4
x + x
2
4
f (x) =
f (x) = x2
1 2 1 4
ˆ
f (x) = x + x ,
2
4
so that
ĝ(x,
0)
= exp
⇣
x2 ⌘
0
Only saddle at origin contributes. We have
s
p
p
1 + 1 + 4 t 1+4 t
0
B 1/2 Ẑ( t, 0 ) =
e
1+4 t
Z
1
0
dt t
1
2
e tB
1/2 Ẑ( t, )
1
⇡
= p e8
4

I
1
4
✓
1
8
Exact result is now completely perturbative
19
1
◆
+ I 14
✓
1
8
◆
Path Integral
These results can be generalized to higher-dimensional integrals
n/2
Z ( )=
Z ( )=e
f (z )/
B ( t) ⌘ ( t)1
= ( t)
n/2
1 n/2
Z
1
Z0
Z
J
⌦
Z
dt tn/2
dz g(z) e
f (z)/
J
1
e
t
B ( t) ,
dz g(z) [f (z)
g(z)
d⌃
|rf (z)|
t
20
f (z )
t]
We can take the continuum for the path integral case.
Consider quantum mechanics where divergencies are all reabsorbed by the
path integral measure (UV limit not expected to be problematic)
Z( ) =
Z
S[x] =
Dx(⌧ ) G[x(⌧ )]e
Z
d⌧

S[x(⌧ )]/
1 2
ẋ + V (x)
2
V (x) analytic function of x
Discrete spectrum
V (x) ! 1 for |x| ! 1
Key result:
If the action S[x(⌧ )] has only one real saddle point x0 (⌧ ), the formal series
expansion of Z( ) around = 0 is Borel resummable to the exact result.
21
Remarks:
1. Number of saddles depends on boundary conditions
and the infinite time limit
2. Depending on the system, some observables are, and
some other are not, Borel reconstructable
If V admits only one minimum, the above subtleties do not arise.
Like in the one-dimensional case, we can define an EPT as follows
Suppose V = V0 +
V such that
1. V0 has a single non-degenerate critical point (minimum);
2. lim|x|!1
Ẑ( ,
0)
=
Z
V /V0 = 0 .
Dx G[x] e
R
d⌧
V
0
e
S0
,
Ẑ( , ) = Z( )
22
S0 ⌘
Z

1 2
d⌧
ẋ + V0
2
Ẑ is guaranteed to be Borel resummable to the exact answer
We have then proved the statement made at the beginning:
All observables in a quantum mechanical system
with a bounded potential can be entirely computed
from a single perturbative series (EPT)
(provided points 1. and 2. above apply)
No need of resurgence, thimble decomposition, trans-series.
Even observables in systems known to have instanton corrections
will be reconstructed by a single perturbative series
Large arbitrariness in the choice of EPT. In principle all choices equally
good, although in numerical studies some choices better than others
We will denote by Standard Perturbation Theory (SPT) the usual
perturbative computation (instantons included)
Time to show explicit results
23
Examples
Numerical analysis to test our results. We compute N orders using package
[Sulejmanpasic, Unsal, 1608.08256]
At fixed order N of perturbative terms, we use Padé approximants
for the Borel function, that is then numerically integrated.
Results are compared with other methods (such as
Rayleigh-Ritz) to get exact answers in QM
p 3
1 2
V (x; ) = x + ↵ x + x4
2
2
p
For |↵| < 2 2/3 V has a unique minimum at x = 0
In principle no need of EPT, being SPT
Borel resummable
Yet, EPT greatly improves the accuracy of result at a
given order in perturbation theory
24
1 2
V0 = x + x 4 ,
2
2
3/2 3
V =↵
x
=5
= 1/20
1
10-21 10-10
SPT
10-31
10
300
400
N
500
10-20
-41
200
SPT
EPT
E0/E0
E0/E0
10-11
EPT
10-30
600
200
300
400
N
500
At weal coupling SPT better than EPT, viceversa at strong
coupling, where SPT is inaccurate
25
600
Symmetric double well
V =
2
✓
x2
1
4
◆2
Prototypical system where instanton occurs. Write
⇣ 1
⌘
⇣
⌘ x2
1
0 2
V0 =
+ x + x4 ,
V =
0+
32
2
2
2 2
Already at small coupling, EPT is able to resolve the
instanton contribution to the ground state energy.
At
E0
⇡ 10
E0
=
0
= 1/25, N = 200
E1
⇡ 4 · 10
E1
8
E0inst
2
p
⇡
e
⇡
26
1
6
⇡ 0.087
14
At larger values of the coupling, instanton computation breaks down
EPT works better and better
No need of many perturbative terms.
With just 2 (using a conformal mapping method) at
= 1 we get
E0
' 3%
E0
We have also computed higher energy states and eigenfunctions.
In all cases EPT give excellent agreement with other methods. No
signal of missing non-perturbative contributions.
27
SUSY double well
Notable tilted double well:
V (x; ) =
2
⇣
x2
1 ⌘2 p
+
x
4
Potential one obtains integrating out fermions in SUSY QM
E0 = 0
to all orders in SPT due to SUSY
Yet, E0 6= 0 by instanton effects that dynamically break supersymmetry
28
We can approach this system in three ways:
1. ordinary instanton methods, SPT
[Jentschura, Zinn-Justin, hep- ph/0405279]
2. Turn the quantum tilt into classical. Alternative
perturbation theory (APT), sort of inverse of EPT
VAPT =
2
⇣
x2
1 ⌘2
0
+p x
4
3. EPT
⇣ 1
⌘
0 2
4
V0 =
+ x + x ,
32
2
2
0
V =p x
29
⇣
1 ⌘ x2
0+
2 2
1
E0
10 -5
+
10 -10
+
×
+
×××
×
+
+×
+
+×
+×
+×
+×
+×
+×
+
+
+
×
×
+
×
+
×
+
×
+
×
+
×
+
×
×
+
×
+
×
+
×
+
×
+×
+×
+×
×
+
×
+
×
+
×
+
10 -15 +
ΔE 0 /E 0
1
10
-10
10
-20
10
-30
+ + + + +××
++++++
10 -40 -2
10
SPT
××
××
++++++++++++++
××
××
0.1
××
××
1
λ
N=200
30
++++++++
× ×EPT
××××××
10
APT
×
×
××
×
×
×
100
False vacuum
1.0
fv
E0
0.5
fv
V (x;1,∞)
V (x)
0.0
V fv (x;1)
- 0.5
fv
E0
fv
- 1.0
E0 (EPT)
- 1.5
- 2.0
-1
0
1
x
31
2
Pure anharmonic potentials:
H = p2 + x2l
Usual perturbation theory breaks down
k
(4)
Ek
(4)
(4)
Ek /Ek
0
1.0603620904
3 · 10
1
3.7996730298
2 · 10
2
7.4556979379
9 · 10
3
11.6447455113
4 · 10
4
16.2618260188
4 · 10
(6)
Ek
(6)
45
0.5724012268
44
2 · 10
2.1692993557
37
3 · 10
4.5365422804
36
9 · 10
7.4675848174
36
7 · 10
10.8570827110
2 · 10
Table 2: Energy eigenvalues Ek(2`) and the corresponding accuracies
(6)
Ek /Ek
(2`)
19
19
17
16
16
(2`)
Ek /Ek of the first five levels
of the pure anharmonic x4 and x6 potentials, computed using EPT with N = 200. Only the first ten
digits after the comma are shown (no rounding on the last digit).
of the pure x4 and x6 oscillators computed comparing EPT to the results from RR methods.
We used N = 200 orders of perturbation theory and in eq.(4.19) we chose c1 = 2 for ` = 2 and
(4)
c1 = 4, c2 = 2 for ` = 3. Notice the accuracy of E0 up to 45 digits! At fixed N , similarly to the
32
Conclusions and future perspectives
We have introduced a new perturbation theory in quantum mechanics
(EPT) that allows us to compute observables perturbatively
No need of trans-series and to worry for possibly ambiguous results
(possible addressed by resurgence or by thimble decomposition)
Is it possible to get the same results directly in Minkowski space?
Generalize to higher dimensional QM systems
33
What about QFT?
Is renormalization an issue?
UV problem
Infinite volume limit, spontaneous symmetry breaking? IR problem
4
Borel summability of 2,3 theories suggest a promising
development in QFT and is the next thing to do
A new approach to strongly coupled physics has begun …
34
Thank You
35