¿Cómo se escala la teoría funcional de densidad con el tamaño del sistema?

11

Teóricamente, ¿cómo escala el tiempo para hacer un cálculo de la teoría funcional de densidad (DFT) con el número de electrones? Estoy interesado en implementaciones DFT "típicas" como VASP, ABINIT, etc., no en códigos O (N).

Max Radin
fuente

Respuestas:

13

La respuesta correcta más simple es que DFT está en . Esto proviene de la idea de que en última instancia está diagonalizando a un hamiltoniano con una dimensión proporcional al número de elecciones y la diagonalización es técnicamente O ( n 3 ) .O(Ne3)O(n3)


O(neMV)MV

El proceso de calcular la acción del hamiltoniano ocurre en un par de pasos:

  • O(nvlnnv)
  • O(nanp)O(nanpnv)
  • O(nanp)O(nanp2)

ne

O(ne2nv)

O(nenvlnnv)

nnvnpnanenvnanenp es un número entero pequeño), pero podría imaginarse aumentando el volumen con un número fijo de electrones (agregando vacío en las geometrías de losas / alambres) o aumentando el número de proyectores con un número fijo de átomos y electrones (utilizando un pseudo potencial más preciso).

O(n2lnn)

Max Hutchinson
fuente
¿Realmente necesita la descomposición propia completa, o solo una pequeña parte del espectro?
Victor Liu
O(ne)ne/2
¡Gracias por tu respuesta! ¿Me puede recomendar algún documento que haya discutido este tema o realizado puntos de referencia?
Max Radin
1
G. Kresse, Computational Materials Science 6, 15 (1996) es la introducción estándar a VASP; es posible que desee comenzar con la sección 6. RM Martin, Estructura electrónica: teoría básica y métodos prácticos (Cambridge Univ Pr, 2004) es una gran introducción a DFT (onda plana y otros), pero probablemente sea menos explícita sobre la complejidad.
Max Hutchinson el