teoría de funcionales de la densidad (dft)química cuántica los métodos dft se basan en la...
Post on 18-Jan-2015
10 Views
Preview:
TRANSCRIPT
Teoría de Funcionales de la Densidad (DFT) Química Cuántica
Los métodos DFT se basan en la premisa de que la energía de una molécula puede determinarse a partir de la densidad electrónica en lugar de usar la función de onda. Lo cual es computacionalmente mucho más eficiente.
Historia:
1927, Thomas y Fermi, consideraciones estadísticas para aproximar la distribución de los electrones en los átomos: Los electrones se encuentran distribuídos uniformemente en el espacio de fase de 6 dimensiones y existe un campo potencial efectivo que queda determinado por la carga nuclear y esta distribución electrónica
Calcularon la energía de un átomo, representando su energía cinética como un funcional de su densidad electrónica, y combinando esto con las expresiones clásicas de las interacciones núcleo-electrón y electrón-electrón (que también se pueden representar en términos de densidad electrónica
Funcional: una función de una función F gf g rf g f r
F g númeroUna función es una regla para ir de una variable a un númeroUn funcional es una regla para ir de una función a un número
Teoría de Funcionales de la Densidad (DFT) Química Cuántica
Aproximación de Thomas y Fermi (TF)
Si se divide el espacio en pequeños cubos (celdas) de lado ℓ, donde cada celda se considera independiente del resto, cada una tendrá un volúmen:
Y cada celda contiene una cantidad fija de electrones Δn, que puede ser diferente para cada celda, entonces la densidad electronica se define como:
Sumando las contribuciones de todas las celdas, la energía cinética total en unidades atómicas será:
Es el funcional de la densidad de TF para la energía cinética
Esta aproximación se corresponde a la LDA, donde la densidad electrónica se considera homogénea
La energía total del sistema atómico se escribe como:
3n
Vn
3V
2 35 3 33 2 871
10//
TF F FT C r dr C . donde
1 21 22
1 2
5 3 1
2/
FTF
rZ dr
r
r rdrE C r dr dr
r r
Y tiene que cumplir con n r dr
nTF eeeVTE V
P. Hohenberg, W. Kong, Phys. Rev. 1964, 864, B136. Demostraron que la E y otras propiedades de cualquier sistema químico están univocamente definidos por la densidad electrónicaNo dan idea de la forma matemática
Teoremas de Hohenberg y Kohn:
1. El potencial electrónico externo, (ri ), de un sistema con N e- está
determinado unívocamente por su (r) + una constante
Teoría de Funcionales de la Densidad (DFT) Química Cuántica
HK
ee
ee
F
K
ne
H
V
V
E
r r dr
r r d
V
F
T
r
T
Si incluimos explicitamente la dependenca de E con el potencial externo
JVee Término no clásico
Repulsión clásica
donde:
Que ellos no definen y que incluye la mayor parte de la E de correlación e intercambio
Independiente del potencial externoCuando se tiene su forma explícita Se puede aplicar a cualquier sistema
Teoría de Funcionales de la Densidad (DFT) Química Cuántica
2. Para toda densidad de prueba, tal que ySe cumple que
ρ r ≥0 ρ r dr =N
0 υE E ρ r Análogo al principio variacional Para la función de onda
Teoremas de Hohenberg y Kohn:
Para encontrar máximos mínimos de una función determinada f(x) sujeta a restricciones del tipo g(x)=0 se usan los multiplicadores de Lagrange.
Se define la función Lagrangiana:
Y encontrar los puntos donde
xgxfxL ,
0
f x
x
equivale a encontrar los puntos donde
0
L x,
x
0L x,
y
Ejemplo: Determinar los extremos de la función x3+y3, si se cumple que x2+y2-1=0
1,
,22
33
yxyxg
yxyxf
2233
,,,,
yxxx
yxgyxfyxL
2
2
2 2
3 2 0
3 2 0
1 0
Lx x
xL
y yy
Lx y
0 xgxf
Teoría de Funcionales de la Densidad (DFT) Química Cuántica
En nuestro caso queremos minimizar E[ρ], y tiene que cumplirse que
o sea: 0r dr n
L x, f x g x
L , E r dr n
f E
g r dr n
0E r dr n
0E r dr n
Potencialquímico
HKE F
rr r
Ecuación fundamental
de la TFD
Teoría de Funcionales de la Densidad (DFT) Química Cuántica
1
1
2
e e en n
iH
j iji j
Ftotal
n
ii
J KE H
eTF eT r r d VrE
Comparemos:
La ρ tiene que ser n-representableTiene que poder obtenerse a partir de alguna función de onda antisimétricaEsto se cumple si:
De este modo la ρ puede representarse en términos de n orbitales ortonormales y por tanto se puede generar a partir de una función de onda monodeterminante:
con
drrndrrr22/10
2
1
M
k kk
x x
1
0 1M
k kk
n M n
Existe una correspondencia entre la densidad electrónica y la función de ondadel estado base
Teoría de Funcionales de la Densidad (DFT) Química Cuántica
Sabemos que se puede obtener la densidad electrónica del estado base a partir de la correspondiente función de onda
Pero ….. ¿Podemos obtener la función de onda a partir de la densidad?
Dificultad: existe un número infinito de funciones de onda antisimétricas (no necesariamente correspondientes al estado base) que pueden dar la misma densidad
¿Cómo distinguir entonces la verdadera función del estado base?
De acuerdo con el principio variacional:
donde
Como υ(r) es un funcional de la densidad:
De modo que de todas las funciones de onda que dan la misma densidad del estado base (ρ), la correcta será aquella que minimice el valor esperado
0 0*r r r dr
000ˆˆ EHH
n
iiee rVTH
1
ˆˆˆ Hamiltoniano del sistema de n electrones
0 0ee eeˆ ˆ ˆ ˆT V r r dr T V r r dr
0 0ee eeˆ ˆ ˆ ˆT V T V
eeVT ˆˆ
Teoría de Funcionales de la Densidad (DFT) Química Cuántica
Como:
Lo que constituye un criterio de búsqueda restringida para F[ρ]
i.e. de todas las funciones de onda antisimétricas que dan la densidad ρ,
F[ρ] entrega solamente la correspondiente al menor valor esperado de
La búsqueda variacional es restringida ya que solamente se incluyen aquellas ψ que dan densidad ρ, en contraste con HF que incluye funciones de onda en todo el espacio de Hilbert de n partículas.
El procedimiento de minimización variacional puede dividirse en 2 pasos:
ee
ee
VT
VTF
ˆˆmin
ˆˆ00
eeVT ˆˆ
01
1
n
ee ii
n
ee ii
ee
ˆ ˆE min T V r
ˆ ˆmin min T V r
ˆ ˆmin min T V r r dr
min F r r dr
Teoría de Funcionales de la Densidad (DFT) Química Cuántica
Recordemos:
El funcional F[ρ] es universal en el sentido de que no depende del campo externo . De modo que la misma expresión de F[ρ] aplica para una molécula de hidrógeno o una proteína. La parte que sí depende la estructura del sistema en cuestión es el campo externo:
Esto sugiere la necesidad de encontrar un método de buena exactitud para tratar F[ρ]
Hasta este punto solo sabemos que la expresión correcta deberá contener términos tanto correspondientes a la energía cinética electrónica como a la repulsión entre electrones:
Este es un funcional muy difícil de aproximar, lo cual es lógico ya que hemos cambiado la incógnita variacionalmente desconocida: de la función de onda de muchas variables, a una función simple de 3 variables: la densidad electrónica
neV r r dr
eeˆ ˆF T V
Teoría de Funcionales de la Densidad (DFT) Química Cuántica
Recordando
Conjuntos estadísticos:
Conjunto microcanónico: sistema aislado con E, N y V constantes
Conjunto canónico: sistema en equilibrio térmico que intercambia energía con el medio, es una mezcla de estados puros todos con = número de partículas
Conjunto gran canónico: sistema en equilibrio térmico que intercambia energía y partículas con el medio, es una mezcla de estados puros todos donde existe una probabilidad ≠ 0 asociada a la existencia de estados con ≠ número de partículas pero de modo tal que el N promedio sea igual a un valor observado
N, T y V ctes
µ, T y V ctes
Teoría de Funcionales de la Densidad (DFT) Química Cuántica
Un sistema en un estado puro puede describirse por medio de una función de onda, un sistema en un estado mezclado NO
Un sistema en un estado mezclado está caracterizado por una distribución de probabilidad sobre todos los estados puros accesibles
El operador densidad se define como:
Generalizando la expresión anterior para un conjunto:
ˆN N N Y contiene la misma información que la función de onda
ˆi i i
i
p donde la sumatoria corre sobre todos los estados puros accesibles del sistema, de modo que:
En un conjunto microcanónico como todos los estados tienen igual energía, existe igual probabilidad de que el sistema se encuentre en cualquiera de ellos
Si η(E) representa el número de estados accesibles con energía E
Pero estamos interesados en los conjuntos canónicos y gran canónicos sistemas con T constante
0 1i ii
p p y
1
ˆ ipE
Teoría de Funcionales de la Densidad (DFT) Química Cuántica
Recordando
La traza (Tr) de un operador es igual a la suma de sus autovalores
Para un estado mezclado el valor esperado del observable Ŷ viene dado por:
ˆ ˆˆi i i
i
Y Tr Y p Y Ejemplo:
Observable: energía
Operador: Hamiltoniano
ˆ ˆˆi i i
i
E Tr H p H
Teoría de Funcionales de la Densidad (DFT) Química Cuántica
Para conjuntos canónicos:
ˆN Ni Ni Ni
i
p El principio variacional queda en función de la Energía Libre de Helmholtz (A):
0N NA A
donde
NA A F r r dr
F T U TS
con
Función de partición
0 1
1
N N
B
HN
A ln Z
k T
Z Tr e
Energía interna
NNB
N
N
TrkS
UTrU
TTrT
ˆlnˆ
ˆˆ
ˆˆ
Teoría de Funcionales de la Densidad (DFT) Química Cuántica
Para conjuntos gran canónicos:
ˆN Ni Ni Ni
N i
p
El principio variacional queda en función del gran potencial (Ω):
0
donde
1 ˆ ˆˆ ˆTr ln H N
con
Función de partición
0 1
1
N
B
ˆ ˆH N
N
ln Z
k T
Z Tr e
Potencialquímico
F r r dr
Teoría de Funcionales de la Densidad (DFT) Química Cuántica
Sentido físico del potencial químico DFT:
Mide la tendencia de un electrón de escapar de una nube electrónica
Es constante a través de todo el espacio, para el estado base de un átomo, molécula o sólido
Es igual a la pendiente de la curva de E vs N, a υ(r) constante
La analogía con el potencial químico de la termodinámica macroscópica clásica es claro y útil
Es igual en magnitud pero de signo contrario a la electronegatividad de Pauling y Mulliken:
totnPTi
clásicoDFT
nGE
,,
2
PI AE
La DFT es una una teoría de estados electrónicos base en equilibrio en la que en el principio variacional la electronegatividad química juega el mismo papel que la energía en ab-initio
Esto hace pensar en DFT como una descripción de la química de los sistemas
0 0GCH E E r dr
Teoría de Funcionales de la Densidad (DFT) Química Cuántica
Potencial químico:
NE
Electronegatividad
Calculables a partir de la DFT o determinables a partir de experimentos
La diferencia en electronegatividades conlleva a la transferencia electrónica
Los electrones tienden a fluir de regiones de baja electronegatividad a regiones de alta electronegatividad (de regiones de alto potencial químico a regiones de bajo potencial químico)
El número de electrones que fluye es proporcional a la diferencia de electronegatividades
La energía de estabilización es proporcional a su cuadrado
En una molécula las electronegatividades de los átomos o grupos que la forman se ecualiza, de modo que todos los fragmentos terminan teniendo igual electronegatividad (Principio de Sanderson)
El potencial químico de una especie microscópica en su estado base es una propiedad del sistema que se comporta forma muy similar al potencial químico termodinámico de sistemas macroscópicos
Teoría de Funcionales de la Densidad (DFT)
Potencial químico:
E
N
Química Cuántica
Recordemos:
E ρ r , r La energía es un funcional de la densidad y del potencial externo
r dr N Y como
E N , r La energía de un sistema se puede describir como función del número de electrones y funcional del potencial externo
DFT E
r
2
2
1 1
2 2
E
N N
Dureza del sistema molecular
Nr
rrf
N
Función de Fukui
Las variaciones de la energía con respecto a N (propiedad global) generan un conjunto de propiedades globales que dan cuenta de la reactividad del sistema.
Las variaciones de la energía con respecto al potencial externo (propiedad local) genera un conjunto de propiedades locales que dan cuenta de la selectividad y que por lo tanto está relacionada con la topología del sistema.
Teoría de Funcionales de la Densidad (DFT) Química Cuántica
Potencial químico: E
N
Para valores enteros de N es necesario usar la aproximación de diferencias finitas ya que la derivada de la energía con respecto a N es discontinua:
donde I = {E(N-1) – E(N)} es el primer potencial de ionización;
A = {E(N+1)-E(N)} es la afinidad electrónica.
1
2I A
Teoría de Funcionales de la Densidad (DFT) Química Cuántica
Dureza y blandura:
La primera derivada del potencial químico se conoce como dureza absoluta de una especie química (η) y su inverso como blandura (S):
2
2
1 1 1
2 2 2
ES
N N
2
IP AE
Desde un punto de vista práctico puede obtenerse como:
Y para aislante o semiconductores la expresión anterior aproxima a:
2LUMO HOMOE E
Desde un punto de vista químico las derivadas primera y segunda de la energía son las que revisten mayor importancia
La Dureza Química es una propiedad global del sistema y mide la resistencia impuesta por este al cambio en su distribución electrónica. En este contexto la dureza es un descriptor de la reactividad.
Teoría de Funcionales de la Densidad (DFT) Química Cuántica
Dureza y Blandura:
Para valores enteros de N es necesario usar la aproximación de diferencias finitas ya que la derivada de la energía con respecto a N es discontinua:
donde I = {E(N-1) – E(N)} es el primer potencial de ionización;
A = {E(N+1)-E(N)} es la afinidad electrónica.
1
2I A
2
2
1 1
2 2
E
N N
1
2S
Teoría de Funcionales de la Densidad (DFT) Química Cuántica
Indice Global de Electrofilicidad
Para un sistema electrofílico definido por {μ, η} y que está inmerso en un mar de electrones con potencial químico igual a cero.
La variación de energía del sistema hasta segundo orden con potencial externo constante, viene dado por:
El punto de saturación ocurre cuando:
entonces:
de donde se obtiene que la cantidad máxima de electrones que puede intercambiar el sistema es:
21
2E N N
0E / N
0E
NN
máxN
Substituyendo Nmax en la ecuación para E se obtiene:
que es el índice de electrofilicidad, una propiedad global del sistema.
El índice de electrofilicidad es entonces una medida de la estabilización energética del sistema cuando se satura de electrones que provienen del medio externo.
2
2
2
2E
Teoría de Funcionales de la Densidad (DFT) Química Cuántica
Los términos de dureza y blandura se utilizan en el vocabulario químico desde 1950, aunque el trabajo que introdujo estos términos se debe a Pearson en 1963. Sin embargo un consenso general en su significado no se alcanzó hasta 193 (Pearson)
Base Blanda: Especie donora de electrones con alta polarizabilidad y baja electronegatividad (alto potencial químico), fácilmente oxidable. Se asocia con HOMO de alta energía
Base Dura: Especie donora de electrones con baja polarizabilidad y alta electronegatividad, difícil de oxidar.
Acido Blando: Especie aceptora de elctrones con pequéña carga positiva. Sus electrones externos son fácilmente exitables.LUMO bajo
Acido Duro: Especie aceptora de electrones cuyos electrones externos no son fácilmente exitables.
Los ácidos duros interactúan preferentemente con bases duras dando lugar a complejos con enlace esencialmente iónico
Los ácidos blandos interactúan preferentemente con bases blandas dando lugar a complejos con enlace esencialmente covalentes (interacciones HOMO LUMO)
Teoría de Funcionales de la Densidad (DFT) Química Cuántica
Para enlaces formados a través de transferencia de electrones en reacciones ácido-base de Lewis
El número de electrones transferidos puede calcularse según:
Y la contribución a la energía debida a la transferencia de electrones según:
Donde rxn implica reacción del tipo: A+B→AB
Que representa una transferencia electrónica de la base de Lewis al ácido de Lewis
BA
ABN
2
2B Arxn
ABA B
E
Teoría de Funcionales de la Densidad (DFT) Química Cuántica
Hasta aquí el análisis de reacciones químicas a través de los conceptos de dureza y blandura, con expresiones en función de y solamente se refería a reacciones del tipo rxn:
Lo cual es muy limitado, sin embargo puede extenderse a esquemas de reacción más generales
Reacciones de intercambio de ácido:
Reacciones de intercambio de base:
Reacciones de doble intercambio:
B B B B
D D D D
A B A B
A B A B
B B D B B D
D D B D D B
B A B A B B
B A B A B B
D B B D B B D DA B A B A B A B
B D B B B D
D B D D D B
A A B A B A
A A B A B A
Teoría de Funcionales de la Densidad (DFT) Química Cuántica
1
2
B D B B B D
D B D D D B
A A B A B A a
A A B A B A a
Reacciones de intercambio de ácido:
2
2B Arxn
ABA B
E
Para encontrar una expresión la contribución a la energía debida a la transferencia de electrones reacciones de transferencia de ácido, reescribimos la ecuación a1 como la resta de las reacciones tipo rxn:
1
2
B B B B
D B D B
A B A B rxn
A B A B rxn
Ya conocemos la expresión de la contribución a la energía debida a la transferencia de electrones para reacciones rxn:
2
1
2B B
B B
B B
B ArxnA B
A B
E
2
2
2B D
D B
D B
B ArxnA B
A B
E
de modo que (rxn1)-(rxn2)=(a1) y la contribución a la E que estamos buscando será:
1 2
1 B B D B
rxn rxn rxnA B A BaE E E
La reacción será exotérmica si:
B B B B
D B B D
A B B A
A B B A
2 2
2 2B D B B
D B B B
B A B A
A B A B
Como el término de la izquierda siempre es <1, si
La reacción será exotérmica
De donde se confirma que la tendencia de los ácidos blandos a combinarse con bases blandas
B DA A
Teoría de Funcionales de la Densidad (DFT) Química Cuántica
2D B D D D BA A B A B A a Reacciones de intercambio de ácido:
Para encontrar una expresión la contribución a la energía debida a la transferencia de electrones reacciones de transferencia de ácido, reescribimos la ecuación a1 como la resta de las reacciones tipo rxn:
3
4
D D D D
B D B D
A B A B rxn
A B A B rxn
2
3
2D D
D D
D D
B ArxnA B
A B
E
2
4
2D B
D B
B D
B ArxnA B
A B
E
de modo que (rxn3)-(rxn4)=(a2) y la contribución a la E que estamos buscando será:
3 4
2 D D B D
rxn rxn rxnA B A BaE E E
La reacción será exotérmica si:
D D D D
B D D B
A B B A
A B B A
2 2
2 2D B D D
B D D D
B A B A
A B A B
Como el término de la izquierda siempre es >1, si
La reacción será exotérmica
Los ácidos duros desplazan a los blandos solo si son significativamente más fuertes que los blandos
B DA A
En los otros casos no cumple el principio HSAB
Teoría de Funcionales de la Densidad (DFT) Química Cuántica
Una lógica similar aplica para las reacciones de intercambio de base:
Serán exotérmicas (favorables) si:
B B D B B D D D B D D BB A B A B B B A B A B B
D D D D
D B B D
A B B A
A B B A
B B B B
B D D B
A B B A
A B B A
<1
Siempre que:
La reacción será exotérmica
La reacción seguirá el principio HSAB
De donde se confirma que la tendencia de las bases blandas a combinarse con acidos blandos
D BB B
>1
Siempre que:
La reacción será exotérmica
La reacción seguirá el principio HSAB
De donde se confirma que la tendencia de las bases blandas a combinarse con acidos blandos depende de sus pot. quimicos
D BB B
Teoría de Funcionales de la Densidad (DFT) Química Cuántica
Reacciones de doble intercambio:
D B B D B B D DA B A B A B A B
1
2
3
4
B B B B
D D D D
D B D B
B D B D
A B A B rxn
A B A B rxn
A B A B rxn
A B A B rxn
1 2 3 4
1 1
B B D D D B B D
rxn rxn rxn rxn rxnA B A B A B A Bab
rxn rxn rxnab a a
E E E E E
E E E
2 2
2 2B D B B
D B B B
B A B Arxnab
A B A B
E
2 2
2 2D B D D
B D D D
B A B A
A B A B
Teoría de Funcionales de la Densidad (DFT) Química Cuántica
La función de Fukui como índice de reactividad:
La función de Fukui es también derivada del potencial químico:
Nr
rrf
N
Involucra la densidad electrónica del átomo en su frontera, en la región de valencia
representa el cambio de la densidad electrónica cuando el número de electrones cambia.
representa la sensibilidad del potencial químico frente a una perturbación externa.
Una especie química que reacciona con un nucleófilo recibe electrones
Cuando la especie reacciona con un electrofilo pierde electrones
Entonces la variación de la densidad electrónica cuando se agregan o se quitan electrones definirá la capacidad electrofílica o nucleofílica de una especie química.
E+NNu-N
Teoría de Funcionales de la Densidad (DFT) Química Cuántica
El problema de la discontinuidad en N se resuelve introduciendo las derivadas por la derecha y por la izquierda. En general el valor de f(r) varía dependiendo de que se calcule la derivada acercandonos por la derecha, por la izquierda o como promedio. De este modo existen 3 índices de Fukui que aproximan como:
Mide la reactividad hacia un agente nucleofílico
Mide la reactividad hacia un agente electrofílico
Mide la reactividad hacia un agente radicálico
0 1
2 HOM
L
O LUM
H
U
O
MO
OMO
f r
f r r
f
r
r r
r
Entonces los puntos de mayor reactividad frente a una perturbación externa son aquellos que presentan un valor grande de f(r)....
...A mayor f(r) mayor reactividad de sitio ...
Puesto que f(r) representa una variación del potencial químico inducida por una perturbación externa, mientras mayor sea esta variación más reactivo será el sistema en el punto r. Se generará un gradiente de potencial químico que inducirá una mayor transferencia de carga.
derecha
izquierda
Teoría de orbitales frontera Cuando un el sistema reacciona con un electrófilo los electrones involucrados en la transferencia provienen del HOMO mientras que frente a un nucleófilo, los electrones entran al orbital LUMO
Son siempre positivas
Teoría de Funcionales de la Densidad (DFT) Química Cuántica
La condensación sobre sitios específicos en la molécula produce los índices de Fukui:
Entonces según la aproximación de diferencias finitas, las funciones de Fukui de un centro atómico x se calculan como diferencias entre las cargas netas en el centro atómico de la molécula neutra y del respectivo ión que se forma cuando se agrega o se retira un electrón.
0
1
1
1 1
2 2
x x x
x x x
x xx xx
f q N q N
f q N q N
q N q Nf ff
Teoría de Funcionales de la Densidad (DFT) Química Cuántica
Otro índice de selectividad es: x x xf f f
Si f es mayor que cero el sitio es electrofílico (rojo)Si f es menor que cero el sitio será nucleofílico (amarillo).
Anilina Fenol
W. Kong, L. J. Sham, Phys. Rev. 1965, 1133, A140.
Mostraron como se puede obtener la a partir de la solución de un conjunto de ecuaciones tipo HF y presentaron la relación matemática entre E y
Sistema de referencia (s) con e- que no interactuan:
1
21
2si
i s i
n
rH
KSih Hamiltoniano Kohn-Sham unielectrónico
Sistema real con e- que sí interactuan:
1
n
i eei
ˆ ˆ ˆH T r V
Teoría de Funcionales de la Densidad (DFT)
xc eeE T V 1
2Total
selectrónica Coulombicae ne ee xcE T V V E
1 21 2
12
1
2
s
ee ee
T T T
r rV V dr dr
r
La calidad de los resultados DFT depende de cómo se defina el término Vxc
xcxc
EV
Teoría de Funcionales de la Densidad (DFT)
CXXC EEE rFEXC
Generalized Gradient Approximation (GGA)Toma en cuenta la densidad electrónica y su gradiente (variación
de la densidad con la posición.Funciona bien para sólidos y para moléculas
drxgrEGGAX
3/4 3/4 x
Local Density Approximation (LDA)La densidad se modela como la de un gas electrónico homogéneo con densidad electrónica Funciona bien para sólidos pero no para moléculas
drrcE XLDAX
3/4
(subestima 10%)3/1
43
23
Xc
Si: 3/1
43
23
xgLDAX
GGAX EE
Funcionales de Intercambio, Ex:
-Slater (LSDA): coeficiente teórico, Local, (< 10%)
-Xa: coeficiente empírico, se usa cuando no se incluye Ec , Local
-Becke88: Local + correcciones que incluyen gradiente
-Perdew-Wang 91: incluye gradiente, mejora a Becke88 para límites de pequeñas y grandes x
...6
sinh61422
1
288
xbbx
xbxbx
xg Becke
21 1lnsinh xxx
-Barone: Modificación de PW91, gradiente (diferentes valores de los parámetros)
....00375.0
/10sinh61
10 42461
466455.12291
2
xxxxbx
xexbbxxg
xPW
3/4 x
-Gill96: Incluye gradiente 2/396 bxxg G
drrE LDA
X3/4
3/1
43
23
drrbE X
X3/4
3/1
43
89
Funcionales de Correlación, Ec:
-VWN: Vosko, Wilk, and Nusair, Local Spin Density (LSD) Correlation-LYP: Lee, Yang, and Parr, Local y Gradiente -Perdew81: Local-Perdew86: Local + correcciones de gradiente-PW91: Perdew-Wang, incluye correcciones de gradiente-Becke96: incluye correcciones de gradiente
Contiene la parte no coulómbica de Vee + la diferencia entre T y T(s)
OpC C C CE E E E
220 01 10
0
2 22
2 2VWNC
x x b xbxx b Q QA ln tan ln tan
X x Q x b X x X x Q x b
0x ,b,c son parámetros 1 22 2 1 24/ /X x x bx c Q c b x
11
1 3
8 3 8 3
2
22 2
2 22 2 22
3
2 2 2
4
147 7 5 1
218 18 2 18
11 2 2 2
9 3 3 3
LYPC
/ /F
/
/
E ad
ab C
a, b, c, d parámetros 1 3 1 3
2 311 3 1 3 21 3 1 3
33
1 1 10
/c /// /
F/ /
e dc C
d d
Funcionales más comunes:
B3LYP: Funcional Híbrido
VWNC
LYPC
BeckeX
HFX
SlaterX ECCEBEEAAE 11 88
A, B y C son parámetros obtenidos empíricamente de modo que los resultados ajusten 56 energies de atomización, 42 potenciales de ionización, 8 afinidades protónicas, y 10 energías atómicas correspondientes a elementos de la 1ra fila A=0.80, B=0.72 y C=0.81
86881 PC
VWNC
BeckeX
HFX
SlaterX CEEBEEAAE
B3P86: Funcional Híbrido
A, B y C = B3LYP, la Ec no-local correspondiente a Perdew86
91881 PWC
VWNC
BeckeX
HFX
SlaterX CEEBEEAAE
B3PW91: Funcional Híbrido
A, B y C = B3LYP, la Ec no-local correspondiente a Perdew-Wang 91
B1B96: Funcional Híbrido, de un solo parámetro
cop= 0.0031 c= c= 0.038
Funcionales más comunes:
BHandHLYP: Funcional Híbrido,
88
0 5 0 5
0 5 0 5 0 5
HF LSDA LYPX X C
HF LSDA Becke LYPX X X C
. E . E E
. E . E . E E
Funcionales definidos por el usuario:
2 1 4 3 6 5HF Slater no local Local no localX X X C CP E P P E P E P E P E
Teoría de Funcionales de la Densidad (DFT) Química Cuántica
Hydrogen molecules - using the LSDA
Teoría de Funcionales de la Densidad (DFT) Química Cuántica
Teoría de Funcionales de la Densidad (DFT) Química Cuántica
Teoría de Funcionales de la Densidad (DFT) Química Cuántica
Water dimer
Fe2SiO4
Teoría de Funcionales de la Densidad (DFT) Química Cuántica
Teoría de Funcionales de la Densidad (DFT) Química Cuántica
Teoría de Funcionales de la Densidad (DFT) Química Cuántica
Teoría de Funcionales de la Densidad (DFT) Química Cuántica
Teoría de Funcionales de la Densidad (DFT) Química Cuántica
Teoría de Funcionales de la Densidad (DFT) Química Cuántica
Teoría de Funcionales de la Densidad (DFT) Química Cuántica
Teoría de Funcionales de la Densidad (DFT) Química Cuántica
Ventajas con respecto a cálculos ab-initio
La función de onda de un sistema de N electrones tiene 3N variablesLa densidad electrónica tiene 3 variables (x,y,z) para cualquier sistema sin importar su tamaño DFT reduce considerablemente el costo computacional (N3 máx, MP2: N4, CI: N6 ó N7).
Se obtienen propiedades de significado químico como electronegatividad (potencial químico), dureza (blandura), funciones de Fukui, etc
Desventajas con respecto a cálculos ab-initio
No existe una forma universal conocida para expresar la E como funcional de la densidad. La estrategia a seguir es emplear aproximaciones, pero no hay un modo sistemático de mejorar los resultados obtenidos ( a diferencia de ab-initio).
DFT está formulada para sistemas en su estado base, por lo que su extensión a estados excitados no es obvia.
Teoría de Funcionales de la Densidad (DFT)
top related