Práctica 2.

Cálculos quı́mico-cuánticos para la

descripción de propiedades termodinámicas
Departamento de Fı́sica y Quı́mica Teórica, UNAM.

Quı́mica Computacional

Avilés Sánchez Alejandro

16 de octubre de 2014

Resumen
En el presente trabajo se estudió la reacción de disociación en fase gaseosa de la molécula
homuclear de F2 : F2 (g)
2F (g). Sobre la base de cálculos quı́mico-cuánticos de estruc-
tura electrónica atómica y molecular ab initio, empleando los niveles de cálculo SCFHF
y CISD, se obtuvieron algunas propiedades termodinámicas del sistema.

1. Introducción
La termodinámica estudia las propiedades macroscópicas de la materia, mientras que la mecáni-
ca cuántica se utiliza para estudiar las propiedades moleculares. El nexo de unión entre la
mecánica cuántica y la termodinámica lo proporciona la mecánica estadı́stica, cuyo propósi-
to es deducir las propiedades macroscópicas de la materia a partir de las propiedades de las
moléculas que componen el sistema. A causa del inmenso número de moléculas existentes en un
sistema macroscópico, se usan métodos estadı́sticos en vez de estudiar el movimiento de cada
molécula del sistema.[3]

Las fórmulas fundamentales de la mecánica estadı́stica relacionan las propiedades termo-

dinámicas con los niveles de energı́a mecano-cuánticos del sistema termodinámico. La manera
de relacionar propiedades termodinámicas con propiedades moleculares de hace del siguiente
modo:

i) Escribimos el operador hamiltoniano H b para el sistema termodinámico. (Esto requiere el

conocimiento de la naturaleza de las energı́as intramoleculares y el conocimiento de las fuerzas
intermoleculares.)

ii) Resolvemos la ecuación de Schrödinger Hψ

b j = Ej ψj para el sistema termodinámico com-
pleto para obtener las energı́as mecano-cuánticas Ej de los posibles estados cuánticos j del
sistema.

iii) Evaluamos la función de partición canónica Z = j e−Ej /kT , donde la suma se extiende
P
sobre todos los estados cuánticos del sistema termodinámico.

iv) Utilizamos lnZ para calcular las propiedades termodinámicas de un sistema a T , V y com-
posición fijos, a partir de las siguientes ecuaciones

1
  
∂lnZ
P = kT (1)
∂V T,NB
 
2 ∂lnZ
U = kT (2)
∂T V,NB

S = U/T + klnZ (3)

A = −kT lnZ (4)

 
∂lnZ
µB = −RT (5)
∂NB T,V,Nc 6=B

donde NB es el número de moléculas de la especie B, Ej es la energı́a del estado cuántico j y

la suma involucra a todos los estados cuánticos posibles del sistema termodinámico.

Siempre que podamos resolver la ecuación de Schrödinger para el sistema termodinámico,

la mecánica estadı́stica nos permite calcular todas las propiedades termodinámicas del sistema.

2. Marco Teórico
2.1. Función de partición canónica de un gas ideal puro
En un gas ideal (donde no hay fuerzas intermoleculares), la energı́a del sistema es la suma de
las energı́as moleculares, y la función de partición canónica es Z = z N /N !, donde z es la función
de partición molecular. Puesto que lo que aparece en las ecuaciones (1) a (5) es lnZ (en vez de
Z) para las propiedades termodinámicas, tomaremos lnZ
lnZ = ln z N /N ! = N lnZ − lnN !


(6)


Es bastante fatigoso calcular un número del orden de ln 1023! . Afortunadamente, existe
una excelente aproximación. La fórmula de Stirling, válida para N elevado, es
 
1/2 N +1/2 −N 1 1
N ! = (2π) N e 1+ + − ... (7)
12N 288N 2
para N grande. Tomando el logaritmo natural de (7), tenemos
   
1 1 1
lnN ! = ln2π + N + lnN − N + ln 1 + + ... (8)
2 2 12N
Puesto que N es del orden de 1023 , tenemos que N + 21 ≈ N , 12 ln2π−N ≈ −N y ln (1 + 1/12N + . . .) ≈
ln1 = 0. Por tanto,
lnN ! ≈ N lnN − N (9)
para N grande.

Consideremos ahora la función de partición molecular del gas ideal

X
z= e−εr /kT (10)
r

2
donde εr es la energı́a del estado cuántico molecular r. Generalmente, es buena aproximación
escribir la energı́a molecular como la suma de las energı́as traslacional, rotacional, vibracional
y electrónica
εr = εtr,s + εrot,t + εvib,v + εelec,u (11)
donde los subı́ndices s, t, v, u indican los estados traslacional, rotacional, vibracional y electróni-
co. El estado molecular r está definido por los estados traslacional, rotacional, vibracional y
electrónico s, t, v y u, por tanto, la suma entorno a r implica una suma entorno a s, t, v y u.
Los números cuánticos de los cuatro tipos de energı́a varı́an de un modo independiente unos de
otros. De esta manera

e−β (εtr,s +εrot,t +εvib,v +εelec,u )

X XXXX
z= e−βεr = (12)
r s t v u

X X X X
= e−βεtr,s e−βεrot,t e−βεvib,v e−βεelec,u (13)
s t v u

z = ztr zrot zvib zelec (14)

lnz = lnztr + lnzrot + lnzvib + lnzelec (15)
Cuando la energı́a molecular es la suma de distintos tipos de energı́a, la función de partición
molecular puede factorizarse en un producto de funciones de partición, cada una de ellas co-
rrespondiente a un tipo de energı́a.[5]

2.2. Termodinámica estadı́stica para gases ideales monoatómicos y

diatómicos
Presentaremos ahora las expresiones de ztr , zrot , zvib , y zelec para un gas ideal puro de moléculas
diatómicas, para expresar las propiedades termodinámicas del gas en términos de sus propie-
dades moleculares.

2.2.1. Función de partición traslacional

La función de partición traslacional de un gas ideal puro formado por moléculas diatómicas es
 3/2
2πmkT
ztr = V (16)
h2

V es el volumen que contiene el gas ideal, m es la masa de la molécula del gas estudiado.

La ecuación (16) muestra que ztr es una función de estado. Debido a que ztr es proporcional
a V , ztr es extensiva.

2.2.2. Función de partición rotacional

La evaluación de zrot extendiendo la suma a todos los estados rotacionales, conduce a la relación
"  2  3 #
T 1 Θrot 1 Θrot 4 Θrot
ztr = 1+ + + + ... (17)
σΘrot 3 T 15 T 315 T

la ecuación (17) es una fórmula utilizada para moléculas diatómicas homo y heteronucleares.
Definimos el número de simetrı́a σ como 2 para una molécula diatómica homonuclear y 1 para

3
una molécula diatómica heteronuclear. Asimismo, se ha definido la temperatura rotacional
caracterı́stica Θrot como
}2 Bhc
e
Θrot = = (18)
2Ik k
I es el momento de inercia de la molécula y B
e la constante rotacional

h
B
e0 = B0 /c = (19)
8π 2 I0 c
donde I0 es el momento de inercia promediado sobre la vibración del punto cero.

2.2.3. Función de partición vibracional

La aproximación
del oscilador armónico da, para las energı́as vibracionales de una molécula
diatómica, v + 12 hν, donde v va de 0 a ∞ y no hay degeneración. La elección del nivel cero
de energı́a es arbitraria, y lo más usual en mecánica estadı́stica molecular es tomarlo en el nivel
de energı́a más bajo permitido de la molécula. Este nivel tiene una energı́a rotacional igual a
cero y una energı́a vibracional igual a 21 hν. Escribiremos por tanto εvib = v + 12 hν − 12 hν =

vhν, donde εvib se mide a partir del nivel v = 0. Para que el razonamiento fuera plenamente
consistente, deberı́amos haber restado de εtr la energı́a traslacional del punto cero, pero ésta es
lo suficientemente pequeña como para ser absolutamente despreciable.

Tenemos que
1 1
ztr = −hν/kT
= −Θ
(20)
1−e 1 − e vib /T
donde Θvib es la temperatura vibracional caracterı́stica. La expresión de εvib obtenida
suponiendo el modelo del oscilador armónico es tan sólo una aproximación, puesto que no
se tiene en cuenta los efectos de anarmonicidad, distorsión centrı́fuga e interacción rotación-
vibración. Afortunadamente, las correcciones que habrı́a que introducir para tener en cuenta la
anarmonicidad, etc., son en general pequeñas (excepto temperaturas extremadamente altas).

2.2.4. Función de partición electrónica

La expresión de zelec como una suma extendida a los niveles de energı́a electrónicos (en lugar
de extendida a los estados) es

zelec = gel,0 + gel,1 e−βεel,1 + gel,2 e−βεel,2 + . . . (21)

donde se incluye la degeneración gel de cada nivel electrónico. Los niveles de energı́a electróni-
cos se numeran 0, 1, 2, ... en orden de energı́a creciente y gel,0 , gel,1 , gel,2 son sus respectivas
degeneraciones. Se toma la energı́a del nivel electrónico fundamental como e−βεel,0 = 0. La
energı́a de cualquier nivel electrónico excitado se mide entonces con relación al nivel electrónico
fundamental.

La contribución electrónica a la función de partición molecular, toma al cero de energı́a

electrónica como la energı́a de los átomos separados en el estado basal:

zelec = gel,1 eβDe + gel,2 e−βεel,2 + . . . (22)

Puesto que no existe una fórmula general para las distintas εel , hay que sumar por separado
cada término de la serie empleando las energı́as electrónicas observadas espectroscópicamente.

4
 Para casi todas las moléculas diatómicas, εel,1 es mucho mayor que kT a temperatura ambien-
te, y todos los términos de (21), salvo el primero, contribuyen con una cantidad insignificante
hasta temperaturas de 5000 ó 10, 000 K. Por tanto,
zelec = gel,0 (23)
para temperaturas no muy altas. En la mayorı́a de las moléculas diatómicas, el nivel electrónico
fundamental no está degenerado: gel,0 = 1. Excepciones importantes son la molécula de O2
(gel,0 = 3) y el N O (gel,0 = 2).

2.3. Constantes de equilibrio

P
Para la reacción entre gases ideales 0 −→ i νi Bi (las magnitudes νi son los coeficientes
estequiométricos de las distintas especies Bi implicadas en la transformación quı́mica):
 ◦  νi
Y νi ∆U0 Y zi
Kc = (Ci ) = exp − (24)
i
RT i
V N A

◦
donde Kc es la constante de equilibrio en función de las concentraciones. La magnitud ∆U0
corresponde a la diferencia U − U0 .

Adicionalmente podemos definir una constante de equilibrio en función de las presiones

parciales KP que puede obtenerse una vez conocida la Kc :
KP = (kT ) Kc ∆ν (25)
con ∆ν = moles de productos − moles de reactivos.

2.3.1. Dependencia de la constante de equilibrio Kp con la temperatura

La constante de equilibrio KP para una reacción entre gases ideales es sólo función de la
temperatura. La dependencia de KP con la temperatura es explı́cita en la ecuación de Van’t
Hoff : ◦
dlnKP ∆H
= (26)
dT RT 2
◦ ◦
En (25) ∆H = ∆HT es la variación de entalpı́a estándar de la reacción entre gases ideales
◦

a la temperatura T . Cuanto mayor sea el valor de ∆H , más rápida será la variación de la

constante de equilibrio KP con la temperatura.

De la ecuación (25) se desprende la relación

◦  
KP (T2 ) ∆H 1 1
ln ≈ − (27)
KP (T1 ) R T1 T2
◦
útil para estimar el valor de ∆H a partir del conocimiento de la constante KP a dos distintas
temperaturas.

3. Sistema de estudio
En el presente trabajo se considera la disociación en fase gaseosa de la molécula de flúor
F2 (g)
2F (g) (28)
◦
Para este sistema, surge el interés de estimar el valor de ∆r H a partir de cálculos quı́mico-
cuánticos.

5
4. Objetivos
4.1. Objetivo General
Familiarizarse con los métodos de estructura electrónica y obtener el valor propiedades termo-
dinámicas del sistema de estudio, a partir de cálculos quı́mico-cuánticos.

4.1.1. Objetivos particulares

Sobre la base de cálculos teóricos:

Obtener las funciones de partición moleculares de gases ideales zF2 y zF .

Determinar la constante de equilibrio KP de la reacción propuesta a 900 y 1200K.

◦
Estimar el valor de ∆r H del proceso de disociación por medio de la ecuación de Van´t
Hoff.

5. Metodologı́a
Se emplearon métodos ab initio a nivel de teorı́a SCFHF y CISD con el conjunto de funciones
de base 6-311++G(d,p)[6]. Las etapas de cálculo fueron las siguientes:

Obtención de la geometrı́a de equilibrio (en este caso la distancia de enlace de equilibrio

Re ) con los distintos niveles de teorı́a.

Cálculo de la constante rotacional B

e y la energı́a de disociación De de la molécula de
flúor.

Estimación de la primera energı́a de excitación εel,1 del átomo de F usando CIS para la
comparación a nivel SCFHF y usando CIS(D) para la comparación a nivel CISD.

Para todos los cálculos efectuados, se utilizó el programa GAMESS-US (sistema de estructura
electrónica atómica y molecular general).

6. Resultados
6.1. Resultados de cálculos
La tabla 1 y 2 resumen los resultados obtenidos correspondientes a la etapa de cálculos quı́mico-
cuánticos propuestos en la metodologı́a.

Método Re (Å) Energı́a de punto cero (J) De (J) e (cm−1 )

B
SCFHF 1.330 7.323x10−20 1.538x10−18 0.567
CISD 1.432 1.753x10−20 1.007x10−19 0.865
Experimental* 1.410 - 2.724x10−19 -

Tabla 1. Diferentes parámetros calculados para la molécula de F2 con los niveles de teorı́a
SCFHF y CISD.*[4]

6
 Método εel,1 (J)
SCFHF 8.010x10−21
CIS 8.026x10−21
CISD 8.026x10−21
Experimental* 8.026x10−21
Tabla 2. Primera energı́a de excitación εel,1 del átomo de F usando diversos métodos ab initio
.*[4]

6.2. Funciones de partición de gases ideales zF2 y zF

6.2.1. Función de partición zF
Las tablas 3 y 4 muestran los valores calculados de las funciones de partición del átomo de flúor
a las temperaturas de 900 y 1200K.

T (K) ztr zelec zF

900 1.327x1037 5.049 6.700x1037
12000 2.044x1037 5.233 1.069x1038
Tabla 3. Valores correspondientes a las funciones de partición traslacional, electrónica y atómi-
ca para el F calculados con el método SCFHF.

T (K) ztr zelec zF

900 1.327x1037 5.048 6.698x1037
12000 2.044x1037 5.231 1.069x1038
Tabla 4. Valores correspondientes a las funciones de partición traslacional, electrónica y atómica
para el F calculados con el método CISD.

6.2.2. Función de partición zF2

Las tablas 5 y 6 resumen los cálculos necesarios para la obtención de las funciones de partición
de la molécula diatómica de flúor zF2 a dos temperaturas.

T (K) ztr zrot zvib zelec zF2

900 3.755x1037 551.367 1.000 6.019x1053 1.259x1094
1200 5.782x1037 735.100 1.000 2.160x1040 9.182x1080
Tabla 5. Funciones de partición traslacional, rotacional, vibracional, electrónica y molecular
para la molécula homonuclear de diflúor, estimadas mediante cálculos de tipo SCFHF.

T (K) ztr zrot zvib zelec zF2

900 3.755x1037 361.380 1.063 3320.466 4.791x1043
1200 5.782x1037 481.785 1.136 437.409 8.997x1042
Tabla 6. Funciones de partición traslacional, rotacional, vibracional, electrónica y molecular
para F2 , estimadas mediante cálculos CISD.

7
6.3. Constantes de equilibrio KP
Las tablas 7 y 8 muestran las constantes de equilibrio en términos de las presiones parciales KP
estimadas mediante las ecuaciones (23) y (24) (véase sección 2.3) sobre la base de los resultados
obtenidos para las funciones de partición de gases ideales zF2 y zF .

T (K) KP
900 4.429x10−39
1200 2.064x10−25

Tabla 7. Constantes KP de equilibrio a dos diferentes temperaturas para el proceso de disociación

en fase gaseosa de la molécula diatómica de flúor. Se han empleado los cálculos quı́mico-cuánti-
cos de tipo SCFHF.

T (K) KP
900 1.163x1012
1200 2.105x1013

Tabla 8. Constantes de equilibrio KP para la reacción de disociación de F2 en fase gaseosa a

dos temperaturas. Se utilizaron los resultados provenientes de cálculos quı́mico-cuánticos CISD.

◦
6.4. Variación de entalpı́a estándar de reacción ∆r H
Finalmente, es la tabla 9 la que recopilan los valores de variación de entalpı́a estándar de reac-
ción hallados mediante la relación (26) (véase la sección 2.3.1.) para los dos niveles de cálculo
considerados en el presente trabajo.

◦
Método ∆r H (kJmol−1 )
SCFHF 941.985
CISD 86.673
Experimental* 157.740

Tabla 9. Valores de variación de entalpı́a estándar de reacción del proceso F2 (g)

2F (g) para
los métodos ab initio SCFHF y CISD.*[2]

7. Análisis y Discusión
7.1. Etapa de cálculos
La tabla 1 muestra que, con el método SCFHF no se obtiene una geometrı́a de equilibrio con-
sistente con el valor informado a nivel experimental. La estimación de la distancia de equilibrio
mejora al manejarse cálculos CISD (aunque el valor de Re estimado con CISD continúa siendo

8
aún ligeramente más grande que el experimental).

La energı́a de disociación de enlace en equilibrio De calculada haciendo uso del método ab

initio SCFHF, resulta ser mucho más grande que el valor hallado experimentalmente (de hecho,
es alrededor de 5.6 veces mayor que la De experimental). En contraste, el valor de De con el
nivel CISD es menor en relación con el informado en tablas experimentales (la De experimental
es prácticamente el triple de la calculada empleando el método CISD).

Las constantes rotacionales y energı́as de punto cero, varı́an de modo poco sustancial entre
ambos métodos.

La tabla 2 muestra los valores de la primera energı́a de excitación εel,1 del átomo de flúor. La
εel,1 obtenida con base en cálculos de tipo SCFHF es sólo ligeramente menor al valor informado
en bases de datos espectroscópicos [8]. Obsérvese también, que los valores εel,1 encontrados con
los niveles CIS y CISD son prácticamente los mismos que el experimental.

7.2. Funciones de partición de gases ideales zF2 y zF

7.2.1. Función de partición zF
Las funciones de partición traslacional vinculadas con el cálculo de zF para ambos métodos ab
initio, son en magnitud iguales, puesto que se ha empleado la ecuación (16) que no contiene
parámetros asociados con cálculos de estructura electrónica.

Las funciones de partición electrónicas no difieren significativamente entre un método y

otro, dado que las energı́as de excitación εel,1 con SCFHF y CISD son muy similares como ya
se ha discutido antes.

7.2.2. Función de partición zF2

Las funciones de partición traslacionales asociadas con el cálculo de zF2 adquieren el mismo
valor para ambos métodos. Esto se explica al considerar que la ecuación (16) (misma que fue
empleada para la obtención de la función de partición traslacional de la molécula), no incluye
variables asociadas con cálculos de estructura electrónica.

Las funciones de partición rotacionales obtenidas con SCFHF son sustancialmente distintas
de las calculadas con CISD en virtud de la diferencia que existe en las geometrı́as de equilibrio
Re y constantes rotacionales B e previstas para cada unos de los niveles de teorı́a.

Para ambos métodos, las funciones de partición vibracionales difieren muy poco en su valor,
porque las energı́a de vibración de punto cero no cambian mucho entre un método y el otro.

Algo que inmediatamente salta a la vista al observar las tablas 5 y 6, es el hecho de que las
funciones de partición electrónicas previstas mediante cálculos SCFHF, son enormes compara-
das con las halladas a través de cálculos CISD. Este hecho tiene su origen en los cálculos de
energı́a de disociación De . La De obtenida con SCFHF es alrededor de 15 veces mayor que la
energı́a de disociación encontrada con CISD. Ésta inmensa diferencia, aunada con el hecho de
que la contribución electrónica a la función de partición molecular, tiene a la De como exponen-
te positivo de una función exponencial [ver ecuación (22)], hace que las funciones de partición
electrónicas de ambos métodos sean totalmente diferentes.

9
 Como consecuencia principalmente de la contribución electrónica, los valores de las funcio-
nes de partición moleculares zF2 calculados con el método SCFHF, son mucho más grandes que
los valores zF2 provenientes de cálculos CISD.

7.3. Constantes de equilibrio KP

Es ampliamente conocido en quı́mica [7] que la facilidad con que la molécula de flúor se disocia es
la causa de su alta reactividad. El argumento anterior nos conduce a pensar en que naturalmente,
la constante de equilibrio quı́mico KP para la reacción de disociación F2 (g)
2F (g) estudiada
a las dos temperaturas consideradas (temperaturas relativamente altas), deberı́a en principio,
adquirir valores superiores a 1 (equilibrio desplazado hacia la formación de productos). Sin
embargo las KP calculadas mediante las ecuaciones (24) y (25) para el método de SCFHF
carecen totalmente de sentido quı́mico, pues su valor es inferior a la unidad. El motivo de que
las constantes de equilibrio en términos de las presiones parciales obtenidas a través de cálculos
tipo SCFHF no sean satisfactorias, es el hecho de que las funciones de partición moleculares zF2
estimadas, tienen un valor demasiado grande, producto (principalmente) de que la contribución
electrónica es enorme (recordemos que esto último sucede en virtud del alto valor de De obtenido
con SCFHF).

Las dos constantes de reacción KP calculadas haciendo uso del método CISD, a distintos
valores de temperatura, adquieren valores razonables por dos motivos; en primera, el valor de
ambas constantes es superior a 1, y en segundo lugar, al aumentar la temperatura, la constante
de reacción también aumenta [esto último, resulta consistente con la ecuación de Van´t Hoff
(26)]

◦
7.4. Entalpı́as estándar de reacción ∆r H
Con base en la tabla de resultados 9, se desprende que el método SCFHF [empleando el conjunto
de funciones de base 6-311++G(d,p)] presenta un error del 497 porciento en el cálculo de la
◦
variación de entalpı́a estándar de reacción ∆r H para el proceso F2 (g)
2F (g). El error tan
grande para esta estimación con SCFHF se vincula directamente con los deficientes valores de
◦
KP que se introdujeron en la ecuación (27) para calcular ∆r H .
◦
Por otra parte, el error asociado con el cálculo de la ∆r H mediante el método CISD es del
45 porciento , y es sustancialmente menor que el error para SCFHF debido a que las constantes
de reacción obtenidas con CISD sı́ son aceptables en contraste con las KP calculadas utilizando
SCFHF, las cuales carecen totalmente de sentido quı́mico.

Bajo esta perspectiva, podemos indicar que CISD constituye un mejor método para el cálculo
◦
de la variación de entalpı́a de reacción ∆r H y de las constantes de reacción KP (a diferentes
valores de T ) que el método SCFHF para la reacción F2 (g)
2F (g) [empleando el conjunto de
funciones de base 6-311++G(d,p)]. Consideramos que un aumento en el tamaño del conjunto
de funciones de base, no mejorarı́a sustancialmente las determinaciones de éstas propiedades
termodinámicas del sistema en estudio, dado que 6-311++G(d,p), ya es una base orbital de
gran tamaño.

8. Conclusiones
Se estudió la reacción de disociación en fase gaseosa F2 (g)
2F (g). Sobre la base de
cálculos quı́mico-cuánticos empleando los niveles de cálculo SCFHF y CISD, se obtuvieron

10
 las funciones de partición moleculares de gases ideales zF2 y zF , lo cual permitió determinar
◦
la constante de equilibrio KP de la reacción a 900 y 1200K y estimar el valor de ∆r H
del proceso.

CISD constituye un mejor método para el cálculo de la variación de entalpı́a de reacción

◦
∆r H y de las constantes de reacción KP (a diferentes valores de T ) que el método SCFHF
para la reacción estudiada.

Se considera que un aumento en el tamaño del conjunto de funciones de base, no mejorarı́a

sustancialmente las determinaciones de éstas propiedades termodinámicas del sistema en
estudio, dado que 6-311++G(d,p), ya es una base orbital de gran tamaño.

Referencias
[1] P. W. Atkins, y R. S. Friedman, Molecular Quantum Mechanics, 3a ed., Oxford University
Press.

[2] B. Darwent, Bond Dissociation Energies in Simple Molecules, National Standard Reference
Data System, 1970

[3] I. N. Levine, Physical Chemistry, 4 a ed. Mc-Graw-Hill, 2000.

[4] I. N. Levine, Molecular Spectroscopy, Wiley, 1975.

[5] D. A. McQuarrie, Statistical Thermodynamics, Harper Row, 1973

[6] M. M. Francl, W. J. Pietro, W. J. Hehre, J. S. Binkley, M. S. Gordon, Self-consistent

molecular orbital methods. XXIII. A polarization-type basis set for second-row elements,
The Journal of Chemical Physics 77 (7), 3654-3665

[7] R. C. Reid, J. M. Pausnitz y B. E. Poling, The Properties of Gases and Liquids, 4a. ed.,
McGraw-Hill, 1987.

[8] NIST ATOMIC SPECTRA DATABASE, http://www.nist.gov/pml/data/asd.cfm

11