A Full Minimal Coupling GW-BSE Framework for Circular Dichroism in Solids: Applications to Chiral 2D Perovskites

Resumen Técnico: Un Marco de Acoplamiento Mínimo Completo GW-BSE para el Dicroísmo Circular en Sólidos

Planteamiento del Problema
El dicroísmo circular (CD) y las respuestas quiropticas sirven como sondas críticas para la quiralidad y la geometría en el espacio de momentos en los sólidos. Sin embargo, los cálculos de primeros principios para el CD en sistemas periódicos con efectos excitónicos fuertes, como los híbridos orgánico-inorgánicos de perovskita 2D (2D-HOIPs), siguen siendo teóricamente desafiantes. Los enfoques tradicionales dependen de expansiones multipolares y técnicas de suma sobre estados (SOS). Estos métodos sufren de limitaciones significativas:

1. Ambigüedades de Calibre y Origen: Las evaluaciones SOS truncadas suelen arrojar resultados que dependen de forma espuria de la elección del origen de la celda unidad o del calibre del potencial vectorial, a menos que se apliquen correcciones explícitas.
2. Problemas de Convergencia: Los enfoques SOS convergen lentamente con el número de estados, particularmente en sistemas con colectores de bandas densos.
3. Inestabilidad Numérica: Las formulaciones estándar de calibre de velocidad exhiben divergencias para transiciones intrabanda o transiciones interbanda entre bandas casi degeneradas.
4. Complejidad Computacional: Las formulaciones alternativas de calibre de longitud (por ejemplo, de Blount) requieren derivadas numéricas de las funciones de onda ($\nabla_k$), lo que requiere mallas $k$ densas y calibres suaves, lo cual resulta computacionalmente prohibitivo para los cálculos de la ecuación de Bethe-Salpeter (BSE) que escalan como $O(N^5)$.
5. Física Incompleta: La teoría del funcional de la densidad (DFT) convencional no logra capturar las correlaciones electrón-hueco de muchos cuerpos, esenciales para las excitaciones ópticas precisas en estos materiales.

Metodología
Los autores introducen un marco de Acoplamiento Mínimo Completo (FMC) para calcular el CD en sólidos, integrando la aproximación $GW$ para las energías de cuasipartícula con la ecuación de Bethe-Salpeter (BSE) para las interacciones electrón-hueco.

• Formulación Teórica: A diferencia de los métodos SOS que expanden el Hamiltoniano de interacción luz-materia en potencias del vector de onda del fotón $q$ (reteniendo los términos de Dipolo Eléctrico (ED), Dipolo Magnético (MD) y Cuadrupolo Eléctrico (EQ) por separado), el método FMC retiene la interacción completa de acoplamiento mínimo $p \to p + \frac{e}{c}A$ sin expandir en $q$.
• Elementos de Matriz: El método evalúa los elementos de matriz ópticos directamente en un vector de onda de fotón finito $q$. El elemento de matriz excitónica se calcula como $\langle 0| e^{iq\cdot r}p_x |S\rangle = \sum_{vck} A^S_{vck,q} \langle vk+q| e^{iq\cdot r}p_x |ck\rangle$.
• Ventajas Clave:
  • Invariancia de Calibre: La formulación es manifiestamente invariante de calibre, eliminando las ambigüedades de origen.
  • Transiciones Intrabanda y Casi-Degeneradas: Al tratar la perturbación en primer orden en $q$ (incluyendo la dispersión espacial), el FMC incluye naturalmente las contribuciones intrabanda y maneja las transiciones entre bandas casi degeneradas sin las divergencias de los denominadores de energía que afectan a la SOS.
  • Igualdad de Condiciones: Los canales ED, MD y EQ se tratan en igualdad de condiciones sin requerir cálculos separados y computacionalmente intensivos de momento angular orbital o derivadas cuadrupolares.
• Aproximaciones: El método aproxima el momento canónico como $-i\hbar\nabla$ para evitar complicaciones numéricas de desplazamientos de momento duales, omitiendo así las contribidades de pseudopotenciales no locales (un efecto explorado en la Información Suplementaria).
• Implementación: El marco se aplica a dos 2D-HOIPs protípicos: (S-NEA)$_2$PbBr$_4$ (S-NPB) y (S-MBA)$_2$PbI$_4$ (S-MPI). Los estados fundamentales se calculan mediante DFT (GGA), las energías de cuasipartícula mediante $G_0W_0$ y los efectos excitónicos mediante BSE en la aproximación de Tamm-Dancoff.

Resultos Clave

• Papel de los Canales MD y EQ: Los cálculos en S-NPB revelan que los canales de Dipolo Magnético (MD) y Cuadrupolo Eléctrico (EQ) contribuyen con la misma magnitud pero con signos opuestos a la señal de CD. Su cancelación parcial resulta en una señal de CD significativamente menor que la predicha por modelos de solo MD.
• Efectos Excitónicos y Forma de Línea: La inclusión de excitones y el término EQ en S-NPB conduce a una estructura de pico asimétrica en el espectro de CD, contrastando con el "efecto Cotton" simétrico predicho por modelos que niegan el EQ o los estados moleculares. La convergencia del espectro de CD en SOS requiere más de 1,000 bandas para capturar la asimetría introducida por los estados moleculares, mientras que el FMC converge de forma más robusta.
• Estabilidad Numérica en Sistemas Degenerados: Para S-MPI, que presenta una estructura de bandas con estados casi degenerados debido al débil acoplamiento entre capas, el enfoque SOS falla completamente, produciendo espectros oscilantes e inestables que dependen de los umbrales de energía. En contraste, el enfoque FMC produce un espectro cualitativamente fiable y estable.
• Contribuciones Intrabanda: El FMC captura explícitamente elementos de matriz de transición diagonal (intrabanda) considerables, que corresponden a términos similares a la conexión de Berry. Estos términos son omitidos en las formulaciones SOS estándar pero se demuestra que son significativos para la respuesta de CD.
• Momento Angular Orbital: El FMC proporciona elementos diagonales suaves y físicamente consistentes para las matrices de momento angular orbital, mientras que los resultados de SOS exhiben fluctuaciones inestables dependientes de parámetros numéricos.

Significancia y Reivindicaciones
El artículo afirma que el marco FMC ofrece una alternativa computacionalmente eficiente y numéricamente robusta a los métodos convencionales de SOS y de expansión multipolar para predecir fenómenos optoelectrónicos quirales en sólidos complejos. Al eliminar los problemas de convergencia lenta, las ambigüedades de calibre y las divergencias asociadas con las degeneraciones de banda, el FMC permite predicciones de primeros principios precisas de CD en materiales con fuertes efectos excitónicos y estructuras de bandas densas. El estudio establece que las transiciones intrabanda y casi-degeneradas, capturadas naturalmente por el FMC, pueden remodelar drásticamente los espectros de CD, resaltando la importancia de la geometría cuántica (términos similares a la conexión de Berry) en la física quiroptica. Este trabajo proporciona una base para comprender el origen del CD en los HOIPs quirales y apoya el desarrollo de dispositivos optoelectrónicos quirales.