Accurate and efficient simulation of photoemission spectroscopy via Kohn-Sham scattering states

Resumen Técnico: Simulación Precisa y Eficiente de la Espectroscopía de Fotoemisión mediante Estados de Dispersión de Kohn-Sham

Planteamiento del Problema
La espectroscopía de fotoemisión con resolución angular (ARPES) es una herramienta primordial para mapear la estructura electrónica de los sólidos; sin embargo, existe una brecha significativa entre las propiedades intrínsecas del material (estados iniciales) y los espectros observados experimentalmente. La intensidad observada depende no solo de los estados de Bloch iniciales, sino también críticamente de los estados de fotoelectrones finales, los cuales están influenciados por la geometría experimental, la polarización de la luz, la dispersión intra- y entre átomos, y la dispersión múltiple dentro del potencial cristalino. Si bien la teoría de "un solo paso" (one-step) de la fotoemisión proporciona un marco riguroso, las implementaciones existentes —particularmente el maduro método de dispersión múltiple de Korringa-Kohn–Rostoker (KKR)— suelen estar ligadas a formalismos de estructura electrónica especializados. Esto limita su integración con códigos de la teoría del funcional de la densidad (DFT) de ondas planas o de espacio real ampliamente utilizados, dificultando la aplicación de funcionales de intercambio y correlación avanzados y de flujos de trabajo modernos. Además, la precisión de los pseudopotenciales (PP) en la descripción de los estados de fotoelectrones de alta energía, esenciales para cálculos de DFT eficientes, sigue siendo una cuestión abierta.

Metodología
Los autores presentan un marco de primeros principios que calcula los estados de fotoelectrones finales resolviendo directamente la ecuación de Kohn-Sham (KS) con condiciones de contorno de dispersión.

• Formalismo: El enfoque asume un sistema periódico bidimensional y emplea una representación de Laue en el plano para las funciones de onda. La ecuación de KS se transforma de un problema de autovalores (para estados ligados) en una ecuación lineal para estados de dispersión mediante la imposición de condiciones de contorno asintóticas: ondas salientes en $z \to +\infty$ y ondas entrantes/estacionarias en $z \to -\infty$.
• Equivalencia con Lippmann-Schwinger: Los autores demuestran que resolver la ecuación de KS con estas condiciones de contorno es formalmente equivalente a la formulación de Lippmann-Schwinger (LS). Sin embargo, el enfoque directo de KS permite el uso de estrategias de precondicionamiento altamente efectivas derivadas de la diagonalización iterativa, lo que resulta en una convergencia significativamente más rápida en comparación con la resolución directa de la ecuación integral de LS.
• Implementación: El método se implementa dentro de códigos estándar de DFT de ondas planas y de espacio real. Utiliza rejillas discretas en $z$ con una discretización de orden superior (hasta $\mathcal{O}(h^5)$) y resolvedores lineales iterativos.
• Análisis de Pseudopotenciales: Un componente clave del estudio es la evaluación sistemática de los pseudopotenciales. Los autores comparan cálculos de todo el electrón (AE) con pseudopotenciales de norma conservada generalizados (PP), investigando específicamente la necesidad de incluir estados semicore y potenciales no locales para reproducir con precisión las propiedades de dispersión de alta energía.

Contribuciones Clave

1. Marco Eficiente: El artículo presenta un método que es totalmente compatible con los flujos de trabajo estándar de DFT (ondas planas y espacio real), permitiendo el uso de funcionales avanzados y simulaciones a gran escala que anteriormente eran difíciles con códigos especializados de un solo paso.
2. Funciones de Onda Explícitas: A diferencia de los métodos que solo producen intensidades, este enfoque proporciona funciones de onda de fotoelectrones explícitas. Esto permite un análisis transparente de los efectos de los elementos de matriz, la dispersión múltiple y la separación de las propiedades intrínsecas del estado inicial de los efectos extrínsecos del estado final.
3. Validación de Pseudopotenciales: El estudio establece que los pseudopotenciales de alta calidad, particularmente aquellos con partes no locales optimizadas y la inclusión de estados semicore, pueden describir con precisión los estados de fotoemisión en el régimen VUV a XUV, haciendo que el método sea práctico para su uso rutinario.

Resultados
El marco fue validado mediante simulaciones de ARPES de dicroísmo circular (CD-ARPES) para dos sistemas:

• Grafeno Monocapa: Los resultados teóricos para el dicroísmo circular integrado en la energía en la distribución angular (CDAD) se compararon con datos experimentales (Ref. [61]). El método reprodujo con éxito las líneas nodales, las texturas generales de CDAD y la reversión de signo de CDAD a 75 eV. Las comparaciones entre los métodos de todo el electrón y de pseudopotencial mostraron una excelente concordancia hasta 100 eV, siempre que el pseudopotencial incluyera contribuciones no locales. Se encontró que la parte no local es crucial para reproducir las propiedades de dispersión, aunque su efecto varió con la energía de los fotones y los canales de momento angular (por ejemplo, insignificante para ondas parciales de tipo $d$ en el carbono a 52 eV).
• Bulk 2H-WSe$_2$: Se compararon nuevos datos experimentales adquiridos en la instalación CELIA con simulaciones teóricas. Se demostró que la inclusión de estados semicore ($5s, 5p$ para W y $3d$ para Se) en el pseudopotencial es esencial. Las simulaciones con estados semicore lograron una concordancia cuantitativa con el CDAD integrado en valle a través de todos los puntos $K$, mientras que las simulaciones sin estados semicore desviaron significativamente, particularmente en los puntos $K$ 1 y 2. Esto confirmó que las fases de dispersión precisas alrededor de los sitios atómicos son necesarias para capturar la naturaleza interferométrica de la fotoemisión.

Significancia y Reivindicaciones
Los autores afirman que este trabajo establece una ruta robusta y accesible hacia el modelado cuantitativo de ARPES. Al cerrar la brecha entre las plataformas ab initio ampliamente utilizadas y la teoría de la fotoemisión, el método abre la puerta a estudios sistemáticos de texturas orbitales, efectos de muchos cuerpos y fenómenos fuera del equilibrio.

• Versatilidad: El cálculo explícito de los elementos de matriz permite el post-procesamiento de geometías experimentales y polarizaciones de luz arbitrarias.
• Extensibilidad: El método puede combinarse con técnicas de función de Green (por ejemplo, $GW$) para tratar efectos de correlación más allá de la DFT estándar. También se posiciona como una base para simular ARPES resuelto en el tiempo (trARPES) en escenarios de bomba-sonda, ofreciendo una alternativa potencial a las implementaciones de trARPES basadas en KKR, que actualmente están limitadas a la DFT estándar.
• Limitaciones: Los autores señalan modestamente que tratar superficies con celdas unitarias muy grandes puede ser computacionalmente desafiante, y que los sistemas verdaderamente semi-infinitos aún podrían requerir formulaciones de función de Green como KKR. Además, a energías de fotones muy altas, las contribuciones de los estados del núcleo pueden requerir la reconstrucción de las funciones de onda de todo el electrón a partir de las pseudofunciones de onda.

En resumen, el artículo proporciona una alternativa computacionalmente eficiente y formalmente equivalente a la formulación de Lippmann-Schwinger que se integra perfectamente con los códigos de DFT modernos, validada por una concordancia de alta precisión con los datos experimentales de CD-ARPES.