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

Resumen Técnico: Diagrama de Fases del Modelo de Gross-Neveu de Masa Nula con Dos Componentes

Planteamiento del Problema
El modelo de Gross-Neveu (GN), un modelo de teoría cuántica de campos de $N$ sabores de fermiones de Dirac sin masa en 1+1 dimensiones con interacciones de cuatro fermiones escalar-escalar, sirve como un modelo de juguete paradigmático para estudiar la ruptura dinámica de la simetría, la libertad asintótica y la integrabilidad. Mientras que el límite de $N$ grande del modelo GN estándar (donde un potencial químico $\mu$ actúa sobre todos los $N$ sabores) ha sido exhaustivamente explorado a temperatura y densidad finitas, desarrollos recientes por Benini et al. introdujeron una generalización donde el potencial químico actúa solo sobre un subconjunto de $a$ sabores ($a = \nu N$, con $0 < \nu \le 1$). Esta configuración de "dos componentes" distingue entre fermiones "cargados" y "neutros", modelando efectivamente "bandas de valencia parcialmente llenas" en materia densa.

A pesar de que existen resultados exactos para este modelo generalizado a temperatura cero ($T=0$), el diagrama de fases a temperatura finita ($T > 0$) permanecía desconocido. Específicamente, no estaba claro cómo la introducción de la fracción de llenado $\nu$ altera la estructura de fases (particularmente respecto a la existencia de fases inhomogéneas o cristalinas) y la naturaleza de las transiciones de fase (primer orden frente a segundo orden) en el plano $(\mu, T)$.

Metodología
Los autores construyen el diagrama de fases completo de $N$ grande sin realizar un cálculo completo de Hartree-Fock (HF) térmico, lo cual sería computacionalmente prohibitivo debido a la necesidad de minimizar el potencial de la gran canónica sobre múltiples parámetros. En su lugar, el estudio se basa en un conjunto de análisis de estabilidad para determinar las fronteras de fase:

1. Aproximación Homogénea: El estudio comienza asumiendo una invariancia de traslación no rota. La densidad del potencial de la gran canónica renormalizado se construye como una combinación lineal de los potenciales de un solo componente: $\Psi^{(2)}_{\text{eff}} = (1-\nu)\Psi^{(1)}_{\text{eff}}(\beta, 0) + \nu\Psi^{(1)}_{\text{eff}}(\beta, \mu)$.

   • Fronteras Perturbativas: La frontera entre la fase quiralmente simétrica ($M=0$) y la fase rota ($M>0$) se determina mediante una expansión de Ginzburg-Landau (GL), específicamente resolviendo $A_2(\nu; T, \mu) = 0$.
   • Puntos Tricríticos: La transición de un comportamiento de segundo orden a uno de primer orden se localiza encontrando dónde el coeficiente cuártico de GL, $A_4$, se anula.
   • Líneas de Primer Orden: Estas se identifican localizando puntos en el plano $(\mu, T)$ donde el potencial de la gran canónica exhibe dos mínimos degenerados.

2. Fases Inhomogéneas (Formación de Cristales): Para mapear el diagrama de fases completo incluyendo los condensados inhomogéneos, los autores emplean análisis de estabilidad específicos:

   • Fronteras Horizontales (Masless $\to$ Cristal): Determinadas mediante el análisis de la inestabilidad de la fase sin masa contra modulaciones espaciales (potenciales periódicos) usando teoría de perturbaciones en el potencial escalar.
   • Fronteras Verticales (Masivo $\to$ Cristal/Solitón): Determinadas mediante el análisis termodinámico de un solo barión (barión DHN). La frontera corresponde al punto donde la fase masiva homogénea se vuelve inestable ante la formación de un solo solitón. Esto implica minimizar el potencial de la gran canónica de un ansatz de un solo barión con respecto a un parámetro variacional $y$.

El estudio se divide en dos regímenes basados en la fracción de llenado $\nu$: una región "conectada" ($\nu \ge \nu_0$) y una región "desconectada" ($\nu < \nu_0$), donde $\nu_0$ es un valor de bifurcación crítico.

Resultados Clave

• Bifurcación en $\nu_0 = 0.8271$: El diagrama de fases homogéneo experimenta un cambio topológico en una fracción de llenado crítica $\nu_0 \approx 0.8271$.

  • Regimen Conectado ($\nu_0 \le \nu \le 1$): La frontera de fase perturbativa de segundo orden (que separa las fases sin masa y masiva) termina en un punto tricrítico, más allá del cual comienza una línea de primer orden. Esta estructura es cualitativamente similar al modelo GN estándar de un solo componente ($\nu=1$).
  • Regimen Desconectado ($\nu < \nu_0$): El punto tricrítico desaparece. La línea perturbativa de segundo orden se extiende indefinidamente sin terminar. En su lugar, una línea de primer orden emerge dentro de la fase masiva, terminando en un punto final crítico que se encuentra por debajo de la frontera perturbativa.

• Diagrama de Fases Completo con Inhomogeneidades:

  • Para $\nu \ge \nu_0$: La inclusión de fases inhomogéneas reemplaza la línea de primer orden de la aproximación homogénea con una estructura de "cúspide" formada por dos líneas de segundo orden. Estas líneas encierran una fase de cristal. Una línea separa la fase sin masa del cristal, y la otra separa la fase masiva del cristal. El punto tricrítico del cálculo homogéneo sigue siendo una característica del diagrama completo.
  • Para $\nu < \nu_0$: La línea de primer orden del cálculo homogéneo es reemplazada por dos líneas de segundo orden que emanan del punto final crítico, encerrando una región de cristal. A medida que $\nu$ disminuye, el tamaño de esta región inhomogénea (cristal) se reduce rápidamente.
  • Nueva Bifurcación (en $\nu = 3/4$): En la región desconectada, surge una característica novedosa donde la línea de inestabilidad "horizontal" (masiva $\to$ cristal) intersecta la frontera homogénea perturbativa en un nuevo punto $Q$. Esto crea una estructura compleja que involucra cuatro fronteras distintas de segundo orden: sin masa-masiva, sin masa-cristal, masiva-cristal y masiva-solitón único.

• Parámetros Críticos: El artículo proporciona coordenadas numéricas específicas para puntos críticos, puntos tricríticos y puntos finales para varios valores de $\nu$ (por ejemplo, para $\nu=1$, el punto tricrítico está en $\mu \approx 0.6082, T \approx 0.3183$; para $\nu_0$, la bifurcación ocurre en $\mu \approx 0.7202, T \approx 0.2240$).

Significancia y Reivindicaciones
El artículo afirma llenar un vacío en la comprensión del diagrama de fases a temperatura finita del modelo GN con llenado parcial de sabores. Al utilizar análisis de estabilidad en lugar de cálculos completos de HF, los autores demuestran que el diagrama de fases completo, incluyendo las fases inhomogéneas, puede construirse de manera eficiente.

La principal significancia radica en la identificación de la bifurcación en $\nu_0 \approx 0.8271$, que marca un cambio cualitativo en la topología del diagrama de fases. El estudio confirma que el concepto de "bandas de valencia parcialmente llenas" (generalizado del barión DHN a $T=0$) conduce a una rica estructura de fases donde la interacción entre la fracción de llenado $\nu$ y la temperatura determina si el sistema exhibe fronteras de fase conectadas o desconectadas. Los resultados resaltan que la libertad adicional proporcionada por el parámetro $\nu$ enriquece significativamente el diagrama de fases, introduciendo nuevos puntos de bifurcación y modificando la estabilidad de las fases cristalinas, de forma análoga pero distinta al modelo GN masivo.