Beth-Uhlenbeck equation for the thermodynamics of fluctuations in a generalised 2+1D Gross-Neveu model

Resumen Técnico: Ecuación de Beth-Uhlenbeck para Fluctuaciones en un Modelo de Gross-Neveu Generalizado en (2+1)D

Planteamiento del Problema
El artículo investiga la termodinámica de un modelo de Gross-Neveu (GN) generalizado en (2+1) dimensiones, un sistema motivado por la física del grafeno y su relación de dispersión lineal cerca de los puntos de Dirac. Aunque la estructura de fase y las propiedades termodinámicas de este modelo han sido estudiadas previamente dentro de la Aproximación de Campo Medio (MFA, por sus siglas en inglés), tales enfoques omiten los efectos de las fluctuaciones de los campos auxiliares. Los intentos previos de ir más allá del campo medio en (3+1) dimensiones utilizaron el enfoque de Beth-Uhlenbeck, pero carecían de implementación numérica y de la dependencia del momento en los desplazamientos de fase. Este trabajo aborda la brecha al explorar numéricamente el modelo en (2+1)D más allá del campo medio, centrándose específicamente en la dependencia del momento de las funciones de polarización y los desplazamientos de fase, los cuales son cruciales para calcular con precisión las contribuciones de las fluctuaciones a las cantidades termodinámicas.

Metodología
Los autores emplean un modelo de GN generalizado con interacciones de cuatro fermiones y acoplamientos iguales ($G_i = G$), respetando la simetría SU(2). La metodología procede en dos etapas:

1. Aproximación de Campo Medio (MFA): Los autores establecen primero la línea base derivando las ecuaciones de brecha utilizando la transformación de Hubbard-Stratonovich. Introducen campos auxiliares para desacoplar la interacción de cuatro fermiones, integran los grados de libertad fermiónicos y renormalizan la teoría utilizando un corte de momento $\Lambda$. La ecuación de brecha resultante determina la masa del condensado de fermiones $m$ como una función de la temperatura $T$ y el potencial químico $\mu$.
2. Más allá del Campo Medio (Fluctuaciones Gaussianas): Para incluir las fluctuaciones, los autores expanden los campos auxiliares alrededor de sus valores de campo medio ($\bar{\Phi}_i \to \bar{\Phi}_i + \delta\Phi_i$). Derivan la contribución de la fluctuación al potencial grandioso ($\Omega_{fl}$) mediante una integración gaussiana, la cual relaciona $\Omega_{fl}$ con el determinante de los propagadores inversos de los campos escalar y pseudoscalar.
   • Función de Polarización: El núcleo del cálculo consiste en evaluar la función de polarización $\Pi_i(q, \omega)$, incluyendo su completa dependencia de momento ($q$) y frecuencia ($\omega$). Los autores derivan expresiones analíticas para la parte imaginaria de $\Pi_i$ utilizando la relación de Sokhotski-Plemelj y funciones delta de Dirac, y la parte real mediante relaciones de Kramers-Kronig.
   • Enfoque de Beth-Uhlenbeck: El propagador inverso se descompone para extraer los desplazamientos de fase ($\phi_i$), separándolos en partes resonantes (estado ligado) y de dispersión (fondo). La presión de fluctuación se calcula entonces utilizando la fórmula generalizada de Beth-Uhlenbeck, que integra los desplazamientos de fase ponderados por la distribución de Bose-Einstein.
   • Dependencia del Momento: Un paso metodológico crítico es la inclusión explícita de los efectos del momento finito en la función de polarización y la brecha de masa, contrastando con las aproximaciones previas de "impulso de Lorentz" donde los desplazamientos de fase a momento finito se derivaban únicamente mediante el impulso de los resultados en el reposo ($q=0$).

Contribuciones Clave y Resultados

• Implementación Numérica: El artículo proporciona la primera exploración numérica del modelo GN en (2+1)D más allá del campo medio, incluyendo el cálculo de masas excitónicas, anchos de decaimiento y desplazamientos de fase con dependencia completa del momento.
• Estados Ligados Excitónicos: Los autores calculan las masas de los estados ligados escalares y pseudoscalares. Confirman que los estados pseudoscalares (bosones de Goldstone en el límite quiral) son estables a bajas temperaturas, pero experimentan una disociación de Mott (se vuelven no ligados) cuando su masa supera el doble de la masa del fermión constituyente ($2m$) en una temperatura de Mott específica ($T_{Mott} \approx 0.955 M$).
• Desplazamientos de Fase Dependientes del Momento:
  • Amortiguamiento de Landau: La inclusión de un momento finito revela desplazamientos de fase no nulos en la región $s = \omega^2 - q^2 < 0$, correspondientes a modos de amortiguamiento de Landau. Estos modos están ausentes en la aproximación de impulso de Lorentz.
  • Resurrección de Estados Ligados: Un hallazgo significativo es la "resurrección" de los estados ligados pseudoscalares a temperaturas por encima de $T_{Mott}$ pero para momentos externos suficientemente grandes ($q > q_{Mott}$). Los autores atribuyen esto a la reducción del bloqueo de Pauli para los pares electrón-hueco que se mueven en relación con el medio. Este fenómeno está ausente en la aproximación estándar de impulso.
• Impacto Termodinámico: El cálculo de la presión de fluctuación ($P_{fl}$) demuestra que el tratamiento dependiente del momento produce resultados significativamente diferentes comparados con la aproximación de impulso de Lorentz. La región de amortiguamiento de Landau y los estados ligados resucitados contribuyen sustancialmente a la presión, lo que conduce a un aumento marcado en la presión de fluctuación respecto a aproximaciones anteriores.
• Renormalizabilidad y Escalamiento: El estudio destaca las diferencias entre las versiones de (2+1)D y (3+1)D del modelo. Los resultados de (2+1)D muestran una menor sensibilidad al corte de momento y exhiben una velocidad del sonido sublumínica que se aproxima correctamente al límite conformal ($c_s^2 = 1/2$) a altas temperaturas.

Significancia y Reivindicaciones
Los autores afirman que su trabajo cierra la brecha entre las formulaciones teóricas del enfoque de Beth-Uhlenbeck y las aplicaciones numéricas en dimensiones (2+1). La principal significancia radica en demostrar que la suposición de invariancia de impulso de Lorentz para los desplazamientos de fase en un medio es inválida para este sistema. El medio define un marco de reposo preferencial, y omitir la dependencia explícita del momento en la función de polarización conduce a predicciones termodinámicas incorrectas.

Específicamente, el artículo asevera que:

1. La presencia de modos de Landau y la resurrección de los estados ligados inducida por el momento son fenómenos físicos críticos que deben incluirse para una descripción precisa de la termodinámica del sistema.
2. Estos efectos resultan en una "diferencia significativa" en la presión de fluctuación calculada en comparación con la literatura previa que dependía de los desplazamientos de fase por impulso.
3. El modelo sirve como una teoría efectiva viable para estudiar las transiciones semimetal-aislante en sistemas tipo grafeno, proporcionando un marco para analizar las fluctuaciones que actualmente son difíciles de observar experimentalmente debido a los efectos de apantallamiento.

El artículo concluye modestamente señalando que, si bien el marco teórico está establecido y explorado numéricamente, la observación experimental directa de estos efectos de fluctuación específicos en el grafeno sigue siendo un desafío debido al desorden y al apantallamiento térmico.