Wigner crystallization in Bernal bilayer graphene

Resumo Técnico: Cristalização de Wigner em grafeno bicamada de Bernal

Enunciado do Problema
O artigo investiga a natureza e a estabilidade do estado de cristal de Wigner (WC) em grafeno bicamada de Bernal (BBG) sob a influência de um campo de deslocamento perpendicular ($U$). Embora experimentos recentes em multicamadas de grafeno torcido e não torcido tenham revelado fases correlacionadas ricas, a interação específica entre a estrutura de bandas única do BBG (especificamente a dispersão "chapéu mexicano" induzida por $U$), a curvatura de Berry e o empenamento trigonal (trigonal warping) permanece uma questão teórica em aberto. Os autores visam determinar as fronteiras de fase do WC em termos de densidade de elétrons ($n$), campo de deslocamento ($U$) e temperatura ($T$), além de caracterizar as propriedades do estado fundamental, incluindo a magnetização orbital e o ordenamento de vales, distinguindo-os de fases de WC convencionais.

Metodologia
Os autores empregam uma descrição de oscilador harmônico (HO) de partícula única para modelar o estado de WC. Neste arcabouço, cada elétron na rede triangular do WC é tratado como uma partícula confinada em um poço de potencial parabólico criado pela repulsão Coulombiana de seus vizinhos. A força desse confinamento é determinada pela densidade de elétrons.

O estudo prossegue em dois regimes teóricos principais:

1. Modelo com Simetria de Rotação (Sem Empenamento Trigonal): Os autores primeiro negligenciam o empenamento trigonal para tratar a dispersão de baixa energia como rotacionalmente simétrica em torno dos pontos $K$ e $K'$. Eles derivam uma equação de Schrödinger de partícula única efetiva no espaço de momento, incorporando a conexão de Berry ($A(p)$) como um potencial vetorial. Isso permite soluções analíticas para estimar espectros de energia, funções de onda e o critério de Lindemann para fusão.
2. Modelo de Empenamento Trigonal: Os autores então incluem os efeitos do empenamento trigonal, que quebra a simetria de rotação e divide a dispersão de "chapéu mexicano" em "mini-vales" discretos (um bolso central e três bolsos laterais). Para este caso, eles realizam a diagonalização numérica da equação de Schrödinger bidimensional em uma grade no espaço de momento para computar a razão de Lindemann e mapear o diagrama de fase.

A estabilidade do WC é avaliada usando o critério de Lindemann, onde a fusão ocorre quando a largura raiz-quadrada do valor médio da função de onda do elétron excede uma fração crítica ($\eta_c \approx 0,23$) da constante de rede. Os autores também calculam integrais de troca para estimar temperaturas de ordenamento e avaliam a magnetização através do valor esperado do operador de momento angular orbital.

Principais Contribuições e Resultados

• Dependência do Campo de Deslocamento e Blindagem Dielétrica: O artigo estabelece que, a zero campo de deslocamento, o estado de WC é precluso no BBG devido à massa de banda evanescente em baixa energia e à forte blindagem dielétrica interbanda. Um estado de WC só emerge quando um campo de deslocamento suficiente é aplicado para abrir um gap de banda e achatar o fundo da banda.
• Magnetização Orbital Espontânea (Sem Empenamento Trigonal): Na ausência de empenamento trigonal, os autores demonstram que a curvatura de Berry impulsiona uma transição no momento angular orbital ($\ell$) dos elétrons do WC. À medida que o campo de deslocamento $U$ aumenta, o fluxo de Berry envolvido pela função de onda aumenta. Quando este fluxo excede um valor crítico ($\tilde{\Phi} > 1/2$), o estado fundamental transita de $\ell=0$ para $\ell=1$ (ou $\ell=-1$ no vale $K'$). Esta transição manifesta-se como um salto descontínuo na magnetização orbital do estado de WC.
• Fusão Duplamente Reentrante (Com Empenamento Trigonal): Quando o empenamento trigonal é incluído, o diagrama de fase exibe um comportamento único de "dupla reentrada" como função da densidade de elétrons em baixo $U$.
  1. Em densidades muito baixas, os elétrons cristalizam nos "bolsos laterais" de menor energia.
  2. À medida que a densidade aumenta, o WC funde-se em um líquido de Fermi (FL) em uma densidade crítica $n_{side}^c$.
  3. Em uma densidade mais alta $n^*$, torna-se energeticamente favorável para os elétrons ocuparem o "bolso central" de massa mais pesada, fazendo com que o sistema recristalize em um estado de WC.
  4. Finalmente, em uma densidade ainda mais alta $n_{central}^c$, este WC do bolso central funde-se novamente.
• Ordenamento de Minivales: A presença de empenamento trigonal quebra a simetria de rotação, tornando o momento angular $\ell$ mal definido. Em vez disso, o estado fundamental é caracterizado por um padrão espacial não trivial de "ordenação de mini-vales". Baseando-se em trabalho recente (Ref. [102]), os autores preveem que os elétrons se polarizam em mini-vales específicos em um padrão de listras alternadas na rede de WC.
• Fronteiras de Fase: Os autores fornecem estimativas numéricas para as fronteiras de fase. A densidade máxima para o estado de WC no BBG é estimada na ordem de $1,2 \times 10^{10} \text{ cm}^{-2}$, ocorrendo em um campo de deslocamento de aproximadamente $0,7 \text{ V/nm}$. A temperatura máxima de fusão é estimada na ordem de $10 \text{ K}$ (diminuindo com a densidade), com a região reentrante do bolso central apresentando uma temperatura máxima inferior de $\sim 1 \text{ K}$.

Significância e Alegações
O artigo afirma fornecer um diagrama de fase teórico abrangente para o estado de WC no BBG, destacando como a geometria quântica não trivial (curvatura de Berry) e a quebra de simetria de rede (empenamento trigonal) alteram fundamentalmente a natureza do cristal de Wigner.

• Efeitos Orbitais: O trabalho identifica um mecanismo para magnetização orbital espontânea em um WC impulsionado puramente pela curvatura de Berry e campos de deslocamento, distinto do ordenamento baseado em spin.
• Comportamento de Fase Complexo: A previsão de fusão "duplamente reentrante" é apresentada como um recurso novel decorrente da competição entre diferentes bolsos de banda (central vs. lateral) com distintas massas efetivas e energias.
• Contexto Experimental: Os autores observam que a faixa de densidade prevista para o WC no BBG ($n \lesssim 10^{10} \text{ cm}^{-2}$) é inferior às densidades onde estados isolantes foram observados em experimentos recentes, sugerindo que as observações experimentais atuais no BBG podem ser dominadas por desordem ou corresponder a regimes diferentes. Eles sugerem que o comportamento reentrante e o ordenamento de mini-vales podem ser mais facilmente observáveis em grafeno de empilhamento rhomboédrico com maior número de camadas ($N$), onde a estrutura de bandas é mais plana e o fluxo de Berry é maior ($N\pi$).
• Complementaridade: O artigo é explicitamente declarado como complementar a um trabalho companheiro (Ref. [1]) que foca no ordenamento de spin; este artigo foca nos efeitos orbitais e no diagrama de fase detalhado incluindo o empenamento trigonal.

Os autores mantêm um tom modesto quanto às limitações da aproximação de HO, reconhecendo que ela pode superestimar a estabilidade da fase de WC em comparação com cálculos de Monte Carlo quântico mais rigorosos, e que o limite "Einstein phonon" usado aqui não pode reproduzir fases com números de Chern não triviais (como o cristal de efeito Hall anômalo quântico).