Disorder-induced symmetry breaking in moiré bands of marginally twisted bilayer MoS$_2$

Enunciado do Problema
Os dicalcogenetos de metais de transição torcidos (tTMDs) abrigam potenciais de moiré altamente ajustáveis e bandas planas que suportam fases eletrônicas fortemente correlacionadas. Embora modelos teóricos predigam fenômenos quânticos emergentes ricos nesses sistemas, o papel específico da desordem intrínseca na moldagem dessas propriedades permanece amplamente não resolvido. Diferente do grafeno bicamada torcido, que requer um "ângulo mágico" específico para bandas planas, os tTMDs de pequeno ângulo exibem larguras de banda continuamente decrescentes próximo ao empilhamento paralelo (0°) e antiparalelo (60°), oferecendo uma ampla plataforma para investigar fases quânticas. No entanto, a desordem intrínseca — como defeitos pontuais e poças de carga — é uma característica inevitável de materiais bidimensionais. A natureza microscópica da desordem nos tTMDs e seu impacto na estrutura eletrônica de superredes de moiré, particularmente em relação à quebra de simetria, recebeu atenção limitada.

Metodologia
Os autores investigaram o MoS$_2$ bicamada torcido (tb-MoS$_2$) marginalmente com um ângulo de torção de $\theta \approx 0,95^\circ$ usando microscopia de tunelamento por varredura (STM) e espectroscopia (STS).

• Experimental: As amostras foram fabricadas usando uma técnica de tear-and-stack para criar bicamadas sobre um substrato de grafite/SiO$_2$. A topografia de STM foi usada para caracterizar a superrede de moiré, enquanto a STS foi empregada para mapear a densidade local de estados (LDOS). Os pesquisadores mediram espectros de condutância diferencial ($dI/dV$) em regiões de empilhamento de alta simetria (MM, MX, XM e paredes de domínio) e geraram mapas de "limiar" (onset maps) espaciais para visualizar a energia na qual a LDOS se eleva acima de um nível de ruído para ambas as bandas de valência e condução.
• Teórico: O estudo utilizou um modelo de contínuo para o MoS$_2$ torcido que incorpora efeitos de relaxação atômica derivados de cálculos de teoria do funcional da densidade (DFT) de primeira ordem. Para modelar as observações experimentais, os autores desenvolveram um modelo de desordem eletrostática tratando o sistema como camadas dielétricas alternadas sobre um substrato metálico de grafite. Este modelo contabilizou defeitos carregados distribuídos aleatoriamente (especificamente vacâncias de enxofre) para calcular o potencial eletrostático resultante e seu efeito nas bandas de moiré.

Resultos Principais

1. Deslocamentos de Limiar de Quebra de Simetria: Medições de STS revelaram um desdobramento inesperado nas energias de limiar das bandas de valência e condução entre regiões de empilhamento MX e XM, que são simétricas em espelho. Especificamente, uma diferença de $15 \pm 4$ meV foi observada para a banda de valência e $15 \pm 6$ meV para a banda de condução. Em um sistema prístino e simétrico, essas regiões deveriam exibir energias de limiar idênticas.
2. Desordem Espacialmente Correlacionada: Os mapas de limiar exibiram variações de energia espacialmente correlacionadas em escalas de comprimento superiores ao período de moiré. A análise de autocorrelação desses mapas mostrou uma dependência radial consistente com um potencial de $1/r^3$, característico de dipolos formados por defeitos carregados e suas cargas de imagem no substrato.
3. Estimativa de Densidade de Defeitos: Ao ajustar os dados de autocorrelação ao modelo de desordem eletrostática, os autores inferiram uma densidade de carga aleatória de fundo de aproximadamente $1 \times 10^{12}$ cm$^{-2}$ a $3 \times 10^{11}$ cm$^{-2}$. Essa faixa alinha-se com as densidades relatadas de vacâncias de enxofre em monocamadas de MoS$_2$ exfoliadas mecanicamente.
4. Concordância Teórica: Cálculos do modelo de contínuo incorporando um campo de deslocamento elétrico entre camadas (simulando o potencial proveniente da desordem inferida) reproduziram os resultados experimentais. Os cálculos mostraram que o campo elétrico induz polarização de camada nas bandas de moiré, causando um desdobramento das energias de limiar entre os sítios MX e XM. A densidade local de estados calculada, quando alargada para considerar a desordem e efeitos térmicos, correspondeu aos mapas de limiar e aos deslocamentos espectrais experimentais quantitativamente.

Significância e Alegações
O artigo afirma que mesmo níveis modestos de desordem eletrostática, inerentes aos TMDs exfoliados, desempenham um papel crítico na determinação da estrutura eletrônica de materiais de moiré em escala nanométrica. A principal significância destas descobertas é a demonstração de que a desordem pode induzir a quebra de simetria nas bandas de moiré, transformando o caráter dos estados eletrônicos. Especificamente, as bandas de valência superiores, que são teoricamente previstas como não polarizadas e do tipo grafeno em um sistema prístino, mostram-se completamente polarizadas em camada e amplamente separadas em energia devido à presença de defeitos carregados.

Os autores concluem que esses efeitos induzidos pela desordem ocorrem em uma escala de energia de dezenas de meV, o que é comparável às larguras de banda das próprias bandas planas. Consequentemente, o artigo afirma que deve-se ter cautela ao interpretar resultados experimentais em tTMDs, pois a desordem nativa pode alterar fundamentalmente a realização de topologias de banda plana e fases correlacionadas específicas. O trabalho destaca que o relaxamento estrutural e a desordem eletrostática devem ser simultaneamente considerados para alcançar concordância quantitativa entre teoria e experimento nestes sistemas.