Hyperstatistical thermodynamics of the one-dimensional Klein-Gordon and Dirac oscillators: a closed-form q-generalized Boltzmann factor and a quantitative comparison with Beck's superstatistics

Resumen Técnico: Termodinámica Hiperestadística de Osciladores de Klein–Gordon y Dirac 1D

Planteamiento del Problema
El artículo aborda la descripción termodinámica de osciladores relativistas unidimensionales —el oscilador de Klein–Gordon (KGO) y el oscilador de Dirac (DO)— cuando están embebidos en entornos fuera del equilibrio donde los parámetros intensivos (como la temperatura inversa $\beta$) fluctúan. En tales escenarios, el factor de Boltzmann estándar $e^{-\beta E}$ es insuficiente. El estudio compara dos marcos estadísticos generalizados:

1. Superestadística de Beck: Un enfoque asintótico donde el factor de Boltzmann generalizado es la transformada de Laplace de una función de densidad de probabilidad (PDF) $f(\beta)$. En la práctica, esto se aproxima frecuentemente mediante una expansión polinómica de baja energía.
2. Hiperestadística: Un marco propuesto por Squillante, Soares, Tsallis y de Souza, que asume una distribución $\gamma$ de los factores de Boltzmann dentro de dominios del sistema, lo que conduce a un peso de exponencial-$q$ de forma cerrada y exacta.

El problema central es determinar qué marco proporciona una descripción más controlada, físicamente consistente y numéricamente estable de la termodinámica de los osciladores relativistas cuando el factor de Boltzmann estándar es reemplazado por pesos estadísticos no boltzmannianos.

Metodología
Los autores realizan una comparación cuantitativa entre los dos marcos utilizando las mismas entradas físicas:

• Modelos: El KGO 1D y el DO 1D. El espectro del DO incluye la degeneración inducida por el espín ($g_0=1, g_n=2$ para $n \ge 1$), mientras que el KGO tiene degeneración unitaria ($g_n=1$).
• Definición de Energía: Para asegurar el cumplimiento de la tercera ley de la termodinámica ($C_v \to 0$ cuando $T \to 0$), los autores utilizan energías de excitación $\varepsilon_n = E_n - E_0$, eliminando el desfase de la energía en reposo.
• Marco de Beck: Los autores reformulan la expansión asintótica de la superestadística de Beck (Ec. 15) utilizando un corchete polinómico $\Pi_q(\langle \beta \rangle \varepsilon)$. Evalúan esto para tres PDFs comunes: Gamma, Log-Normal y F-distribución, con coeficientes $g(q)$ específicos para cada una.
• Marco de la Hiperestadística: Los autores aplican el factor de Boltzmann efectivo de forma cerrada $B_q(\varepsilon) = \exp_q(-\langle \beta \rangle \varepsilon)$, donde $\exp_q$ es la función exponencial-$q$. Esta forma se deriva de una construcción promediada en el dominio donde la transformación $\gamma$ interna produce exactamente la exponencial-$q$, independientemente de la PDF específica utilizada para el ensamble de dominios (siempre que el $\langle \beta \rangle$ esté fijo).
• Cálculo Termodinámico: Las funciones de partición ($Z$), la entropía ($S$) y el calor específico ($C_v$) se computan numéricamente mediante la suma de 8000 niveles de energía. Para la hiperestadística, se utilizan derivadas analíticas de la exponencial-$q$ para calcular $U$ y $C_v$ sin recurrir a la diferenciación numérica. Para el marco de Beck, las derivadas se toman de la expansión polinómica truncada, aplicando cortes de positividad donde el polinomio se vuelve negativo.

Contribuciones Clave

1. Peso de Forma Cerrada de la Hiperestadística: El artículo demuestra que la hiperestadística produce un peso exponencial-$q$ ($B_q$) universal, de forma cerrada, que permanece positivo y decreciente monótonamente para todas las energías, a diferencia de la expansión polinómica asintótica de la superestadística de Beck.
2. Comparación Cuantitativa: El estudio proporciona una comparación numérica directa de la entropía y el calor específico tanto para el KGO como para el DO bajo ambos marcos a través de un rango de índices entrópicos ($q$) y temperaturas.
3. Identificación de la Ruptura Asintótica: Los autores identifican rigurosamente los regímenes de parámetros donde la expansión polinómica de Beck falla. Muestran que para $\langle \beta \rangle \varepsilon$ moderados y $q \gtrsim 1.2$, los corchetes polinómicos para las distribuciones Gamma y Log-Normal se vuelven negativos, mientras que el corchete de la distribución F crece de forma infundada, violando la interpretación de $B(E)$ como un peso estadístico.
4. Efectos de Degeneración: El análisis aísla el impacto termodinámico de la degeneración inducida por el espín del oscilador de Dirac, mostrando que esta eleva la entropía por $\ln 2$ en el límite de alta temperatura y modifica el perfil del calor específico comparado con el KGO.

Resultos

• Regímenes de Validez: Ambos marcos coinciden cuantitativamente en el límite de baja energía y alta temperatura ($q-1 \ll 1$ y $\langle \beta \rangle \varepsilon \lesssim 2$).
• Estabilidad: La hiperestadística produce curvas suaves y monótonas para la entropía y el calor específico en todo el rango numérico considerado ($q \le 1.40$). En contraste, la expansión asintótica de Beck exhibe oscilaciones espurias, valores negativos o crecimiento infundado cuando el corchete polinómico pierde la positividad.
• Límites Termodinámicos: El uso de energías de excitación asegura que $C_v \to 0$ cuando $T \to 0$ para ambos modelos en ambos marcos (dentro del dominio de convergencia $q < 3/2$ para la hiperestadística).
• Diferencias entre Osciladores: El DO exhibe una entropía mayor que el KGO debido a la degeneración doble de los estados excitados. El calor específico del DO se aproxima a $2k_B$ a altas temperaturas, consistente con la estadística de Boltzmann estándar, pero las curvas hiperestadísticas muestran una transición más suave que la expansión truncada de Beck.

Significancia y Reivindicaciones
El artículo sostiene que la hiperestadística es el marco de trabajo más robusto para las sumas termodinámicas de osciladores relativistas en ausencia de una distribución microscópica conocida de las variables intensivas.

• Ventajas: Ofrece una solución de forma cerrada e independiente de la PDF que evita las inestabilidades numéricas (pesos negativos, divergencia) inherentes a la expansión polinómica truncada de la superestadística de Beck.
• Limitaciones del Enfoque de Beck: Los autores argumentan que, si bien la superestadística de Beck es físicamente directa cuando la distribución completa $f(\beta)$ es conocida y se puede computar la transformada de Laplace completa, la expansión asintótica truncada comúnmente utilizada en aplicaciones de osciladores está estructuralmente limitada. Es fiable solo para $\langle \beta \rangle \varepsilon$ pequeños y $q-1$ pequeños.
• Conclusión: Para el problema específico de la termodinámica de los osciladores 1D KGO y DO, la hiperestadística proporciona una descripción unificada, analíticamente tratable y numéricamente estable. El artículo no pretende que la hiperestadística sea universalmente superior en todos los contextos físicos, sino que es la herramienta práctica superior para estas sumas termodinámicas específicas donde la distribución microscópica no se proporciona explícitamente.