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

Résumé Technique : Diagramme de Phase du Modèle de Gross-Neveu Sans Masse à Deux Composantes

Énoncé du Problème
Le modèle de Gross-Neveu (GN), une théorie quantique des champs de $N$ saveurs de fermions de Dirac sans masse en 1+1 dimensions avec des interactions de type quatre-fermions scalaire-scalaire, sert de modèle de base paradigmatique pour l'étude de la brisure de symétrie dynamique, de l'asymptotique liberté et de l'intégrabilité. Bien que la limite de grand $N$ du modèle GN standard (où un potentiel chimique $\mu$ agit sur toutes les $N$ saveurs) ait été approfondie à température et densité finies, des développements récents par Benini et al. ont introduit une généralisation où le potentiel chimique n'agit que sur un sous-ensemble de $a$ saveurs ($a = \nu N$, avec $0 < \nu \le 1$). Cette configuration à « deux composantes » distingue les fermions « chargés » des fermions « neutres », modélisant efficacement les « bandes de valence partiellement remplies » dans la matière dense.

Malgré l'existence de résultats exacts pour ce modèle généralisé à température nulle ($T=0$), le diagramme de phase à température finie ($T > 0$) restait inconnu. Plus précisément, il n'était pas clair comment l'introduction de la fraction de remplissage $\nu$ modifie la structure de phase (notamment l'existence de phases inhomogènes ou cristallines) et la nature des transitions de phase (premier ordre vs second ordre) dans le plan $(\mu, T)$.

Méthodologie
Les auteurs construisent le diagramme de phase complet à grand $N$ sans effectuer un calcul complet de Hartree-Fock (HF) thermique, ce qui serait computationnellement prohibitif en raison de la nécessité de minimiser le potentiel grand canonique sur plusieurs paramètres. Au lieu de cela, l'étude repose sur une série d'analyses de stabilité pour déterminer les frontières de phase :

1. Approximation Homogène : L'étude commence en supposant une invariance par translation non brisée. La densité du potentiel grand canonique renormalisé est construite comme une combinaison linéaire des potentiels à une composante : $\Psi^{(2)}_{\text{eff}} = (1-\nu)\Psi^{(1)}_{\text{eff}}(\beta, 0) + \nu\Psi^{(1)}_{\text{eff}}(\beta, \mu)$.

   • Frontières Perturbatives : La frontière entre la phase chirale symétrique ($M=0$) et la phase brisée ($M>0$) est déterminée via une expansion de Ginzburg-Landau (GL), spécifiquement en résolvant $A_2(\nu; T, \mu) = 0$.
   • Points Tricritiques : La transition d'un comportement de second ordre vers un comportement de premier ordre est localisée en trouvant où le coefficient quartique de GL, $A_4$, s'annule.
   • Lignes de Premier Ordre : Celles-ci sont identifiées en localisant les points dans le plan $(\mu, T)$ où le potentiel grand canonique présente deux minima dégénérés.

2. Phases Inhomogènes (Formation de Cristal) : Pour cartographier le diagramme de phase complet incluant les condensats inhomogènes, les auteurs emploient des analyses de stabilité spécifiques :

   • Frontières Horizontales (Sans Masse $\to$ Cristal) : Déterminées par l'analyse de l'instabilité de la phase sans masse contre les modulations spatiales (potentiels périodiques) en utilisant la théorie des perturbations dans le potentiel scalaire.
   • Frontières Verticales (Massive $\to$ Cristal/Soliton) : Déterminées par l'analyse thermodynamique d'un baryon unique (baryon DHN). La frontière correspond au point où la phase massive homogène devient instable face à la formation d'un soliton unique. Cela implique la minimisation du potentiel grand canonique d'une ansatz de baryon unique par rapport à un paramètre variationnel $y$.

L'étude est divisée en deux régimes basés sur la fraction de remplissage $\nu$ : une région « connectée » ($\nu \ge \nu_0$) et une région « déconnectée » ($\nu < \nu_0$), où $\nu_0$ est une valeur de bifurcation critique.

Résultats Clés

• Bifurcation à $\nu_0 = 0,8271$ : Le diagramme de phase homogène subit un changement topologique à une fraction de remplissage critique $\nu_0 \approx 0,8271$.

  • Régime Connecté ($\nu_0 \le \nu \le 1$) : La frontière de phase de second ordre perturbative (séparant les phases sans masse et massive) se termine par un point tricritique, au-delà duquel une ligne de premier ordre commence. Cette structure est qualitativement similaire au modèle GN standard à une composante ($\nu=1$).
  • Régime Déconnecté ($\nu < \nu_0$) : Le point tricritique disparaît. La ligne de second ordre perturbative s'étend indéfiniment sans se terminer. À la place, une ligne de premier ordre émerge au sein de la phase massive, se terminant à un point critique final qui se situe en dessous de la frontière perturbative.

• Diagramme de Phase Complet avec Inhomogénéités :

  • Pour $\nu \ge \nu_0$ : L'inclusion des phases inhomogènes remplace la ligne de premier ordre de l'approximation homogène par une structure en « cusp » (pointe) formée par deux lignes de second ordre. Ces lignes entourent une phase cristalline. Une ligne sépare la phase sans masse du cristal, et l'autre sépare la phase massive du cristal. Le point tricritique du calcul homogène reste une caractéristique du diagramme complet.
  • Pour $\nu < \nu_0$ : La ligne de premier ordre de la phase homogène est remplacée par deux lignes de second ordre émanant du point critique final, entourant une région cristalline. À mesure que $\nu$ diminue, la taille de cette région inhomogène (cristal) rétrécit rapidement.
  • Nouvelle Bifurcation (à $\nu = 3/4$) : Dans la région déconnectée, une caractéristique inédite apparaît où la ligne d'instabilité horizontale (massive $\to$ cristal) intersecte la frontière homogène perturbative en un nouveau point $Q$. Cela crée une structure complexe impliquant quatre frontières de second ordre distinctes : sans masse-massive, sans masse-cristal, massive-cristal, et massive-soliton unique.

• Paramètres Critiques : L'article fournit les coordonnées numériques spécifiques des points critiques, des points tricritiques et des points finaux pour diverses valeurs de $\nu$ (par exemple, pour $\nu=1$, le point tricritique est à $\mu \approx 0,6082, T \approx 0,3183$ ; à $\nu_0$, la bifurcation se produit à $\mu \approx 0,7202, T \approx 0,2240$).

Signification et Revendications
L'article prétend combler une lacune dans la compréhension du diagramme de phase à température finie du modèle GN avec remplissage partiel de saveur. En utilisant des analyses de stabilité plutôt que des calculs HF complets, les auteurs démontrent que le diagramme de phase complet, incluant les phases inhomogènes, peut être construit efficacement.

La principale signification réside dans l'identification de la bifurcation à $\nu_0 \approx 0,8271$, qui marque un changement qualitatif de la topologie du diagramme de phase. L'étude confirme que le concept de « bandes de valence partiellement remplies » (généralisé à partir du baryon DHN à $T=0$) conduit à une structure de phase riche où l'interaction entre la fraction de remplissage $\nu$ et la température détermine si le système présente des frontières de phase connectées ou déconnectées. Les résultats soulignent que la liberté supplémentaire fournie par le paramètre $\nu$ enrichit considérablement le diagramme de phase, introduisant de nouveaux points de bifurcation et modifiant la stabilité des phases cristallines, de manière analogue mais distincte du modèle GN massif.