Common integrals in quantum field theory

Common integrals in quantum field theory

There are common integrals in quantum field theory that appear repeatedly.[1] These integrals are all variations and generalizations of gaussian integrals to the complex plane and to multiple dimensions. Other integrals can be approximated by versions of the gaussian integral. Fourier integrals are also considered.

Contents

Variations on a simple gaussian integral

Gaussian integral

The first integral, with broad application outside of quantum field theory, is the gaussian integral.

 G \equiv \int_{-\infty}^{\infty} e^{-{1 \over 2} x^2}\,dx

In physics the factor of 1/2 in the argument of the exponential is common.

Note:

 G^2 = \left ( \int_{-\infty}^{\infty} e^{-{1 \over 2} x^2}\,dx \right ) \cdot \left ( \int_{-\infty}^{\infty} e^{-{1 \over 2} y^2}\,dy \right ) = 2\pi \int_{0}^{\infty} r e^{-{1 \over 2} r^2}\,dr = 2\pi \int_{0}^{\infty} e^{- w}\,dw = 2 \pi

Thus we obtain

 \int_{-\infty}^{\infty} e^{-{1 \over 2} x^2}\,dx = \sqrt{2\pi}

Slight generalization of the gaussian integral

 \int_{-\infty}^{\infty} e^{-{1 \over 2} a x^2}\,dx = \sqrt{2\pi \over a}

where we have scaled

 x \rightarrow {x \over \sqrt{a}} .

Integrals of exponents and even powers of x

 \int_{-\infty}^{\infty} x^2 e^{-{1 \over 2} a x^2}\,dx = -2{d\over da} \int_{-\infty}^{\infty} e^{-{1 \over 2} a x^2}\,dx = -2{d\over da} \left ( {2\pi \over a } \right ) ^{1\over 2} = \left ( {2\pi \over a } \right ) ^{1\over 2} {1\over a}

and

 \int_{-\infty}^{\infty} x^4 e^{-{1 \over 2} a x^2}\,dx = \left ( -2{d\over da} \right) \left ( -2{d\over da} \right) \int_{-\infty}^{\infty} e^{-{1 \over 2} a x^2}\,dx = \left ( -2{d\over da} \right) \left ( -2{d\over da} \right) \left ( {2\pi \over a } \right ) ^{1\over 2} = \left ( {2\pi \over a } \right ) ^{1\over 2} {3\over a^2}

In general

 \int_{-\infty}^{\infty} x^{2n} e^{-{1 \over 2} a x^2}\,dx = \left ( {2\pi \over a } \right ) ^{1\over {2}} {1\over a^{n}} \left ( 2n -1 \right ) \left ( 2n -3 \right ) \cdots 5 \cdot 3 \cdot 1 = \left ( {2\pi \over a } \right ) ^{1\over {2}} {1\over a^{n}} \left ( 2n -1 \right )!!

Note that the integrals of exponents and odd powers of x are 0, due to odd symmetry.

Integrals with a linear term in the argument of the exponent

 \int_{-\infty}^{\infty} \exp\left( -{1 \over 2} a x^2 + Jx\right ) dx

This integral can be performed by completing the square.

  \left( -{1 \over 2} a x^2 + Jx\right ) = -{1 \over 2} a \left ( x^2 - { 2 Jx \over a } + { J^2 \over a^2 } - { J^2 \over a^2 } \right ) = -{1 \over 2} a \left ( x -  { J \over a } \right )^2 + { J^2 \over 2a }
 \int_{-\infty}^{\infty} \exp\left( -{1 \over 2} a x^2 + Jx\right ) dx = \exp\left( { J^2 \over 2a } \right ) \int_{-\infty}^{\infty} \exp \left [ -{1 \over 2} a \left ( x -  { J \over a } \right )^2 \right ] dx = \exp\left( { J^2 \over 2a } \right )\int_{-\infty}^{\infty} \exp\left( -{1 \over 2} a x^2 \right ) dx =  \left ( {2\pi \over a } \right ) ^{1\over 2} \exp\left( { J^2 \over 2a }\right )

Integrals with an imaginary linear term in the argument of the exponent

The integral

 \int_{-\infty}^{\infty} \exp\left( -{1 \over 2} a x^2 + iJx\right ) dx  =  \left ( {2\pi \over a } \right ) ^{1\over 2} \exp\left( -{ J^2 \over 2a }\right )

is proportional to the Fourier transform of the gaussian where J is the conjugate variable of x.

By again completing the square we see that the Fourier transform of a gaussian is also a gaussian, but in the conjugate variable. The larger a is, the narrower the gaussian in x and the wider the gaussian in J. This is a demonstration of the uncertainty principle.

Integrals with a complex argument of the exponent

The integral of interest is (for an example of an application see Relation between Schrödinger's equation and the path integral formulation of quantum mechanics)

 \int_{-\infty}^{\infty} \exp\left( {1 \over 2} i a x^2 + iJx\right ) dx.

We assume that a_{ }^{ } may be complex. All other quantities are real.

Completing the square

 \left( {1 \over 2} i a x^2 + iJx\right )  =  {1\over 2} ia \left ( x^2 + {2Jx \over a} + \left ( { J \over a} \right )^2 - \left ( { J \over a} \right )^2 \right ) = -{1\over 2} {a \over i} \left ( x + {J\over a} \right )^2 - { iJ^2 \over a}

By analogy with the previous integrals

 \int_{-\infty}^{\infty} \exp\left( {1 \over 2} i a x^2 + iJx\right ) dx  =  \left ( {2\pi i \over a } \right ) ^{1\over 2} \exp\left( { -iJ^2 \over 2a }\right )

This result is valid as an integration in the complex plane as long as a_{ }^{ } has a positive imaginary part.

Gaussian integrals in higher dimensions

The one-dimensional integrals can be generalized to multiple dimensions.[2]


\int \exp\left( - \frac 1 2 x \cdot A \cdot x +J \cdot x \right) d^nx
=
\sqrt{\frac{(2\pi)^n}{\det A}} \exp \left( {1\over 2} J \cdot A^{-1} \cdot J \right)

Here A is a real symmetric matrix.

This integral is performed by diagonalization of A with an orthogonal transformation


D_{ }^{ } = O^{-1}  A  O = O^T  A  O

where D is a diagonal matrix and O is an orthogonal matrix. This decouples the variables and allows the integration to be performed as n_{ }^{ } one-dimensional integrations.

This is best illustrated with a two-dimensional example.

Example: Simple gaussian integration in two dimensions

The gaussian integral in two dimensions is


\int \exp\left( - \frac 1 2 A_{ij} x^i x^j \right) d^2x
=
\sqrt{\frac{(2\pi)^2}{\det A}}

where A is a two-dimensional symmetric matrix with components specified as

 A = \bigl[ \begin{smallmatrix}
  a&c\\ c&b
\end{smallmatrix} \bigr]

and we have used the Einstein summation convention.

Diagonalize the matrix

The first step is to diagonalize the matrix.[3] Note that

 A_{ij} x^i x^j \equiv x^T  A  x = x^T  \left( O  O^T\right)  A  \left( O  O^T\right)  x  = \left( x^T  O \right) \left(  O^T  A  O \right) \left(  O^T  x \right)

where, since A is a real symmetric matrix, we can choose O to be an orthogonal matrix, and hence also a unitary matrix.

We choose O such that

 D \equiv   O^T  A  O

is diagonal.

O can be obtained from the eigenvectors of A.

Eigenvalues of A

To find the eigenvectors of A one first finds the eigenvalues λ of A given by

  \bigl[ \begin{smallmatrix}
  a&c\\ c&b
\end{smallmatrix} \bigr] 
\bigl[ \begin{smallmatrix}
  u\\ v
\end{smallmatrix} \bigr] =
\lambda \bigl[ \begin{smallmatrix}
  u\\ v
\end{smallmatrix} \bigr]
    .

The eigenvalues are solutions of the characteristic polynomial

 \left( a - \lambda \right) \left( b-\lambda\right) -c^2 = 0

which are

  \lambda_{\pm} = {1\over 2}\left( a+b\right) \pm {1\over 2}\sqrt{ \left(a-b\right)^2+4c^2} .
Eigenvectors of A

Substitution of the eigenvalues back into the eigenvector equation yields

 v = -{ \left( a - \lambda_{\pm} \right)u \over c }

or

 v = -{cu \over  \left( b - \lambda_{\pm}  \right)}.

From the characteristic equation we know

 { \left( a - \lambda_{\pm} \right)\over c } = {c \over  \left( b - \lambda_{\pm}  \right)}.

Also note

 { \left( a - \lambda_{\pm} \right)\over c } = -{\left( b - \lambda_{\mp}  \right) \over  c   } .

The eigenvectors can be written

  
   \begin{bmatrix}
  {1\over \eta}\\ - \left( {  a - \lambda_{-}  \over c\eta } \right)
\end{bmatrix}

and

  
   \begin{bmatrix}
  - \left( {  b - \lambda_{+}  \over c\eta } \right) \\ { 1\over \eta} 
\end{bmatrix}

for the two eigenvectors. Here η is a normalizing factor given by

  
   \eta = \sqrt{ 1 + \left( {  a - \lambda_{-}  \over c } \right)^2 } = \sqrt{ 1 + \left( {  b - \lambda_{+}  \over c } \right)^2 }

    .

It is easily verified that the two eigenvectors are orthogonal to each other.

Construction of the orthogonal matrix

The orthogonal matrix is constructed by assigning the normalized eigenvectors as columns in the orthogonal matrix

  
O =  \begin{bmatrix}
  {1\over \eta } & -\left( {  b - \lambda_{+}  \over c \eta }\right) \\ -\left( {  a - \lambda_{-}  \over c \eta }\right) & {1\over \eta }
\end{bmatrix} 
    .

Note that the detrminant of O is equal to one.

If we define

 \sin \left( \theta \right) \equiv  -\left( {  a - \lambda_{-}  \over c \eta }\right)

then the orthogonal matrix can be written

  
O =  \begin{bmatrix}
  \cos \left( \theta \right) & -\sin \left( \theta \right) \\ \sin \left( \theta \right) & \cos \left( \theta \right)
\end{bmatrix}

which is simply a rotation of the eigenvectors.

The inverse is

  
O^{-1} = O^T =  \begin{bmatrix}
  \cos \left( \theta \right) & \sin \left( \theta \right) \\ -\sin \left( \theta \right) & \cos \left( \theta \right)
\end{bmatrix}
Diagonal matrix

The diagonal matrix becomes

  
D = O^T A O =  \begin{bmatrix}
 \lambda_{-}&0\\ 0 & \lambda_{+}
\end{bmatrix}

with eigenvectors

  
  \begin{bmatrix}
 1\\ 0
\end{bmatrix}

and

  
  \begin{bmatrix}
 0\\ 1
\end{bmatrix}
Numerical example
  
A =  \begin{bmatrix}
 2&1\\ 1 & 1
\end{bmatrix}

The eigenvalues are

  
\lambda_{\pm}  =  {3\over 2} \pm {\sqrt{ 5} \over 2}
 
    .

The eigenvectors are

  {1\over \eta}
  \begin{bmatrix}
  1\\ { -{1\over 2} - {\sqrt{ 5} \over 2}   }
\end{bmatrix}

and

  {1\over \eta}
  \begin{bmatrix}
  { {1\over 2} + {\sqrt{ 5} \over 2 } } \\  1  
\end{bmatrix}

where

  \eta = \sqrt{ {5\over 2} + {\sqrt{5}\over 2} }
    .

The orthogonal vector is

  O =
  \begin{bmatrix}
  {1\over \eta} & {1\over \eta} \left({ {1\over 2} + {\sqrt{ 5} \over 2}   }\right) \\ {1\over \eta} \left({ -{1\over 2} - {\sqrt{ 5} \over 2}   }\right) & {1\over \eta}
\end{bmatrix} 
    .

It is easily verified that the detrminant of O is 1.

The inverse of O is

  O =
  \begin{bmatrix}
  {1\over \eta} & {1\over \eta} \left({ -{1\over 2} - {\sqrt{ 5} \over 2}   }\right) \\ {1\over \eta} \left({ {1\over 2} + {\sqrt{ 5} \over 2}   }\right) & {1\over \eta}
\end{bmatrix} 
    .


The diagonal matrix becomes

  D = O^TAO = 
  \begin{bmatrix}
  \lambda_{-} &0\\ 0 & \lambda_{+}
\end{bmatrix} =

  \begin{bmatrix}
  \left({ {3\over 2} - {\sqrt{ 5} \over 2}   }\right)& 0\\ 0 & \left({ {3\over 2} + {\sqrt{ 5} \over 2}   }\right)
\end{bmatrix}

with the eigenvectors

  
  \begin{bmatrix}
  1\\ 0 
\end{bmatrix}

and

   \begin{bmatrix}
  0\\ 1 
\end{bmatrix}  
    .

Rescale the variables and integrate

With the diagonalization the integral can be written

   
\int \exp\left( - \frac 1 2 x^T A x \right) d^2x
= \int \exp\left( - \frac 1 2 \sum_{j=1}^2 \lambda_{j} y_j^2 \right) d^2x

where

   y_{ }^{ } = O^T x .

Since the coordinate transformation is simply a rotation of coordinates the Jacobian determinant of the transformation is one yielding

   dy_{ }^2  = dx^2 .

The integrations can now be performed.

   
\int \exp\left( - \frac 1 2 x^T A x \right) d^2x
= \int \exp\left( - \frac 1 2 \sum_{j=1}^2 \lambda_{j} y_j^2 \right) d^2y
= \prod_{j=1}^2 \left( { 2\pi \over \lambda_j } \right)^{1\over 2}
= \left( { ( 2\pi )^2  \over \prod_{j=1}^2 \lambda_j } \right)^{1\over 2}
= \left( { ( 2\pi )^2  \over \det{ \left( O^{-1}AO \right)}  } \right)^{1\over 2}
= \left( { ( 2\pi )^2  \over \det{ \left( A \right)}  } \right)^{1\over 2}

which is the advertised solution.

Integrals with complex and linear terms in multiple dimensions

With the two-dimensional example it is now easy to see the generalization to the complex plane and to multiple dimensions.

Integrals with a linear term in the argument


\int \exp\left( - \frac 1 2 x \cdot A \cdot x +J \cdot x \right) d^nx
=
\sqrt{\frac{(2\pi)^n}{\det A}} \exp \left( {1\over 2} J \cdot A^{-1} \cdot J \right)


Integrals with an imaginary linear term


\int \exp\left( - \frac 1 2 x \cdot A \cdot x +iJ \cdot x \right) d^nx
=
\sqrt{\frac{(2\pi)^n}{\det A}} \exp \left( -{1\over 2} J \cdot A^{-1} \cdot J \right)

Integrals with a complex quadratic term


\int \exp\left(  \frac i 2 x \cdot A \cdot x +iJ \cdot x \right) d^nx
=
\sqrt{\frac{(2\pi i)^n}{\det A}} \exp \left( -{i\over 2} J \cdot A^{-1} \cdot J \right)

Integrals with differential operators in the argument

As an example consider the integral[4]


\int \exp\left[  \int d^4x \left ( -\frac 1 2 \varphi  \hat A  \varphi +  J  \varphi \right) \right ] D\varphi

where  \hat A is a differential operator with  \varphi and J functions of spacetime, and Dφ indicates integration over all possible paths. In analogy with the matrix version of this integral the solution is


\int \exp\left( - \frac 1 2 \varphi  \hat A  \varphi +J  \varphi \right) D\varphi \; \propto \;
\exp \left( {1\over 2} \int d^4x \; d^4y J\left ( x \right ) D\left ( x - y \right )  J\left( y \right )  \right)

where


\hat A D\left ( x - y \right ) = \delta^4 \left ( x - y \right )

and  D\left ( x - y \right ) , called the propagator, is the inverse of  \hat A , and  \delta^4 \left ( x - y \right ) is the Dirac delta function.

Similar arguments yield


\int \exp\left[  \int d^4x \left ( -\frac 1 2 \varphi  \hat A  \varphi +  i J  \varphi \right) \right ] D\varphi \; \propto \;
\exp \left( - { 1\over 2} \int d^4x \; d^4y J\left ( x \right ) D\left ( x - y \right )  J\left( y \right )  \right)

and


\int \exp\left[ i \int d^4x \left ( \frac 1 2 \varphi  \hat A  \varphi +  J  \varphi \right) \right ] D\varphi \; \propto \;
\exp \left( { i\over 2} \int d^4x \; d^4y J\left ( x \right ) D\left ( x - y \right )  J\left( y \right )  \right)
.

See Path-integral formulation of virtual-particle exchange for an application of this integral.

Integrals that can be approximated by the method of steepest descent

In quantum field theory n-dimensional integrals of the form

  \int_{-\infty}^{\infty} \exp\left( -{1 \over \hbar} f\left( q \right) \right ) d^nq

appear often. Here \hbar is the reduced Planck's constant and f is a function with a positive minimum at q = q0. These integrals can be approximated by the method of steepest descent.

For small values of Planck's constant, f can be expanded about its minimum

  \int_{-\infty}^{\infty} \exp\left[ -{1 \over \hbar} \left( f\left( q_0 \right)  + {1\over 2} \left( q-q_0\right)^2f^{\prime \prime} \left( q-q_0\right) + \cdots \right ) \right] d^nq.

Here  f^{\prime \prime} is the n by n matrix of second derivatives evaluated at the minimum of the function.

If we neglect higher order terms this integral can be integrated explicitly.

  \int_{-\infty}^{\infty} \exp\left[ -{1 \over \hbar} \left( f\left( q \right)   \right ) \right] d^nq \approx
\exp\left[ -{1 \over \hbar} \left( f\left( q_0 \right)   \right ) \right] \sqrt{ (2 \pi \hbar )^n \over \det f^{\prime \prime} } .

Integrals that can be approximated by the method of stationary phase

A common integral is a path integral of the form

  \int \exp\left( {i \over \hbar}  S\left( q, \dot q \right) \right ) Dq

where  S\left( q, \dot q \right) is the classical action and the integral is over all possible paths that a particle may take. In the limit of small  \hbar the integral can be evaluated in the stationary phase approximation. In this approximation the integral is over the path in which the action is a minimum. Therefore, this approximation recovers the classical limit of mechanics.

Fourier integrals

Dirac delta function

The Dirac delta function in spacetime can be written as a Fourier transform[5]

  \int { d^4 k \over \left ( 2 \pi \right ) ^4 } \; {\exp \left ( ik \left ( x-y \right) \right ) = \delta^4 \left ( x-y \right )  } .

In general, for any dimension N

  \int { d^N k \over \left ( 2 \pi \right ) ^N } \; {\exp \left ( ik \left ( x-y \right) \right ) = \delta^N \left ( x-y \right )  } .

Fourier integrals of forms of the Coulomb potential

Laplacian of 1/r

While not an integral, the identity in three-dimensional Euclidean space

 
 -{1 \over 4\pi} \nabla^2 \left( {1 \over r} \right) 
= \delta \left( \mathbf r \right)

where

 
 r^2 =  \mathbf r \cdot \mathbf r

is a consequence of Gauss's theorem and can be used to derive integral identities. For an example see Longitudinal and transverse vector fields.

This identity implies that the Fourier integral representation of 1/r is

  \int { d^3 k \over \left ( 2 \pi \right ) ^3 } { \exp \left ( i\mathbf k \cdot \mathbf r \right) \over k^2  } = {1 \over 4 \pi r }   .

Yukawa Potential: The Coulomb potential with mass

The Yukawa potential in three dimensions can be represented as an integral over a Fourier transform[6]

  \int { d^3 k \over \left ( 2 \pi \right ) ^3 } { \exp \left ( i\mathbf k \cdot \mathbf r \right) \over k^2 +m^2 } = {e^{  - m r } \over 4 \pi r }

where

  r^2 =  \mathbf r \cdot \mathbf r

and

  k^2 =  \mathbf k \cdot \mathbf k .

See Static forces and virtual-particle exchange for an application of this integral.

In the small m limit the integral reduces to

   {1 \over 4 \pi r }   .

To derive this result note:

  \int { d^3 k \over \left ( 2 \pi \right ) ^3 } { \exp \left ( i \mathbf k \cdot \mathbf r  \right ) \over k^2 +m^2 } = 
\int_0^{\infty} {k^2 dk \over \left ( 2 \pi \right )^2 } \int_{-1}^{1} du {\exp\left( ikru \right) \over k^2 + m^2}

= 
{2\over r} \int_0^{\infty} {k dk \over \left ( 2 \pi \right )^2 }  {\sin\left( kr \right) \over k^2 + m^2} =
{1\over i r} \int_{-\infty}^{\infty} {k dk \over \left ( 2 \pi \right )^2 }  {\exp\left( ikr \right) \over k^2 + m^2}

= 
{1\over i r} \int_{-\infty}^{\infty} {k dk \over \left ( 2 \pi \right )^2 }  {\exp\left( ikr \right) \over \left(k + i m \right)\left(k - i m \right)}
=
{1\over i r} { 2\pi i \over \left( 2 \pi \right)^2 } {im \over 2 i m} \exp \left( -m r \right)

Modified Coulomb potential with mass

  \int { d^3 k \over \left ( 2 \pi \right ) ^3 } 
\left( \mathbf{\hat k}\cdot \mathbf{\hat r}\right)^2
{ \exp \left ( i\mathbf \mathbf k \cdot \mathbf r \right ) \over k^2 +m^2 } = {e^{  - m r } \over 4 \pi r }  \left\{ 1+  {2\over mr}  
-  {2 \over \left( mr \right)^2 } \left( e^{mr} -1 \right)  \right \}

where the hat indicates a unit vector in three dimensional space.

In the small m limit the integral goes to zero.

To derive this result note:

  \int { d^3 k \over \left ( 2 \pi \right ) ^3 } \left( \mathbf{\hat k}\cdot \mathbf{\hat r}\right)^2
 { \exp \left ( i\mathbf \mathbf k \cdot \mathbf r \right ) \over k^2 +m^2 } = 
\int_0^{\infty} {k^2 dk \over \left ( 2 \pi \right )^2 } \int_{-1}^{1} du u^2 {\exp\left( ikru \right) \over k^2 + m^2}

= 
{2} \int_0^{\infty} {k^2 dk \over \left ( 2 \pi \right )^2 }  {1 \over k^2 + m^2} 
 
\left\{   {1\over kr }  \sin\left( kr \right) + 2 {1\over \left(kr\right)^2 } \cos\left( kr \right) 
- 2 {1\over \left(kr\right)^3 } \sin\left( kr \right) \right \}

=
{\exp \left( -m r \right) \over 4\pi r}

\left\{   g\left( mr\right)  
  \right \}

where


 g\left( mr\right) 
= 1+{2\over mr}-{2 \over \left( mr \right)^2 } \left( e^{mr} -1 \right)

.

Note that in the limit that


 m\rightarrow 0

that


 g\left( mr\right)  \rightarrow 1
.

Longitudinal potential with mass

  \int { d^3 k \over \left ( 2 \pi \right ) ^3 } 
 \; \mathbf{\hat k} \mathbf{\hat k} \;
{ \exp \left ( i\mathbf \mathbf k \cdot \mathbf r \right ) \over k^2 +m^2 } 
= 
{1\over 2} {e^{  - m r } \over 4 \pi r }  \left[ \mathbf 1 - \mathbf{\hat r} \mathbf{\hat r} \right]
+
 {1\over 2} {e^{  - m r } \over 4 \pi r }  \left\{ 1+  {2\over mr}  
-  {2 \over \left( mr \right)^2 } \left( e^{mr} -1 \right)  \right \}
  \left[\mathbf 1 + \mathbf{\hat r} \mathbf{\hat r}\right]

where the hat indicates a unit vector in three dimensional space.

In the small m limit the integral reduces to

   {1\over 2} {1 \over 4 \pi r }  \left[ \mathbf 1 - \mathbf{\hat r} \mathbf{\hat r} \right] .

To derive this result note:

  \int { d^3 k \over \left ( 2 \pi \right ) ^3 }  \mathbf{\hat k} \mathbf{\hat k}
 { \exp \left ( i\mathbf k \cdot \mathbf r \right ) \over k^2 +m^2 } 
=
\int { d^3 k \over \left ( 2 \pi \right ) ^3 } 
\left[
\left( \mathbf{\hat k}\cdot \mathbf{\hat r}\right)^2\mathbf{\hat r} \mathbf{\hat r}
+ \left( \mathbf{\hat k}\cdot \mathbf{\hat \theta}\right)^2\mathbf{\hat \theta} \mathbf{\hat \theta}
+ \left( \mathbf{\hat k}\cdot \mathbf{\hat \phi}\right)^2\mathbf{\hat \phi} \mathbf{\hat \phi}
\right]
 { \exp \left ( i\mathbf k \cdot \mathbf r \right ) \over k^2 +m^2 }

where the cross terms vanish. The integral can be written


= 
 {1 \over 4 \pi r } \exp \left ( - m r \right ) \left\{ 1+   {2\over mr }  
-  {2\over \left(mr\right)^2 } \left( e^{mr} -1 \right) \right \}
\left\{\mathbf 1  - {1\over 2} \left[\mathbf 1 - \mathbf{\hat r} \mathbf{\hat r}\right] \right\}
+
\int_0^{\infty} {k^2 dk \over \left ( 2 \pi \right )^2 } \int_{-1}^{1} du  {\exp\left( ikru \right) \over k^2 + m^2} 
{1\over 2}  \left[ \mathbf 1  - \mathbf{\hat r} \mathbf{\hat r} \right]

= 
{1\over 2} {e^{  - m r } \over 4 \pi r }  \left[ \mathbf 1 - \mathbf{\hat r} \mathbf{\hat r} \right]
+
 {e^{  - m r } \over 4 \pi r }  \left\{ 1+   {2\over mr }  
-  {2\over \left(mr\right)^2 } \left( e^{mr} -1 \right)  \right \}
\left\{ {1\over 2} \left[\mathbf 1 + \mathbf{\hat r} \mathbf{\hat r}\right] \right\}

.

Transverse potential with mass

  \int { d^3 k \over \left ( 2 \pi \right ) ^3 } 
 \; \left[ \mathbf 1 - \mathbf{\hat k} \mathbf{\hat k} \right] \;
{ \exp \left ( i\mathbf \mathbf k \cdot \mathbf r \right ) \over k^2 +m^2 } 
= 

 {1\over 2} {e^{  - m r } \over 4 \pi r }  \left\{    
  {2 \over \left( mr \right)^2  } \left( e^{mr} -1 \right) -  {2\over mr} \right \}
\left[\mathbf 1 + \mathbf{\hat r} \mathbf{\hat r}\right]

In the small mr limit the integral goes to

  

 {1\over 2} {1 \over 4 \pi r }  
\left[\mathbf 1 + \mathbf{\hat r} \mathbf{\hat r}\right]
  .

For large distance, the integral falls off as the inverse cube of r

  

  {1\over 4 \pi m^2r^3 }  
\left[\mathbf 1 + \mathbf{\hat r} \mathbf{\hat r}\right]
  .

For applications of this integral see Darwin Lagrangian and Darwin interaction in a vacuum.

Angular integration in cylindrical coordinates

There are two important integrals. The angular integration of an exponential in cylindrical coordinates can be written in terms of Bessel functions of the first kind[7][8]

  
 \int_0^{2 \pi} {d\varphi \over 2 \pi}  \exp\left( i p \cos\left( \varphi \right) \right)
=
 J_0 \left( p \right)

and

  
 \int_0^{2 \pi} {d\varphi \over 2 \pi} \cos\left( \varphi \right) \exp\left( i p \cos\left( \varphi \right) \right)
=
i J_1 \left( p \right)
  .

For applications of these integrals see Magnetic interaction between current loops in a simple plasma or electron gas.

Bessel functions

Integration of the cylindrical propagator with mass

First power of a Bessel function

  
\int_0^{\infty} {k\; dk \over k^2 +m^2} J_0 \left( kr \right)
=
K_0 \left( mr \right)
  .

See Abramowitz and Stegun.[9]

For mr < < 1, we have[10]


K_0 \left( mr \right) \rightarrow -\ln \left( {mr \over 2}\right) + 0.5772
.

For an application of this integral see Two line charges embedded in a plasma or electron gas.

Square of Bessel Functions

The integration of the propagator in cylindrical coordinates is[11]

  
\int_0^{\infty} {k\; dk \over k^2 +m^2} J_1^2 \left( kr \right)
=
I_1 \left( mr \right)K_1 \left( mr \right)
  .

For small mr the integral becomes

  
\int_o^{\infty} {k\; dk \over k^2 +m^2} J_1^2 \left( kr \right)
\rightarrow
{1\over 2 }\left[ 1 - {1\over 8} \left(mr\right)^2 \right]
  .

For large mr the integral becomes

  
\int_o^{\infty} {k\; dk \over k^2 +m^2} J_1^2 \left( kr \right)
\rightarrow
{1\over 2}\;\left( {1\over mr}\right)
  .

For applications of this integral see Magnetic interaction between current loops in a simple plasma or electron gas.

In general

  
\int_0^{\infty} {k\; dk \over k^2 +m^2} J_{\nu}^2 \left( kr \right)
=
I_{\nu} \left( mr \right)K_{\nu} \left( mr \right)
\;\;\; \Re \;{\nu} > -1
  .

Integration over a magnetic wave function

The two-dimensional integral over a magnetic wave function is[12]

  
{2 a^{2n+2}\over n!}
\int_0^{\infty} { dr }\;r^{2n+1}\exp\left( -a^2 r^2\right) J_{0} \left( kr \right)
=
M\left( n+1, 1, -{k^2 \over 4a^2}\right)
  .

Here, M is a confluent hypergeometric function. For an application of this integral see Charge density spread over a wave function.

References

  1. ^ A. Zee (2003). Quantum Field Theory in a Nutshell. Princeton University. ISBN 0-691-01019-6.  pp. 13-15
  2. ^ Frederick W. Byron and Robert W. Fuller (1969). Mathematics of Classical and Quantum Physics. Addison-Wesley. ISBN 0-201-00746-2. 
  3. ^ Herbert S. Wilf (1978). Mathematics for the Physical Sciences. Dover. ISBN 0-486-63635-6. 
  4. ^ Zee, pp. 21-22.
  5. ^ Zee, p. 23.
  6. ^ Zee, p. 26, 29.
  7. ^ I. S. Gradshteyn and I. M. Ryzhik (1965). Tables of Integrals, Seies, and Products. Academic Press. ISBN 65-29097.  p. 402
  8. ^ Jackson, John D. (1998). Classical Electrodynamics (3rd ed.). Wiley. ISBN 047130932X.  p. 113
  9. ^ M. Abramowitz and I. Stegun (1965). Handbook of Mathematical Functions. Dover. ISBN 486-61272-4.  Section 11.4.44
  10. ^ Jackson, p. 116
  11. ^ I. S. Gradshteyn and I. M. Ryzhik, p. 679
  12. ^ Abramowitz and Stegun, Section 11.4.28

See also

  • Relation between Schrödinger's equation and the path integral formulation of quantum mechanics

Wikimedia Foundation. 2010.

Игры ⚽ Поможем сделать НИР

Look at other dictionaries:

  • Scalar field theory — In theoretical physics, scalar field theory can refer to a classical or quantum theory of scalar fields. A field which is invariant under any Lorentz transformation is called a scalar , in contrast to a vector or tensor field. The quanta of the… …   Wikipedia

  • Quantum mechanics — For a generally accessible and less technical introduction to the topic, see Introduction to quantum mechanics. Quantum mechanics …   Wikipedia

  • Field electron emission — It is requested that a diagram or diagrams be included in this article to improve its quality. For more information, refer to discussion on this page and/or the listing at Wikipedia:Requested images. Field emission (FE) (also known as field… …   Wikipedia

  • Superstring theory — is an attempt to explain all of the particles and fundamental forces of nature in one theory by modelling them as vibrations of tiny supersymmetric strings. It is considered one of the most promising candidate theories of quantum gravity.… …   Wikipedia

  • Quantum Monte Carlo — is a large class of computer algorithms that simulate quantum systems with the idea of solving the many body problem. They use, in one way or another, the Monte Carlo method to handle the many dimensional integrals that arise. Quantum Monte Carlo …   Wikipedia

  • Stationary phase approximation — In mathematics, the stationary phase approximation is a basic principle of asymptotic analysis, applying to oscillatory integrals taken over n dimensional space Rn where i is the imaginary unit. Here f and g are real valued smooth functions. The… …   Wikipedia

  • Darwin Lagrangian — The Darwin Lagrangian (named after Charles Galton Darwin, grandson of the biologist) describes the interaction to order between two charged particles in a vacuum and is given by[1] where the free particle Lagrangian is …   Wikipedia

  • Renormalization — Quantum field theory (Feynman diagram) …   Wikipedia

  • Path integral formulation — This article is about a formulation of quantum mechanics. For integrals along a path, also known as line or contour integrals, see line integral. The path integral formulation of quantum mechanics is a description of quantum theory which… …   Wikipedia

  • Complex number — A complex number can be visually represented as a pair of numbers forming a vector on a diagram called an Argand diagram, representing the complex plane. Re is the real axis, Im is the imaginary axis, and i is the square root of –1. A complex… …   Wikipedia

Share the article and excerpts

Direct link
Do a right-click on the link above
and select “Copy Link”