La informática se superponen, cinética, electrón-núcleo de atracción y la repulsión electrón-electrón integrales es posible, incluso para Cartesiano gaussianas de mayor momento angular. Hay fórmulas analíticas que le da la solución de estas integrales. En particular, usted puede encontrarlos en Cocinar el libro [1] (prestar atención al hecho de que el electrón-electrón fórmula está mal!).
Estas fórmulas son muy sencillos de implementar (aparte de en el cálculo de los Chicos de la función, que puede ser escrito como una función Gamma incompleta, por lo que usted puede retransmitir en alguna biblioteca), pero son muy ineficientes. Sólo la cuenta de cuantas nasted sumatorias que tiene (que se traduce en nasted bucles en una aplicación sencilla). Este es el principal problema!
Esta es la razón por la alternativa (recursivo) se han desarrollado métodos. En Helgaker del libro [2] se puede encontrar una buena descripción de estos métodos alternativos
- Obara-Saika esquema de
- McMurchie-Davidson esquema de
- Rys cuadratura
Si sólo quieres escribir un trabajo de Hartree-Fock programa puede usar la solución analítica. Sin embargo, si usted necesita eficiencia que usted tendrá que mirar a estos métodos.
Nota: personalmente me encanta Python, pero es realmente unappropriate para computacionalmente exigente cálculos! Mi Hartree-Fock programa en Python que lleva al menos 2 minutos para calcular dos electrones integrales en el STO-3G de la base establecida para $\ce{H_2O}$, mientras que un Fortran versión del mismo programa tiene un montón de segundos.
Para la fórmula correcta para dos electrones integrales usted puede mirar en esta discusión.
[1] D. Cook, Manual de la Química Computacional, Oxford University Press, 1998.
[2] T. Helgaker, P. Jørgensen y J. Olsen, Molecular Electrónica-Teoría de la Estructura, Wiley, 2000.
