Cell Natural Orbitals in Quantum Materials

Resumen Técnico: Orbitales Naturales de Celda en Materiales Cuánticos

Planteamiento del Problema
El modelado preciso de la materia cuántica correlacionada requiere una descripción mínima de los grados de libertad de una sola partícula de baja energía que sean relevantes. Mientras que las descripciones en el espacio de momentos (funciones de onda de Bloch) capturan la dispersión y la simetría, los modelos en el espacio real son esenciales para tratar las interacciones electrónicas locales. Sin embargo, identificar un conjunto mínimo de orbitales locales que reproduzca fielmente la dispersión, la simetría y la geometría cuántica (específicamente la densidad de carga intracelular y los elementos de matriz de dipolo) de la banda objetivo es una tarea no trivial. Los métodos existentes, como las Funciones de Wannier Máximamente Localizadas (MLWF) o las técnicas de desentrelazamiento, a menudo requieren preespecificar el número de orbitales o dependen de proyecciones iniciales motivadas físicamente, fallando al optimizar el espacio de Hilbert local mínimo dictado por las propias bandas. Además, las obstrucciones topológicas (por ejemplo, en bandas de Chern) o los límites atómicos obstruidos pueden impedir la existencia de una base simétrica y localizada que abarque el espacio de la variedad objetivo utilizando enfoques estándar.

Metodología
Los autores introducen una construcción sistemática de orbitales locales, denominada Orbitales Naturales de Celda (CNO, por sus siglas en inglés), derivados únicamente del proyector sobre un conjunto de bandas de Bloch objetivo. El objeto matemático central es la Matriz de Densidad Reducida de Una Partícula de la Celda Unidad (UC-1pRDM), denotada como $L$.

1. Construcción de la UC-1pRDM: Sea $P$ el proyector sobre las bandas de Bloch objetivo y $\Pi$ el proyector que restringe el espacio de Hilbert a una sola celda unidad de referencia. La UC-1pRDM se define como $L = \Pi P \Pi$. En la base de orbitales, los elementos de matriz de $L$ están dados por $L_{\alpha\beta} = \langle c^\dagger_\beta c_\alpha \rangle$, calculados mediante la suma sobre las bandas objetivo e incluyendo una matriz de incrustación $V_k$ que da cuenta de la posición real de los orbitales dentro de la celda unidad.
2. Diagonalización: Los autoestados de $L$, denotados como $|\tau_a\rangle$, son los CNO. Sus correspondientes autovalores, $\lambda_a$, representan la ocupación de cada CNO dentro de las bandas objetivo.
3. Interpretación Física:
   • Peso Espectral: $\lambda_a$ cuantifica la importancia de cada orbital local. Un $\lambda_1 \approx 1$ implica que la banda está bien descrita por un solo orbital; las ocupaciones fraccionarias indican un carácter multi-orbital.
   • Entrelazamiento: Los autovalores $\lambda_a$ son los valores singulares del mapa entre los orbitales locales y el subespacio de la banda. Definen el espectro de entrelazamiento a través de la partición de la celda unidad, donde la entropía de von Neumann $S_{cell}$ mide el fallo de la subvariedad ocupada al factorizarse a través de la frontera de la celda.
   • Simetría y Topología: Al elegir una partición de la celda unidad centrada en una posición Wyckoff específica $q$, los CNO se transforman como representaciones irreducibles (irreps) del grupo de simetría del sitio $G_q$. Esto permite que los CNO identifiquen candidatos a Representaciones de Bandas Elementales (EBR).
   • Ceros: Las funciones de envolvente $s_{a,k}$ (coeficientes de Bloch reescalados) revelan ceros en la zona de Brillouin. Las obstrucciones topológicas (por ejemplo, números de Chern no nulos) imponen ceros irremovibles, mientras que los límites atómicos obstruidos o las elecciones de Wyckoff incompatibles conducen a ceros forzados por la simetría que pueden eliminarse desplazando el centro de la partición.

Contribuciones Clave y Resultados

• Selección Sistemática de la Base: El artículo proporciona un método para determinar la base de orbitales locales óptima sin supuestos previos sobre el número de orbitales o su carácter atómico específico. Los CNO están determinados de forma única por las funciones de onda de Bloch y la elección de la partición en el espacio real.
• Aplicación a Aislantes Atómicos Obstruidos: Utilizando un modelo de monocapa de 1H-TMD, los autores demuestran que elegir una partición de la celda unidad centrada en el sitio atómico (posición Wyckoff 1a) produce ocupaciones fraccionarias y ceros forzados por la simetría en las funciones de envolvente, reflejando la naturaleza obstruida de la banda. Desplazar la partición al centro de Wannier (posición Wyckoff 1c) resulta en un CNO dominante con $\lambda_1 \approx 1$ y sin ceros, identificando correctamente la descripción local mínima.
• Aplicación a Bandas Topológicas: En el modelo Bernevig-Hughes-Zhang (BHZ) sin espín, las ocupaciones de los CNO varían suavemente a través de la transición de fase topológica, pero las funciones de envolvente exhiben ceros irremovibles en la fase topológica ($C=1$), manifestando la obstrucción de Wannier en el espacio real.
• Construcción de Modelos Locales para tWSe$_2$ de Bicapa Retorcida: Los autores aplican el método al WSe$_2$ de bicapa retorcida (tWSe$_2$) a través de un rango de ángulos de rotación ($1^\circ$ a $5^\circ$).
  • Identifican que las dos bandas superiores de moiré requieren descripciones locales diferentes dependiendo del ángulo de rotación y de la elección de la partición (sitios de apilamiento MM, MX o XM).
  • En el régimen topológico ($\theta > 1.51^\circ$), ningún orbital local único abarca las bandas debido al número de Chern total finito.
  • Al seleccionar subconjuntos de CNO (por ejemplo, un conjunto multicéntrico que combina orbitales MM, MX y XM) que maximizan el peso espectral total y satisfacen la condición de completitud (determinante no nulo de la matriz de solapamiento $S_k$), construyen modelos de red de Tight-Binding.
  • Wannierización de un solo paso (Single-Shot): Estos subconjuntos de CNO sirven como estados de prueba para un procedimiento de Wannierización de un solo paso. Los modelos de red de Tight-Binding resultantes reproducen no solo la dispersión de las bandas sino también la geometría cuántica (curvatura de Berry) del modelo de continuo.
  • Compromisos (Trade-offs): Los autores comparan un modelo centrado en MM (3 orbitales en un sitio) con un modelo multicéntrico (orbitales en diferentes sitios). El modelo multicéntrico produce funciones de Wannier más localizadas (hoppings de corto alcance) pero requiere un clúster más grande para tratamientos de muchos cuerpos que preserven la simetría (como el DMFT de clúster), mientras que el modelo centrado en MM preserva la simetría $C_3$ más fácilmente pero tiene interacciones de más largo alcance.

Significancia y Reivindicaciones
El artículo afirma que los CNO ofrecen una medida directa, derivada de la función de onda, del carácter multi-orbital de las bandas objetivo. Al revertir el enfoque convencional de inferir orbitales a partir de la química o los autovalores de simetría, los CNO permiten que los grados de libertad locales se infieran directamente desde la propia variedad de bandas.

La significancia radica en:

1. Optimalidad: Los CNO proporcionan los estados de prueba óptimos para la Wannierización, minimizando el número de orbitales requeridos para abarcar las bandas objetivo mientras se preserva la geometría cuántica.
2. Poder Diagnóstico: El espectro de ocupación y los ceros de la función de envolvente sirven como diagnósticos para las obstrucciones topológicas y los límites atómicos obstruidos, distinguiendo entre ceros removibles (simetría) e irremovibles (topológicos).
3. Base para Modelos Interactuantes: Dado que los CNO diagonalizan la matriz de correlación de la celda unidad, proporcionan una base natural para las interacciones proyectadas sobre la banda. La jerarquía de ocupación $\lambda_a$ ordena las escalas de interacción, ofreciendo una ruta para construir modelos de Hubbard multi-orbitales y definir espacios de impureza para métodos de incrustación cuántica (como DMFT) sin elecciones arbitrarias.
4. Generalidad: El método se aplica tanto a modelos de red como de continuo (mediante una rejilla en el espacio real derivada de la base de ondas planas) y es compatible con cálculos de primeros principios, ya que solo requiere el proyector de la banda como entrada.

Los autores concluyen que esta perspectiva tiende un puente entre las funciones de onda de primeros principios y los modelos de continuo con los grados de libertad locales mínimos necesarios para las teorías efectivas de materiales cuánticos topológicos y correlacionados.