Barreras de potencial
Consideremos el potencial representado en la figura que consta de dos escalones y que se denomina barrera de potencial de altura V0 y anchura a.
El caso más interesante se da, cuando la energía de las partículas sea menor que la de la barrera. La Mecánica Clásica requiere que una partícula proveniente de la izquierda con E<V0 se refleje en el origen x=0, ya que en la región (0, a) la energía cinética de la partícula es negativa.
Las partículas que hayan penetrado una distancia mayor o igual que a, tendrían una energía cinética igual a su energía total (la energía potencial vuelve a ser cero) y por tanto, se moverán hacia la derecha con igual velocidad que las incidentes. Estas partículas que han atravesado la barrera se denominan transmitidas y han pasado de la primera a la tercera región de potencial a través de la región intermedia, clásicamente prohibida (la energía cinética de la partícula es negativa).
Una barrera de potencial
Discutiremos ahora el problema desde el punto de vista de la Mecánica Cuántica, resolviendo la ecuación de Schrödinger en las tres regiones y aplicando las condiciones de continuidad de la función de onda y de su derivada primera en los puntos x=0, y x=a.
Resolveremos primero, el caso en el que la energía de las partículas E es menor que la del escalón V0, el caso más interesante desde el punto de vista físico. Posteriormente, estudiamos el caso en el que la energía de la partícula E es mayor que la del escalón E>V0.
E<V0
Región I, x<0
Región II, 0<x<a, aquí E<V0
Región III, x>a
La función de onda Ψ1(x) contiene las partículas incidentes y reflejadas, Ψ2(x) decrece exponencialmente, la exponencial positiva no está excluida ya que la región clásicamente prohibida no es indefinida como en el caso del escalón de potencial. Debido a que Ψ2(x) no ha alcanzado el valor cero en x=a, la función de onda continúa a la derecha de dicho punto, con amplitud F. La función de onda Ψ3(x) representa las partículas transmitidas.
Desde el punto de vista de la Mecánica Cuántica, es posible que una partícula atraviese la barrera de potencial aún cuando su energía cinética sea menor que la altura de la barrera.
Aplicamos las condiciones de continuidad de la función de onda y de su derivada primera en los puntos x=0
Aplicamos las condiciones de continuidad de la función de onda y de su derivada primera en los puntos x=a
Para calcular el coeficiente de transmisión T=|F/A|2 tenemos que relacionar el coeficiente F con A
A partir, de las dos primeras ecuaciones despejamos C y D
A partir, de las otras dos ecuaciones despejamos C y D
Llegamos al sistema de dos ecuaciones y despejamos F
El cuadrado de su módulo es
Se denomina coeficiente de transmisión a la proporción de partículas incidentes que son transmitidas
Aproximación
Cuando la barrera es alta (V0 grande) y estrecha (a pequeño), y la energía de las partículas E es pequeña, aproximamos sinh(qa)≈qa
Un resultado similar se obtiene para una barrera de potencial delta de Dirac
El coeficiente de transmisión disminuye rápidamente a medida que se incrementa la anchura de la barrera de potencial.
Representamos el coeficiente de transmisión T para una barrera de anchura a=0.5 y de altura V0=2
a=0.5; %anchura barrera V0=2; %altura ee=linspace(0,V0,200); T=zeros(1,length(ee)); i=1; for E=ee q=sqrt(V0-E); T(i)=4*E*(V0-E)/(4*E*(V0-E)+V0^2*sinh(q*a)^2); i=i+1; end plot(ee,T); grid on ylim([0,1.1]) xlabel('E') ylabel('T') title('Coeficiente de transmisión')
Como ejercicio adicional, se despeja el coeficiente B y se calcula el coeficiente de reflexión R, verificando que R+T=1
Función de onda
Representamos la función de onda para una barrera de potencial de altura E0=10 y anchura a=1. La energía E=5
a=1; %anchura barrera V0=10; %altura E=5; %energía menor que V0 k=sqrt(E); q=sqrt(V0-E); A=1; F=2*1i*k*q*exp(-1i*k*a)*A/(2*1i*k*q*cosh(q*a)-(q^2-k^2)*sinh(q*a)); C=(1+1i*k/q)*F*exp(1i*k*a)*exp(-q*a)/2; D=(1-1i*k/q)*F*exp(1i*k*a)*exp(q*a)/2; B=((1-1i*k/q)*F*exp(1i*k*a)*exp(q*a)-(1-1i*k/q)*A)/(1+1i*k/q); hold on fplot(@(x) real(A*exp(1i*k*x)+B*exp(-1i*k*x)),[-10,0]) %incidente+reflejada %fplot(@(x) real(A*exp(1i*k*x)),[-10,0]) %incidente %fplot(@(x) real(B*exp(-1i*k*x)),[-10,0]) %reflejada fplot(@(x) real(D*exp(-q*x)+C*exp(q*x)),[0,a]) fplot(@(x) real(F*exp(1i*k*x)),[a,10]) %transmitida hold off grid on xlabel('x') ylabel('\Phi(x)') title('Función de onda')
Se muestra la función de onda, en la región I (color azul), II (color rojo) y III (color amarillo)
En la barrera de potencial la función de onda tiene la forma Ψ2(x)=Cexp(qx)+Dexp(-qx). Comprobamos que el coeficiente C es pequeño comparado con D, por lo que la función de onda en esta región es casi una exponencial decreciente con x
>> C,D C = 0.0113 + 0.0113i D = 0.9887 - 0.9887i
Separamos la parte incidente de la función de onda de la reflejada, en región I
E>V0
Para E<V0, el coeficiente de transmisión T es menor que la unidad. Sin embargo, para E>V0, T alcanza el valor máximo, para valores concretos del cociente E/V0.
Resolviendo de nuevo, la ecuación de Schrödinger en las tres regiones
Región I, x<0
Región II, 0<x<a, ahora E>V0.
Región III, x>a
Las ecuaciones de continuidad de la función de onda y de su derivada primera en x=a, relacionan C y D con F, y en x=0, relacionan A y B con C y D, y por tanto, con F
Aplicamos las condiciones de continuidad de la función de onda y de su derivada primera en los puntos x=0
Aplicamos las condiciones de continuidad de la función de onda y de su derivada primera en los puntos x=a
A partir, de las dos primeras ecuaciones despejamos C y D
A partir, de las otras dos ecuaciones despejamos C y D
Llegamos al sistema de dos ecuaciones y despejamos F
Finalmente, obtenemos la siguiente expresión para el coeficiente de transmisión
T toma el valor máximo 1, cuando qa=nπ, siendo n un número entero. Como q es el número de onda, q=2π/λ, se obtiene que
que relaciona la longitud de onda λ de la partícula en la barrera de potencial con la anchura a de la misma, para que se obtenga el máximo en el coeficiente de transmisión. Los valores de la energía E, o mejor del cociente E/V0 ,para los cuales hay un máximo del coeficiente de transmisión se denominan resonancias.
Representamos el coeficiente de transmisión T para una barrera de anchura a=1 y de altura V0=10
a=1; %anchura barrera V0=10; %altura ee=linspace(V0,50,200); T=zeros(1,length(ee)); i=1; for E=ee k=sqrt(E); q=sqrt(E-V0); T(i)=4*k^2*q^2/((2*k*q*cos(q*a))^2+(k^2+q^2)^2*sin(q*a)^2); i=i+1; end plot(ee,T); E1=(pi/a)^2+V0; line([E1,E1],[0,1],'lineStyle','--') E1=(2*pi/a)^2+V0; line([E1,E1],[0,1],'lineStyle','--') grid on ylim([0,1.1]) xlabel('E') ylabel('T') title('Coeficiente de transmisión')
Función de onda
Representamos la función de onda para una barrera de potencial de altura E0=10 y anchura a=1. La energía E=5
a=1; %anchura barrera V0=10; %altura E=15; %energía mayor que V0 k=sqrt(E); q=sqrt(E-V0); A=1; F=2*k*q*exp(-1i*k*a)*A/(2*k*q*cos(q*a)-(q^2+k^2)*1i*sin(q*a)); C=(1-k/q)*F*exp(1i*k*a)*exp(1i*q*a)/2; D=(1+k/q)*F*exp(1i*k*a)*exp(-1i*q*a)/2; B=((1+k/q)*F*exp(1i*k*a)*exp(-1i*q*a)-(1+k/q)*A)/(1-k/q); hold on fplot(@(x) real(A*exp(1i*k*x)+B*exp(-1i*k*x)),[-10,0]) fplot(@(x) real(C*exp(-1i*q*x)+D*exp(1i*q*x)),[0,a]) fplot(@(x) real(F*exp(1i*k*x)),[a,10]) hold off grid on xlabel('x') ylabel('\Phi(x)') title('Función de onda')
Se muestra la función de onda, en la región I (color azul), II (color rojo) y III (color amarillo)
Separamos la parte incidente de la función de onda en región I de la reflejada
a=1; %anchura barrera V0=10; %altura E=15; %energía mayor que V0 k=sqrt(E); q=sqrt(E-V0); A=1; F=2*k*q*exp(-1i*k*a)*A/(2*k*q*cos(q*a)-(q^2+k^2)*1i*sin(q*a)); C=(1-k/q)*F*exp(1i*k*a)*exp(1i*q*a)/2; D=(1+k/q)*F*exp(1i*k*a)*exp(-1i*q*a)/2; B=((1+k/q)*F*exp(1i*k*a)*exp(-1i*q*a)-(1+k/q)*A)/(1-k/q); hold on fplot(@(x) real(A*exp(1i*k*x)),[-10,0]) fplot(@(x) real(B*exp(-1i*k*x)),[-10,0]) fplot(@(x) real(C*exp(-1i*q*x)+D*exp(1i*q*x)),[0,a]) fplot(@(x) real(F*exp(1i*k*x)),[a,10]) hold off grid on xlabel('x') ylabel('\Phi(x)') title('Función de onda')
Dos barreras de potencial
Consideremos el potencial representado en la figura que consta de dos barreras de potencial de anchura a y altura V0 separadas una distancia b.
Resolveremos el caso en el que la energía de las partículas E es menor que la del escalón V0, el caso más interesante desde el punto de vista físico.
E<V0
Región I, x<0
Región II, 0<x<a, aquí E<V0
Región III, a<x<a+b
Región IV, a+b<x<2a+b
Región V, x>2a+b
Aplicamos las condiciones de continuidad de la función de onda y de su derivada primera en los puntos:
x=0
x=a
x=a+b
x=2a+b
Para calcular el coeficiente de transmisión T=|A5/A1|2 tenemos que relacionar el coeficiente A5 con A1
El coeficiente de transmisión T
Representamos el coeficiente de transmisión T para una barrera doble de potencial de anchura a=1, separación b=2 y de altura V0=5
b=2; %separación entre barreras a=1; %anchura barrera V0=5.0; %altura barrera E=linspace(0,5,200); %energía trans=zeros(length(E),1); for j=1:length(E) k=sqrt(E(j)); q=sqrt(V0-E(j)); beta=q/k; gamma=(1+beta^2)/(2*beta); eta=(beta^2-1)/(2*beta); u=(gamma^2*cos(2*k*b)-eta^2)*sinh(q*a)^2+cosh(q*a)^2; v=(gamma^2*sin(2*k*b)*sinh(q*a)+2*eta*cosh(q*a))*sinh(q*a); trans(j)=1/(u^2+v^2); end plot(E,trans); grid on ylim([0 1.1]) xlabel('E'); ylabel('T_r') title('Coeficiente de transmisión, N=2')
N barreras de potencial
Estudiaremos ahora el caso en el que hay N barreras de potencial de la misma anchura a y separadas unas de otras la misma cantidad b tal como se aprecia en la figura. Observaremos que se producen picos de resonancia adicionales, dando lugar a un comportamiento complejo del coeficiente de transmisión.
Resolvemos la ecuación de Schrödinger para cada una de las regiones, 1, 2...2N+1
Donde qj es un número complejo real o imaginario dependiendo de que E>Vj o E<Vj.
En las fronteras entre las regiones, aplicamos las condiciones de continuidad. Sea la frontera entre las regiones j-1 y j, cuya abscisa es xj.
que relaciona los coeficientes Aj y Bj con Aj+1 y Bj+1
Teniendo en cuenta que solamente hay partículas trasmitidas en la región 2N+1, resulta que B2N+1=0. Obtenemos los valores de todos los coeficientes Aj y Bj en términos de A2N+1 que actúa como factor de escala
El coeficiente de transmisión se define como la proporción de partículas incidentes que se transmiten y se obtiene mediante el cociente.
En la figura, vemos el comportamiento del coeficiente de transmisión, para N=2 barreras de potencial, de anchura a=1, separación b=2 y altura V0=5.
b=2; %separación entre barreras a=1; %anchura barrera N=2; %número de barreras V0=5.0; %altura barrera x=zeros(2*N,1); V=zeros(2*N+1,1); q=zeros(2*N+1,1); x(1)=-N*(a+b)/2+b/2; for i=2:2:2*N x(i)=x(i-1)+a; x(i+1)=x(i)+b; end %potenciales for i=1:2:2*N+1 V(i)=0.0; if i==2*N+1 break; end V(i+1)=V0; end f=@(z,q) [exp(1i*q*z) exp(-1i*q*z); q*exp(1i*q*z) -q*exp(-1i*q*z)]; g=@(z,q) [exp(-1i*q*z) exp(-1i*q*z)/q; exp(1i*q*z) -exp(1i*q*z)/q]/2; E=linspace(0,V0,200); trans=zeros(length(E),1); for j=1:length(E) for i=1:2*N+1 q(i)=sqrt(E(j)-V(i)); end A=[1;0]; for i=2*N:-1:1 A=g(x(i),q(i))*f(x(i),q(i+1))*A; end trans(j)=1/abs(A(1))^2; end plot(E,trans); grid on ylim([0 1.1]) xlabel('E'); ylabel('T_r') title('Coeficiente de transmisión, N=2')
La misma figura que hemos obtenido anteriormente, a partir de un cálculo explícito del coeficiente de transmisión T
En la figura, vemos el comportamiento del coeficiente de transmisión, para N=3 barreras de potencial, de anchura a=1, separación b=2 y altura V0=5.
Referencias
Zbigniew Ficek. Problems and Solutions in Quantum Physics. CRC Press. Problem 9.1, pp. 75-79