ads:
Universidade de São Paulo
Instituto de Física de São Carlos
Introdução a Química Quântica
τ
τ
τ
dH
d
dH
E
∫
∫
∫
ΨΨ=
ΨΨ
ΨΨ
=
*
*
*
André Luis Bonfim Bathista e Silva
IFSC/USP 2003
ads:
Livros Grátis
http://www.livrosgratis.com.br
Milhares de livros grátis para download.
Autor: André Luis Bonfim Bathista e Silva – IFSC/USP 2003
Cap. 1: Operadores em Química Quântica
Cap. 2: Equação de Schrödinger e Suas Aplicações
Cap. 3: Método de Hückel – Método Semi-Empírico
Cap. 4: Método Variacional
Cap. 5: Teoria do Campo Autoconsistente (SCF)
Cap. 6: Método de Hartree-Fock-Roothaan
Cap. 7: Teoria da Perturbação de Rayleigh-Schrödinger.
ads:
1 Edição
André Luis Bonfim Bathista e Silva
concluiu o Mestrado em
Ciência e Engenharia de Materiais no Instituto de Física de
São Carlos IFSC – USP em 2004, graduou em Licenciatura
Plena Em Física pela Universidade Federal de Mato Grosso
(UFMT) em 2001, conclui o curso técnico em Química
Industrial na Escola Técnica Federal de Mato Grosso (ETF-MT)
em 2000. Atualmente está cursando Doutorado em Ciência e
Engenharia de Materiais no Grupo de LEAR-IFSC-USP e o
curso Técnico em Eletrônica pela Escola Técnica Estadual –
Paulino Botelho. Publicou. Publicou 5 artigos em periódicos
especializados e 52 trabalhos em anais de eventos, sendo 6
em eventos Internacionais. Participou do desenvolvimento de
2 produtos tecnológicos. Recebeu 2 prêmio na Iniciação
Científica da Associação de Usuários de Ressonância
Magnética Nuclear (AUREMN) e uma Menção Honrosa da
Universidade Federal do Rio de Janeiro (UFRJ). Atua na área
de Física, com ênfase em Ressonância Magnética Nuclear. Em
suas atividades profissionais interagiu com 28 colaboradores
em co-autorias de trabalhos científicos.
Capítulo 1: Operadores em Química Quântica
1.1 Introdução
1.2 Conceitos Preliminares
1.3 O Formalismo Matemático da Mecânica Quântica
1.4 Operadores
1.5 Comutadores
1.6 Observáveis
1.7 Valor Esperado
1.8 Fluxo de Probabilidade
1.9 Postulados da Mecânica Quântica da Função de Onda
1.10 Bibliografia:
Capítulo 2: Equação de Schrödinger
2.1 Introdução
2.2 Construção da Equação de Schrödinger
2.3 Números Quânticos
2.4 Função de Onda na Concepção de Max Born
2.5 Valor Esperado da Função de Onda
2.6 Aplicação da Equação de Schrödinger
2.6.1 Poço de Potencial não relativístico
2.7 Bibliografia:
Capítulo 3: Método de Hückel.
3.1 Método Semi-Empírico
3.2 Elaboração
3.3 Considerações do Método de Hückel
3.4 Aplicação do Método de Hückel no Etileno
3.5 Aplicação do Método de Hückel no Butadieno
3.6 Aplicação do Método de Hückel para o Ciclobutadieno
3.7 Bibliografia:
Capítulo 4: Método Variacional.
4.1 Introdução
4.2 Extensão do método Variacional
4.3 Função Variacional Linear
4.4 Aplicação do Método Variacional
4.5 Método Variacional para o Átomo de Hidrogênio
4.6 Método Variacional para o Átomo de Hélio
4.7 Método Variacional para Sistema de dois Spins
4.8 Referências:
Capítulo 5: Teoria do Campo Autoconsistente (SCF)
5.1 Introdução
5.2 Valor médio da Energia
5.3 O Valor médio de H
1
5.4 O Valor médio de H
2
5.5 Cálculo da Energia Média
5.6 As Equações de Hatree-Fock.
5.7 Referência:
Capítulo 6: Método de Hartree-Fock-Roothaan.
6.1 Introdução
Capítulo 7: Teoria da Perturbação
7.1 Introdução
7.2 Correção de ordem um para a energia
7.3 Correção de ordem um para a função de onda
7.4 Tratamento de um sistema de spins pela Teoria da Perturbação
7.5 Referência:
Capítulo 1: Operadores em Química Quântica
1.1 Introdução
A Mecânica Quântica foi construída na década de 20 para lidar com o átomo, os
físicos consideravam-na como um instrumento provisório destinado ao fracasso fora dos
domínios do átomo, entretanto a mecânica quântica prosperou e ultrapassou os mais
fantásticos sonhos de seus inventores. Atualmente, devido a aceitação dos Postulados
Quânticos e suas extensivas aplicações em Centros Científicos, vem proporcionando maior
divulgação de experimento clássicos da Mecânica Quântica, bem como seus inventores. A
internet, onde é um meio de informação de mais alto padrão, nos favorece diversas fontes
de trabalhos diretos ou correlacionados à Mecânica Quântica, sendo estes: revistas,
fascículos e até mesmo jornais científicos.
1.2 Conceitos Preliminares
Para explicar a mecânica dos átomos e dos sistemas nucleares, foram desenvolvidas,
a partir de 1900, diversas teorias. Estas teorias foram se enquadrando a uma denominação
comum: “Mecânica Quântica”(MQ).
Assim, MQ é um conjunto de teoria dos sistemas atômicos e nucleares. Este conjunto
não é unitário, mas sim uma sucessão de diversas teorias, umas complementares das outras.
Estas teorias surgiram e desenvolveram-se da Física Clássica, particularmente da mecânica
newtoniana e da teoria eletromagnética de Maxwell. Demonstrando-se que os conceitos
clássico e eletromagnético eram suficientes para explicação do mundo macroscópico, mas,
incapazes para uma explicação adequada e coerente do mundo microscópico, i.e., atômico-
nuclear e também molecular, nasceram à chamada “Antiga teoria dos quanta” ou Mecânica
Quântica Antiga. Esta mecânica compreende as teorias de Planck (1900), Einstein (1905),
Bohr (1913) e De Broglie (1924).
A antiga MQ embora explicasse muitos fenômenos até então incompreendidos,
falhava em sua base lógica para outros fenômenos. E assim Heisenberg (1925),
desenvolveu a “Mecânica das Matrizes”, que não implica em nenhum modelo atômico, pois
é um edifício puramente matemático. Quase simultaneamente, Schröndiger (1926)
baseando-se nas idéias de De Broglie, desenvolveu a “Mecânica ondulatória”. Esta
mecânica tem um aspecto mais físico, embora algo impreciso. Esta imprecisão foi
desaparecendo com a aplicação das teorias de Pauli (1924) e Dirac (1926). As teorias de
Heisenberg até Dirac, são conhecidas como “Mecânica Quântica Moderna”.
Esquematicamente, temos:
MQ
MQA
MQM
teoria dos quantas de Planck
teoria atômica de Bohr
dualismo de De Broglie
mecânica das matrizes de Hiesenberg
mecânica ondulatória de Schöndiger
teoria do Spin de Pauli
teoria da energia positiva e negativa de
Dirac
Juntamente com estes físicos, representantes da MQ, figura outro grande número de
cientista que, experimentalmente, deram confirmação a cada um destas teorias, assim como
aplicaram teoricamente a muitos fenômenos Físicos-Químicos.
Primeiramente estaremos considerando a visão matemática, que a mecânica
quântica possui com uma linguagem acessível para estudantes de ciências e assim uma
breve discussão da MQA (Planck→ de Broglie) e MQM (Pauli→ Dirac). tratando
historicamente os fatos ocorridos e descrevendo alguns passos tomados pelos pesquisadores
da MQ, explanando o formalismo matemático.
1.3 O Formalismo Matemático da Mecânica Quântica
Na mecânica quântica as formas de se estabelecer relações numéricas e
representação física de eventos são de modelos não convencionais como o de costumeiro da
física e química. Às vezes tendo tratamento até filosófico de problemas.
1.4 Operadores
Em muitos processos matemáticos um elemento tem que se submeter a um conjunto
de operações quase sempre complicadas, mas que se repetem muitas vezes, a fim de aplicar
a outros elementos diferentes para correlacioná-los. O conjunto de operações e
transformações métricas são representadas por símbolos (Â, ∇
2
, , Ψ) que se denominam
operadores de transformação. Este conjunto de operações de transformação pode ser
representado também por uma matriz, a qual tem propriedades intrínsecas da teoria de
matrizes. Veja abaixo a representação de operadores e matrizes respectivamente.
hAf
ˆ
=
[1]
=
333231
232221
131211
ˆ
aaa
aaa
aaa
A
[2]
Estes operadores são denotados usando letras com acento circunflexo, como por
exemplo
Ô. Existe em princípio, um operador para cada informação que se deseja sobre o
sistema, sendo os mais importantes os operadores que fornecem a energia, chamados de
operadores Hamiltonianos que são denotados pelo símbolo
H. Os operadores, por serem
ordens ou instruções matemáticas, devem ser aplicadas pela direita nas funções desejadas
por exemplo a função de onda Ψ ou usando a chamada representação de Dirac
Ψ
. Se Ô =
d/dx, ou seja a ordem matemática é “derive com relação a x”, e
Ψ = f(x),podemos ver
claramente que
Ô
Ψ
(df(x)/dx) não é o mesmo que
Ψ
Ô , ou f(x)d/dx (4). Como em
alguns casos é necessário trabalhar com combinações lineares de funções de onda, os
operadores também devem ser lineares, i.e.:
()
2121
ˆˆˆ
Ψ+Ψ=Ψ+Ψ OOO
[3]
Quando aplicamos o operador na função de onda pode acontecer dois eventos. No
caso mais simples a operação sobre Ψ fornece a própria função Ψ multiplicada por uma
constante
a, ou seja:
Ψ=Ψ aO
ˆ
[4]
Neste caso a função de onda Ψ é dita ser autofunção do operador
Ô com autovalor
a. O autovalor obtido corresponde à propriedade que se deseja medir, de forma que, quando
se fizer uma medição desta propriedade em um conjunto de sistemas idênticos, sempre será
obtido o mesmo valor
a. Mas é óbvio que as medições realizadas usando os operadores
devem fornecer valores reais e não números complexos. Apresentando esta condição o
operador é dito Hermitiano, o que matematicamente corresponde a expressão:
(
)
*
ˆˆ
ijji
OO ΨΨ=ΨΨ
[5]
A equação [5] garante que se
Ψ=Ψ aO
ˆ
e que
***
ˆ
Ψ=Ψ aO
, então a = a
*
, o
que só acontece se
a for um número real, ou seja, se não contiver o número i na sua
expressão.
O segundo tipo de resultado que pode ser esperado quando se aplica um operador
em uma função de onda é que a função de onda seja alterada pelo operador, como mostrado
na equação 6:
21
ˆ
Ψ=Ψ bO
[6]
Agora a função não é autofunção do operador, e o valor
b corresponde a uma única
medição da probabilidade desejada em um só sistema dentro de um universo de sistemas
idênticos. Neste caso , cada vez que se fizer a medição experimental da propriedade seria
obtida uma resposta diferente, sendo portanto necessário calcular o valor médio das
medições, também chamado de valor esperado, <
b>, para ter um resultado confiável da
medição. O valor médio é calculado exatamente como se calculam os valores médios em
estatística:
ΨΨ
ΨΨ
=
O
b
ˆ
[7]
O principal problema agora é definir matematicamente os operadores e as funções
de onda. Os operadores são obtidos através dos comutadores dos operadores.
Se tivermos uma função de onda descrita por
)(
2
Etxp
h
i
x
Ae
−
=Ψ
π
[8]
que pode ser aplicada a descrever qualquer partícula de massa m e velocidade V
x
que
apresente comportamento ondulatório, sendo portanto uma equação geral e não só referente
ao fóton. Agora utilizando-se dum operador para derivá-la em relação à
x, e considerando
que a energia
E e o momento P
x
da partícula associada à onda são constantes, temos:
Ψ=
∂
Ψ
∂
h
ip
x
x
π
2
[9]
rearranjando temos:
Ψ=
∂
Ψ
∂
−
x
p
x
ih
π
2
[10]
analisando a equação [10] vemos que a função de onda Ψ é uma autofunção do operador
x
ih
∂
∂
−
π
2
com autovalor
p
x
, ou seja, que o operador extrai o valor de momento da função
de onda
x
ih
∂
∂
−
π
2
. Isto significa que, para obter operadores podemos escrever a expressão
clássica da quantidade que nos interessa, e substituir os termos referentes a componente em
x do momento,
p
x
, por
x
ih
∂
∂
−
π
2
. Por analogia, para as componentes em y e z a substituição
deve ser por
y
ih
∂
∂
−
π
2
e
z
ih
∂
∂
−
π
2
respectivamente. Já onde existem coordenadas (x,y,z) ou
tempo (t) não é necessário fazer nenhuma modificação para converter a expressão em um
operador. Esta metodologia pode ser expressa em forma de tabela. Veja Tabela 1
Tabela 1: substituições para a construção de operadores.
Termo na
expressão clássica
Termo no
operador
Símbolo do
operador
X X
X
ˆ
Y Y
Y
ˆ
Z Z
Z
ˆ
P
x
x
ih
∂
∂
−
π
2
x
p
ˆ
P
y
y
ih
∂
∂
−
π
2
y
p
ˆ
P
z
z
ih
∂
∂
−
π
2
z
p
ˆ
T t
t
ˆ
Outros exemplos que expressa estes operadores pode serem vistos nas Tabelas 2 e 3
Tabela 2: Exemplos do operador Hamiltoniano para o movimento de uma
partícula de massa m em diferentes campos de força definidos pela função
(operador) potencial V.
Operador
(a) Partícula livre
V=0
x
2
22
2
ˆ
dx
d
m
H
=
−=
(b) barreira de
potencial
V=V
0
x
0 a
2
22
2
ˆ
dx
d
m
H
=
−=
(x<0; x>0)
0
2
22
2
ˆ
V
dx
d
m
H +−=
=
(0<x<a)
(c) Oscilador
harmônico
V=(1/2)Kx
2
x
K
2
2
22
2
1
2
ˆ
Kx
dx
d
m
H +−=
=
Tabela 3: Exemplos do operador para átomos e moléculas.
Operador H
(a) Átomos de
um elétron
r
Ze
m
H
2
2
2
2
ˆ
−∇−=
=
(b) Átomos de
muitos elétrons
∑∑∑∑
>==
+−∇−=
iij
ij
n
j
j
n
i
r
e
r
Ze
m
H
2
1
2
1
2
2
2
ˆ
=
(c) moléculas
∑∑ ∑∑∑∑∑∑
>>==
++−∇−∇−=
iij
ij
n
ii
i
i
n
i
r
eZZ
r
e
r
eZ
mm
H
ααβ
αβ
αα
β
α
α
α
α
2
2
1
2
2
2
1
2
2
2
1
2
ˆ
==
1.5 Comutadores
Em 1926, Born, Heisenberg e Jordan apresentaram na Zeitschrift für Physik a
mecânica quântica, desenvolvida por Born em 1924, e Heisenberg, em 1925, sob o ponto de
vista matricial. Nesse trabalho, eles apresentaram as relações de comutações para o
momento angular
L
G
para um sistema de muitas partículas:
[
]
[]
[]
.
2
,
,
2
,
,
2
,
xzy
yxz
zyx
L
i
h
LL
L
i
h
LL
L
i
h
LL
π
π
π
=
=
=
[11]
Hoje em alguns cursos de MQ, podemos ver que esses comutadores podem ser visto
da forma apresentada acima.
Mas qual será a idéia de comutação ?, comutar significa ‘
Realizar diversas
funções
’. Imaginemos uma situação na qual temos dois operadores Hermitianos Ô e Ê que
comutam entre si, ou seja
ÔÊ – ÊÔ = 0 (o comutador também é representado como [Ô,Ê]
= 0
). Suponhamos ainda que existe um conjunto de funções de onda Ψ
i
que são autofunções
do operador
Ô, com autovalores diferentes entre a
i
, i.e., sem degenerescência
.
iii
aO Ψ=Ψ
ˆ
neste caso qual seria a relação do operador
Ê com este conjunto de funções?. Se operarmos
pela esquerda com
Ê teremos que
(
)
iiiii
ÊaaÊOÊ Ψ=Ψ=Ψ )(
ˆ
[12]
e pela condição de comutação:
)(
ii
ÊÔÊÔÔÊÊÔ Ψ=Ψ→=
[13]
de (12) e (13) temos que
)((
iii
ÊaÊÔ Ψ=Ψ
[14]
A equação [14] nos mostra explicitamente que
Ê
i
Ψ
é uma autofunção do operador
Ô, com o autovalor a
i
. Como Ψ
i
é a única função com autovalor a
i
, devido ‘a não
degenerescência dos sistema, a única forma da equação [14] ser verdadeira é quanto a
função
Ê
i
Ψ
seja um múltiplo de Ψ
i
, i. e.,
iii
dÊ Ψ=Ψ
[15]
Em outras palavras, a função
i
Ψ
também é autofunção do operador Ê. Este
resultado na forma de teoria seria: ‘Se dois operadores Hermitianos Ô e Ê comutam, então
existe um conjunto de funções que são autofunções dos dois operadores’. Este teorema é
extremamente importante pois mostra que as únicas propriedades de um sistema que
podemos medir simultaneamente e com precisão são aquelas para as quais a função de onda
representando o estado do sistema é a autofunção dos seus respectivos operadores, ou seja,
para medir precisa e simultaneamente duas ou mais propriedades de um sistema, seus
Degenerescência: significa os mesmos autovalores de energia.
respectivos operadores devem comutar. De fato, este resultado é uma forma de expressão
do Princípio de Incerteza de Heisenberg, pois indica que só é possível conhecer
simultaneamente com precisão as propriedades de um sistema cujos operadores comutam.
Claramente, de acordo com a formulação mais conhecida deste Princípio, os operadores
referentes à posição e ao momento de um elétron em um átomo não comutam.
1.6 Observáveis
As propriedades (E, P, X
i
) que formam um observável têm duas características
importantes:
1) Podem excluir-se mutuamente, ou 2) Constituem um conjunto completo, i.e., incluem
entre si todos os possíveis resultados da medição do observável (E, P, X
i
).
Em MQ, pode-se distinguir dois tipos importantes de relações entre observáveis.
Analisemos dois observáveis, q e r são compatíveis se as propriedades (E, P, X
i
) de todos
os pares (q
i
, r
j
) – q
i
representa qualquer propriedade do observável q, r
j
qualquer
propriedade de r – forem compatíveis ou de mútua exclusão. Sendo o último caso, o par de
propriedades que não correspondem a qualquer situação fisicamente real ou previsível.
No entanto, para que dois observáveis q e s sejam incompatíveis basta que, entre
todos os pares de propriedades (q
i
, s
j
) exista um que seja incompatível.
Um exemplo de incompatibilidade é o momento linear (P) e a posição (x
i
). A energia de
uma partícula livre (energia cinética) e a sua quantidade de movimento são também
observáveis compatíveis. Neste caso, os pares de propriedades formados por uma
propriedade de cada um destes observáveis ou são compatíveis ou excluem mutuamente; a
exclusão múua aplica-se a situações fisicamente absurdas, por exemplo quando, uma das
prorpiedades for nula e a outra diferente de zero porque a energia cinética e a quantidade de
movimento estão relacionadas entre si através de uma dependência funcional, K = p
2
/(2m).
1.7 Valor Esperado
Se considerarmos, uma partícula e onda associada, a função
Ψ
(r,t)
e se essa função
não se anula num intervalo entre r e r + dr, na medida de sua posição há uma probabilidade
finita dessa partícula ser encontrada. Não podemos, atribuir a coordenada um valor bem
definido, no entanto, é possível especificarmos uma posição média da partícula.
Imaginemos a medida da posição da partícula no instante t, a probabilidade de
encontrá-la entre r e r + dr é dada pela equação.
),(
),(
),( tr
tr
tr
P
ψψ
∗
=
[16]
onde P
(r,t)
é a probabilidade de encontramos a partícula. Repetindo essa experiência a uma
certa freqüência no mesmo instante e registrando os valores de P
(r,t)
, podemos usar a média
dos valores observados para caracterizar a posição da partícula no instante t.
Este valor é representado por <r>, valor esperado da coordenada r. Veja abaixo
como podemos demonstrar matematicamente este cálculo.
∫
+∞
∞−
>=< drrPr
tr ),(
[17]
como:
),(
),(
),( tr
tr
tr
P
ψψ
∗
=
,
substituímos [16] em [ 17] temos que:
∫
+∞
∞−
∗
>=< drrr
tr
tr
),(
),(
ψψ
[18]
1.8 Fluxo de Probabilidade
A probabilidade de encontrar uma partícula a uma certa posição
x e num tempo t é
dado por
),(
),(
),( tr
tr
tr
P
ψψ
∗
=
Sendo o fluxo qualquer grandeza que varia o tempo.
][),(
),(
),(
tx
tx
t
tx
t
P
∗
∂
∂
=
∂
∂
ψψ
[19]
),(),(),(
),(
),(
tx
t
tx
t
tx
t
P
tx
tx
∂
∗
∂
+
∂
∂
=
∂
∂
∗
ψ
ψ
ψ
ψ
[20]
como:
t
iV
x
m
txx
∂
∂
=+
∂
∂
−
ψ
ψ
=
=
),()(
2
22
2
consideramos que ela esteja livre de V(x).
Parte real:
t
x
m
i
∂
∂
=
∂
∂
ψψ
2
2
2
=
[21]
Parte imaginária:
t
x
m
i
∂
∂
=
∂
∂
∗∗
ψψ
2
2
2
=
[22]
Substituindo [21] e [22] em [20]
∂
∂
−+
∂
∂
=
∂
∂
∗
2
22
),(
2
2
),(
22
),(
xm
i
xm
i
tx
t
P
tx
tx
ψ
ψ
ψ
ψ
==
∂
∂
−
∂
∂
∂
∂
−=
∂
∂
∗
∗
xxmix
tx
t
P
tx
tx
ψ
ψ
ψ
ψ
),(
),(
2
),(
=
∂
∂
−
∂
∂
=
∗
∗
ψ
ψψ
ψ
xxmi
S
tx
2
),(
=
[23]
Sendo S o fluxo de probabilidade
),(),( tx
x
S
tx
t
P
∂
∂
−=
∂
∂
ou
0),(),( =
∂
∂
+
∂
∂
tx
x
S
tx
t
P
[24]
Sendo esta a equação da continuidade.
1.9 Postulados da Mecânica Quântica da Função de Onda
I) A todo sistema corresponde uma função de onda Ψ (x,t).
II) A seguinte correspondência existe entre as variáveis dinâmicas e os seguintes:
Variáveis Operadores
x, y, z, t x, y, z, t
F(x,t) F(x,t)
P
x
, P
y
, ou P
z
zi
ou
yixi ∂
∂
∂
∂
∂
∂
===
,
E
ti
∂
∂
−
=
III) Ψ (x,t) (onde x representa x, y, z) e suas derivadas primeiras, devem ser finitas e
unívocas.
IV) Ψ (x,t) deve ser normalizada
, isto é.
dxdydzdv
dv
tx
tx
=
=ΨΨ
∫∫∫
+∞
∞−
+∞
∞−
+∞
∞−
1
),(
),(
*
ou simplesmente
quando tem-se uma função = 1, a probabilidade é 100% de encontrar a particula numa região, isso significa
que esta função está normalizada.
1
*
=ΨΨ
∫
+∞
∞−
dv
ou mais simplesmente
1>=ΨΨ<
V) A cada observável corresponde um operador Â. O valor médio <Â> ou valor
esperado deste operador e definido como.
∫∫∫
+∞
∞−
+∞
∞−
+∞
∞−
ΨΨ>=< dvAA
ˆˆ
*
∫
ΨΨ>=< dvAA
*
ˆ
ou
>ΨΨ>=<< AA
*
ˆ
porém, se Ψ não for normalizada, então <A> é definido como
>ΨΨ<
>ΨΨ<
>=<
A
A
*
1.10 Bibliografia:
(1) Coelho, J. V. Mecânica Quântica (
2000), In: 100 Anos de Mecânica Quântica,
Departamento de Física, UFMT.
(2) Rey, A. B. Mecânica Quântica e Ondulatória, In: Fisica/Química Modernas. V. 3, SP,
Ed. Fortaleza,
1970
(3) Bathista, A. L. B. S., Nogueira, J. S. Uma Breve Discussão da Mecânica Quântica
,Cuiabá, UFMT. 2001.
(4) Almeida, W. B., Santos, Química Nova,
(5) Bassalo, J. M. F. Nascimentos da Física (1901-1950), Belém: EDUFPA,
2000. 503 p.
(6) Teixeira Dias, J. J. C. Química Quântica fundamentos e métodos. 1982. Ed. Fundação
Calouste gulbenkian. 448 p.
(7) Millomen, W. C. C. Notas de Aula, 2001. Departamento de Física – Universidade
Federal de Mato Grosso.
(8) Eisberg, R. Resnick, R. Física Quântica, Ed. Campus, Rio de Janeiro,
1979. 15
ª
Edição.
(9) Bathista, A. L. B. S., Nogueira, J. S. Da Physis à Física, O Livre Contexto Sobre a
essência da matéria: O Átomo. IX Encontro de Iniciação Científica, Cuiabá, UFMT.
2001.
(10) Born, M., Física Atômica,
1962, 4
ª
Edição. Ed: Fundação Calouste Gulbenkian,
Lisboa.
(11) Coelho, J. V. Mecânica Quântica, Cuiabá: EdUFMT, 215 p.1993.
(12) Millomen, W. C. Born e a Teoria Quântica (
2000), In: 100 anos de Mecânica
Quântica, Departamento de Física, UFMT
(17) De Broglie, L. Recherches sur la téorie des quanta. 1963, Paris: Masson et C
ie
editeurs.
128p. Réédition du text de
1924
(18) Landau, L. Lifshitz, E. Mecânica Quântica, Teoria não relativista,
1985. Mir, Moscou.
Vol.1.
Capítulo 2: Equação de Schrödinger.
2.1 Introdução
Em 1926, o físico austríaco Erwin Schrödinger (1887-1961) publicou quatro
trabalhos nos Annales de Physique Leipzig nos quais desenvolveu a sua famosa Mecânica
Quântica Ondulatória, cujo resultado principal é a equação para as órbitas estacionárias dos
elétrons atômicos, a igualmente famosa equação de Schrödinger:
() ()
[]
()
0
2
,,,,
2
,,
2
=Ψ−+Ψ∇
zyxzyxzyx
VE
m
=
[1]
Em relação ao trabalho de Bohr, o trabalho de Schrödinger foi bem mais completo.
Uma vez que prevê também o seguinte:
•
As autofunções são correspondentes a cada autovalor.
•
Prevê o cálculo da probabilidade de um determinado estado.
•
Prevê o cálculo da probabilidade de transição de um estado para outro.
•
Calcula os momentos angulares orbitais.
A equação de Schrödinger nada mais é que uma equação diferencial de segunda
ordem,a qual podemos aplicar para um sistema como o átomo de
1
H e calcularmos os seus
níveis de energias correspondentes. Historicamente foi o primeiro sistema que Schrödinger
tratou, onde os autovalores de energia são os mesmos que previstos por Bohr.
2.2 Construção da Equação de Schrödinger
Primeiramente temos que resolver um problema de dois corpos. Neste caso,
podemos reduzir o sistema de dois corpos a um sistema de um corpo, considerando a
massa reduzida do sistema:
21
21
mm
mm
+
=
µ
[2]
Este é o termo que é introduzido na equação Schrödinger e podemos adquiri-lo através do
cálculo do centro de massa e para acharmos o centro de massa temos que fazer a seguinte
igualdade m
1
r
1
=m
2
r
2
e se variarmos r logo teremos o centro de massa. Considerando
m
2
>>m
1
2211
rmrm
=
[3]
21
rrr
+
=
[4]
21
rrr
−
=
;
12
rrr
−
=
[4.1]
Substituindo [4.1] em [3] obtemos,
r
mm
m
r
12
2
1
+
=
e r
mm
m
r
21
1
2
+
=
[5]
Como o próton e o elétron estão translacionando e girando com velocidades próprias.
Podemos obter a expressão da energia cinética total do sistema e a velocidade angular.
22
2
22
2
11
vmvm
E
K
+=
[6]
1
1
r
v
=
ω
e
2
2
r
v
=
ω
(
)
22
22
22
11
2
1
ωω
rmrmE
K
+= [7]
(
)
22
22
2
11
2
1
ω
rmrmE
K
+=
sendo
(
)
2
22
2
11
2
1
rmrmI += [8]
Substituindo [5] em [8]
2
21
21
2
1
r
mm
mm
I
+
=
[9]
21
21
mm
mm
+
=
µ
[10]
22
2
1
ωµ
rE
K
=
[11]
No átomo de
1
H o núcleo é massivo e o elétron tem massa reduzida
µ
dada pela
equação [10] de tal maneira que gira em torno do núcleo estacionário.
elétron
Figura 1: representação do átomo de hidrogênio
O potencial de interação (coulombiano) do elétron-próton é dado pela relação
()
2/1222
0
2
,,
)(4
.
zyx
eZ
V
zyx
++
−=
πε
[12]
Onde
e
-
= carga do elétron, Z = carga do núcleo ( para o
1
H, Z = 1)
A eq. de Schroedinger dependente do tempo para este sistema é a seguinte relação,
conforme mostrado abaixo:
t
iV
zyx
tzyx
tzyxzyxtzyx
∂
Ψ∂
=Ψ+Ψ
∂
+
∂
+
∂
−
),,,(
),,,(),,(),,,(
2
2
2
2
2
22
2
=
=
µ
[13]
ou numa notação mais compacta: onde se consideramos o operador laplaciano
2
2
2
2
2
2
2
zyx
∂
+
∂
+
∂
=∇
[14]
Assim a eq. de Schroedinger pode ser reescrita assim:
2
2
2
t
iV
∂
Ψ∂
=Ψ+Ψ∇ =
=
µ
[15]
Cuja solução é dada:
e
-iEt/
z)y,(x,t)z,y,(x,
=
ϕ
=Ψ
[16]
Onde ψ
(x,y,z)
é a função onda independente do tempo dada pela resolução da eq. de
Schroedinger independente do tempo a seguir:
),,(),,(),,(),,(
2
2
2
zyxzyxzyxzyx
EV Ψ=Ψ+Ψ∇−
µ
=
[17]
ou
Ψ=Ψ+Ψ
∂
+
∂
+
∂
− EV
zyx
zyxzyxzyx ),,,(),,(),,,(
2
2
2
2
2
22
2
µ
=
(x,y,z)
[18]
Para resolvermos a equação acima é melhor escreve-la em coordenadas esféricas e usar o
método de separação de variáveis
elétron
φ
θ
Próton
y
x
z
Figura 2: representação da interação coulômbiana entre o
próton e o elétron, aplicando as notações de coordenadas
esféricas para resolver a equação de Schrödinger.
De acordo com a Figura 2, r é o raio vetor posição do elétron,
222
zyxr ++=
[19]
θ é o angulo polar, o qual cresce a partir do eixo z para o plano xy
222
yx
arccos
z
z
++
=
θ
[20]
e φ é o ângulo azimutal, o qual cresce de x para y.
x
y
arctg=
ϕ
[21]
onde temos
=
=
=
θ
φθ
φθ
rcosz
senrseny
cossenrx
E o laplaciano em coordenadas esféricas é dado assim:
2
2
222
2
2
2
sen
1
)(sen
sen
1
)(
1
ϕθ
θ
θ
θ
θ
∂
∂
+
∂
∂
∂
∂
+
∂
∂
∂
∂
=∇
rr
r
r
r
r
[22]
A equação de Schroedinger em coordenadas esféricas fica assim após aplicarmos o
operador laplaciano esférico: eq. [18]
),,(),,(),,(
2
2
222
2
2
2
sen
1
sen
sen
11
2
φθφθϕθ
ϕθ
θ
θ
θ
θ
µ
rrr
EV
rr
r
r
r
r
Ψ=Ψ+
∂
Ψ∂
+
∂
Ψ∂
∂
∂
+
∂
Ψ∂
∂
∂
−
=
[23]
Nota importante: Podemos escrever a função de onda
),,(
φθ
r
Ψ
Como um produto de três
funções tal que
(
)
)()()(,,
θ
φ
θ
φ
Θ
Φ
=
Ψ
rRr
1.
A primeira R(r) = dependência em r ( radial afastamento do elétron)
2.
A segunda Θ(θ) = dependência em θ (mostra a posição polar do elétron)
3.
A outra função Φ(ϕ) = mede a posição azimutal do elétron na sua trajetória
Podemos realizar esta representação porque há um campo de força central na equação.
Aplicando agora a separação de variáveis. A equação [23] pode ser reescrita (veja notação
da derivada e o que deriva parte constante fica fora da derivada)
ΘΦ=ΘΦ+
ΦΘ
+
ΘΦ
+
ΘΦ
− ER
sen
1
sen
sen
11
2
),,(
2
2
222
2
2
2
RV
d
d
r
R
d
d
d
d
r
R
dr
dR
r
dr
d
r
r
ϕθ
ϕθ
θ
θ
θ
θ
µ
=
ΘΦ=ΘΦ+
ΦΘ
+
ΘΦ
+
ΘΦ
−
ER )(
sen
1
sen
sen
11
2
2
2
222
2
2
2
RrV
d
d
r
R
d
d
d
d
r
R
dr
dR
r
dr
d
r
ϕθ
θ
θ
θ
θ
µ
=
[24]
Substituindo o potencial V(r )=-
ΘΦRpor divindo e
4
0
2
r
e
πε
a equação [23] fica:
2
2
2
0
2
22
2
2
2
2
1
4
sen
2
sen
sensen
m
d
d
r
e
Er
d
d
d
d
dr
dR
r
dr
d
R
−=
Φ
Φ
−=
++
Θ
Θ
+
ϕ
πε
θ
µ
θ
θ
θ
θθ
=
[25]
)(constante
1
2
2
2
m
d
d
−=
Φ
Φ
−
ϕ
2
0
2
22
2
2
2
2
4
sen
2
sen
sensen
m
r
e
Er
d
d
d
d
dr
dR
r
dr
d
R
−=
++
Θ
Θ
+
πε
θ
µ
θ
θ
θ
θθ
=
[26]
Multiplicando por r
2
sen
2
θ todos os membros fica
0
4
sen
2
1
sen
sen1sen
0
2
22
22
22
2
2
=
++
Φ
Φ
+
Θ
Θ
+
r
e
Er
d
d
d
d
d
d
dr
dR
r
dr
d
R
πε
θ
µ
ϕ
θ
θ
θ
θθ
=
Ou
0
4
sen
2
sen
sen
sen
sen
sen
sen
0
2
22
22
2
22
22
2
22
2
2
22
=
++
Φ
Φ
+
Θ
Θ
+
r
e
Er
d
d
r
r
d
d
d
d
r
r
dr
dR
r
dr
d
rR
r
πε
θ
µ
ϕθ
θ
θ
θ
θ
θ
θθ
=
Isolando os termos em
ϕ (passando ao segundo membro)
Note : do lado direito a equação depende só de
ϕ e do lado esquerdo a dependência é em θ
e r. Como temos uma igualdade a constante de separação (-m
2
) deve ser a mesma para
ambas as equações.
2
2
2
0
2
22
2
2
2
2
1
4
sen
2
sen
sensen
m
d
d
r
e
Er
d
d
d
d
dr
dR
r
dr
d
R
−=
Φ
Φ
−=
++
Θ
Θ
+
ϕ
πε
θ
µ
θ
θ
θ
θθ
=
)(constante
1
2
2
2
m
d
d
−=
Φ
Φ
−
ϕ
2
0
2
22
2
2
2
2
4
sen
2
sen
sensen
m
r
e
Er
d
d
d
d
dr
dR
r
dr
d
R
−=
++
Θ
Θ
+
πε
θ
µ
θ
θ
θ
θθ
=
[27]
imraizes cujas 0
d
ou
1
2
2
2
2
2
2
±==+
Φ
−=
Φ
Φ
λ
ϕϕ
m
d
m
d
d
logo a solução é
Ne)(
im
ϕ
ϕ
±
=Φ
[28]
que é a solução particular, sendo
m número quântico magnético. Φ(ϕ) é a função onda
átomo de
1
H que contém toda a dependência em função do Potencial ϕ. Onde ϕ varia de 0 a
2
π. Normalizando a função Φ(ϕ)
1=ΦΦ
∫
==ΦΦ
π
ϕϕ
ϕ
2
0
im-
1Ne * deN
im
∫
==ΦΦ
π
ϕϕ
ϕ
2
0
im-2
1e deN
im
[29]
2
2
0
2
22
2
2
2
2
1
4
sen
2
sen
sen1sen
ϕ
πε
θ
µ
θ
θ
θ
θθ
d
d
r
e
Er
d
d
d
d
dr
dR
r
dr
d
R
Φ
Φ
−=
++
Θ
Θ
+
=
ou
π
π
2
1
N 1]2[
2
=⇒=N
logo a função vale
ϕ
π
ϕ
im
e
2
1
)(
=Φ
[30]
como
única é )2()( ela⇒+Φ=Φ
π
ϕ
ϕ
mimeeee
imimimim
ππ
ππϕϕ
2sen2cos11
22
+=⇒=⇒=
magnético) quântico número ( inteirom 3...210 m se satisfaz
=
±
±±=queveja
escolhe a solução
λ = im (movimento de giro elétron sentido horário).
Assim temos até o momento
)()(
2
1
),,(
ϕ
π
ϕθ
ϕ
Φ=Ψ rRer
im
[31]
Tomando a equação [27],
2
0
2
22
2
2
2
2
4
sen
2
sen
sensen
m
r
e
Er
d
d
d
d
dr
dR
r
dr
d
R
−=
++
Θ
Θ
+
πε
θ
µ
θ
θ
θ
θθ
=
Dividindo por
θ
2
sen
fica
sen
4
2
sen
11
2
2
0
2
2
2
2
θ
πε
µ
θ
θ
θ
m
r
e
Er
d
d
d
d
dr
dR
r
dr
d
R
−=
++
Θ
Θ
+
=
[32]
Separando os termos da R(r ) da função
Θ(θ) fica
sen
1
sen
4
2
1
2
2
0
2
2
2
2
Θ
Θ
−−=
++
θ
θ
θ
θ
πε
µ
d
d
d
dm
r
e
Er
dr
dR
r
dr
d
R
=
[33]
mesma a é separação de constante a igualdade a devido
)1(sen
1
sen
4
2
1
2
2
0
2
2
2
2
+=
Θ
Θ
−−=
++
ll
d
d
d
dm
r
e
Er
dr
dR
r
dr
d
R
θ
θ
θ
θ
πε
µ
=
[34]
)1(
4
2
1
0
2
2
2
2
+=
++
ll
r
e
Er
dr
dR
r
dr
d
R
πε
µ
=
[35]
Ou
)1(sen
1
sen
2
2
+=
Θ
Θ
−−
ll
d
d
d
dm
θ
θ
θ
θ
[36]
Divide a [35] por
r
2
e arrumando fica
0)1(
4
2
1
)1(
4
2
1
0
2
2
2
2
2
0
2
2
2
2
2
=
++
++
+=
++
Rll
r
e
E
dr
dR
r
dr
d
r
ou
RllR
r
e
E
dr
dR
r
dr
d
r
πε
µ
πε
µ
=
=
[37]
Esta é a equação Radial – cujas soluções são as funções de Laguerre. Apresentaremos a
solução que é: Solução particular:
)(
2
)(
12
1
2
1
2/3
ξξ
ξ
+
+
−
−=
l
n
l
n
Le
na
rR
[38]
Sendo
[
]
∑
−−
=
++
+
++−−−
+
−
+
−−
=
1
0
2
1
2/1
3
12
!)!12()!1(
)!(
)1(.
])
As soluções desta equação são as funções associadas de Legendre. Os harmônicos esféricos
– são soluções das equações diferenciais: )().(
θ
ϕ
Θ
Φ
representadas por
[]
θξ
ξ
ξ
ξ
θ
ξ
θθ
θ
π
ϕθ
ϕ
cos
)1(
!2
1
)(
)(cos
)(
sen)(
)(
)!(4
)!(12
),(
2
2/1
=
−
=
=
+
−+
=
onde
d
d
l
P
d
Pd
Ponde
eP
ml
mll
Y
l
ll
l
l
m
l
m
m
m
l
imm
llm
A solução geral para o átomo de Hidrogênio é :
),()(),,(
ϕθϕθ
lmnnlm
YrRr =Ψ
e
-(i/
=
)Et
Que pode ser escrita assim:
Com
na
r2
=
ξ
n =1,2,3... número quântico principal
l =0,1,2,3... número quântico de momento angular
m = l,-l+1,...l-1,l número quântico magnético
2.3 Números Quânticos
Cada conjunto de n, l, m define uma função de onda – que é um estado eletrônico do
átomo.
n = 1 camada k
n = 2 camada l
n = 3 camada m
n = 4 camada n
[
]
t
n
Ei
lm
ln
p
p
p
l
ln
l
nlm
eY
pplpln
ln
lnn
ln
Le
na
r
)/(
1
0
2
1
2/1
12
2
1
2/3
),(.
!)!12()!1(
)!(
)1(.
! ][2
!)1(
)()
2
(),,(
=−
−−
=
+
+
+
−
∑
++−−−
+
−
+
−−
−=Ψ
ϕθ
ξ
ξξϕθ
ξ
n = 5 camada o
As funções de onda de cada camada são chamadas de orbitais.
a) para cada valor de n, há n-1 valor de l
n = 1 l = 0
n = 2 l = 1
n = 3 l = 2
n define o estado de energia
b) para cada valor de l há (2l+1) valor de m. -l< m < l
l = 0 m = 0
l = 1
=
=
−=
1
0
1
m
m
m
c) para cada valor de n há n
2
autofunções. Exemplo, se n = 3, há 9 autofunções
l = 2 -l<m<l
=
=
=
−=
−=
2
1
0
1
2
m
m
m
m
m
l = 1 -l<m<l
=
=
−=
1
0
1
m
m
m
l = 0 m = 0
Concluindo, quando resolvemos a equação de Schrödinger para o átomo de
hidrogênio, obtemos as energias para os níveis energéticos (E
n
) e os orbitais atômicos
(Ψ
n,l,m
). Temos então um diagrama de níveis energéticos no interior do átomo, i. e.,
dizemos que um átomo se constitui de um sistema energético quantizado. Na Figura 3
demonstramos soluções da equação de Schroedinger radial.
E foi no Annales 80, que Schrödinger publicou o seu estudo do efeito Stark através de sua
Mecânica Quântica Ondulatória (1). Um dos seus primeiros trabalhos utilizando a nova
teoria.
Figura 3: Soluções da equação de Schroedinger radial.
2.4 Função de Onda na Concepção de Max Born
Born (1928) deu um passo a esta dificuldade propondo uma interpretação estatística
das funções de ondas do elétron, à qual, devido às inúmeras vantagens apresentadas, tem
sido amplamente aceita. Born supôs que as ondas não têm existência real, e assim, as define
como ondas de probabilidade.
O produto Ψ
*
Ψ ou
2
Ψ
em um ponto representa a densidade de probabilidade de
encontrar o elétron, ou um outro corpúsculo qualquer, em um ponto x, y, z, num dado
instante t e igualdade
dvdxdydz
22
Ψ=Ψ
Representando a densidade de probabilidade de encontrar o mesmo elétron em um
elemento de volume dv, e também o número de elétrons dentro do mesmo volume. Esta
interpretação teve um pleno acordo com as condições de Schrödinger (5).
Em processos vibratórios o conhecimento da amplitude é importante como o
conhecimento da freqüência própria; analogamente, é de se esperar que, em mecânica
ondulatória, esteja ligado um importante significado físico à função de onda Ψ ou antes, ao
quadrado do seu módulo, visto ser evidente que o valor instantâneo da própria função
oscilatória não pode desempenhar qualquer papel em virtude da sua alta freqüência. O
motivo por que se toma o quadrado do módulo é que a própria função de onda (devido ao
coeficiente imaginário da derivada em ordem ao tempo da equação diferencial) é uma
quantidade complexa, enquanto as grandezas suscetíveis de interpretação física devem
evidentemente ser reais (6).
Suponhamos que no estado caracterizado pela função de onda Ψ
1
se efetua uma
medição que conduz com certeza a um determinado resultado, e que o mesmo fazendo o
estado no estado Ψ
2
, conduz ao resultado 2. Admite-se então a combinação linear de Ψ
1
e
Ψ
2
, o que significa que toda função de forma CΨ
1
+ CΨ
2
(C
1
e C
2
, constantes) representa
um estado em que a mesma medição pode dar um resultado 1 ou o resultado 2. Podendo
afirmar, que se conhecemos a dependência dos estados com respeito ao tempo, dependência
a qual é dada pela função Ψ
1
(x,t) e em outro, por Ψ
2
(x,t), pode-se notar que qualquer
combinação linear destas dá também a possível dependência de um estado do tempo. Estas
afirmações constituem o conteúdo do princípio de superposição dos estados – um princípio
positivo fundamental de mecânica quântica (7).
2.5 Valor Esperado da Função de Onda
Consideramos uma partícula e onda associada, a função Ψ
(r,t)
e se essa função não se
anula num intervalo entre r e r + dr, na medida de sua posição há uma probabilidade finita
dessa partícula ser encontrada. Não podemos, atribuir a coordenada um valor bem definido,
no entanto, é possível especificarmos uma posição média da partícula.
Imaginemos a medida da posição da partícula no instante t, a probabilidade de
encontrá-la entre r e r + dr é dada pela equação.
),(
),(
),( tr
tr
tr
P
ψψ
∗
=
[16]
onde P
(r,t)
é a probabilidade de encontramos a partícula. Repetindo essa experiência a uma
certa freqüência no mesmo instante e registrando os valores de P
(r,t)
, podemos usar a média
dos valores observados para caracterizar a posição da partícula no instante t.
Este valor é representado por <r>, valor esperado da coordenada r. Veja abaixo
como podemos demonstrar matematicamente este cálculo.
∫
+∞
∞−
>=< drrPr
tr ),(
[17]
como:
),(
),(
),( tr
tr
tr
P
ψψ
∗
=
,
substituímos (16) em (17) temos que:
∫
+∞
∞−
∗
>=< drrr
tr
tr
),(
),(
ψψ
[18]
2.6 Aplicação da Equação de Schrödinger
Tabela 1: Exemplos do operador Hamiltoniano para o movimento de uma
partícula de massa m em diferentes campos de força definidos pela função
(operador) potencial V.
Operador
(a) Partícula livre
V=0
x
2
22
2
ˆ
dx
d
m
H
=
−=
(b) barreira de
potencial
V=V
0
x
0 a
2
22
2
ˆ
dx
d
m
H
=
−=
(x<0; x>0)
0
2
22
2
ˆ
V
dx
d
m
H +−=
=
(0<x<a)
(c) Oscilador
harmônico
V=(1/2)Kx
2
x
K
2
2
22
2
1
2
ˆ
Kx
dx
d
m
H +−=
=
Tabela 2: Exemplos do operador para átomos e moléculas.
Operador H
(a) Átomos de
um elétron
r
Ze
m
H
2
2
2
2
ˆ
−∇−=
=
(b) Átomos de
muitos elétrons
∑∑∑∑
>==
+−∇−=
iij
ij
n
j
j
n
i
r
e
r
Ze
m
H
2
1
2
1
2
2
2
ˆ
=
(c) moléculas
∑∑ ∑∑∑∑∑∑
>>==
++−∇−∇−=
iij
ij
n
ii
i
i
n
i
r
eZZ
r
e
r
eZ
mm
H
ααβ
αβ
αα
β
α
α
α
α
2
2
1
2
2
2
1
2
2
2
1
2
ˆ
==
2.6.1 Poço de Potencial não relativístico
Considerando um potencial degrau unidimensional definido por
(
)
0,0
<
=
xxV
(
)
0,
0
>
=
xVxV
Ψ
Ι
Ψ
ΙΙ
V
0
x
(0)
= 0
x
Supor que uma energia incidente da esquerda para a direita tem energia E = 4V
0
. Neste
exemplo podemos calcular a probabilidade de que a onda será refletida (coeficiente de
reflexão). Analisando o contorno do problema, para
0)(
=
xV quando x<0, o potencial na
região I é nulo e consideramos a partícula livre. A equação de Schrödinger se reduz à.
I
I
E
x
m
Ψ=
∂
Ψ∂
−
2
2
2
2
=
I
I
E
m
x
Ψ=
∂
Ψ∂
−
22
2
2
=
0
2
22
2
=Ψ+
∂
Ψ∂
I
I
E
m
x
=
sendo
2
2
1
2
=
mE
k =
0
2
1
2
2
=Ψ+
∂
Ψ∂
I
I
Ek
x
As raízes da equação:
0
2
1
2
=+ k
λ
,
2
1
2
k−=
λ
,
1
ik
±
=
λ
xikxik
I
BeAe
11
−
+=Ψ
sendo
→
→
−
refletidaondaBe
incidenteondaAe
xik
xik
1
1
Solução geral
=/iEt
II
e
−
Ψ=Ψ
(
)
=/
11
iEt
I
xikxik
I
eBeAe
−
−
Ψ+=Ψ
(
)
(
)
== /
1
/
1
EtxkiiEtxki
I
BeAe
+−−
+=Ψ
Na região II há um potencial )(
xV , sendo assim a partícula sofre uma ação do potencial.
IIIIII
EV
m
Ψ=Ψ+Ψ∇−
0
2
2
2
=
[]
IIII
VE
m
Ψ−=Ψ∇−
0
2
2
2
=
[]
0
2
0
2
2
=Ψ−+Ψ∇
IIII
VE
m
=
Considerando
[]
0
2
2
2
VE
m
k
II
−=
=
0
22
=Ψ+Ψ∇
IIIIII
k
Achando as raízes da equação
0
2
2
2
=+ k
λ
,
2
2
2
k−=
λ
,
2
ik±=
λ
(
)
xikxik
II
DeCe
22
−
+=Ψ
A solução aceitável é
xik
II
Ce
2
=Ψ
Solução geral
=/iEt
IIII
eΨ=Ψ
=/
2
iEt
xik
II
eCe=Ψ
(
)
=/
2
Etxki
II
Ce
−
=Ψ
logo temos as duas funções para as duas regiões I e II
xikxik
I
BeAe
11
−
+=Ψ
para
0
≤
x
xik
II
Ce
2
=Ψ
para
0
≥x
Analisando as condições de contorno
a)
00
||
==
Ψ=Ψ
xIIxI
xikxikxik
CeBeAe
211
=+
−
A+B=C deve ser contínua no ponto x=0
b) A sua derivada também
00
||
==
∂
Ψ∂
=
∂
Ψ∂
x
II
x
I
x
x
xikxikxik
CeikBeikAeik
2
2
1
1
1
1
=−
−
(
)
CikBAik
21
=
−
()
C
k
k
BA
1
2
=−
Montando o sistema
Achando B Achando C
=−
=+
+
C
k
k
BA
CBA
1
2
CC
k
k
A +=
1
2
2
+=
1
2
1
2
k
k
C
A
−⋅=−
=+
+
)1(
1
2
C
k
k
BA
CBA
−=
1
2
12
k
k
CB
−
=
1
21
2 k
kk
C
B
CB
k
kk
C
=+
+
1
21
2
0
2
1
1
21
=
+
+
k
kk
C
B
+
−=
1
21
2
k
kk
BC
Substituindo os valores de A, B e C em
I
Ψ
e
II
Ψ
() ()
== /
1
1
21
/
1
1
21
22
EtxkiEtxki
I
e
k
kk
C
e
k
kk
C
−−−
−
+
+
=Ψ
x<0
() ()
−
+
+
=Ψ
−−− == /
1
1
21
/
1
1
21
2
EtxkiEtxki
I
e
k
kk
e
k
kkC
x<0
()
=/
2
1
21
2
Etxki
II
e
k
kk
B
−
−
−=Ψ
Obedecendo ao princípio da complementaridade temos a combinação linear do dois estados
da função de onda
IIITotal
Ψ
+
Ψ
=
Ψ
() () ()
=== /
2
1
21
/
1
1
21
/
1
1
21
2
2
EtxkiEtxkiEtxki
Total
e
k
kk
Be
k
kk
e
k
kk
C
−−−−
−
−
−
+
+
=Ψ
2.7 BIBLIOGRAFIA:
(1) Bassalo, J. M. F. Nascimentos da Física (1901-1950), Belém: EDUFPA,
2000. 503 p
(2) Peixoto, E.M.A. Química Nova, (
1978)
(3) Eisberg, R. Resnick, R. Física Quântica, Ed. Campus, Rio de Janeiro,
1979. 15
ª
Edição.
(4) Bathista, A.L.B.S., Nogueira, J. S. Uma Breve Discussão da Mecânica Quântica IX
Encontro de Iniciação Científica, Cuiabá, UFMT.
2001.
(5) Rey, A.B. Mecânica Quântica e Ondulatória, In: Fisica/Química Modernas. V. 3, SP,
Ed. Fortaleza,
1970.
(6) Born, M., Física Atômica,
1962, 4
ª
Edição. Ed: Fundação Calouste Gulbenkian, Lisboa.
(7) Landau, L. Lifshitz, E. Mecânica Quântica, Teoria não relativista,
1985. Mir, Moscou.
Vol.1.
Capítulo 3: Método de Hückel.
3.1 Método Semi-Empírico
A equação de Schrödinger para moléculas não pode ser exatamente resolvida. Isto
ocorre porque a equação exata não é separável. Uma estratégia é de assumirmos que
podemos permitir de escrever uma forma aproximativa de equação de Schrödinger para
moléculas e principalmente que esta nova equação de Schrödinger se torne separável. E há
então de procurar como resolveremos esta equação separável. O método de Hückel é uma
possibilidade e de campo autoconsistente é um outro.
3.2 Elaboração
Separação
π
σ
−
Para o tratamento de moléculas com elétrons π é possível escrever o Hamiltoniano total
como:
πσσπ
HHHH
+
+
=
[1]
o método de Hückel extendido é
)()()()(
σσσσ
Ψ=Ψ EH
[2]
e o método de Hückel simples é
)()()()(
ππππ
Ψ=Ψ EH
[3]
O Hamiltoniano que descreve a energia do elétron π é
)()(
2
1
)(
2
iViiH
ππ
+∇−= [4]
Sendo o Hamiltoniano total
∑
=
n
i
Total
iHiH )()(
ππ
[5]
∑
+∇−=
n
i
Total
iViiH )()(
2
1
)(
2
ππ
[6]
A função
π
Ψ
é antissimétrica de n órbitas moleculares π. A função antissimétrica
π
Ψ descreve partículas que obedecem a estatística de Fermi-Dirac, para Férmions. O
método de Hückel é um método classificado como semiempírico. Este pode ser
simplificado no máximo a resolução das equações
∫
ΨΨ=
τ
dHE
*
, ou seja ele utiliza do
método variacional para solucionar problemas de moléculas com orbitais π, e para isso ele
faz algumas considerações.
3.3 Considerações do Método de Hückel
i)
ijij
S
δ
=
onde δ é o delta de Krönecker.
ii)
α
=
ii
H
onde α é um parâmetro coulombiano
iii)
β
=
ij
H
se i e j pertencem a átomos ligados entre si e
Estas considerações são realizadas sobre a determinante secular
0=−
ijij
ESH
[7]
3.4 Aplicação do Método de Hückel no Etileno
CH
2
CH
2
A determinante secular
π
H
para o etileno
)1(1
2
z
p
=
Ψ
e
)2(2
2
z
p
=
Ψ
será 2x2
0
22212121
12121111
=
−−
−−
ESHESH
ESHESH
Utilizando-se as condições de Hückel
0
1
=
−
−
÷
β
αβ
βα
E
E
0
1
1
=
−
−
β
α
β
α
E
E
,
β
α
E
x
−
=
0
1
1
=
x
x
Esta é conhecida como determinante de Hückel.
01
2
=−x
,
1
±
=
x
para x=1
β
α
−=E e para x=-1
β
α
+
=
E
=
⋅
0
0
1
1
2
1
c
c
x
x
=+
=+
0
0
21
21
xcc
cxc
, c
1
=c
2
Quando normalizamos as funções
∫
=ΨΨ 1
*
τ
d
nós achamos os valores de c
1
(
)
(
)
1
2
2
1
1
*
22
*
11
=Ψ+ΨΨ+Ψ
∫
τ
dcccc
1
2
2
2
1
0
*
2
1
2
1
0
2
*
1
2
1
2
1
2
1
=Ψ+ΨΨ+ΨΨ+Ψ
∫∫∫∫
ττττ
dcdcdcdc
2
1
1
=c
logo achamos a funções
1
Ψ e
2
Ψ
)(707.0)(707.0
211 zz
pp
Ψ
+
Ψ
=Ψ
)(707.0)(707.0
212 zz
pp
Ψ
−
Ψ
=Ψ
-2 -1 0 1 2
-1
0
1
2
3
Etileno
Figura 1: comportamento da função de segundo grau do Etileno
3.5 Aplicação do Método de Hückel no Butadieno
A molécula 1,3 Butadieno é plana e tem preferência como isômero trans e tem
fórmula
0
4444434342424141
3434333332323131
2424232322222121
1414131312121111
=
−−−−
−−−−
−−−−
−−−−
ESHESHESHESH
ESHESHESHESH
ESHESHESHESH
ESHESHESHESH
Fazendo as considerações do método de Hückel
0
00
0
0
00
1
=
−
−
−
−
÷
β
αβ
βαβ
βαβ
βα
E
E
E
E
0
100
110
011
001
=
−
−
−
−
β
α
β
α
β
α
β
α
E
E
E
E
sendo
β
α
E
x
−
=
0
100
110
011
001
=
x
x
x
x
quando se tem uma determinante diagonal simétrica, utilizamos a fórmula abaixo
()
∏
=
+
⋅
−
n
j
n
j
bca
1
2/1
1
cos2
π
sendo j = 1,2,..n
0
100
110
011
001
=
ac
bac
bac
ba
x
x
x
x
0134
24
=+− xx
as raízes são
62.0±=x
e
62.1±=x
βα
βα
βα
β
α
62.1
62.0
62.0
62.1
4
3
2
1
−=
−=
+=
+
=
E
E
E
E
α − 1.62 β
α − 0.62 β
α + 0.62 β
α + 0.62 β
C
2p
4π
∗
3π
∗
2π
1π
=
⋅
0
0
0
0
100
110
011
001
4
3
2
1
c
c
c
c
x
x
x
x
=+
=++
=++
=+
0
0
0
0
43
432
321
21
xcc
cxcc
ccc
cxc
Agora que temos o sistema montado para o Butadieno, podemos calcular todas as funções
de onda relacionadas com as ligações do Butadieno. Iniciando a calcular
1
Ψ
o orbital de
menor energia, utilizamos x= -1.62, vem:
=
=
=
=+−=
=
41
32
43
1213
12
62.1
62.162.1
62.1
cc
cc
cxc
cccc
cc
Agora com a condição de normalização de
1
Ψ
(para S
ij
=0)
1))((
44332211
*
4
*
4
*
3
*
3
*
2
*
2
*
1
*
1
=Ψ+Ψ+Ψ+ΨΨ+Ψ+Ψ+Ψ
∫
τ
dcccccccc
43
3
2
2
1
62.162.1 cffcfccf
cc
i
+
+
+
=Ψ
1
2
4
2
3
2
2
2
1
=+++ cccc
(
)
(
)
161.162.1
2
4
22
2
1
=+++ cccc
372.0124,7
2
=→= cc
substituindo c em
1
Ψ , temos o resultado da primeira combinação
43211
372.0602.0602.0372.0
ψ
ψ
ψ
ψ
+
+
+
=Ψ
e assim para acharmos
1
Ψ ,
2
Ψ e
3
Ψ
devemos substituir no sistema acima apresentado.
-2-1012
0,1
1
10
100
Y
Coeficiente C
i
Representação
de 4x
4
-3x
2
+1
Figura 3: comportamento do determinante do butadieno
3.6 Aplicação do Método de Hückel para o Ciclobutadieno
A molécula do Ciclobutadieno difere da molécula do Butadieno porque esta possui
os carbonos ligados 2 em 2. lembrando que a determinante secular é
0
4444434342424141
3434333332323131
2424232322222121
1414131312121111
=
−−−−
−−−−
−−−−
−−−−
ESHESHESHESH
ESHESHESHESH
ESHESHESHESH
ESHESHESHESH
0
0
0
0
0
1
=
−
−
−
−
÷
β
αββ
βαβ
βαβ
ββα
E
E
E
E
0
101
110
011
101
=
−
−
−
−
β
α
β
α
β
α
β
α
E
E
E
E
sendo
β
α
E
x
−
=
0
101
110
011
101
=
x
x
x
x
;
24
4xx −
20 ±== xex
βα
α
α
β
α
2
2
4
3
2
1
−=
=
=
+=
E
E
E
E
onde E
2
e E
3
são degenerados
α−2β
α
α+2β
=
⋅
0
0
0
0
101
110
011
101
4
3
2
1
c
c
c
c
x
x
x
x
=++
=++
=++
=++
0
0
0
0
431
432
321
321
xccc
cxcc
cxcc
ccxc
-2 -1 0 1 2
-4
-3
-2
-1
0
Y
X
Figura 5: comportamento do determinante do ciclobutadieno
3.7 Referência:
(1) Vidal, B. Chimie Quantique: de l’atome à la théorie de Hückel. Ed. Masson: Paris, 1993
(2) Gil, V. M. S., Orbitais em Átomos e Moléculas, Fundação Calouste Gulbenkian, 1996
(3) Rogers, D. W. Computational Chemistry Using the PC, 2º Ed. VCH: New York, 1994
(4) Bernath, P. F. Spectra of Atoms and Molecules, Oxford Univesity Press, New York, 1995
(5) Levine, I. N. Quantum Chemistry, 4º Ed. Prentice Hall, New Jersey, 1991
(6) Hanna, M. W. Quantum Mechanics in Chemistry, 2º Ed. Benjamin, Colorado, 1969
(7) Anderson, P. W. Phys. Rev. 181, 25 (1969)
Capítulo 4: Método Variacional.
4.1 Introdução
O método variacional trata-se de um princípio importante em Química Quântica,
não só porque permite obter soluções aproximadas para a equação de Schrödinger como
também porque fornece um critério, o valor mínimo de energia. Considerando a equação de
Schrödinger.
Ψ
=
Ψ
E
H
[0]
onde
Ψ
e E são autofunção e autovalor do Hamiltoniano, são soluções exatas desta
equação. Normalmente o Hamiltoniano é de tal forma que a equação [0] não é separável,
qualquer que seja o sistema de coordenadas. Por essa razão não é possível pensarmos em
obter uma expressão analítica exata para
Ψ
. Passemos então a examinar as propriedades
da funcional )(
ψ
E , definida como
τ
τ
τ
dH
d
dH
E
∫
∫
∫
ΨΨ=
ΨΨ
ΨΨ
=
*
*
*
[1]
Esta é a integral variacional. Uma função qualquer, que seja normalizada, bem comportada
e que satisfaz as condições de contorno do problema, podemos mostrar que esta função
admite que o valor esperado do operador Hamiltoniano para uma função de onda
aproximada é sempre maior que a energia exata
E
0
correspondente a esse Hamiltoniano.
∫
≥ΨΨ
0
*
EdH
τ
[2]
Em outras palavras, a função de onda aproximada conduz à menor energia. Assim sendo, os
coeficientes
k
c
nas orbitais moleculares
Ψ
são determinados de modo a conduzirem ao
valor mínimo de energia
E associada com essa orbital molecular. Em relação à obtenção de
soluções aproximadas, o que se faz na prática consiste em partir de uma função aceitável,
com vários parâmetros ajustáveis, e variar esses parâmetros de modo minimizar
E pois,
segundo o princípio variacional, conseguirá obter-se uma melhor aproximação para a
energia do estado fundamental
E
0
∑
∑
∑∑
∫
∑∑
∫
∫
∑∑
∫
∑∑
=
ΨΨ
ΨΨ
=
ΨΨ
ΨΨ
=
kj
kjjk
kj
kjjk
kj
jkjk
kj
jkjk
k
k
k
k
k
k
k
k
Scc
Hcc
cc
Hcc
dcc
dcHc
E
*
*
*
*
**
**
ˆˆ
τ
τ
Provando a afirmação [2],
k
k
k
a Ψ=Ψ
∑
[3]
Substituindo [3] no lado esquerdo de [2]
∫
∑
∫
∑
∑
∫
∑
ΨΨ=ΨΨ=ΨΨ
j
jj
k
k
k
j
jj
k
k
k
dHaadaHadH
τττ
ˆˆ
*****
∫
∑
∑
∫
∑
∑
∑
∫
∑
Ψ=ΨΨ=ΨΨ=ΨΨ
kj
kjjjjkj
k
k
jj
j
k
j
jjj
k
k
k
EaadEaadEaadH
δτττ
******
∫
∑
∑
Ψ=ΨΨ
kj
kjjjjk
EaadH
δτ
**
[4]
onde
kj
δ
é o delta de Kroenicker,
≠
=
kse
jkse
0
1
como temos
kj
δ
podemos conseguir o termo
∑
k
kkk
Eaa
*
∑
∫
∑
==ΨΨ
k
kk
k
kkk
EaEaadH
2
**
τ
[5]
sendo que o nível k é sempre maior
0
EE
k
>
, onde E
0
é o estado fundamental do sistema.
∑∑
∫
∑
=≥=ΨΨ
k
k
k
k
k
kkk
aEEaEaadH
2
00
2
**
τ
[6]
Através da normalização
∫
=ΨΨ 1
*
τ
dH [7]
agora substituindo
Ψ
∫
∑
∑
∑
∑
∫
ΨΨ=ΨΨ=
kj kj
jkjkjjkk
daadaa
ττ
****
1
1
**
=ΨΨ
∑∑
∫
kj
jkjk
daa
τ
[8]
abrindo a somatória em
j
∑∑∑
=
k
kkj
kj
jk
aaa
2
*
δ
[9]
Substituindo [9] em [6]
∑
∑
∫
=ΨΨ
kj
kjjk
aaEdH
δτ
*
0
*
[10]
∫
≥ΨΨ
0
*
EdH
τ
[11]
Isso mostra a prova. Por isso temos que ajustar bem a solução de
Ψ
muito bem para
atingirmos o estado
0
E
(valor real) que é o estado fundamental. Ajustar a função
Ψ
é
ajustar a função de base computacionalmente.
4.2 Extensão do método Variacional
O método variacional apresentado na seção passada tem duas limitações. Primeira,
ele provê informação de energia somente sobre o estado fundamental e da função de onda.
Segundo, ele provê um estado acima do estado fundamental somente. Nesta seção
gostaríamos de estender o método variacional.
Considerando que nós desejamos estender o método variacional para obter um valor
de energia estimada para o primeiro estado excitado.
Demonstrado que a função variacional normalizada [5] e [8]
∑
∫
=
=ΨΨ
1
2
*
k
kk
EadH
τ
e
1
1
2
*
==ΨΨ
∑
∫
=k
kjk
ad
τ
onde
k
a
é a expansão do coeficiente em
k
k
k
a Ψ=Ψ
∑
. Tendo
ΨΦ=
kk
a
. A função
Ψ
que é ortogonal para o estado fundamental real da função de onda
0
Φ
. Então nós temos
0
00
=ΨΦ=a
e
e
∑
∫
=
=ΨΨ
2
2
*
k
kk
EadH
τ
e
1
2
2
*
==ΨΨ
∑
∫
=k
kjk
ad
τ
[12]
para 1
≥k , temos
1
EE
k
≥
e
1
22
EaEa
kkk
≥
∑∑∑
∞
=
∞
=
∞
=
==≥
1
1
2
1
1
1
2
1
2
k
k
k
k
k
kk
EaEEaEa [13]
Combinando [12] e [13], temos o resultado:
∫
≥ΨΨ
1
*
ˆ
EdH
τ
, se
∫
=ΨΨ 0
ˆ
*
τ
dH
e
∫
=ΨΨ 1
ˆ
*
τ
dH
[14]
a extensão do método variacional para altos estados excitados acima de E
1
.
1
*
*
ˆ
+
≥
ΨΨ
ΨΨ
∫
∫
k
E
d
dH
τ
τ
, se 0
**
1
*
0
∫
∫
∫
=ΦΨ=⋅⋅⋅=ΦΨ=ΦΨ
τττ
ddd
k
[15]
4.3 Função Variacional Linear
Um tipo especial de função variacional é usado largamente no estudo de moléculas,
é a
Função Variacional Linear ou Método Variacional de Raleigh-Ritz. Uma função
variacional linear é uma combinação de
n funções linearmente independente
n
fff ,...,
21
:
∑
=
=++=Ψ
n
j
jjnn
fcfcfcfc
1
2211
... [16]
onde
Ψ
é a função trial variacional e os coeficientes
j
c são parâmetros para ser
determinados pela minimização da integral variacional. A função
j
f deve satisfazer a
condição de contorno do problema. Nós apresentamos que
Ψ
,
j
c
e
j
f são reais.
Aplicando o teorema variacional para a função variacional linear, nós temos
∫
∑
∫
∑∑∑
∫
====
==ΨΨ
n
j
kj
n
k
kj
n
k
kkj
n
j
j
dffccdfcfcd
1
*
1
*
1
*
1
**
τττ
[17]
Agora devemos definir a integral de superposição
jk
S
como
∫
=
τ
dffS
kjjk
*
[18]
sendo
jk os elementos de matriz. Substituindo [18] em [17] temos
∑∑
∫
==
=ΨΨ
n
j
n
k
jkkj
Sccd
11
**
τ
[20]
Note que
jk
S
não é necessariamente igual a
jk
δ
, desde não há razão para supor que a
função
j
f
são mutuamente ortogonais. O numerador da integral variacional é
∫∫
∑∑ ∑∑
∫
== ==
==ΨΨ
n
j
n
k
n
j
n
k
kjkjkkjj
dHffccdfcHfcdH
11 11
*****
τττ
e usando a abreviação
∫
≡
τ
dHffH
kjjk
*
[21]
Reescrevendo o numerador
∫
∑∑
==
=ΨΨ
n
j
n
k
kj
dccdH
11
**
ττ
[22]
Substituindo [22] e [20] na integral variacional W,
∑∑
∑∑
∫
∫
==
==
=
ΨΨ
ΨΨ
=
n
j
n
k
jkkj
n
j
n
k
jkkj
Scc
Hcc
d
dH
W
11
11
*
*
τ
τ
[23]
Arranjando [23]
∑∑∑∑
====
=
n
j
n
k
jkkj
n
j
n
k
jkkj
HccSccW
11
*
11
*
[24]
agora devemos minimizar W para aproximar tão perto de E
0
)(
0
EW ≥
. A integral
variacional W é uma função de
n variáveis independentes
n
ccc ,...,
21
),...,(
21 n
cccWW
=
Uma condição necessária para minimizar uma função W de muitas variáveis é que a
sua derivada parcial a respeito a cada uma de sua variável seja zero no ponto mínimo.
,0=
∂
∂
i
c
W
ni ,...2,1
=
[25]
agora nos diferenciamos [24] parcialmente em cada c
i
para obter n equações:
∑∑ ∑∑ ∑∑
== == ==
∂
∂
=
∂
∂
+
∂
∂
n
j
n
k
n
j
n
k
n
ij
n
ik
jkkj
i
jkkj
i
jkkj
i
Hcc
c
Scc
c
WScc
c
W
11 11
[26]
i= 1,2,...,n
Agora
()
jk
n
j
n
j
n
k
n
j
n
ki
k
j
i
j
kjkkj
i
jk
n
k
kj
i
S
c
c
c
c
c
cScc
c
Scc
c
∑∑∑ ∑∑∑
=== ===
∂
∂
+
∂
∂
=
∂
∂
=
∂
∂
111 111
os c
j
são variáveis independentes, e portanto
0=
∂
∂
i
j
c
c
se i
j
≠
1=
∂
∂
i
j
c
c
se i
j
=
ij
i
j
c
c
δ
=
∂
∂
Nós temos então
∑∑∑∑∑∑∑∑
========
+=+=
∂
∂
n
k
n
j
ji
jikkjk
n
j
n
j
n
k
n
j
n
k
ikjjkijkjk
n
k
kj
i
cScScScScc
c
11111111
δδ
∑∑∑∑
====
+=
∂
∂
n
k
n
j
jijikk
n
j
jk
n
k
kj
i
ScScScc
c
1111
[27]
Onde nós temos desenvolvido
ijijji
SSS ==
*
[28]
Resolvendo a equação [27] a partir desta relação, conseguimos compacta-la:
∑∑∑∑∑∑
======
+=+=
∂
∂
n
k
n
k
ikkikk
n
k
n
j
ijjikk
n
j
jk
n
k
kj
i
ScScScScScc
c
111111
∑∑∑
===
=
∂
∂
n
k
ikk
n
j
jk
n
k
kj
i
ScScc
c
111
2
[29]
onde o fato de
j ser uma variável aleatória utilizada e recolocando
jk
S
por
jk
H
em cada
uma destas manipulações, nós temos o numerador da equação vai
∑∑∑
===
=
∂
∂
n
k
ikk
n
j
jk
n
k
kj
i
HcHcc
c
111
2
[30]
o resultado depende de
ijijji
HHH ==
*
Porque H é Hermitiano, agora substituindo [25], [29] e [30] em [26]
niHcScW
n
k
ikk
n
k
ikk
,...2,1,22
11
==
∑∑
==
Rearranjando
ou
()
[]
nicWSH
n
k
kikik
,..2,1,0
1
==−
∑
=
[31]
Sendo esta conhecida como Sistema de Equações Seculares. A equação [31] é um
conjunto de
n equações simultâneas, lineares e homogêneas em n de c
1
, c
2
,...c
n
desconhecidos. E pode ser representada como
0)( =−=
ikikm
WSHWD
, onde m é a linha
da determinante.
Uma solução trivial das equações seculares é todos os coeficientes
k
c
serem nulos.
Qualquer solução fisicamente relevante exige que seja nulo o chamado determinante
secular, constituído pelos elementos
ikik
ESH
−
0=−
ikik
ESH
[32]
ou
0)( ==
ikik
ESHED
da equação [31] podemos obter os valores reais de
mmmm
WWWW ,...,,
21
=
os quais
correspondem
mmmm
ΨΨ
Ψ
=Ψ ,...,,
21
. Onde o valor esperado da energia
Ψ
é facilmente
visto por
∫
∫
ΨΨ
ΨΨ
=
τ
τ
d
dH
E
*
*
ou
SCC
CC
+
+
=
H
E
4.4 Aplicação do Método Variacional
Sabemos que agora temos a parte secular da integral variacional e podemos a partir
daqui aplicar o método variacional no estudo de energias de orbitais atômicos e
moleculares.
0
0
*
=−→≥ΨΨ=
∫
iKik
SecularEquação
ESHEdHE
τ
Para átomo as funções de base são os elétrons e para molécula as funções de base
são os orbitais moleculares. Ex: para o Benzeno há 36 orbitais, ou seja 36 combinações
lineares. Logo
φ
Ψ será
∑
=
=Ψ
36
1
i
ii
c
φ
φ
e a sua determinante será de 36x36. há 42 elétrons no benzeno e 21 orbitais estão ocupados
e outros 21 são virtuais. Sendo que o último orbital ocupado é chamado de HOMO
(Highest Occupied molecular Orbital) e o primeiro orbital não ocupado é o LUMO (Lowest
Unccupied Molecular Orbital). Se soubermos a diferença de energia LUMO e HOMO,
podemos achar a energia de ionização da molécula.
A energia total dos orbitais ocupados é da seguinte forma
()
2
1
⋅++=
∑
=
n
i
ExCiT
EEEE
Onde n é o número de orbitais ocupados, 2 o par de elétrons e E
i
número de orbitais, E
C
repulsão coulombiana e E
Ex
é a energia de troca.
4.5 Método Variacional para o Átomo de Hidrogênio
A partir da função tentativa
r
e
α
α
−
=Ψ
para o estado fundamental do átomo de
Hidrogênio, calcule o valor de
α
que minimiza a energia.
αα
αα
ΨΨ
ΨΨ
=
H
E
[1]
Sendo o Hamiltoniano H,
r
Z
r
rr
r
r
r
r
r
H −
∂
∂
+
∂
∂
∂
∂
+
∂
∂
∂
∂
−=
2
2
222
2
2
sen
1
sen
sen
11
2
1
φθ
θ
θ
[2]
φθθτ
ddrdrd sen
2
=
)0(
!
1
≥=
+
−
∫
a
a
n
drre
n
nr
α
Como o átomo de Hidrogênio depende somente de r, o operador se reduz à
r
Z
r
r
r
r
H −
∂
∂
∂
∂
−=
2
2
1
2
1
[3]
∂
∂
∂
∂
≡∇
r
r
r
r
2
2
2
1
+=∇
dr
d
dr
d
r
r
dr
d
r
r
2
22
2
1
2
1
+=∇
2
2
2
2
dr
d
dr
d
r
[4]
Substituindo [4] em [3]
r
Z
H −∇−=
2
2
1
[5]
Reduzindo a equação [1]
αα
ΨΨ= HE
[6]
e aplicando o operador em [6]
αα
Ψ−∇−Ψ=
r
Z
E
2
2
1
αααα
ΨΨ−Ψ∇Ψ−=
r
ZE
1
2
1
2
αααα
ΨΨ−Ψ+Ψ−=
r
Z
dr
d
dr
d
r
E
12
2
1
2
2
αααααα
ΨΨ−ΨΨ−ΨΨ−=
r
Z
dr
d
dr
d
r
E
1
2
11
2
2
[7]
para que a variação funcional,
)(
α
δ
ΨE
, seja mínima →
0)( ==Ψ HE
δδ
α
sendo que
α
αα
α
αα
Ψ−=−==Ψ
−− rr
ee
dr
d
dr
d
e que
αα
α
Ψ=Ψ
2
2
2
dr
d
,
teremos:
αααααα
α
α
ΨΨ−ΨΨ−ΨΨ=
r
Z
r
E
1
2
1
2
r
Z
r
E −−=
2
2
αα
2
1
)(
2
α
α
−−=
r
ZE
[8]
Normalizando
α
Ψ
1=ΨΨ
αα
lembrando que
)0(
!
1
≥=
+
−
∫
a
a
n
drre
n
nr
α
1
22
=
∫
−
drre
r
α
1
)4(2
!2
3
2
=
α
N
→
3
4
α
=N
[9]
calculando o valor médio de
r
1
∫
−
= dr
r
r
eN
r
r
2
2
1
α
∫
−
= rdre
r
r
α
α
2
3
4
1
⋅=
2
3
)2(
1
4
1
α
α
r
α
=
r
1
[10]
Substituindo [10] em [8]
2
)(
2
α
αα
−−= ZE
α
αα
αα
ZZE −=−−=
22
22
2
[11]
como
0)( ==Ψ HE
δδ
α
, o mínimo de energia será obtido se derivarmos E em relação a
α
e igualando-se a zero,
ZZ
d
dE
=→=−=
α
α
α
0
2
2 [12]
Substituindo [12] em [11], e lembrando que Z=1 para o átomo de hidrogênio.
..
2
1
)1(
2
)1(
2
auE −=−=
..
2
1
auE −= [13]
4.6 Método Variacional para o Átomo de Hélio
Aplicando o método variacional para o átomo de Hélio e encontrando a energia do
estado fundamental a partir da função tentativa.
1
3
)1(1
r
es
α
π
α
=
e
2
3
)2(1
r
es
α
π
α
=
[1]
Um Hamiltoniano de N elétrons e M núcleos é
nnrepulsão
M
k
M
kl
kl
lk
N
i
M
k
eerepulsão
N
i
N
j
ik
nerepulsão
N
i
M
k
ik
k
k
k
i
R
ZZ
rr
Z
M
H
−
=>==
−−
−
=>
−
−
==
∑∑∑∑ ∑∑∑∑
++−∇−∇−=
111 1111
22
1
2
1
2
1
ˆ
[2]
No átomo de Hélio, consideremos os núcleos fixos, podemos omitir o termo de
energia cinética dos núcleos para obter o Hamiltoniano puramente eletrônico, separando-o
do nuclear. Sendo que há 5 termos no Hamiltoniano do Hélio
[
]
ijjkik
jieletrônico
rr
Z
r
Z
H
1
2
1
ˆ
22
−−−∇−∇−= [3]
spinespacialap
Φ×Φ=Ψ
)2()1()2(1)1(1
α
α
×=Ψ ss
ap
, há 4 possibilidades de função de spin desde que
negligenciamos o efeito magnético.
A)
)2(1)1(1
α
α
ss
ap
=
Ψ
B)
)2(1)1(1
β
β
ss
ap
=
Ψ
C)
)2(1)1(1
β
α
ss
ap
=
Ψ
D)
)2(1)1(1
α
α
ss
ap
=
Ψ
Onde A) e B) são simétricas e C) e D) são antissimétricas.
C)
0
)2(1)2(1
)1(1)1(1
!2
1
)2(1)1(1)2(1)1(1 ==Ψ→−=Ψ
βα
βα
αββα
ss
ss
ssss
apap
D)
0
)2(1)2(1
)1(1)1(1
!2
1
)2(1)1(1)2(1)1(1 ==Ψ→−=Ψ
αβ
αβ
βααβ
ss
ss
ssss
apap
Utilizando C)
)2(1)1(1)2(1)1(1
α
β
β
α
ssss
ap
−=Ψ
() ()
)2()1()2(1)1(1)2()1()2(1)1(1
α
β
β
α
ssss
ap
−=Ψ
[]
)2()1()2()1()2(1)1(1
α
β
β
α
−=Ψ ss
ap
[4]
A parte
espacial é simétrica e lembrando que a função de spin já está normalizada
[]
spin
espacial
ap
ss )2()1()2()1(
!2
1
)2(1)1(1
αββα
−=Ψ
[5]
Verificando a ortonormalização da função de spin
[][]
1)2()1(
2
1
)2()2()1()1(
2
1
)2()2()1()1(
2
1
)2()1(
2
1
)2()1()2()1(
!2
1
)2()1()2()1(
!2
1
2
2
2
1
2
2
*
1
2
2
*
1
2
2
2
1
*
12
=+
+−
=−−
∑∑∑∑
∑∑∑∑
∑∑
αββααβ
αββαβα
αββααββα
[6]
∑
−=
=
2/1
2/1
2
1)(
s
m
s
m
α
e
∑
−=
=
2/1
2/1
2
1)(
s
m
s
m
β
[7]
Desde que as funções αe β correspondem a diferentes autovalores do operador Hermitiano
z
S
ˆ
, eles são ortogonais.
∑
−=
=
2/1
2/1
*
0)()(
s
m
ss
mm
βα
[8]
Aplicando estas notações na verificação de ortonormalização
1)1()1(
2
1
)0()0(
2
1
)0()0(
2
1
)1()1(
2
1
=⋅+⋅+⋅−⋅
1=1 [10]
A integral variacional é
()
[]
()
[]
21
*
*
12
*
)2()1()2()1(
2
1
)2(1)1(1)2()1()2()1(
2
1
)2(1)1(1 dvdvssHssH
αββααββα
−−=ΨΨ
∫∫∫
∑∑
Neste caso
H não atua na função de spin e a integral variacional torna-se
()()
∫∫∫
∑∑
−=ΨΨ
1
2
12
21
*
*
)2()1()2()1(
2
1
)2(1)1(1)2(1)1(1
αββα
dvdvssHssH
[11]
Como a função de spin é normalizada, a integral variacional se reduz a
(
)
(
)
∫∫∫
=ΨΨ
21
*
*
)2(1)1(1)2(1)1(1 dvdvssHssH
[12]
Substituindo o Hamiltoniano [3] em [12]
()
[
]
()
∫∫∫
−−−∇−∇−=ΨΨ
21
22
*
*
)2(1)1(1
1
2
1
)2(1)1(1 dvdvss
rr
Z
r
Z
ssH
ijjkik
ji
[13]
()()
21
*
2
*
1
*
2
2*
1
*
1
2*
2
*
)2(1)1(1
1
)2(1)1(1)2(1)2(1)1(1)1(1
)2(1
2
1
)2(1)1(1)1(1)1(1
2
1
)1(1)2(1)1(1
dvdvss
r
ssdvs
r
Z
sdvs
r
Z
s
dvssdvssdvssdvssE
ijji
ji
∫∫∫∫
∫∫∫∫
+++
∇−+
∇−=
Aplicando os operadores nas funções espaciais
Ψ===Ψ
α
π
α
α
π
α
αα
i
r
i
r
ii
ee
dr
d
dr
d
33
,
Ψ−=−==Ψ
−−
α
π
α
α
π
α
αα
j
r
j
r
jj
ee
dr
d
dr
d
33
e que
Ψ=Ψ
2
2
2
α
i
dr
d
,
Ψ=Ψ
2
2
2
α
j
dr
d
Logo a expressão de Energia na notação de Dirac torna-se
ΨΨΨΨΨΨ−ΨΨ−ΨΨ−ΨΨ−ΨΨ−=
ijji
j
jj
i
ii
r
Z
r
Z
r
Z
dr
d
dr
d
r
dr
d
dr
d
r
E
111
2
11
2
11
2
2
2
2
αα
α
α
α
α
α
α
α
)8/5(
1
2
1
2
111
2
1
2
1
ΨΨ+ΨΨ−ΨΨ−ΨΨ−ΨΨ+ΨΨ−ΨΨ−=
ijjiji
r
Z
r
Z
r
Z
rr
E
ααα
8
5
2
2
+−−= zE
como
0)( ==Ψ HE
δδ
α
, o mínimo de energia será obtido se derivarmos E em relação a
α
e igualando-se a zero,
0
8
5
22 =+−−= z
d
dE
α
α
8
5
)(2 −=+− z
α
16
5
+−= z
α
6875.1
=
α
Substituindo o valor de
α na expressão de energia
)6875.1(
8
5
)6875.1)(2(2)6875.1(
2
+−=E
HartreeE 8477.2−=
4.7 Método Variacional para Sistema de dois Spins
ΨΨ
ΨΨ
=
H
E
[1]
Sendo o Hamiltoniano H,
212211
ˆˆˆˆ
IIII ⋅+−−= JH
νν
[2]
Agora analisemos a interação spin-spin e podemos ver que esta atua como uma
perturbação no sistema de spins. E por isso consideramos esta como uma perturbação, logo
podemos escrever o Hamiltoniano [2] da forma:
)1()0(
H
H
H
+
=
[3]
onde
)0(
2
)0(
1
)0(
HHH +=
e
11
)0(
1
ˆ
I
ν
−=H
e
22
)0(
2
ˆ
I
ν
−=H
, o qual a sua função
aproximada é
spinespacialap
Φ
×
Φ
=
Ψ
[4]
)2()1()2(1)1(1
α
α
×=Ψ ss
ap
, há 4 possibilidades de função de spin desde que
negligenciamos o efeito magnético.
a)
)2(1)1(1
α
α
ss
ap
=
Ψ
b) )2(1)1(1
β
β
ss
ap
=
Ψ
c) )2(1)1(1
β
α
ss
ap
=
Ψ
d)
)2(1)1(1
α
α
ss
ap
=
Ψ
Onde a)e b) são simétricas e c) e d) são antissimétricas.
c)
0
)2(1)2(1
)1(1)1(1
!2
1
)2(1)1(1)2(1)1(1
==Ψ→−=Ψ
βα
βα
αββα
ss
ss
ssss
apap
d)
0
)2(1)2(1
)1(1)1(1
!2
1
)2(1)1(1)2(1)1(1
==Ψ→−=Ψ
αβ
αβ
βααβ
ss
ss
ssss
apap
Utilizando c)
)2(1)1(1)2(1)1(1
α
β
β
α
ssss
ap
−
=
Ψ
(
)
(
)
)2()1()2(1)1(1)2()1()2(1)1(1
α
β
β
α
ssss
ap
−
=
Ψ
[
]
)2()1()2()1()2(1)1(1
α
β
β
α
−
=
Ψ ss
ap
[5]
A parte
espacial é simétrica e lembrando que a função de spin já está normalizada
[]
spin
espacial
ap
ss )2()1()2()1(
!2
1
)2(1)1(1
αββα
−=Ψ
[6]
Verificando a ortonormalização da função de spin
[][]
1)2()1(
2
1
)2()2()1()1(
2
1
)2()2()1()1(
2
1
)2()1(
2
1
)2()1()2()1(
!2
1
)2()1()2()1(
!2
1
2
2
2
1
2
2
*
1
2
2
*
1
2
2
2
1
*
12
=+
+−
=−−
∑∑∑∑
∑∑∑∑
∑∑
αββααβ
αββαβα
αββααββα
[7]
∑
−=
=
2/1
2/1
2
1)(
s
m
s
m
α
e
∑
−=
=
2/1
2/1
2
1)(
s
m
s
m
β
[8]
Desde que as funções α e β correspondem a diferentes autovalores do operador Hermitiano
z
S
ˆ
, eles são ortogonais.
∑
−=
=
2/1
2/1
*
0)()(
s
m
ss
mm
βα
[9]
Aplicando estas notações na verificação de ortonormalização
1)1()1(
2
1
)0()0(
2
1
)0()0(
2
1
)1()1(
2
1
=⋅+⋅+⋅−⋅
1=1 [10]
A integral variacional é
()
[]
()
[]
21
*
*
12
*
)2()1()2()1(
2
1
)2(1)1(1)2()1()2()1(
2
1
)2(1)1(1 dvdvssHssH
αββααββα
−−=ΨΨ
∫∫∫
∑∑
Neste caso H não atua na função de spin e a integral variacional torna-se
()()
∫∫∫
∑∑
−=ΨΨ
1
2
12
21
*
*
)2()1()2()1(
2
1
)2(1)1(1)2(1)1(1
αββα
dvdvssHssH
[11]
Como a função de spin é normalizada, a integral variacional se reduz a
(
)
(
)
∫∫∫
=ΨΨ
21
*
*
)2(1)1(1)2(1)1(1 dvdvssHssH
[12]
Substituindo o Hamiltoniano [3] em [12]
()
[
]
()
∫∫∫
++=ΨΨ
21
)1()0(
2
)0(
1
*
*
)2(1)1(1)2(1)1(1 dvdvssHHHssH
[13]
[]
[
]
()
[]
()
2121
*
222
*
1
*
111
*
2
*
)2(1)1(1
ˆˆ
)2(1)1(1
)2(1
ˆ
)2(1)1(1)1(1)1(1
ˆ
)1(1)2(1)1(1
dvdvssJss
dvssdvssdvssdvssE
∫∫
∫∫∫∫
⋅+
−+−=
II
II
νν
4.8 Referências:
(1) Teixeira, J. J. C. D.,
Química Quântica: Fundamentos e Métodos, Fundação Calouste
Gulbenkian,
1982
(2) Gil, V. M. S.,
Orbitais em Átomos e Moléculas, Fundação Calouste Gulbenkian, 1996
(3) Levine, I. N.
Quantum Chemistry, Prentice Hall, 1991
(4) Peixoto, E. M. A.
Teoria Quântica, São Paulo, 1988
(5) Christoffersen, R. E. Basic Principles and Techniques of Molecular Quantum
Mechanics
, Springer-Verlag, 1989
()
Hameka, H. F. Advanced Quantum Chemistry, Addison-Wesley, 1965
(6) MacDonald, J. K. L.
Phys. Rev. 43, 830 (1933)
(7) Löwdin, P. O.
Phys. Rev. 110, 1466 (1958)
(8) You, S. K., et. al.
Eur. J. Phys., 19, 179 (1998)
(9) Mur-Petit, J. et al.
Am. J. Phys., 70, 808 (2002)
(10) Rath, B.
Eur. J. Phys., 11, 184 (1990)
(11) Garrod, C. et. al.
J. Math. Phys. 16, 868 (1975)
Capítulo 5: Teoria do Campo Autoconsistente (SCF)
5.1 Introdução
A teoria do campo autoconsistente (self-consistent field, SCF) para sistemas
atômicos e moleculares, na formulação de Slater (1), se baseia no modelo seguinte: em vez
de considerar todas as interações instantâneas entre os pares de elétrons do sistema, supõe-
se que o potencial atuando sobre cada elétron é devido 1) aos núcleos e 2) à distribuição de
carga média dos outros elétrons.
Para os átomos com camadas fechadas esta distribuição é esférica, e o modelo é de
campo central (1/r), logo um campo atuando sobre um elétron em
i
r
G
, depende então da
distribuição de carga média dos outros elétrons, que tem a seguinte forma:
NN
ij
N
irepulsão
dd
r
rV
ττ
∫
ΨΨ= ...
1
)(
)2,...,2,1(
)2,...,2,1(
*
G
[1]
ou seja, que está definido em termos da autofunção que é justamente a variável. O
problema se resolve utilizando um método iterativo: postula-se uma certa distribuição
eletrônica que permite calcular o potencial aproximado, constroe-se então o Hamiltoniano
aproximado e acham-se as suas soluções. A nova Ψ é supostamente melhor que a primeira;
o campo construído a partir dela dá origem a um novo H cujas soluções dão origem à um
novo campo, e o processo se repete até que o campo seja autoconsistente, ou seja, até que o
potencial calculado em duas interações consecutivas seja constante, dentro da margem de
precisão pré estabelecida.
A teoria SCF, nesta formulação conduz às equações de Hartree-Fock-Slater, que são
a base de métodos usadas especialmente em física do estado sólido, e popularizadas em
química.
Existem porém várias teorias SCF: a que os químicos usam geralmente é a de
Hatree-Fock-Roothaan (HFR) (2), que é uma forma aproximada da teoria de Hartree-Fock
(HF). Ambas as teorias se baseiam num modelo puramente matemático: para o caso de
sistemas de camadas fechadas, consiste em representar o estado do sistema por uma função
Ψ escrita como um único determinante de Slater em termos de um conjunto de funções de
base monoeletrônicas
{}
i
Ψ
, sendo que as funções
i
Ψ
são otimizadas utilizando o princípio
variacional. No modelo de HF, as
i
Ψ
são variadas sem restrições, as equações são
resolvidas numericamente, e as soluções são os orbitais de Hartree-Fock,
)(HF
i
Ψ
: a
energia calculada formado um determinante de Slater com estes orbitais, é a energia
mínima que é possível conseguir com uma função de onda deste tipo.
Na prática, as variações efetuadas em
{
}
i
Ψ
estão em vários sentidos, e os orbitais
)(HF
i
Ψ
não são estritamente calculados. No modelo de HFR, as equações não são
resolvidas numericamente, e os orbitais são expandidos como combinações lineares de
funções de certos conjuntos finitos,
{
}
µ
Ψ
, que podem ser funções gaussianas, orbitais
atômicas, etc.:
∑
Ψ=Ψ
m
iHFR
i
c
µ
µµ
)()(
[2]
onde m é a dimensão do conjunto de base
{
}
µ
Ψ
. Os orbitais
)(HFR
i
Ψ
são aproximações aos
)(HF
i
Ψ
:
)()(
lim
HF
i
HFR
i
m
ite
Ψ=Ψ
∞→
[3]
No caso de moléculas, os orbitais moleculares (MO) são expandidos como combinação
linear de orbitais atômicos (LCAO), e o modelo chama-se de SCF-LCAO-MO. A
interpretação física das equações resultantes do modelo matemático acima não é tão
simples quanto a do modelo de Slater: além da interação de Coulomb entre os elétrons, o
potencial contém termos, chamados de troca ou de permuta (exchange) que surgem devido
à antissimetria da função de onda.
5.2 Valor médio da Energia
Dada uma função de onda
Total
Ψ
eletrônica para um sistema de 2N elétrons e M
núcleos, escrita como determinante de Slater em termos de um conjunto de 2N spin-
orbitais,
{
}
Spinespacial
Total
 Ψ×Ψ=Ψ [4]
Onde  é o operador antissimetrizador. A expressão para o valor médio da energia é:
τ
τ
d
dH
E
∫
∫
ΨΨ
ΨΨ
=
*
*
[5]
Sabendo que a função de base é antissimétrica e o operador Hamiltoniano
∑∑∑∑
<
+
−∇−=
N
i
N
i
N
j
ij
M
A
Ai
A
rr
z
H
222
2
1
2
1
[6]
lembrando que este é um Hamiltoniano especial (a interação de um elétron com o restante
do caroço), a Primeira parte entre colchete está constituída por operador de um elétron, os
termos de energia cinética dos elétrons e a porção da energia potencial relativa à interação
dos elétrons com os núcleos:
∑
=
N
i
ihH
2
1
)( [7]
onde
∑
−∇−=
M
A
Ai
A
i
r
z
h
2
2
1
)1(
[8]
Freqüentemente se utilizam aproximações nas quais os elétrons das camadas
internas não são considerados individualmente; o efeito dos elétrons internos é introduzido
no Hamiltoniano h(1) o qual contém a energia potencial do caroço molecular:
∑
=
N
i
caroço
ihH
2
1
)( [9]
a segunda parte do Hamiltoniano:
∑∑
<
=
N
i
N
j
ij
r
H
22
2
1
[10]
e depende de pares de elétrons.
Para calcular <E> é conveniente escrever o operador antissimetrizador  em termos
do operador permutação,
P
ˆ
(umas vez que as funções de base são antissimétricas):
∑
−=
P
P
P
N
Â
ˆ
)1(
!2
1
[11]
O fator (-1)
P
vale 1± segundo que a permutação seja par ou ímpar, e a somatória é sobre
todas as possíveis permutações de 2N elementos, ou seja, que tem 2N! termos. O operador
P
ˆ
permuta as coordenadas dos elétrons na função sobre a qual atua, assim, por exemplo:
{
}
{}
)2()2()4()4()3()3()1()1()1(
)4()4()3()3()2()2()1()1(
ˆ
2211
2211342
βαβα
βαβα
ΨΨΨΨ−=
=ΨΨΨΨ
P
P
[12]
para determinar se a permutação é par ou ímpar, compara-se os subíndices, neste caso 342,
com a ordem crescente normal, 234, e obtém-se o número de transposições necessárias para
transformar uma seqüência na outra. No caso considerado, transpondo o número dois em
duas casas para a esquerda obtemos a seqüência 234; como são necessárias duas
transposições, a permutação é par, e (-1)
P
= (-1)
2
= +1.
Como a função Ψ [4] está normalizada, o valor médio <E> pode ser escrito como:
{}
[]
{}
N
N
N
PP
PP
dddNNP
HHNNP
N
E
221
1
21
**
1
'
'
...)2()2()...1()1('
ˆ
)2()2()...1()1(
ˆ
)1(
!2
1
τττβα
βα
ΨΨ
+ΨΨ−=
∫∫
∑∑
+
[13]
5.3 O Valor médio de H
1
Consideramos primeiro o integral sobre H
1
, que chamamos de
1
H
. Como os
elétrons são indistinguíveis, e como Ψ é antissimétrica, as contribuições de todos os h(1)
são iguais; assim:
∫∫
∑
ΨΨ=Ψ
Ψ=
ττ
dihNdihH
N
i
)(2)(
*
2
*
1
[14]
Substituindo Ψ em [14]:
{
}
{}
NN
P
N
P
PP
ddddNNPh
NNP
N
H
2321
*
211
**
2
*
1
*
1
'
'
1
...)2()2()...3()3()2()2()1()1(')1(
)2()2()...3()3()2()2()1()1(
ˆ
...)1(
)!12(
1
ττττβαβα
βαβα
ΨΨΨΨ
ΨΨΨΨ−
−
=
∑
∫∫
∑
+
[
15]
a expressão [15] é uma soma de
2
)1!2)(!2(
−
NN
integrais múltiplas. É preciso analisar
separadamente dois casos:
'
ˆˆ
)1( PP ≠
Suponhamos por exemplo que
P
ˆ
é a identidade, e
12
ˆ
'
ˆ
PP = , ou seja que
'
ˆ
P
permuta
os elétrons 1 e 2. A integral múltipla correspondente é:
N
N
N
N
N
N
dNNNNddh
dddNNhNN
2
*
2
1
*
11
1
*
1
221
11
**
1
*
1
)2()2()2()2(...)2()2()2()2()1()1()1()1()1(
...)2()2()...1()1()2()2()1()2()2()...2()2()1()1(...
τββταβτβα
τττββαββα
ΨΨΨΨΨΨ=
ΨΨΨΨΨΨ
∫
∫
∫∫
∫∫
[16]
O operador h(1) só depende das coordenadas espaciais do elétron 1, e não das suas
coordenadas de spin; assim, a primeira integral do membro da direita da equação [16] pode
ser fatorizada ainda em:
1
*
11
*
1
)1()1()1()1()1(
ξβα
ddvh
∫
∫
ΨΨ
e, como as funções de spin são ortogonais, a segunda integral é zero. O mesmo acontece
com a integral sobre elétron 2. quanto às outras 2N-2 integrais, são todas de normalização,
ou seja que são iguais a 1.
Um caso levemente diferente ocorre se permutamos os elétrons 1 e 3 por exemplo.
Temos:
N
N
N
N
N
N
dNNNNd
dvdddh
ddddNN
hNN
2
*
3
31
*
22
1
*
1
1
1
2
*
1
23212
11
**
2
*
1
*
1
)2()2()2()2(...)3()3(
)3()3()2()2()2()2()1()1()1()1()1(
...)2()2()...1()1()2()2()3()3(
)1()2()2()...3()3()2()2()1()1(...
τββαα
τββαατ
ττττβαβα
βαβα
ξ
ξ
ΨΨ
ΨΨΨΨΨΨ=
ΨΨΨΨ
ΨΨΨΨ
∫∫∫
∫∫∫∫
∫∫
Observamos que mesmo se as integrais de spin não são zero, mas a integral
∫
ΨΨ
31
*
2
)3()3( dv
se anula pois
1
Ψ e
2
Ψ são ortogonais.
Assim, se
'
ˆˆ
P
P
≠
, todas as integrais sobre o operador h(1) (ou, em geral, sobre
qualquer operador de um elétron) são iguais a zero.
'
ˆˆ
)2( PP =
Obtém-se 2N! integrais com P
ˆ
e '
ˆ
P iguais; para todas elas, o fator
'
)1(
PP+
− é
1)1(
2
+=−
P
e diferença reside no spin-orbital associado ao elétron 1. são (2N-1)! Termos
para os quais 1 está associado a um certo spin-orbital
α
i
Ψ
ou
β
i
Ψ
:
∫
∫∫∫
ΨΨ=
ΨΨΨΨ
1
*
2
*
1
*
)1()1()1(
)2()2()2()2(...)1({)1()1({
dvh
dNNNNdh
ii
N
N
N
i
i
τββτ
α
β
α
β
pois é possível fazer (2N-1)! Permutações dos 2N-1 elétrons, deixando o elétron 1
associado ao spin orbital
α
β
{
i
Ψ
.
∫
∑
ΨΨ−
−
=
1
*
1
)1()1()1(2)!12(
)!12(
1
dvhN
N
H
ii
N
i
ii
N
i
hH
∑
= 2
1
[17]
o fator 2 aparece devido a que o mesmo orbital
i
Ψ
está associado às duas funções de spin
α
e
β
.
5.4 O Valor médio de H
2
O valor médio de
2
H
pode ser calculado de maneira análoga. Neste caso, para a
função Ψ antissimétrica, a dupla soma:
∑∑
<
N
ij
ij
r
2
1
pode ser substituída por
12
1
2
)12(2
r
NN
⋅
−
e
{
}
{}
N
N
N
PP
PP
dddNNP
r
NNP
N
H
221
1
12
**
1
'
'
2
...)2()2()...1()1('
ˆ
1
)2()2()...1()1(
ˆ
...)1(
)!22(2
1
τττβα
βα
ΨΨ
ΨΨ−
−
=
∫∫
∑∑
+
[18]
Consideremos separadamente os três casos seguintes '
ˆˆ
,'
ˆˆ
PPPP ≠= em mais de
dois elétrons, e
'
ˆˆ
P
P
≠
em dois elétrons somente.
'
ˆˆ
)1( PP =
Quando
'
ˆˆ
P
P
=
, podemos fazer (2N-2)! Permutações nas quais os elétrons 1 e 2
estão associados aos mesmos spin-orbitais
α
β
{
i
Ψ
e
α
β
{
j
Ψ
. Se os dois orbitais
i
Ψ
e
j
Ψ
são diferentes
)( ji ≠
, estes podem estar associados a quatro possíveis combinações das
funções de spin:
)2()1(
α
α
;
)2()1(
β
β
;
)2()1(
β
α
;
)2()1(
α
β
, e para cada integral
21
12
**
ij
)2()1(
1
)2()1(
ττ
dd
r
ji
ji
ΨΨΨΨ≡
∫
J
[19]
Há 4 (2N-2)! Termos. O integral
ij
J é chamada de integral de Coulomb. Portanto, se i=j,
ou seja, se as funções espaciais do dois elétrons são as mesmas, as funções de spin devem
ser diferentes, sendo que as duas possibilidades são:
)2()1(
β
α
;
)2()1(
α
β
.
21
12
**
ii
)2()1(
1
)2()1(
ττ
dd
r
ii
ii
ΨΨΨΨ≡
∫
J
[20]
(2)
'
ˆˆ
PP ≠ em mais de dois elétrons
Se
P
ˆ
difere de
'
ˆ
P
em mais de dois elétrons, é fácil demonstrar que as integrais valem zero.
(3)
'
ˆˆ
P
P
≠
em dois elétrons
Se
P
ˆ
difere de
'
ˆ
P
em só dois elétrons, porém as funções têm apenas dois spin-orbitais
invertidos (ou seja duas fileiras no determinante de Slater de uma das ψ). Novamente é
preciso considerar vários casos:
(
i) se ambos spin-orbitais tem a mesma parte espacial,
ji
Ψ
=
Ψ
, a integral é do tipo
∫
ΨΨΨΨ
21
12
**
)1()2()1()1(
1
)2()2()1()1(
τταββα
dd
r
iiii
que é igual a zero porque as funções de spin são ortogonais.
(
ii) se as partes espaciais são diferentes, e os spin também,
∫
ΨΨΨΨ
21
12
**
)2()2()1()1(
1
)2()2()1()1(
τταββα
dd
r
ijji
a qual também é igual a zero por causa da ortogonalidade das funções de spin.
(
iii) se as partes espaciais são diferentes mas as de spin são iguais, a integral é, por exemplo
∫
ΨΨΨΨ
21
12
**
)2()2()1()1(
1
)2()2()1()1(
τταααα
dd
r
ijji
[21]
que se reduz a
∫
ΨΨΨΨ≡
21
12
**
)2()1(
1
)2()2()1(
ττα
dd
r
ijjiij
K
[22]
a integral
ij
K se denomina integral de troca (exchange). Se aos orbitais
i
Ψ
e
j
Ψ da
equação [22] se associam funções de spin β, o resultado será o mesmo. Ainda para cada
para de orbitais
),(
ji
ΨΨ há 2(2N-2)! permutações dos outros spin-orbitais. Assim, a
contribuição de cada par ),(
ji
ΨΨ à energia total devido às integrais de troca é de
ij
N K)!22(2 − . Finalmente, voltando à equação [18] temos:
∑∑∑
≠
−+=
N
i
N
i
N
ij
ijijii
H )2(
2
KJJ
[23]
Os
ij
K
aparecem negativos porque correspondem à uma permutação ímpar.
5.5 Cálculo da Energia Média
Juntando os dois valores médios [17] e [23] temos:
∑∑∑∑
≠
−+=
N
i
N
i
N
ij
ijijii
N
i
ii
hE )2(2 KJJ
[24]
esta expressão pode ser simplificada observando que
∫
ΨΨΨΨ==
21
12
**
)2()1(
1
)2()1(
ττ
dd
r
iiiiiiii
KJ
de modo que
∑∑∑
−+=
N
i
N
j
ijij
N
i
ii
hE )2(2 KJ
[25]
Cada termo da primeira somatória na equação [25], 2h
ii
, corresponde à integração entre um
elétron no orbital
i
Ψ
e o núcleo; o fator 2 é devido a que dois elétrons ocupam cada orbital
i
Ψ . Os termos
∑
−
N
j
ijij
)2( KJ
representam a repulsão entre um elétron no orbital
i
Ψ
e os outros elétrons. É interessante
definir um conjunto de energias orbitais E
i
para um elétron ocupando o orbital
i
Ψ
∑
−+=
N
j
ijijiii
hE )2( KJ [26]
Se fizermos a soma destas orbitais sobre todos os elétrons do sistema, teremos:
∑∑∑∑
−+=
N
i
N
j
ijij
N
i
ii
N
i
i
hE )2(222 KJ
[27]
comparando esta expressão co a da energia total <E> da equação [25], observamos que
∑∑∑
−+=
N
i
N
j
ijij
N
i
i
EE )2(2 KJ
[28]
o qual é facilmente justificado já que na energia total, a repulsão entre cada par de elétrons
em
i
Ψ
e
j
Ψ deve ser contada só uma vez.
5.6 As Equações de Hatree-Fock.
Na seção anterior vimos que, dada uma função de onda aproximada, escrita como
um determinante de Slater em termos de um conjunto de funções de um elétron
{}
i
Ψ , a
expressão geral para a energia do sistema é:
∑∑∑
−+=
N
i
N
j
ijij
N
i
ii
hE )2(22 KJ
podemos aplicar o método variacional para achar as melhores funções de base
i
Ψ
, isto é,
aquelas tais que se formamos Ψ como um único determinante de Slater na base das
i
Ψ , a
energia é mínima:
0
=
E
δ
A variação do conjunto
{}
i
Ψ
porém está sujeita à restrição de que continue sendo um
conjunto ortonormal
=
≠
=
jise
jise
S
ij
1
0
ou seja, que:
0
=
ij
S
δ
O método que permite variar
E mantendo a ortonormalidade da base é um método
matemático bem conhecido, o dos
multiplicadores de Lagrande. Formamos uma função G:
∑
∑
−=
ij
ijij
SEG
λ
onde os
ij
λ
são os multiplicadores de Lagrange, por enquanto, parâmetros indeterminados,
e procuramos o mínimo na função G:
∑
∑
=−=
ij
ijij
SEG 0
δλδδ
Substituindo
E por seu valor (equação [25])
∑∑∑∑∑
=−−+
ij
ijij
i
N
i
N
j
ijijii
Sh 0)2(2
δλδδδ
KJ
[29]
Os termos da primeira soma na equação anterior:
∫
∫∫
ΨΨ+ΨΨ=ΨΨ=
1
*
1
*
1
*
))1()(1()1()1()1()1(()1()1()1(
τδτδτδδ
dhdhdhh
iiiiiiii
podem ser reescritos como:
conjugadocomplexodhh
iiii
+ΨΨ=
∫
1
*
)1()1()1((
τδδ
Analogamente:
..)2()1(
1
))2()(1()2()1(
1
)2())1(
21
i
12
*
i21
i
12
*
i
ccdd
r
dd
r
j
j
j
jij
+ΨΨΨΨ+ΨΨΨΨ=
∫∫
ττδττδδ
**
(J
ou utilizando o
operador integral de Coulomb:
∫
ΨΨ≡
2
12
*
)2(
1
)2()1(
ˆ
τ
d
r
jjj
J [30]
..)1(
ˆ
))1(()1(
ˆ
))1(
1i
*
1
i
ji
ccdd
j
jij
+ΨΨ+ΨΨ=
∫
∫
τδτδδ
JJ(J
*
Para as integrais de troca
..)1()1(
ˆ
))1()1()1(
ˆ
))1(
1j1i
ccdd
jiijij
+ΨΨ+ΨΨ=
∫
∫
τδτδδ
K(K(K
**
e definimos o operador integral de troca, )1(
ˆ
j
K , como um operador que tem a seguinte
propriedade:
)1()2(
1
)2()1()1(
ˆ
2
12
*
jijij
d
r
Ψ
ΨΨ=Ψ
∫
τ
K
[31]
Finalmente, a variação das integrais de recobrimento é dada por
∫
+ΨΨ= ..)1())1((
1
*
ccdS
jiij
τδδ
Juntando todos os termos e pondo em evidência as diferenciais:
0..)1()1()1(
ˆ
)1(
ˆ
2(2)1()1(2)1((
1
*
=+
Ψ−Ψ−+ΨΨ
∑
∫
∑∑
ccdh
i
N
j
iijij
N
j
jii
τλδ
KJ
O fator 2 que aparece na frente da primeira soma sobre j é devido à que tanto J
ij
como J
ji
contribuem com um termo idêntico,e analogamente K
ij
e K
ji
. Como as variações das
funções e dos seus complexos conjugados são independentes, cada termo da soma em i
deve ser identicamente igual a zero. Assim:
0)1()1())1(
ˆ
)1(
ˆ
2(2)1(2 =Ψ−Ψ
−+
∑∑
j
N
j
ijij
N
j
j
h
λ
KJ , i=1,2,...N
substituindo, por conveniência,
ijij
E2
≡
λ
e dividindo tudo por dois, temos:
0)1()1())1(
ˆ
)1(
ˆ
2()1(
=Ψ=Ψ
−+
∑∑
j
N
j
ijij
N
j
j
Eh KJ
, i=1,2,...N [32]
A Expressão contida entre os colchetes se denomina operador de Fock,
F
ˆ
:
))1(
ˆ
)1(
ˆ
2()1(
ˆ
j
N
j
j
hF KJ −+=
∑
[33]
temos então um conjunto de equações diferenciais de um elétron:
)1()1(
ˆ
i
N
j
iji
EF Ψ=Ψ
∑
, i=1,2,...N [34]
Vemos imediatamente que se todos os
ij
E com
j
i
≠
fossem iguais a zero, a equação [34]
seria uma equação de autovalores (devemos diagonalizar
ij
E para adquirirmos
ii
E ). Isto
pode ser conseguido mudando-se o conjunto de base
{
}
i
Ψ
. Com efeito, os orbitais
i
Ψ
não
são únicos. Qualquer combinação linear,
∑
Ψ=Ψ
j
iiji
T'
também é autofunção do operador
F
ˆ
, se os coeficientes
ij
T (transformação unitária) são
tais que:
ij
k
kjik
δ
=
∑
TT
*
[35]
uma transformação que obedece o requisito [35] é uma transformação unitária. Em forma
diferencial é
1
=
+
TT [36]
onde
+
T
é a matriz adjunta de T, ou seja, a matriz transposta conjugada. A equação [34]
pode ser escrita em forma matricial:
Ψ
=
Ψ
E
F
ˆ
[37]
Multipliquemos ambos membros desta equação por T pela esquerda, e insiramos a matriz
1, ou, o que é equivalente, o produto
TT
+
:
Ψ
=
Ψ
++
TTTTTT
E
F
ˆ
Agrupando os fatores, temos:
Ψ=Ψ
++
)TTT)TT(T EF (
ˆ
ou
''''
ˆ
Ψ
=
Ψ
E
F
que é análoga à equação [37]. Em particular podemos escolher uma matriz
T tal que E’seja
diagonal. Neste caso, teremos uma série de equações de autovalores.
''''
ˆ
iiii
EF Ψ=Ψ
i=1,2,...N [38]
para o operador
'
ˆ
F
. É possível demonstrar que o operador
'
ˆ
F
na base
{}
i
'Ψ
é exatamente
análogo ao operador
F
ˆ
na base
{}
i
Ψ
de modo que as equações [38] representam o mesmo
problema que as equações [34]. As equações [38] são as equações de Hartree-Fock, e os
orbitais
{}
i
'Ψ são os orbitais canônicos de Hartree-Fock.
Podemos escrever a equação [38] diretamente
iii
EF Ψ=Ψ'
ˆ
i=1,2,...N
com
)
ˆˆ
2()1(
ˆ
j
j
j
hF KJ −+=
∑
Notamos o Valor E
i
:
iiiii
FdE =ΨΨ=
∫
τ
F
ˆ
)
ˆˆ
2()
ˆˆ
2()1(
ij
j
ijiiij
j
ji
hdh KJKJ −+=Ψ
−+Ψ=
∑
∫
∑
τ
é exatamente igual a energia total E
i
definida pela equação [26]. Estas são chamadas de
potencial de Koopmans ou potenciais verticais.
A resolução das equações de Hartree-Fock se faz da maneira seguinte. Como o
operador de Fock,
F
ˆ
, contém os operadores integrais
j
J
ˆ
e
j
K
ˆ
definidos em termos de um
conjunto de orbitais de base
{}
i
Ψ
, inicia-se o processo postulando uma forma para o
conjunto
{}
i
Ψ
: por exemplo, os
)0(
i
Ψ
podem ser autofunções do operador )1(
ˆ
h
)0()0()0(
)1()1(
ˆ
iii
Eh Ψ=Ψ
com essas funções calcula-se o operador de Fock e resolvem-se as equações de Hartree-
Fock. As soluções
{
}
)1(
i
Ψ
formam um novo conjunto de base no qual pode se definir
novamente o operador de Fock. O processo continua até que os orbitais permaneçam
invariantes sob novas iterações, dentro da margem de precisão requerida. Se os
)1(
)0(
i
Ψ
estão bem escolhidos geralmente o processo converge.
É interessante discutir o significado físico do operador de Fock: ele representa um
Hamiltoniano efetivo de um elétron tal que a energia potencial do campo no qual o elétron
se encontra é constituída das seguintes partes:
(1) A energia de interação com os núcleos
(2) A energia de repulsão com todos os elétrons de spin oposto ao do elétron
considerado.
∑
j
j
J
ˆ
(3) A energia de interação com todos os elétrons de igual spin ao do elétron
considerado:
(
)
∑
−
j
jj
KJ
ˆˆ
que é a interação (2) devido a que, para uma função antissimétrica, dois elétrons com o
mesmo spin não podem ocupar o mesmo orbital e conseqüentemente nunca estão muito
perto um do outro; assim a sua energia de repulsão é menor
5.7 Referência:
(1) Slater, J. C.
Phys. Rev. (1953)
(2) Roothaan, C. C. J.
Rev. Mod. Phys. 32, 239, 245. (1960)
(3) Teixeira, J. J. C. D.,
Química Quântica: Fundamentos e Métodos, Fundação Calouste
Gulbenkian,
1982
(4) Gil, V. M. S.,
Orbitais em Átomos e Moléculas, Fundação Calouste Gulbenkian, 1996
(5) Bunge, A.
Introdução a Química Quântica, Edgar Blücher
(6) Levine, I. N.
Quantum Chemistry, Prentice Hall, 1991
(7) Peixoto, E. M. A.
Teoria Quântica, São Paulo, 1988
(8) Christoffersen, R. E.
Basic Principles and Techniques of Molecular Quantum
Mechanics
, Springer-Verlag, 1989
(9) Bazley, N. W.
Phys. Rev., 124, 483 (1961)
(10) Hameka, H. F.
Advanced Quantum Chemistry, Addison-Wesley, 1965
Capítulo 6: Método de Hartree-Fock-Roothaan.
6.1 Introdução
Para sistemas com muitos elétrons a resolução das equações de Hatree-Fock pode
ser feita de maneira aproximada pelo método de Hatree-Fock-Roothaan no qual os orbitais
de Hartree-Fock
{}
i
Ψ
, são aproximados como combinações lineares de funções de base
{
}
µ
χ
:
∑
=Ψ
µ
µµ
χ
)(i
i
c
[1]
para moléculas, os
i
Ψ
são orbitais moleculares e os
µ
χ
são orbitais atômicos. O método
variacional é aplicado aos parâmetros lineares
)(i
c
µ
.
Calculemos primeiro as integrais
ii
h
,
ij
J e
ij
K da expressão abaixo, em termos das
funções
{
}
µ
χ
∑∑∑∑
−+=
N
i
N
j
ijij
N
i
ii
N
i
i
hE )2(222 KJ
[2]
sendo assim obtemos,
∑
∫
∑∑
∫
==ΨΨ=
µν
µννµ
µν
νµνµ
τχχτ
hccdhccdhh
iiii
iiii
)(*)(
1
*
)(*)(
1
*
)1(
ˆ
)1()1(
ˆ
)1(
[3]
∑
∑∑∑∑
∫∫
=
==ΨΨΨΨ=
λσµ
σ
λ
µ
µν λσ
σλµσ
λ
µ
σµλ
ττχχχχ
v
jij
v
i
v
jjj
v
i
ji
jiij
vcccc
dd
r
cccc
r
)|(
)2()1(
1
)2()1()2()1(
1
)2()1(
)()(*)(*)(
21
12
**
)(*)(*)(*)(
12
**
J
[4]
onde o símbolo (µλ|νσ) se refere à integral sobre
12
1
r
tal que as funções da esquerda estão
associadas ao elétron 1 e as da direita, ao elétron 2. Analogamente:
∑
=
λσµ
σ
λ
µ
λµσ
v
ijj
v
i
ij
vccccK )|(
)()(*)(*)(
[5]
Substituindo na equação abaixo
∑∑∑∑
−+=
N
i
N
j
ijij
N
i
ii
N
i
i
hE )2(222 KJ
temos:
[]
∑∑ ∑∑∑
−+=
N
i
N
jv
jij
v
i
N
iv
v
i
v
i
vvcccchccE
λσµ
σ
λ
µ
µ
µµ
λµσσµλ
)|()|(22
)()(*)(*)()(*)(
[6]
onde é conveniente definir:
∑
≡
N
i
i
v
i
v
ccP
)(*)(
2
µµ
Substituindo na expressão de E:
∑∑
−+=
λσµ
σµλ
µ
µµ
λµσσµλ
v
v
v
vv
vvPPhPE )|(
2
1
)|(
2
1
procurando agora os melhores valores das constantes
)(i
c
µ
, variando a energia com respeito
a essas constantes. O procedimento e inteiramente análogo ao da secção anterior, e o
resultado e,
[]
∑∑∑∑
=−+
v
j
N
j
ij
N
j
j
i
j
v
v
v
i
v
Scvvccchc
µννσ
λ
µ
ελµσσµλ
)|()|(2
)(
)(
*)()(
onde ,...2,1
=
µ
fazendo uma transforma~c”ao unitária da base para diagonalizar a matriz
ε, e definido
[]
∑
−+≡
νσ
νσµµ
λµσσµλ
)|()|(2 vvPhF
vv
obtemos
0)(
)(
=−
∑
i
viv
cSF
ν
ν
µµ
ε
[7]
que são as equações de
Roothaan. Em notação matricial, torna-se
FC = SCε [8]
A grande vantagem pratica destas equações e que elas são simplesmente equações
algébricas e não mais diferencias.
Capítulo 7: Teoria da Perturbação.
7.1 Introdução
Agora nós iremos discutir um outro método de aproximação em Química Quântica,
a
teoria da Perturbação ou Teoria da Perturbação de Rayleigh-Schrödinger.
Consideramos um sistema o qual a equação de Schrödinger é
iii
EH Ψ=Ψ
[1]
e para um estado estacionário
)0()0()0(
)0(
iii
EH Ψ=Ψ
[2]
onde
)0(
H
é Hamiltoniano de um sistema não perturbado, ou seja um sistema na ausência
de radiação eletromagnética ou correlação interna. Quando há presença de uma perturbação
o Hamiltoniano que descreve a perturbação é representado por
H’. O Hamiltoniano
Perturbado
H e um não perturbado H
(0)
é
)0(
'
H
H
H
−
≡
[3]
'
)0(
H
H
H
+
=
[4]
onde
'
)0(
H
H
>
. Quando dizemos que
'
)0(
H
H
>
, queremos dizer que '
)0(
HH > ,
isto é, estamos comparando os valores esperados dos dois operadores
)0(
H
e '
H
.
O objetivo é relacionar as autofunções e autovalores do sistema perturbado às
autofunções e autovalores conhecidos do sistema não perturbado. Para isto, imaginamos
que a perturbação é aplicada de maneira a variar continuamente o sistema do estado não
perturbado ao estado perturbado. Matematicamente, isto corresponde à introdução de um
parâmetro
λ
no Hamiltoniano, modo que
'
)0(
HHH
λ
+=
[5]
com
λ
variando entre 0 e 1. Quando
λ
=0, o sistema é não perturbado e quando
λ
=1 o
sistema é perturbado. Assim temos:
iii
EHH Ψ=Ψ+ )'(
)0(
λ
[6]
desde que o Hamiltoniano depende do parâmetro
λ
, tanto a autofunção
i
Ψ
, como o
autovalor
i
E
, devem depender também de
λ
e podemos expandi-lo em séries de Taylor de
potencias de
λ
.
)()2(
2
)1()0(
...
k
i
k
iii
i
Ψ++Ψ+Ψ+Ψ=Ψ
λλλ
[7]
)()2(
2
)1()0(
...
k
i
k
iii
i
EEEEE
λλλ
++++= [8]
Onde
k é a ordem da função e energia. Substituindo [7] e [8] em [6], temos:
)...)(...(
)...)('(
)()2(
2
)1()0()()2(
2
)1()0(
)()2(
2
)1()0(
)0(
k
i
k
iii
k
i
k
iii
k
i
k
iii
EEEE
HH
Ψ++Ψ+Ψ+Ψ++++
=Ψ++Ψ+Ψ+Ψ+
λλλλλλ
λλλλ
[9]
...()(
)...'()'(
)2()0()1()1()0()2(
2
)1()0()0()1()0()0(
)1()2(
)0(2
)1(
)0(
)0()0(
)0(
+Ψ+Ψ+Ψ+Ψ+Ψ+Ψ
=Ψ+Ψ+Ψ+Ψ+Ψ
iiiiiiiiiiii
iiiii
EEEEEE
HHHHH
λλ
λλ
[10]
Igualando as famílias de mesma potência
)0()0()0(
)0(
iii
EH Ψ=Ψ
[11]
)()'(
)1()0()0()1()1(
)0(
)0(
iiiiii
EEHH Ψ+Ψ=Ψ+Ψ
λλ
0'
)1()0()0()1()1(
)0(
)0(
=Ψ+Ψ=Ψ+Ψ
iiiiii
EEHH
[12]
)2()0()1()1()0()2(
2
)1()2(
)0(2
()'(
iiiiiiii
EEEHH Ψ+Ψ+Ψ=Ψ+Ψ
λλ
0'
)2()0()1()1()0()2()1()2(
)0(
=Ψ+Ψ+Ψ=Ψ+Ψ
iiiiiiii
EEEHH
[13]
Agora devemos agrupar as expressões [11], [12] e [13]
(
)
(
)
0......)()'()(
2
)1()0(
)0(
)0()1()0()0()0(
)0(
=++Ψ−+Ψ−+Ψ−Ψ
λλ
iiiiiii
EHEHEH
[15]
se
i
Ψ
e
i
E
são funções contínuas de
λ
, para que a equação [15] seja verdadeira para
qualquer valor de
λ
, os coeficientes de cada potência de
λ
devem ser iguais a zero. As
equações [11] e [12] são equações de ordem zero e ordem um respectivamente. E como as
soluções de ordem zero são conhecidas elas podem ser utilizadas se obter as equações de
ordem um.
7.2 Correção de ordem um para a energia
Para obter a energia reescrevemos a equação [12]
)0()1()1()0(
)0(
)0(
)('
iiiii
EEHH Ψ=Ψ−+Ψ
multiplicando por
*)0(
i
Ψ
ambos os termos e integrando
)0()1(*)0()1()0(
)0(
*)0()0(*)0(
)('
iiiiiiii
EEHH ΨΨ=Ψ−Ψ+ΨΨ
)1()0()0()1()0(
)0(
)0(
'
iiiiii
EHEH =ΨΨ+Ψ−Ψ
O primeiro termo é zero porque
H é um operador Hermitiano. Então, a correção de ordem
um para a energia do nível
i é:
ii
iii
HdHE '
)0(
)1(
*)0()1(
≡ΨΨ=
∫
τ
[16]
ou
ii
iii
HHE '
)0(
)1(
)0()1(
≡ΨΨ=
.
,
)1(
)1
)0(
1
)(
)0(
)2(
)0()1(
)1(
)0()2(
−
=
ΨΨ=
ΨΨ+ΨΨ=
∑
m
i
m
i
n
m
n
i
iiiii
HE
HHE
#
[17]
7.3 Correção de ordem um para a função de onda
Para obtermos a correção de ordem um para uma função de onda, expandimos
)1(
i
Ψ
em um conjunto ortonormal completo de funções de ordem zero,
{
}
)0(
j
Ψ
:
∑
Ψ=Ψ
j
j
ji
i
a
)0()1(
, onde
)1()0(
ij
j
a ΨΨ= [18]
Substituindo [18] em [16] temos:
)0()1()0()0()0(
)0(
)'()(
iiji
j
jji
HEEHa Ψ−=Ψ−Ψ
∑
ou
)0()1()0()0()0(
)'()(
iiji
j
jji
HEEEa Ψ−=Ψ−
∑
[19]
multiplicando a equação [19] por
)*0(
k
Ψ
e integrando, temos:
)0()0()0()0()1()0()0()0()0(
')(
i
k
i
k
ij
k
i
j
jji
HEEEa ΨΨ−ΨΨ=ΨΨ−
∑
[20]
Efetuando a somatória, somente os termos com k = j permanecem, devido à
ortonormalidade:
)0()0()0()0()1()0()0(
')(
i
k
i
k
ii
k
ki
HEEEa ΨΨ−ΨΨ=− [21]
temos agora dois casos a considerar: para k=i, obtém-se a equação [16].
)0()0()1(
0
)0()0(
')(
i
k
ii
k
ki
HEEEa ΨΨ−=−
)0()0()0()0()1(
''
iii
k
i
HHE ΨΨ=ΨΨ=
[16]
Para
ik ≠
, o primeiro termo à direita da igualdade è zero e a equação [21] torna-se:
)0()0()0()0(
')(
i
k
i
k
ki
HEEa ΨΨ−=−
arranjando
)0()0(
)0()0(
'
k
i
i
k
ki
EE
H
a
−
ΨΨ
=
[22]
a equação [22] determina todos os coeficientes
ki
a , exceto
ii
a , o coeficiente de
)0(
i
Ψ .
Aplicando a condição de ortonormalização, se determina
0
=
ii
a
. Os resultados finais para
i
E
e
i
Ψ
corretos até a primeira ordem são os seguintes:
iiii
HEE '
)0()1(
λ
+=
[23]
∑
≠
Ψ
−
ΨΨ
+Ψ=Ψ
ik
i
k
i
i
k
ii
EE
H
)0(
)0()0(
)0()0(
)0()1(
'
λ
[24]
7.4 Tratamento de um sistema de spins pela Teoria da Perturbação
Quando estudamos um ensemble de spins através de perturbações, realizamos um
deslocamento de equilíbrio do sistema de spins e podemos fazer transições entre níveis de
energias que ocasionam mudanças na magnetização na magnetização líquida. Esta por sua
vez nos dá o sianl de RMN.
Considerando um sistema de dois níveis representados por
a
e
b
, com energias
Ea e Eb, sujeito a uma perturbação V(t)=Vf(t), onde V é um operador independente do
tempo e f(t) é um fator numérico flutuante. O operador V, tem com efeito a modulação nos
níveis de energia Ea e Eb, o qual induz transições entre os níveis.
aaba
EtfVEV )(=
bbab
EtfVEV )(=
Temos para esse sistema o Hamiltoniano
H
EVft
Vft E
aab
ba b
=
()
()
e a função de onda que é solução da equação de Schrödinger:
Ψ= −
+−
ct i
Et
act i
Et
b
a
a
b
b
()exp ()exp
==
Resolvendo o sistema abaixo:
i
t
ct i
Et
ct i
Et
EVft
Vft E
ct i
Et
ct i
Et
a
a
b
b
aab
ba b
a
a
b
b
=
=
=
=
=
∂
∂
()exp
()exp
()
()
()exp
()exp
−
−
=
−
−
temos:
i
t
c t Vc tft i
EEt
i
t
c t Vc tft i
EEt
aabb
ab
bbaa
ba
=
=
=
=
∂
∂
∂
∂
() () ()exp
()
() () ()exp
()
=
−
=
−
Supondo que no instante t=0 o sistema esteja no equilíbrio, ou seja, no seu estado
a
, em
que c
a
(0)=1 e c
b
(0)=0. A integração da equação acima fornece o valor da correção em c
b
num tempo t posterior:
dt
tEE
itf
V
itc
dt
tEE
itf
i
V
tc
t
abba
b
t
abba
b
∫
∫
−
−=
−
=
0
0
)(
exp)()(
)(
exp)()(
==
==
A probabilidade de que ocorra uma transição do estado
a
para o estado
b
no tempo t é
dada por:
Pab c t
b
(,) ()=
2
Um caso em que P(a,b) cresce linearmente com o tempo é aquele em que f(t) é uma
função periódica do tipo:
ft t() cos( )
=
2
ω
dt
tEE
it
V
itc
t
abba
b
∫
−
−=
0
)(
exp)cos(2)(
==
ω
escrevendo
()
ω
ba
Eb Ea
=
−
=
, (ou seja, =ω
ba
é a diferença de energia entre os níveis)
()
dttit
V
itc
t
ab
ba
b
∫
−=
0
exp)cos(2)(
ωω
=
e sendo:
cos( )
exp( ) exp( )
ω
ω
ω
t
it it
=
+
−
2
()
dttititi
V
itc
t
ba
ba
b
∫
−+−=
0
)exp()exp()exp()(
ωωω
=
()
dttititi
V
itc
t
ba
ba
b
∫
−+−=
0
)exp()exp()exp()(
ωωω
=
()
−++−=
∫
dttiti
V
itc
t
baba
ba
b
0
])(exp[])(exp[)(
ωωωω
=
−
−−
+
+
−+
−=
ωω
ωω
ωω
ωω
ba
ba
ba
baba
b
titiV
tc
1])(exp[1])(exp[
)(
=
Para ω próximo de ω
ba
podemos desprezar o primeiro termo da equação acima com
relação ao segundo termo.A probabilidade
2
)( b)P(a, = tc
b
fica então igual a:
−
−−−
⋅
−
−−
==
ωω
ωω
ωω
ωω
ba
ba
ba
ba
ba
b
titi
V
TcbaP
1])(exp[1])(exp[
)(),(
2
2
2
=
()
−
+−−−−−
=
22
2
1])(exp[])(exp[1
),(
ωω
ωωωω
ba
baba
ba
titi
V
baP
=
()
()
()
t
t
t
V
baP
t
V
baP
t
V
baP
t
V
baP
ba
ba
ba
ba
ba
ba
ba
ba
ba
ba
ba
ba
−
−
=
−
−
=
−
−−
=
−
−−
=
2
2
2
2
2
2
2
2
2
2
2
22
2
2
2
2
)(
sen
2
),(
2
1
2
)(
sen2
),(
2
1
)cos(1
),(
)cos(22
),(
ωω
π
ωω
π
ωω
ωω
ωω
ωω
ωω
ωω
=
=
=
=
P(a,b)
ω
ba
ω
probabilidade de transição em função da frequencia da perturbação (fig.3)
Na figura 2 temos uma simulação da função probabilidade de transição P(a,b) entre dois níveis
com energias Ea e Eb. Aproximando essa função por uma função delta de Dirac:
δ
()
;
;
xx
xx
xx
−=
=
≠
0
0
0
1
0
t
V
baP
ba
ba
)(
2
),(
2
2
ωωδ
π
−=
=
A primeira condição a ser satisfeita pela perturbação é que sua atuação forneça ao sistema
uma energia
=ω=hν, onde ν= (Eb-Ea)/h, ou seja, a energia fornecida ao sistema deve ser
igual à diferença de energia entre os níveis.
O termo independente do tempo da perturbação deve ser tal que V
ba
2
≠0.
VaVb
ba
2
2
0=≠
7.5 Referência:
(1) Bunge, A.
Introdução a Química Quântica, Edgar Blücher
(2) Levine, I. N.
Quantum Chemistry, Prentice Hall, 1991
(3) Peixoto, E. M. A.
Teoria Quântica, São Paulo, 1988
(4) Christoffersen, R. E.
Basic Principles and Techniques of Molecular Quantum
Mechanics
, Springer-Verlag, 1989
(5) Bazley, N. W.
Phys. Rev., 124, 483 (1961)
(6) Teixeira, J. J. C. D.,
Química Quântica: Fundamentos e Métodos, Fundação Calouste
Gulbenkian,
1982
(7) Hanna, M. W.
Quantum Mechanics in Chemistry, The Benjamin/Cummings
Publisingh
1969
(8) Hameka, H. F.
Advanced Quantum Chemistry, Addison-Wesley, 1965
Livros Grátis
( http://www.livrosgratis.com.br )
Milhares de Livros para Download:
Baixar livros de Administração
Baixar livros de Agronomia
Baixar livros de Arquitetura
Baixar livros de Artes
Baixar livros de Astronomia
Baixar livros de Biologia Geral
Baixar livros de Ciência da Computação
Baixar livros de Ciência da Informação
Baixar livros de Ciência Política
Baixar livros de Ciências da Saúde
Baixar livros de Comunicação
Baixar livros do Conselho Nacional de Educação - CNE
Baixar livros de Defesa civil
Baixar livros de Direito
Baixar livros de Direitos humanos
Baixar livros de Economia
Baixar livros de Economia Doméstica
Baixar livros de Educação
Baixar livros de Educação - Trânsito
Baixar livros de Educação Física
Baixar livros de Engenharia Aeroespacial
Baixar livros de Farmácia
Baixar livros de Filosofia
Baixar livros de Física
Baixar livros de Geociências
Baixar livros de Geografia
Baixar livros de História
Baixar livros de Línguas
Baixar livros de Literatura
Baixar livros de Literatura de Cordel
Baixar livros de Literatura Infantil
Baixar livros de Matemática
Baixar livros de Medicina
Baixar livros de Medicina Veterinária
Baixar livros de Meio Ambiente
Baixar livros de Meteorologia
Baixar Monografias e TCC
Baixar livros Multidisciplinar
Baixar livros de Música
Baixar livros de Psicologia
Baixar livros de Química
Baixar livros de Saúde Coletiva
Baixar livros de Serviço Social
Baixar livros de Sociologia
Baixar livros de Teologia
Baixar livros de Trabalho
Baixar livros de Turismo
