Cell Natural Orbitals in Interacting Topological Bands

Resumen Técnico: Orbitales Naturales de Celda en Bandas Topológicas con Interacciones

Planteamiento del Problema
Un desafío fundamental en la física de la materia condensada es reconciliar la descripción en el espacio de momentos de las bandas topológicas con el lenguaje del espacio real requerido para describir las interacciones electrón-electrón. La teoría de bandas describe naturalmente las fases topológicas, pero la construcción de funciones de Wannier simétricas y exponencialmente localizadas se ve obstruida por las propiedades globales de las funciones de onda de Bloch en sistemas topológicos. Esta obstrucción complica la representación de las interacciones proyectadas, que suelen modelarse mediante orbitales locales de rango finito. En las heteroestructuras de moiré, donde las bandas planas portan una geometría cuántica no trivial y albergan fases fuertemente correlacionadas, esta tensión es aguda. Los enfoques existentes, como la analogía de los niveles de Landau o los modelos fenomenológicos de fermiones pesados, a menudo dependen de suposiciones basadas en ansatz para los orbitales localizados o no logran proporcionar un marco sistemático para clasificar la importancia de las diferentes contribuciones orbitales. Existe una falta de un método para generar orbitales locales directamente desde las bandas de baja energía que cuantifique simultáneamente la complejidad orbital mínima requerida y establezca una jerarquía de fuerzas de interacción.

Metodología
Los autores introducen un marco basado en Orbitales Naturales de Celda (CNO, por sus siglas en inglés), definidos como los autoestados de la matriz de densidad de una partícula reducida de la celda unidad (uc-RDM). La metodología procede de la siguiente manera:

1. Definición de los CNO: Para un conjunto de bandas seleccionadas con un proyector $P$, los autores definen una matriz $L$ que representa el promedio del proyector de la banda restringido a una sola celda unidad. Al diagonalizar esta matriz de $N_\alpha \times N_\alpha$ (donde $N_\alpha$ es el número de grados de libertad internos por celda), se obtienen los autovalores $\lambda_a$ (ocupaciones naturales) y los autovectores $|\tau_a\rangle$ (los CNO).
2. Descomposición en Valores Singulares (SVD): Los autores establecen una conexión entre la uc-RDM y los factores de forma de la densidad. Al definir un operador $U$ que mapea el subespacio de baja energía al espacio de Hilbert local de la celda unidad, demuestran que la SVD de $U$ proporciona una descomposición del núcleo de solapamiento $\Lambda$. Los autovalores $\lambda_a$ corresponden a los valores singulares al cuadrado, y los CNO sirven como los vectores singulares izquierdos. Los vectores singulares derechos son las funciones envolventes de espacio de momentos $s_{a,n,k}$.
3. Jerarquía y Truncamiento: El espectro de ocupación $\{\lambda_a\}$ proporciona una medida del entrelazamiento de la celda unidad y la complejidad orbital. Los autores proponen un esquema de truncamiento basado en la geometría, donde el espacio de Hilbert se reduce reteniendo solo los $N_\tau$ CNO principales. Esto minimiza el error de Hilbert-Schmidt en el núcleo de solapamiento.
4. Interacciones Proyectadas: El formalismo se aplica a las interacciones de densidad-densidad proyectadas. Los autores sustituyen la expansión truncada de los CNO en el Hamiltoniano de interacción para derivar elementos de matriz de interacción efectivos en la base de orbitales locales.

Resultados Clave

• Jerarquía Intrínseca de Interacciones: La descomposición CNO revela una jerarquía sistemática en las interacciones proyectadas de la banda. El CNO dominante (mayor $\lambda_a$) captura la contribución primaria a la densidad de carga local, mientras que los canales subdominantes portan elementos de matriz de interacción progresivamente más débiles. Esto permite una aproximación controlada donde los canales dominantes pueden requerir un tratamiento dinámico (por ejemplo, autoenergía), mientras que los canales subdominantes pueden tratarse a un nivel de campo medio estático.
• Ceros Topológicamente Protegidos: Para bandas topológicas (por ejemplo, bandas de Chern), las funciones envolventes $s_{a,n,k}$ de los CNO principales exhiben ceros topológicamente protegidos en la zona de Brillouin. Un solo CNO no puede representar el estado de Bloch en todo el espacio debido a estos ceros. En consecuencia, un truncamiento de un solo orbital conduce a singularidades (colas de ley de potencia) en las interacciones del espacio real. La inclusión de CNO subdominantes regulariza la matriz de Gram, eliminando estas singularidades y restaurando el decaimiento exponencial en las interacciones del espacio real.
• Aplicación al Grafeno de Bicapa Retorcido de Ángulo Mágico (TBG): Aplicado al TBG quiral, se encuentra que el CNO dominante está centrado en el sitio AA, asemejándose al fermión-$f$ del modelo de fermiones pesados. Los CNO subdominantes son más deslocalizados. La jerarquía explica por qué la interacción en el sitio AA es dominante.
• Correlaciones Dependientes del Momento: La presencia de ceros topológicos en la envolvente del CNO principal conduce a escalas de correlación dependientes del momento. En los momentos donde la envolvente principal se anula (por ejemplo, el punto $\Gamma$ en el TBG quiral), el canal de interacción dominante se desacopla y las correlaciones están gobernadas por los canales subdominantes. Esto resulta en una función espectral dependiente del momento donde las brechas pueden cerrarse en puntos específicos de alta simetría.
• Medidas de Complejidad: Los autores introducen el Ratio de Participación Inversa (IPR) y la entropía de Shannon normalizada del espectro de ocupación como medidas cuantitativas de la complejidad orbital. Estas métricas distinguen entre bandas que pueden describirse con un solo orbital y aquellas que requieren una descripción multiorbital debido a una compleja geometría cuántica.

Significancia y Reivindicaciones
El artículo afirma establecer los CNO como un "puente informado por la geometría" entre la topología de banda y las correlaciones en el espacio real. La significancia de este trabajo radica en:

1. Marco Sistemático: Proporciona un método sistemático y libre de ansatz para generar orbitales locales a partir de bandas de baja energía, determinando el conjunto mínimo de orbitales necesarios para representar fielmente los factores de forma y las interacciones.
2. Resolución de la Obstrucción Topológica: Demuestra cómo la obstrucción topológica a las funciones de Wannier localizadas se manifiesta como una jerarquía de canales CNO. La obstrucción no es una barrera para la representación local, sino que dicta la necesidad de múltiples orbitales para capturar con precisión la física de corto alcance.
3. Aproximaciones Controladas: La jerarquía de las ocupaciones CNO ofrece una forma fundamentada de truncar el espacio de Hilbert para materiales topológicos interactuantes, permitiendo aproximaciones controladas para métodos como la Teoría de Campo Medio Dinámico (DMFT) o esquemas de incrustación cuántica.
4. Perspectiva Física: El formalismo aclara cómo las variaciones de la densidad de carga dentro de una celda unidad generan la dependencia del momento en las envolventes de interacción y en las funciones espectrales de partícula única, vinculando la geometría cuántica directamente con los efectos de correlación observables.

Los autores enfatizan que, si bien el marco ofrece utilidades conceptuales y numéricas, no pretende resolver todos los problemas de muchos cuerpos en materiales topológicos, sino proporcionar un punto de partida robusto para analizar interacciones y fases emergentes respetando la geometría subyacente de las funciones de onda.