Trimers in the Extended Hubbard Model

Sintesi Tecnica: Instabilità del Trimero nel Modello di Hubbard Esteso

Enunciato del Problema
Questo lavoro investiga la stabilità della fase isolante ferrimagnetica di Lieb all'interno del modello di Hubbard esteso su una catena di trimeri accoppiati bipartiti al semipieno (un elettrone per sito). Mentre lo standard modello di Hubbard su reticoli bipartiti con hopping a primo vicino e repulsione on-site $U$ è noto per esibire uno stato fondamentale ferrimagnetico robusto predetto dal teorema di Lieb, l'introduzione di una repulsione coulombiana tra primi vicini $V$ dovrebbe indurre instabilità di carica. Il problema centrale affrontato è l'evoluzione di questo stato ferrimagnetico all'aumentare di $V$, identificando specificamente le condizioni in cui la fase isolante omogenea diventa instabile rispetto alla separazione di fase (PS) e caratterizzando la natura delle fasi coesistenti risultanti.

Metodologia
Gli autori utilizzano il metodo Density Matrix Renormalization Group (DMRG) per risolvere l'Hamiltoniana di Hubbard estesa per catene di trimeri accoppiati con condizioni al contorno aperte. Il sistema consiste di $L$ trimeri, dove ogni trimero contiene tre siti etichettati come $A$, $B_1$ e $B_2$. L'Hamiltoniana include:

• Ampiezza di hopping intra-trimero $t$ (fissata a $t=1$).
• Ampiezza di hopping inter-trimero $t'$ (fissata a $t'=0.8$).
• Repulsione coulombiana on-site $U$.
• Repulsione coulombiana tra primi vicini $V$.

I calcoli sono stati eseguiti su catene con dimensioni fino a $L=60$ (totale $N=180$ siti), utilizzando dimensioni di legame massime tra 512 e 2048, con un errore di troncamento di $\sim 1 \times 10^{-6}$. Per sondare le proprietà dello stato fondamentale, gli autori hanno analizzato:

1. Medie spaziali e distribuzioni locali dell'occupazione di particelle ($\langle n_i \rangle$) e dell'occupazione di doublon ($\langle d_i \rangle = \langle n_{i\uparrow} n_{i\downarrow} \rangle$).
2. Spin totale per trimero ($S/L$), calcolato tramite funzioni di correlazione di spin nel settore $S_z=0$.
3. Compressione ($\kappa$), derivata dalle differenze di energia a dimensione finita per identificare regimi incomprimibili (gapati) rispetto a regimi comprimibili (metallici/separati di fase).
4. Metriche di omogeneità, specificamente $\Delta_{max}$, definita come la deviazione massima del numero di particelle in qualsiasi trimero dalla densità omogenea di 3 elettroni per trimero.

Contributi Chiave e Risultati

• Analisi del Limite Atomico ($t=t'=0$): Gli autori hanno prima analizzato le configurazioni statiche nel limite atomico. Hanno identificato che l'aumento di $V$ da zero favorisce una transizione dallo standard stato semipieno (un elettrone per sito) a uno stato in cui i siti $A$ sono vuoti e i siti $B$ sono doppiamente occupati. Questa transizione avviene a un valore critico $V = U/4$, dove l'energia dello stato semipieno ($E/L = 4V$) è uguale all'energia dello stato di trasferimento di carica ($E/L = U$).

• Il Caso $U=0$: Per una repulsione on-site nulla, qualsiasi $V$ finito induce separazione di fase. Il sistema si separa in due fasi distinte: una fase ricca di doublon sul sottoreticolo $B$ (isolante), e l'altra ricca di doublon sul sottoreticolo $A$ (metallica). Lo spin totale per trimero è zero ($S/L = 0$) in questo regime.

• Stabilità della Fase Ferrimagnetica di Lieb ($U \neq 0$):

  • Per $V=0$, il sistema esibisce l'atteso stato isolante ferrimagnetico di Lieb con $S/L = 1/2$.
  • All'aumentare di $V$, il sistema rimane in una fase isolante omogenea fino a un punto critico. Tuttavia, si verifica un significativo trasferimento di carica dai siti $A$ ai siti $B$, aumentando la densità di doublon sul sottoreticolo $B$.
  • Gli autori identificano una linea di crossover vicino a $V \approx U/4$. Al di sotto di questa linea, il sistema è omogeneo. Tra la linea di crossover e la linea di separazione di fase, il sistema rimane omogeneo ma esibisce un carattere "ricco di doublon" con un aumento del trasferimento di carica, mantenendo l'ordine ferrimagnetico ($S/L = 1/2$).

• Separazione di Fase e Coesistenza:

  • Al di sopra di una linea critica di separazione di fase (PS) (osservata a $V \gtrsim U/4$), lo stato omogeneo diventa instabile. Il sistema entra in una regione di coesistenza dove esistono spazialmente due fasi con diverse occupazioni di trimero.
  • Per $U \gtrsim 1$: La regione di coesistenza consiste in:
    1. Una fase isolante dove i siti $B$ sono prevalentemente doppiamente occupati (circa 4 elettroni per trimero), corrispondente a uno stato ricco di doublon.
    2. Una fase ferromagnetica insatura metallica dove la distribuzione di carica è più uniforme ma manca la stretta condizione di semipieno. Questa fase esibisce un momento magnetico netto ($S/L > 0$), derivante dall'interazione tra l'itineranza di carica e l'allineamento ferrimagnetico dei momenti sui siti $A$ e $B$.
  • Per $U \lesssim 1$: La regione di coesistenza somiglia al caso $U=0$, consistendo in due fasi paramagnetiche ($S/L \approx 0$), una ricca di doublon sui siti $B$ e l'altra sui siti $A$.

• Comportamento ad Alto $U$: Gli autori stimano che la linea di crossover e la linea PS si fondano approssimativamente a $U \sim 16$. Per $U > 16$, il sistema transita direttamente dalla fase ferrimagnetica omogenea alla regione di separazione di fase, come indicato da un salto nelle medie spaziali piuttosto che da un punto di inflessione.

Significato e Rivendicazioni
L'articolo sostiene che la fase isolante ferrimagnetica di Lieb è straordinariamente robusta contro la repulsione tra primi vicini $V$, persistendo ben oltre il punto in cui la densità di doublon aumenta significativamente. Il contributo primario è la mappatura del diagramma di fase per catene di trimeri accoppiati, distinguendo tra un regime ferrimagnetico "ricco di doublon" omogeneo e il successivo regime di separazione di fase.

Gli autori affermano che i loro risultati chiariscono la natura dell'instabilità dello stato fondamentale in questa specifica geometria di reticolo. Evidenziano che sopra la linea PS, il sistema supporta una coesistenza unica di una fase paramagnetica isolante e una fase ferromagnetica insatura metallica per $U$ intermedio o grande. Il lavoro fornisce una stima quantitativa delle linee di PS e di crossover usando DMRG, confermando che la transizione è guidata dalla competizione tra il costo energetico della doppia occupazione ($U$) e il guadagno energetico dall'evitare la repulsione tra primi vicini ($V$). Gli autori concludono che questi risultati forniscono una base per future investigazioni su modelli drogati e sugli aspetti topologici di tali transizioni di fase, senza proporre specifiche nuove realizzazioni sperimentali oltre il quadro teorico.