En física, la aproximación WKB es un método para encontrar soluciones aproximadas a ecuaciones diferenciales lineales con coeficientes variables. Se usa especialmente para cálculos semiclásicos en mecánica cuántica en los que la función de onda se escribe como una exponencial cuya amplitud o fase
varían lentamente.
El nombre de este método es un acrónimo de aproximación Wentzel-Kramers-Brillouin. Otros acrónimos usualmente usados son aproximación JWKB y aproximación WKBJ, donde la "J" representa a Jeffreys.
Este método lleva el nombre de los físicos Wentzel, Kramers, y Brillouin, quienes lo desarrollaron en 1926. En 1923, el matemático Harold Jeffreys desarrolló un método general de aproximación a soluciones lineales de ecuaciones diferenciales de segundo orden que incluía la ecuación de Schrödinger. A pesar de que la ecuación de Schrödinger fue propuesta dos años después, Wentzel, Kramers y Brillouin parece que no estaban al tanto del trabajo previo de Jeffreys por lo que a veces se excluye a éste del reconocimiento. Los primeros textos de mecánica cuántica contienen combinaciones de sus iniciales, que incluyen WBK, BWK, WKBJ, JWKB y BWKJ.
Referencias anteriores al método son:
Carlini en 1817, Liouville en 1837, Green en 1837, Rayleigh en 1912 y Gans en 1915. Liouville y Green pueden ser llamados los fundadores del método, en 1837, y esto es también comúnmente llamado como "Liouville-Green" o "método LG".
La importante contribución de Jeffreys, Wentzel, Kramers y Brillouin al método fue la inclusión del tratamiento momentos, conectando la Fugacidad y oscillatorio soluciones en ningún lado del momento. Por ejemplo, esto puede ocurrir en la ecuación de Schrödinger, ambos para un pico potencial de energía.
En general, la teoría WKB es un método para aproximar la solución de una ecuación diferencial cuya más alta derivada va multiplicada por un pequeño parámetro ε. El método de aproximación es como sigue. Dada una ecuación diferencial
![{\displaystyle \epsilon {\frac {\mathrm {d} ^{n}y}{\mathrm {d} x^{n}}}+a(x){\frac {\mathrm {d} ^{n-1}y}{\mathrm {d} x^{n-1}}}+\cdots +k(x){\frac {\mathrm {d} y}{\mathrm {d} x}}+m(x)y=0}](//wikimedia.org/api/rest_v1/media/math/render/svg/9f754fbf639145c1b4938cdeb23e6aa791f860a3)
Asume una solución de la forma de una expansión de serie asintótica
![{\displaystyle y(x)\sim \exp \left[{\frac {1}{\delta }}\sum _{i=0}^{\infty }\delta ^{i}S_{i}(x)\right]}](//wikimedia.org/api/rest_v1/media/math/render/svg/1f9e57404679fd2dfc266416dad956732681478b)
En el límite
, la sustitución del ansatz anterior dentro de la ecuación diferencial y la cancelación de los términos exponenciales permite obtener la solución de un número arbitrario de términos
de la expansión. La teoría WKB es un caso especial de Análisis de escala múltiple.
Considere la ecuación diferencial lineal homogénea de segundo orden.
![{\displaystyle \epsilon ^{2}{\frac {d^{2}y}{dx^{2}}}=Q(x)y}](//wikimedia.org/api/rest_v1/media/math/render/svg/db8e603b55b7f04deefe0a0542b38da440135abd)
donde
. Reemplezando con
![{\displaystyle y(x)=\exp \left[{\frac {1}{\delta }}\sum _{n=0}^{\infty }\delta ^{n}S_{n}(x)\right]}](//wikimedia.org/api/rest_v1/media/math/render/svg/8d7ab03767fae0bcbd9465105922ea1878725028)
resulta en la ecuación
![{\displaystyle \epsilon ^{2}\left[{\frac {1}{\delta ^{2}}}\left(\sum _{n=0}^{\infty }\delta ^{n}S_{n}'\right)^{2}+{\frac {1}{\delta }}\sum _{n=0}^{\infty }\delta ^{n}S_{n}''\right]=Q(x)}](//wikimedia.org/api/rest_v1/media/math/render/svg/cb495cfee3874faedb93e401165effac6585ea5a)
De primer orden, (suponiendo, por el momento, la serie será asintóticamente consistente) sobre él se puede aproximar como:
![{\displaystyle {\frac {\epsilon ^{2}}{\delta ^{2}}}S_{0}'^{2}+{\frac {2\epsilon ^{2}}{\delta }}S_{0}'S_{1}'+{\frac {\epsilon ^{2}}{\delta }}S_{0}''=Q(x)}](//wikimedia.org/api/rest_v1/media/math/render/svg/d3934d29131f7aeb70a78faea46b2e7bd7b2e5a8)
En el límite
, el equilibrio dominante está dado por:
![{\displaystyle {\frac {\epsilon ^{2}}{\delta ^{2}}}S_{0}'^{2}\sim Q(x)}](//wikimedia.org/api/rest_v1/media/math/render/svg/4e695cba9c3fceb3266c38d58cdc9876e733549f)
Por lo tanto δ es proporcional a ε. Valorizando a ellos igual y comparando podemos hacer:
![{\displaystyle \epsilon ^{0}:\;\;\;S_{0}'^{2}=Q(x)}](//wikimedia.org/api/rest_v1/media/math/render/svg/475ec63145f9a6d87e788f9cd9143877ceef8140)
La cual puede ser reconocida como la ecuación Eikonal, con solución
![{\displaystyle S_{0}(x)=\pm \int _{x_{0}}^{x}{\sqrt {Q(t)}}\,dt}](//wikimedia.org/api/rest_v1/media/math/render/svg/55b5467e7114e18f78238ddc41a8b8215aee0714)
En cuanto a las potencias de primer orden de
da
![{\displaystyle \epsilon ^{1}:\;\;\;2S_{0}'S_{1}'+S_{0}''=0}](//wikimedia.org/api/rest_v1/media/math/render/svg/b2f3f137cf861aded9d98000f0c6eb2042c99a4c)
La cual es ecuación de transporte unidimensional, el cual tiene la solución
![{\displaystyle S_{1}(x)=-{\frac {1}{4}}\log \left(Q(x)\right)+k_{1}.\,}](//wikimedia.org/api/rest_v1/media/math/render/svg/ed7479ac44dc0c0006ef147ca821bbaa998e5ea8)
Y
es una constante arbitraria. Nosotros ahora tenemos un par de aproximaciones para el sistema(un par porque
puede tomar dos signos); la aproximación-WKB de primer orden será una combinación lineal de las dos:
![{\displaystyle y(x)\approx c_{1}Q^{-{\frac {1}{4}}}(x)\exp \left[{\frac {1}{\epsilon }}\int _{x_{0}}^{x}{\sqrt {Q(t)}}dt\right]+c_{2}Q^{-{\frac {1}{4}}}(x)\exp \left[-{\frac {1}{\epsilon }}\int _{x_{0}}^{x}{\sqrt {Q(t)}}dt\right]}](//wikimedia.org/api/rest_v1/media/math/render/svg/ecf8d08bab3d084d20ac38886939df9beb809781)
Términos de mayor orden pueden ser obtenidos por examinar en ecuaciones para más altas potencia de ε. Explícitamente
![{\displaystyle 2S_{0}'S_{n}'+S''_{n-1}+\sum _{j=1}^{n-1}S'_{j}S'_{n-j}=0}](//wikimedia.org/api/rest_v1/media/math/render/svg/3e354b3114e4c6f7bfc2215768a88bf1e31befe5)
para
. Este ejemplo viene de los libros de texto Bender y Orszag (ver referencias).
La ecuación de Schrödinger independiente del tiempo en una dimensión viene dada por:
,
La cual puede reescribirse como:
.
La función de onda puede reescribirse como la exponencial de otra función Φ (La cual esta estrechamente relacionada con la acción):
![{\displaystyle \Psi (x)=e^{\Phi (x)},\!}](//wikimedia.org/api/rest_v1/media/math/render/svg/857736b04359299f7326b0639ae5935d8b3ad7fc)
Así que:
![{\displaystyle \Phi ''(x)+\left[\Phi '(x)\right]^{2}={\frac {2m}{\hbar ^{2}}}\left(V(x)-E\right),}](//wikimedia.org/api/rest_v1/media/math/render/svg/8890bfffc1739c238e0acc62f3755e5224dd61f6)
Donde Φ' indica la derivada de Φ con respecto a x. La derivada
puede separarse en parte real e imaginaria introducciendo las funciones reales A y B:
![{\displaystyle \Phi '(x)=A(x)+iB(x).\;}](//wikimedia.org/api/rest_v1/media/math/render/svg/4d3eecc389cf2d77bab807d924ad29a790e7dea0)
La amplitud de la función de onda es entonces
, mientras que la fase es
. Las partes real e imaginaria de la ecuación de Schrödinger entonces son:
![{\displaystyle A'(x)+A(x)^{2}-B(x)^{2}={\frac {2m}{\hbar ^{2}}}\left(V(x)-E\right),}](//wikimedia.org/api/rest_v1/media/math/render/svg/4cdbdd9a7939455a7c9b1b54dffeb78a9479cc06)
![{\displaystyle B'(x)+2A(x)B(x)=0.\;}](//wikimedia.org/api/rest_v1/media/math/render/svg/57aefbd5b5468e207dddd63819bea397eaaa9a11)
Tras esto, y usando la aproximación semiclásica, podemos escribir cada función como una serie de potencias en
. Desde la ecuación, esto puede ser visto como que la serie de potencias puede comenzar con al menos un orden de
para satisfacer la parte real de la ecuación. Con el fin de alcanzar un límite clásico bueno, es necesario comenzar con tan alta potencia de la constante de Planck como sea posible.
![{\displaystyle A(x)={\frac {1}{\hbar }}\sum _{n=0}^{\infty }\hbar ^{n}A_{n}(x)}](//wikimedia.org/api/rest_v1/media/math/render/svg/45d5e18baa97dc59b984ff6fce461c768f9f07b3)
![{\displaystyle B(x)={\frac {1}{\hbar }}\sum _{n=0}^{\infty }\hbar ^{n}B_{n}(x)}](//wikimedia.org/api/rest_v1/media/math/render/svg/343a26352f3c6edd0e67dbe769a5cc7cd682a97f)
De primer orden en esta expansión, las condiciones sobre A y B pueden ser escritas.
![{\displaystyle A_{0}(x)^{2}-B_{0}(x)^{2}=2m\left(V(x)-E\right)}](//wikimedia.org/api/rest_v1/media/math/render/svg/8d579d4c3e26ebe95d9e8eba484bb2e9b7c828fe)
![{\displaystyle A_{0}(x)B_{0}(x)=0\;}](//wikimedia.org/api/rest_v1/media/math/render/svg/6a2cbc5d764554198754fc0ca74b76b4c4b18ae3)
Si la amplitud varía con la suficiente lentitud en comparación con la fase (
), se deduce que:
![{\displaystyle B_{0}(x)=\pm {\sqrt {2m\left(E-V(x)\right)}},}](//wikimedia.org/api/rest_v1/media/math/render/svg/fb768050074c2836fb55be5248bb2c8fbe75d0ec)
que sólo es válido cuando la energía total es mayor que la energía potencial, como es siempre el caso en movimiento clásico. Después que el mismo procedimiento sobre el siguiente orden de la expansión se deduce que:
![{\displaystyle \Psi (x)\approx C_{0}{\frac {e^{i\int \mathrm {d} x{\sqrt {{\frac {2m}{\hbar ^{2}}}\left(E-V(x)\right)}}+\theta }}{\sqrt[{4}]{{\frac {2m}{\hbar ^{2}}}\left(E-V(x)\right)}}}}](//wikimedia.org/api/rest_v1/media/math/render/svg/436a908df742a9449e6dd7d2e5d355286cce4a9d)
De otro lado, si se varía la fase que varía lentamente (en comparación con la amplitud), (
) entonces
![{\displaystyle A_{0}(x)=\pm {\sqrt {2m\left(V(x)-E\right)}}}](//wikimedia.org/api/rest_v1/media/math/render/svg/eef3dd016a5c710587f5bade739bc5bf291ffb0c)
El cual es solo válido cuando la energía potencial es más grande que la energía total (el régimen en el que túnel cuántico ocurre). rectificando el siguiente orden en el campo de expansión.
![{\displaystyle \Psi (x)\approx {\frac {C_{+}e^{+\int \mathrm {d} x{\sqrt {{\frac {2m}{\hbar ^{2}}}\left(V(x)-E\right)}}}+C_{-}e^{-\int \mathrm {d} x{\sqrt {{\frac {2m}{\hbar ^{2}}}\left(V(x)-E\right)}}}}{\sqrt[{4}]{{\frac {2m}{\hbar ^{2}}}\left(V(x)-E\right)}}}.}](//wikimedia.org/api/rest_v1/media/math/render/svg/1da57426d8f8f2e28d6b65b1bca1c922608c9e47)
Se desprende del denominador, que ambas de las soluciones aproximadas divergen cerca del punto de inflexión clásico donde
y no puede ser válida. Estas son las soluciones aproximadas lejos del valle del potencial y por debajo del valle del potencial. Lejos del valle del potencial, la partícula actúa similarmente a la fase-onde libre que está oscilando. Por debajo del valle de potencial, la partícula sufre cambios exponenciales en amplitud.
Para completar la derivación, las soluciones aproximadas que se encuentran en todas partes y sus coeficientes emparejados se pueden utilizar para obtener una solución aproximada global. La solución aproximada cerca de los puntos de inflexión clásicos
esta aún por ser encontrada.
Para un punto clásico de inflexión
y cerca a
,
puede ser expandida en series de potencias:
![{\displaystyle {\frac {2m}{\hbar ^{2}}}\left(V(x)-E\right)=U_{1}(x-x_{1})+U_{2}(x-x_{1})^{2}+\cdots }](//wikimedia.org/api/rest_v1/media/math/render/svg/21c14e75c73b4b245d6f5e41b4abe4fe23c7d576)
De primer orden se encuentra:
![{\displaystyle {\frac {\mathrm {d} ^{2}}{\mathrm {d} x^{2}}}\Psi (x)=U_{1}(x-x_{1})\Psi (x).}](//wikimedia.org/api/rest_v1/media/math/render/svg/6c0b3e4e19bc22f94512c61d6eb38a7e6ecca66a)
Esta ecuación diferencial es conocida como la ecuación de Airy, y la solución puede ser escrita en términos de la función de Airy.
![{\displaystyle \Psi (x)=C_{A}{\textrm {Ai}}\left({\sqrt[{3}]{U_{1}}}(x-x_{1})\right)+C_{B}{\textrm {Bi}}\left({\sqrt[{3}]{U_{1}}}(x-x_{1})\right).}](//wikimedia.org/api/rest_v1/media/math/render/svg/01cf8de7fdb5aa57e32a7a223ac9125cfecb72c5)
Esta solución debe conectar lejos y por debajo del potencial con las soluciones. Dados los dos coeficientes en un lado del punto de inflexión clásico, los dos coeficientes en el otro lado del punto de inflexión clásico se pueden determinar mediante el uso de esta solución local para conectarlos. Por lo tanto, una relación entre
y
puede ser encontrada.
Afortunadamente, las funciones de Airy divergen en seno, coseno y funciones exponenciales en los límites adecuados. La relación puede ser encontrada para ser de la siguiente manera (a menudo referido como "fórmulas de conexión"):
![{\displaystyle {\begin{aligned}C_{+}&=+{\frac {1}{2}}C_{0}\cos {\left(\theta -{\frac {\pi }{4}}\right)},\\C_{-}&=-{\frac {1}{2}}C_{0}\sin {\left(\theta -{\frac {\pi }{4}}\right)}.\end{aligned}}}](//wikimedia.org/api/rest_v1/media/math/render/svg/6f97b07e26ec7e5998ad8c037ac7aad13b43c8e5)
Ahora las soluciones globales (aproximadas) pueden ser construidas.
Precisión de la serie asintótica
La serie asintótica para
es usualmente una serie divergente cuyos términos generales
comienzan a aumentar después de un cierto valor
. Por lo tanto, el más mínimo error obtenido por el método WKB es de orden menor al último término incluido. Para la ecuación:
con
una función analítica, el valor
y la magnitud del último término puede ser estimada como sigue (ver Winitzki 2005),
![{\displaystyle n_{\max }\approx 2\epsilon ^{-1}\left|\int _{x_{0}}^{x_{*}}dz{\sqrt {-Q(z)}}\right|,}](//wikimedia.org/api/rest_v1/media/math/render/svg/89fe16412e1f6304ed9ee9b28c45827ea2c7eee4)
![{\displaystyle \delta ^{n_{\max }}S_{n_{\max }}(x_{0})\approx {\sqrt {\frac {2\pi }{n_{\max }}}}\exp[-n_{\max }],}](//wikimedia.org/api/rest_v1/media/math/render/svg/3c6d5f3c3f38727f520b9ff7281655890c27ca91)
donde
es el punto en el cual
necesita ser evaluado y
es el punto de inflexión (complejo) donde
, más cerca a
. El número
puede ser interpretado como el número de oscilaciones entre
y el punto de inflexión más cercano. Si
es una función que cambia lentamente,
![{\displaystyle \epsilon \left|{\frac {dQ}{dx}}\right|\ll Q^{2},}](//wikimedia.org/api/rest_v1/media/math/render/svg/5cd7e807c53615065e8c4987bdc3c036d0b9b4b5)
el número
será grande, y el error mínimo de la serie asintótica será exponencialmente pequeño.
Referencias históricas
- Carlini, Francesco (1817). Ricerche sulla convergenza della serie che serva alla soluzione del problema di Keplero. Milano.
- Liouville, Joseph (1837). «Sur le développement des fonctions et séries..». Journal de Mathématiques Pures et Appliquées 1: 16-35.
- Green, George (1837). «On the motion of waves in a variable canal of small depth and width». Transactions of the Cambridge Philosophical Society 6: 457-462.
- Rayleigh, Lord (John William Strutt) (1912). «On the propagation of waves through a stratified medium, with special reference to the question of reflection». Proceedings of the Royal Society London, Series A 86: 207-226. doi:10.1098/rspa.1912.0014.
- Gans, Richard (1915). «Fortplantzung des Lichts durch ein inhomogenes Medium». Annalen der Physik 47: 709-736.
- Jeffreys, Harold (1924). «On certain approximate solutions of linear differential equations of the second order». Proceedings of the London Mathematical Society 23: 428-436. doi:10.1112/plms/s2-23.1.428.
- Brillouin, Léon (1926). «La mécanique ondulatoire de Schrödinger: une méthode générale de resolution par approximations successives». Comptes Rendus de l'Academie des Sciences 183: 24-26.
- Kramers, Hendrik A. (1926). «Wellenmechanik und halbzählige Quantisierung». Zeitschrift der Physik 39: 828-840. doi:10.1007/BF01451751.
- Wentzel, Gregor (1926). «Eine Verallgemeinerung der Quantenbedingungen für die Zwecke der Wellenmechanik». Zeitschrift der Physik 38: 518-529. doi:10.1007/BF01397171.