430 likes | 970 Views
QUIMICA FISICA AVANZADA. QUIMICA CUANTICA – METODOLOGIA Supongamos una molécula formada por 3 átomos, cada uno aportando un electrón descrito por un orbital hidrogenoide: át.1 + át2 + át.3 Mol La teoría LCAO-MO propone:
E N D
QUIMICA FISICA AVANZADA QUIMICA CUANTICA – METODOLOGIA Supongamos una molécula formada por 3 átomos, cada uno aportando un electrón descrito por un orbital hidrogenoide: át.1 + át2 + át.3 Mol La teoría LCAO-MO propone: O sea, la función molecular de prueba i, llamada Orbital Molecular, se puede expresar como una Combinación Lineal de funciones átomicas, orbitales atómicos, j conocidas como “funciones base”: Siguiendo la metodología de la Teoría de la Variación Lineal para nuestro caso de 3 electrones, cada uno sobre un átomo distinto: 1) Función de prueba
QUIMICA FISICA AVANZADA QUIMICA CUANTICA – METODOLOGIA 2) La energía: Donde: Las H son las integrales de energía promedio y la S integrales de solapamiento. Sustituyendo: 3) Condición de mínimo: Derivando:
QUIMICA FISICA AVANZADA QUIMICA CUANTICA – METODOLOGIA Derivando: Se trata de un sistema de ecuaciones homogéneas con una solución trivial (ci,j=0) y otra cuando el determinante de los coeficientes, “determinante secular”, es igual a cero: Resolviendo el determinante secular se obtiene una ecuación en E3 con tres raices (E1, E2, E3). Para ello es necesario conocer el valor de las integrales Hij y Sij. Si la base de orbitales atómicos está normalizada entonces Sij = ij , lo cual simplifica mucho la resolución del determinante. Luego, con la relación entre los coeficientes y la condición de normalización se puede obtener el valor de tres juegos de coeficientes, uno para cada valor de Ei.
QUIMICA FISICA AVANZADA QUIMICA CUANTICA – METODOLOGIA Con los valores de los coeficientes individuales se obtienen tres orbitales moleculares i del tipo: Procedimiento matricial: El sistema de ecuaciones homogéneas (1) se puede reescribir en forma matricial: A su vez:
QUIMICA FISICA AVANZADA QUIMICA CUANTICA – METODOLOGIA
QUIMICA FISICA AVANZADA QUIMICA CUANTICA – METODOLOGIA Que también se puede escribir: H c = E S c (2) Donde H es la matriz de las integrales del Hamiltoniano, c es el vector de los coeficientes del orbital molecular i , S es la matriz solapamiento y E es la energía que corresponde al orbital i. Con respecto a la matriz S se debe tener presente que los orbitales atómicos están centrados sobre distintos átomos. Si la base de orbitales atómicos está ortonormalizada entonces Sij = ij , o sea, la matriz S es la matriz unidad U (Todos los elementos diagonales iguales a 1 y los demás a 0).. Entonces: H c = E c (3) Si la base no estuviera ortonormalizada se puede ortogonalizar S diagonalizándola, o sea, mediante una transformación de similitud que lleve a un sistema donde S sea diagonal y luego normalizarla por el procedimiento habitual. El sistema (3) representa el problema de Autovalores y Autovectores que se resuelve diagonalizando H.
QUIMICA FISICA AVANZADA QUIMICA CUANTICA – METODOLOGIA Donde: C es la matriz de los coeficientes de los autovectores (orbitales moleculares i). Cada columna es un autovector. E es la matriz diagonal de los autovalores. Cada elemento diagonal Eii es un autovalor, o sea, la energía del correspondiente orbital molecular i. Por ejemplo, en MatLab esta operación se codifica: (C,E) = EIG(H) que significa: diagonalice H devolviendo la matriz ortogonal C y la matriz diagonal E. Finalmente, con las columnas de la matriz C se construyen los OM:
QUIMICA FISICA AVANZADA QUIMICA CUANTICA – METODOLOGIA METODO DE HUCKEL SIMPLE Se trata de un método semiempírico donde se utiliza la metodología teórica junto con parámetros empíricos. No se toma en cuenta la geometría molecular sino la conectividad. 1. Intervienen únicamente orbitales pz centrados en distintos átomos formando la base de orbitales atómicos. 2. Los OM son únicamente del tipo π. 3. La base de orbitales atómicos está ortonormalizada, <ij> = Sij= ij. 4. Hii= (integrales de coulomb); 5. Hij (integrales de resonancia). Hij= para átomos vecinos y Hij=0 para átomos no vecinos. Dado que la base es ortogonal, se cumple el sistema 3: H c = E c Ejemplo: Propeno, H3C1-HC2-C3H2 . Consideremos =0 y =1. Diagonalizar H y formular los OM con sus respectivas energías.
QUIMICA FISICA AVANZADA QUIMICA CUANTICA – METODOLOGIA METODO DE HOFFMAN (EXTEDED HUCKEL) 1983 Hoffman y Fukuy Se tiene en cuenta la geometría molecular. 1. La base de orbitales atómicos incluye todos los orbitales de valencia (s,p) y no está ortonormalizada: <ij> = Sijij. No intervienen los electrones del “core”. 2. Los OM son de tipo y π. 3. Las integrales Hij y Hij se expresan como funciones de las integrales de solapamiento: Hij = f(Sij) Debido a que la base no es ortogonal se cumple el sistema (2): H c = E S c Por ello se procede: a. Se diagonaliza S para obtener una nueva base ortogonalizada ’. Se normaliza. b. Se transforma el sistema (2) al (3) aplicando la transformación a la nueva base ’ (en la cual S’ es diagonal H’ c’ = E c’ (3’) c. Se diagonaliza H’ como en Huckel simple.
QUIMICA FISICA AVANZADA QUIMICA CUANTICA – METODOLOGIA UNIDADES ATOMICAS Permiten simplificar las ecuaciones mediante la eliminación de las constantes fundamentales. La unidad atómica de longitud es el radio de Bohr: La energía se expresa en Hartrees y la masa en unidades de masa electrónica: Para el átomo de hidrógeno:
QUIMICA FISICA AVANZADA QUIMICA CUANTICA – METODOLOGIA METODO DE HARTREE (Átomo de He) Aparece el término de correlación electrónica (repulsión). No tiene solución exacta. Vimos que la aproximación de partículas independientes, ignorando este término, no conducía a buenos resultados. Consideremos ahora que dicho término se pueda dividir en el aporte individual de cada electrón, o sea, en dos componentes monoelectrónicas. Entonces la solución de la Ecuación de Schrodinger estará dada por un conjunto de dos ecuaciones, monoelectrónicas, una para cada electrón: Sistema de ecuaciones diferenciales donde la función de onda total y la energía electrónica total son:
QUIMICA FISICA AVANZADA QUIMICA CUANTICA – METODOLOGIA Esta sugerencia fue hecha por Douglas R. Hartree. Se trata de funciones de onda monoelectrónicas: Además propuso una forma de resolver aproximadamente el término Vi(ee): si se coloca al electrón 1 en una posición del espacio, la energía de repulsión de este electrón en relación al campo promedio producido por el otro electrón 2 será Esta expresión se lleva al sistema de ecuaciones anterior y se obtiene un sistema de ecuaciones integrodiferenciales: Solo se puede resolver mediante un procedimiento iterativo hasta lograr autoconsistencia. Por esta razón a este procedimiento se lo llama Campo Auto Consistente, Self Consistent Field, SCF.
QUIMICA FISICA AVANZADA QUIMICA CUANTICA – METODOLOGIA El procedimiento (algoritmo) es el siguiente: 1) COMIENZO 2) Se proponen las funciones de prueba º1(1), º2(2). 3) Se calculan los potenciales de repulsión interelectrónica que experimenta cada electrón (integrales) y el resultado se lleva al sistema de ecuaciones diferenciales. 4) Se resuelve el sistema de ecuaciones diferenciales y como resultado se obtienen las funciones 1(1), 2(2). 5) Se comparan las funciones º1(1), º2(2) con las obtenidas 1(1), 2(2). SI 1(1), 2(2) º1(1), º2(2) ENTONCES se reasignan º1(1) 1(1) º2(2) 2(2) y se continúa en el paso 3). SINO se continúa en el paso 6). 6) Se imprimen los resultados (1(1), 2(2)) y (ε1, ε2). 7) FIN El método de Hartree se puede extender a átomos polielectrónicos y también, una vez aplicada la aproximación de Born-Oppenheimer (núcleos fijos), a moléculas.
QUIMICA CUANTICA METODOLOGIA
QUIMICA FISICA AVANZADA QUIMICA CUANTICA – METODOLOGIA PRINCIPIO DE EXCLUSIÓN DE PAULI: 1925 Wolfwang Pauli. Considera la indistinguibilidad de los electrones. “Las funciones de onda de los electrones deben ser antisimétricas respecto al intercambio de electrones”. Los electrones pueden tener igual función espacial (orbital) pero distinta función de espín: espín-orbitales. 1929 J. C. Slater. Señaló que para representar las funciones de onda antisimétricas se pueden usar determinantes. Los conceptos anteriores, sumados a la aproximación de Born-Oppenheimer, definen el método de Hartree-Fock para la resolución de la ecuación de Scrhodinger molecular. Consiste en resolver un conjunto de N ecuaciones integro-diferenciales acopladas (N es el número de OM). Para ello se tiene en cuenta el principio variacional. “Cualquier función de onda que uno proponga para el sistema tendrá una energía mayor, o a lo sumo, igual que la energía verdadera del estado fundamental de dicho sistema”.
QUIMICA FISICA AVANZADA QUIMICA CUANTICA – METODOLOGIA METODO DE HARTREE-FOCK La ecuaciones de Hartree-Fock resultan de minimizar la energía respecto a los OM. Donde Fi es el operador de Fock: Los primeros dos términos corresponden al Hamiltoniano de Core que representa la energía de cada electrón en el campo promedio desnudo de los núcleos. Jij es el operador de Coulomb y Kij el operador de Intercambio. Dado que se trata de operadores integro-diferenciales los mismos dependen de sus propias soluciones. Por ello las ecuaciones de Hartree-Fock se deben resolver en forma autoconsistente (SCF).
QUIMICA FISICA AVANZADA QUIMICA CUANTICA – METODOLOGIA Cuando se utiliza la teoría LCAO-OM los OM están desarrollados como combinación lineal de orbitales atómicos: La minimización de la energía total se realiza entonces con respecto a los coeficientes de la expansión. Este procedimiento conduce a las ecuaciones deducidas por Roothaan y Hall que se pueden expresar matricialmente: FC = SCε Donde F es la matriz formada por el operador de Fock y S es la matriz de las integrales de solapamiento: La matriz de Fock incluye otras matrices como: la matriz Hcore, la matriz de Coulomb, la matriz de Intercambio y la matriz Densidad. Esta última es importante porque sobre ella se realiza el proceso iterativo para llegar a autoconsistencia, comparando sus elementos entre dos iteraciones sucesivas.
QUIMICA FISICA AVANZADA QUIMICA CUANTICA – METODOLOGIA El algoritmo de Hartree-Fock, en forma resumida, es: 1. Comienzo 2. Se calculan todas las integrales necesarias. 3. Se construye una matriz densidad aproximada P0. 4. Se construye la matriz del operador de Fock F. 5. Se resuelve FC = SCε y se determinan C. 6. Se construye una nueva matriz densidad P1. 7. Si P1- P0 test, ENTONCES Fin del ciclo de autoconsistencia. SINO P0 P1 y se continúa en el paso 4. Para resolver el sistema FC = SCε es necesario ortonormalizar la base (SU) y luego diagonalizar F en la nueva base. Una mejor aproximación a nuestro problema molecular dependerá de la base que utilicemos para expandir o desarrollar los orbitales moleculares. Si la base que se utiliza tuviera infinitos términos se dice que se ha alcanzado el Límite de Hartree-Fock. Las funciones gaussianas facilitan el cálculo de las integrales analíticamente pero aumentan la dimensión de las matrices y con ello la cantidad de cálculo. La situación resulta ventajosa para resolver con computadoras.
QUIMICA CUANTICA METODOLOGIA
QUIMICA FISICA AVANZADA QUIMICA CUANTICA – METODOLOGIA CAPA ABIERTA Hasta aquí hemos considerado Hartree-Fock para electrones apareados. Para el caso de electrones desapareadotes necesario definir dos conjuntos distintos de coeficientes para los electrones alfa y beta: METODOS POST HARTREE-FOCK Al incluir la antisimetría de la función de onda el método de Hartree-Fock está considerando una parte de la correlación electrónica que se llama de Intercambio. Sin embargo no tiene en cuenta la correlación dinámica debida a la cual los electrones tratan de permanecer alejados por repulsión. Los métodos que tratan de resolver este problema se llaman métodos Post Hartree-Fock
QUIMICA FISICA AVANZADA QUIMICA CUANTICA – METODOLOGIA INTERACCION DE CONFIGURACIONES El método de HF produce un conjunto {i} de 2K espín orbitales (K es el número de funciones espaciales i y 2K el número de espín orbitales: i y iβ). Para el estado fundamental se formula el determinante (ver Fig. a): El método CI incluye además del determinante de Hartree-Fock para la función de onda otros determinantes que se obtienen al reemplazar un orbital ocupado por otro virtual. Por ejemplo para una simple sustitución (ai) el determinante quedaría (ver Fig. b): Equivale a una excitación. Para una doble sustitución (ar) y (bs) (ver Fig. c) : Equivale a doble excitación.
2K 2K 2K s r s r s r N+2 N+1 N+2 N+1 N+2 N+1 N N N b a b a b a 2 1 2 1 2 1 QUIMICA FISICA AVANZADA QUIMICA CUANTICA – METODOLOGIA INTERACCION DE CONFIGURACIONES
QUIMICA FISICA AVANZADA QUIMICA CUANTICA – METODOLOGIA INTERACCION DE CONFIGURACIONES FULL CI En esta aproximación la función de onda es una combinación lineal de determinantes de Hartree-Fock: Donde 0 indica nivel Hartree-Fock y s corre sobre todas las posibles sustituciones. Full CI proporciona la mejor corrección a la correlación electrónica. Para una base infinita se dice que se ha alcanzado la solución exacta de la Ecuación de Schrodinger independiente del tiempo. Sin mbargo Full CI resulta impracticable. Por ello se utiliza la expansión Ci truncada a algún nivel de sustitución: CIS, CID, CISD, CISDT, etc.
QUIMICA FISICA AVANZADA QUIMICA CUANTICA – METODOLOGIA TEORIA DE PERTURBACIONES Donde Ĥ0 tiene solución exacta y ĥ es una corrección que se incluye como una pequeña perturbación. Por lo tanto se puede desarrollar como una serie de potencias. De esta manera la función de onda y la energía quedan: Estas expresiones se llevan a la Ecuación de Schrodinger: Se distribuyen los productos y se igualan aquellos que corresponden la misma potencia de :
QUIMICA FISICA AVANZADA QUIMICA CUANTICA – METODOLOGIA TEORIA DE MOLLER-PLESSET La teoría de Moller-Plesset define Ĥ0 como la suma de los operadores monoelectrónicos de Fock: Este hamiltoniano se aplica a las aproximaciones anteriores resultando los distintos órdenes de perturbación: 2, 3, 4, etc. conocidos como MP2, MP3, MP4, etc. Por tratarse de un método de perturbación en este caso no existe una cota inferior a la energía del sistema.
QUIMICA FISICA AVANZADA QUIMICA CUANTICA – METODOLOGIA METODOS SEMIEMPIRICOS Los distintos términos del operador de Fock producen un número muy grande de integrales bi-, tri- y tetracéntricas, y el cálculo de las mismas supone un costo computacional muy alto. El número de integrales es igual al número de orbitales atómicos elevados a la curta potencia). El propósito de los métodos semiempíricos es reducir el número de integrales sin afectar el resultado final. CNDO: Anulación completa del solapamiento diferencial (ZDO). Desprecia algunas integrales bieletrónicas y reemplaza algunas monoelectrónicas por parámetros ajustables empíricamente. Autor, J. A. Pople (1998 Kohn y Pople ) INDO: Anulación parcial del solapamiento diferencial. Incluye en el cálculo integrales bielectrónicas cuando son monocéntricas.
QUIMICA FISICA AVANZADA QUIMICA CUANTICA – METODOLOGIA METODOS SEMIEMPIRICOS Debido a los pobres resultados que se obtenían con estos métodos, Dewar y colaboradores propusieron parametrizar masivamente, reemplazando las integrales por parámetros que dependen principalmente de la naturaleza de cada átomo. Dewar: MINDO, MNDO y AM1 (Austin Method 1), este último orientado a sistemas con uniones hidrógeno. Stewart: PM3 (Parametric Method 3), similar al AM1 pero con un conjunto distinto de parámetros.
QUIMICA FISICA AVANZADA QUIMICA CUANTICA – METODOLOGIA TEORIA DEL FUNCIONAL DE LA DENSIDAD Thomas, Fermi y Dirac en la década del 1920 y Slater en los 50. La idea es usar funcionales de la densidad electrónica para atender el problema de la correlación. Se busca presentar a la energía total de un sistema dado como un funcional de la densidad electrónica () en lugar de la función de onda (). 1964 Hohenberg y Kohn demuestran un teorema que prueba la existencia de un único funcional que determina la energía y la densidad del estado fundamental en forma exacta. (1998 Kohn y Pople ) 1965 Kohn y Sham definieron un sistema de N electrones sometidos a un potencial externo vext y con una densidad: Obtuvieron como resultado N ecuaciones monoelectrónicas del tipo de HF.
QUIMICA FISICA AVANZADA QUIMICA CUANTICA – METODOLOGIA TEORIA DEL FUNCIONAL DE LA DENSIDAD De acuerdo con Khon y Sham la teoría DFT particiona la energía electrónica en varios téminos: E = EK() + EV() + EJ() + EXC() EK: energía cinética debida al movimiento de los electrones. EV: energía potencial de atracción núcleo electrón y repulsión entre pares de núcleos. Esta última interacción no depende de . EJ: término de repulsión electrón-electrón de tipo coulombica. EXC: término de intercambio y correlación que incluye el resto de la interacción electrón-electrón: Energía de intercambio debida a la antisimetría de la función de onda. Correlación dinámica debida al movimiento de los electrones. EK, EV y EJ corresponden a la energía clásica de la distribución de carga
QUIMICA FISICA AVANZADA QUIMICA CUANTICA – METODOLOGIA TEORIA DEL FUNCIONAL DE LA DENSIDAD EXC() = EX() + EC() EX(): Término de Intercambio (interacciones del mismo espín): Funcional de Intercambio. EC():Término de Correlación (interacciones del mezcla de espín): Funcional de Correlación. Estos funcionales se llaman “locales” (LDA: Local Density Aproximation) cuando solo dependen de y “corregidos por gradiente” (NLDA) cuando dependen además del gradiente de la densidad .
QUIMICA FISICA AVANZADA QUIMICA CUANTICA – METODOLOGIA Para un gas homogéneo de electrones (no muy apropiado para describir un sistema molecular) el funcional local de intercambio es: En 1988 Becke formuló el siguiente funcional basado en la LDA corregido por introducción del : De la misma manera, en 1991 Perdew y Wang formularon un funcional de correlación corregido por donde aparecen los términos que corresponden a las densidades de espín y . BLYP: es un funcional muy conocido que complementa el funcional de intercambio de Becke corregido por gradiente con el funcional de correlación de Lee, Yang y Parr, también corregido por gradiente.
QUIMICA FISICA AVANZADA QUIMICA CUANTICA – METODOLOGIA En 1993 Becke propuso mezclar el intercambio que proviene de la DFT con el que surge de HF, calculando este último con las ecuaciones de Kohn y Sham. Este tipo de funcional se llama híbrido: Donde los términos NLDA son términos corregidos por gradiente y corresponden al funcional de Becke. Los parámetros c0, cX y cC se determinan por ajuste de datos. Este funcional es muy conocido y se llama B3LYP. Los métodos DFT tienen un costo computacional similar al de HF pero con la ventaja de que incluyen la correlación electrónica.
QUIMICA FISICA AVANZADA QUIMICA CUANTICA – METODOLOGIA MECÁNICA MOLECULAR (MM) Se basa en formular una función de energía potencial (FEP) analítica que describe la interacción entre los átomos. El objetivo es optimizar la FEP en función de longitudes, ángulos y diedros de enlace. Comprende: 1. La expresión de la FEP. 2. Los Tipos de átomos. 3. Los parámetros estándar.
QUIMICA FISICA AVANZADA QUIMICA CUANTICA – METODOLOGIA MECÁNICA MOLECULAR (MM) La FEP: Función de Energía Potencial Incluye los siguientes términos: a. Átomos enlazados: b. Átomos segundos vecinos: c. Átomos terceros vecinos:
QUIMICA FISICA AVANZADA QUIMICA CUANTICA – METODOLOGIA MECÁNICA MOLECULAR (MM) d. Átomos no enlazados: e. Enlaces hidrógeno:
QUIMICA FISICA AVANZADA QUIMICA CUANTICA – METODOLOGIA MECÁNICA MOLECULAR (MM) Tipos de átomos Los cálculos de la MM se efectúan sobre tipos de átomos y no sobre átomos. El tipo de átomo tiene en cuenta la naturaleza del átomo, su hibridación, carga formal y el entorno inmediato. Así un átomo de C sp3 en metilo, por ejemplo, tendrá un tipo de átomo distinto de un C aromático en benceno. Parámetros Son componentes de la expresión de la FEP que se determinan por ajuste de datos experimentales dentro de una familia de compuestos. Por eso cuando se aplica la MM se debe tener en cuenta sobre cual familia de compuestos se ha ajustado los parámetros.
QUIMICA FISICA AVANZADA QUIMICA CUANTICA – METODOLOGIA DINÁMICA MOLECULAR (DM) (2013Karplus, Levitt y Warshel) Consiste en resolver numéricamente ecuaciones de movimiento clásicas (DM Clásica) con el objeto de construir trayectorias y de calcular propiedades como promedio de dichas trayectorias. El sistema es un conjunto de partículas que interactúan entre sí mediante un potencial V, aditivo de a pares. El punto de partida es el conjunto de partículas en su estado inicial ro,vo, posición y velocidad de las partículas. De acuerdo con la mecánica newtoniana: Donde Fij es la fuerza ejercida sobre la partícula “i” por la partícula “j”. Dicha fuerza se calcula a partir del potencial:
QUIMICA FISICA AVANZADA QUIMICA CUANTICA – METODOLOGIA DINÁMICA MOLECULAR (DM) De esta forma las ecuaciones de Newton se resuelven en forma numérica para un intervalo de tiempo dado [to,tn] obteniéndose la trayectoria de las partículas {r(t),v(t)}. Una propiedad A se calcula: Etapas de la DM: a. Inicialización: se asignan las coordenadas y velocidades iniciales de las partículas. b. Equilibración: cantidades como la energía cinética y la energía potencial, durante un determinado número de pasos, son controladas para que alcancen valores constantes. Al comienzo de la simulación se observan fuertes oscilaciones, que indican la falta de equilibrio del sistema. c. Adquisición de datos: esta etapa comienza cuando se alcanza el equilibrio, trayectoria estable. Se calculan las propiedades de interés a lo largo del resto de la trayectoria.
QUIMICA FISICA AVANZADA QUIMICA CUANTICA – METODOLOGIA DINÁMICA MOLECULAR (DM) En DM se deben especificar las condiciones en que se llevará acabo la simulación. Pueden ser a NVE (Número de partículas, volumen y energía constantes) o NVT (Número de partículas, volumen y temperatura constantes). Para el primero (DM-NVE) se establece una condición adicional, el momento lineal total debe ser nulo. Un ejemplo simplificado del algoritmo correspondiente podría ser: 1. Se especifican posiciones y velocidades iniciales para cada partícula ro,vo, 2. Se calculan las fuerzas {Fo}. 3. Se realiza una iteración durante el número de pasos preestablecidos y se calculan ri,vi y {Fi} en cada paso. 4. Se verifica que la energía total es igual a la deseada. En caso contrario se aplica vinueva = factor x vivieja el factor tiene en cuenta la desviación de la energía.