Tab<strong>la</strong> 3.3: Energías y principales parámetros geométricos (distancias en Å y ángulos en grados) para los puntos estacionariosinvolucrados en <strong>la</strong>s transposición sigmatrópica [1,7] <strong>de</strong> hidrógeno en el 7-metilocta-1,3(Z),5(Z)-trieno calcu<strong>la</strong>dos al nivel AM1-SRPEstructura Energía φ a 1 φ a 2 φ a 3 α1 a α2 a α3 a d(C 7 − H 17 ) d(C 1 − C 7 ) d(C 1 − H 17 )R1 0.00 174.07 170.0 -32.3 108.5 114.4 111.0 1.113 6.373 5.990R2/R2* 1.82 -179.8 -141.8 166.5 107.0 114.4 111.7 1.115 6.177 7.040R3 2.21 -180. 180. -180.0 106.9 114.4 112.1 1.115 6.527 7.495R4/R4* 1.99 -32.6 145.2 -166.8 107.0 114.5 111.7 1.115 4.921 5.607R6/R6* 0.13 -35.9 153.5 48.7 107.6 114.5 111.0 1.114 5.075 5.395R7/R7* 0.22 -36.2 153.4 -28.7 108.6 114.5 111.1 1.112 5.061 4.822R8/R8* 2.23 -178.6 91.9 154.6 107.1 114.4 111.1 1.113 5.028 6.097R9/R9* 2.03 -178.2 116.0 123.3 106.0 114.4 110.7 1.116 5.698 6.741R10/R10* 0.52 -179.4 63.5 -11.2 108.7 114.4 111.0 1.114 4.271 3.458R11/R11* 2.43 -25.6 125.0 122.1 106.1 114.4 110.6 1.115 4.511 5.517R12/R12* 1.78 14.8 86.3 15.9 108.5 114.4 111.3 1.113 3.836 3.899TS12/TS12* 24.40 9.1 27.5 -72.4 100.0 113.6 114.0 1.355 2.568 1.342P12/P12* -6.32 -2.3 2.5 -134.6 113.6 108.3 115.3 2.506 3.151 1.103TS13/TS13* 41.70 -126.3 31.5 -85.0 112.4 112.4 114.7 1.448 2.683 1.336P13/P13* -7.24 178.5 2.1 -69.1 104.6 107.8 115.2 4.034 4.170 1.100a φ1 = φ(C 1 − C 2 − C 3 − C 4 ), φ 2 = φ(C 3 − C 4 − C 5 − C 6 ), φ 3 = φ(C 5 − C 6 − C 7 − C 1 0), α 1 = α(C 6 − C 7 − H 17 ), α 2 = α(H 10 − C 1 − H 11 ), α 3 = α(C 8 − C 7 − C 9 )92 3.1. Superficie <strong>de</strong> energía potencial
3. Resultados y discusión 933.2. Cálculos <strong>de</strong> dinámica directaPara llevar a cabo simu<strong>la</strong>ciones <strong>de</strong> <strong>la</strong>s reacciones químicas se necesita informaciónacerca <strong>de</strong> <strong>la</strong> superficie <strong>de</strong> energía potencial <strong><strong>de</strong>l</strong> sistema que es objeto <strong>de</strong>estudio. La superficie <strong>de</strong> energía potencial se pue<strong>de</strong> obtener como una funciónanalítica ajustada total o parcialmente a energías calcu<strong>la</strong>das mediante métodosteóricos <strong>de</strong> <strong>la</strong> estructura electrónica o a propieda<strong>de</strong>s empíricas ajustables.Sin embargo, con <strong>la</strong> actual potencia <strong>de</strong> cálculo es posible usar directamente <strong>la</strong>steorías <strong>de</strong> <strong>la</strong> estructura electrónica en <strong>la</strong>s simu<strong>la</strong>ciones, prescindiendo <strong>de</strong> <strong>la</strong>senergías <strong>de</strong> potencial analíticas o empíricas, es <strong>de</strong>cir, que <strong>la</strong> energía potencialy sus <strong>de</strong>rivadas se calcu<strong>la</strong>n “al vuelo”. Este tipo particu<strong>la</strong>r <strong>de</strong> simu<strong>la</strong>ciones basadasen <strong>la</strong> aproximación <strong>de</strong> Born-Oppenheimer se <strong>de</strong>nominan, simu<strong>la</strong>ciones<strong>de</strong> dinámica directa Born-Oppenheimer.En <strong>la</strong> sección 2.4 <strong>de</strong> esta memoria se presenta <strong>la</strong> teoría <strong>variacional</strong> <strong><strong>de</strong>l</strong><strong>estado</strong> <strong>de</strong> <strong>transición</strong>. Esta teoría requiere cierta información <strong>de</strong> <strong>la</strong> superficie<strong>de</strong> energía potencial, que pue<strong>de</strong> ser calcu<strong>la</strong>da al vuelo. Este es el método quehemos utilizado para obtener <strong>la</strong>s constantes <strong>de</strong> velocidad <strong><strong>de</strong>l</strong> <strong>de</strong>sp<strong>la</strong>zamientosigmatrópico [1,7] <strong>de</strong> hidrógeno en el 7-metilocta-1,3(Z),5(Z)-trieno .3.2.1. Detalles computacionalesLa superficie <strong>de</strong> energía potencial se ha obtenido con el programaMOPAC [65] , mientras que para el cálculo <strong>de</strong> <strong>la</strong>s constantes <strong>de</strong> velocidad hemosutilizado el programa <strong>de</strong>sarrol<strong>la</strong>do en el grupo <strong><strong>de</strong>l</strong> Profesor D. G. Truh<strong>la</strong>r,l<strong>la</strong>mado POLYRATE [66] .La teoría <strong>variacional</strong> <strong><strong>de</strong>l</strong> <strong>estado</strong> <strong>de</strong> <strong>transición</strong> requiere, al menos, el cálculo<strong><strong>de</strong>l</strong> camino <strong>de</strong> mínima energía para obtener <strong>la</strong>s constantes velocidad canónicas<strong>variacional</strong>es k CVT/tun . En el caso <strong>de</strong> que, lo que se <strong>de</strong>see es evaluar elefecto túnel <strong>de</strong> gran curvatura, es necesario calcu<strong>la</strong>r una región más amplia<strong>de</strong> <strong>la</strong> superficie <strong>de</strong> energía potencial. En cálculos <strong>de</strong> dinámica directa el programaPOLYRATE necesita valores <strong>de</strong> <strong>la</strong> superficie <strong>de</strong> energía potencial y sus<strong>de</strong>rivadas en estas regiones, concretamente en <strong>de</strong>terminados puntos <strong>de</strong> <strong>la</strong>s mismas.El programa <strong>de</strong> <strong>la</strong> estructura electrónica, MOPAC, <strong>de</strong>be calcu<strong>la</strong>r para<strong>la</strong> geometría nuclear correspondiente a esos puntos, <strong>la</strong> energía potencial y en<strong>de</strong>terminados casos el gradiente y <strong>la</strong> matriz <strong>de</strong> constantes <strong>de</strong> fuerza. Un tercerprograma <strong>de</strong>be servir <strong>de</strong> en<strong>la</strong>ce entre ambos programas.A pesar <strong>de</strong> que el grupo <strong>de</strong> D. G. Truh<strong>la</strong>r <strong>de</strong>sarrolló una interfaz <strong>de</strong> unión[65] J. J. P. Stewart, “MOPAC: A semiempirical molecu<strong>la</strong>r orbital program,” Journal ofComputer-Ai<strong>de</strong>d Molecu<strong>la</strong>r Design, vol. 4, no. 1, pp. 1–103, 1990.[66] J. C. Corchado, Y.-Y. Chuang, P. L. Fast, W.-P. Hu, Y.-P. Liu, G. C. Lynch, K. A.Nguyen, C. F. Jackels, A. Fernan<strong>de</strong>z Ramos, B. A. Ellingson, B. J. Lynch, V. S. Melissas,J. Villà, I. Rossi, E. L. Coitiño, J. Pu, T. V. Albu, R. Steckler, B. C. Garrett, A. D.Isaacson, and D. G. Truh<strong>la</strong>r , “POLYRATE 9.7: Computer Program for the Calcu<strong>la</strong>tionof Chemical Reaction Rates for Polyatomics, version 9.3.1.” Department of Chemistryand Supercomputing Institute University of Minnesota, Minneapolis, MN, June 2005.
- Page 5:
VANTONIO FEZNÁNDEZ RAMOS Y SAULO A
- Page 9 and 10:
Índice general1. Introducción y o
- Page 12 and 13:
XIIÍNDICE DE FIGURAS3.4. Confórme
- Page 15 and 16:
Capítulo 1Introducción y objetivo
- Page 17 and 18:
1. Introducción y objetivos 3su ef
- Page 19 and 20:
1. Introducción y objetivos 5Tabla
- Page 21 and 22:
1. Introducción y objetivos 7algun
- Page 23 and 24:
1. Introducción y objetivos 91.3.
- Page 25 and 26:
1. Introducción y objetivos 11Tabl
- Page 27 and 28:
1. Introducción y objetivos 13la p
- Page 29 and 30:
1. Introducción y objetivos 15C 8
- Page 31 and 32:
1. Introducción y objetivos 17Tabl
- Page 33 and 34:
1. Introducción y objetivos 19C 8
- Page 35 and 36:
1. Introducción y objetivos 21Figu
- Page 37 and 38:
1. Introducción y objetivos 23medi
- Page 39 and 40:
Capítulo 2Métodos teóricosPara a
- Page 41 and 42:
2. Métodos teóricos 27Sin embargo
- Page 43 and 44:
2. Métodos teóricos 29mínimo. Mi
- Page 45 and 46:
2. Métodos teóricos 31podrán omi
- Page 47 and 48:
2. Métodos teóricos 33Para poder
- Page 49 and 50:
2. Métodos teóricos 35esférica d
- Page 51 and 52:
2. Métodos teóricos 37recibe el n
- Page 53 and 54:
2. Métodos teóricos 39Las integra
- Page 55 and 56: 2. Métodos teóricos 41electrón l
- Page 57 and 58: 2. Métodos teóricos 43son de áto
- Page 59 and 60: 2. Métodos teóricos 45e INDO. Sin
- Page 61 and 62: 2. Métodos teóricos 47centros y c
- Page 63 and 64: 2. Métodos teóricos 492.3.3. Mét
- Page 65 and 66: 2. Métodos teóricos 51cierto tama
- Page 67 and 68: 2. Métodos teóricos 53Teorema var
- Page 69 and 70: 2. Métodos teóricos 55el funciona
- Page 71 and 72: 2. Métodos teóricos 57• Por úl
- Page 73 and 74: 2. Métodos teóricos 59este nivel
- Page 75 and 76: 2. Métodos teóricos 61transición
- Page 77 and 78: 2. Métodos teóricos 63generalizad
- Page 79 and 80: 2. Métodos teóricos 65expresión
- Page 81 and 82: 2. Métodos teóricos 67(o 2) rotac
- Page 83 and 84: 2. Métodos teóricos 69partición
- Page 85 and 86: 2. Métodos teóricos 71en un incre
- Page 87: 2. Métodos teóricos 73La segunda
- Page 90 and 91: 76 3.1. Superficie de energía pote
- Page 92 and 93: 78 3.1. Superficie de energía pote
- Page 94 and 95: 80 3.1. Superficie de energía pote
- Page 96 and 97: 82 3.1. Superficie de energía pote
- Page 98 and 99: Tabla 3.2: Energías en kcal/mol y
- Page 100 and 101: 86 3.1. Superficie de energía pote
- Page 102 and 103: 88 3.1. Superficie de energía pote
- Page 104 and 105: Figura 3.9: Estructuras conformacio
- Page 108 and 109: 94 3.2. Cálculos de dinámica dire
- Page 110 and 111: 96 3.2. Cálculos de dinámica dire
- Page 112 and 113: 98 3.2. Cálculos de dinámica dire
- Page 114 and 115: 100 3.2. Cálculos de dinámica dir
- Page 117 and 118: Apéndice AMétodo simplex descende
- Page 119: A. Método simplex descendente mult
- Page 122 and 123: 108el que se encuentra en favor de
- Page 124 and 125: 110 C.1. MORATE2end programCsubrout
- Page 126 and 127: 112 C.1. MORATE2if (i.eq.icurt.or.i
- Page 128 and 129: 114 C.2. CONFORATE& ,wts(nconfs,nfr
- Page 130 and 131: 116 C.2. CONFORATEC TSif (nat_ts.gt
- Page 132 and 133: 118 C.2. CONFORATECCCCif (iop6.eq.1
- Page 134 and 135: 120 C.2. CONFORATEakp_vtst(n,i)=ct(
- Page 136 and 137: 122 C.2. CONFORATECCwrite(6,30) tem
- Page 138 and 139: 124 C.2. CONFORATEC Compute traslat
- Page 140 and 141: 126 C.2. CONFORATEC Read input data
- Page 142 and 143: 128 C.2. CONFORATECCCCenddoelsewrit
- Page 144 and 145: 130 C.2. CONFORATEstop1011 write(6,
- Page 146 and 147: 132 C.2. CONFORATElin=0C write(6,*)
- Page 148 and 149: 134 C.2. CONFORATE& ,numsta,numstb,
- Page 150 and 151: 136 C.2. CONFORATEwrite(6,33) (dmir
- Page 152 and 153: 138 C.2. CONFORATEendC Compute the
- Page 154 and 155: 140 C.2. CONFORATEC Full F matrixdo
- Page 156 and 157:
142 C.2. CONFORATE* NAGOYA UNIVERSI
- Page 158 and 159:
144 C.2. CONFORATE* * * * * * * * *
- Page 160 and 161:
146 C.2. CONFORATEDO 400 J=1,N400 V
- Page 162 and 163:
148 C.2. CONFORATE298. 63.0843 # Co
- Page 165 and 166:
Bibliografía[1] C. Funk, “On the
- Page 167 and 168:
BIBLIOGRAFÍA 153[24] A. Verloop, A
- Page 169 and 170:
BIBLIOGRAFÍA 155[47] D. Hartree,
- Page 171:
BIBLIOGRAFÍA 157G. Liu, A. Liashen