Phase diagram of the massless Gross-Neveu model with two components

Sintesi Tecnica: Diagramma di Fase del Modello di Gross-Neveu Senza Massa con Due Componenti

Problematica
Il modello di Gross-Neveu (GN), una teoria quantistica dei campi di $N$ sapori di fermioni di Dirac senza massa in 1+1 dimensioni con interazioni quattro-fermione scalare-scalare, funge da modello paradigmatico per lo studio della rottura dinamica della simmetria, dell'asintotica libertà e dell'integrabilità. Mentre il limite per $N$ grande del modello GN standard (dove un potenziale chimico $\mu$ agisce su tutti i $N$ sapori) è stato ampiamente esplorato a temperatura e densità finite, recenti sviluppi di Benini et al. hanno introdotto una generalizzazione in cui il potenziale chimico agisce solo su un sottoinsieme di $a$ sapori ($a = \nu N$, con $0 < \nu \le 1$). Questa configurazione a "due componenti" distingue tra fermioni "carichi" e "neutri", modellando efficacemente le "bande di valenza parzialmente riempite" nella materia densa.

Nonostante esistano risultati esatti per questo modello generalizzato a temperatura zero ($T=0$), il diagramma di fase a temperatura finita ($T > 0$) rimaneva ignoto. Nello specifico, non era chiaro come l'introduzione della frazione di riempimento $\nu$ alteri la struttura di fase, in particolare riguardo all'esistenza di fasi disomogenee (cristalline) e alla natura delle transizioni di fase (del primo ordine rispetto al secondo ordine) nel piano $(\mu, T)$.

Metodologia
Gli autori costruiscono l'intero diagramma di fase nel limite $N$ grande senza eseguire un calcolo Hartree-Fock (HF) termico completo, che sarebbe computazionalmente proibitivo a causa della necessità di minimizzare il potenziale grand canonical rispetto a molteplici parametri. Inve al, lo studio si basa su una serie di analisi di stabilità per determinare i confini di fase:

1. Approssimazione Omogenea: Lo studio inizia assumendo l'invarianza traslazionale non rotta. La densità del potenziale grand canonical rinormalizzato è costruita come una combinazione lineare dei potenziali a una componente: $\Psi^{(2)}_{\text{eff}} = (1-\nu)\Psi^{(1)}_{\text{eff}}(\beta, 0) + \nu\Psi^{(1)}_{\text{eff}}(\beta, \mu)$.

   • Confini Perturbativi: Il confine tra la fase chiralmente simmetrica ($M=0$) e quella rotta ($M>0$) è determinato tramite l'espansione di Ginzburg-Landau (GL), specificamente risolvendo $A_2(\nu; T, \mu) = 0$.
   • Punti Tricritici: La transizione dal comportamento del secondo ordine a quello del primo ordine è individuata trovando dove il coefficiente quartico GL $A_4$ si annulla.
   • Linee del Primo Ordine: Queste sono identificate individuando i punti nel piano $(\mu, T)$ in cui il potenziale grand canonical presenta due minimi degeneri.

2. Fasi Disomogenee (Formazione di Cristalli): Per mappare il diagramma di fase completo includendo i condensati disomogenei, gli autori impiegano specifiche analisi di stabilità:

   • Confini Orizzontali (Senza Massa $\to$ Cristallo): Determinati analizzando l'instabilità della fase senza massa contro modulazioni spaziali (potenziali periodici) usando la teoria delle perturbazioni nel potenziale scalare.
   • Confini Verticali (Massivo $\to$ Cristallo/Solitone): Determinati analizzando la termodinamica di un singolo barione (barione DHN). Il confine corrisponde al punto in cui la fase massiva omogenea diventa instabile alla formazione di un singolo solitone. Ciò comporta la minimizzazione del potenziale grand canonical di un ansatz di singolo barione rispetto a un parametro variazionale $y$.

Lo studio è diviso in due regimi basati sulla frazione di riempimento $\nu$: una regione "connessa" ($\nu \ge \nu_0$) e una regione "disconnessa" ($\nu < \nu_0$), dove $\nu_0$ è un valore di biforcazione critico.

Risultati Chiave

• Biforcazione a $\nu_0 = 0.8271$: Il diagramma di fase omogeneo subisce un cambiamento topologico a una frazione di riempimento critica $\nu_0 \approx 0.8271$.

  • Regime Connesso ($\nu_0 \le \nu \le 1$): Il confine di fase perturbativo del secondo ordine (che separa le fasi senza massa e massiva) termina in un punto tricritico, oltre il quale inizia una linea del primo ordine. Questa struttura è qualitativamente simile al modello GN standard a una componente ($\nu=1$).
  • Regime Disconnesso ($\nu < \nu_0$): Il punto tricritico scompare. La linea perturbativa del secondo ordine si estende indefinitamente senza terminare. Inveve, una linea del primo ordine emerge all'interno della fase massiva, terminando in un endpoint critico che si trova al di sotto del confine perturbativo.

• Diagramma di Fase Completo con Disomogeneità:

  • Per $\nu \ge \nu_0$: L'inclusione delle fasi disomogenee sostituisce la linea del primo ordine dell'approssimazione omogenea con una struttura a "cuspide" formata da due linee del secondo ordine. Queste linee racchiudono una fase cristallina. Una linea separa la fase senza massa dal cristallo, e l'altra separa la fase massiva dal cristallo. Il punto tricritico del calcolo omogeneo rimane una caratteristica del diagramma completo.
  • Per $\nu < \nu_0$: La linea del primo ordine del calcolo omogeneo è sostituita da due linee del secondo ordine che emanano dall'endpoint critico, racchiudendo una regione cristallina. Al diminuire di $\nu$, la dimensione di questa regione disomogenea (cristallina) si contrae rapidamente.
  • Nuova Biforcazione (a $\nu = 3/4$): Nel regime disconnesso, emerge una nuova caratteristica in cui la linea di instabilità "orizzontale" (massivo $\to$ cristallo) interseca il confine omogeneo perturbativo in un nuovo punto $Q$. Ciò crea una struttura complessa che coinvolge quattro distinte frontiere del secondo ordine: senza massa-massiva, senza massa-cristallo, massiva-cristallo e massiva-solitone singolo.

• Parametri Critici: Il documento fornisce coordinate numeriche specifiche per i punti critici, i punti tricritici e gli endpoint per vari valori di $\nu$ (ad esempio, per $\nu=1$, il punto tricritico è a $\mu \approx 0.6082, T \approx 0.3183$; per $\nu_0$, la biforcazione avviene a $\mu \approx 0.7202, T \approx 0.2240$).

Significato e Rivendicazioni
Il lavoro rivendica di aver colmato una lacuna nella comprensione del diagramma di fase a temperatura finita del modello GN con riempimento parziale dei sapori. Utilizzando analisi di stabilità anziché calcoli HF completi, gli autori dimostrano che il diagramma di fase completo, incluse le fasi disomogenee, può essere costruito efficientemente.

La principale significatività risiede nell'identificazione della biforcazione a $\nu_0 \approx 0.8271$, che segna un cambiamento qualitativo nella topologia del diagramma di fase. Lo studio conferma che il concetto di "bande di valenza parzialmente riempite" (generalizzato dal barione DHN a $T=0$) porta a una ricca struttura di fase dove l'interazione tra la frazione di riempimento $\nu$ e la temperatura determina se il sistema esibisce confini di fase connessi o disconnessi. I risultati evidenziano che la libertà extra fornita dal parametro $\nu$ arricchisce significativamente il diagramma di fase, introducendo nuovi punti di biforcazione e modificando la stabilità delle fasi cristalline, in modo analogo ma distinto rispetto al modello GN massivo.