Capítulo 9 O POÇO RETANGULAR FINITO

9.1 Introdução

O poço retangular finito humaniza o poço retangular infinito do [ capítulo 6 ]. Uma parede de potencial infinita é uma idealização matemática. Em alguns casos, pode até ser usada como aproximação. Mas, para ser mais realístico, devemos considerar uma parede finita.

Podemos imaginar o perfil de um poço retangular finito pensando em uma barreira retangular invertida [ capítulo 8 ]. Se uma barreira é invertida, um poço é formado.

O poço retangular finito tem, em suas margens, que chamamos de barreiras, um valor de potencial maior que o valor do potencial que cobre a extensão do próprio poço. Uma maneira de construirmos seu perfil, é igualar à zero a região das barreiras, esquerda e direita, e fazer a região do poço igual a (V0), sendo V0 uma quantidade positiva:

V(x)=0,se  x<0,=V0,se  0xa,=0,se  x>a.(9.1)

Visualizamos o poço retangular finito na Fig. 9.1. Observe que a região do poço se estende de x=0 até x=a.

O poço retangular finito.

Figura 9.1: O poço retangular finito.

Tomando como referencial a energia potencial das barreiras, no caso, V=0, a energia da partícula que interage com o poço retangular finito, da Fig. 9.1, pode assumir valores positivos ou negativos. Nosso primeiro objetivo será estudar a partícula com energia maior que as barreiras, E>0 [ seção 9.2 ]. Depois, vamos tratar da partícula com energia menor que as barreiras e maior que o fundo do poço, V0<E<0 [ seção 9.3 ].

9.2 A transmissão do estado não-ligado

A partícula que vamos estudar nessa seção possui energia maior que a energia potencial das barreiras (B) do poço (P) retangular finito, quer dizer, a partícula interage com energia positiva". No regime E>0, o estado quântico dessa partícula é do tipo não-ligado, o que proporciona condições para ela ser refletida ou transmitida pelo poço de potencial.

A solução da equação de Schrödinger segue a mesma estrutura da [ seção 8.3 ]:

ψB=Aeikx+Beikx,x<0,ψP=CeikPx+DeikPx,0xa,ψB=Feikx,x>a.(9.2)(9.3)(9.4)

Em cada região, o número de onda apresenta uma peculiar expressão:

k=2mE,x<0,kP=2m(E+V0),0xa,k=2mE,x>a.(9.5)(9.6)(9.7)

Ao compararmos, as soluções e os números de onda acima, com os correspondentes desenvolvidos na [ seção 8.3 ], fica evidente a perfeita identidade estrutural. Por essa razão, o coeficiente de transmissão do poço retangular finito compartilha a mesma expressão do coeficiente de transmissão da barreira retangular no regime de alta energia:

1T=1+(k2k2P)24k2k2Psen2(kPa).(9.8)

Nota: O que precisa ser destacado é que o número de onda da barreira retangular, na região 0xa, é kD=2m(EV0)/, enquanto que o número de onda do poço é dado por (9.6).

Investigando a expressão do coeficiente (9.8), transmissões com o valor mávimo T=1 ocorrem para sen(kPa)=0, ou seja:

kP=nπa,     n=1,2,3...(9.9)

Substituindo (9.9) em (9.6), encontramos a expressão das energias que geram 100% de transmissão:

En=n2π222ma2V0.(9.10)

Restrição aos valores de n: Por causa da condição a partícula interage com energia positiva" e também por (9.10) ser uma subtração de dois termos, os inteiros válidos devem conferir ao primeiro termo de (9.10) um resultado maior que V0, quer dizer, todo n que resulte em E<0 deve ser descartado.

  Um exemplo numérico.

Considere um elétron interagindo com um poço retangular finito, de largura a=100 Å e profundidade V0=100 meV. Os gráficos dos coeficientes de transmissão e reflexão, em função da energia do elétron, são apresentados na Fig. 9.2.

Os coeficientes de transmissão e reflexão do poço retangular finito discutido no exemplo.

Figura 9.2: Os coeficientes de transmissão e reflexão do poço retangular finito discutido no exemplo.

Conforme previsto pela equação de ressonância (9.10), há transmissões máximas no decorrer da curva do coeficiente de transmissão. As energias de ressonância foram impressas na própria figura, ao longo da linha pontilhada T=1. O coeficiente de reflexão foi determinado fazendo R=1T. Quando T é máximo, R é zero, confirmando que não há reflexão.

Há pontos intermediários, por exemplo, quando E=56 meV, temos T=0,78 e R=0,22. Isso significa que há 78% de chance do elétron atravessar a região do poço, e 22%, do elétron não atravessar a região do poço. Se um físico clássico" fosse explicar a reflexão, diria: Um elétron incidindo com energia maior que a energia das barreiras de um poço, não pode sofrer reflexão por esse poço". Mas um físico quântico" garante que há 22% de chance do elétron ser refletido pelo poço. No ponto E=19 meV, ocorre a intersecção das curvas, daí T=R=50%.

9.3 O confinamento do estado ligado

Na [ seção 9.2 ], a partícula interage com energia positiva". Nessa seção, ela interage com energia menor que as barreiras e maior que o fundo do poço, V0<E<0, quer dizer, agora a partícula interage com energia negativa".

A solução da equação de Schrödinger segue a mesma estrutura da [ seção 9.2 ]:

ψB=Aeikx+Beikx,x<0,ψP=CeikPx+DeikPx,0xa,ψB=Feikx,x>a.(9.11)(9.12)(9.13)

Em cada região, o número de onda apresenta uma peculiar expressão:

k=2mE,x<0,kP=2m(E+V0),0xa,k=2mE,x>a.(9.14)(9.15)(9.16)

Visto que E<0, os números de onda nas regiões (9.14) e (9.16) são imaginários. É favorável escrevê-los na forma: k=iq, onde q é um número real e positivo:

q=2m(0E).(9.17)

Com o auxílio de (9.17), as soluções (9.11) e (9.13) se transformam em exponenciais reais. Ademais, fazendo x, percebe-se que a exponencial A é divergente. Daí vem a necessidade de ajustar A=0. Após essas considerações, as soluções mudam para:

ψB=Beqx,x<0,ψP=CeikPx+DeikPx,0xa,ψB=Feqx,x>a.(9.18)(9.19)(9.20)

Nas interfaces, há necessidade de se aplicar as condições de continuidade da função de onda:

Interface x=0 ψB(0)=ψP(0)dψB(0)dx=dψP(0)dx(9.21)

Interface x=a ψP(a)=ψB(a)dψP(a)dx=dψB(a)dx(9.22)

Resolvendo (9.21) e (9.22), encontramos:

Interface x=0 {B=C+DqB=ikPCikPD(9.23)

Interface x=a {CeikPa+DeikPa=FeqaikPCeikPaikPDeikPa=qFeqa(9.24)

Ao resolver o sistema de equações (9.23):

2C=(ikP+qikP)B,2D=(ikPqikP)B.(9.25)

E ao resolver o sistema de equações (9.24):

2CeikPa=(ikPqikP)Feqa,2DeikPa=(ikP+qikP)Feqa.(9.26)

A utilização de 2C e 2D, impressos em (9.25), nas relações (9.26), resulta em:

(ikP+qikP)BeikPa=(ikPqikP)Feqa,(ikPqikP)BeikPa=(ikP+qikP)Feqa.(9.27)

Com um pouco de matemática em (9.27), podemos isolar o termo F:

(ikP+qikPq)BeikPa=Feqa,(ikPqikP+q)BeikPa=Feqa.(9.28)

Enfim, dividindo as equações (9.28), encontramos:

(ikPqikP+q)2=e2ikPa.(9.29)

A igualdade (9.29) é verdadeira se:

ikPqikP+q=±eikPa.(9.30)

A igualdade (9.30) nos oferece duas possibilidades, de acordo com a escolha do sinal ±. Vamos, então, analisar os dois casos.

9.3.1 O caso tangente, transcendental

Vamos analisar o caso (9.30) com o sinal (+).

Primeiro, dividimos as equações (9.25):

DC=ikPqikP+q.(9.31)

Agora usamos (9.30) em (9.31):

D=CeikPa.(9.32)

E agora usamos (9.32) em (9.24), referente à interface x=a:

CeikPa+CeikPaeikPa=Feqa,ikPCeikPaikPCeikPaeikPa=qFeqa.(9.33)

Efetuando a matemática de (9.33):

C(eikPa+1)=Feqa,C(eikPa1)=qikPFeqa.(9.34)

Dividimos as equações (9.34):

eikPa1eikPa+1=qikP.(9.35)

Identificamos o lado esquerdo de (9.35) como sendo a seguinte relação trigonométrica:

eikPa1eikPa+1=itg(kPa2).(9.36)

Por fim, levando (9.36) ao resultado (9.35):

tg(kPa2)=qkP.(9.37)

Observe que a variável kP se encontra em ambos lados da equação (9.37). Uma equação desse tipo é chamada transcendental. Equações transcendentais vão além dos limites da álgebra e somente podem ser resolvidas numericamente. Aqui temos uma transcendental envolvendo a função tangente.

9.3.2 O caso cotangente, transcendental

Agora vamos analisar o caso (9.30) com o sinal ().

Seguindo o passo a passo da [ seção 9.3.1 ], temos:

D=CeikPa.(9.38)

eikPa+1eikPa1=qikP.(9.39)

Identificamos o lado esquerdo de (9.39) como:

eikPa+1eikPa1=icotg(kPa2).(9.40)

Por fim, levando (9.40) ao resultado (9.39):

cotg(kPa2)=qkP.(9.41)

Então terminamos com uma equação transcendental envolvendo a função cotangente.

9.3.3 Os níveis de energia

Conforme deduzido nas [ seções 9.3.1 e 9.3.2 ], o poço retangular finito gera duas equações transcendentais, as quais necessitam de métodos numéricos para serem resolvidas. A energia envolvida é negativa, V0<E<0. É oportuno trabalhar com variáveis positivas. Isso pode ser feito seguindo o raciocínio ilustrado na Fig. 9.3. Ao somar V0 nas energias E e V do perfil do lado esquerdo de 9.3, chega-se ao perfil do lado direito de 9.3, que toma como referência o fundo do poço.

Definição de nível de energia.

Figura 9.3: Definição de nível de energia.

A nova variável (E+V0) é chamada de nível de energia, e será simbolizada por:

ε=E+V0.(9.42)

Logo, o número de onda (9.15) muda para:

kP=2mε.(9.43)

E o coeficiente (9.17) passa para:

q=2m(V0ε).(9.44)

Também, vale ressaltar que as energias percorrem o intervalo V0<E<0, e os valores dos níveis de energia, 0<ε<V0.

9.3.4 Um exemplo numérico

Considere o mesmo poço retangular finito utilizado na [ seção 9.2 ], aquele de largura a=100 Å e profundidade V0=100 meV. Nele, um elétron está confinado, por possuir energia menor que as barreiras.

Sabemos, do estudo do poço infinito [ capítulo 6 ], que as energias permitidas para um elétron confinado são discretas, quantizadas. Inclusive, na [ seção 6.2 ], deduzimos uma fórmula para que possam ser calculadas. Mas, no caso atual, do poço finito, não há fórmula, então, como determinar as energias? O cálculo é numérico, substituiremos os dados do nosso poço nas equações transcendentais (9.37) e (9.41):

tg(kPa2)=qkP,cotg(kPa2)=qkP.(9.45)

A Fig. 9.4 apresenta as curvas das equações (9.45). Segundo (9.43) e (9.44), kP e q dependem de ε, então, fizemos 0<ε<V0. Como se vê, as curvas tangente e cotangente cruzam a curva q/kP. São nas intersecções que aparecem os valores dos níveis de energia, sinalizados por ε1 até ε6. Nota-se que ε1, ε3 e ε5 são determinados pela curva tangente, e que ε2, ε4 e ε6, pela cotangente. Os valores dos níveis de energia εn e das energias de confinamento En são mostrados logo abaixo da figura.

Curvas das equações transcendentais do poço retangular finito discutido no exemplo numérico.

Figura 9.4: Curvas das equações transcendentais do poço retangular finito discutido no exemplo numérico.

Nível de energia (meV) Energia (meV)
ε1=3,0 E1=97,0
ε2=11,9 E2=88,1
ε3=26,6 E3=73,4
ε4=46,8 E4=53,2
ε5=71,6 E5=28,4
ε6=97,9 E6=2,1

É de se esperar que os valores dos níveis de energia do poço finito tendam aos valores dos níveis de energia do poço infinito, quando V0. Para comprovar essa hipótese, determinamos os níveis de energia de um poço finito, de largura a=100 Å e profundidade V0=106 meV, e comparamos com os valores dos níveis de energia de um poço infinito, de mesma largura. A Fig. 9.5 mostra a comparação. Fica claro que os valores do poço finito 106 meV se aproximam muito dos valores do poço infinito. Na mesma figura, acrescentamos os valores dos níveis de energia do poço finito da Fig. 9.4.

Comparação energética entre 3 poços quânticos.

Figura 9.5: Comparação energética entre 3 poços quânticos.

9.4 A densidade de probabilidade da posição

O poço retangular finito, de largura a, divide o espaço em três regiões: há a região do próprio poço, 0xa, a região da barreira da esquerda, x<0, e a região da barreira da direita, x>a. A função de onda nessas regiões tem o aspecto:

ψB=Beqx,x<0,ψP=CeikPx+DeikPx,0xa,ψB=Feqx,x>a.

Nosso objetivo é determinar uma expressão para a densidade de probabilidade da posição. Para isso, vamos utilizar alguns resultados da [ seção 9.3 ]. Antes, porém, precisamos de algumas relações trigonométricas:

2cosz=eiz+eiz,(2cosz)(2cosz)=(eiz+eiz)(eiz+eiz),4cos2z=2+e2iz+e2iz.(9.46)

2isenz=eizeiz,(2isenz)(2isenz)=(eizeiz)(eizeiz),4sen2z=2e2ize2iz.(9.47)

9.4.1 Dentro do poço

A função de onda dentro do poço é:

ψP=CeikPx+DeikPx.(9.48)

Se ela for vinculada ao caso tangente, D=CeikPa [ seção 9.3.1 ]:

ψP=C(eikPx+eikPaeikPx).(9.49)

A partir de (9.49), determina-se a densidade de probabilidade da posição:

|ψP|2=|C|2(2+e2ikPxeikPa+e2ikPxeikPa).(9.50)

Então, recorrendo à trigonometria (9.46):

|ψP|2=4|C|2cos2(kPx12kPa)    ε1,ε3,ε5.(9.51)

  Todavia, se a função de onda (9.48) for vinculada ao caso cotangente, D=CeikPa [ seção 9.3.2 ]:

ψP=C(eikPxeikPaeikPx).(9.52)

E a densidade de probabilidade fica:

|ψP|2=|C|2(2e2ikPxeikPae2ikPxeikPa).(9.53)

Então, recorrendo à trigonometria (9.47):

|ψP|2=4|C|2sen2(kPx12kPa)    ε2,ε4,ε6.(9.54)

Queremos destacar o sinal de () que aparece em (9.52), ele esconde uma fase:

ψP=C(eikPx+eiπeikPaeikPx).(9.55)

Logo, no caso dos estados pares, a partícula ganha uma fase adicional de 1800 ao passar pelo processo de reflexão na interface.

9.4.2 Na barreira da esquerda

Na barreira da esquerda, a função de onda é:

ψB=Beqx.(9.56)

Ela deve respeitar a condição de continuidade em x=0, B=C+D, além do mais, se for vinculada ao caso tangente, D=CeikPa:

ψB=C(1+eikPa)eqx.(9.57)

Daí,

|ψB|2=|C|2(2+eikPa+eikPa)e2qx.(9.58)

Ou seja,

|ψB|2=4|C|2cos2(12kPa)e2qx    ε1,ε3,ε5.(9.59)

  Contudo, se a função de onda for vinculada ao caso cotangente, D=CeikPa:

ψB=C(1eikPa)eqx.(9.60)

Daí,

|ψB|2=|C|2(2eikPaeikPa)e2qx.(9.61)

Ou seja,

|ψB|2=4|C|2sen2(12kPa)e2qx    ε2,ε4,ε6.(9.62)

9.4.3 Na barreira da direira

Na barreira da direita, a função de onda é:

ψB=Feqx.(9.63)

Ela deve respeitar a condição de continuidade em x=a:

CeikPa+DeikPa=Feqa.(9.64)

Então,

F=(CeikPa+DeikPa)eqa.(9.65)

Podemos usar em (9.65), o resultado do caso tangente: D=CeikPa:

F=C(eikPa+1)eqa.(9.66)

Surge uma nova expressão para (9.63):

ψB=C(eikPa+1)eq(xa).(9.67)

A partir de (9.67), determina-se a densidade de probabilidade da posição:

|ψB|2=|C|2(2+eikPa+eikPa)e2q(xa).(9.68)

Então, recorrendo à trigonometria (9.46):

|ψB|2=4|C|2cos2(12kPa)e2q(xa)    ε1,ε3,ε5.(9.69)

  Entretanto, podemos usar em (9.65), o resultado do caso cotangente, D=CeikPa:

F=C(eikPa1)eqa.(9.70)

E isso modifica a função de onda (9.63):

ψB=C(eikPa1)eq(xa).(9.71)

Daí,

|ψB|2=|C|2(2eikPaeikPa)e2q(xa).(9.72)

Então, recorrendo à trigonometria (9.47):

|ψB|2=4|C|2sen2(12kPa)e2q(xa)    ε2,ε4,ε6.(9.73)

9.4.4 Um exemplo numérico

Depois de tantas deduções, as tabelas abaixo nos ajudam a ter um visão geral do poço retangular finito.

DENSIDADE DE PROBABILIDADE PARA ε1,ε3,ε5
Barreira (x<0) Poço (0xa) Barreira (x>0)
4|C|2cos2(12kPa)e2qx 4|C|2cos2(kPx12kPa) 4|C|2cos2(12kPa)e2q(xa)
DENSIDADE DE PROBABILIDADE PARA ε2,ε4,ε6
Barreira (x<0) Poço (0xa) Barreira (x>0)
4|C|2sen2(12kPa)e2qx 4|C|2sen2(kPx12kPa) 4|C|2sen2(12kPa)e2q(xa)
NÚMEROS DE ONDA
Poço Barreiras
kP(εn)=2mεn/ q(εn)=2m(V0εn)/

Agora vomos considerar o mesmo poço retangular finito utilizado nas [ seções 9.2 e 9.3.4 ], aquele de largura a=100 Å e profundidade V0=100 meV. Já calculamos numericamente os valores dos níveis de energia.

NÍVEIS DE ENERGIA
Valores (meV)
ε1=3,0
ε2=11,9
ε3=26,6
ε4=46,8
ε5=71,6
ε6=97,9

As Figs. 9.6, 9.7 e 9.8, apresentam os contornos das densidades de probabilidade referentes aos estados n=1, n=2 e n=6. Na construção dessas figuras, o termo 4|C|2 foi desconsiderado.

Densidades de probabilidade para $n=1$.

Figura 9.6: Densidades de probabilidade para n=1.

Densidades de probabilidade para $n=2$.

Figura 9.7: Densidades de probabilidade para n=2.

Densidades de probabilidade para $n=6$.

Figura 9.8: Densidades de probabilidade para n=6.

Alguém poderá estar se perguntando como fazer para normalizar as curvas acima, pois o termo 4|C|2 foi desprezado. Um dos caminhos é tentar determinar o valor desse termo por meio de integração numérica. Sabemos que o conceito de normalização implica em: 4|C|2(área da curva não normalizada)=1.

Há vários softwares que podem determinar a área debaixo da curva não normalizada. No caso em questão, a integração feita com certo software livre resultou nos seguintes valores:

4|C|2=1,78×108 m1    n=1,4|C|2=1,12×108 m1    n=6.

Analisando as figuras acima, o que chama atenção é a penetração da densidade de probabilidade nas barreiras, comportamento puramente quântico. E a penetração cresce à medida que subimos de estado. Na interface x=100 Å, quando n=1, a densidade de probabilidade marca 3%, mas, quando n=6, sobe para 98% nesse caso, há grande chance da partícula ser encontrado dentro da barreira. É claro, a penetração está vinculada à energia da partícula. Uma partícula no estado n=6 tem muito mais energia que no estado n=1.

As barreias possuem comprimento infinito. Não há propagação de onda dentro delas. Ali, a função de onda se transforma em exponencial decrescente. Portanto, não há fluxo de probabilidade dentro das barreiras. Por não haver transmissão, as barreiras são 100% reflexivas. Isso significa que toda partícula que “bate” nas paredes do poço, “volta.” O fluxo se estabelece somente dentro do poço. Mesmo assim, há probabilidade de se encontrar uma partícula dentro das barreiras. Não há fluxo de partículas, mas uma experiência para “colher” uma partícula dentro da barreira pode registrar a presença dela.