Potencial de Morse

Fazer download em pptx, pdf ou txt
Fazer download em pptx, pdf ou txt
Você está na página 1de 19

Química Computacional

Método de Euler aplicado para solução da equação de Schrodinger

Professor: Dr. Alejandro Lopez Castilho

Felipe Ferri Hilario RA 640077


Guilherme Cazagrande Felix RA 728599
Equação de Schrodinger unidimensional independente do tempo:

HΨ=EΨ
• H = Hamiltoniano do sistema, representa a
energia total do sistema, ele não é definido
pela equação e deve ser determinado pelas
ħ2 𝑑 2 propriedades físicas do sistema.
H=hamiltoniano do sistema 𝐻=− + 𝑉(𝑥)
2𝑚 𝑑𝑥 2 ħ2
E • termo constante aonde ħ é a constante de
2𝑚
planck e m a massa do elétron.
• 𝑉(𝑥) energia potencial do Sistema
• E
Manipulação algébrica:

ħ2 𝑑 2
HΨ(x) = − +𝑉 𝑥 Ψ = EΨ(x)
2𝑚 𝑑𝑥 2
• A equação de Shrodinger admite diversas soluções
ħ2 𝑑 2 Ψ • Em todos os casos os valores de energia potencial e
− + 𝑉 𝑥 Ψ=EΨ total devem ser tais que:
2𝑚 𝑑𝑥 2 ∞
න Ψ 𝑥 2 𝑑(𝑥) =1
ħ2 𝑑2Ψ −∞
− 2 =EΨ − 𝑉 𝑥 Ψ = (E−V(x))Ψ(x) • A equação acima representa a probabilidade de
2𝑚 𝑑𝑥
encontrar a partícula em um espaço infinito

ħ2 𝑑 2 Ψ
= 2 (V(x)−E) Ψ(x)
𝑚 𝑑𝑥 2
Método de Euler:

𝑑𝑦 Δ𝑦 𝑦𝑛+1 − 𝑦𝑛 𝑦𝑛+1 − 𝑦𝑛
= 𝐹(𝑥, 𝑦) = 𝐹(𝑥, 𝑦) = 𝐹(𝑥, 𝑦) = 𝐹(𝑥, 𝑦)
𝑑𝑥 Δ𝑥 𝑥𝑛+1 − 𝑥𝑛 ℎ
𝑥𝑛+1 − 𝑥𝑛 = ℎ

• O desenho gráfico da função de interesse é obtido calculando-se o


valor desta ponto a ponto
𝑌𝑛+1 = 𝐹 𝑥, 𝑦 ℎ + 𝑌𝑛 • h=passo dos cálculos
• Yn= valor de Y no n esimo passo
• Yn+1= valor da função no próximo passo
• F(x,y) = função a ser fornecida (determina o formato do gráfico),
também é calculado em cada passo
Aplicação do método de Euler para solução da equação de Shrodinger:

ħ2
ħ2 𝑑 2 Ψ • é aproximado para 1
𝑚
= 2 (V(x)−E) Ψ(x) • Ψ é a variável dependente “Y”
𝑚 𝑑𝑥 2
• 2 (V(x)−E) Y(x) = F(x,y)

dY dY
=𝑧 =𝑧
𝑑𝑥 𝑑𝑥
ħ2 𝑑 2 Ψ 𝑑2Y
= 2 (V(x)−E) Ψ(x) = F(x,y)
𝑑 2 Y 𝑑𝑧 𝑚 𝑑𝑥 2 𝑑𝑥 2 𝑑𝑧
= = F(x,y)
𝑑𝑥 2 𝑑𝑥 𝑑𝑥
Aplicação do método de Euler para solução da equação de Shrodinger:
dY
=𝑧 𝑌𝑛+1 = 𝑧𝑛 ℎ + 𝑌𝑛
ħ2 𝑑 2 Ψ 𝑑𝑥
= 2 (V(x)−E) Ψ(x)
𝑚 𝑑𝑥 2 𝑑𝑧
= F(x,y) 𝑧𝑛+1 = 𝐹 𝑥𝑛 , 𝑦𝑛 , 𝑧𝑛 ℎ + 𝑧𝑛
𝑑𝑥

Dessa forma o valor da função de onda pode ser


calculado ponto a ponto, de forma numérica, se
forem fornecidos determinados parâmetros.
• F(x,y)
• z
• Y0
• h:determinado a partir dos valores de xi e xf
Soluções da equação de Shrodinger unidimensional independente do tempo:

4 soluções estudadas: Cada solução fornece


• Oscilador harmônico valores diferentes de
• Potencial de Morse energia cinética, V(x), e
• Átomo de hidrogênio energia potencial
• Poço duplo

1
OH: 𝑉 𝑥 = 2 𝑘𝑥 2
Potencial de morse:
𝑧 𝑙(𝑙+1)
H: 𝑉 𝑥 = − 𝑟 + 𝑉 𝑥 = 𝐷𝑒(1 − 𝑒 −𝑎1 𝑥−𝑎0
)2
2𝑟 2

Poço Duplo: 𝑉 𝑥 = −𝑥 2 + 0,2𝑥 4


Potencial de morse:

𝑉 𝑥 = 𝐷𝑒(1 − 𝑒 −𝑎1 𝑥−𝑎0 )2

2
𝑑2Ψ 𝑑𝑧
ħ 2 = = 2m (V(x) − E) Ψ(x) = 2m (𝐷𝑒(1 − 𝑒 −𝑎 𝑟−𝑟𝑒
)2 −E) Ψ(x) = F(x,y)
𝑑𝑥 𝑑𝑥

• 𝑎1 = 1,02 𝑢. 𝑎 E: deve ser tal que:



• 𝑎0 = 1,4 𝑢. 𝑎 2 𝑑(𝑥)
න Ψ 𝑥 =1
• 𝐷𝑒 = 0,176 𝑢. 𝑎 −∞
• 𝑚 = 916,57
Primeira energia: E= 0.0098525553

𝑑2 Ψ 𝑑𝑧
2
= = −2m (𝐷𝑒(1 − 𝑒 −𝑎 𝑟−𝑟𝑒
)2 −E) Ψ(x) = F(x,y,z) 𝑧𝑛+1 = 𝐹. ℎ + 𝑧𝑛
𝑑𝑥 𝑑𝑥

𝑑Ψ
=𝑧 Ψ𝑛+1 = 𝑧ℎ + Ψ𝑛
𝑑𝑥
E0=0,00985255
E1=0,0287065
E2=0.0464235
E3=0.063027
E3=0.063027
E2=0.046423
E1=0,0287065
E0=0,00985255
Sentido físico:

V(r)

O potencial de Morse descreve de forma


simplificada a interação de um sistema diatômico.
Na curva de potencial ao lado r representa a
distância entre os núcleos dos átomos e V(r) a
energia potencial da ligação

Potencial de morse:

𝑉 𝑥 = 𝐷𝑒(1 − 𝑒 −𝑎 𝑟−𝑟0 )2
r
As energias encontradas (resultantes da equacão de
Shrodinger) são interpretadas como modos vibracionais.

A absorção de energia (por interação com radiação


por exemplo) excita a molécula a um modo
vibracional mais energético

Estados excitados

Estado fundamental
Ao retornar para o estado fundamental (ou
passar para um estado menos energético)
a molécula libera o excesso de energia na
forma de radiação.

ℎ𝑐
𝐸= A diferença de energia entre
ƛ
os níveis fornece o
comprimento de onda.
hv
Referências:

1- MORSE, P.M. Diatomic molecules according to the wave mechanics II. Vibrational levels. Phys. Ver. N 34, p.57-64, 1929

2- Daniel D. Knowalow, Joseph O. Hirschfelfer, “Morse Potential Parameters For O-O, N-N, and N-O
Interations”, The Physics of Fluids, 1961, volume 4

Você também pode gostar