Flat-band projected versus fully atomistic twisted bilayer graphene

Resumen Técnico: Proyección de banda plana frente a grafeno de bicapa retorcida totalmente atomístico

Planteamiento del problema
El grafeno de bicapa retorcida de ángulo mágico (MATBG) sirve como una plataforma primaria para el estudio de la topología y las correlaciones electrónicas fuertes. Sin embargo, existe una tensión entre el costo computacional de los modelos totalmente atomísticos (que requieren miles de átomos por celda unidad de moiré) y la necesidad de capturar con precisión los efectos de muchos cuerpos. Un trabajo previo introdujo un método de proyección para construir teorías efectivas de baja energía integrando las bandas remotas, de forma análoga a la renormalización wilsoniana. El problema central abordado en este trabajo es evaluar la precisión de este enfoque de "proyección de banda plana" frente al modelo "totalmente atomístico" (tight-binding completo) a través de varias fases de ruptura de simetría en la neutralidad de carga. Específicamente, los autores buscan determinar si el esquema de proyección, que asume que las bandas remotas están efectivamente congeladas, puede reproducir las estructuras de bandas, las energías totales y las características de ruptura de simetría de la teoría completa sin ambigüedades de doble conteo.

Metodología
Los autores realizan cálculos de Hartree-Fock auto-consistentes en el punto de neutralidad de carga para un sistema MATBG con un ángulo de giro de $1.05^\circ$ (que contiene 11,908 átomos en la celda unidad de moiré).

1. Modelos: Se comparan dos modelos distintos:
   • Totalmente Atomístico: Un modelo de tight-binding definido por una parametrización de Slater-Koster, que incluye relajación de red y un potencial de Coulomb de doble compuerta regularizado por una energía de Hubbard en el sitio ($U=4$ eV).
   • Proyección de Banda Plana: Una teoría efectiva construida mediante la proyección del estado normal renormalizado sobre las bandas planas de baja energía, siguiendo el procedimiento establecido en la Ref. [1].
2. Parámetros: La constante dieléctrica se establece en $\epsilon_r = 10$, y la distancia de la compuerta es $\xi = 10$ nm.
3. Ruptura de Simetría: El estudio investiga varias fases correlacionadas, incluyendo el semimetal nemático (NSM), el efecto Hall anómalo cuántico (QAH), la coherencia de intervalo de Kramers (KIVC), la polarización orbital (OP), la coherencia de intervalo de inversión temporal (TIVC) y la polarización de valle (VP).
4. Herramientas de Análisis: Los autores emplean un nuevo conjunto de parámetros de orden locales derivados de los solapamientos de funciones de onda en el espacio real. Estos generalizan el operador de valle microscópico para visualizar las funciones de onda localmente y cuantificar la ruptura de simetría, complementando los análisis en el espacio de momentos.

Resultados Clave

1. Acuerdo en la Estructura de Bandas y Energía: Las soluciones de la proyección de banda plana muestran una excelente concordancia con los resultados totalmente atomísticos.
   • Estructuras de Bandas: Las estructuras de bandas de los estados de ruptura de simetría difieren por solo unos pocos meV. Las fases totalmente atomísticas exhiben brechas ligeramente mayores en el punto $K$ (oscilando entre 3 meV en el estado VP y 6 meV en el estado TIVC) en comparación con el modelo proyectado.
   • Energías Totales: Los cálculos totalmente atomísticos producen energías totales más bajas que los cálculos de banda plana por aproximadamente 3 meV en todas las fases. Las energías del estado normal simétrico son idénticas por construcción, con discrepancias menores atribuidas a las mallas de momento más finas utilizadas en los cálculos proyectados.
2. Justificación de la Proyección: La precisión de la proyección se justifica por la separación de escalas de energía. En las fases de ruptura de simetría, el desdoblamiento de la banda plana ($\pm 15\text{--}20$ meV) es menor que la mitad de la energía de los bordes de las bandas remotas en el estado normal ($\approx +60$ meV y $-50$ meV). Esto indica que las bandas remotas están efectivamente congeladas, a diferencia del estado no interactuante donde las bandas remotas aparecen a energías comparables al desdoblamiento de la banda plana.
3. Parámetros de Orden Locales: Los autores logran visualizar con éxito las funciones de onda localmente en el espacio real.
   • Los parámetros de orden se concentran en el centro de la celda unidad (región AA), lo cual es consistente con el peso espectral de los orbitales activos.
   • Se identifican simetrías específicas para cada fase: KIVC rompe la inversión temporal ($T$) y el $U(1)_v$ de valle pero preserva $C_{2z}T$; QAH rompe las simetrías de espejo y $T$; OP y NSM rompen $C_{2z}$ y $C_{3z}$ respectivamente; VP rompe $C_{2z}$ y $T$; y TIVC rompe $U(1)_v$ mientras preserva todas las simetrías de punto y $T$.
   • La magnitud de los parámetros de orden integrados en las bandas planas es comparable al máximo teórico (cercano a 1), lo que indica una fuerte ruptura de simetría en las bandas planas. Incluir las bandas remotas no altera significativamente estos valores numéricos.

Significancia y Reivindicaciones
El artículo afirma validar la fiabilidad del esquema de proyección de banda plana para MATBG. Al demostrar que el modelo proyectado reproduce los resultados de la teoría de tight-binding completa con alta fidelidad (diferencias de solo unos pocos meV), los autores confirman que la proyección de muchos cuerpos captura correctamente los efectos de interacción de las bandas integradas, incluyendo la renormalización del ancho de banda de la banda plana.

Este trabajo también introduce un conjunto de herramientas prácticas para caracterizar sistemas de MATBG y sistemas de panal genéricos. Los novedosos parámetros de orden locales permiten el cálculo directo de los solapamientos de las funciones de onda en modelos atomísticos sin necesidad de construir la base de Chern de las descripciones en el espacio de momentos. Esto proporciona una perspectiva complementaria a las matrices de densidad en el espacio de momentos y los análisis de base de fermiones pesados, permitiendo la visualización de la ruptura de simetría directamente en el espacio real. Los autores señalan que, si bien el tratamiento preciso de la función dieléctrica (por ejemplo, mediante la teoría GW) y el acoplamiento electrón-fonón están fuera del alcance de este benchmark específico, los resultados actuales establecen la robustez del método de proyección para estudiar fases correlacionadas en MATBG.