INTRODUCCIÓN A LOS CÁLCULOS DE

ESTRUCTURA ELECTRÓNICA
Universidad de Ciencias Aplicadas y Ambientales
Quı́mica

Bermudez Abril Luisa Fernanda, Tovar Bahamon Angie Paola

25 de marzo de 2020

Resumen
La tabla periódica nos permite observar las propiedades de los elementos a lo largo de
los periodos y grupos. Se estudiaron las energı́as de ionización y afinidad electrónica para
los primeros 10 elementos de la tabla periódica, utilizando cálculos computacionales con la
ayuda de los programas Gabedit y Gamess empleando el modelo de Hartree-Fock. Para las
moléculas H2 , N2 y N H3 , se realizó la optimización de energı́a, con diferentes funciones base,
3-21G, 6-31G, 6-31G (d,p), 6-31++G (2df,2dp), y el tiempo de cálculo para cada una. Tam-
bien se realizó un cálculo de punto sencillo con la geometrı́a optimizada del agua, etileno y
benceno para observar los orbitales moleculares HOMO y LUMO.

Palabras claves: energı́a de ionización, afinidad electrónica, Gamess, Hartree-Fock

Abstract
The periodic table allows observing the properties of the elements through the periods and
groups. Ionization energies and electron affinity for the first 10 elements of the periodic ta-
ble were calculated using computational programs, Gabedit and Gamess, with Hartree-Fock
model. For H2 , N2 and N H3 molecules, energy optimization were performed, with different
basic sets, 3-21G, 6-31G, 6-31G (d, p), 6-31 ++ G (2df, 2dp), and the calculation time for
each one. Finally, a simple point calculation was performed with the optimized geometry of
water, ethylene and benzene to observe the HOMO and LUMO molecular orbitals.

Keywords: ionization energy, electron affinity, Gamess, Hartree-Fock

Introducción las y átomos.[1] Ciertas propiedades fı́sicas

de los elementos reflejna en una correspon-
La idoneidad de los métodos ab initio diente regularidad de comportamiento quı́mi-
para una determinación precisa del poten- co. Existen propiedades fı́sicas de naturaleza
cial de ionización y afinidad electrónica es electrónica como lo es la energı́a de ionización
objeto de análisis sistemático para molécu- y la afinidad electrónica [2].

1
 de φi [4].
La energı́a de ionización es uno de los paráme-
tros más importantes de un átomo, esta es el El teorema de Koopmans proporciona cierto
trabajo realizado para eliminar el electrón de apoyo para esto, al garantizar la estabilidad
un átomo en estado gaseoso y se mide conve- de φi en sı́. La ecuación 4 se beneficia de
nientemente en electronvoltios [3]. Para esta una cancelación parcial fortuita de errores:
energı́a se emplea el simbolo I. Su definición el descuido de Hartree-Fock de la correla-
formal serı́a la siguiente [2]: ción electrónica hace que E(X + i ) y E(X) sean
demasiado altos, pero más aún esto último
en el proceso porque X tiene más electrones. Sin embargo,
esto se compensa en cierta medida al no tener
A(g) → A+ + e− ; ∆E = I (1) en cuenta la relajación de X+ , la cual aumen-
i
ta E(X + i ). Ası́, la ecuación 4 podrı́a ser una
donde el electrón que sale es el menos ligado
aproximación razonable. Cuando está escri-
al átomo A, la energı́a desarrollada se define
to para el orbital de mayor energı́a en una
como energı́a de ionización, I [2].
molécula, se convierte en [4],
Ası́, llamando E+ a la energı́a total de un Ii ≈| εHOMO | (5)
◦
ion positivo y E a la del átomo neutro en es-
en el que Ii es la primera energı́a de ioniza-
tado gaseoso, la primera energı́a de ionización
ción y εHOMO se refiere a la energı́a del orbital
puede expresarse como [2]
molecular más ocupado (HOMO) [4].
I = E+ − E◦ (2)
La afinidad electrónica es un parámetro que
Este concepto es generalizable para siste- esta relacionado con la energı́a del proceso de
mas moleculares, donde A representaria una aceptación de un electrón por un átomo da-
molécula neutra. [2]. do. En este sentido, complementa al dato de
energı́a de ionización, pues con ambos se tie-
Dentro del marco de Hartree-Fock, la energı́a ne información sobre la facilidad relativa con
de ionización Ii de un electrón en un orbital la que un átomo acepta o pierde un electrón.
φi en un átomo o molécula X puede deter- Formalmente la afinidad electrónica se define
minarse fácilmente simplemente insertando como la diferencia de energı́as totales del áto-
las expresiones de energı́a de Hartree-Fock mo A y de su ion negativo A− en sus estados
apropiadas [4]: basales [2]:

Ii = E(X + − AEA = Etotal (A) − Etotal (A− ) (6)

i ) + E(e ) − E(X) (3)
la que corresponde a la energı́a liberada en la
donde E(X) y E(X + i ) son las energı́as de X y reacción de ionización de A− :
el ion positivo formado por la pérdida de un
electrón de φi , el i-ésimo orbital canónico de A− −
(g) → A(g) + e ; ∆E = AEA (7)
Hartree-Fock. Esto lleva a [4],
Las afinidades electrónicas son generalmente
Ii ≈| εi | (4) positivas. El nombre de esta propiedad se de-
be a la reacción inversa [2]
donde εi es la energı́a de un electrón en φi . La
e− + A(g) → A− (g) ; ∆E = −AEA (8)
ecuación 4 se basa una suposición principal:
que los orbitales ocupados de X no se ven generalmente es exotérmica e indica la mayor
afectados por la eliminación de un electrón o menor facilidad con la que A acepta un

2
electrón [2]. y Orbitales de tipo gaussiano (GTO) [6].

El teorema de Koopmans tambien permite Un procedimiento que se ha utilizado am-

mencionar la afinidad electrónica como [4]: pliamente es ajustar un orbital de tipo Slater
(STO) a una combinación lineal de n = 1,
A ≈ εLUMO (9)
2, 3,. . ., funciones Gaussianas primitivas.
donde A es la afinidad electrónica de la Este es el procedimiento STO-nG. En parti-
molécula y εLUMO es la energı́a del orbital cular, los conjuntos bases STO-3G a menudo
molecular no ocupado más bajo (LUMO) [4]. se usan en cálculos poliatómicos, en lugar de
evaluar integrales con funciones Slater [6].
Uno de los programas que permite estudiar
las propiedades atómicas con métodos de es- Una vez decidido el tipo de función base
tructura electrónica es el programa GAMESS (STO/GTO) y su ubicación (núcleos), el fac-
(General Atomic Molecular and Electronic tor más importante es la cantidad de fun-
Structure System) que es un programa esta- ciones que se utilizarán. El menor número
blecido de simulación de quı́mica cuántica ab posible de funciones es un conjunto mı́nimo
initio [5]. de bases. Solo se emplean suficientes funcio-
nes orbitales atómicas para contener todos
Una de las aproximaciones inherentes a esen- los electrones de los átomos neutros. Los cin-
cialmente todos los métodos ab initio es la juntos base se clasifican en [6]:
introducción de un conjunto de bases finitas. Conjunto de base mı́nima: contiene
Las funciones centradas en el átomo utiliza- el mı́nimo número de funciones por áto-
das para describir los orbitales atómicos se mo para describir los orbitales atómicos
conocen como funciones base y forman colec- ocupados [7]. Para hidrógeno (y helio)
tivamente un conjunto base [6]. esto significa una sola función s. Para la
primera fila de la tabla periódica, signi-
Expandir una función desconocida, como un fica dos funciones s (1s y 2s) y un con-
orbital molecular, en un conjunto de funcio- junto de funciones p (2px, 2py y 2pz ).
nes conocidas no es una aproximación, si la Para los elementos de la segunda fila, se
base está completa. Sin embargo, una ba- utilizan tres funciones s (1s, 2s y 3s) y
se completa significa que se debe utilizar un dos conjuntos de funciones p (2p y 3p)
número infinito de funciones, lo cual es impo- [6].
sible en los cálculos reales. Cuando se utiliza
una base finita, solo se pueden representar los Bases tipo STO-nG: propuesta por
componentes del MO (Molecular Orbital) a John Pople en 1969, la más usada es
lo largo de los ejes de coordenadas correspon- la STO-3G. Son bases mı́nimas en las
dientes a la base seleccionada. Cuanto más que se intenta simular la correspondien-
pequeña es la base, más pobre es la represen- te base STO empleando para ello n fun-
tación y cuanto mejor una función de base es ciones GTO por cada STO. Ası́ para
capaz de reproducir la función desconocida, cada orbital se tienen 3 funciones. En
menos funciones base son necesarias para lo- estas bases el GTO individual es llama-
grar un determinado nivel de precisión [6]. do orbital primitivo y la combinación de
funciones es llamada función contraı́da
Hay dos tipos de funciones base comúnmen- o contraccion [8].
te utilizadas en los cálculos de estructuras Bases tipo N-21G y N-31G: de es-
electrónicas: Orbitales de tipo Slater (STO) tas funciones las más populares son las

3
 3-21G, y las 6-31G en ellas se indica OF GAMESS TERMINATED NORMALLY,
que para el core (centro) se ha empleado se tomó el valor de energı́a total donde se
una GTO contraı́da con N GTO’s pri- observó TOTAL ENERGY, este cálculo se
mitivas, y la capa de valencia se repre- realizó para el catión y el anión del mis-
senta con dos funciones base contraı́das, mo átomo. Mediante los datos obtenidos se
la primera con 2 o 3 primitivas respec- cálculo la energı́a de ionización y la afinidad
tivamente, y la segunda con tan sólo 1 electrónica los cuales fueron comparados con
función GTO primitiva. Los parámetros el teorema de Koopmans (energı́a de ioni-
de estas funciones (coeficientes y expo- zación es igual a energı́a HOMO del átomo
nentes) se optimizan variacionalmente neutro y la afinidad electrónica es igual a
empleando el método UHF [8]. energı́a LUMO del átomo neutro).
Bases extendidas: también en este
Posteriormente se calculó la energia de ioni-
grupo se tiene un notable conjunto de-
zación vertical para la molécula de agua usan-
bido a Pople, y se suele denominar
do la base 6-31G(d,p) y también la afinidad
añadiendo ** si se extiende el conjunto
electrónica para la molécula de dióxido de
de funciones base con funciones de po-
nitrógeno (N O2 , usando las bases 6-31G(d,p)
larización (incrementan una unidad el
y 6-31+G(d,p) y se compararon los valores
momento angular), y ++ si se añaden
de esta dos bases.
funciones difusas (para estados exita-
dos, aniones, moléculas de longitudes de
Para el cálculo de la geometrı́a de equilibrio
enlace grandes) [8].
de las moléculas de H2 , N2 y N H3 , seguido
A continuación se describen lo procedimien- se calcularon las energı́as para cada molécu-
tos realizados durante la práctica para la de- las con las funciones base 3-21G, 6-31G, 6-
terminación de energı́a de ionización y adini- 31G(d,p) y 6-31++G(2df, 2pd) donde se tuvo
dad electrónica de átomos y moléculas. en cuenta la energı́a total y el tiempo que tar-
do el programa en realizar el cálculo.

Detalles computacionales Por último se dibujo la molécula de agua y se

hizo un cálculo para optimizar la geometrı́a
Se realizó el cálculo de la energı́a total
utilizando un nivel de teória HF/STO-3G,
de los primeros 10 elementos de la tabla
una vez optimizada la molécula se guardo y
periódica con el método Hartree-Fock y 6-
posteriormente se realizó un cálculo de punto
31G y 6 bases de cálculo distintas (STO-
sencillo con la geometrı́a optimizada, final-
3G, 3-21G, 6-31G, 6-31++G, 6-31G(d,p) y
mente se visualizaron los orbitales molecula-
6-31++G(d,p)), se tuvo en cuenta la car-
res HOMO y LUMO. Se repitió el procedi-
ga y la multiplicidad, para asegurar que el
miento para la molécula de etileno y benceno
calculo terminara satisfactoriamente se ob-
en avogadro.
servo en el archivo de salida EXECUTION

Resultados y discusión
1. Cálculo de la energia de ionización y afi- los primeros 10 átomos de la tabla periodi-
nidad electrónica de átomos. ca, posteriormente se calcuó la energı́a de
ionización y la afinidad electrónica utilizan-
Se realizó el cálculo de la energı́a total de

4
do las ecuaciones 2 y 6, respectivamente, y
ción (Tabla 1 - Tabla 10) se pueden observar
finalmente se comparó con la teorı́a de Koop-
los valores de energı́a asociados a los cálculos
mans teniendo en cuenta el valor del orbital
realizados, esto se realizó con las 6 funciones
HOMO y LUMO del átomo neutro. base registradas y teniendo en cuenta los va-
lores de energı́a total reportado en el archivo
En las tablas que se muestran a continua- de salida.

Atomo Afinidad Electronica Energı́a Ionización

H 0,0277 0,4997
He 0 0,9035
Li 0,0227 0,1981
Be 0 0,3426
B 0,01017 0,1535
C 0,05862 0,4138
N 0,00270 0,5341
O 0,0337 0,5004
F 0,1249 0,6403
Ne 0 0,7924

Tabla 1: Afinidad Electrónica y Energı́a de Ionización teórica, unidades en Hartrees

Hidrogeno Energia Ionización E.I Koopmans Afinidad Electrónica A.E Koopmans

STO-3G 0,46658185 -0,4666 0,30802409 /
3-21G 0,49619864 -0,4962 0,09577795 1,069
6-31++G 0,4988011 -0,4988 0,01172851 0,0968
6-31G(d,p) 0,49823291 -0,4982 0,07579093 0,9162
6-31G 0,49823291 -0,4982 0,07579093 0,9162
6-31++G(d,p) 0,4988011 -0,4988 0,01172851 0,0968

Tabla 2: Resultados obtenidos con el programa Gabedit-Gamess, unidades en Hartrees

Helio Energia Ionización E.I Koopmans Afinidad Electrónica A.E Koopmans

STO-3G 0,87603551 -0,876 2,80778396 /
3-21G 0,86029772 -0,9036 2,83567987 2,0817
6-31++G 0,86113604 -0,917 2,85565054 0,2041
6-31G(d,p) 0,86154265 -0,9141 2,85516043 1,3999
6-31G 0,86154265 -0,9141 2,85516043 1,3999
6.31++G(d,p) 0,88026839 -0,9036 2,85565054 2,0817

Tabla 3: Resultados obtenidos con el programa Gabedit-Gamess, unidades en Hartrees

5
Litio Energia Ionización E.I Koopmans Afinidad Electrónica A.E Koopmans
STO-3G 0,180078379 -0,1801 0,10225259 0,1301
3-21G 0,194418681 -0,1944 0,01475276 0,0264
6-31++G 0,195987622 -0,196 0,00578175 0,0118
6-31G(d,p) 0,195735863 -0,1958 0,02424214 0,0361
6-31G 0,195755791 -0,1958 0,02584869 0,0373
6.31++G(d,p) 0,196026622 -0,196 0,0057201 0,0118

Tabla 4: Resultados obtenidos con el programa Gabedit-Gamess, unidades en Hartrees

Berilio Energia Ionización E.I Koopmans Afinidad Electrónica A.E Koopmans

STO-3G 0,25403762 -0,254 0,10225259 0,2211
3-21G 0,29684835 -0,3056 0,01475276 0,0714
6-31++G 0,29417707 -0,3063 0,00578175 0,0725
6-31G(d,p) 0,29142203 -0,3072 0,02424214 0,03
6-31G 0,29135336 -0,3015 0,02584869 0,0822
6.31++G(d,p) 0,29420186 -0,3073 0,0057201 0,0299

Boro Energia Ionización E.I Koopmans Afinidad Electrónica A.E Koopmans

STO-3G 0,20051828 -0,2005 -0,3862544 0,2914
3-21G 0,2933854000 -0,3072 -0,1099228 0,0634
6-31++G 0,28966706 -0,3078 -0,0543396 0,0277
6-31G(d,p) 0,28840729 -0,3062 -0,1161981 0,0744
6-31G 0,28545018 -0,3008 -0,1222658 0,0759
6.31++G(d,p) 0,292201198 -0,3128 -0,129057 0,0291

Tabla 5: Resultados obtenidos con el programa Gabedit-Gamess, unidades en Hartrees

Carbono Energia Ionización E.I Koopmans Afinidad Electrónica A.E Koopmans

STO-3G 0,219221101 -0,2192 -0,2722684 0,2723
3-21G 0,309057792 -0,3405 -0,045313 0,0726
6-31++G 0,310838746 -0,355 0,02411662 0,0208
6-31G(d,p) 0,30731253 -0,3457 -0,027468 0,0634
6-31G 0,301128052 -0,3422 -0,0345309 0,0669
6.31++G(d,p) 0,308485896 -0,3566 0,02492779 0,0209

Tabla 6: Resultados obtenidos con el programa Gabedit-Gamess, unidades en Hartrees

Nitrogeno Energia Ionización E.I Koopmans Afinidad Electrónica A.E Koopmans

STO-3G 0,435252932 -0,353 -0,3275516 0,2042
3-21G 0,489020279 -0,4536 -0,1176798 0,0548
6-31++G 0,497628275 -0,483 -0,02338 0,1458
6-31G(d,p) 0,499250861 -0,4715 -0,0688739 0,0257
6-31G 0,495651016 -0,4730 -0,0715273 0,0245
6.31++G(d,p) 0,500642614 -0,4796 -0,0265242 0,1422
6
Tabla 7: Resultados obtenidos con el programa Gabedit-Gamess, unidades en Hartrees
 Oxigeno Energia Ionización E.I Koopmans Afinidad Electrónica A.E Koopmans
STO-3G 0,408010743 -0,2192 -0,2722684 0,2349
3-21G 0,4925755 -0,3405 -0,045313 0,0454
6-31++G 0,503420586 -0,355 0,02411662 0,1906
6-31G(d,p) 0,498309588 -0,3457 -0,027468 0,005
6-31G 0,501828666 -0,3422 -0,0345309 0,0048
6.31++G(d,p) 0,500320636 -0,3566 0,02492779 0,1823

Tabla 8: Resultados obtenidos con el programa Gabedit-Gamess, unidades en Hartrees

Fluor Energia Ionización E.I Koopmans Afinidad Electrónica A.E Koopmans

STO-3G 0,624323355 -0,4626 0,37319028 0,3732
3-21G 0,721343211 -0,5943 0,072874024 0,1345
6-31++G 0,736843503 -0,6482 -0,049854522 0.0367
6-31G(d,p) 1,365892246 -0,6236 0,645730186 0.1964
6-31G 1,350555834 -0,6254 0,62925749 0,1887
6.31++G(d,p) 1,377906566 -0,6449 0,593447322 0,0391

Tabla 9: Resultados obtenidos con el programa Gabedit-Gamess, unidades en Hartrees

Neon Energia Ionización E.I Koopmans Afinidad Electrónica A.E Koopmans

STO-3G 0,543052789 -0,5431 126,604525 0
3-21G 0,708434028 -0,7903 127,803825 2,6873
6-31++G 0,7252445 -0,8536 128,473877 0,2928
6-31G(d,p) 0,723726953 -0,8307 128,474464 1,6385
6-31G 0,722164571 -0,8308 128,473877 1,7558
6.31++G(d,p) 0,731083708 -0,8536 128,484025 0,2919

Tabla 10: Resultados obtenidos con el programa Gabedit-Gamess, unidades en Hartrees

Comparando las energı́as de ionización La energı́a I teórica del Helio es 0,90356986

(I ) del Hidrógeno (Tabla 2) con la energı́a Ha, observando las energı́as I calculadas que
I teórica, la cual es 0,49973175 Ha, se ob- se encuentran en la Tabla 3 se observa que
serva que realizando el cálculo con GAMESS son muy cercanas, sobre todo la calculada
la energı́a I obtenida es cercana a la teórica. con la funcion base 6-31++G(d,p), aunque la
En este caso la más cercana es la cálcula- energı́a proporcionada por el orbital HOMO
da con la base 6-31G++(d,p), con esto se es igual a la teórica. Con esto se evidencia de
evidencia que aumentando el número de fun- nuevo la tendencia que tiene la energı́a I de
ciones base se obtiene una energı́a I con un ser más precisa cuando se aumenta la funcion
valor casi igual al teórico. Teniendo encuen- base.
ta el teorema de Koopmans (ecuación 5) se
evidencia que la energı́a I obtenida por el or- Con la energia I calculada (Tabla 1) y de
bital HOMO es igual a la energı́a I calculada. Koopmans de los atomos de Li, Be, N y O
de igual manera que con el H y He coini-

7
ciden con la energı́a teórica, obteniendose
tambien el dato más cercano con la base 6- En general se puede evidenciar que para la
31++G(d,p) afinidad electrónica calculada y el cálculo
del orbital LUMO teniendo en cuenta el teo-
La energia I calculada de los átomos de B, rema de Koopmans estas no se acercan al
C, F y Ne no se acerca al valor teórico con las valor teórico reportado en la Tabla 1, para
diferentes funciones bases utilizadas del pro- ningún átomo, lo cual nos indicarı́a que no es
grama Gamess, en algunos de ellos el valor confiable el dato calculado para la afinidad
esta muy por debajo del teórico y en otros electrónica.
por encima. Por otro lado teniendo en cuenta
el teorema de Koopmans podemos evidenciar Con las energı́as de ionización calculadas se
que la energı́a I proporcionada por el orbi- realizó la gráfica representada en la Figura 1,
tal HOMO del átomo neutro del F Y Ne se en la cual se observa que la energı́a claulada
va acercando a la teórica a medida que se es más precisa utilizando más funciones base
aunmenta la función base, evidenciando la que con la base STO-3G, siendo esta la ba-
tendencia que tiene la energı́a I de ser más se que proporciona un valor más alejado al
precisa cuando se aumenta el número de fun- teórico.
ciones base.

Figura 1: Energı́a de Ionización

2. Cálculo de la energı́a de ionización y afinidad electrónica para móleculas

Molécula Energı́a (Ha) E.I Koopmans E.I Calculada

Agua -76,0160947857 -0,5001 -0,498303998
Agua Ionizada -76,5143987842 /

Tabla 11: Energı́a de Ionizacion H2 O, unidades en Hartrees

8
 Base Afinidad electrónica N O2 (Ha)
6-31G(d,p) 0,0995
6-31+G(d,p) 0,0896
Tabla 12: Afinidad Electrónica N O2

Realizando el cálculo de energı́a total de teórica de la molécula de N O2 , 0.0835312

la molécula de agua y la molécula de agua io- Ha, la base que porporciona un resultado más
nizada, aplicando la ecuación 2 se obtuvó la cercano al teorico es la base 6-31+G(d,p).
energı́a de ionización, en la Tabla 11 se en-
cunetran las energı́as obtenidas, donde se 3. Cálculo de geometrı́a de equilibrio con di-
evidencia que la energı́a calculada se acerca ferentes bases.
mucho al valor arrojado por el orbital HOMO
(E.I Koopmans). Comparando las energı́as Para el cálculo de la geometrı́a de equili-
obtenidas con la energı́a de ionización del brio de las moléculas de H2 , N2 , Y N H3 se
agua, la cual es de 0.463813 Ha, se obtuvó un utilizaron diferentes funciones base (3-21 G,
valor muy cercano al teórico. 6-31 G, 6-31G (d,p) y 6-31++G(2df,2pd)).
Con estas base se encontró la enrgı́a total de
En la Tabla 12 se encuentra la afinidad cada molécula teniendo en cuenta la función
electrónica del dióxido de nitrógeno (N O2 ) aplicada y el tiempo que tardaba el progra-
calculada con dos funciones base distintas. ma Gabbedit en arrojar el archivo de salida
Teniendo en cuenta la afinidad electrónica correspondiente.

Teoria Energia total (KJ/mol) Numero de funciones Tiempo de calculo(s)

3-21 G -108.3543 6 0.36
6-31 G -108.7276 10 0.28
6-31G (d,p) -109.1623 12 0.36
6-31++G(2df,2pd) -109.2319 18 0.36

Tabla 13: Resultados obtenidos con el programa Gabedit-Gamess para la molécula H2

Teoria Energia total (KJ/mol) Numero de funciones Tiempo de calculo(s)

3-21 -10445.6269 6 0.4
6-31 G -10500.3190 10 0.36
6-31G (d,p) -10512.0016 12 0.64
6-31++G(2df,2pd) 10513.5895 18 0.63

Tabla 14: Resultados obtenidos con el programa Gabedit-Gamess para la molécula N2

Teoria Energia total (KJ/mol) Numero de funciones Tiempo de calculo

3-21 -5391.1089 6 0.36
6-31 G -5419.4111 10 0.36
6-31G (d,p) -5422.3081 12 0.64
6-31++G(2df,2pd) -5423.2838 18 0.448

Tabla 15: Resultados obtenidos con el programa Gabedit-Gamess para la molécula NH3

9
 En la Tabla 13 se consignaron los datos
obtenidos de estos cálculos para la molécula
de H2 donde se muestra la tendencia que tie-
nen estos datos, al observar los datos de la
energı́a total se puede concluir que a medida
que aumentan las funciones base aplicadas en
el cálculo, mayor es la energia total del sis-
tema. En los valores del tiempo para llevar a Figura 2: Orbital molecular LUMO del ben-
cabo el calculo no se observa mucha variación ceno
entre ellos.

En la Tabla 14 se encuentran los valores de

los cálculos correspondientes para la molécu-
la N2 con las diferentes funciones base uti-
lizadas, donde se tiene un comportamiento
parecido con la molécula de H2 , a medida
que aumenta el numero de funiones que se
aplican a la molécula la energı́a del sistema Figura 3: Orbital molecular HOMO del ben-
también aumenta obteniendo un valor mas ceno
cercano al teórico. En este caso los valores
del tiempo empleados para realizar el calculo En la Figura 2 y 3 se pueden observar los
si varian, donde los resultados obtenidos que orbitales LUMO y HOMO del benceno. En
tienen mas funciones tienen un tiempo mayor el orbital HOMO del benceno la energÍa se
de cálculo. encuentra sobre los carbonos y los dobles en-
laces de este.
Por último, en la Tabla 15 se encuentran
los resultados asociados al cálculo de energı́a
de la molécula de N H3 , donde el comporta-
miento que se evidenció anteriormente en el
H2 y N2 es el mismo, a medida que la canti-
dad de funciones aplicadas sobre la molécula
de aumenta, el valor de energı́a reportado se
acerca mas al encontrado en la teorı́a. La va-
Figura 4: Orbital molcular HOMO del etileno
riación del tiempo aplicado en el cálculo no
varı́a mucho entre sı́, aunque si se observa que
los de más funciones base tieneden a tener un
tiempo más largo de cálculo.

4. Cálculo de la energia con el método

Hartree-Fock.

Para realizar el cálculo de la energı́a de la

molécula de agua y etileno se utilizó la base Figura 5: Orbital molcular LUMO del etileno
STO-3, donde se evidenciaron los orbitales
HOMO y LUMO en la molecula de cada Para el caso de la molécula de etileno se
compuesto. puede observar en la Figura 4 y 5 los orbitales

10
HOMO y LUMO de la molécula. En el orbi- métodos existentes para la aplicación de
tal LUMO se puede observar que la energÍa se los mismos.
encuentra sobre los carbonos en la dirección
opuesta a los enlaces que posee el átomo y en Se puede concluir que al aunmentar el
el orbital HOMO de la molécula la energı́a número de funciones base utilizadas, el
presente se encuentra posada sobre el doble cálculo obtenido en el programa Ga-
enlace entre los carbonos de la molécula. mess se acerca más al teórico.

Llevar acabó primero el cálculo con es-

tructura en equilibrio nos permite ob-
tener una mejor geometria optimizada
para posteriormente realizar un cálculo
de punto sencillo lo que permite obser-
var los orbitales LUMO Y HOMO de la
molecula.
Figura 6: Orbital LUMO del agua
El tiempo de cálculo es pequeño, por
esto es util ya que no se invierte tanto
tiempo esperando el cáculo,corriendo el
riesgo de que salga abnormally.

El teórema de koopmans se acerca más

a los valores teóricos que la energı́a de
ionización calculada.

Figura 7: Orbital HOMO del agua

Referencias
La Figura 6 y 7 muestra los orbitales LU-
MO y HOMO del agua. El orbital LUMO se [1] R. Shankar, K. Senthilkumar, P. Kolan-
encuentra ubicado sobre los hidrógenos de la daivel, Int. J. Quantum Chem., 2008,
molécula de agua, por el contrario en el orbi- 109, 764.
tal HOMO se puede observar que está ubica-
[2] El comportamiento periódi-
do sobre el oxı́geno.
co de los elementos, en lı́nea,
https://andoni.garritz.com/
Conclusiones documentos/Capitulo_9.pdf, (consul-
tado en marzo de 2020)
El programa Gamess permite analizar
[3] P. Lang, B. C. Smith, J. Chem. Educ.,
caracteristicas de una molécula o áto-
2003, 80, 938.
mo realizando algunas variaciones en el
sistema. En este caso se realizaron op- [4] P. Politzer, J. S. Murray, F. A. Bulat, J.
timizaciones de geometrı́a para obser- Mol. Model., 2010, 16, 1731-1732
var el orbital LUMO y HOMO de cier-
tas móleculas. Tambien permite cálcu- [5] K. K. Baldridge, J. P. Greenberg,
lar la energı́a de ionización y afinidad S. T. Elbert, S. Mock, P. Papa-
electrónica de los primeros 10 elemen- dopoulos, IEEE Xplore., 2002, DOI:
tos de la tabla periódica, utilizando los 10.1109/SC.2002.10014.

11
[6] Basis functions in quantum chemistry, [9] Chemical elements listed by
en lı́nea, https://documents.epfl.ch/ ionization energy, en lı́nea,
users/m/mk/mkilic/public/chapter_ https://www.lenntech.com/
3.pdf, (consultado en marzo de 2020). periodic-chart-elements/
ionization-energy.htm, (consul-
[7] Quı́mica Cuántica I: Introducción
tado en marzo de 2020).
a la quı́mica computacional, en
lı́nea, http://depa.fquim.unam.
[10] Electron Affinity of the elements, en
mx/jesusht/qci_quim_comp.pdf,
lı́nea, https://periodictable.com/
(consultado en marzo de 2020).
Properties/A/ElectronAffinity.v.
[8] E. F. Maroto, Universidad de Alicante, log.html, (consultado en marzo de
2017, 59-60. 2020).

12