Métodos compuestos de química cuántica - Quantum chemistry composite methods

Los métodos compuestos de química cuántica (también conocidos como recetas termoquímicas) son métodos de química computacional que apuntan a una alta precisión al combinar los resultados de varios cálculos. Combinan métodos con un alto nivel de teoría y un pequeño conjunto de bases con métodos que emplean niveles más bajos de teoría con conjuntos de bases más grandes. Se utilizan comúnmente para calcular cantidades termodinámicas como entalpías de formación, energías de atomización, energías de ionización y afinidades electrónicas. Su objetivo es la precisión química, que generalmente se define dentro de 1 kcal / mol del valor experimental. El primer modelo químico sistemático de este tipo con amplia aplicabilidad se denominó Gaussian-1 (G1) introducido por John Pople . Este fue reemplazado rápidamente por el Gaussian-2 (G2) que se ha utilizado ampliamente. El Gaussian-3 (G3) se introdujo más tarde.

Teorías de Gauss-n

Gaussiano-2 (G2)

El G2 utiliza siete cálculos:

  1. la geometría molecular se obtiene mediante una optimización de MP2 utilizando el conjunto de bases 6-31G (d) y todos los electrones incluidos en la perturbación. Esta geometría se utiliza para todos los cálculos posteriores.
  2. El nivel más alto de teoría es un cálculo de interacción de configuración cuadrática con excitaciones simples y dobles y una contribución de excitación triple (QCISD (T)) con el conjunto de bases 6-311G (d). Dicho cálculo en los programas Gaussiano y Spartan también proporciona las energías MP2 y MP4 que también se utilizan.
  3. El efecto de las funciones de polarización se evalúa mediante un cálculo MP4 con el conjunto de bases 6-311G (2df, p).
  4. El efecto de las funciones difusas se evalúa mediante un cálculo MP4 con el conjunto de bases 6-311 + G (d, p).
  5. El conjunto de bases más grande es 6-311 + G (3df, 2p) utilizado en el nivel de teoría MP2.
  6. Una optimización de geometría de Hartree-Fock con el conjunto de bases 6-31G (d) utilizado para dar una geometría para:
  7. Un cálculo de frecuencia con la base 6-31G (d) establecida para obtener la energía vibratoria de punto cero (ZPVE)

Se supone que los diversos cambios de energía son aditivos, por lo que la energía combinada viene dada por:

EQCISD (T) de 2 + [EMP4 de 3 - EMP4 de 2] + [EMP4 de 4 - EMP4 de 2] + [EMP2 de 5 + EMP2 de 2 - EMP2 de 3 - EMP2 de 4]

El segundo término corrige el efecto de sumar las funciones de polarización. El tercer término corrige las funciones difusas. El término final corrige el conjunto de bases más amplio con los términos de los pasos 2, 3 y 4, lo que evita que las contribuciones se cuenten dos veces. Se hacen dos correcciones finales a esta energía. La ZPVE tiene una escala de 0,8929. Luego se agrega una corrección empírica para tener en cuenta los factores no considerados anteriormente. Esto se denomina corrección de nivel superior (HC) y viene dado por -0,00481 x (número de electrones de valencia) -0,00019 x (número de electrones de valencia no apareados). Los dos números se obtienen calibrando los resultados contra los resultados experimentales para un conjunto de moléculas. Se añaden la ZPVE a escala y la HLC para dar la energía final. Para algunas moléculas que contienen uno de los elementos de la tercera fila Ga – Xe, se agrega un término adicional para explicar el acoplamiento de la órbita de espín .

Se han utilizado varias variantes de este procedimiento. Eliminar los pasos 3 y 4 y confiar solo en el resultado de MP2 del paso 5 es significativamente más barato y solo un poco menos preciso. Este es el método G2MP2. A veces, la geometría se obtiene utilizando un método de teoría funcional de densidad como B3LYP y, a veces, el método QCISD (T) en el paso 2 se reemplaza por el método de clúster acoplado CCSD (T).

La variante G2 (+), donde el símbolo "+" se refiere a funciones difusas agregadas, describe mejor los aniones que la teoría G2 convencional. El conjunto de bases 6-31 + G (d) se utiliza en lugar del conjunto de bases 6-31G (d) tanto para la optimización de la geometría inicial como para la optimización de la segunda geometría y el cálculo de frecuencia. Además, la aproximación de núcleo congelado se realiza para la optimización inicial de MP2, mientras que G2 generalmente usa el cálculo completo.

Gaussiano-3 (G3)

El G3 es muy similar al G2 pero aprende de la experiencia con la teoría G2. El conjunto de bases 6-311G se reemplaza por la base 6-31G más pequeña. Los cálculos finales de MP2 utilizan un conjunto de bases más grande, generalmente llamado G3large, y correlacionando todos los electrones no solo los electrones de valencia como en la teoría G2, además se introduce un término de corrección de órbita de espín y una corrección empírica para los electrones de valencia. Esto da algunas contribuciones de correlación central a la energía final. El HLC adopta la misma forma pero con diferentes parámetros empíricos.

Gaussiano-4 (G4)

G4 es un método compuesto en el espíritu de las otras teorías gaussianas e intenta llevar la precisión lograda con G3X un pequeño paso más allá. Esto implica la introducción de un esquema de extrapolación para obtener las energías Hartree-Fock límite establecidas de base, el uso de geometrías y correcciones termoquímicas calculadas en el nivel B3LYP / 6-31G (2df, p), un cálculo de punto único de nivel más alto en CCSD (T ) en lugar del nivel QCISD (T), y la adición de funciones de polarización extra en los cálculos MP2 del conjunto de base más grande. Por lo tanto, la teoría gaussiana 4 (G4) es un enfoque para el cálculo de energías de especies moleculares que contienen elementos del grupo principal de la primera, segunda y tercera fila. La teoría G4 es una modificación mejorada del enfoque anterior de la teoría G3 . Las modificaciones a la teoría G3 son el cambio en una estimación del límite de energía de Hartree-Fock , un conjunto de polarización expandido para el cálculo del conjunto de bases grandes, el uso de energías CCSD (T), el uso de geometrías de la teoría funcional de densidad y el punto cero. energías y dos parámetros de corrección de nivel superior añadidos. Según los desarrolladores, esta teoría ofrece una mejora significativa sobre la teoría G3. Los métodos G4 y G4MP2 relacionados se han ampliado para cubrir los metales de transición. Se ha desarrollado una variante de G4MP2, denominada G4 (MP2) -6X, con el objetivo de mejorar la precisión con componentes de química cuántica esencialmente idénticos. Aplica escalamiento a los componentes energéticos además de utilizar el HLC. En el método G4 (MP2) -XK que está relacionado con G4 (MP2) -6X, los conjuntos de bases de tipo Pople se sustituyen por conjuntos de bases de tipo Karlsruhe personalizados. En comparación con G4 (MP2) -6X, que cubre elementos del grupo principal hasta el criptón, G4 (MP2) -XK se aplica a los elementos del grupo principal hasta el radón.

Enfoque Feller-Peterson-Dixon (FPD)

A diferencia de las "fórmulas químicas modelo" de receta fija, el enfoque FPD consiste en una secuencia flexible de (hasta) 13 componentes que varían con la naturaleza del sistema químico en estudio y la precisión deseada en los resultados finales. En la mayoría de los casos, el componente principal se basa en la teoría de conglomerados acoplados, como CCSD (T), o la teoría de la interacción de la configuración combinada con grandes conjuntos de bases gaussianas (hasta aug-cc-pV8Z, en algunos casos) y la extrapolación al conjunto de bases completo límite. Al igual que con algunos otros enfoques, generalmente se incluyen correcciones aditivas para efectos de núcleo / valencia, relativistas escalares y de correlación de orden superior. Se presta atención a las incertidumbres asociadas con cada uno de los componentes para permitir una estimación aproximada de la incertidumbre en los resultados generales. Los parámetros estructurales precisos y las frecuencias vibratorias son un subproducto natural del método. Si bien las propiedades moleculares calculadas pueden ser muy precisas, la naturaleza computacionalmente intensiva del enfoque FPD limita el tamaño del sistema químico al que se puede aplicar a aproximadamente 10 o menos átomos de la primera / segunda fila.

El enfoque FPD se ha comparado en gran medida con el experimento. Cuando se aplica al nivel más alto posible, FDP es capaz de producir una desviación de la raíz cuadrada media (RMS) con respecto al experimento de 0.30 kcal / mol (311 comparaciones que cubren energías de atomización, potenciales de ionización, afinidades electrónicas y afinidades protónicas). En términos de equilibrio, estructuras del fondo del pozo, FPD da una desviación RMS de 0.0020 Å (114 comparaciones que no involucran hidrógenos) y 0.0034 Å (54 comparaciones que involucran hidrógeno). Se encontró una buena concordancia similar para las frecuencias vibratorias.

T1

El calor de formación T1 calculado (eje y) comparado con el calor de formación experimental (eje x) para un conjunto de> 1800 moléculas orgánicas diversas de la base de datos termoquímica del NIST con errores medios absolutos y RMS de 8.5 y 11.5 kJ / mol, respectivamente. .

El método T1. es un enfoque computacional eficiente desarrollado para calcular calores precisos de formación de moléculas de capa cerrada sin carga que comprenden H, C, N, O, F, Si, P, S, Cl y Br, dentro del error experimental. Es práctico para moléculas de hasta 500 uma de peso molecular.

El método T1 incorporado en Spartan consiste en:

  1. Optimización HF / 6-31G *.
  2. RI-MP2 / 6-311 + G (2d, p) [6-311G *] energía de un solo punto con conjunto de base dual.
  3. Una corrección empírica utilizando recuentos de átomos, órdenes de enlace de Mulliken , energías HF / 6-31G * y RI-MP2 como variables.

Sin embargo, T1 sigue la receta G3 (MP2) al sustituir una geometría HF / 6-31G * por MP2 / 6-31G *, eliminando tanto la frecuencia HF / 6-31G * como QCISD (T) / 6-31G * energía y aproximando la energía MP2 / G3MP2large usando técnicas RI-MP2 de base dual, el método T1 reduce el tiempo de cálculo hasta en 3 órdenes de magnitud. Los recuentos de átomos, los órdenes de enlace de Mulliken y las energías HF / 6-31G * y RI-MP2 se introducen como variables en un ajuste de regresión lineal a un conjunto de 1126 G3 (MP2) calores de formación. El procedimiento T1 reproduce estos valores con errores medios absolutos y RMS de 1,8 y 2,5 kJ / mol, respectivamente. T1 reproduce los calores experimentales de formación para un conjunto de 1805 moléculas orgánicas diversas de la base de datos termoquímica del NIST con errores medios absolutos y RMS de 8,5 y 11,5 kJ / mol, respectivamente.

Enfoque compuesto consistente de correlación (ccCA)

Este enfoque, desarrollado en la Universidad del Norte de Texas por el grupo de investigación de Angela K. Wilson , utiliza los conjuntos de bases consistentes de correlación desarrollados por Dunning y sus colaboradores. A diferencia de los métodos de Gauss-n, ccCA no contiene ningún término ajustado empíricamente. El método funcional de densidad B3LYP con el conjunto de bases cc-pVTZ y cc-pV (T + d) Z para los elementos de la tercera fila (Na - Ar), se utilizan para determinar la geometría de equilibrio. Luego, se utilizan cálculos de un solo punto para encontrar la energía de referencia y contribuciones adicionales a la energía. La energía ccCA total para el grupo principal se calcula mediante:

E ccCA = E MP2 / CBS + ΔE CC + ΔE CV + ΔE SR + ΔE ZPE + ΔE SO

El E energía de referencia MP2 / CBS es la MP2 / ago-cc-pVnZ (donde n = D, T, Q) energías extrapolados en el conjunto completo de base límite por el Peterson mezcla esquema extrapolación exponencial gaussiano. CCSD (T) / cc-pVTZ se utiliza para dar cuenta de la correlación más allá de la teoría MP2:

ΔE CC = E CCSD (T) / cc-pVTZ - E MP2 / cc-pVTZ

Las interacciones núcleo-núcleo y núcleo-valencia se contabilizan mediante MP2 (FC1) / aug-cc-pCVTZ:

ΔE CV = E MP2 (FC1) / aug-cc-pCVTZ - E MP2 / aug-cc-pVTZ

Los efectos relativistas escalares también se tienen en cuenta con un conjunto de bases de Douglass Kroll Hess hamiltoniano y recontratado de una partícula :

ΔE SR = E MP2-DK / cc-pVTZ-DK - E MP2 / cc-pVTZ

Los dos últimos términos son correcciones de energía de punto cero escaladas con un factor de 0,989 para tener en cuenta las deficiencias en la aproximación armónica y las correcciones de órbita de espín consideradas solo para átomos.

El enfoque compuesto coherente de correlación está disponible como palabra clave en NWChem y GAMESS (ccCA-S4 y ccCA-CC (2,3))

Métodos completos de conjuntos de bases (CBS)

Los métodos Complete Basis Set (CBS) son una familia de métodos compuestos, cuyos miembros son: CBS-4M, CBS-QB3 y CBS-APNO, en orden creciente de precisión. Estos métodos ofrecen errores de 2.5, 1.1 y 0.7 kcal / mol cuando se prueban con el equipo de prueba G2. Los métodos CBS fueron desarrollados por George Petersson y colaboradores, y hacen extrapolar varias energías de un solo punto a la energía "exacta". En comparación, los métodos de Gauss-n realizan su aproximación usando correcciones aditivas. Similar al método G2 (+) modificado, CBS-QB3 ha sido modificado mediante la inclusión de funciones difusas en el paso de optimización de geometría para dar CBS-QB3 (+). La familia de métodos CBS está disponible mediante palabras clave en el conjunto de programas Gaussian 09 .

Teorías de Weizmann- n

El Weizmann- n métodos ab initio (W n , n = 1-4) son altamente teorías compuestos precisos desprovistos de parámetros empíricos. Estas teorías son capaces de precisiones inferiores a kJ / mol en la predicción de cantidades termoquímicas fundamentales, como los calores de formación y las energías de atomización, y precisiones sin precedentes en la predicción de constantes espectroscópicas. La capacidad de estas teorías para reproducir con éxito las energías CCSD (T) / CBS (W1 y W2), CCSDT (Q) / CBS (W3) y CCSDTQ5 / CBS (W4) se basa en una combinación juiciosa de conjuntos de bases gaussianas muy grandes con técnicas de extrapolación basadas en el conjunto de bases. Por lo tanto, la alta precisión de las teorías W n viene con el precio de un costo computacional significativo. En la práctica, para sistemas que constan de más de ~ 9 átomos distintos de hidrógeno (con simetría C1), incluso la teoría W1, más económica desde el punto de vista computacional, se vuelve prohibitivamente cara con el hardware de servidor convencional actual.

En un intento por extender la aplicabilidad de los métodos termoquímicos W n ab initio , se han desarrollado versiones explícitamente correlacionadas de estas teorías: W n- F12 ( n = 1-3) y más recientemente incluso una teoría W4-F12. W1-F12 se aplicó con éxito a grandes hidrocarburos (p. Ej., Dodecaedrano , así como a sistemas de relevancia biológica (p. Ej., Bases de ADN ). La teoría W4-F12 se ha aplicado a sistemas tan grandes como el benceno . De manera similar, la W n Los protocolos X que se han desarrollado de forma independiente reducen aún más los requisitos de recursos computacionales mediante el uso de conjuntos de bases más eficientes y, para los componentes menores, métodos de correlación electrónica que son computacionalmente menos exigentes.

Referencias

  1. ^ a b c Ohlinger, William S .; Philip E. Klunzinger; Bernard J. Deppmeier; Warren J. Hehre (enero de 2009). "Cálculo eficiente de calores de formación". El Journal of Physical Chemistry A . Publicaciones ACS. 113 (10): 2165–2175. Código Bibliográfico : 2009JPCA..113.2165O . doi : 10.1021 / jp810144q . PMID  19222177 .
  2. a b A. Karton (2016). "Guía de un químico computacional para la termoquímica precisa de moléculas orgánicas" (PDF) . Revisiones interdisciplinarias de Wiley: ciencia molecular computacional . 6 (3): 292–310. doi : 10.1002 / wcms.1249 .
  3. ^ Mikhail N. Glukhovtsev; Addy Pross; Leo Radom (1996). "Reacciones SN2 de no identidad en fase gaseosa de aniones haluro con haluros de metilo: un estudio computacional de alto nivel". Mermelada. Chem. Soc . 118 (26): 6273–6284. doi : 10.1021 / ja953665n .
  4. ^ Curtiss, Larry A .; Paul C. Redfern; Krishan Raghavachari (2007). "Teoría de Gauss-4" . La Revista de Física Química . 126 (8): 084108. Código Bibliográfico : 2007JChPh.126h4108C . doi : 10.1063 / 1.2436888 . PMID  17343441 .
  5. ^ Mayhall, Nicholas J .; Raghavachari, Krishnan; Redfern, Paul C .; Curtiss, Larry A. (30 de abril de 2009). "Investigación de la teoría de Gaussian4 para la termoquímica de metales de transición". El Journal of Physical Chemistry A . 113 (17): 5170–5175. doi : 10.1021 / jp809179q . ISSN  1089-5639 . PMID  19341257 .
  6. ^ Chan, Bun; Deng, Jia; Radom, Leo (11 de enero de 2011). "G4 (MP2) -6X: una mejora rentable a G4 (MP2)". Revista de teoría química y computación . 7 (1): 112–120. doi : 10.1021 / ct100542x . ISSN  1549-9618 . PMID  26606224 .
  7. ^ Chan, Bun; Karton, Amir; Raghavachari, Krishnan (13 de agosto de 2019). "G4 (MP2) -XK: una variante del método compuesto G4 (MP2) -6X con aplicabilidad ampliada para elementos del grupo principal hasta radón". Revista de Teoría y Computación Química . 15 (8): 4478–4484. doi : 10.1021 / acs.jctc.9b00449 . ISSN  1549-9618 . PMID  31287695 .
  8. ^ a b "Conjunto de bases (química)" , Wikipedia , 3 de diciembre de 2019 , consultado el 6 de diciembre de 2019
  9. ^ David Feller, Kirk A Peterson y David A Dixon (2008). "Un estudio de los factores que contribuyen a las predicciones teóricas precisas de las energías de atomización y estructuras moleculares". Revista de física química . 129 (20): 204105–1–204105–30. Código Bibliográfico : 2008JChPh.129t4105F . doi : 10.1063 / 1.3008061 . PMID  19045850 .
  10. ^ David A Dixon, David Feller y Kirk A Peterson (2012). "Una guía práctica de primeros principios fiables predicciones termoquímicas computacionales a través de la tabla periódica". Informes anuales en Química Computacional Volumen 8 . Informes anuales en química computacional. 8 . págs. 1–28. doi : 10.1016 / B978-0-444-59440-2.00001-6 . ISBN 9780444594402.
  11. ^ David Feller, Kirk A Peterson y David A Dixon (2012). "Otros puntos de referencia de un enfoque basado en clúster acoplado compuesto, convergente, calibrado estadísticamente para estudios termoquímicos y espectroscópicos". Física molecular . 110 (19–20): 2381–2399. Código bibliográfico : 2012MolPh.110.2381F . doi : 10.1080 / 00268976.2012.684897 .
  12. ^ Kirk A Peterson, David Feller y David A Dixon (2012). "Precisión química en termoquímica y espectroscopia ab initio: estrategias actuales y retos de futuro". Cuentas de Química Teórica . 131 : 1079–5. doi : 10.1007 / s00214-011-1079-5 .
  13. ^ David Feller, Kirk A Peterson y Branko Ruscic (2014). "Puntos de referencia de precisión mejorados de moléculas pequeñas utilizando conjuntos de bases consistentes de correlación". Cuentas de Química Teórica . 133 : 1407-16. doi : 10.1007 / s00214-013-1407-z .
  14. ^ a b [1] Libro web de química del NIST
  15. ^ Mulliken, RS (1955). "Análisis de población electrónico en funciones de onda molecular LCAO-MO. I". La Revista de Física Química . 23 (10): 1833–1840. Código Bibliográfico : 1955JChPh..23.1833M . doi : 10.1063 / 1.1740588 .
  16. ^ Deyonker, Nathan J .; Cundari, Thomas R .; Wilson, Angela K. (2006). "El enfoque compuesto consistente de correlación (ccCA): una alternativa a los métodos de Gauss-n" . J. Chem. Phys . 124 (11): 114104. Código bibliográfico : 2006JChPh.124k4104D . doi : 10.1063 / 1.2173988 . PMID  16555871 .
  17. ^ Fabián, Walter MF (2008). "¿Termoquímica precisa a partir de cálculos químicos cuánticos?". Monatshefte für Chemie . 139 (4): 309–318. doi : 10.1007 / s00706-007-0798-8 .
  18. ^ "Enfoque compuesto consistente de correlación (ccCA)" . NWChem.
  19. ^ "Enfoque compuesto consistente de correlación de GAMESS (ccCA)" (PDF) . Gordon Group.
  20. ^ Petersson, G. (2002). "Modelos de conjuntos básicos completos para la reactividad química: desde el átomo de helio a la cinética enzimática". En Cioslowski, J. (ed.). Predicción cuántico-mecánica de datos termoquímicos . Comprensión de la reactividad química. 22 . Springer Holanda. págs. 99–130. doi : 10.1007 / 0-306-47632-0_4 . ISBN 0-7923-7077-5.
  21. ^ Srinivasan Parthiban; Glênisson de Oliveira; Jan ML Martin (2001). "Benchmark ab Initio Energy Profiles for the Gas-Phase SN2 Reactions Y- + CH3X → CH3Y + X- (X, Y = F, Cl, Br). Validation of Hybrid DFT Methods". J. Phys. Chem. Una . 105 (5): 895–904. arXiv : física / 0011029 . Código Bibliográfico : 2001JPCA..105..895P . doi : 10.1021 / jp0031000 .
  22. ^ "Métodos CBS" . Referencia del usuario de Gaussian 09 . Gaussian, Inc.
  23. ^ JML Martin y G. de Oliveira (1999). "Hacia métodos estándar para la termoquímica ab initio de calidad de referencia: teoría W1 y W2". Revista de física química . 111 (5): 1843–1856. arXiv : física / 9904038 . Código Bibliográfico : 1999JChPh.111.1843M . doi : 10.1063 / 1.479454 .
  24. ^ AD Boese; M. Oren; O. Atasoylu; JML Martin; M. Kállay y J. Gauss (2004). "Teoría W3: termoquímica computacional robusta en el rango de precisión kJ / mol". Revista de física química . 120 (9): 4129–4141. arXiv : física / 0311067 . Código bibliográfico : 2004JChPh.120.4129B . doi : 10.1063 / 1.1638736 . PMID  15268579 .
  25. ^ A. Karton; E. Rabinovich; JML Martin; B. Ruscic (2006). "Teoría de W4 para la termoquímica computacional: en la búsqueda de predicciones confiables sub-kJ / mol". Revista de física química . 125 (14): 144108. arXiv : física / 0608123 . Código Bibliográfico : 2006JChPh.125n4108K . doi : 10.1063 / 1.2348881 . PMID  17042580 .
  26. ^ A. Karton, S. Daon y JML Martin (2011). "W4-11: un conjunto de datos de alta confianza para la termoquímica computacional derivado de los datos ab initio de W4" (PDF) . Letras de física química . 510 (4–6): 165. Bibcode : 2011CPL ... 510..165K . doi : 10.1016 / j.cplett.2011.05.007 .
  27. ^ A. Karton y JML Martin (2010). "Rendimiento de la teoría W4 para constantes espectroscópicas y propiedades eléctricas de moléculas pequeñas". Revista de física química . 133 (14): 144102. arXiv : 1008.4163 . Código Bibliográfico : 2010JChPh.133n4102K . doi : 10.1063 / 1.3489113 . PMID  20949982 .
  28. a b A. Karton y JML Martin (2012). "Teoría Wn explícitamente correlacionada: W1-F12 y W2-F12" . Revista de Física Química . 136 (12): 124114. Código bibliográfico : 2012JChPh.136l4114K . doi : 10.1063 / 1.3697678 . PMID  22462842 .
  29. a b N. Sylvetsky; KA Peterson; A. Karton y JML Martin (2016). "Hacia un enfoque W4-F12: ¿Se pueden reconciliar los límites CCSD (T) ab initio explícitamente correlacionados y basados ​​en orbitales?". Revista de Física Química . 144 (21): 214101. arXiv : 1605.03398 . Código bibliográfico : 2016JChPh.144u4101S . doi : 10.1063 / 1.4952410 . PMID  27276939 .
  30. ^ A. Karton; PR Schreiner y JML Martin (2016). "Calores de formación de jaulas de hidrocarburos platónicos mediante procedimientos termoquímicos de alto nivel" (PDF) . Revista de Química Computacional . 37 (1): 49–58. doi : 10.1002 / jcc.23963 . PMID  26096132 .
  31. ^ Chan, Bun; Radom, Leo (13 de noviembre de 2012). "W1X-1 y W1X-2: precisión de calidad W1 con un orden de reducción de magnitud en el costo computacional". Revista de teoría química y computación . 8 (11): 4259–4269. doi : 10.1021 / ct300632p . ISSN  1549-9618 . PMID  26605589 .
  32. ^ Chan, Bun; Radom, Leo (12 de noviembre de 2013). "W3X: un procedimiento compuesto de Post-CCSD (T) rentable". Revista de teoría química y computación . 9 (11): 4769–4778. doi : 10.1021 / ct4005323 . ISSN  1549-9618 . PMID  26583395 .
  33. ^ Chan, Bun; Radom, Leo (12 de mayo de 2015). "W2X y W3X-L: Aproximaciones rentables a W2 y W4 con precisión kJ mol – 1". Revista de teoría química y computación . 11 (5): 2109–2119. doi : 10.1021 / acs.jctc.5b00135 . ISSN  1549-9618 . PMID  26574414 .
  • Cramer, Christopher J. (2002). Fundamentos de la química computacional . Chichester: John Wiley and Sons. págs. 224–228. ISBN 0-471-48552-7.
  • Jensen, Frank (2007). Introducción a la Química Computacional . Chichester, Inglaterra: John Wiley and Sons. págs. 164-169. ISBN 978-0-470-01187-4.