TEORÍA DE FUNCIONALES DE LA DENSIDAD

1. TEORÍA DE FUNCIONALES DE LA DENSIDAD Los métodos de Funcionales de la Densidad (DF: Density Functional) proveen una aproximación alternativa al tratamiento de la correlación electrónica en sistemas con varios electrones. La base de la Teoría de Funcionales de la Densidad (DFT: Density Functional Theory) es la prueba por Hohenber y Kohn que la energía electrónica del estado fundamental, la función de onda y todas las demás propiedades electrónicas, están determinadas completamente por la densidad electrónica r. Existe una correspondencia unívoca entre la densidad de electrones del sistema y la energía. DENSIDAD ELECTRÓNICA Representa el número de electrones por unidad de volumen en un punto del espacio. ni: número de ocupación, número de electrones en el orbital

2. Orbital molecular Espín-orbital molecular ni= 0,1 ó 2 ni= 0 ó 1 Orbital molecular espacial MOLÉCULA DE CAPA CERRADA Podemos escribir la densidad de electrones en términos de las funciones de base.

3. Definimos MATRIZ DENSIDAD Existe una correspondencia entre la densidad de electrones del sistema y la energía, esto supone ciertas ventajas prácticas. Una función de onda para un sistema de n electrones contiene 3n coordenadas (3 para cada electrón ó 4 si se incluye el espín). La densidad de electrones sólo depende de 3 coordenadas independientemente del número de electrones. Mientras la complejidad de una función de onda se incrementa con el número de electrones, la densidad electrónica tiene el mismo número de variables, independientemente del tamaño del sistema.

4. El “único” inconveniente es que aunque se ha probado que cada densidad diferente genera un estado de energía diferente, el funcional que conecta estas dos cantidades es desconocido. Función: f(x) Instrucciones o reglas para producir un número a partir de un conjunto de variables x. Es una regla que asocia un número a cada valor de la variable x para el que la función f está definida. Funcional: F[f] Instrucciones o reglas para producir un número a partir de una función f, la cual depende a su vez de variables. Es una regla que asocia un número a cada función f. Así, un funcional puede considerarse como una función de una función.

5. TEOREMA DE HOHENBERG-KOHN (HK) En 1964, Pierre Hohenberg y Walter Kohn probaron que la energía molecular del estado fundamental, la función de onda y todas las demás propiedades electrónicas, están determinadas unívocamente por la densidad de probabilidad electrónica del estado fundamental, una función de solamente tres variables. Se dice que la energía electrónica es un funcional de la densidad electrónica. La función de onda electrónica del estado fundamental de una molécula de n electrones es una función propia del Hamiltoniano electrónico. Potencial externo que actúa sobre el electrón i

6. HK demostraron que la densidad de probabilidad electrónica del estado fundamental, r0(r), determina el potencial externo y determina el número de electrones. r0(r), determina el número de electrones: r0(r), determina el potencial externo: Para demostrarlo suponemos que es falso y que hay dos potenciales externos va y vb que dan lugar a la misma densidad electrónica del estado fundamental r0. Sean y0,a y y0,b, E0,a y E0,b las funciones de onda del estado fundamental normalizadas y las energías para estos Hamiltonianos. Si el estado fundamental es no-degenerado, solamente hay una función de onda normalizada, la función de onda del estado fundamental exacta y0, que da la energía del estado fundamental exacta, E0.

7. Usamos y0,b como función de prueba variacional con el Hamiltoniano Ha Los Hamiltonianos difieren en sus potenciales externos Haciendo el mismo razonamiento intercambiando a por b

8. Por hipótesis las dos funciones de onda diferentes dan la misma densidad electrónica, r0,a = r0,b Haciendo r0,a = r0,b y sumando las desigualdades absurdo La densidad de probabilidad electrónica del estado fundamental r0 determina el potencial externo y el número de electrones  r0 determina el Hamiltoniano electrónico molecular y, de esta forma, determina la función de onda del estado fundamental, la energía y otras propiedades. La energía electrónica del estado fundamental E0 es un funcional de la función r0(r)

9. Si y0 es la función de onda del estado fundamental Cada uno de estos promedios es un funcional r

10. se conoce, pero los funcionales y son desconocidos es independiente del potencial externo La ecuación no proporciona la vía práctica para calcular E0 ya que el funcional F[r0] es desconocido TEOREMA VARIACIONAL DE HK HK demostraron que para toda función densidad de prueba rpr(r) que satisface para todo r, es válida la siguiente desigualdad La verdadera densidad electrónica del estado fundamental minimiza el funcional energía

11. EL MÉTODO DE KOHN-SHAM (KS) KS consideraron un sistema ficticio (denotado mediante el subíndice s y llamado sistema no-interactuante) de n electrones no interactuantes que experimentan todos ellos la misma función de energía potencial vs(ri), donde vs(ri) es tal que hace que la densidad electrónica del sistema de referencia sea igual a la densidad electrónica exacta. Los electrones no interactúan entre sí en el sistema de referencia, así, el Hamiltoniano del sistema de referencia Hamiltoniano de un electrón de KS La función de onda del estado fundamental ys,0 del sistema de referencia es un determinante de Slater de los espín-orbitales de KS Espín-orbitales de KS

12. parte espacial función de espín energías orbitales de KS KS reescribieron la ecuación de HK de la siguiente forma. Definimos

13. FUNCIONALES DESCONOCIDOS Se define el funcional de la energía de correlación-intercambio Eci[r] como La clave para efectuar con precisión un cálculo KS-DFT es tener una buena aproximación para la energía de correlación-intercambio. Necesitamos obtener la densidad electrónica del estado fundamental. El sistema ficticio de electrones se define por tener la misma r que el estado fundamental de la molécula: rs = r0.

14. Evaluamos los términos de la ecuación Podemos obtener E0 a partir de la densidad electrónica si podemos obtener los orbitales de KS y si conocemos cuál es el funcional Eci[r].

15. El teorema variacional de HK nos dice que podemos obtener la energía del estado fundamental variando r de forma que minimice el funcional Ev[r]. De forma equivalente, en lugar de variar r podemos variar los orbitales de KS que determinan r mediante Los orbitales de KS que minimizan E satisfacen Función energía potencial del sistema de referencia ECUACIONES DE KOHN-SHAM

16. Potencial de correlación-intercambio. Incluye los efectos de las correlaciones de intercambio (antisimetría) y electrónica. REALIZACIÓN DE LOS CÁLCULOS DF DE KS Se comienza con una elección inicial para r A partir de la elección inicial de r se obtiene una estimación inicial de vci(r) dependiendo de la forma funcional aproximada elegida (aquí es donde se diferencian los métodos DFT). Esta vci(r) se usa en las ecuaciones de KS Que se resuelven para dar la estimación inicial de los orbitales de KS, que se desarrollan en términos de funciones de base.

17. Para dar lugar a una serie de ecuaciones que se parecen a las ecuaciones de Roothaan Se resuelven estas ecuaciones, se obtienen las energías orbitales y los coeficientes. Como los orbitales se expresan como combinacion lineal de las funciones de base y los coeficientes, entonces a partir de los coeficientes calculados, se obtienen orbitales mejorados. Los orbitales mejorados se usan para calcular una densidad electrónica mejorada, que permite obtener un potencial mejorado, que se usa para obtener orbitales mejorados y así sucesivamente. La diferencia entre los métodos DFT radica en la elección de la forma funcional de la energía de correlación-intercambio (o del potencial de correlación-intercambio)

18. LA APROXIMACIÓN DENSIDAD LOCAL (LDA: Local Density Approximation) Son los modelos más simples. Asumen que la densidad local puede ser tratada como un gas de electrones uniforme, o lo que es lo mismo, que la densidad es una función que varía lentamente. HK mostraron que si r es una función que varía de forma extremadamente lenta con la posición, Eci[r] está dada con precisión por eci es la energía de intercambio más la correlación en un gas de electrones homogéneo con densidad electrónica r. Jellium: sistema hipotético eléctricamente neutro de volumen infinito, consistente en un número infinito de electrones interactuantes moviéndose en un espacio a través del cual la carga positiva está distribuida de forma continua y uniforme; el número de electrones por unidad de volumen tiene un valos constante no-nulo, r. Los electrones del Jellium constituyen un gas de electrones uniforme.

19. Vosko Wilk Nusair (VWN) Si tenemos ei y ec, tenemos E

20. EL MÉTODO Xa (X: exchange) Método desarrolado por Slater en 1951. Aquí se omite la contribución de la correlación a la Eci (es sustancialmente menor que la contribución de intercambio) y la contribución de intercambio se toma como a es un parámetro ajustable (en general entre 2/3 y 1) LA APROXIMACIÓN DENSIDAD DE ESPÍN LOCAL (LSDA: Local Spin Dentisty Approximation) Para moléculas de capa abierta y geometrías cercanas a la disociación, la LSDA da mejores resultados que la LDA. Mientras que en LDA los electrones apareados, de espín opuesto, tienen el mismo orbital espacial de KS, la LSDA permite (análogo al método UHF) que tales electrones tengan diferentes orbitales espaciales de KS.

21. FUNCIONALES DEL GRADIENTE CORREGIDO Buscan corregir el LSDA para la variación de la densidad electrónica con la posición. Esto se hace incluyendo las derivadas de la densidad en el integrando. Se hace que le energía de correlación-intercambio dependa no solamente de la densidad electrónica, sino también de sus gradientes. GGA: Aproximación del gradiente generalizado (Generalized Gradient Approximation) Normalmente las partes de intercambio y correlación se modelan separadamente

22. Funcionales de intercambio corregidos por gradiente PW86 (Perdew y Wang 1986) PWx91 (Perdew y Wang 1991) B88 (Becke 1988) Funcionales de correlación corregidos por gradiente LYP (Lee, Yang y Parr) P86 (Perdew 1986) PW91 (Perdew y Wang 1991) B86 (Becke 1986) Cualquier funcional de intercambio se puede combinar con cualquier funcional de correlación Ej: BLYP/6-31G* Cálculo DF con el funcional de intercambio B88 y el funcional de correlación LYP, con los orbitales KS expandidos en la base 6-31G*.

23. FUNCIONALES HÍBRIDOS Un funcional híbrido mezcla la fórmula para Ei de LDA con las fórmulas de Ei y Ec corregidas por gradiente. Por ejemplo el funcional híbrido B3LYP (donde el 3 indica un funcional de tres parámetros) viene definido por: Se usa una definición HF de la Ei a0= 0.20 ai= 0.72 ac= 0.81 Estos parámetros se eligieron de forma que se obtuviera un buen ajuste de las energías de disociación molecular experimentales. B3LYP funcional de tres parámetros con intercambio de B88 y correlación de LYP B3PW91 funcional de tres parámetros con intercambio de B88 y correlación de PW91 B3P86 funcional de tres parámetros con intercambio de B88 y correlación de P86

24. ALGUNAS CONSIDERACIONES • el trabajo de HK fue publicado en 1964-1965. Poco después, los físicos aplicaron DFT usando LSDA para estudiar los sólidos con considerable éxito. • en la mitad de los años ’80 aparecieron los funcionales de gradiente corregido, mejorando los resultados para el estudio de moléculas. DFT se aplicó a la química. • en 1993 se añadió la posibilidad de cálculos DF en los programas Gaussian. En la mitad de los años ’90, los cálculos DFT moleculares experimentaban un crecimiento explosivo. • DFT tiene la ventaja de permitir que se incluyan los efectos de la correlación en los cálculos, que consumen aproximadamente el mismo tiempo que los cálculos HF que no incluyen correlación.

25. las interacciones débiles debidas a la dispersión (interacciones de van der Waals) son pobremente descritas por los funcionales comunes. • ¿los métodos DFT son ab initio o semiempíricos? • Si ab initio se toma como ausencia de parámetros, los métodos LSDA son ab initio, pero los métodos de gradiente corregido pueden no serlo y los métodos híbridos no los son. Sin embargo el número de parámetros es significativamente menor que para los métodos empíricos. • El método semiempírico PM3 tiene 18 parámetros por cada átomo, mientras que el funcional de intercambio B88 sólo tiene una constante de ajuste, válida para toda la tabla periódica. • DFT también se basa en la Ecuación de Schrödinger. Los métodos ab initio usan el Hamiltoniano exacto y hacen aproximaciones sobre la función de onda, mientras que DFT hace aproximaciones sobre el Hamiltoniano.