Existen una serie de integrales típicas en teoría cuántica de campos que aparecen continuamente al aplicarla a los problemas concretos de la teoría.[ 1] Estas integrales son todas variaciones y generalizaciones de integrales gaussianas en el plano complejo y en varias (muchas, o incluso infinitas) dimensiones. Otro conjunto de integrales pueden ser aproximadas por integrales gaussianas . También se presentan transformaciones de Fourier .
Integral gaussiana en una dimensión
editar
Integral Gaussiana
editar
La primera integral aquí presentada, de uso común también en aplicaciones diferentes a la teoría cuántica de campos , es la integral gaussiana , que es la integral sobre todo su dominio de la campana de Gauss .
G
≡
∫
−
∞
∞
e
−
1
2
x
2
d
x
{\displaystyle G\equiv \int _{-\infty }^{\infty }e^{-{1 \over 2}x^{2}}\,dx}
En Física, es habitual usar un factor 1/2 en el argumento de la exponencial.
La forma usual de resolver esta integral es intentando resolver su versión bidimensional haciendo uso de la propiedad de factorización de la integral respecto a la suma de argumentos y de un cambio de variables de cartesianas a polares:
G
2
=
(
∫
−
∞
∞
e
−
1
2
x
2
d
x
)
⋅
(
∫
−
∞
∞
e
−
1
2
y
2
d
y
)
=
2
π
∫
0
∞
r
e
−
1
2
r
2
d
r
=
2
π
∫
0
∞
e
−
w
d
w
=
2
π
.
{\displaystyle 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 }re^{-{1 \over 2}r^{2}}\,dr=2\pi \int _{0}^{\infty }e^{-w}\,dw=2\pi .}
Con lo que se obtiene el resultado buscado
∫
−
∞
∞
e
−
1
2
x
2
d
x
=
2
π
.
{\displaystyle \int _{-\infty }^{\infty }e^{-{1 \over 2}x^{2}}\,dx={\sqrt {2\pi }}.}
Una minúscula generalización de la integral gaussiana
editar
∫
−
∞
∞
e
−
1
2
a
x
2
d
x
=
2
π
a
{\displaystyle \int _{-\infty }^{\infty }e^{-{1 \over 2}ax^{2}}\,dx={\sqrt {2\pi \over a}}}
donde hemos reescalado la variable
x
→
x
a
{\displaystyle x\rightarrow {x \over {\sqrt {a}}}}
.
Integrales de gaussianas con potencias pares de x
editar
∫
−
∞
∞
x
2
e
−
1
2
a
x
2
d
x
=
−
2
d
d
a
∫
−
∞
∞
e
−
1
2
a
x
2
d
x
=
−
2
d
d
a
(
2
π
a
)
1
2
=
(
2
π
a
)
1
2
1
a
{\displaystyle \int _{-\infty }^{\infty }x^{2}e^{-{1 \over 2}ax^{2}}\,dx=-2{d \over da}\int _{-\infty }^{\infty }e^{-{1 \over 2}ax^{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}}
y
∫
−
∞
∞
x
4
e
−
1
2
a
x
2
d
x
=
(
−
2
d
d
a
)
(
−
2
d
d
a
)
∫
−
∞
∞
e
−
1
2
a
x
2
d
x
=
(
−
2
d
d
a
)
(
−
2
d
d
a
)
(
2
π
a
)
1
2
=
(
2
π
a
)
1
2
3
a
2
{\displaystyle \int _{-\infty }^{\infty }x^{4}e^{-{1 \over 2}ax^{2}}\,dx=\left(-2{d \over da}\right)\left(-2{d \over da}\right)\int _{-\infty }^{\infty }e^{-{1 \over 2}ax^{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}}}
En general:
∫
−
∞
∞
x
2
n
e
−
1
2
a
x
2
d
x
=
(
2
π
a
)
1
2
1
a
n
(
2
n
−
1
)
(
2
n
−
3
)
⋯
5
⋅
3
⋅
1
=
(
2
π
a
)
1
2
1
a
n
(
2
n
−
1
)
!
!
{\displaystyle \int _{-\infty }^{\infty }x^{2n}e^{-{1 \over 2}ax^{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)!!}
Nótese que, en el caso en el que tengamos potencias impares de x , la integral se anula debido a la simetría impar del argumento, por lo que
∫
−
∞
∞
x
2
n
+
1
e
−
1
2
a
x
2
d
x
=
0
{\displaystyle \int _{-\infty }^{\infty }x^{2n+1}e^{-{1 \over 2}ax^{2}}\,dx=0}
Integrales con un término lineal en el argumento de la exponencial
editar
∫
−
∞
∞
exp
(
−
1
2
a
x
2
+
J
x
)
d
x
{\displaystyle \int _{-\infty }^{\infty }\exp \left(-{1 \over 2}ax^{2}+Jx\right)dx}
Esta integral se puede resolver completando el cuadrado.
(
−
1
2
a
x
2
+
J
x
)
=
−
1
2
a
(
x
2
−
2
J
x
a
+
J
2
a
2
−
J
2
a
2
)
=
−
1
2
a
(
x
−
J
a
)
2
+
J
2
2
a
{\displaystyle \left(-{1 \over 2}ax^{2}+Jx\right)=-{1 \over 2}a\left(x^{2}-{2Jx \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}}
∫
−
∞
∞
exp
(
−
1
2
a
x
2
+
J
x
)
d
x
=
exp
(
J
2
2
a
)
∫
−
∞
∞
exp
[
−
1
2
a
(
x
−
J
a
)
2
]
d
x
=
exp
(
J
2
2
a
)
∫
−
∞
∞
exp
(
−
1
2
a
x
2
)
d
x
=
(
2
π
a
)
1
2
exp
(
J
2
2
a
)
{\displaystyle {\begin{aligned}&\int _{-\infty }^{\infty }\exp \left(-{1 \over 2}ax^{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\\[8pt]&=\exp \left({J^{2} \over 2a}\right)\int _{-\infty }^{\infty }\exp \left(-{1 \over 2}ax^{2}\right)\,dx=\left({2\pi \over a}\right)^{1 \over 2}\exp \left({J^{2} \over 2a}\right)\end{aligned}}}
Integrales con un término lineal imaginario en el argumento de la exponencial
editar
La integral
∫
−
∞
∞
exp
(
−
1
2
a
x
2
+
i
J
x
)
d
x
=
(
2
π
a
)
1
2
exp
(
−
J
2
2
a
)
{\displaystyle \int _{-\infty }^{\infty }\exp \left(-{1 \over 2}ax^{2}+iJx\right)dx=\left({2\pi \over a}\right)^{1 \over 2}\exp \left(-{J^{2} \over 2a}\right)}
es proporcional a la transformada de Fourier de la campana de Gauss, donde
J
{\displaystyle J}
es la variable conjugada de
x
{\displaystyle x}
.
Se puede resolver como el caso anterior completando cuadrados, viendo así que la transformada de Fourier de una gaussiana es también otra gaussiana, pero de la variable conjugada. Cuando mayor es el valor de
a
{\displaystyle a}
, más estrecha es la gaussiana en
x
{\displaystyle x}
y más amplia es la gaussiana en
J
{\displaystyle J}
. Este es un caso particular del principio de incertidumbre .
Integrales con un argumento complejo en el exponente
editar
La integral de interés es (para ver un ejemplo de uso, véase Relación entre la ecuación de Schrödinger y la formulación de la integral de camino de la mecánica cuántica ).
∫
−
∞
∞
exp
(
1
2
i
a
x
2
+
i
J
x
)
d
x
.
{\displaystyle \int _{-\infty }^{\infty }\exp \left({1 \over 2}iax^{2}+iJx\right)dx.}
Ahora asumimos que
a
{\displaystyle a_{}^{}}
y
J
{\displaystyle J_{}^{}}
pueden ser complejos.
Completando el cuadrado
(
1
2
i
a
x
2
+
i
J
x
)
=
1
2
i
a
(
x
2
+
2
J
x
a
+
(
J
a
)
2
−
(
J
a
)
2
)
=
−
1
2
a
i
(
x
+
J
a
)
2
−
i
J
2
2
a
.
{\displaystyle \left({1 \over 2}iax^{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 2a}.}
Por analogía con las integrales previas
∫
−
∞
∞
exp
(
1
2
i
a
x
2
+
i
J
x
)
d
x
=
(
2
π
i
a
)
1
2
exp
(
−
i
J
2
2
a
)
.
{\displaystyle \int _{-\infty }^{\infty }\exp \left({1 \over 2}iax^{2}+iJx\right)dx=\left({2\pi i \over a}\right)^{1 \over 2}\exp \left({-iJ^{2} \over 2a}\right).}
Este resultado es válido como una integración en el plano complejo siempre que
a
{\displaystyle a_{}^{}}
tenga una parte imaginaria positiva (por muy pequeña que esta sea).
Integrales gaussianas en varias dimensiones
editar
La integral gaussiana unidimensional puede ser generalizada a
n
{\displaystyle n_{}^{}}
dimensiones.[ 2]
∫
exp
(
−
1
2
x
⋅
A
⋅
x
+
J
⋅
x
)
d
n
x
=
(
2
π
)
n
det
A
exp
(
1
2
J
⋅
A
−
1
⋅
J
)
{\displaystyle \int \exp \left(-{\frac {1}{2}}x\cdot A\cdot x+J\cdot x\right)d^{n}x={\sqrt {\frac {(2\pi )^{n}}{\det A}}}\exp \left({1 \over 2}J\cdot A^{-1}\cdot J\right)}
Donde
A
{\displaystyle A}
es una matriz simétrica real.
Esta integral se resuelve por medio de una diagonalización de la matriz
A
{\displaystyle A}
a través de una transformación ortogonal
D
=
O
−
1
A
O
=
O
T
A
O
{\displaystyle D_{}^{}=O^{-1}AO=O^{T}AO}
donde
D
{\displaystyle D}
es una matriz diagonal y
O
{\displaystyle O}
es una matriz ortogonal . Esta transformación desacopla las variables de integración, por lo que tenemos
n
{\displaystyle n_{}^{}}
integrales unidimensionales que podemos resolver una a una independientemente.
Esto quedará seguramente más claro con un ejemplo en el caso bidimensional.
Ejemplo: Integración gaussiana en dos dimensiones
editar
La integral gaussiana en 2 dimensiones es
∫
exp
(
−
1
2
A
i
j
x
i
x
j
)
d
2
x
=
(
2
π
)
2
det
A
{\displaystyle \int \exp \left(-{\frac {1}{2}}A_{ij}x^{i}x^{j}\right)d^{2}x={\sqrt {\frac {(2\pi )^{2}}{\det A}}}}
donde
A
{\displaystyle A}
es una matriz simétrica bidimensional definida como
A
=
[
a
c
c
b
]
{\displaystyle A={\bigl [}{\begin{smallmatrix}a&c\\c&b\end{smallmatrix}}{\bigr ]}}
y donde hemos usado el convenio de sumación de Einstein .
Diagonalización de la matriz
editar
El primer paso es diagonalizar la matriz.[ 3] Nótese que
A
i
j
x
i
x
j
≡
x
T
A
x
=
x
T
(
O
O
T
)
A
(
O
O
T
)
x
=
(
x
T
O
)
(
O
T
A
O
)
(
O
T
x
)
{\displaystyle A_{ij}x^{i}x^{j}\equiv x^{T}Ax=x^{T}\left(OO^{T}\right)A\left(OO^{T}\right)x=\left(x^{T}O\right)\left(O^{T}AO\right)\left(O^{T}x\right)}
donde, como
A
{\displaystyle A}
es una matriz simétrica real, podemos elegir una matriz ortogonal (y por tanto también unitaria )
O
{\displaystyle O}
.
Elegimos la matriz
O
{\displaystyle O}
de tal manera que
D
≡
O
T
A
O
{\displaystyle D\equiv O^{T}AO}
sea diagonal.
O
{\displaystyle O}
puede obtenerse a partir de los autovalores de
A
{\displaystyle A}
.
Autovalores de A
editar
Para obtener los autovectores de
A
{\displaystyle A}
, primero tenemos que conocer los autovalores
λ
{\displaystyle \lambda }
de
A
{\displaystyle A}
,
[
a
c
c
b
]
[
u
v
]
=
λ
[
u
v
]
.
{\displaystyle {\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 ]}.}
Los autovalores de
A
{\displaystyle A}
son las soluciones de su polinomio característico
(
a
−
λ
)
(
b
−
λ
)
−
c
2
=
0
{\displaystyle \left(a-\lambda \right)\left(b-\lambda \right)-c^{2}=0}
que son
λ
±
=
1
2
(
a
+
b
)
±
1
2
(
a
−
b
)
2
+
4
c
2
.
{\displaystyle \lambda _{\pm }={1 \over 2}\left(a+b\right)\pm {1 \over 2}{\sqrt {\left(a-b\right)^{2}+4c^{2}}}.}
Autovectores de A
editar
Sustituyendo los autovalores otra vez en la ecuación que define los autovectores, tenemos
v
=
−
(
a
−
λ
±
)
u
c
{\displaystyle v=-{\left(a-\lambda _{\pm }\right)u \over c}}
o
v
=
−
c
u
(
b
−
λ
±
)
.
{\displaystyle v=-{cu \over \left(b-\lambda _{\pm }\right)}.}
De la ecuación característica tenemos que
(
a
−
λ
±
)
c
=
c
(
b
−
λ
±
)
.
{\displaystyle {\left(a-\lambda _{\pm }\right) \over c}={c \over \left(b-\lambda _{\pm }\right)}.}
Nótese también que
(
a
−
λ
±
)
c
=
−
(
b
−
λ
∓
)
c
.
{\displaystyle {\left(a-\lambda _{\pm }\right) \over c}=-{\left(b-\lambda _{\mp }\right) \over c}.}
Los autovectores pueden ser escritos como
[
1
η
−
(
a
−
λ
−
c
η
)
]
{\displaystyle {\begin{bmatrix}{1 \over \eta }\\-\left({a-\lambda _{-} \over c\eta }\right)\end{bmatrix}}}
y
[
−
(
b
−
λ
+
c
η
)
1
η
]
{\displaystyle {\begin{bmatrix}-\left({b-\lambda _{+} \over c\eta }\right)\\{1 \over \eta }\end{bmatrix}}}
Aquí,
η
{\displaystyle \eta }
es una constante de normalización dada por
η
=
1
+
(
a
−
λ
−
c
)
2
=
1
+
(
b
−
λ
+
c
)
2
.
{\displaystyle \eta ={\sqrt {1+\left({a-\lambda _{-} \over c}\right)^{2}}}={\sqrt {1+\left({b-\lambda _{+} \over c}\right)^{2}}}.}
Es fácil verificar que los 2 autovectores son ortogonales entre sí.
Construcción de la matriz ortogonal
editar
La matriz ortogonal
O
{\displaystyle O}
se construye usando los autovectores normalizados de la matriz
A
{\displaystyle A}
como columnas de la matriz ortogonal.
O
=
[
1
η
−
(
b
−
λ
+
c
η
)
−
(
a
−
λ
−
c
η
)
1
η
]
.
{\displaystyle 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}}.}
Nótese que el determinante de
O
{\displaystyle O}
, al ser una matriz unitaria, es igual a uno.
Si definimos
sin
(
θ
)
≡
−
(
a
−
λ
−
c
η
)
{\displaystyle \sin \left(\theta \right)\equiv -\left({a-\lambda _{-} \over c\eta }\right)}
entonces, la matriz ortogonal puede escribirse como
O
=
[
cos
(
θ
)
−
sin
(
θ
)
sin
(
θ
)
cos
(
θ
)
]
{\displaystyle O={\begin{bmatrix}\cos \left(\theta \right)&-\sin \left(\theta \right)\\\sin \left(\theta \right)&\cos \left(\theta \right)\end{bmatrix}}}
lo que representa una rotación de los autovectores.
La inversa de esta matriz es
O
−
1
=
O
T
=
[
cos
(
θ
)
sin
(
θ
)
−
sin
(
θ
)
cos
(
θ
)
]
.
{\displaystyle O^{-1}=O^{T}={\begin{bmatrix}\cos \left(\theta \right)&\sin \left(\theta \right)\\-\sin \left(\theta \right)&\cos \left(\theta \right)\end{bmatrix}}.}
Matriz diagonal
editar
La matriz diagonal es
D
=
O
T
A
O
=
[
λ
−
0
0
λ
+
]
{\displaystyle D=O^{T}AO={\begin{bmatrix}\lambda _{-}&0\\0&\lambda _{+}\end{bmatrix}}}
con autovectores
[
1
0
]
{\displaystyle {\begin{bmatrix}1\\0\end{bmatrix}}}
y
[
0
1
]
{\displaystyle {\begin{bmatrix}0\\1\end{bmatrix}}}
Ejemplo numérico
editar
A
=
[
2
1
1
1
]
{\displaystyle A={\begin{bmatrix}2&1\\1&1\end{bmatrix}}}
Los autovalores son
λ
±
=
3
2
±
5
2
.
{\displaystyle \lambda _{\pm }={3 \over 2}\pm {{\sqrt {5}} \over 2}.}
Los autovectores son
1
η
[
1
−
1
2
−
5
2
]
{\displaystyle {1 \over \eta }{\begin{bmatrix}1\\{-{1 \over 2}-{{\sqrt {5}} \over 2}}\end{bmatrix}}}
y
1
η
[
1
2
+
5
2
1
]
{\displaystyle {1 \over \eta }{\begin{bmatrix}{{1 \over 2}+{{\sqrt {5}} \over 2}}\\1\end{bmatrix}}}
donde
η
=
5
2
+
5
2
{\displaystyle \eta ={\sqrt {{5 \over 2}+{{\sqrt {5}} \over 2}}}}
.
La matriz ortogonal es
O
=
[
1
η
1
η
(
1
2
+
5
2
)
1
η
(
−
1
2
−
5
2
)
1
η
]
.
{\displaystyle 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}}.}
Es fácil comprobar que el determinante de
O
{\displaystyle O}
es 1.
La matriz inversa de
O
{\displaystyle O}
es
O
=
[
1
η
1
η
(
−
1
2
−
5
2
)
1
η
(
1
2
+
5
2
)
1
η
]
.
{\displaystyle 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}}.}
Por lo que la matriz diagonal
D
{\displaystyle D}
es
D
=
O
T
A
O
=
[
λ
−
0
0
λ
+
]
=
[
(
3
2
−
5
2
)
0
0
(
3
2
+
5
2
)
]
{\displaystyle D=O^{T}AO={\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}}}
con autovectores
[
1
0
]
{\displaystyle {\begin{bmatrix}1\\0\end{bmatrix}}}
y
[
0
1
]
.
{\displaystyle {\begin{bmatrix}0\\1\end{bmatrix}}.}
Reescalado de las variables e integración
editar
Tras la diagonalización, la integral a resolver se puede escribir como
∫
exp
(
−
1
2
x
T
A
x
)
d
2
x
=
∫
exp
(
−
1
2
∑
j
=
1
2
λ
j
y
j
2
)
d
2
y
,
{\displaystyle \int \exp \left(-{\frac {1}{2}}x^{T}Ax\right)d^{2}x=\int \exp \left(-{\frac {1}{2}}\sum _{j=1}^{2}\lambda _{j}y_{j}^{2}\right)\,d^{2}y,}
donde
y
=
O
T
x
{\displaystyle y_{}^{}=O^{T}x}
.
Como la transformación de coordenadas aplicada es una rotación, el determinante jacobiano de la transformación es igual a la unidad, por lo que
d
y
2
=
d
x
2
{\displaystyle dy_{}^{2}=dx^{2}}
.
Con lo que ya podemos resolver la integral
∫
exp
(
−
1
2
x
T
A
x
)
d
2
x
=
∫
exp
(
−
1
2
∑
j
=
1
2
λ
j
y
j
2
)
d
2
y
=
∏
j
=
1
2
(
2
π
λ
j
)
1
2
=
(
(
2
π
)
2
∏
j
=
1
2
λ
j
)
1
2
=
(
(
2
π
)
2
det
(
O
−
1
A
O
)
)
1
2
=
(
(
2
π
)
2
det
(
A
)
)
1
2
,
{\displaystyle \int \exp \left(-{\frac {1}{2}}x^{T}Ax\right)d^{2}x=\int \exp \left(-{\frac {1}{2}}\sum _{j=1}^{2}\lambda _{j}y_{j}^{2}\right)d^{2}y=\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},}
que es la solución buscada.
Integrales con términos complejos lineales en varias dimensiones
editar
Una vez que tenemos resuelto el ejemplo bidimensional, es fácil generalizar el resultado obtenido al plano complejo y a varias dimensiones.
Integrales con un término lineal real en el argumento
editar
∫
exp
(
−
1
2
x
⋅
A
⋅
x
+
J
⋅
x
)
d
n
x
=
(
2
π
)
n
det
A
exp
(
1
2
J
⋅
A
−
1
⋅
J
)
{\displaystyle \int \exp \left(-{\frac {1}{2}}x\cdot A\cdot x+J\cdot x\right)d^{n}x={\sqrt {\frac {(2\pi )^{n}}{\det A}}}\exp \left({1 \over 2}J\cdot A^{-1}\cdot J\right)}
Integrales con un término lineal imaginario en el argumento
editar
∫
exp
(
−
1
2
x
⋅
A
⋅
x
+
i
J
⋅
x
)
d
n
x
=
(
2
π
)
n
det
A
exp
(
−
1
2
J
⋅
A
−
1
⋅
J
)
{\displaystyle \int \exp \left(-{\frac {1}{2}}x\cdot A\cdot x+iJ\cdot x\right)d^{n}x={\sqrt {\frac {(2\pi )^{n}}{\det A}}}\exp \left(-{1 \over 2}J\cdot A^{-1}\cdot J\right)}
Integrales con un término complejo cuadrático
editar
∫
exp
(
i
2
x
⋅
A
⋅
x
+
i
J
⋅
x
)
d
n
x
=
(
2
π
i
)
n
det
A
exp
(
−
i
2
J
⋅
A
−
1
⋅
J
)
{\displaystyle \int \exp \left({\frac {i}{2}}x\cdot A\cdot x+iJ\cdot x\right)d^{n}x={\sqrt {\frac {(2\pi i)^{n}}{\det A}}}\exp \left(-{i \over 2}J\cdot A^{-1}\cdot J\right)}
Integrales con operadores diferenciales en el argumento
editar
Por ejemplo, véase la integral[ 4]
∫
exp
[
∫
d
4
x
(
−
1
2
φ
A
^
φ
+
J
φ
)
]
D
φ
{\displaystyle \int \exp \left[\int d^{4}x\left(-{\frac {1}{2}}\varphi {\hat {A}}\varphi +J\varphi \right)\right]D\varphi }
donde
A
^
{\displaystyle {\hat {A}}}
es un operador diferencial con
φ
{\displaystyle \varphi }
y
J
{\displaystyle J}
funciones definidas en el espacio-tiempo .
D
φ
{\displaystyle D\varphi }
representa una integración sobre todos los posibles caminos . Análogamente a la versión matricial de esta integral, la solución es
∫
exp
(
−
1
2
φ
A
^
φ
+
J
φ
)
D
φ
∝
exp
(
1
2
∫
d
4
x
d
4
y
J
(
x
)
D
(
x
−
y
)
J
(
y
)
)
{\displaystyle \int \exp \left(-{\frac {1}{2}}\varphi {\hat {A}}\varphi +J\varphi \right)D\varphi \;\propto \;\exp \left({1 \over 2}\int d^{4}x\;d^{4}y\,J\left(x\right)D\left(x-y\right)J\left(y\right)\right)}
donde
A
^
D
(
x
−
y
)
=
δ
4
(
x
−
y
)
,
{\displaystyle {\hat {A}}D\left(x-y\right)=\delta ^{4}\left(x-y\right),}
D
(
x
−
y
)
{\displaystyle D\left(x-y\right)}
, es el propagador , que es la inversa de
A
^
{\displaystyle {\hat {A}}}
, y
δ
4
(
x
−
y
)
{\displaystyle \delta ^{4}\left(x-y\right)}
es la delta de Dirac .
Análogamente se obtienen los siguientes resultados
∫
exp
[
∫
d
4
x
(
−
1
2
φ
A
^
φ
+
i
J
φ
)
]
D
φ
∝
exp
(
−
1
2
∫
d
4
x
d
4
y
J
(
x
)
D
(
x
−
y
)
J
(
y
)
)
{\displaystyle \int \exp \left[\int d^{4}x\left(-{\frac {1}{2}}\varphi {\hat {A}}\varphi +iJ\varphi \right)\right]D\varphi \;\propto \;\exp \left(-{1 \over 2}\int d^{4}x\;d^{4}y\,J\left(x\right)D\left(x-y\right)J\left(y\right)\right)}
∫
exp
[
i
∫
d
4
x
(
1
2
φ
A
^
φ
+
J
φ
)
]
D
φ
∝
exp
(
−
i
2
∫
d
4
x
d
4
y
J
(
x
)
D
(
x
−
y
)
J
(
y
)
)
{\displaystyle \int \exp \left[i\int d^{4}x\left({\frac {1}{2}}\varphi {\hat {A}}\varphi +J\varphi \right)\right]D\varphi \;\propto \;\exp \left(-{i \over 2}\int d^{4}x\;d^{4}y\,J\left(x\right)D\left(x-y\right)J\left(y\right)\right)}
.
Integrales que pueden ser aproximados por el método de descenso más pronunciado
editar
Integrales que pueden ser aproximadas por el método de la fase estacionaria
editar
Otra integral común es la integral de caminos de la forma
∫
exp
(
i
ℏ
S
(
q
,
q
˙
)
)
D
q
{\displaystyle \int \exp \left({i \over \hbar }S\left(q,{\dot {q}}\right)\right)Dq}
donde
S
(
q
,
q
˙
)
{\displaystyle S\left(q,{\dot {q}}\right)}
es la acción clásica y el dominio de integración son todos los posibles caminos que una partícula podría tomar. En el límite de
ℏ
{\displaystyle \hbar }
pequeña, la integral puede aproximarse por el método de la fase estacionaria . En esta aproximación, la integral queda definida solamente sobre el camino en el cual la acción es un punto estacionario (un mínimo, un máximo o un punto de silla ). Por tanto, a través de esta aproximación se recupera la mecánica clásica como el límite clásico de la mecánica cuántica.
Integrales de Fourier
editar
Función delta de Dirac
editar
La función delta de Dirac puede representarse de la siguiente manera como la siguiente transformada de Fourier [ 5]
∫
d
4
k
(
2
π
)
4
exp
(
i
k
(
x
−
y
)
)
=
δ
4
(
x
−
y
)
{\displaystyle \int {d^{4}k \over \left(2\pi \right)^{4}}\;{\exp \left(ik\left(x-y\right)\right)=\delta ^{4}\left(x-y\right)}}
.
En el caso
N
{\displaystyle N}
dimensional tenemos la generalización inmediata
∫
d
N
k
(
2
π
)
N
exp
(
i
k
(
x
−
y
)
)
=
δ
N
(
x
−
y
)
{\displaystyle \int {d^{N}k \over \left(2\pi \right)^{N}}\;{\exp \left(ik\left(x-y\right)\right)=\delta ^{N}\left(x-y\right)}}
.
Integrales de Fourier de generalizaciones del potencial de Coulomb
editar
Laplaciano de 1/r
editar
Aplicando el teorema de Gauss se tiene que, en el espacio tridimensional
−
1
4
π
∇
2
(
1
r
)
=
δ
(
r
)
{\displaystyle -{1 \over 4\pi }\nabla ^{2}\left({1 \over r}\right)=\delta \left(\mathbf {r} \right)}
donde
r
2
=
r
⋅
r
{\displaystyle r^{2}=\mathbf {r} \cdot \mathbf {r} }
Esta identidad se puede usar para derivar igualdades entre integrales. Por ejemplo, véase la descomposición en campos longitudinal y transversal .
A partir de la igualdad de arriba, se obtiene inmediatamente que la representación de 1/r como integral de Fourier es
∫
d
3
k
(
2
π
)
3
exp
(
i
k
⋅
r
)
k
2
=
1
4
π
r
.
{\displaystyle \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}.}
Potencial de Yukawa: El potencial de Coulomb con masa
editar
El potencial de Yukawa en 3 dimensiones puede representarse como una integral de su transformada de Fourier [ 6]
∫
d
3
k
(
2
π
)
3
exp
(
i
k
⋅
r
)
k
2
+
m
2
=
e
−
m
r
4
π
r
{\displaystyle \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^{-mr} \over 4\pi r}}
donde
r
2
=
r
⋅
r
{\displaystyle r^{2}=\mathbf {r} \cdot \mathbf {r} }
y
k
2
=
k
⋅
k
.
{\displaystyle k^{2}=\mathbf {k} \cdot \mathbf {k} .}
En el límite de masa
m
{\displaystyle m}
pequeña, esta integral se reduce a
1
4
π
r
{\displaystyle {1 \over 4\pi r}}
.
Para obtener este límite de la forma correcta, nótese que:
∫
d
3
k
(
2
π
)
3
exp
(
i
k
⋅
r
)
k
2
+
m
2
=
∫
0
∞
k
2
d
k
(
2
π
)
2
∫
−
1
1
d
u
exp
(
i
k
r
u
)
k
2
+
m
2
{\displaystyle \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
r
∫
0
∞
k
d
k
(
2
π
)
2
sin
(
k
r
)
k
2
+
m
2
=
1
i
r
∫
−
∞
∞
k
d
k
(
2
π
)
2
exp
(
i
k
r
)
k
2
+
m
2
{\displaystyle ={2 \over r}\int _{0}^{\infty }{kdk \over \left(2\pi \right)^{2}}{\sin \left(kr\right) \over k^{2}+m^{2}}={1 \over ir}\int _{-\infty }^{\infty }{kdk \over \left(2\pi \right)^{2}}{\exp \left(ikr\right) \over k^{2}+m^{2}}}
=
1
i
r
∫
−
∞
∞
k
d
k
(
2
π
)
2
exp
(
i
k
r
)
(
k
+
i
m
)
(
k
−
i
m
)
=
1
i
r
2
π
i
(
2
π
)
2
i
m
2
i
m
exp
(
−
m
r
)
{\displaystyle ={1 \over ir}\int _{-\infty }^{\infty }{kdk \over \left(2\pi \right)^{2}}{\exp \left(ikr\right) \over \left(k+im\right)\left(k-im\right)}={1 \over ir}{2\pi i \over \left(2\pi \right)^{2}}{im \over 2im}\exp \left(-mr\right)}
Potencial de Coulomb modificado con masa
editar
∫
d
3
k
(
2
π
)
3
(
k
^
⋅
r
^
)
2
exp
(
i
k
⋅
r
)
k
2
+
m
2
=
e
−
m
r
4
π
r
{
1
+
2
m
r
−
2
(
m
r
)
2
(
e
m
r
−
1
)
}
{\displaystyle \int {d^{3}k \over \left(2\pi \right)^{3}}\left(\mathbf {\hat {k}} \cdot \mathbf {\hat {r}} \right)^{2}{\exp \left(i\mathbf {k} \cdot \mathbf {r} \right) \over k^{2}+m^{2}}={e^{-mr} \over 4\pi r}\left\{1+{2 \over mr}-{2 \over \left(mr\right)^{2}}\left(e^{mr}-1\right)\right\}}
donde el gorro indica que el vector es de módulo unidad en el espacio tridimensional.
En el límite de masa
m
{\displaystyle m}
pequeña, esta integral se anula.
Para obtener este límite de la forma correcta, nótese que:
∫
d
3
k
(
2
π
)
3
(
k
^
⋅
r
^
)
2
exp
(
i
k
⋅
r
)
k
2
+
m
2
=
∫
0
∞
k
2
d
k
(
2
π
)
2
∫
−
1
1
d
u
u
2
exp
(
i
k
r
u
)
k
2
+
m
2
{\displaystyle \int {d^{3}k \over \left(2\pi \right)^{3}}\left(\mathbf {\hat {k}} \cdot \mathbf {\hat {r}} \right)^{2}{\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}duu^{2}{\exp \left(ikru\right) \over k^{2}+m^{2}}}
=
2
∫
0
∞
k
2
d
k
(
2
π
)
2
1
k
2
+
m
2
{
1
k
r
sin
(
k
r
)
+
2
1
(
k
r
)
2
cos
(
k
r
)
−
2
1
(
k
r
)
3
sin
(
k
r
)
}
{\displaystyle ={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
(
−
m
r
)
4
π
r
{
g
(
m
r
)
}
{\displaystyle ={\exp \left(-mr\right) \over 4\pi r}\left\{g\left(mr\right)\right\}}
donde
g
(
m
r
)
=
1
+
2
m
r
−
2
(
m
r
)
2
(
e
m
r
−
1
)
.
{\displaystyle g\left(mr\right)=1+{2 \over mr}-{2 \over \left(mr\right)^{2}}\left(e^{mr}-1\right).}
Nótese que, en este límite, cuando
m
→
0
,
{\displaystyle m\rightarrow 0,}
g
(
m
r
)
→
−
m
r
3
,
{\displaystyle g\left(mr\right)\rightarrow -{\frac {mr}{3}},}
por lo que el lado derecho de la igualdad es proporcional a
m
{\displaystyle m}
, que en este límite es cero.
Potencial longitudinal con masa
editar
∫
d
3
k
(
2
π
)
3
k
^
k
^
exp
(
i
k
⋅
r
)
k
2
+
m
2
=
1
2
e
−
m
r
4
π
r
[
1
−
r
^
r
^
]
+
1
2
e
−
m
r
4
π
r
{
1
+
2
m
r
−
2
(
m
r
)
2
(
e
m
r
−
1
)
}
[
1
+
r
^
r
^
]
{\displaystyle \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}}={1 \over 2}{e^{-mr} \over 4\pi r}\left[\mathbf {1} -\mathbf {\hat {r}} \mathbf {\hat {r}} \right]+{1 \over 2}{e^{-mr} \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]}
donde el gorro indica que el vector es de módulo unidad en el espacio tridimensional.
En el límite de masa
m
{\displaystyle m}
pequeña, esta integral se aproxima a
1
2
1
4
π
r
[
1
−
r
^
r
^
]
.
{\displaystyle {1 \over 2}{1 \over 4\pi r}\left[\mathbf {1} -\mathbf {\hat {r}} \mathbf {\hat {r}} \right].}
La integral se resuelve en coordenadas esféricas:
∫
d
3
k
(
2
π
)
3
k
^
k
^
exp
(
i
k
⋅
r
)
k
2
+
m
2
=
∫
d
3
k
(
2
π
)
3
[
(
k
^
⋅
r
^
)
2
r
^
r
^
+
(
k
^
⋅
θ
^
)
2
θ
^
θ
^
+
(
k
^
⋅
ϕ
^
)
2
ϕ
^
ϕ
^
]
exp
(
i
k
⋅
r
)
k
2
+
m
2
{\displaystyle \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}}}
Usando los resultados anteriores, se obtiene
=
1
4
π
r
exp
(
−
m
r
)
{
1
+
2
m
r
−
2
(
m
r
)
2
(
e
m
r
−
1
)
}
{
1
−
1
2
[
1
−
r
^
r
^
]
}
+
∫
0
∞
k
2
d
k
(
2
π
)
2
∫
−
1
1
d
u
exp
(
i
k
r
u
)
k
2
+
m
2
1
2
[
1
−
r
^
r
^
]
{\displaystyle ={1 \over 4\pi r}\exp \left(-mr\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
2
e
−
m
r
4
π
r
[
1
−
r
^
r
^
]
+
e
−
m
r
4
π
r
{
1
+
2
m
r
−
2
(
m
r
)
2
(
e
m
r
−
1
)
}
{
1
2
[
1
+
r
^
r
^
]
}
{\displaystyle ={1 \over 2}{e^{-mr} \over 4\pi r}\left[\mathbf {1} -\mathbf {\hat {r}} \mathbf {\hat {r}} \right]+{e^{-mr} \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\}}
.
Potencial transversal con masa
editar
∫
d
3
k
(
2
π
)
3
[
1
−
k
^
k
^
]
exp
(
i
k
⋅
r
)
k
2
+
m
2
=
1
2
e
−
m
r
4
π
r
{
2
(
m
r
)
2
(
e
m
r
−
1
)
−
2
m
r
}
[
1
+
r
^
r
^
]
{\displaystyle \int {d^{3}k \over \left(2\pi \right)^{3}}\;\left[\mathbf {1} -\mathbf {\hat {k}} \mathbf {\hat {k}} \right]\;{\exp \left(i\mathbf {k} \cdot \mathbf {r} \right) \over k^{2}+m^{2}}={1 \over 2}{e^{-mr} \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]}
En el límite de
m
r
{\displaystyle mr}
pequeño, esta integral se aproxima a
1
2
1
4
π
r
[
1
+
r
^
r
^
]
.
{\displaystyle {1 \over 2}{1 \over 4\pi r}\left[\mathbf {1} +\mathbf {\hat {r}} \mathbf {\hat {r}} \right].}
Para grandes distancias, la integral decae cúbicamente con
r
{\displaystyle r}
1
4
π
m
2
r
3
[
1
+
r
^
r
^
]
.
{\displaystyle {1 \over 4\pi m^{2}r^{3}}\left[\mathbf {1} +\mathbf {\hat {r}} \mathbf {\hat {r}} \right].}
Integración angular en coordenadas cilíndricas
editar
La integración angular de una exponencial en coordenadas cilíndricas se puede escribir en función de funciones de Bessel de primer tipo[ 7] [ 8]
∫
0
2
π
d
φ
2
π
exp
(
i
p
cos
(
φ
)
)
=
J
0
(
p
)
,
{\displaystyle \int _{0}^{2\pi }{d\varphi \over 2\pi }\exp \left(ip\cos \left(\varphi \right)\right)=J_{0}\left(p\right),}
∫
0
2
π
d
φ
2
π
cos
(
φ
)
exp
(
i
p
cos
(
φ
)
)
=
i
J
1
(
p
)
.
{\displaystyle \int _{0}^{2\pi }{d\varphi \over 2\pi }\cos \left(\varphi \right)\exp \left(ip\cos \left(\varphi \right)\right)=iJ_{1}\left(p\right).}
Funciones de Bessel
editar
Integración del propagador del campo con masa en coordenadas cilíndricas
editar
Función de Bessel lineal
editar
∫
0
∞
k
d
k
k
2
+
m
2
J
0
(
k
r
)
=
K
0
(
m
r
)
.
{\displaystyle \int _{0}^{\infty }{k\;dk \over k^{2}+m^{2}}J_{0}\left(kr\right)=K_{0}\left(mr\right).}
Véase Abramowitz and Stegun.[ 9]
Para
m
r
<<
1
{\displaystyle mr<<1}
, tenemos[ 10]
K
0
(
m
r
)
→
−
ln
(
m
r
2
)
+
γ
,
{\displaystyle K_{0}\left(mr\right)\rightarrow -\ln \left({mr \over 2}\right)+\gamma ,}
donde
γ
≈
0.5772
{\displaystyle \gamma \approx 0.5772}
es la constante de Euler-Mascheroni .
Cuadrado de funciones de Bessel
editar
∫
0
∞
k
d
k
k
2
+
m
2
J
1
2
(
k
r
)
=
I
1
(
m
r
)
K
1
(
m
r
)
.
{\displaystyle \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).}
Véase[ 11]
Para
m
r
<<
1
{\displaystyle mr<<1}
la integral se reduce a
∫
o
∞
k
d
k
k
2
+
m
2
J
1
2
(
k
r
)
→
1
2
[
1
−
1
8
(
m
r
)
2
]
,
{\displaystyle \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],}
mientras que para
m
r
>>
1
{\displaystyle mr>>1}
se tiene
∫
o
∞
k
d
k
k
2
+
m
2
J
1
2
(
k
r
)
→
1
2
(
1
m
r
)
.
{\displaystyle \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).}
En general
∫
0
∞
k
d
k
k
2
+
m
2
J
ν
2
(
k
r
)
=
I
ν
(
m
r
)
K
ν
(
m
r
)
ℜ
ν
>
−
1.
{\displaystyle \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.}
Integración de una función de onda magnética
editar
La integral bidimensional sobre una función de onda magnética es[ 12]
2
a
2
n
+
2
n
!
∫
0
∞
d
r
r
2
n
+
1
exp
(
−
a
2
r
2
)
J
0
(
k
r
)
=
M
(
n
+
1
,
1
,
−
k
2
4
a
2
)
,
{\displaystyle {2a^{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),}
donde
M
{\displaystyle M}
es una función hipergeométrica confluente .
Véase también
editar
↑ A. Zee (2003). Quantum Field Theory in a Nutshell . Princeton University. ISBN 0-691-01019-6 . pp. 13-15
↑ Frederick W. Byron and Robert W. Fuller (1969). Mathematics of Classical and Quantum Physics . Addison-Wesley. ISBN 0-201-00746-2 .
↑ Herbert S. Wilf (1978). Mathematics for the Physical Sciences . Dover. ISBN 0-486-63635-6 .
↑ Zee, pp. 21-22.
↑ Zee, p. 23.
↑ Zee, p. 26, 29.
↑ I. S. Gradshteyn and I. M. Ryzhik (1965). Tables of Integrals, Series, and Products . Academic Press. ISBN 65-29097 . p. 402
↑ Jackson, John D. (1998). Classical Electrodynamics (3rd ed.) . Wiley. ISBN 047130932X . p. 113
↑ M. Abramowitz and I. Stegun (1965). Handbook of Mathematical Functions . Dover. ISBN 486-61272-4 . Section 11.4.44
↑ Jackson, p. 116
↑ I. S. Gradshteyn and I. M. Ryzhik, p. 679
↑ Abramowitz and Stegun, Section 11.4.28