New Dirac Delta function based methods - Kempf, Jackson, Morales - - PowerPoint PPT Presentation

new dirac delta function based methods kempf jackson
SMART_READER_LITE
LIVE PREVIEW

New Dirac Delta function based methods - Kempf, Jackson, Morales - - PowerPoint PPT Presentation

New Dirac Delta function based methods - Kempf, Jackson, Morales Marco Knipfer Goethe Universitaet Frankfurt knipfer@fias.uni-frankfurt.de July 30, 2014 Marco Knipfer (FIAS) New Dirac Delta July 30, 2014 1 / 23 The paper I am talking about


slide-1
SLIDE 1

New Dirac Delta function based methods - Kempf, Jackson, Morales

Marco Knipfer

Goethe Universitaet Frankfurt knipfer@fias.uni-frankfurt.de

July 30, 2014

Marco Knipfer (FIAS) New Dirac Delta July 30, 2014 1 / 23

slide-2
SLIDE 2

The paper I am talking about

Achim Kempf, David M. Jackson, Alejandro H. Morales New Dirac Delta function based methods with applications to perturbative expansions in quantum field theory and beyond http://arxiv.org/abs/1404.0747

Marco Knipfer (FIAS) New Dirac Delta July 30, 2014 2 / 23

slide-3
SLIDE 3

Overview

1

Introduction New way of writing the Dirac Delta New way of Fourier transforming New integration method

2

Application Quantum Field Theory Application to our work?

3

New delta, blurring, deblurring and QFT

Marco Knipfer (FIAS) New Dirac Delta July 30, 2014 3 / 23

slide-4
SLIDE 4

Fourier Transform Convention

  • g(y) :=

1 √ 2π

  • g(x)eixy dx

(1)

Marco Knipfer (FIAS) New Dirac Delta July 30, 2014 4 / 23

slide-5
SLIDE 5

Fourier Transform Convention

  • g(y) :=

1 √ 2π

  • g(x)eixy dx

(1) g(x) := 1 √ 2π

  • g(y)e−ixy dx

(2)

Marco Knipfer (FIAS) New Dirac Delta July 30, 2014 4 / 23

slide-6
SLIDE 6

New Dirac Delta

δ(x) = 1 2π

  • eixy dy

(3) = 1 2π

  • 1

g(y)g(y)eixy dy (4) = 1 √ 2π 1 g(−i∂x) 1 √ 2π g(y)eixy dy

  • g(x)

(5)

Marco Knipfer (FIAS) New Dirac Delta July 30, 2014 5 / 23

slide-7
SLIDE 7

New Dirac Delta

δ(x) = 1 2π

  • eixy dy

(3) = 1 2π

  • 1

g(y)g(y)eixy dy (4) = 1 √ 2π 1 g(−i∂x) 1 √ 2π g(y)eixy dy

  • g(x)

(5) δ(x) = 1 √ 2π 1 g(−i∂x) g(x), (6)

Marco Knipfer (FIAS) New Dirac Delta July 30, 2014 5 / 23

slide-8
SLIDE 8

New Dirac Delta

δ(x) = 1 2π

  • eixy dy

(3) = 1 2π

  • 1

g(y)g(y)eixy dy (4) = 1 √ 2π 1 g(−i∂x) 1 √ 2π g(y)eixy dy

  • g(x)

(5) δ(x) = 1 √ 2π 1 g(−i∂x) g(x), (6) where g must be a “sufficiently well-behaved” function, e.g. a Gaussian, then

Marco Knipfer (FIAS) New Dirac Delta July 30, 2014 5 / 23

slide-9
SLIDE 9

New Dirac Delta

Example: g is a Gaussian, then δ(x) = 1 √ 2πσ

  • n=0

( −σ

2 ∂2 x)n

n! e−x2/2σ, (7)

Marco Knipfer (FIAS) New Dirac Delta July 30, 2014 6 / 23

slide-10
SLIDE 10

New Dirac Delta

Example: g is a Gaussian, then δ(x) = 1 √ 2πσ

  • n=0

( −σ

2 ∂2 x)n

n! e−x2/2σ, (7) can be truncated → approximation of δ(x).

Marco Knipfer (FIAS) New Dirac Delta July 30, 2014 6 / 23

slide-11
SLIDE 11

New way of Fourier transforming

δ(x) = 1 √ 2π 1 g(−i∂x) g(x) , (8)

Marco Knipfer (FIAS) New Dirac Delta July 30, 2014 7 / 23

slide-12
SLIDE 12

New way of Fourier transforming

δ(x) = 1 √ 2π 1 g(−i∂x) g(x) , (8) can be inverted

  • g(x) =

√ 2πg(−i∂x)δ(x) . (9)

Marco Knipfer (FIAS) New Dirac Delta July 30, 2014 7 / 23

slide-13
SLIDE 13

New way of Fourier transforming

δ(x) = 1 √ 2π 1 g(−i∂x) g(x) , (8) can be inverted

  • g(x) =

√ 2πg(−i∂x)δ(x) . (9)

Benefits?

Maybe easier to compute

Marco Knipfer (FIAS) New Dirac Delta July 30, 2014 7 / 23

slide-14
SLIDE 14

New way of Fourier transforming

δ(x) = 1 √ 2π 1 g(−i∂x) g(x) , (8) can be inverted

  • g(x) =

√ 2πg(−i∂x)δ(x) . (9)

Benefits?

Maybe easier to compute Put in approximation of δ(x) → approximative FT

Marco Knipfer (FIAS) New Dirac Delta July 30, 2014 7 / 23

slide-15
SLIDE 15

Second version of the FT

First one:

  • g(x) =

√ 2πg(−i∂x)δ(x) (10)

Marco Knipfer (FIAS) New Dirac Delta July 30, 2014 8 / 23

slide-16
SLIDE 16

Second version of the FT

First one:

  • g(x) =

√ 2πg(−i∂x)δ(x) (10) Second one (proof easy, see paper):

  • g(y) =

√ 2πeixyδ(i∂x − y)g(x) (11)

Marco Knipfer (FIAS) New Dirac Delta July 30, 2014 8 / 23

slide-17
SLIDE 17

Second version of the FT

First one:

  • g(x) =

√ 2πg(−i∂x)δ(x) (10) Second one (proof easy, see paper):

  • g(y) =

√ 2πeixyδ(i∂x − y)g(x) (11) How can this be used to integrate?

Marco Knipfer (FIAS) New Dirac Delta July 30, 2014 8 / 23

slide-18
SLIDE 18

Second version of the FT

First one:

  • g(x) =

√ 2πg(−i∂x)δ(x) (10) Second one (proof easy, see paper):

  • g(y) =

√ 2πeixyδ(i∂x − y)g(x) (11) How can this be used to integrate? ∞

−∞

g(x) dx = 2π g(0) (12)

Marco Knipfer (FIAS) New Dirac Delta July 30, 2014 8 / 23

slide-19
SLIDE 19

New integration method using the new Dirac Delta

−∞

g(x) dx = 2π g(0) , (13)

Marco Knipfer (FIAS) New Dirac Delta July 30, 2014 9 / 23

slide-20
SLIDE 20

New integration method using the new Dirac Delta

−∞

g(x) dx = 2π g(0) , (13) so by using the above representations of g

  • g(x) dx = 2πg(−i∂x)δ(x)|x=0

(14)

Marco Knipfer (FIAS) New Dirac Delta July 30, 2014 9 / 23

slide-21
SLIDE 21

New integration method using the new Dirac Delta

−∞

g(x) dx = 2π g(0) , (13) so by using the above representations of g

  • g(x) dx = 2πg(−i∂x)δ(x)|x=0

(14)

  • g(x) dx = 2πδ(i∂x)g(x)

(15)

Marco Knipfer (FIAS) New Dirac Delta July 30, 2014 9 / 23

slide-22
SLIDE 22

What do we got so far?

Two representations of FT

  • g(x) =

√ 2πg(−i∂x)δ(x) (16)

  • g(y) =

√ 2πeixyδ(i∂x − y)g(x) (17)

Marco Knipfer (FIAS) New Dirac Delta July 30, 2014 10 / 23

slide-23
SLIDE 23

What do we got so far?

Two representations of FT

  • g(x) =

√ 2πg(−i∂x)δ(x) (16)

  • g(y) =

√ 2πeixyδ(i∂x − y)g(x) (17)

Two representations of integrals from 0 to ∞

  • g(x) dx = 2πg(−i∂x)δ(x)|x=0

(18)

  • g(x) dx = 2πδ(i∂x)g(x)

(19)

Marco Knipfer (FIAS) New Dirac Delta July 30, 2014 10 / 23

slide-24
SLIDE 24

Overview

1

Introduction New way of writing the Dirac Delta New way of Fourier transforming New integration method

2

Application Quantum Field Theory Application to our work?

3

New delta, blurring, deblurring and QFT

Marco Knipfer (FIAS) New Dirac Delta July 30, 2014 11 / 23

slide-25
SLIDE 25

Quantum Field Theory I

Generating functional (applying eqs. for FT) Z[J] =

  • exp
  • iS[φ] + i
  • Jφ dnx
  • Dφ ,

(20)

Marco Knipfer (FIAS) New Dirac Delta July 30, 2014 12 / 23

slide-26
SLIDE 26

Quantum Field Theory I

Generating functional (applying eqs. for FT) Z[J] =

  • exp
  • iS[φ] + i
  • Jφ dnx
  • Dφ ,

(20) translates to Z[J] = NeiS[−iδ/δJ]δ[J] (21) = Nei

  • φJ dnxδ(iδ/δφ − J)eiS[φ]

(22)

Marco Knipfer (FIAS) New Dirac Delta July 30, 2014 12 / 23

slide-27
SLIDE 27

Quantum Field Theory I

Generating functional (applying eqs. for FT) Z[J] =

  • exp
  • iS[φ] + i
  • Jφ dnx
  • Dφ ,

(20) translates to Z[J] = NeiS[−iδ/δJ]δ[J] (21) = Nei

  • φJ dnxδ(iδ/δφ − J)eiS[φ]

(22) Approximation by using approximative δ[J]

Marco Knipfer (FIAS) New Dirac Delta July 30, 2014 12 / 23

slide-28
SLIDE 28

Quantum Field Theory I

Generating functional (applying eqs. for FT) Z[J] =

  • exp
  • iS[φ] + i
  • Jφ dnx
  • Dφ ,

(20) translates to Z[J] = NeiS[−iδ/δJ]δ[J] (21) = Nei

  • φJ dnxδ(iδ/δφ − J)eiS[φ]

(22) Approximation by using approximative δ[J] “Usual pertubative expansion of Z[J] is not convergent and is at best asymptotic”

Marco Knipfer (FIAS) New Dirac Delta July 30, 2014 12 / 23

slide-29
SLIDE 29

Quantum Field Theory I

Generating functional (applying eqs. for FT) Z[J] =

  • exp
  • iS[φ] + i
  • Jφ dnx
  • Dφ ,

(20) translates to Z[J] = NeiS[−iδ/δJ]δ[J] (21) = Nei

  • φJ dnxδ(iδ/δφ − J)eiS[φ]

(22) Approximation by using approximative δ[J] “Usual pertubative expansion of Z[J] is not convergent and is at best asymptotic” This could lead to better convergence

Marco Knipfer (FIAS) New Dirac Delta July 30, 2014 12 / 23

slide-30
SLIDE 30

Quantum Field Theory II

Generating functional (applying eqs. for integration) Z[J] = NeiS[−iδ/δφ]+

  • Jδ/δφ dnxδ[φ]|φ=0

(23) = Nδ[iδ/δφ]eiS[φ]+i

  • Jφ dnx

(24)

Marco Knipfer (FIAS) New Dirac Delta July 30, 2014 13 / 23

slide-31
SLIDE 31

Application to our work

To solve

  • dnp

1 + βpn eixp =

  • f (p)eixp dnp

(25)

Marco Knipfer (FIAS) New Dirac Delta July 30, 2014 14 / 23

slide-32
SLIDE 32

Application to our work

To solve

  • dnp

1 + βpn eixp =

  • f (p)eixp dnp

(25)

The possibilities

1

  • g(x) =

√ 2πg(−i∂x)δ(x)

Marco Knipfer (FIAS) New Dirac Delta July 30, 2014 14 / 23

slide-33
SLIDE 33

Application to our work

To solve

  • dnp

1 + βpn eixp =

  • f (p)eixp dnp

(25)

The possibilities

1

  • g(x) =

√ 2πg(−i∂x)δ(x)

2

  • g(y) =

√ 2πeixpδ(i∂x − y)g(x)

Marco Knipfer (FIAS) New Dirac Delta July 30, 2014 14 / 23

slide-34
SLIDE 34

Application to our work

To solve

  • dnp

1 + βpn eixp =

  • f (p)eixp dnp

(25)

The possibilities

1

  • g(x) =

√ 2πg(−i∂x)δ(x)

2

  • g(y) =

√ 2πeixpδ(i∂x − y)g(x)

3

g(x) dx = 2πg(−i∂x)δ(x)|x=0

Marco Knipfer (FIAS) New Dirac Delta July 30, 2014 14 / 23

slide-35
SLIDE 35

Application to our work

To solve

  • dnp

1 + βpn eixp =

  • f (p)eixp dnp

(25)

The possibilities

1

  • g(x) =

√ 2πg(−i∂x)δ(x)

2

  • g(y) =

√ 2πeixpδ(i∂x − y)g(x)

3

g(x) dx = 2πg(−i∂x)δ(x)|x=0

4

g(x) dx = 2πδ(i∂x)g(x)

Marco Knipfer (FIAS) New Dirac Delta July 30, 2014 14 / 23

slide-36
SLIDE 36

Application to our work

To solve

  • dnp

1 + βpn eixp =

  • f (p)eixp dnp

(25)

The possibilities

1

  • g(x) =

√ 2πg(−i∂x)δ(x)

2

  • g(y) =

√ 2πeixpδ(i∂x − y)g(x)

3

g(x) dx = 2πg(−i∂x)δ(x)|x=0

4

g(x) dx = 2πδ(i∂x)g(x) Since 3-4 follow from 1-2, they should not give any new method for solving eq. (25). Right?

Marco Knipfer (FIAS) New Dirac Delta July 30, 2014 14 / 23

slide-37
SLIDE 37

#1

The Fourier Transformation #1

  • dnp

1 + βpn eixp =

  • f (p)eixp dnp

(26)

Marco Knipfer (FIAS) New Dirac Delta July 30, 2014 15 / 23

slide-38
SLIDE 38

#1

The Fourier Transformation #1

  • dnp

1 + βpn eixp =

  • f (p)eixp dnp

(26) The equation g(x) = √ 2πg(−i∂x)δ(x) has so be extended to n dimensions

  • f (x) = (2π)n/2f (−i∇x)δ(x)

(27) = (2π)n/2 1 1 + β(−i∇x)n δ(x), (28)

Marco Knipfer (FIAS) New Dirac Delta July 30, 2014 15 / 23

slide-39
SLIDE 39

#1

The Fourier Transformation #1

  • dnp

1 + βpn eixp =

  • f (p)eixp dnp

(26) The equation g(x) = √ 2πg(−i∂x)δ(x) has so be extended to n dimensions

  • f (x) = (2π)n/2f (−i∇x)δ(x)

(27) = (2π)n/2 1 1 + β(−i∇x)n δ(x), (28) That’s where we started some time ago . . .

Marco Knipfer (FIAS) New Dirac Delta July 30, 2014 15 / 23

slide-40
SLIDE 40

#2

The Fourier Transformation #2

  • dnp

1 + βpn eixp =

  • f (p)eixp dnp

(29)

Marco Knipfer (FIAS) New Dirac Delta July 30, 2014 16 / 23

slide-41
SLIDE 41

#2

The Fourier Transformation #2

  • dnp

1 + βpn eixp =

  • f (p)eixp dnp

(29)

  • f (x) = (2π)n/2eixpδn(i∇p − x)f (p)

(30) = (2π)n/2eixpδn(i∇p − x) 1 1 + βpn (31) = (2π)n/2eixp ∞ δn (i(nβs) − x) e−(1−βpn) ds , (32)

Marco Knipfer (FIAS) New Dirac Delta July 30, 2014 16 / 23

slide-42
SLIDE 42

#2

The Fourier Transformation #2

  • dnp

1 + βpn eixp =

  • f (p)eixp dnp

(29)

  • f (x) = (2π)n/2eixpδn(i∇p − x)f (p)

(30) = (2π)n/2eixpδn(i∇p − x) 1 1 + βpn (31) = (2π)n/2eixp ∞ δn (i(nβs) − x) e−(1−βpn) ds , (32) and now?

Marco Knipfer (FIAS) New Dirac Delta July 30, 2014 16 / 23

slide-43
SLIDE 43

#2

The Fourier Transformation #2

  • dnp

1 + βpn eixp =

  • f (p)eixp dnp

(29)

  • f (x) = (2π)n/2eixpδn(i∇p − x)f (p)

(30) = (2π)n/2eixpδn(i∇p − x) 1 1 + βpn (31) = (2π)n/2eixp ∞ δn (i(nβs) − x) e−(1−βpn) ds , (32) and now? In eq. (31) use approximation for the δ-function → approx. solution

Marco Knipfer (FIAS) New Dirac Delta July 30, 2014 16 / 23

slide-44
SLIDE 44

#2

The Fourier Transformation #2

  • dnp

1 + βpn eixp =

  • f (p)eixp dnp

(29)

  • f (x) = (2π)n/2eixpδn(i∇p − x)f (p)

(30) = (2π)n/2eixpδn(i∇p − x) 1 1 + βpn (31) = (2π)n/2eixp ∞ δn (i(nβs) − x) e−(1−βpn) ds , (32) and now? In eq. (31) use approximation for the δ-function → approx. solution Maybe in eq. (32) something more can be done. Reminds me of that Fadeev Popov trick.

Marco Knipfer (FIAS) New Dirac Delta July 30, 2014 16 / 23

slide-45
SLIDE 45

Another idea

Take the equation for the Dirac Delta δ(x) = 1 √ 2π 1 g(−i∂x) g(x) , (33) and use g =

1 1+βpn

δn(x) = 1 (2π)n/2 1 g(−i∇x) g(x) (34) = 1 (2π)n/2 (1 + β(−i∇x)n)

  • 1

(2π)n/2 eixp 1 1 + βpn dnp . (35) Rearranging the last line gives

  • eixp

1 + βpn dnp = 1 (2π)n δn(x) − (−i)nβ

  • ∇n

x

  • (ip)n

eixp 1 1 + βpn dnp (36)

Marco Knipfer (FIAS) New Dirac Delta July 30, 2014 17 / 23

slide-46
SLIDE 46

Another idea

  • eixp

1 + βpn dnp = 1 (2π)n δn(x) − (−i)nβ

  • (ip)neixp

1 1 + βpn dnp (37)

Marco Knipfer (FIAS) New Dirac Delta July 30, 2014 18 / 23

slide-47
SLIDE 47

Another idea

  • eixp

1 + βpn dnp = 1 (2π)n δn(x) − (−i)nβ

  • (ip)neixp

1 1 + βpn dnp (37) The integral only gives the energy density, has to be integrated again from 0 to r (spher. coords.).

Marco Knipfer (FIAS) New Dirac Delta July 30, 2014 18 / 23

slide-48
SLIDE 48

Another idea

  • eixp

1 + βpn dnp = 1 (2π)n δn(x) − (−i)nβ

  • (ip)neixp

1 1 + βpn dnp (37) The integral only gives the energy density, has to be integrated again from 0 to r (spher. coords.). Maybe the last integral is somehow easier to solve, but I didn’t go further.

Marco Knipfer (FIAS) New Dirac Delta July 30, 2014 18 / 23

slide-49
SLIDE 49

Overview

1

Introduction New way of writing the Dirac Delta New way of Fourier transforming New integration method

2

Application Quantum Field Theory Application to our work?

3

New delta, blurring, deblurring and QFT

Marco Knipfer (FIAS) New Dirac Delta July 30, 2014 19 / 23

slide-50
SLIDE 50

Blurring a signal

Blurred signal fB fB(y) :=

  • f (x)g(x − y) dx ,

(38) e.g. Gaussian blurring.

Marco Knipfer (FIAS) New Dirac Delta July 30, 2014 20 / 23

slide-51
SLIDE 51

Blurring a signal

Blurred signal fB fB(y) :=

  • f (x)g(x − y) dx ,

(38) e.g. Gaussian blurring. Also possible: Deblurring the signal. Deblurring operator Dg DgfB(y) = 1 √ 2π g(i∂y) fB(y) = f (y) , (39) proof see paper. For Gaussian:

Marco Knipfer (FIAS) New Dirac Delta July 30, 2014 20 / 23

slide-52
SLIDE 52

Blurring a signal

Blurred signal fB fB(y) :=

  • f (x)g(x − y) dx ,

(38) e.g. Gaussian blurring. Also possible: Deblurring the signal. Deblurring operator Dg DgfB(y) = 1 √ 2π g(i∂y) fB(y) = f (y) , (39) proof see paper. For Gaussian: DGauss = e− 1

2a ∂2 y =

  • n=0

1 n!

  • −∂2

y

2a n , (40) can be truncated to get an approximation.

Marco Knipfer (FIAS) New Dirac Delta July 30, 2014 20 / 23

slide-53
SLIDE 53

Euclidean φ4 theory

Generating functional: Z[J] =

  • e
  • dnx(− 1

2 φ(x)(−∆+m2)φ(x)−λφ4(x)+J(x)φ(x))Dφ

(41) in quadratic term substitute φ by derivatives. Z[J] = e− 1

2

  • dnx

δ δJ(x) (−∆+m2) δ δJ(x)

  • Gaussian deblurring operator
  • e
  • dnx(−λφ4(x)+J(x)φ(x))Dφ

(42)

Marco Knipfer (FIAS) New Dirac Delta July 30, 2014 21 / 23

slide-54
SLIDE 54

Euclidean φ4 theory

Generating functional: Z[J] =

  • e
  • dnx(− 1

2 φ(x)(−∆+m2)φ(x)−λφ4(x)+J(x)φ(x))Dφ

(41) in quadratic term substitute φ by derivatives. Z[J] = e− 1

2

  • dnx

δ δJ(x) (−∆+m2) δ δJ(x)

  • Gaussian deblurring operator
  • e
  • dnx(−λφ4(x)+J(x)φ(x))Dφ

(42) Can be expanded in powers of

1 √ λ (strong coupling).

Marco Knipfer (FIAS) New Dirac Delta July 30, 2014 21 / 23

slide-55
SLIDE 55

Euclidean φ4 theory

Generating functional: Z[J] =

  • e
  • dnx(− 1

2 φ(x)(−∆+m2)φ(x)−λφ4(x)+J(x)φ(x))Dφ

(41) in quadratic term substitute φ by derivatives. Z[J] = e− 1

2

  • dnx

δ δJ(x) (−∆+m2) δ δJ(x)

  • Gaussian deblurring operator
  • e
  • dnx(−λφ4(x)+J(x)φ(x))Dφ

(42) Can be expanded in powers of

1 √ λ (strong coupling).

The typical way of pulling φ4 in front of the integral is analytically not justified (see ref [10,11] in paper).

Marco Knipfer (FIAS) New Dirac Delta July 30, 2014 21 / 23

slide-56
SLIDE 56

Marco Knipfer (FIAS) New Dirac Delta July 30, 2014 22 / 23

slide-57
SLIDE 57

The End

Marco Knipfer (FIAS) New Dirac Delta July 30, 2014 23 / 23