Trimers in the Extended Hubbard Model

Resumo Técnico: Instabilidade de Trímero no Modelo de Hubbard Estendido

Definição do Problema
Este trabalho investiga a estabilidade da fase isolante ferrimagnética de Lieb no modelo de Hubbard estendido em uma cadeia de trímeros acoplados bipartida em meio preenchimento (um elétron por sítio). Embora o modelo de Hubbard padrão em redes bipartidas com salto de vizinho mais próximo e repulsão on-site $U$ seja conhecido por exibir um estado fundamental ferrimagnético robusto previsto pelo teorema de Lieb, a introdução de uma interação Coulombiana de vizinho mais próximo $V$ deve induzir instabilidades de carga. O problema central abordado é a evolução deste estado ferrimagnético conforme $V$ aumenta, identificando especificamente as condições sob as quais a fase isolante homogênea torna-se instável à separação de fases (PS) e caracterizando a natureza das fases coexistentes resultantes.

Metodologia
Os autores utilizam o método de Grupo de Renormalização de Matriz de Densidade (DMRG) para resolver o Hamiltoniano de Hubbard estendido para cadeias de trímeros acoplados com condições de contorno abertas. O sistema consiste em $L$ trímeros, onde cada trímero contém três sítios rotulados como $A$, $B_1$ e $B_2$. O Hamiltoniano inclui:

• Amplitude de salto intra-trímero $t$ (definida como $t=1$).
• Amplitude de salto inter-trímero $t'$ (definida como $t'=0.8$).
• Repulsão Coulombiana on-site $U$.
• Repulsão Coulombiana de vizinho mais próximo $V$.

Os cálculos foram realizados em cadeias com tamanhos de até $L=60$ (total de $N=180$ sítios), utilizando dimensões de ligação máximas entre 512 e 2048, com um erro de truncamento de $\sim 1 \times 10^{-6}$. Para investigar as propriedades do estado fundamental, os autores analisaram:

1. Médias espaciais e distribuições locais da ocupação de partículas ($\langle n_i \rangle$) e ocupação de doublons ($\langle d_i \rangle = \langle n_{i\uparrow} n_{i\downarrow} \rangle$).
2. Spin total por trímero ($S/L$), calculado via funções de correlação de spin no setor $S_z=0$.
3. Compressibilidade ($\kappa$), derivada de diferenças de energia de tamanho finito para identificar regimes incompressíveis (gapados) versus compressíveis (metálicos/separados de fase).
4. Métricas de homogeneidade, especificamente $\Delta_{max}$, definido como o desvio máximo do número de partículas em qualquer trímero em relação à densidade homogênea de 3 elétrons por trímero.

Contribuições Principais e Resultados

• Análise do Limite Atômico ($t=t'=0$): Os autores primeiro analisaram as configurações estáticas no limite atômico. Eles identificaram que o aumento de $V$ a partir de zero favorece uma transição de um estado padrão de meio preenchimento (um elétron por sítio) para um estado onde os sítios $A$ estão vazios e os sítios $B$ estão duplamente ocupados. Esta transição ocorre em um valor crítico $V = U/4$, onde a energia do estado de meio preenchimento ($E/L = 4V$) é igual à do estado de transferência de carga ($E/L = U$).

• O Caso $U=0$: Para repulsão on-site nula, qualquer $V$ finito induz separação de fases. O sistema separa-se em duas fases distintas: uma fase rica em doublons na sub-rede $B$ (isolante), e outra rica em doublons na sub-rede $A$ (metálica). O spin total por trímero é zero ($S/L = 0$) neste regime.

• Estabilidade da Fase Ferrimagnética de Lieb ($U \neq 0$):

  • Para $V=0$, o sistema exibe o esperado estado isolante ferrimagnético de Lieb com $S/L = 1/2$.
  • À medida que $V$ aumenta, o sistema permanece em uma fase isolante homogênea até um ponto crítico. No entanto, ocorre uma transferência significativa de carga dos sítios $A$ para os sítios $B$, aumentando a densidade de doublons na sub-rede $B$.
  • Os autores identificam uma linha de crossover próxima a $V \approx U/4$. Abaixo desta linha, o sistema é homogêneo. Entre a linha de crossover e a linha de separação de fases, o sistema permanece homogêneo, mas exibe um caráter "rico em doublons" com aumento da transferência de carga, mantendo a ordem ferrimagnética ($S/L = 1/2$).

• Separação de Fases e Coexistência:

  • Acima de uma linha crítica de separação de fases (PS) (observada em $V \gtrsim U/4$), o estado homogêneo torna-se instável. O sistema entra em uma região de coexistência onde duas fases com distintas ocupações de trímero existem espacialmente.
  • Para $U \gtrsim 1$: A região de coexistência consiste em:
    1. Uma fase isolante onde os sítios $B$ são predominantemente duplamente ocupados (aproximadamente 4 elétrons por trímero), correspondendo a um estado rico em doublons.
    2. Uma fase ferromagnética insaturada metálica onde a distribuição de carga é mais uniforme, mas carece da condição estrita de meio preenchimento. Esta fase exibe um momento magnético líquido ($S/L > 0$), resultante da interação entre a itinerância de carga e o alinhamento ferrimagnético dos momentos nos sítios $A$ e $B$.
  • Para $U \lesssim 1$: A região de coexistência assemelha-se ao caso $U=0$, consistindo em duas fases paramagnéticas ($S/L \approx 0$), uma rica em doublons nos sítios $B$ e outra nos sítios $A$.

• Comportamento de Alto $U$: Os autores estimam que a linha de crossover e a linha de PS se fundem em aproximadamente $U \sim 16$. Para $U > 16$, o sistema transita diretamente da fase ferrimagnética homogênea para a região de separação de fases, indicado por um salto nas médias espaciais em vez de um ponto de inflexão.

Significância e Alegações
O artigo afirma demonstrar que a fase isolante ferrimagnética de Lieb é notavelmente robusta contra a repulsão de vizinho mais próximo $V$, persistindo bem além do ponto onde a densidade de doublons aumenta significativamente. A principal contribuição é o mapeamento do diagrama de fases para cadeias de trímeros acoplados, distinguindo entre um regime ferrimagnético homogêneo "rico em doublons" e a subsequente região de separação de fases.

Os autores afirmam que seus resultados esclarecem a natureza da instabilidade do estado fundamental neste tipo específico de geometria de rede. Eles destacam que, acima da linha de PS, o sistema suporta uma coexistência única de uma fase paramagnética isolante e uma fase ferromagnética insaturada metálica para $U$ intermediário a grande. O trabalho fornece uma estimativa quantitativa das linhas de PS e de crossover usando DMRG, confirmando que a transição é impulsionada pela competição entre o custo energético da dupla ocupação ($U$) e o ganho de energia ao evitar a repulsão de vizinho mais próximo ($V$). Os autores concluem que estas descobertas fornecem uma base para investigações futuras em modelos dopados e aspectos topológicos de tais transições de fase, sem propor novas realizações experimentais específicas além do arcabouço teórico.