Práctica 1: Átomos de dos electrones
Fı́sica Atómica y Molecular
Curso 2018-19
1 Unidades
Compile el programa unidades.f y ejecútelo. Responda a las siguientes
preguntas:
• Anote los factores de conversión entre unidades atómicas, electron-
voltios, julios, y cm−1 .
• ¿En qué unidades está expresada la constante de estructura fina?
• ¿Cuál es el error relativo (en %) cometido al considerar la masa del
electrón en lugar de la masa reducida (en el átomo de He)?
• Estime el orden de magnitud del cociente entre la intensidad de la
interacción gravitatoria y la electrostática.
2 Método perturbativo
Compile el programa he1s2s2p.f y ejecútelo. El programa incluye la re-
pulsión electrónica 1/r12 usando teorı́a de perturbaciones a primer orden
como una corrección al Hamiltoniano en el que los dos electrones solamente
sienten el potencial culombiano nuclear. Responda a las siguientes pregun-
tas:
• Para el átomo de Helio (Z=2), ¿cuál es la energı́a de la configuración
electrónica 1s2s (términos 21 S y 23 S) para el átomo de He? Relacione
este valor con la energı́a sin incluir el término 1/r12 y las integrales Jnl
y Knl cuyo resultado se muestra en pantalla. ¿Cuál es la diferencia
relativa de energı́a (en %) entre los estados singlete y triplete?
• Repetir el apartado anterior para la configuración electrónica 1s2p
(términos 21 P y 23 P ). Notar que la repulsión electrónica produce un
desdoblamiento con respecto a los niveles sin perturbar (las integrales
Jnl y Knl dependen de l) ¿Cuál es la diferencia relativa de energı́a (en
%) entre ambas configuraciones electrónicas?
1
• Ejecute el programa de nuevo para el ión H − . ¿Cuál es la diferencia
relativa de energı́a (en %) entre ambas configuraciones electrónicas en
este caso? Compare con el caso del Helio y justifique el resultado.
• Ejecute el programa de nuevo para el ión C +4 . ¿Cuál es la diferencia
relativa de energı́a (en %) entre ambas configuraciones electrónicas en
este caso? Compare con el caso del Helio y justifique el resultado.
3 Cálculo variacional: estado fundamental (1)
Compile el programa dosexp.f y ejecútelo. El programa utiliza el método
variacional usando funciones de onda prueba del tipo ψ(r1 , r2 ) = [A exp(−αr1 )+
B exp(−βr1 )][A exp(−αr2 )+B exp(−βr2 )] que son una generalización de las
funciones de onda hidrogenoides, y las variables A1 , A2 , α1 , α2 son parámetros
variacionales. El algoritmo necesita una estimación inicial de dichos parámetros
para realizar la minimización, pero puede comprobar que ligeras variaciones
en dicha estimación dan lugar al mismo mı́nimo. Responda a las siguientes
preguntas:
• Ejecute el programa para el átomo de Helio y anote los valores óptimos
de los parámetros ası́ como la energı́a del estado fundamental tras
la minimización. Compare con el resultado del método perturbativo
(visto en teorı́a).
• Repita el cálculo para los iones H− y C+4 y anote las energı́as del estado
fundamental. Compare con el resultado del método perturbativo (visto
en teorı́a). ¿Se obtiene un estado ligado para el ión H− con esta familia
de funciones de onda prueba ?
4 Cálculo variacional: estado fundamental (2)
Compile el programa correlacionado.f y ejecútelo. El programa utiliza el
método variacional usando funciones de onda prueba del tipo ψ(r1 , r2 , r12 ) =
exp(−α(r1 +r2 ))[1+b exp(−βr12 )]. A diferencia del caso anterior, la función
de onda prueba depende de r12 con lo que los electrones en esta función de
onda están correlacionados. El algoritmo necesita una estimación inicial de
dichos parámetros para realizar la minimización, pero puede comprobar que
ligeras variaciones en dicha estimación dan lugar al mismo mı́nimo.
• Ejecute el programa para el átomo de Helio y anote los valores óptimos
de los parámetros ası́ como la energı́a del estado fundamental tras la
minimización. Compare con el resultado que se obtiene usando la
función de onda prueba del apartado anterior. ¿Se obtiene una mejora
al incluir la correlación ?
2
• Repita el cálculo para los iones H− y C+4 y anote las energı́as del
estado fundamental. Compare con el resultado que se obtiene usando
la función de onda prueba del apartado anterior. ¿Se obtiene una
mejora al incluir la correlación ?
• Mencione otra aplicación que pueda tener dicha función de onda prueba,
aparte de mejorar la cota superior de la energı́a del estado fundamen-
tal.
5 Cálculo variacional: estado fundamental (3)
Compile el programa diago.f y ejecútelo. El programa utiliza el método
variacional usando funciones de onda prueba del tipo ψ(r1 , r2 , r12 ) como
suma de monomios del tipo exp(−α1 r1 − α2 r2 )r1n1 r2n2 r12m . A diferencia del
caso anterior, la función de onda prueba es una suma de funciones cuyo
grado máximo está fijado. El algoritmo requiere especificar dicho grado
máximo, ası́ como una estimación del parámetro α y si queremos incluir la
dependencia en la variable r12 (cuando el programa pregunte, introducir 0 si
no se desea dependencia en r12 ). A diferencia de los programas anteriores,
el usuario tiene que indicar el intervalo de α donde se espera encontrar el
mı́nimo de energı́a, de forma que hay que especificar el salto en α y el número
de iteraciones. El usuario deberá ejecutar el programa repetidas veces hasta
encontrar un mı́nimo aceptable.
• Ejecute el programa para el átomo de Helio y use un valor de n máximo
de 1, y compare con el resultado del apartado anterior. Repita para
un valor de n máximo de 10, y compare con lo que se obtiene
• Con la potencia de cálculo que le permite su ordenador, compruebe
cuántos términos de este tipo puede incluir para que el cálculo dure
menos de 30 segundos. Compare el resultado con el valor exacto (in-
finitos términos) que viene en el libro.
6 Cálculo variacional: estado fundamental (4)
Compile el programa hlog.f y ejecútelo. El programa utiliza el método
variacional usando funciones de onda prueba de tipo Hylleraas:
max mmax
lX X/2 nX
max
ψ(s, t, u) = exp(−ks) (sl t2m un ) (1)
l=0 m=0 n=0
donde s = r1 + r2 , t = r1 − r2 , y u = r12 . Vemos que la función de
onda se expresa como suma de monomios. El algoritmo requiere que el
usuario especifique los valores máximos de los exponentes (l, m, n), teniendo
3
en cuenta que m ha de ser un número par (para que la función de onda
espacial del estado fundamental sea simétrica, es decir un estado para).
También se ha de especificar un valor para el parámetro variacional k.
• Ejecute el programa para el átomo de Helio y use un valor de n máximo
de 1, y compare con el resultado del apartado anterior. Repita para
un valor de n máximo de 10, y compare con lo obtenido en el apartado
anterior.
7 Cálculo variacional: estados excitados
Compile el programa excit-he.f y ejecútelo. El programa utiliza el método
variacional usando funciones de onda prueba del tipo ψ(r1 , r2 ) = [A1 exp(−α1 r1 )][exp(−α2 r2 ) Ai r2i ],
P
donde las variables α1 , α2 y los coeficientes polinómicos son parámetros
variacionales. Responda a las siguientes preguntas:
• Ejecute el programa para la configuración electrónica 1s2s. Compare
con el resultado obtenido con el método perturbativo. Razone por
qué se obtiene una energı́a mayor (menos ligada) con este método.
Mencione una interpretación fı́sica para los valores óptimos de α1 y
α2 .
• Según la página del NIST, el átomo de Helio presenta las siguientes
transiciones en el espectro visible: 3888.65Å, 5875.61Å, 6678.15Å,
7065.18Å. Deduzca, a partir de las energı́as de los distintos estados
excitados obtenidos con el programa, a qué tipo de transiciones corre-
sponden.