On the Lefschetz thimbles structure of the Thirring model F. Di - - PowerPoint PPT Presentation

on the lefschetz thimbles structure of the thirring model
SMART_READER_LITE
LIVE PREVIEW

On the Lefschetz thimbles structure of the Thirring model F. Di - - PowerPoint PPT Presentation

On the Lefschetz thimbles structure of the Thirring model F. Di Renzo 1 K. Zambello 1 , 2 francesco.direnzo@pr.infn.it kevin.zambello@pr.infn.it Wuhan, 18/06/2019 1 University of Parma and INFN - Gruppo collegato di Parma 2 Speaker Outline


slide-1
SLIDE 1

On the Lefschetz thimbles structure of the Thirring model

  • F. Di Renzo1
  • K. Zambello1,2

francesco.direnzo@pr.infn.it kevin.zambello@pr.infn.it

Wuhan, 18/06/2019

1 University of Parma and INFN - Gruppo collegato di Parma 2 Speaker

slide-2
SLIDE 2

Outline

Introduction The sign problem Lefschetz thimbles regularization On the Lefschetz thimbles structure of the Thirring model The theory Lefschetz thimbles regularization Conclusions

slide-3
SLIDE 3

Introduction

The sign problem

◮ We would like to compute expectation values

O ≡ 1 Z

  • Dψ O[ψ]e−S[ψ]

by importance sampling, sampling fields configurations from P ∝ e−S[ψ].

◮ However, there are cases of physical interests where S is

complex → sign problem

◮ Some interesting approaches to attack the sign problem are

those in which one complexifies the field variables (these are the subject of today’s and tomorrow’s parallel talks).

slide-4
SLIDE 4

Introduction

Lefschetz thimbles regularization

◮ One of such approaches is thimble regularization. Idea:

complexify the degrees of freedom of the theory and deform the integration paths. Picard-Lefschetz theory: attached to each critical point pσ exists a manifold Jσ s.t.

  • C

dzn O(z) e−S(z) =

  • σ

nσe−iSI

σ

dzn O(z) e−SR

σ

◮ The thimble Jσ attached to a critical point pσ is the union of

the steepest ascent paths leaving the critical points dzi dt = ∂ ¯ S ∂ ¯ zi , with i.c. zi(−∞) = zσ,i Along the flow, the imaginary part of the action is constant.

◮ The tangent space at pσ is spanned by the Takagi vectors,

which can be found by diagonalizing the Hessian at the critical point H(pσ)v(i) = λσ

i ¯

v(i)

slide-5
SLIDE 5

Introduction

Lefschetz thimbles regularization

◮ A natural parametrization for a point on the thimble is

z ∈ Jσ ↔ (ˆ n, t), where ˆ n defines the direction along which the path leaves the critical point and t is the integration time.

◮ Using this parametrization, the thimbles decomposition of an

expectation value O takes the form O =

  • σ nσ
  • D ˆ

n 2

i λσ i n2 i

  • dt e−Seff (ˆ

n,t) O eiω(ˆ n,t)

  • σ nσ
  • D ˆ

n 2

i λσ i n2 i

  • dt e−Seff (ˆ

n,t) eiω(ˆ n,t)

where V (ˆ n, t) is the parallel transported basis, Seff (ˆ n, t) = SR(ˆ n, t) − log |detV (ˆ n, t)| and ω(ˆ n, t) = arg(detV (ˆ n, t))).

slide-6
SLIDE 6

Introduction

Lefschetz thimbles regularization

◮ Observe that, when only one thimble contributes, one can

rewrite O = Oeiωσ

eiωσ , having defined

f σ =

  • D ˆ

n Zˆ

n

Zσ fˆ n

Zσ =

  • D ˆ

n Zˆ

n, Zˆ n = (2 i λσ i n2 i )

  • dt e−Seff (ˆ

n,t)

n = 1 Zˆ

n (2

i λσ i n2 i )

  • dt f (ˆ

n, t)e−Seff (ˆ

n,t)

→ importance sampling, Pacc(ˆ n′ ← ˆ n) = min

  • 1, Zˆ

n′

n

  • .

◮ Can be generalized to more than one thimble:

O =

  • σ nσZσOeiωσ
  • σ nσZσeiωσ
slide-7
SLIDE 7

On the Lefschetz thimbles structure of the Thirring model

The theory

◮ Let’s consider the 0 + 1-dimensional Thirring model

S = β

  • (1 − cos(zn)) + log detD

detD = 1 2L−1 (cosh(Lˆ µ + i

  • zn) + cosh(L ˆ

m)) , ˆ µ ≡ aµ , ˆ m = asinh(am)

◮ It has been shown before that one thimble is not enough to

capture the full content of the theory (i.e. JHEP1511(2015)078; but see also

JHEP 1605(2016)053)

◮ We wanted to try collecting contributions from different

thimbles with our approach

slide-8
SLIDE 8

On the Lefschetz thimbles structure of the Thirring model

Critical points of the theory (L = 2, β = 1)

◮ To test if we are able to recover the results of the full theory,

we first start from the simple case of L = 2, with a coupling β = 1; the critical points are determined by imposing

∂S ∂zn = βsin(zn) − i sinh(Lµ + i zn) cosh(Lµ + i zn) + cosh(Lm) = 0

The second term depends on the fields only through the sum s ≡ zn, then it must be sin(zn) = sin(z) ∀n and zn can be either z or π − z.

◮ Critical points in the n− = 0 sector for µ = 0 . . . 2

  • 4
  • 2

2 4

  • 1

1 2 3 Real(z) Imag(z)

slide-9
SLIDE 9

On the Lefschetz thimbles structure of the Thirring model

Stokes phenomenon (L = 2, β = 1)

◮ We look for changes in the intersection numbers after Stokes

phenomena

0.5 1 1.5 2

  • 20
  • 10

10 20 mu Imag(S)

(a) SI

0.5 1 1.5 2

  • 4
  • 2

2 4 6 mu Real(S)

(b) SR

Note: due to a simmetry, the two thimbles that enter the thimbles decomposition at µ 0.30 give conjugate contributions.

slide-10
SLIDE 10

On the Lefschetz thimbles structure of the Thirring model

Results from MC (L = 2, β = 1)

◮ How to collect the contributions from more than one thimble? ◮ Last year we proposed to calculate the corrections to the

semiclassical weight as Z G

ˆ n

n −1, it turns out this doesn’t work

very well for this theory.

◮ Nevertheless, we can resort to the method we proposed

previously: O = n0Z0e−iSI

0Oeiω00 + n12Z12e−iSI 12Oeiω1212

n0Z0e−iSI

0eiω00 + n12Z12e−iSI 12eiω1212

≡ Oeiω00 + αOeiω1212 eiω00 + αeiω1212 .

slide-11
SLIDE 11

On the Lefschetz thimbles structure of the Thirring model

Results from MC (L = 2, β = 1)

◮ Fermion number density

0.5 1 1.5 2 0.2 0.4 0.6 0.8 µ n

(c) Results from σ0

0.5 1 1.5 2 0.2 0.4 0.6 0.8 µ n

(d) Results from σ0,i

◮ Condensate

0.5 1 1.5 2 0.2 0.3 0.4 0.5 0.6 0.7 0.8 µ ¯ χχ

(e) Results from σ0

0.5 1 1.5 2 0.2 0.3 0.4 0.5 0.6 0.7 0.8 µ ¯ χχ

(f) Results from σ0,i

slide-12
SLIDE 12

On the Lefschetz thimbles structure of the Thirring model

Numerical integration

◮ At least in this case, where we only have 2 degrees of

freedom, one may also compute expectation values by direct numerical integration.

  • 1e-09 -5e-10

5e-10 1e-09 10 10.5 11 n0 log(Zn)

(g) LogZn vs n0

  • 1e-09 -5e-10

5e-10 1e-09

  • 12
  • 10
  • 8
  • 6
  • 4
  • 2

2 n0 n

(h) n vs n0

◮ Strong dependence of Zn on n0, that is the component of the

initial displacement along the tangent space at the critical point associated to the larger Takagi value. Sharp peaks in some regions of n0. These observations appear to hold also for higher L.

slide-13
SLIDE 13

On the Lefschetz thimbles structure of the Thirring model

Numerical integration (L = 2, β = 1)

◮ Fermion number density

0.5 1 1.5 2 0.2 0.4 0.6 0.8 µ n

(i) Results from σ0

0.5 1 1.5 2 0.2 0.4 0.6 0.8 µ n

(j) Results from σ0,i

◮ Condensate

0.5 1 1.5 2 0.2 0.3 0.4 0.5 0.6 0.7 0.8 µ ¯ χχ

(k) Results from σ0

0.5 1 1.5 2 0.2 0.3 0.4 0.5 0.6 0.7 0.8 µ ¯ χχ

(l) Results from σ0,i

slide-14
SLIDE 14

On the Lefschetz thimbles structure of the Thirring model

Results from MC (L = 4, β = 1)

◮ Fermion number density

0.5 1 1.5 2

  • 0.2

0.2 0.4 0.6 0.8 1 µ n

(m) Results from σ0

0.5 1 1.5 2

  • 0.2

0.2 0.4 0.6 0.8 1 µ n

(n) Results from σ0,i

◮ Condensate

0.5 1 1.5 2 0.2 0.4 0.6 0.8 1 µ ¯ χχ

(o) Results from σ0

0.5 1 1.5 2 0.2 0.4 0.6 0.8 1 µ ¯ χχ

(p) Results from σ0,i

slide-15
SLIDE 15

On the Lefschetz thimbles structure of the Thirring model

Reweighting and Taylor expansion for thimbles

◮ Reweighting:

  • dx O[x] e−β (1−cos(zn))−SF =

=

  • dx O[x]e−(β−β′) (1−cos(zn)) e−β′ (1−cos(zn))−SF

◮ Taylor expansion:

O(µ) = O(µ0) + ∂O ∂µ

  • µ0

(µ − µ0) + . . .

0.5 1 1.5 2 0.2 0.4 0.6 0.8 µ n

slide-16
SLIDE 16

Conclusions

Conclusions

◮ We have studied the (0 + 1)-dimensional Thirring model for

L = 2 and L = 4 with a strong coupling β = 1.

◮ Discrepancies between the analytical solution and the results

from one-thimble simulations seem to disappear after collecting the contribution from the sub-leading thimble.

◮ Very preliminary, we explored reweighting and Taylor

expansion for thimbles.