Aproximación WKB

En Física, la aproximación WKB es el ejemplo más familiar de un cálculo semi-clásico (ver antigua teoría cuántica) en el cual la función de onda es redactada como una función exponencial, semiclásicamente expandida, y entonces o la amplitud o la fase es tomada para estar cambiando lentamente.

El nombre de este método es un acrónimo para la aproximación Wentzel–Kramers–Brillouin. Otro acrónimo usualmente usado para el método incluido Aproximación JWKB and aproximación WKBJ, donde la "J" es debido a Jeffreys.

Contenido

1 Breve historia
2 Método WKB
3 Un ejemplo
4 Aplicación a la ecuación de Schrödinger
- 4.1 Precisión de la serie asintótica
5 Véase también
6 References
- 6.1 Modern references
- 6.2 Historical references
7 Enlaces externos

Breve historia

Este método fue nombrado por los físicos Wentzel, Kramers, y Brillouin, los que desarrollaron esto en 1926. En 1923, el matemático Harold Jeffreys ha desarrollado un método general de aproximación a soluciones lineales, ecuaciones diferenciales de segundo orden, los cuales incluido la ecuación de Schrödinger. Pero desde que la Ecuación de Schrödinger fue desarrollada dos años después, y Wentzel, Kramers, y Brillouin estuvieron aparentemente inconscientes de su trabajo muy adelantado, a Jeffrey algunas veces se le negó el crédito. Textos tempranos en mecánica cuántica contienen algún número de combinaciones de sus iniciales, incluyendo WBK, BWK, WKBJ, JWKB y BWKJ.

Tempranas referencias para el 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.

Método WKB

Generalmente, la teoria WKB es un metodo para aproximar la solución de una ecuación diferencial la cuya mas alta derivada es multiplicada por un pequeño parámetro ε. El método de aproximación es como sigue:

Para una ecuación diferencial:

$\epsilon \frac{\mathrm{d}^ny}{\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$

Asume una solución de la forma de una expansión de serie asintótica

$y(x) \sim \exp\left[\frac{1}{\delta}\sum_{n=0}^{\infty}\delta^nS_n(x)\right]$

En el límite $\delta \rightarrow 0$ . La sustitución de anteriores ansatz dentro de la ecuación diferencial y la cancelación de los tereminos exponenciales permite la solución de un número arbitrario de términos $S n (x)$ en la expansión. La teoría WKB es un caso espefcial de Análisis de escala múltiple.

Un ejemplo

Considerar la ecuación diferencial lineal homogénea de segundo orden.

$\epsilon^2 \frac{d^2 y}{dx^2} = Q(x) y$

donde $Q(x) \neq 0$ . Reemplezando con

$y(x) = \exp\left[\frac{1}{\delta}\sum_{n=0}^{\infty}\delta^nS_n(x)\right]$

resulta en la ecuación

$\epsilon^2\left[\frac{1}{\delta^2}\left(\sum_{n=0}^{\infty}\delta^nS_n'\right)^2 + \frac{1}{\delta}\sum_{n=0}^{\infty}\delta^nS_n''\right] = Q(x)$

De primer orden, (suponiendo, por el momento, la serie será asintóticamente consistente) sobre el se puede aproximar como:

$\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)$

En el límite $\delta \rightarrow 0$ , el equilibrio dominante está dado por:

$\frac{\epsilon^2}{\delta^2}S_0'^2 \sim Q(x)$

Por lo tanto δ es proporcional a ε. Valorizando a ellos igual y comparando podemos hacer:

$\epsilon^0: \; \; \; S_0'^2 = Q(x)$

La cual puede ser reconocida como la ecuación Eikonal, con solución

$S_0(x) = \pm \int_{x_0}^{x}\sqrt{Q(t)}\,dt$

En cuanto a las potencias de primer orden de $\epsilon$ da

$\epsilon^1: \; \; \; 2S_0'S_1' + S_0'' = 0$

La cual es ecuación de transporte unidimensional, el cual tiene la solución

$S_1(x) = -\frac{1}{4}\log\left(Q(x)\right) + k_1.\,$

Y $k 1$ es una constante arbitraria. Nosotros ahora tenemos un par de aproximaciones para el sistema(un par porque $S 0$ puede tomar dos signos); la aproximación-WKB de primer orden será una combinación lineal de las dos:

$y(x) \approx c_1Q^{-\frac{1}{4}}(x)\exp\left[\frac{1}{\epsilon}\int_{x_0}^x\sqrt{Q(t)}dt\right] + c_2Q^{-\frac{1}{4}}(x)\exp\left[-\frac{1}{\epsilon}\int_{x_0}^x\sqrt{Q(t)}dt\right]$

Términos de mayor orden pueden ser obtenidos por examinar en ecuaciones para mas altas potencia de ε. Explicitamente

$2S_0'S_n' + S''_{n-1} + \sum_{j=1}^{n-1}S'_jS'_{n-j} = 0$

para $n > 2$ . Este ejemplo viene de los libros de texto Bender y Orszag (ver referencias).

Aplicación a la ecuación de Schrödinger

La unidimensional, ecuación de Schrödinger independiente del tiempo es: $-\frac{\hbar^2}{2m} \frac{\mathrm{d}^2}{\mathrm{d}x^2} \Psi(x) + V(x) \Psi(x) = E \Psi(x)$ ,

La cual puede ser reescrita como: $\frac{\mathrm{d}^2}{\mathrm{d}x^2} \Psi(x) = \frac{2m}{\hbar^2} \left( V(x) - E \right) \Psi(x)$ .

La función de onda puede ser reescrita como el exponencial de otra función Φ (La cual esta estrechamente relacionada a la acción):

$\Psi(x) = e^{\Phi(x)}, \!$

Así que:

$\Phi''(x) + \left[\Phi'(x)\right]^2 = \frac{2m}{\hbar^2} \left( V(x) - E \right),$

Donde Φ' indica la derivada de Φ con respecto a x. La derivada $Φ'(x)$ puede ser separada dentro de partes reales e imaginarias por la introducción de funciones reales A y B:

$\Phi'(x) = A(x) + i B(x). \;$

La amplitud de la función de onda es entonces $\exp\left[\int^x A(x')dx'\right]\,\!$ , mientras que la fase es $\int^x B(x')dx'\,\!$ . Las partes reales e imaginarias de la ecuación de Schrödinger entonces llegan a ser:

$A'(x) + A(x)^2 - B(x)^2 = \frac{2m}{\hbar^2} \left( V(x) - E \right),$

$B'(x) + 2 A(x) B(x) = 0. \;$

Seguido, la aproximación semiclásica es invocada. Esto significa que cada función como una potente serie en $\hbar$ . Desde la ecuación esto puede ser visto que la serie de potencias puede comenzar con al menos un orden de $\hbar^{-1}$ 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.

$A(x) = \frac{1}{\hbar} \sum_{n=0}^\infty \hbar^n A_n(x)$

$B(x) = \frac{1}{\hbar} \sum_{n=0}^\infty \hbar^n B_n(x)$

De primer orden en esta expansión, las condiciones sobre A y B pueden ser escritas.

$A_0(x)^2 - B_0(x)^2 = 2m \left( V(x) - E \right)$

$A_0(x) B_0(x) = 0 \;$

Si la amplitud varía con la suficiente lentitud en comparación con la fase ( $A 0 (x) = 0$ ), se deduce que:

$B_0(x) = \pm \sqrt{ 2m \left( E - V(x) \right) },$

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:

$\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)}}$

De otro lado, si se varía la fase que varía lentamente (en comparación con la amplitud), ( $B 0 (x) = 0$ ) entonces

$A_0(x) = \pm \sqrt{ 2m \left( V(x) - E \right) }$

El cuál es solo válido cuando la energía potenciañ es mas grande que la energía total (el régimen en el que túnel cuántico ocurre). rectificando el siguiente orden en el campoi de expansión.

$\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)}}.$

Se desprende del denominador, que ambas de las soluciones aproximadas 'explotan' cerca del punto de inflexión clásico donde $E = V (x)$ y no puede ser válida. Estas son las soluciones aproximadas lejos de el valle del potencial y por debajo del valle del potencial. Lejos del valle del potencial, la partícula actua similarmente a la fase-onde libre que está oscilando. por debajo de la villa 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 para obtener una solución aproximada global. La solución aproximada cerca de los puntos de inflexión clásicos $E = V (x)$ estan aón por ser encontrados.

Para un punto clásico de inflexión $x 1$ y cerca a $E = V (x 1)$ , $\frac{2m}{\hbar^2}\left(V(x)-E\right)$ can be expanded in a power series.

$\frac{2m}{\hbar^2}\left(V(x)-E\right) = U_1 (x - x_1) + U_2 (x - x_1)^2 + \cdots$

De primer orden se encuentra:

$\frac{\mathrm{d}^2}{\mathrm{d}x^2} \Psi(x) = U_1 (x - x_1) \Psi(x).$

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.

$\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).$

Esta solución debe conectar el lejos y por debajo de soluciones. Dadas los dos coeficientes en un lado del punto de inflexión clásico, los 2 coeficientes en el otro lado del punto de inflexión clásico se puede determinar mediante el uso de esta solución local para conectarlos. Por lo tanto, una relación entre $C 0,θ$ y $C +, C -$ puede ser encotnrada.

Afortunadamente, las funciones de Airy se asíntota en seno, coseno y funciones exponenciales en los límites adecuados. La relación puede ser encontrado para ser de la siguiente manera (a menudo referido como "fórmulas de conexión"):

$\begin{align} 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{align}$

Ahora las soluciones globales (aproximadas) pueden ser construidas.

Precisión de la serie asintótica

La serie asintótica para $y (x)$ es usualmente una serie divergente cuyos término general $δ n S n (x)$ comienza a aumentar después de un cierto valor $n = n max$ . Por lo tanto el más mínimo error obtenidos por el método WKB es en el mejor de la orden del último término incluido. Para la ecuación: $\epsilon^2 \frac{d^2 y}{dx^2} = Q(x) y$ con $Q (x) < 0$ una función analítica, el valor $n max$ y la magnitud del último termino puede ser estimado como sigue (ver Winitzki 2005),

$n_\max \approx 2\epsilon^{-1} \left| \int_{x_0}^{x_{*}} dz\sqrt{-Q(z)} \right| ,$

$\delta^{n_\max}S_{n_\max}(x_0) \approx \sqrt{\frac{2\pi}{n_\max}} \exp[-n_\max],$

donde $x 0$ es el punto en el cual $y (x 0)$ necesita ser evaluado y $x *$ es el punto de inflexión (complejo) donde $Q (x *) = 0$ , mas cerca a $x = x 0$ . El número $n max$ puede ser interpretado como el número de oscilaciones entre $x 0$ y el punto de inflexión mas cercano. Si $\epsilon^{-1}Q(x)$ es una función que cambia lentamente,

$\epsilon\left| \frac{dQ}{dx} \right| \ll Q^2 ,$

el número $n max$ será largo, y el error mínimo de la serie asintótica será exponencialmente pequeño.

Véase también

Función de Airy
Field electron emission
Langer correction
Method of steepest descent / Laplace Method
Old quantum theory
Perturbation methods
Quantum tunneling
Slowly varying envelope approximation

References

Modern references

Razavy, Moshen (2003). Quantum Theory of Tunneling. World Scientific. ISBN 981-238-019-1.
Griffiths, David J. (2004). Introduction to Quantum Mechanics (2nd ed.). Prentice Hall. ISBN 0-13-111892-7.
Liboff, Richard L. (2003). Introductory Quantum Mechanics (4th ed.). Addison-Wesley. ISBN 0-8053-8714-5.
Sakurai, J. J. (1993). Modern Quantum Mechanics. Addison-Wesley. ISBN 0-201-53929-2.
Bender, Carl; Orszag, Steven (1978). Advanced Mathematical Methods for Scientists and Engineers. McGraw-Hill. ISBN 0-07-004452-X.
Olver, Frank J. W. (1974). Asymptotics and Special Functions. Academic Press. ISBN 0-12-525850-X.
Winitzki, Sergei (2005). «Cosmological particle production and the precision of the WKB approximation». Physical Review D 72: pp. 104011. doi:10.1103/PhysRevD.72.104011.

Historical references

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: pp. 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: pp. 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: pp. 207–226. doi:10.1098/rspa.1912.0014.
Gans, Richard (1915). «Fortplantzung des Lichts durch ein inhomogenes Medium». Annalen der Physik 47: pp. 709–736.
Jeffreys, Harold (1924). «On certain approximate solutions of linear differential equations of the second order». Proceedings of the London Mathematical Society 23: pp. 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: pp. 24–26.
Kramers, Hendrik A. (1926). «Wellenmechanik und halbzählige Quantisierung». Zeitschrift der Physik 39: pp. 828–840. doi:10.1007/BF01451751.
Wentzel, Gregor (1926). «Eine Verallgemeinerung der Quantenbedingungen für die Zwecke der Wellenmechanik». Zeitschrift der Physik 38: pp. 518–529. doi:10.1007/BF01397171.

Enlaces externos

Richard Fitzpatrick, The W.K.B. Approximation (2002). (An application of the WKB approximation to the scattering of radio waves from the ionosphere.)
Free WKB library for Microsoft Visual C v6 for some special functions

Categoría:

Mecánica cuántica

Wikimedia foundation. 2010.

Игры ⚽ Поможем решить контрольную работу

Mira otros diccionarios:

Teoría perturbacional — En mecánica cuántica, la teoría perturbacional o teoría de perturbaciones es un conjunto de esquemas aproximados para describir sistemas cuánticos complicados en términos de otros más sencillos. La idea es empezar con un sistema simple y… … Wikipedia Español
Efecto túnel — Saltar a navegación, búsqueda Reflexión y tunelado de un electrón dirigido hacia una barrera potencial. El punto resplandeciente moviéndose de derecha a izquierda es la sección reflejada del wavepacket. Un vislumbre puede observarse a la derecha… … Wikipedia Español
George Green — Portada de un ensayo publicado en 1828 por George Green. George Green (julio de 1793, 31 de mayo de 1841) fue un matemático británico cuyo trabajo influenció notablemente el desarrollo de importantes conceptos en física. Entre sus obras más… … Wikipedia Español

Los diccionarios y las enciclopedias sobre el Académico

Aproximación WKB

Contenido

Breve historia

Método WKB

Un ejemplo