Quimica Fisica de La Atmosfera
Quimica Fisica de La Atmosfera
Quimica Fisica de La Atmosfera
UNIVERSITAT DE VALENCIA
Servei de Publicacions
2005
Aquesta Tesi Doctoral va ser presentada a Valencia el dia 19 de
Novembre de 2004 davant un tribunal format per:
Depòsit legal:
I.S.B.N.:84-370-6138-5
Edita: Universitat de València
Servei de Publicacions
C/ Artes Gráficas, 13 bajo
46010 València
Spain
Telèfon: 963864115
DEPARTAMENTO DE QUÍMICA FÍSICA
CERTIFICA:
Que la memoria con título «Estudio teórico de reacciones del ozono con
diferentes radicales de interés en química atmosférica», ha sido realizada,
bajo mi dirección, por Julio José Peiró García, Licenciado en Ciencias
Químicas por la Universitat de València, para optar al grado de Doctor en
Química.
Garcilaso de la Vega
ÍNDICE DE CONTENIDOS
I. INTRODUCCIÓN
9
5.7. El método Quadratic Configuration Interaction (QCI) .............. 82
5.8. Técnicas computacionales G2 y G3........................................... 84
5.9. Métodos multiconfiguracionales................................................ 87
5.9.1. El método CASSCF ........................................................... 88
5.9.2. El método CASPT2............................................................ 91
6. Materiales y programas empleados.................................................... 93
6.1. Estaciones de trabajo.................................................................. 93
6.2. Programas utilizados .................................................................. 95
7. Objetivos ............................................................................................ 97
II. RESULTADOS
10
10.3. Entalpía de reacción y energía de activación ....................... 145
10.3.1. Energía de activación ....................................................... 145
10.3.2. Entalpía de reacción ......................................................... 151
10.3.3. El diagnóstico T1 .............................................................. 151
10.4. Conclusiones ........................................................................ 152
11. La reacción F + O3 → FO + O2 .................................................... 155
11.1. Reactivos y productos .......................................................... 157
11.2. El mecanismo de reacción.................................................... 158
11.3. Energía de activación ........................................................... 165
11.4. El diagnóstico T1 .................................................................. 168
11.5. Entalpía de reacción ............................................................. 169
11.6. Estudio comparativo de la reactividad con ozono a lo largo de
la serie isoelectrónica de radicales NH2, OH y F..................................... 170
11.7. Conclusiones ........................................................................ 173
12. La reacción NO2 + O3 → NO3 + O2 ............................................. 177
12.1. Reactivos y productos .......................................................... 180
12.2. El mecanismo de reacción.................................................... 183
12.3. Energía de activación y entalpía de reacción ....................... 187
12.4. El diagnóstico T1 .................................................................. 192
12.5. La entalpía de reacción ........................................................ 193
12.6. Cálculo de la constante de velocidad ................................... 194
12.7. Conclusiones ........................................................................ 195
12.8. Perspectivas de futuro .......................................................... 196
13. La reacción NO + O3 → NO2 + O2 .............................................. 199
13.1. Consideraciones generales ................................................... 203
11
13.2. Reactivos y productos .......................................................... 206
13.3. El mecanismo de reacción.................................................... 208
13.4. Energía de activación y entalpía de reacción ....................... 212
13.5. El diagnóstico T1 .................................................................. 217
13.6. Otros canales de reacción..................................................... 218
13.7. Conclusiones ........................................................................ 220
III. CONCLUSIONES
IV. BIBLIOGRAFÍA
V. APÉNDICES
12
AGRADECIMIENTOS
…por estos y por los otros, por los de aquí y de allá, […],
sin distinción, para todos…
13
A José Luis Fayos, compañero de colegio primero, de instituto después y
amigo desde siempre, por haber compartido conmigo la pasión por la
química, ¿o sería por el Na? Gracias también por tu compañía en tantos
crepúsculos en la inmensidad de los bosques de España y las siempre
apasionantes travesías en tu velero. ¿Hasta las Baleares será la próxima?
¿Italia? ¿Qué vendrá…?
14
A Maria Àngels Lafont, Josep Maria Lluch y Xavier Prat, por su calurosa
acogida en la Universitat Autònoma de Barcelona y por la ayuda
incondicional que me prestaron las semanas que compartí con ellos.
Al resto de mis amigos, Mabel Garrido, Ana Solbes, José Ramón Barreiro,
Yolanda Falcó, Cristina Toledo, Nieves Toledo, José Urizar, Eva Gómez y
Guillem Brocal, por las interminables partidas de frontón durante las
noches de los miércoles y las no menos interesantes cenas de los fines de
semana.
Por último, a Juanjo Serrano, Daniel Roca y Gloria Olaso, los becarios
nuevos, por ser la savia nueva y por acogerme con los brazos abiertos aun
cuando hayamos coincido en tan contadas ocasiones.
15
Prólogo
PRÓLOGO
17
Prólogo
finales del siglo XVIII, con Lavoisier, fiel defensor del método científico,
cuando la alquimia inició un rápido declive, para transformarse en la ciencia
que hoy conocemos, razón por la cual se considera a Lavoisier el padre
fundador de la química.
Como ejemplo del rigor con que Lavoisier realizaba sus experimentos,
tenemos uno que posteriormente le llevó a enunciar, junto a otros tantos, su
famosa ley de conservación de la masa.
18
Prólogo
19
Prólogo
Ojalá esta Tesis Doctoral tenga sólo una pequeña porción del rigor, la
tenacidad y la genialidad que a la historia de la química ha aportado el genial
y célebre químico francés Gay-Lussac.
20
I. INTRODUCCIÓN
Los albores de la química atmosférica
La preocupación por el medio ambiente y la calidad del aire que nos rodea
no es sólo un fenómeno de nuestros días. De hecho, ya hacia finales del siglo
XII, en plena Edad Media, el filósofo, científico y médico judío de origen
hispano, nacido en Córdoba, Moses ben Maimon Maimónides (1135-1204),
escribió:
«Comparar el aire de las ciudades con el aire de las tierras áridas y los
desiertos, es como comparar aguas sucias y túrbidas con otras límpidas y
cristalinas. En la ciudad, a causa de la altura de sus edificios, la estrechez de
sus calles y todo cuanto es desechado por sus habitantes, el aire se enrarece
y estanca haciéndose turbio, denso, brumoso y nebuloso…».
«Si el aire es alterado, aunque sólo sea ligeramente, el estado del Espíritu
Psíquico será afectado de forma perceptible. En consecuencia, encontrarás
muchos hombres en los que podrán percibirse defectos en las acciones de su
psíquico a causa de la corrupción del aire. Es decir, encontrarás que dichos
hombres adolecen de falta de entendimiento, fallos de inteligencia y ausencia
de memoria…».
23
Los albores de la química atmosférica
Figura 1.1: Portada del tratado pionero publicado por John Evelyn en el siglo XVII a
cerca de la contaminación del aire en Londres.
«Is this horrid Smoake which obscures our Church and makes our Palaces
look old, which fouls our Cloth and corrupts the Waters, so as the very Rain,
24
Los albores de la química atmosférica
and refreshing Dews which fall in the several Seasons, precipitate to impure
vapour, which, with its black and tenacious quality, spots and contaminates
whatever is exposed to it».
Fue también John Evelyn quien encontró una relación que hoy en día
nadie pone en duda entre la calidad del aire y sus efectos en la salud humana,
sentenciando lo siguiente:
25
Química atmosférica
2. QUÍMICA ATMOSFÉRICA
Hace tan sólo unas décadas, la química atmosférica era objeto de estudio
minoritario por parte de muy pocos grupos de investigación, que además se
encontraban aislados entre sí. La enseñanza de esta rama de la química en
universidades, así como la publicación de libros centrados exclusivamente en
esta materia, era así mismo muy escasa y superficial. El despertar de una
conciencia sobre nuestro entorno, primero entre los especialistas en la
materia, y seguidamente en el público en general, ha llevado sin embargo a
un cambio notable en la importancia relativa que tiene la química atmosférica
respecto a otros campos de la química ya consolidados. De hecho, en
nuestros días, la química atmosférica es impartida de forma integrada dentro
de estudios superiores como pueden ser las ciencias del medio ambiente y
química general, al menos de forma opcional. Además, el número de grupos
de investigación que centran su actividad exclusivamente en esta área han
multiplicado su número e importancia, paralelamente al desarrollo de nuevos
proyectos que tratan de dar respuesta a problemas de nivel cada vez más
local.
Como consecuencia de la concienciación general que ha ido generándose,
sobre todo entre los ciudadanos de los países más industrializados, incluso las
industrias comienzan a asumir el reto de minimizar el impacto ambiental de
sus emisiones, aunque sea a través de las directivas acordadas en la
Conferencia sobre el Cambio Climático que tuvo lugar en noviembre de 1997
en la ciudad japonesa de Kioto, y que se cerró con el acuerdo de 39 de los
países más industrializados del mundo (entre los que se cuenta España) de
27
Química atmosférica
28
Química atmosférica
80
70
60
50
Altitud (km)
40
30
20
10
0
0 10 20 30 40 50 60 70 80 90 100
% masa total de aire
Figura 2.1: Porcentaje de la masa total de la atmósfera que se encuentra por debajo
de los primeros 80 km de altitud.
29
Química atmosférica
120
100 TERMOSFERA
90
Mesopausa
80
Altitud (km)
70
MESOSFERA
60
Estratopausa
50
40 ESTRATOSFERA
30
Capa de ozono
20
Tropopausa
10 TROPOSFERA
0
-100 -80 -60 -40 -20 0 20 40 60
Temperatura (ºC)
30
Química atmosférica
2.1.1. La exosfera
2.1.2. La termosfera
31
Química atmosférica
32
Química atmosférica
dando lugar a especies excitadas cuya emisión espontánea genera una débil
luminiscencia capaz incluso de iluminar débilmente la superficie terrestre
durante las noches sin luna.
O+ + O2 → O + O2+
O+ + N2 → N + NO+
En la región E, por el contrario, son los iones O2+ y N2+ los productos
primarios de la ionización. La proporción de iones N2+ es, sin embargo, muy
escasa, debido a que se consumen rápidamente en reacciones secundarias
como:
N2+ + O → N* + NO+
N2+ + O2 → N2 + O2+
O2+ + NO → O2 + NO+
33
Química atmosférica
2.1.3. La mesosfera
O3 + O (3P) → 2O2
34
Química atmosférica
150
140
130
120
110 O2
100
90
Altitud (km)
O
80
70
60
50
40 O3
30
20
10
0
10 11 10 12 10 13 10 14 10 15 10 16 10 17 10 18 10 19 10 20 10 21 10 22 10 23 10 24 10 25 10 26
-3
Densidad numérica (m )
35
Química atmosférica
2.1.4. La estratosfera
36
Química atmosférica
[O3] (moléculas/m3)
40
Unidades relativas
35 (ppmv)
30
Altitud (km)
25
20
Tropopausa
15
10
Unidades absolutas
3
(moléculas/m )
5
TROPOSFERA
0
0 1 2 3 4 5 6 7 8
[O3] (ppmv)
37
Química atmosférica
70
65
60
OH
55
Altitud (km)
Clx
50
45
NOx
40
35
30
0 10 20 30 40 50 60 70 80 90 100
% de contribución en la destrucción de O3
Figura 2.6: Contribución de los radicales OH, Clx y NOx a la destrucción de ozono
en la media y alta estratosfera.
38
Química atmosférica
39
Química atmosférica
BrO + ClO → Br + Cl + O2
Br + O3 → BrO + O2
Cl + O3 → ClO + O2
2O3 → 3O2
2.1.5. La troposfera
40
Química atmosférica
41
Química atmosférica
La troposfera libre
42
Química atmosférica
La capa límite
Los veloces vientos reinantes en la troposfera libre son moderados en su
parte más baja a causa de la fricción con la superficie terrestre, dando lugar
como resultado a una capa parcialmente estática de aproximadamente 1 km
de profundidad.
Además de los gases presentes en la parte más alta de la troposfera, la
capa límite contiene gran cantidad de partículas sólidas en suspensión (PM10
y PM2.5)a, además de cientos de compuestos de origen antropogénico que
pueden encontrarse en concentraciones localmente muy elevadas. Sus
concentraciones vienen determinadas por el balance que se establece entre la
emisión, su transformación química y la pérdida por transferencia a la
troposfera libre.
La química de esta delgada capa viene determinada principalmente por
fenómenos meteorológicos como la temperatura, la intensidad lumínica y el
contenido de agua presente en forma de vapor o niebla. Dicha química es tan
compleja que le dedicaremos todo un apartado.
a
Se trata de partículas sólidas en suspensión con diámetros menores a 10 y 2.5 µm
respectivamente. Las de menor tamaño, las PM2.5, son responsables en un alto grado de
muchas de las afecciones pulmonares que afectan a los seres humanos, incluidas las alergias.
43
Química atmosférica
O (1D) + H2O → 2 OH
Otra importante fuente del radical hidroxilo es la reacción del radical HO2
con NO:
HO2 + NO → OH + NO2
⋅CH3 + O2 → CH3OO⋅
44
Química atmosférica
M
OH + NO2 → HNO3
M
HO2 + HO2 → H2O2
El ozono presente en la troposfera es también un importante oxidante,
especialmente para compuestos insaturados como los terpenos. El ozono se
adiciona al doble enlace de estos compuestos formando un anillo de cinco
miembros (-COOOC-) llamado ozonólido, el cual se fragmenta muy
rápidamente dando lugar a radicales de tipo peroxídico, que acaban
evolucionando hacia productos oxigenados estables de una forma similar a
las reacciones representadas más arriba para el caso del metano.
Hacia el final de la década de los 70 se descubrió una nueva faceta de la
química atmosférica, no sólo en las áreas más contaminadas, sino también en
aquellas más alejadas de los núcleos urbanos: la química nocturna del radical
NO3 y del HONO. Fue debido a la aplicación del espectrofotómetro
diferencial de absorción ultravioleta-visible (DOAS) por Platt y
16
colaboradores . Esta novedosa técnica experimental reveló la existencia de
dos importantes especies nitrogenadas durante la noche: el ácido nitroso17 y
el radical NO318 que, de otra manera, en presencia de luz, se fotodisociarían
rápidamente.
45
Química atmosférica
Es éste un proceso que tiene lugar en las áreas urbanas del mundo
industrializado como consecuencia de una compleja serie de reacciones
fotoquímicas.
Se observa cuando las masas de VOCs y óxidos de nitrógeno (NOx = NO
+ NO2) emitidas a la troposfera son atrapadas por una inversión térmica e
irradiadas por la luz solar, condiciones que suceden muy frecuentemente en
la cuenca mediterránea19.
Durante cualquier episodio de contaminación fotoquímica se observan las
siguientes etapas:
46
Química atmosférica
0,50
0,45
0,40 O3
0,35
Concentración (ppm)
0,30
0,25
0,20
0,15
NO2
0,10
0,05 NO
0,00
0:00 2:00 4:00 6:00 8:00 10:00 12:00 14:00 16:00 18:00 20:00 22:00
Tiempo (horas)
O + O2 + M → O3 + M [2]
47
Química atmosférica
NO + O3 → NO2 + O2 [3]
J NO2 [NO 2 ]
[O3 ] =
k [NO ]
48
Química atmosférica
OH + O3 → HO2 + O2
HO2 + O3 → OH + 2O2
NO2 + O3 → NO3 + O2
VOCs + NOx + hν → O3 + …
49
Química atmosférica
La oxidación del SO2 es iniciada por los radicales OH en fase gas, aunque
la oxidación en fase acuosa, presente en la atmósfera en forma de aerosoles,
partículas, nubes y nieblas; es si cabe más importante.
50
Química atmosférica
El NO2 puede dar lugar al ácido nítrico a través de tres mecanismos bien
diferenciados:
OH + NO2 → HNO3
NO3 + RH → R + HNO3
3. La reacción del radical NO3 con NO2 para formar N2O5, el cual
reacciona a continuación con agua en forma de vapor o en estado
líquido generando el ácido nítrico:
51
Química atmosférica
NO2 + O3 → NO3 + O2
52
Química teórica y computacional
53
Química teórica y computacional
54
Química teórica y computacional
“La química cuántica es hoy una ciencia madura que penetra en otras
muchas ramas de la química, como por ejemplo, en el conocimiento de
mecanismos de reacciones enzimáticas (bioquímica) o explicando la relación
existente entre la estructura de un determinado compuesto y su actividad
farmacológica o sus propiedades espectroscópicas”.
55
Química teórica y mecanismos de reacción
57
Química teórica y mecanismos de reacción
Figura 4.1: Diagrama esquemático de un modelo tipo caja mostrando los elementos
básicos que se tienen en cuenta en el mismo.
58
Química teórica y mecanismos de reacción
59
Química teórica y mecanismos de reacción
60
Química cuántica: metodología
61
Química cuántica: metodología
Η̂Φ = ΕΦ
62
Química cuántica: metodología
ˆ = −∑ 1∇ i2 − ∑ 1 ∇ 2A − ∑ ∑ Z A + ∑ ∑ 1 + ∑ ∑ Z A Z B
N M N M N N M M
Η
i =1 2 A =1 2 M A i =1 A =1 riA i =1 j > i rij A =1 B > A R AB
Ε elec = Ε elec (R A )
63
Química cuántica: metodología
Φ( x1 ) = ∑ a i χ i ( x1 )
i
b
Los métodos variacionales son aquellos para los que cualquier función de onda proporciona
un valor esperado de la energía mayor o igual que el valor exacto.
64
Química cuántica: metodología
Φ = ∑ c I ΨI
I
Φ ≅ Ψ0
ˆ Ψ0
Ε 0 = Ψ0 Η
δΕ 0 = 2 δΕ 0 Ηˆ Ψ0
δΕ 0 = 2∑ ∑ c ar Ψar Ηˆ Ψ0 = 0
a r
65
Química cuántica: metodología
Ψar = χ 1 K χ r K χ N
ˆ Ψ0 = 0
Ψar Η
( )
N
fˆ (1) = hˆ(1) + ∑ Jˆ b (1) − Kˆ b (1)
b
ˆ Ψ0 = f ar = 0
Ψar Η
fˆχ a = ε a χ a
66
Química cuántica: metodología
67
Química cuántica: metodología
Ello nos permite resolver las ecuaciones de una forma algebraica, que es
mucho más sencilla a efectos computacionales.
Una base exacta de este tipo debería tener una dimensión i = ∞ , lo que es
imposible a efectos prácticos, por lo que es preciso recurrir a un número
finito de funciones de base. La elección de este tipo de bases constituye uno
de los puntos críticos que pueden condicionar drásticamente la calidad de los
resultados obtenidos. Dos son los principales tipos de funciones de base
utilizados:
68
Química cuántica: metodología
Los conjuntos de base Doble Zeta (DZ) y Triple Zeta (TZ), en los
que el tamaño de la base se aumenta reemplazando cada función por
dos o tres funciones respectivamente.
69
Química cuántica: metodología
70
Química cuántica: metodología
Φ = c 0 Ψ0 + ∑ c ar Ψar + ∑ c ab
rs
Ψabrs + L
ar abrs
71
Química cuántica: metodología
72
Química cuántica: metodología
Ε corr = Ε FCI − Ε 0
O O O O O O
O O O
73
Química cuántica: metodología
enlace y di-radical), nos mostraría que la función de onda del ozono contiene
hasta un 40% de carácter di-radicalario, mientras que la función RHF sólo
incluye un 12% de dicho carácter45. Los métodos que incorporan la
correlación no dinámica se denominan multiconfiguracionales (MCSCF), y
se verán más adelante.
La correlación dinámica, por el contrario, proviene de sustituir el término
N N
1
interelectrónico ∑∑ r
i =1 j > i
del Hamiltoniano de la ecuación de Schrödinger
ij
N N
1
independiente del tiempo por un campo promedio. El término ∑∑ r
i =1 j > i ij
74
Química cuántica: metodología
ΨCISD = c 0 Ψ0 + ∑ c ar Ψar + ∑ c ab
rs
Ψabrs
ar abrs
Ε AB (r → ∞ ) = Ε A + Ε B
75
Química cuántica: metodología
H2 → 2H
Ε CISD (H 2 ) ≠ 2Ε CISD (H )
Ejemplo del cual puede concluirse que el método CISD no debe utilizarse
si antes no se corrige de alguna manera el error de consistencia con el
tamaño. La forma más común de hacerlo es por medio de la corrección de
Davidson, basada en métodos perturbacionales y que presenta la ventaja de
que se añade a posteriori. La expresión de esta corrección es muy simple, y
precisamente por esto su uso está muy generalizado.
Ε corr = Ε CISD
corr + ∆Ε Davidson
( )
donde: ∆Ε Davidson = 1 − c02 Ε CISD
corr
ˆ + Vˆ
ˆ =Η
Η 0
76
Química cuántica: metodología
En dicha teoría, las funciones de orden cero son los determinantes Hartree-
Fock, en tanto que las energías de orden cero son la suma de energías
monoelectrónicas.
La energía de segundo orden, o método MP2 (método que ha sido
ampliamente utilizado en esta tesis para obtener geometrías de partida), viene
dada por medio de la expresión:
2
ˆ Ψabrs
Ψ0 Η
Ε ( 2) = ∑
a <b εa +εb −εr −ε s
r<s
77
Química cuántica: metodología
Ψk = exp Tˆ Ψ0 ()
En dicha exponencial, Ψ0 es el determinante Hartree-Fock y T̂ es el
“operador de cluster”, que puede expandirse en función de los denominados
operadores de excitación.
Siendo:
Tˆ1 = ∑ t ar a r+ a a
a,r
1
Tˆ2 = ∑ t ab rs + +
a r a s ab a a L
2! a ,b
r,s
()
La expansión de la exponencial exp T̂ en serie de Taylor da lugar a la
ecuación:
() 1 1
exp Tˆ = 1 + Tˆ + Tˆ 2 + Tˆ 3 + L
2! 3!
78
Química cuántica: metodología
Tˆ ≅ Tˆ1 + Tˆ2
79
Química cuántica: metodología
t1
T1 =
N
80
Química cuántica: metodología
81
Química cuántica: metodología
ˆ T 2 Ψ0 = ε
Ψ0 Η
ˆ 1 + T1 + T2 + T1T2 Ψ0 = εt ar
Ψar Η
82
Química cuántica: metodología
ˆ 1 + T1 + T2 + 2 −1 T22 Ψ0 = εt ab
Ψabrs Η rs
ˆ T2 Ψ0 + 2 −1 T12 = ε
Ψ0 Η
ˆ 1 + T + T + T T + 2 −1 T 2 + 3!−1 T 3 Ψ = εt r
Ψar Η 1 2 1 2 1 1 0 a
Como conclusión, grosso modo, podemos afirmar que, puesto que todas
las diferencias entre ambos métodos se reducen a los términos en T1 , las
83
Química cuántica: metodología
84
Química cuántica: metodología
c
La opción Full que figura en el programa de cálculo GAUSSIAN implica que se utilizan
también los electrones de core en la correlación. Dicho programa tiene definido, por defecto,
la opción FC (Frozen Core) que no los incluye.
85
Química cuántica: metodología
∆ = ∆( A) + ∆( B)
Donde:
86
Química cuántica: metodología
87
Química cuántica: metodología
88
Química cuántica: metodología
Fˆij = ∑ aˆ i+σ aˆ jσ
σ
aniquilación con σ = α , β .
En función de estos operadores el Hamiltoniano puede expresarse de la
siguiente forma:
ij ij
2 ijkl
ij kl (
ˆ = ∑ h Fˆ + 1 ∑ (ij kl )⋅ Εˆ Εˆ − δ Εˆ
Η jk il )
ij
ˆ Ψ = ∑ h D + ∑ (ij kl )⋅ P
Ε0 = Ψ Η ij ij ijkl
ij ijkl
Dij = Ψ Εˆ ij Ψ = ∑ c m* c n Dijmn
m, n
89
Química cuántica: metodología
90
Química cuántica: metodología
Ψ = c 0 Ψ (0 ) + c1 Ψ (1)
91
Materiales y programas empleados
93
Materiales y programas empleados
94
Materiales y programas empleados
95
Materiales y programas empleados
96
Objetivos
7. OBJETIVOS
Entre los objetivos que nos proponemos abordar los hay unos de carácter
más general y otros más particulares.
97
Objetivos
Estudio de las reacciones del radical H2NO con ozono con el fin de
racionalizar las observaciones registradas durante los experimentos
llevados a cabo sobre la cinética de la reacción NH2 + O3.
98
II. RESULTADOS
Resultados
Fernando de Ávalos
Marqués de Pescara
101
Resultados
Como puede verse en el apéndice de esta tesis, cada una de las reacciones
estudiadas se corresponde en líneas generales con uno de los artículos (I-VI)
que se muestran a continuación:
F + O3 → FO + O2 (artículo IV)
NO + O3 → NO2 + O2 (artículo V)
102
La reacción NH2 + O3
d
An Ab Initio Study on the Mechanism of the Atmospheric Reaction NH2 + O3 → H2NO +
O2, Julio Peiró-García, Ignacio Nebot-Gil y Manuela Merchán, ChemPhysChem, 4, 2003,
366-372.
103
La reacción NH2 + O3
104
La reacción NH2 + O3
105
La reacción NH2 + O3
106
La reacción NH2 + O3
107
La reacción NH2 + O3
NH2
r(N-H) 1.024 1.019 1.029 1.024
∠(H-N-H) 102.8 103.5 102.5 103.3
O3
r(O-O) 1.281 1.274 1.256 1.272
∠(O-O-O) 117.1 116.9 117.9 116.8
H2NOd
r(N-H) 1.014 1.011 1.015 -
r(N-O) 1.262 1.255 1.281 -
∠(H-N-O) 120.1 120.2 118.4 -
O2
r(O-O) 1.223 1.215 1.203 1.207
a
Las distancias de enlace se dan en Ángstroms, los ángulos en grados.
b
Las optimizaciones a nivel UMP2 han sido efectuadas teniendo en cuenta todos los
electrones en la correlación.
c
Datos recopilados de las referencias92-95.
d
La geometría del radical H2NO tiene simetría C2v a nivel UMP2, sin embargo, a
nivel QCISD se encuentra que el ángulo diedro ∠(H-H-N-O) es de 155º.
108
La reacción NH2 + O3
→ NO3 + H2 (8.9)
→ HNO2 + OH (8.10)
109
La reacción NH2 + O3
Para que esto suceda así, al menos una de las dos reacciones, (8.12) ó
(8.13), debe ser suficientemente rápida.
Si bien hasta la fecha no hay ninguna evidencia experimental que
confirme la existencia de alguna de estas dos reacciones, nos limitaremos a
suponer, por el momento, que el producto único de la reacción NH2 + O3 es
el radical H2NO. Esta hipótesis viene respaldada por el hecho de que no se
110
La reacción NH2 + O3
han encontrado otras vías de reacción, aunque esto no implica que no puedan
existir. Dejaremos, en cualquier caso, para más adelante (artículo II, capítulo
siguiente), la constatación de la hipótesis que acabamos de exponer.
Como puede observarse en la tabla 8.1, la simetría del radical H2NO a
nivel MP2 es C2v, mientras que a nivel QCISD nos encontramos con que el
diedro ∠(H-H-N-O) es de 155º y, por tanto, que la simetría es Cs. Dicha
contradicción aparente ha suscitado que la estructura del radical H2NO haya
sido objeto, hasta fechas muy recientes, de varios estudios teóricos que han
dictaminado estructuras con simetría C2v o Cs dependiendo del nivel de
cálculo empleado96-98.
El estudio sistemático llevado a cabo por Komaromi y Tronchet99 en 1993
utilizando todas las opciones posibles del paquete GAUSSIAN, sin embargo,
concluye que la búsqueda de una estructura exacta carece de importancia, ya
que la barrera existente entre ambas conformaciones es menor que la energía
correspondiente al primer nivel vibracional de dicha molécula.
En la tabla 8.2 se detallan los modos normales de vibración de los
reactivos y productos de la reacción (8.1), tanto a nivel MP2 como QCISD.
La concordancia con los resultados experimentales es en general
razonablemente correcta. Sin embargo, podemos observar que el método
MP2 no describe adecuadamente el modo normal de vibración b2
correspondiente a la molécula de ozono con los 2285 cm-1 (1135 cm-1
medidos experimentalmente) que predice para la base 6-311+G(d,p). Esta
discrepancia se puede atribuir a las limitaciones intrínsecas del método MP2,
de tipo monorreferencial, cuando se aplica a un sistema como el ozono, que
precisa de hasta tres configuraciones para describir correctamente su función
111
La reacción NH2 + O3
NH2
MP2/6-311+G(d,p) 1520.2, 3461.7, 3564.5
MP2/6-311+G(3df,2p) 1542.6, 3480.7, 3581.6
QCISD/6-311+G(d,p) 1527.5, 3389.8, 3481.5
Exp.a 1497, 3221, 3300
O3
MP2/6-311+G(d,p) 751.0, 1165.0, 2285.3
MP2/6-311+G(3df,2p) 750.8, 1176.9, 2195.9
QCISD/6-311+G(d,p) 749.3, 935.2, 1240.9
Exp. 716, 1089, 1135
H2NO
MP2/6-311+G(d,p) 229.5, 1300.9, 1515.9, 1712.5, 3504.7, 3652.4
MP2/6-311+G(3df,2p) 288.5, 1288.8, 1533.7, 1702.0, 3509.3, 3656.3
QCISD/6-311+G(d,p) 441.3, 1310.5, 1369.5, 1699.5, 3489.7, 3624.3
O2
MP2/6-311+G(d,p) 1461.0
MP2/6-311+G(3df,2p) 1502.3
QCISD/6-311+G(d,p) 1649.2
Exp. 1580
a
Valores experimentales extraídos de las referencias92,93,95.
112
La reacción NH2 + O3
113
La reacción NH2 + O3
114
La reacción NH2 + O3
Los resultados obtenidos son muy similares a los que Hwang y Mebel
calcularon para la reacción Cl + O3102. De hecho, el estado de transición
determinante de dicha reacción viene caracterizado por una frecuencia
vibracional de 242i cm-1 a nivel MP2/6-311G(d) y de 188i cm-1 con esta
misma base pero a nivel QCISD. Aunque en dicho trabajo se optimiza este
mismo estado de transición análogo al TS1 nuestro a nivel CASSCF(17,12),
éste parece no haber sido caracterizado y, por tanto, no se da la frecuencia
correspondiente a la transición.
Por último, indicaremos que el valor obtenido para 〈 S 2 〉 a nivel QCISD
es de 0.833, lo que queda muy lejos del límite de S + 1 que Cremer y Yuan64
tomaron como valor de diagnóstico indicativo de la incapacidad de un
método basado en la teoría Coupled Cluster para describir adecuadamente un
sistema concreto debido a su carácter excesivamente multiconfiguracional.
115
La reacción NH2 + O3
e
El IRC es un método capaz de seguir el camino de reacción dada la estructura de un estado
de transición cualquiera. Por lo tanto, nos puede conducir fácilmente a los reactivos y
productos de la reacción en que interviene.
f
En realidad en este caso a través de dos intermedios, uno de los cuales, con una barrera
insignificante y que no observamos nosotros, corresponde al tránsito del intermedio
[ClOOO] de simetría C1 a una simetría Cs.
116
La reacción NH2 + O3
dos etapas, distinto del mecanismo directo que indican los métodos QCISD y
CASSCF.
Lógicamente, el método MP2 nos ha conducido a un resultado erróneo,
debido a que la función de referencia de orden cero (UHF) en la que se basa
dicho método no debe ser una aproximación “suficientemente buena”. La
razón por la que el método QCI, que también se basa en una función de
referencia UHF, sí proporciona resultados satisfactorios (en concordancia con
el método CASSCF) se debe a que introduce correcciones a orden infinito,
por lo que tolera mejor que los métodos Møller-Plesset una función de orden
cero inadecuada52.
El buen comportamiento de las metodologías basadas en la teoría Coupled
Cluster contrastadas con una metodología de tipo multiconfiguracional (la
única que en principio puede describir correctamente el carácter
multirreferencial de la molécula de ozono), resulta muy interesante, ya que
nos permitirá comparar esta reacción con el resto de reacciones que
constituyen el conjunto de la presente tesis.
QCISD(T)/6-311+G(d,p)
QCISD(TQ)/6-311+G(d,p)
CCSD(T)/6-311+G(3df,2p)
117
La reacción NH2 + O3
CCSD(Full,T)/6-311+G(3df,2p)
CASPT2/4s3p1d|2s1p y 6-311+G(d,p)
Tabla 8.4: Entalpía de reacción y barreras energéticas calculadas para los estados de
transición TS1 y TS2 y el intermedio Int de la reacción (8.1).
118
La reacción NH2 + O3
119
La reacción NH2 + O3
incluso menor que cero cuando se calcula por medio de otra metodología más
sofisticada, en este caso la QCISD(T).
∆E (MP2) > 0
∆E (QCISD) < 0
Energía
MP2
QCISD(T)
120
La reacción NH2 + O3
121
La reacción NH2 + O3
122
La reacción NH2 + O3
8.3.2. El diagnóstico T1
Tabla 8.5: Valores del diagnóstico T1 para todos los puntos estacionarios de la
reacción (8.1).
CCSD/6-311G(3df,2p)
Estructura MP2/6-311+G(3df,2p) QCISD/6-311+G(d,p)
NH2 0.0079 0.0081
O3 0.0272 0.0256
TS1 0.0476 0.0472
Int 0.0338 -
TS2 0.0494 -
H2NO 0.0200 0.0228
O2 0.0147 0.0144
123
La reacción NH2 + O3
una función de referencia HF de capa abierta para el que los valores de corte,
según Jayatilaka y Lee55, pueden ser superiores incluso a 0.045. Es decir, que
nos encontramos precisamente en el límite de calidad aceptable. De hecho,
Rienstra-Kiracofe y colaboradores56 encuentran aceptables valores de 0.045,
incluso superiores, para sistemas similares. Para llegar a dicha conclusión
comparan, como nosotros, los resultados Coupled Cluster con los obtenidos
por medio de métodos multiconfiguracionales.
La inclusión de las triples excitaciones, así como la capacidad que tienen
los métodos Coupled Cluster de incluir todas las correcciones de un tipo
hasta orden infinito, hacen, por tanto, de esta metodología una alternativa
muy útil a las de tipo CAS, especialmente en sistemas de mayor tamaño
molecular.
124
La reacción NH2 + O3
QCISD(TQ) son francamente precisos, lo que nos lleva a pensar que el valor
propuesto en base a métodos ab initio en este trabajo puede ser una
estimación bastante razonable de la entalpía de reacción, ∆H r0 .
A modo de resumen final, en la figura 8.3, se representa el perfil de
energía CCSD(T) correspondiente a la reacción (8.1).
5.9
NH2 + O3
TS1
H2NO + O2 -65.2
8.4. Conclusiones
125
La reacción NH2 + O3
126
La reacción H2NO + O3
9. LA REACCIÓN H2NO + O3
127
La reacción H2NO + O3
g
A Theoretical Ab Initio Study on the H2NO + O3 Reaction, Julio Peiró-García, Víctor M.
Ramírez-Ramírez e Ignacio Nebot-Gil, Journal of Computational Chemistry, 24, 2003, 1321-
1327.
128
La reacción H2NO + O3
h
No se han optimizado estructuras a nivel MP2/6-311+G(3df,2p) por dos razones: porque no
aporta ninguna mejora sustancial, y porque el tamaño molecular mayor de este nuevo sistema
implica tiempos de cálculo fuera de nuestras posibilidades.
129
La reacción H2NO + O3
UMP2b QCISD
Parámetrosa Experimentalc
6-311+G(d,p) 6-311+G(d,p)
HNO
r(N-H) 1.053 1.058 1.062
r(N-O) 1.220 1.207 1.221
∠(H-N-O) 107.9 108.4 108.5
ν (cm-1) 1504, 1583, 3039 1580, 1650, 2963 1501, 1565, 2684
OH
r(O-H) 0.968 0.971 0.970
ν (cm-1) 3338.5 3776.0 3739.9
a
Las distancias de enlace se dan en Ángstroms, los ángulos en grados.
b
Las optimizaciones a nivel UMP2 han sido efectuadas teniendo en cuenta todos
los electrones en la correlación (opción Full).
c
Datos recopilados de las referencias107-109.
130
La reacción H2NO + O3
131
La reacción H2NO + O3
132
La reacción H2NO + O3
133
La reacción H2NO + O3
TS2
828i, 75, 138, 344, 435, 580, 688, 727, 1148, 1260,
MP2a
1350, 1361, 1708, 3498, 3650
H2NOO
QCISD 150, 563, 843, 1066, 1247, 1327, 1650, 3464, 3563
MP2 146, 580, 840, 1211, 1251, 1312, 1636, 3476, 3587
TS4
QCISD 583i, 89, 368, 815, 956, 1299, 1509, 3408, 3511
MP2 1068i, 82, 428, 945, 1076, 1507, 1592, 3461, 3578
B
QCISD 132, 519, 563, 876, 1102, 1421, 1480, 3410, 3785
MP2 136, 555, 575, 957, 1156, 1401, 1468, 3456, 3787
TS5
QCISD 2436i, 562, 807, 955, 1019, 1333, 1718, 2241, 3497
MP2 2605i, 637, 789, 966, 1051, 1435, 1615, 2257, 3523
TS6
QCISD 795i, 178, 460, 689, 1094, 1183, 1445, 3302, 3812
MP2 1222i, 90, 483, 744, 1109, 1123, 1472, 3316, 3857
a
El estado de transición TS2 no se ha caracterizado a nivel QCISD debido al
elevado coste computacional de este cálculo, que lo deja fuera de la capacidad
de las mejores máquinas de que disponíamos a finales de 2002. Todos lo
cálculos han sido realizados con la base 6-311+G(d,p).
134
La reacción H2NO + O3
135
La reacción H2NO + O3
QCISD(T)/6-311+G(d,p)//MP2/6-311+G(d,p)
QCISD(T)/6-311+G(2d,2p)//MP2/6-311+G(d,p)
CCSD(T)/6-311+G(3df,2p)//QCISD/6-311+G(d,p)
Gaussian-3 (G3)
136
La reacción H2NO + O3
137
La reacción H2NO + O3
TS 3 TS 4
H2NO + O3 → H2NOO + O2 → NH2 + 2O2
138
La reacción H2NO + O3
Figura 9.4: Perfil de energías QCISD(T) para las reacciones (8.1) y (9.4) a (9.6). Los
puntos de llegada de las reacciones (8.1) y (9.4) se han tomado como los puntos de
partida de las reacciones (9.4) y (9.5) y (9.6) respectivamente.
139
La reacción H2NO + O3
9.3. Conclusiones
140
La reacción OH + O3
OH + O3 → HO2 + O2 (10.2)
i
Ab Initio Study on the Mechanism of the Atmospheric Reaction OH + O3 → HO2 + O2,
Julio Peiró-García e Ignacio Nebot-Gil, ChemPhysChem, 4, 2003, 843-847.
141
La reacción OH + O3
HO2
r(O-H) 0.972 0.971 0.971 0.971
r(O-O) 1.310 1.303 1.335 1.331
∠(H-O-O) 105.3 105.1 104.7 104.3
a
Las distancias de enlace se dan en Ángstroms, los ángulos en grados.
b
Las optimizaciones a nivel UMP2 han sido efectuadas con la opción Full.
c
Datos recopilados de las referencia124.
142
La reacción OH + O3
143
La reacción OH + O3
como de la entalpía de reacción, con los que contrastar los valores calculados
de forma teórica.
El mecanismo de reacción observado es análogo al de la reacción
isoelectrónica NH2 + O3, tal como queda representado en la figura 10.1.
1.281
1.274 1.328
1.256 1.307
1.321
1.223
1.215
1.830 1.203
1.888
1.971
0.968
0.967
0.971
0.968 1.310
0.966 1.303
0.971 1.335
0.972
0.971
0.971
144
La reacción OH + O3
QCISD(T)/6-311+G(d,p)
QCISD(TQ)/6-311+G(d,p)
CCSD(T)/6-311+G(3df,2p)
Como podemos ver en la tabla 10.2, estamos ante un sistema con una
componente multiconfiguracional importante que el método MP2 no puede
resolver correctamente, por lo que a este nivel de cálculo se obtienen barreras
de transición exageradamente altas, en torno a las 39 kcal mol-1. La
145
La reacción OH + O3
Tabla 10.2: Entalpía de reacción y barreras energéticas calculadas para los estados
de transición TS1 y TS2 y el intermedio Int correspondientes a la reacción (10.1).
146
La reacción OH + O3
TS1
Exp: 2.0
OH + O3 4.7
(2.2)
-43.7
(-45.0)
HO2 + O2
Como podemos ver en la figura 10.1, las longitudes de los enlaces O(1)-
O(2), O(2)-O(3) y O(4)-H(5) correspondientes al estado de transición TS1,
147
La reacción OH + O3
148
La reacción OH + O3
149
La reacción OH + O3
Hasta aquí, por tanto, tenemos que los resultados obtenidos concuerdan
razonablemente bien; por una parte, con los obtenidos por Back y Yokota106 y
Hack y colaboradores89 en cuanto a la reactividad comparada entre las
reacciones (8.1) y (10.1); y por la otra, con los estudios experimentales
realizados sobre el sistema OH + O3, que proporcionan energías de activación
tan sólo 2 – 3 kcal mol-1 por debajo de nuestras predicciones.
150
La reacción OH + O3
10.3.3. El diagnóstico T1
Por último, nos queda cerciorarnos de hasta qué punto las especies
involucradas en el mecanismo de la reacción (10.1) están o no bien descritas
por medio de una función Coupled Cluster basada en una función de
151
La reacción OH + O3
referencia de tipo UHF. Con tal fin, como se hizo previamente con el sistema
NH2 + O3, se ha calculado el valor del diagnóstico T1 para cada uno de los
puntos estacionarios presentes en la superficie de energía potencial. Los
resultados obtenidos se muestran en la tabla 10.4.
Tabla 10.4: Valores del diagnóstico T1 para todos los puntos estacionarios
involucrados en el mecanismo de la reacción (10.1).
CCSD/6-311G(3df,2p)
Estructura MP2/6-311+G(3df,2p) QCISD/6-311+G(d,p)
OH 0.0066 0.0067
O3 0.0272 0.0256
TS1 0.0481 0.0475
HO2 0.0265 0.0286
O2 0.0147 0.0144
10.4. Conclusiones
152
La reacción OH + O3
153
La reacción OH + O3
154
La reacción F + O3
11. LA REACCIÓN F + O3 → FO + O2
155
La reacción F + O3
reacción de los átomos de flúor con el ozono aislado sobre una matriz sólida
de argón.
El interés suscitado por este reciente descubrimiento llevó a Li y
colaboradores137, tan sólo tres años después, a efectuar un estudio teórico ab
initio con el método Gaussian-3 (G3) sobre la reacción F + O3 con el que
confirmaron en parte los resultados obtenidos por Misochko y
136
colaboradores . Desgraciadamente, este interesante trabajo obvia por
completo la naturaleza multiconfiguracional de alguna de las especies
involucradas en la reacción que denominaremos (11.1) y las estructuras sólo
se han optimizado al bajo nivel de cálculo MP2/6-31G*, el cual hemos
demostrado con los tres capítulos anteriores que, en ocasiones, cuando el
carácter multiconfiguracional de las estructuras involucradas en el
mecanismo de reacción es importante, puede conducir a resultados
insatisfactorios.
j
An ab initio study on the mechanism of the F + O3 → FO + O2: comparative reactivity study
along the isoelectronic NH2, OH, and F radical series, Julio Peiró-García e Ignacio Nebot-
Gil, Chemical Physic Letters, 391, 2004, 195-199.
156
La reacción F + O3
UMP2b QCISD
Parámetrosa Experimentalc
6-311+G(d,p) 6-311+G(d,p)
FO
r(F-O) 1.331 1.360 1.358
a
Las distancias de enlace se dan en Ángstroms.
b
Las optimizaciones a nivel UMP2 han sido efectuadas teniendo en cuenta
todos los electrones en la correlación.
c
Datos recopilados de la referencia139.
157
La reacción F + O3
158
La reacción F + O3
activo seleccionado para la reacción (11.1) elimina los electrones de core del
sistema, así como los orbitales y electrones 2s correspondientes a los átomos
de oxígeno y flúor, quedando, en consecuencia, 17 electrones y 12 orbitales
moleculares de tipo p.
MP2/6-311+G(d,p) 1.281
QCISD/6-311+G(d,p) 1.296
CASSCF(17,12)/4s3p1d 1.306
TS1
1.317
1.281 1.310
1.256 1.360
1.223
1.203
1.678
1.880
1.840
F
1.331
1.360
159
La reacción F + O3
160
La reacción F + O3
161
La reacción F + O3
162
La reacción F + O3
QCISD/6-311+G** 1.269
FO3
1.468
1.404
163
La reacción F + O3
Tabla 11.3: Frecuencias de vibración del intermedio de reacción FO3 y del estado de
transición TS2 correspondiente a la reacción (11.1).
FO3
MP2/6-311+G(d,p) 127.4, 486.7, 669.3, 1144.8, 1706.5, 2200.3
QCISD/6-311+G(d,p) 96.8, 402.9, 588.5, 716.3, 1069.3, 1154.6
TS2
MP2/6-311+G(d,p) 1379.1i, 103.4, 307.4, 759.4, 1583.5, 1685.2
QCISD/6-311+G(d,p) 443.1i, 81.6, 376.3, 620.0, 1106.0, 1237.1
164
La reacción F + O3
165
La reacción F + O3
166
La reacción F + O3
167
La reacción F + O3
11.4. El diagnóstico T1
Como en casos anteriores, y con el fin de analizar hasta qué grado los
resultados Coupled Cluster obtenidos son fiables, se ha calculado el
diagnóstico T1 para cada una de las especies moleculares involucradas en el
mecanismo de la reacción (11.1). Los resultados obtenidos se muestran en la
tabla 11.5.
Tabla 11.5: Valores del diagnóstico T1 para todos los puntos estacionarios
involucrados en el mecanismo QCISD de la reacción (11.1).
CCSD/6-311G(3df,2p)
Estructura QCISD/6-311+G(d,p)
F 0.0040
O3 0.0256
TS1 0.0497
Int 0.0399
TS2 0.0531
FO 0.0266
O2 0.0144
168
La reacción F + O3
169
La reacción F + O3
Este valor, además, está de acuerdo con los valores sugeridos por Li y
colaboradores137 a los niveles de cálculo G3 y G3MP2 (con -26.6 y -27.0 kcal
mol-1 respectivamente) que proporcionan buenos resultados en este caso.
170
La reacción F + O3
171
La reacción F + O3
172
La reacción F + O3
11.7. Conclusiones
173
La reacción F + O3
174
La reacción F + O3
175
La reacción NO2 + O3
O + O2 →
M
O3 (12.2)
NO + O3 → NO2 + O2 (12.3)
k
Ab Initio Study of the Mechanism of the Atmospheric Reaction: NO2 + O3 → NO3 + O2,
Julio Peiró-García e Ignacio Nebot-Gil, Journal of Computational Chemistry, 24, 2003,
1657-1663.
177
La reacción NO2 + O3
178
La reacción NO2 + O3
NO + O3 → NO2 + O2 (12.11)
179
La reacción NO2 + O3
180
La reacción NO2 + O3
NO2
r(N-O) 1.216 1.202 1.209 1.193
∠(O-N-O) 133.7 134.0 134.2 134.1
NO3
r(N-O) 1.255 1.243 1.241 1.240
a
Las distancias de enlace se dan en Ángstroms, los ángulos en grados.
b
Las optimizaciones a nivel UMP2 han sido efectuadas teniendo en cuenta
todos los electrones en la correlación.
c
Datos recopilados de las referencias163,164.
181
La reacción NO2 + O3
energéticos entre los dos tipos de estructuras son mínimas, lo que está de
acuerdo con que sólo ha de encontrarse una única solución para el estado
fundamental electrónico del radical nitrato.
En la tabla 12.2, además, se detallan las frecuencias calculadas a nivel
MP2 y QCISD para las moléculas de NO2 y NO3.
En dicha tabla podemos encontrar dos discrepancias a nivel MP2 respecto
de los valores determinados experimentalmente. Se trata del “stretching”
asimétrico correspondiente al radical NO3 y el modo normal de vibración b2
de la molécula de NO2. Ambas anomalías se deben a la contaminación de
espín que presenta la función de referencia UHF173 y que no resuelve el
método MP2. Estas discrepancias sí que son resueltas a nivel QCISD debido
a que este método introduce todas las correcciones de un determinado tipo
hasta orden infinito.
Tabla 12.2: Modos normales de vibración de los radicales NO2 y NO3 a nivel MP2 y
QCISD.
NO2
MP2/6-31G(d) 752.7, 1382.4, 2290.9
MP2/6-311G(d) 775.9, 1390.3, 2510.1
QCISD/6-31G(d) 747.1, 1380.0, 1632.3
Experimentala 750, 1320, 1617
NO3
MP2/6-31G(d) 719.5, 745.2, 1084.3, 2366.8
MP2/6-311G(d) 743.1, 768.5, 1096.7, 2422.1
QCISD/6-31G(d) 125.1, 787.5, 1134.6, 1163.0
Experimental 704, 765, 940, 1158
a
Valores experimentales extraídos de las referencias174,175.
182
La reacción NO2 + O3
TS1
1.296 MP2/6-31G(d)
1.311 QCISD/6-31G(d)
1.351
1.369
1.207 1.940
1.200 2.015
1.210
1.202
183
La reacción NO2 + O3
184
La reacción NO2 + O3
185
La reacción NO2 + O3
Superficie QCISD(T)
10,0
TS1
5,0
0,0
Región
Energía (kcal mol )
-1
correspondiente Región
-5,0
al intermedio A correspondiente
(MP2) al TS2 (MP2)
-10,0
-15,0
-20,0
-25,0
-30,0
Coordenada de reacción
186
La reacción NO2 + O3
CCSD(T)/6-311G(d)
CCSD(T)/6-311G(2d)
CASSCF/6-31G(d)l
CASPT2/6-31G(d)
G2 y G2Qm
l
Se han empleado concretamente los siguientes espacios activos: (15,11), (17,12) y (19,13),
este último en el límite de la capacidad del programa MOLCAS 5.0.
187
La reacción NO2 + O3
m
La técnica G2Q es idéntica a la G2 (Gaussian-2), excepto que emplea una estructura de
partida optimizada al nivel de cálculo QCISD/6-31G(d).
188
La reacción NO2 + O3
189
La reacción NO2 + O3
190
La reacción NO2 + O3
191
La reacción NO2 + O3
MP2 son unas 2 kcal mol-1 menores, debido precisamente al error que
introduce el modo normal de vibración b2 calculado a este nivel para la
molécula de ozono.
A modo de resumen gráfico final, en la figura 12.3, se representa el perfil
de energía CCSD(T)/6-311G(2d) basado en estructuras optimizadas a nivel
MP2/6-31G(d) correspondiente a la reacción (12.5).
12.4. El diagnóstico T1
Tabla 12.5: Valores del diagnóstico T1 para todos los puntos estacionarios
involucrados en el mecanismo de la reacción (12.5).
CCSD/6-311G(2d)
Estructura MP2/6-31G(d) MP2/6-311G(d)
192
La reacción NO2 + O3
193
La reacción NO2 + O3
k B T (∆G ≠ 0 / RT )
k r (T ) = e
hc 0
194
La reacción NO2 + O3
12.7. Conclusiones
Las conclusiones que pueden deducirse del presente capítulo son muy
similares a las expuestas en los capítulos 8 al 11, por lo que las resumiremos
lo más brevemente posible a continuación:
195
La reacción NO2 + O3
196
La reacción NO2 + O3
O N O N O N
O O O
197
La reacción NO + O3
14
Ab Initio Study of the Mechanism and Thermochemistry of the Atmospheric Reaction: NO
+ O3 → NO2 + O2, Julio Peiró-García e Ignacio Nebot-Gil, Journal of Physical Chemistry A,
106, 2002, 10302-10310.
199
La reacción NO + O3
O + O2 →
M
O3 (13.3)
200
La reacción NO + O3
NO + O3 → NO2 + O2 (13.1)
NO + O3 → NO2 + O2 (13.1)
NO2 + O → NO + O2 (13.5)
O3 + O → 2O2 (13.6)
201
La reacción NO + O3
202
La reacción NO + O3
203
La reacción NO + O3
NO + O3 → ONOO + O (13.9)
ONOO → NO + O2 (13.10)
15
Los métodos semiempíricos se caracterizan por el uso de parámetros obtenidos a partir de
datos experimentales con el fin de simplificar el cálculo. Resuelven una forma aproximada
de la ecuación de Schrödinger que depende de la idoneidad de los parámetros seleccionados
para el sistema químico en estudio.
204
La reacción NO + O3
205
La reacción NO + O3
16
Como en el caso de la metodología MP2, se ha utilizado la opción full que implica incluir
también los electrones de core en la correlación.
206
La reacción NO + O3
NO
r(N-O) 1.134 1.137 1.175 1.151
O3
r(O-O) 1.281 1.250 1.276 1.278
∠(O-O-O) 116.9 117.8 117.4 116.8
NO2
r(N-O) 1.202 1.194 1.209 1.193
∠(O-N-O) 134.0 134.5 134.2 134.1
O2
r(O-O) 1.223 1.204 1.222 1.207
a
Las distancias de enlace se dan en Ángstroms, los ángulos en grados.
b
Las optimizaciones a nivel MP2 y MP4 han sido efectuadas teniendo en cuenta
todos los electrones en la correlación.
c
Datos nuevos recopilados de la referencia213.
Tabla 13.2: Modos normales de vibración del radical NO a nivel MP2, MP4 y
QCISD.
MP2/6-311G(d) 3534.0
MP4/6-31G(d) 3428.0
QCISD/6-31G(d) 1849.8
Experimentala 1876
a
Valor experimental extraído de la referencia213.
207
La reacción NO + O3
208
La reacción NO + O3
1.134 MP2/6-31G*
TS1 1.159 QCISD/6-31G*
1.996
2.096
1.370
1.364 1.290
1.310
209
La reacción NO + O3
210
La reacción NO + O3
cis-TS1 1.137
1.163
2.041 MP2/6-31G*
2.156 QCISD/6-31G*
1.296
1.316
1.355
1.359
cis-TS1
MP2/6-31G(d) 333.7i, 118.0, 154.8, 228.5, 449.1, 667.0, 1060.3,
1138.6, 3643.1
QCISD/6-31G(d) 303.3i, 102.7, 129.7, 213.8, 379.3, 632.5, 818.7,
1008.1, 1831.0
trans-TS1
MP2/6-31G(d) 352.9i, 73.6, 125.4, 227.9, 345.2, 667.7, 1101.3,
1171.4, 3774.8
QCISD/6-31G(d) 363.8i, 55.3, 99.2, 185.9, 329.2, 632.4, 787.7,
1039.1, 1861.4
211
La reacción NO + O3
212
La reacción NO + O3
213
La reacción NO + O3
214
La reacción NO + O3
TS1
cis-TS1
NO + O3
3.02 y 3.79
Experimental: 2.72
NO2 + O2 -54.13
215
La reacción NO + O3
216
La reacción NO + O3
13.5. El diagnóstico T1
Tabla 13.5: Valores del diagnóstico T1 para todos los puntos estacionarios
involucrados en el mecanismo de la reacción (13.1).
CCSD(T)/6-311G(d)
Estructura MP2/6-31G(d) MP2/6-311G(df)
NO 0.0261 0.0204
O3 0.0349 0.0274
TS1 0.0501 0.0525
NO2 0.0275 0.0241
O2 0.0152 0.0148
Este valor, eso sí, supera con creces a los obtenidos en el resto de
sistemas, lo que, unido a la quimiluminiscencia que presenta la reacción
(13.1), revela el interés que podría tener un futuro estudio CASSCF/CASPT2
de este sistema.
217
La reacción NO + O3
Como queda plasmado en el artículo VI, se han estudiado además las vías
de reacción (13.8) a (13.12), cuyos mecanismos de reacción se representan en
las figuras 13.4 y 13.5.
Estos cálculos se basan exclusivamente en la metodología MP2, que,
posteriormente a la publicación del artículo VI, comenzó a revelar las
carencias que ya se han puesto de manifiesto previamente. Todas las vías de
reacción, además, parten de la estructura del intermedio de reacción MP2
[ONOOO], denominado A, que no aparece en la superficie de energía
potencial QCISD, lo que nos conduce a plantearnos la existencia real de las
mismas.
218
La reacción NO + O3
219
La reacción NO + O3
que especies como el radical ONOO puedan existir el tiempo suficiente como
para poder ser observados.
13.7. Conclusiones
Las conclusiones que pueden deducirse del presente capítulo son en parte
idénticas a las encontradas en los capítulos 8 al 12, pero también las hay
diferentes debido al distinto enfoque que se ha dado al sistema NO + O3. En
cualquier caso, las resumimos todas lo más brevemente posible a
continuación:
220
La reacción NO + O3
221
III. CONCLUSIONES
Conclusiones generales
225
Conclusiones generales
226
Conclusiones generales
227
Conclusiones generales
228
Conclusiones generales
229
Conclusiones generales
conduce a que sean muy bajas. Esto se debe a que la estructura MP2
del TS1 en la que se basan dichos métodos está sensiblemente
desplazada hacia los productos, pero también al error que introduce la
ZPC calculada a nivel MP2 (aproximadamente -2 kcal mol-1).
230
Conclusiones generales
231
IV. BIBLIOGRAFÍA
Bibliografía
15. BIBLIOGRAFÍA
235
Bibliografía
236
Bibliografía
237
Bibliografía
238
Bibliografía
239
Bibliografía
240
Bibliografía
241
Bibliografía
242
Bibliografía
243
Bibliografía
244
Bibliografía
[105] Evaluated Kinetic and Photochemical Data, J. Phys. Chem. Ref. Data,
26 (1997) 1396.
[114] J. G. Anderson and F. Kaufman, Chem. Phys. Lett., 19, (1973) 483.
245
Bibliografía
246
Bibliografía
247
Bibliografía
[137] L-C. Li, J. Wang, X. Wang, A-M. Tian, N-B. Wong, Int. J. Quantum
Chem. 2002, 87, 288.
[139] Handbook of Chemistry and Physics, CRC Press LLC (D. R. Lide
Editor-in-chief), (1999-2000).
248
Bibliografía
249
Bibliografía
[163] D. R. Cide (Ed.) Handbook of Chemistry and Physics, 74th edn. CRL
Press, Boca Raton, FL (1993), reprinted from: Kagaku Benran, 3rd
edn, vol. II, (1984) 649.
[169] D. Wang, P. Jiang, X. Qian, and G. Hong, J. Chem. Phys., 106 (1997)
3003.
250
Bibliografía
251
Bibliografía
252
Bibliografía
253
Bibliografía
254
V. APÉNDICES
Acrónimos
257
Acrónimos
258
Acrónimos
259
Artículos publicados
Lista de publicaciones
261
Artículos publicados
262
Artículos publicados
263
ARTÍCULO I