stringtranslate.com

Teorema de Koopmans

El teorema de Koopmans establece que en la teoría de Hartree-Fock de capa cerrada (HF), la primera energía de ionización de un sistema molecular es igual al negativo de la energía orbital del orbital molecular más alto ocupado ( HOMO ). Este teorema recibe su nombre en honor a Tjalling Koopmans , quien publicó este resultado en 1934. [1]

El teorema de Koopmans es exacto en el contexto de la teoría restringida de Hartree-Fock si se supone que los orbitales del ion son idénticos a los de la molécula neutra (la aproximación de orbital congelado [2] ). Las energías de ionización calculadas de esta manera están en concordancia cualitativa con el experimento: la primera energía de ionización de moléculas pequeñas a menudo se calcula con un error de menos de dos  electronvoltios . [3] [4] [5] Por lo tanto, la validez del teorema de Koopmans está íntimamente ligada a la precisión de la función de onda subyacente de Hartree-Fock . [ cita requerida ] Las dos principales fuentes de error son la relajación orbital, que se refiere a los cambios en el operador de Fock y los orbitales de Hartree-Fock al cambiar el número de electrones en el sistema, y ​​la correlación electrónica , que se refiere a la validez de representar toda la función de onda de muchos cuerpos utilizando la función de onda de Hartree-Fock, es decir, un único determinante de Slater compuesto de orbitales que son las funciones propias del operador de Fock autoconsistente correspondiente.

Las comparaciones empíricas con valores experimentales y cálculos ab initio de mayor calidad sugieren que en muchos casos, pero no en todos, las correcciones energéticas debidas a los efectos de relajación casi cancelan las correcciones debidas a la correlación electrónica. [6] [7]

Existe un teorema similar (teorema de Janak) en la teoría del funcional de la densidad (DFT) para relacionar la primera energía de ionización vertical exacta y la afinidad electrónica con las energías HOMO y LUMO , aunque tanto la derivación como la declaración precisa difieren de las del teorema de Koopmans. [8] Las energías de ionización calculadas a partir de las energías orbitales DFT son generalmente más pobres que las del teorema de Koopmans, con errores mucho mayores que dos electronvoltios posibles dependiendo de la aproximación de correlación de intercambio empleada. [3] [4] La energía LUMO muestra poca correlación con la afinidad electrónica con aproximaciones típicas. [9] El error en la contraparte DFT del teorema de Koopmans es un resultado de la aproximación empleada para el funcional de energía de correlación de intercambio de modo que, a diferencia de la teoría HF, existe la posibilidad de mejores resultados con el desarrollo de mejores aproximaciones.

Generalizaciones

Si bien el teorema de Koopmans se estableció originalmente para calcular energías de ionización a partir de funciones de onda Hartree-Fock restringidas (de capa cerrada), el término desde entonces adquirió un significado más generalizado como una forma de usar energías orbitales para calcular cambios de energía debidos a cambios en la cantidad de electrones en un sistema.

Iones en estado fundamental y estado excitado

El teorema de Koopmans se aplica a la eliminación de un electrón de cualquier orbital molecular ocupado para formar un ion positivo. La eliminación del electrón de diferentes orbitales moleculares ocupados hace que el ion se encuentre en diferentes estados electrónicos. El más bajo de estos estados es el estado fundamental y, a menudo, aunque no siempre, esto surge de la eliminación del electrón del HOMO. Los otros estados son estados electrónicos excitados.

Por ejemplo, la configuración electrónica de la molécula de H 2 O es (1a 1 ) 2 (2a 1 ) 2 (1b 2 ) 2 (3a 1 ) 2 (1b 1 ) 2 , [10] donde los símbolos a 1 , b 2 y b 1 son etiquetas orbitales basadas en la simetría molecular . Del teorema de Koopmans la energía del HOMO 1b 1 corresponde a la energía de ionización para formar el ion H 2 O + en su estado fundamental (1a 1 ) 2 (2a 1 ) 2 (1b 2 ) 2 (3a 1 ) 2 (1b 1 ) 1 . La energía del segundo MO más alto 3a 1 se refiere al ion en el estado excitado (1a 1 ) 2 (2a 1 ) 2 (1b 2 ) 2 (3a 1 ) 1 (1b 1 ) 2 , y así sucesivamente. En este caso, el orden de los estados electrónicos de los iones corresponde al orden de las energías orbitales. Las energías de ionización en estado excitado se pueden medir mediante espectroscopia fotoelectrónica .

Para H 2 O, las energías orbitales cercanas a Hartree-Fock (con signo cambiado) de estos orbitales son 1a 1 559,5, 2a 1 36,7 1b 2 19,5, 3a 1 15,9 y 1b 1 13,8 eV . Las energías de ionización correspondientes son 539,7, 32,2, 18,5, 14,7 y 12,6 eV. [10] Como se explicó anteriormente, las desviaciones se deben a los efectos de la relajación orbital, así como a las diferencias en la energía de correlación electrónica entre el estado molecular y los diversos estados ionizados.

En cambio, en el caso del N 2 , el orden de las energías orbitales no es idéntico al orden de las energías de ionización. Los cálculos de Near-Hartree-Fock con un conjunto de bases amplio indican que el orbital de enlace 1π u es el HOMO. Sin embargo, la energía de ionización más baja corresponde a la eliminación de un electrón del orbital de enlace 3σ g . En este caso, la desviación se atribuye principalmente a la diferencia en la energía de correlación entre los dos orbitales. [11]

Para afinidades electrónicas

A veces se afirma [12] que el teorema de Koopmans también permite el cálculo de las afinidades electrónicas como la energía de los orbitales moleculares desocupados más bajos ( LUMO ) de los respectivos sistemas. Sin embargo, el artículo original de Koopmans no hace ninguna afirmación con respecto a la importancia de los valores propios del operador de Fock distintos del correspondiente al HOMO . No obstante, es sencillo generalizar la afirmación original de Koopmans para calcular la afinidad electrónica en este sentido.

Los cálculos de afinidades electrónicas que utilizan este enunciado del teorema de Koopmans han sido criticados [13] sobre la base de que los orbitales virtuales (desocupados) no tienen interpretaciones físicas bien fundamentadas y que sus energías orbitales son muy sensibles a la elección del conjunto de bases utilizado en el cálculo. A medida que el conjunto de bases se vuelve más completo, aparecerán cada vez más orbitales "moleculares" que no están realmente en la molécula de interés y se debe tener cuidado de no utilizar estos orbitales para estimar las afinidades electrónicas.

Las comparaciones con experimentos y cálculos de mayor calidad muestran que las afinidades electrónicas predichas de esta manera son generalmente bastante pobres.

Para sistemas de capa abierta

El teorema de Koopmans también es aplicable a sistemas de capa abierta, sin embargo, las energías orbitales (valores propios de las ecuaciones de Roothaan) deben corregirse, como se mostró en la década de 1970. [14] [15] A pesar de este trabajo temprano, la aplicación del teorema de Koopmans a sistemas de capa abierta continuó causando confusión, por ejemplo, se afirmó que el teorema de Koopmans solo se puede aplicar para eliminar el electrón desapareado. [16] Más tarde, se revisó la validez del teorema de Koopmans para ROHF y se informaron varios procedimientos para obtener energías orbitales significativas. [17] [18] [19] [20] Las energías orbitales de espín hacia arriba (alfa) y espín hacia abajo (beta) no necesariamente tienen que ser las mismas. [21]

Contraparte en la teoría del funcional de la densidad

La teoría funcional de la densidad de Kohn-Sham (KS) (KS-DFT) admite su propia versión del teorema de Koopmans (a veces llamado teorema de Koopmans-DFT ) muy similar en espíritu a la teoría de Hartree-Fock. El teorema iguala la primera energía de ionización (vertical) de un sistema de electrones al negativo de la energía HOMO KS correspondiente . De manera más general, esta relación es verdadera incluso cuando el sistema KS describe un conjunto de temperatura cero con un número no entero de electrones para números enteros y . Al considerar los electrones, la carga excedente infinitesimal ingresa al LUMO KS del sistema electrónico N pero luego el potencial KS exacto salta por una constante conocida como la "discontinuidad derivada". [22] Se puede argumentar que la afinidad electrónica vertical es exactamente igual al negativo de la suma de la energía LUMO y la discontinuidad derivada. [22] [23] [24] [25]

A diferencia del estado aproximado del teorema de Koopmans en la teoría de Hartree Fock (debido a la negligencia de la relajación orbital), en la aplicación KS exacta el teorema es exacto, incluido el efecto de la relajación orbital. Una prueba esquemática de esta relación exacta se realiza en tres etapas. Primero, para cualquier sistema finito determina la forma asintótica de la densidad, que decae como . [22] [26] Luego, como corolario (ya que el sistema que interactúa físicamente tiene la misma densidad que el sistema KS), ambos deben tener la misma energía de ionización. Finalmente, dado que el potencial KS es cero en el infinito, la energía de ionización del sistema KS es, por definición, el negativo de su energía HOMO, es decir, . [27] [28]

Si bien estas son afirmaciones exactas en el formalismo de DFT, el uso de potenciales de intercambio-correlación aproximados hace que las energías calculadas sean aproximadas y, a menudo, las energías orbitales son muy diferentes de las energías de ionización correspondientes (¡incluso en varios eV!). [29]

Un procedimiento de ajuste es capaz de "imponer" el teorema de Koopmans en aproximaciones DFT, mejorando así muchas de sus predicciones relacionadas en aplicaciones reales. [29] [30] En las DFT aproximadas se puede estimar con un alto grado de precisión la desviación del teorema de Koopmans utilizando el concepto de curvatura de energía. [31] Proporciona energías de excitación de orden cero [32] y . [8] [33]

Imagen orbital dentro de los formalismos de muchos cuerpos

El concepto de orbitales moleculares y una imagen similar a la de Koopmans de los procesos de ionización o de unión de electrones se pueden extender a funciones de onda de muchos cuerpos correlacionadas mediante la introducción de orbitales de Dyson. [34] [35] Los orbitales de Dyson se definen como la superposición generalizada entre una función de onda molecular de un electrón y la función de onda del electrón del sistema ionizado (o la función de onda del electrón de un sistema unido a un electrón):

Los orbitales canónicos de Hartree-Fock son orbitales de Dyson calculados para la función de onda de Hartree-Fock del sistema de electrones y la aproximación de Koopmans del sistema de electrones. Cuando se utilizan funciones de onda correlacionadas, los orbitales de Dyson incluyen efectos de correlación y relajación orbital. Los orbitales de Dyson contienen toda la información sobre los estados inicial y final del sistema necesaria para calcular cantidades observables experimentalmente, como las secciones eficaces de fotoionización/fotodesprendimiento total y diferencial.

Referencias

  1. ^ Koopmans, Tjalling (1934). "Über die Zuordnung von Wellenfunktionen und Eigenwerten zu den einzelnen Elektronen eines Atoms". Física . 1 (1–6): 104–113. Código Bib : 1934Phy.....1..104K. doi :10.1016/S0031-8914(34)90011-2.
  2. ^ Szabo, Attila ; Ostlund, Neil S. (1996). Química cuántica moderna. Dover Publications. pág. 128. ISBN 0-486-69186-1.OCLC 34357385  .
  3. ^ ab Politzer, Peter; Abu-Awwad, Fakher (1998). "Un análisis comparativo de las energías orbitales de Hartree-Fock y Kohn-Sham". Teoría de la química teórica: teoría, computación y modelado . 99 (2): 83–87. doi :10.1007/s002140050307. S2CID  96583645.
  4. ^ ab Hamel, Sebastien; Duffy, Patrick; Casida, Mark E.; Salahub, Dennis R. (2002). "Orbitales Kohn-Sham y energías orbitales: construcciones ficticias pero buenas aproximaciones de todos modos". Journal of Electron Spectroscopy and Related Phenomena . 123 (2–3): 345–363. doi :10.1016/S0368-2048(02)00032-4.
  5. ^ Véase, por ejemplo, Szabo, A.; Ostlund, NS (1982). "Capítulo 3". Química cuántica moderna . ISBN 978-0-02-949710-4.
  6. ^ Michl, Josef; Bonačić-Koutecký, Vlasta (1990). Aspectos electrónicos de la fotoquímica orgánica . Wiley. pag. 35.ISBN 978-0-471-89626-5.
  7. ^ Aquí, Warren J.; Radom, Leo; Schleyer, Paul vR; Pople, John A. (1986). Teoría de los orbitales moleculares ab initio . Wiley. pag. 24.ISBN 978-0-471-81241-8.
  8. ^ ab Janak, J. F (1978). "Prueba de que ∂E∂ni=ε en la teoría de los funcionales de la densidad". Physical Review B . 18 (12): 7165–7168. Bibcode :1978PhRvB..18.7165J. doi :10.1103/PhysRevB.18.7165.
  9. ^ Zhang, Gang; Musgrave, Charles B. (2007). "Comparación de métodos DFT para cálculos de valores propios de orbitales moleculares". The Journal of Physical Chemistry A . 111 (8): 1554–1561. Bibcode :2007JPCA..111.1554Z. doi :10.1021/jp061633o. PMID  17279730. S2CID  1516019.
  10. ^ ab Levine, IN (1991). Química cuántica (4.ª ed.). Prentice-Hall. pág. 475. ISBN 978-0-7923-1421-9.
  11. ^ Duke, Brian J.; O'Leary, Brian (1995). "Moléculas que no son de Koopman". Revista de educación química . 72 (6): 501. Bibcode :1995JChEd..72..501D. doi :10.1021/ed072p501.
  12. ^ Véase, por ejemplo, Szabo, A.; Ostlund, NS (1982). Química cuántica moderna . pág. 127. ISBN 978-0-02-949710-4.
  13. ^ Jensen, Frank (1990). Introducción a la química computacional . Wiley. págs. 64-65. ISBN 978-0-471-98425-2.
  14. ^ Sauer, Joaquín; Jung, Christoph (1975). "Consecuencias de los teoremas de Koopmansschen en el método restringido de Hartree Fock para sistemas de capa abierta". Teor. Chim. Acta . 40 (2): 129-141. doi :10.1007/BF01135884. S2CID  92462090.
  15. ^ Sauer, J.; Jung, Ch.; Jaffé, HH; Singerman, J. (1978-07-01). "Energías orbitales en sistemas de capas abiertas". The Journal of Chemical Physics . 69 (1): 495–496. Bibcode :1978JChPh..69..495S. doi :10.1063/1.436380. ISSN  0021-9606.
  16. ^ Fulde, Peter (1995). Correlaciones electrónicas en moléculas y sólidos . Springer. pp. 25-26. ISBN. 978-3-540-59364-5.
  17. ^ Plakhutin, BN; Gorelik, EV; Breslavskaya, NN (2006). "Teorema de Koopmans en el método ROHF: forma canónica para el hamiltoniano de Hartree-Fock". The Journal of Chemical Physics . 125 (20): 204110. Bibcode :2006JChPh.125t4110P. doi :10.1063/1.2393223. PMID  17144693.
  18. ^ Davidson, Ernest R. ; Plakhutin, Boris N. (2010). "Teorema de Koopmans en el método Hartree-Fock de capa abierta restringida. II. El segundo conjunto canónico para orbitales y energías orbitales". The Journal of Chemical Physics . 132 (18): 184110. Bibcode :2010JChPh.132r4110D. doi :10.1063/1.3418615. Archivado desde el original el 14 de enero de 2021 . Consultado el 29 de abril de 2018 .
  19. ^ Plakhutin, Boris N.; Davidson, Ernest R. (2009). "Teorema de Koopmans en el método Hartree−Fock de capa abierta restringida. 1. Un enfoque variacional†". The Journal of Physical Chemistry A . 113 (45): 12386–12395. Bibcode :2009JPCA..11312386P. doi :10.1021/jp9002593. PMID  19459641.
  20. ^ Glaesemann, Kurt R.; Schmidt, Michael W. (2010). "Sobre el ordenamiento de las energías orbitales en ROHF† de alto espín". The Journal of Physical Chemistry A . 114 (33): 8772–8777. Bibcode :2010JPCA..114.8772G. doi :10.1021/jp101758y. PMID  20443582.
  21. ^ Tsuchimochi, Takashi; Scuseria, Gustavo E. (2010). "Comunicación: teoría ROHF simplificada". The Journal of Chemical Physics . 133 (14): 141102. arXiv : 1008.1607 . Código Bibliográfico :2010JChPh.133n1102T. doi :10.1063/1.3503173. PMID  20949979. S2CID  31648260.
  22. ^ abc Perdew, John P.; Parr, Robert G.; Levy, Mel; Balduz, Jr., Jose L. (1982). "Teoría funcional de la densidad para el número fraccionario de partículas: discontinuidades derivadas de la energía". Physical Review Letters . 49 (23): 1691–1694. Código Bibliográfico :1982PhRvL..49.1691P. doi :10.1103/PhysRevLett.49.1691.
  23. ^ Perdew, John P.; Levy, Mel (1997). "Comentario sobre la "Importancia del valor propio de Kohn-Sham más alto ocupado"". Physical Review B . 56 (24): 16021–16028. Código Bibliográfico :1997PhRvB..5616021P. doi :10.1103/PhysRevB.56.16021.
  24. ^ Kleinman, Leonard (1997). "Significación del valor propio de Kohn-Sham ocupado más alto". Physical Review B . 56 (19): 12042. Bibcode :1997PhRvB..5612042K. doi :10.1103/PhysRevB.56.12042.
  25. ^ Kleinman, Leonard (1997). "Respuesta al "Comentario sobre 'Importancia del valor propio de Kohn-Sham más alto ocupado'"". Physical Review B . 56 (24): 16029. Código Bibliográfico :1997PhRvB..5616029K. doi :10.1103/PhysRevB.56.16029.
  26. ^ Almbladh, C. -O.; von Barth, U. (1985). "Resultados exactos para las densidades de carga y espín, potenciales de correlación de intercambio y valores propios funcionales de la densidad". Physical Review B . 31 (6): 3231–3244. Bibcode :1985PhRvB..31.3231A. doi :10.1103/PhysRevB.31.3231. PMID  9936207.
  27. ^ Savin, A; Umrigar, CJ; Gonze, Xavier (1998). "Relación de los valores propios de Kohn-Sham con las energías de excitación". Chemical Physics Letters . 288 (2–4): 391. Bibcode :1998CPL...288..391S. doi : 10.1016/S0009-2614(98)00316-9 .
  28. ^ Levy, Mel; Perdew, John P; Sahni, Viraht (1984). "Ecuación diferencial exacta para la densidad y la energía de ionización de un sistema de muchas partículas". Physical Review A . 30 (5): 2745. Bibcode :1984PhRvA..30.2745L. doi :10.1103/PhysRevA.30.2745.
  29. ^ ab Salzner, U.; Baer, ​​R. (2009). "Koopmans' springs to life" (PDF) . The Journal of Chemical Physics . 131 (23): 231101–4. Bibcode :2009JChPh.131w1101S. doi :10.1063/1.3269030. hdl : 11693/11792 . PMID  20025305. Archivado desde el original el 14 de enero de 2021 . Consultado el 24 de septiembre de 2019 .
  30. ^ Baer, ​​R.; Livshits, E.; Salzner, U. (2010). "Híbridos separados por rangos "ajustados" en la teoría del funcional de la densidad". Revisión anual de química física . 61 : 85–109. doi :10.1146/annurev.physchem.012809.103321. hdl : 11693/22326 . PMID  20055678.Stein, T.; Eisenberg, H.; Kronik, L.; Baer, ​​R. (2010). "Brechas fundamentales de sistemas finitos a partir de los valores propios de un método generalizado de Kohn-Sham". Physical Review Letters . 105 (26): 266802. arXiv : 1006.5420 . Bibcode :2010PhRvL.105z6802S. doi :10.1103/PhysRevLett.105.266802. PMID  21231698. S2CID  42592180. Kornik, L.; Stein, T.; Refaely-Abramson, S.; Baer, ​​R. (2012). "Brechas de excitación de sistemas de tamaño finito a partir de funcionales híbridos separados por rangos ajustados de forma óptima". Journal of Chemical Theory and Computation . 8 (5): 1515–31. doi : 10.1021/ct2009363 . PMID  26593646.
  31. ^ Stein, Tamar; Autschbach, Jochen; Govind, Niranjan; Kronik, Leeor; Baer, ​​Roi (2012). "Curvatura y energías orbitales de frontera en la teoría del funcional de la densidad". The Journal of Physical Chemistry Letters . 3 (24): 3740–4. arXiv : 1208.1496 . doi :10.1021/jz3015937. PMID  26291104. S2CID  22495102.
  32. ^ Andrejkovics, I; Nagy, Á (1998). "Energías de excitación en la teoría del funcional de la densidad: Comparación de varios métodos para las moléculas de H2O, N2, CO y C2H4". Chemical Physics Letters . 296 (5–6): 489. Bibcode :1998CPL...296..489A. doi :10.1016/S0009-2614(98)01075-6.
  33. ^ Perdew, John P; Levy, Mel (1983). "Contenido físico de las energías orbitales exactas de Kohn-Sham: brechas de banda y discontinuidades derivadas". Physical Review Letters . 51 (20): 1884. Bibcode :1983PhRvL..51.1884P. doi :10.1103/PhysRevLett.51.1884.
  34. ^ Krylov, Anna I. (24 de agosto de 2020). "De los orbitales a los observables y viceversa". The Journal of Chemical Physics . 153 (8): 080901. doi : 10.1063/5.0018597 . ISSN  0021-9606. PMID  32872858.
  35. ^ Ortiz, JV (18 de agosto de 2020). "Conceptos de orbitales de Dyson para la descripción de electrones en moléculas". The Journal of Chemical Physics . 153 (7): 070902. doi : 10.1063/5.0016472 . ISSN  0021-9606. PMID  32828082.

Enlaces externos