Cell Natural Orbitals in Quantum Materials

Resumo Técnico: Orbitais Naturais de Célula em Materiais Quânticos

Definição do Problema
A modelagem precisa da matéria quântica correlacionada requer uma descrição mínima dos graus de liberdade de partícula única relevantes de baixa energia. Embora as descrições no espaço de momento (funções de onda de Bloch) capturem a dispersão e a simetria, modelos no espaço real são essenciais para tratar interações eletrônicas locais. No entanto, identificar um conjunto mínimo de orbitais locais que reproduza fielmente a dispersão, a simetria e a geometria quântica (especificamente a densidade de carga intracelular e os elementos de matriz de dipolo) da banda alvo é uma tarefa não trivial. Métodos existentes, como Funções de Wannier Maximamente Localizadas (MLWFs) ou técnicas de desenredamento (disentanglement), frequentemente exigem a especificação prévia do número de orbitais ou dependem de projeções iniciais motivadas fisicamente, falhando em otimizar o espaço de Hilbert local mínimo ditado pelas próprias bandas. Além disso, obstruções topológicas (por exemplo, em bandas de Chern) ou limites atômicos obstruídos podem impedir a existência de uma base localizada e simétrica que abranja o manifold alvo usando abordagens padrão.

Metodologia
Os autores introduzem uma construção sistemática de orbitais locais, denominada Orbitais Naturais de Célula (CNOs), derivados exclusivamente do projetor sobre um conjunto de bandas de Bloch alvo. O objeto matemático central é a Matriz de Densidade Reduzida de Uma Partícula da Unidade de Célula (UC-1pRDM), denotada por $L$.

1. Construção da UC-1pRDM: Seja $P$ o projetor sobre as bandas de Bloch alvo e $\Pi$ o projetor que restringe o espaço de Hilbert a uma única célula de referência. A UC-1pRDM é definida como $L = \Pi P \Pi$. Na base orbital, os elementos da matriz são dados por $L_{\alpha\beta} = \langle c^\dagger_\beta c_\alpha \rangle$, calculados somando-se as bandas alvo e incluindo uma matriz de embutimento $V_k$ que contabiliza a posição real dos orbitais dentro da célula unitária.
2. Diagonalização: Os autovetores de $L$, denominados $|\tau_a\rangle$, são os CNOs. Seus respectivos autovalores, $\lambda_a$, representam a ocupação de cada CNO dentro das bandas alvo.
3. Interpretação Física:
   • Peso Espectral: $\lambda_a$ quantifica a importância de cada orbital local. Um $\lambda_1 \approx 1$ implica que a banda é bem descrita por um único orbital; ocupações fracionárias indicam um caráter multi-orbital.
   • Emaranhamento: Os autovalores $\lambda_a$ são os valores singulares do mapeamento entre orbitais locais e o subspace da banda. Eles definem o espectro de emaranhamento através da partição da célula unitária, onde a entropia de von Neumann $S_{cell}$ mede a falha do subspace ocupado em fatorar através da fronteira da célula.
   • Simetria e Topologia: Ao escolher uma partição de célula unitária centrada em uma posição de Wyckoff específica $q$, os CNOs transformam-se como representações irredutíveis (irreps) do grupo de simetria de sítio $G_q$. Isso permite que os CNOs identifiquem candidatos a Representações de Bandas Elementares (EBRs).
   • Zeros: As funções de envelope $s_{a,k}$ (coeficientes de Bloch reescalonados) revelam zeros na zona de Brillouin. Obstruções topológicas (por exemplo, números de Chern não nulos) impõem zeros irremovíveis, enquanto limites atômicos obstruídos ou escolhas de Wyckoff incompatíveis levam a zeros forçados pela simetria, que podem ser removidos deslocando o centro da partição.

Contribuições Principais e Resultados

• Seleção Sistemática de Base: O artigo fornece um método para determinar a base de orbitais locais ideal sem pressuposições sobre o número de orbitais ou seu caráter atômico específico. Os CNOs são determinados unicamente pelas funções de onda de Bloch e pela escolha da partição no espaço real.
• Aplicação a Isolantes Atômicos Obstruídos: Usando um modelo de monocamada de 1H-TMD, os autores demonstram que escolher uma partição de célula unitária centrada no sítio atômico (posição de Wyckoff 1a) produz ocupações fracionárias e zeros forçados pela simetria nas funções de envelope, refletindo a natureza obstruída da banda. Deslocar a partição para o centro de Wannier (posição de Wyckoff 1c) resulta em um CNO dominante com $\lambda_1 \approx 1$ e sem zeros, identificando corretamente a descrição local mínima.
• Aplicação a Bandas Topológicas: No modelo Bernevig-Hughes-Zhang (BHZ) sem spin, as ocupações dos CNOs variam suavemente através da transição de fase topológica, mas as funções de envelope exibem zeros irremovíveis na fase topológica ($C=1$), manifestando a obstrução de Wannier no espaço real.
• Construção de Modelos Locais para WSe$_2$ de Camadas Torcidas: Os autores aplicam o método ao WSe$_2$ de camadas torcidas (tWSe$_2$) através de uma gama de ângulos de torção ($1^\circ$ a $5^\circ$).
  • Eles identificam que as duas bandas superiores requerem descrições locais diferentes dependendo do ângulo de torção e da escolha da partição (sítios de empilhamento MM, MX ou XM).
  • No regime topológico ($\theta > 1.51^\circ$), nenhum orbital local único abrange as bandas devido ao número de Chern total finito.
  • Ao selecionar subconjuntos de CNOs (por exemplo, um conjunto multicêntrico combinando orbitais MM, MX e XM) que maximizam o peso espectral total e satisfazem a condição de completude (determinante não nulo da matriz de sobreposição $S_k$), eles constroem modelos de tight-binding.
  • Wannierização de Passo Único (Single-Shot): Esses subconjuntos de CNOs servem como estados de teste para um procedimento de Wannierização de passo único. Os modelos de tight-binding resultantes reproduzem não apenas a dispersão de banda, mas também a geometria quântica (curvatura de Berry) do modelo contínuo.
  • Compromissos (Trade-offs): Os autores comparam um modelo centrado em MM (3 orbitais em um sítio) com um modelo multicêntrico (orbitais em diferentes sítios). O modelo multicêntrico produz funções de Wannier mais localizadas (saltos de curto alcance) mas requer um cluster maior para tratamentos de muitos corpos que preservam a simetria (como o cluster DMFT), enquanto o modelo centrado em MM preserva a simetria $C_3$ mais facilmente, mas possui interações de longo alcance.

Significância e Alegações
O artigo afirma que os CNOs oferecem uma medida direta, derivada da função de onda, do caráter multi-orbital das bandas alvo. Ao inverter a abordagem convencional de inferir orbitais a partir da química ou dos autovalores de simetria, os CNOs permitem que os graus de liberdade locais sejam inferidos diretamente do próprio manifold de bandas.

A significância reside em:

1. Otimalidade: Os CNOs fornecem os estados de teste ideais para a Wannierização, minimizando o número de orbitais necessários para abranger as bandas alvo enquanto preservam a geometria quântica.
2. Poder Diagnóstico: O espectro de ocupação e os zeros das funções de envelope servem como diagnósticos para obstruções topológicas e limites atômicos obstruídos, distinguindo entre zeros removíveis (simetria) e irremovíveis (topológicos).
3. Base para Modelos Interagentes: Como os CNOs diagonalizam a matriz de correlação da célula unitária, eles fornecem uma base natural para interações projetadas de banda. A hierarquia de ocupação $\lambda_a$ ordena as escalas de interação, oferecendo uma rota para construir modelos de Hubbard multi-orbitais e definir espaços de impureza para métodos de embutimento quântico (como DMFT) sem escolhas arbitrárias.
4. Generalidade: O método aplica-se tanto a modelos de rede quanto de contínuo (via uma grade no espaço real derivada da base de ondas planas) e é compatível com cálculos de primeiros princípios, pois requer apenas o projetor de banda como entrada.

Os autores concluem que esta perspectiva une as funções de onda de primeiros princípios e os modelos contínuos com os graus de liberdade locais mínimos necessários para teorias efetivas de materiais quânticos topológicos e correlacionados.