Introdução à Mecânica Quântica

Jaime E. Villate.
Universidade do Porto, Portugal, 2026.

☰

5. Oscilador harmónico quântico

5.1. Hamiltoniano do oscilador harmónico

Na mecânica clássica, a função hamiltoniana de um objeto de massa m que oscila no eixo x devido a uma força elástica com constante elástica k, é:

H(x,p)=12mp2+k2x2
(5.1)

O respetivo operador hamiltoniano será então:

H^=12mp^2+k2x^2
(5.2)

onde x^ e p^ são os operadores de posição e quantidade de movimento.

Os valores próprios desse operador são os possíveis valores da energia mecânica do sistema, E. Para determinar esses possíveis valores, teremos de resolver o problema de valores próprios do hamiltoniano:

H^|E⟩=E|E⟩
(5.3)

onde |E⟩ é um vetor próprio, normalizado.

Aparentemente o hamiltoniano depende de dois parámetros, m e k, mas na realidade é necessário apenas um parámetro ω. Isso pode ser demonstrado introduzindo novos operadores de posição e quantidade de movimento, definidos da forma seguinte:

X^=(mk)1/4x^P^=1(mk)1/4p^
(5.4)

e a expressão do hamiltoniano em função desses novos operadores é,

H^=ω2(P^2+X^2)
(5.5)

onde a frequência angular ω é definida por:

ω=km
(5.6)

Observe-se que ω tem unidades de frequência (inverso do tempo) e os operadores P^2 e X^2 têm ambos unidades de energia vezes tempo. Os valores próprios de P^2 e X^2 são os valores próprios de P^ e X^, ao quadrado, que são números positivos ou zero. Como tal, os valores próprios do hamiltoniano não poderão ser valores negativos: E≥0.

Usando a definição (5.4) dos operadores X^ e P^ o seu comutador é:

[X^,P^]=X^P^−P^X^=x^p^−p^x^=[x^,p^]
(5.7)

e, como vimos no capítulo anterior, o comutador dos operadores de posição e quantidade de movimento é igual iℏ:

[X^,P^]=iℏ
(5.8)

5.2. Operadores de escada

A forma do hamiltoniano (5.5) faz lembrar o produto entre um número complexo e o seu conjugado; neste caso:

(X^−iP^)(X^+iP^)=X^2+P^2+i(X^P^−P^X^)
(5.9)

o ultimo termo, entre parênteses, é o comutador de X^ e P^, que é igual a iℏ. Como tal,

(X^−iP^)(X^+iP^)=X^2+P^2−ℏ
(5.10)

Vamos definir um operador, designado por operador de abaixamento (o fator ℏ no denominador garante que seja um operador sem unidades):

a^=12ℏ(X^+iP^)
(5.11)

e o seu operador adjunto, designado por operador de levantamento, é:

a^†=12ℏ(X^−iP^)
(5.12)

Conjuntamente, estes dois operadores designam-se operadores de escada e, em função deles, a expressão (5.10) escreve-se:

X^2+P^2−ℏ=2ℏa^†a^
(5.13)

E o hamiltoniano (5.5) é então:

H^=ℏω(a^†a^+12)
(5.14)

Os operadores de escada não são hermíticos, mas o operador a^†a^ sim é hermítico, e sem unidades.

O comutador dos operadores de escada é igual a 1, como pode mostrar-se substituindo as suas definições e usando o valor do comutador dos operadores de posição e quantidade de movimento:

[a^,a^†]=12ℏ[X^+iP^,X^−iP^]=−iℏ[X^,P^]=1
(5.15)

E o comutador do hamiltoniano (5.14) com os operadores de escada é:

[H^,a^]=ℏω[a^†a^+12,a^]=ℏω(a^†[a^,a^]+[a^†,a^]a^)=−ℏωa^
(5.16)
[H^,a^†]=ℏω[a^†a^+12,a^†]=ℏω(a^†[a^,a^†]+[a^†,a^†]a^)=ℏωa^†
(5.17)

5.3. Espetro de energia do oscilador harmónico

Se |E⟩ for vetor próprio do hamiltoniano, correspondente ao valor próprio E de energia,

H^|E⟩=E|E⟩
(5.18)

Consideremos agora o ket a^|E⟩. O hamiltoniano atuando nesse ket é igual a:

H^(a^|E⟩)=([H^,a^]+a^H^)|E⟩
(5.19)

e, usando a expressão (5.16) obtém-se:

H^(a^|E⟩)=(a^H^−ℏωa^)|E⟩=(E−ℏω)a^|E⟩
(5.20)

Isto mostra que a^|E⟩ é também vetor próprio do hamiltoniano, mas com energia E−ℏω. A norma do ket a^|E⟩ é:

‖a^|E⟩‖=⟨E|a^†a^|E⟩=⟨E|(H^ℏω−12)|E⟩=Eℏω−12
(5.21)

Aplicando sucessivamente o operador de abaixamento ao ket a^|E⟩, conclui-se que E−2ℏω, E−3ℏω, … também são valores próprios da energia. Como esses valores próprios da energia não podem ser negativos, deverá existir um valor mínimo da energia, Eo, que quando aplicarmos no respetivo vetor próprio o operador de abaixamento, o resultado seja o ket nulo, com norma igual a zero. Igualando a zero a expressão (5.21) do quadrado da norma de a^|E0⟩, obtemos o valor mínimo da energia:

E0=ℏω2
(5.22)

O hamiltoniano atuando no ket a^†|E0⟩ dá:

H^(a^†|E0⟩)=([H^,a^†]+a^†H^)|E0⟩
(5.23)

e, usando a expressão (5.17) obtém-se:

H^(a^†|E0⟩)=(a^†H^+ℏωa^†)|E0⟩=3ℏω2a^†|E0⟩
(5.24)

Isto mostra que a^†|E0⟩ é também vetor próprio do hamiltoniano, mas com energia 3ℏω/2. Aplicando sucessivamente o operador de elevamento, a^†, obtemos que 5ℏω/2, 7ℏω/2, … também são valores próprios. Como tal, os valores próprios do hamiltoniano do oscilador harmónico simples são:

En=ℏω(n+12)(n=0,1,2,…)
(5.25)

Aplicando o operador de abaixamento num dos estados próprios de energia, a^|En⟩, o resultado é o estado |En−1⟩, multiplicado por uma constante; constante essa que é a norma de a^|En⟩, que é obtida substituindo E por En=ℏω(n+1/2) na equação (5.21). Temos assim que,

a^|En⟩=n|En−1⟩
(5.26)

Já |Φn⟩=a^†|En⟩ é proporcional a |En+1⟩, e como,

⟨En|a^a^†|En⟩=⟨En|(H^ℏω+12)|En⟩=Enℏω+12=n+1
(5.27)

O efeito do operador de levantamento num estado próprio da energia é:

a^†|En⟩=n+1|En+1⟩
(5.28)

Exercícios

5.1 Encontre as expressões dos operadores de posição, x^, e quantidade de movimento, p^, em função de a^ e a^†.

Resolução. substituindo a definição (5.4) dos operadores X^ e P^, na expressão (5.11) do operador de abaixamento, e na expressão (5.12), e tendo em conta que k/m=ω, obtemos:

a^=mω2ℏ(x^+imωp^)a^†=mω2ℏ(x^−imωp^)

A soma e a diferença desses operadores são,

a^+a^†=2mωℏx^a^−a^†=i2mℏωp^

E as expressões dos operadores de posição e quantidade de movimento são:

x^=ℏ2mω(a^+a^†)p^=−imℏω2(a^−a^†)

5.2 Se o estado do sistema for um dos estados próprios do hamiltoniano do oscilador harmónico, |En⟩, de energia En, determine os valores esperados de: ⟨x^⟩, ⟨p^⟩, ⟨x^2⟩ e ⟨p^2⟩, e mostre que o valor esperado do hamiltoniano,

⟨H^⟩=12m⟨p^2⟩+mω22⟨x^2⟩

é igual à energia En.

Resolução. Usando a expressão de x^ obtida no exercício 5.1, e as equações (5.26) e (5.28),

x^|En⟩=ℏ2mω(a^+a^†)
=ℏ2mω(n|En−1⟩+n+1|En+1⟩)

e como os kets |En−1⟩ e |En+1⟩ são ortogonais a |En⟩, o valor esperado e x^ é zero:

⟨x^⟩=⟨En|x^|En⟩=0

De forma análoga, a partir da expressão de p^ obtida no exercício 5.1, e as equações (5.26) e (5.28), obtemos:

p^|En⟩=−imℏω2(a^−a^†)
=−imℏω2(n|En−1⟩−n+1|En+1⟩)

e o valor esperado e p^ também é zero:

⟨p^⟩=⟨En|p^|En⟩=0

A expressão de x^2 é:

x^2=ℏ2mω(a^+a^†)2=ℏ2mω(a^2+a^†2+a^a^†+a^†a^)

e usando as equações (5.26) e (5.28), obtemos:

x^2|En⟩=ℏ2mω(n(n−1)|En−2⟩+(n+1)(n+2)|En+2⟩
+(n+1)|En⟩+n|En⟩)

O produto de ⟨En| com os dois primeiros kets é zero, e com os dois últimos é 1. Como tal, o valor esperado do operador x^2 é:

⟨x^2⟩=⟨En|x^2|En⟩=ℏ(2n+1)2mω

A expressão de p^2 é:

p^2=−mℏω2(a^−a^†)2=−mℏω2(a^2+a^†2−a^a^†−a^†a^)

e usando as equações (5.26) e (5.28), obtemos:

p^2|En⟩=−mℏω2(n(n−1)|En−2⟩+(n+1)(n+2)|En+2⟩
−(n+1)|En⟩−n|En⟩)

O produto de ⟨En| com os dois primeiros kets é zero, e com os dois últimos é 1. Como tal, o valor esperado do operador p^2 é:

⟨p^2⟩=⟨En|x^2|En⟩=mℏω(2n+1)2

Com os valores obtidos para ⟨x^2⟩ e ⟨p^2⟩, temos que:

12m⟨p^2⟩+mω22⟨x^2⟩=ℏω(2n+1)4+ℏω(2n+1)4=ℏω(n+12)

que é exatamente o valor de En.