a variacional del estado de transición a la - Páxinas persoais - USC ...

88 3.1. Superficie de energía potencial1. Con los parámetros por defecto disponibles en el programa MOPAC [65]se obtienen los valores de la variables que hemos elegido.2. Se calcula la función de los errores cuadrados∑(Y i,final − Y i ) 2 W i (3.1)idonde Y i,final son los valores óptimos de las variables, Y i son los valoresque toman esas mismas variables después de cada bucle del ciclo y W ison los factores de peso que tiene cada variable.3. Se minimiza la función mediante el método simplex, que busca mínimoslocales de la función 3.1 y el método del temple simulado que intentalocalizar mínimos globales. En cada paso de la optimización se deben generarunos nuevos valores de los parámetros AM1 y reevaluar la funciónde los errores cuadráticos medios y recalcular las variables definidas enla ecuación anterior.4. El proceso se repite hasta que los errores están por debajo de un valorumbral elegido por el usuario.Con este procedimiento podemos elegir los parámetros más adecuados quemejoren las características de nuestro modelos teórico. Así hemos sido capacesde generar los parámetros para nuestra superficie a nivel AM1-SRP en base alas propiedades físicas y químicas que nos interesan.De manera análoga a la obtención de la superficie de energía potencialMPWB1K y partiendo de los confórmeros optimizados a este nivel hemos obtenidolos nuevos confórmeros a nivel AM1-SRP. Geométricamente y energéticamente(Tabla 3.3) son muy similares a los obtenidos mediante el métodoMPWB1K con la única salvedad de que una de las estructuras, la R5 desapareceen esta superficie. De igual modo hemos calculado el estado de transiciónpara la migración de hidrógeno. En este caso la barrera energética es de 22.62kcal/mol, que es algo mayor que a nivel DFT. También existen diferencias enla frecuencia imaginaria del estado de transición; a nivel AM1-SRP la frecuenciaes mayor (1800i cm −1 frente a los 1400i cm −1 del método DFT). Estose traduce en un perfil de la energía de reacción más afilado, lo que conllevaimportantes diferencias en la contribución de túnel a estos dos niveles.[65] J. J. P. Stewart, “MOPAC: A semiempirical molecular orbital program,” Journal ofComputer-Aided Molecular Design, vol. 4, no. 1, pp. 1–103, 1990.