⚛️ phenomenology

Order-reversed Kubo formulas in relativistic kinetic theory

This paper verifies that the stress-energy tensor response functions derived from the massive Anderson-Witting kinetic model satisfy necessary analytic conditions and are fully consistent with recently established order-reversed Kubo formulas.

Original authors: Sangyong Jeon, Juhee Hong, Alina Czajka

Published 2026-07-20

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Sangyong Jeon, Juhee Hong, Alina Czajka

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Imagine the universe as a giant, invisible soup. Sometimes this soup is so hot and dense that it behaves like a perfect fluid, flowing without any friction, like water sliding over ice. Other times, it's a bit more like honey, thick and sticky, resisting that flow. Scientists call this resistance "viscosity," and understanding exactly how much the universe resists moving is crucial for figuring out how the Big Bang evolved into the stars and galaxies we see today. To measure this, physicists use special mathematical recipes called "Kubo formulas." Think of these formulas as a set of instructions for a chef: if you know how the soup reacts when you poke it, you can calculate exactly how sticky it is. For a long time, scientists had one main way to poke the soup: they would wait until the poke was completely still (zero frequency) and then check how the soup settled (zero wavenumber). But recently, a new group of scientists discovered that if you change the order of these steps—checking the stillness first, then the settling—you get a different, more accurate recipe. This new method is tricky because it requires the mathematical "soup" to have specific, hidden features that only appear when you look at it very closely.

In this paper, the authors, Sangyong Jeon, Juhee Hong, and Alina Czajka, decide to test these new, tricky recipes. They don't use a real soup or a supercomputer simulation of the entire universe; instead, they use a simplified model called "kinetic theory," which treats the particles in the soup like a crowd of billiard balls bouncing around. They specifically look at a scenario where these "balls" have a little bit of weight (mass), which makes the math much harder than if they were weightless. Their goal was to see if the mathematical results from their billiard-ball model actually fit the new, complex rules of the reversed-order recipes. They found that, indeed, the model works perfectly. The "billiard balls" behave exactly as the new formulas predict, showing that these new mathematical tools are reliable and consistent. This is a big deal because it proves that the new recipes aren't just theoretical guesses; they hold up even when the physics gets complicated with mass. It's like verifying that a new, complex map of a city works perfectly even when you add traffic jams and construction zones, giving scientists confidence that they can use these new tools to explore the deepest secrets of the universe's fluid behavior.

The core of the paper involves checking eight specific mathematical formulas. In the old way of doing things, scientists would take a limit where the frequency of a disturbance goes to zero first, and then the wavelength goes to zero. The new formulas flip this order: they take the wavelength to zero first, then the frequency. The authors showed that their kinetic theory model, which includes particles with mass, produces correlation functions (mathematical descriptions of how the soup reacts) that have the exact right "poles" or peaks in the right places to satisfy these new formulas. They calculated the shear viscosity (how the soup resists sliding layers past each other) and the bulk viscosity (how it resists being squeezed) using these new methods.

The results were consistent. When they applied the new formulas to their model, they got the same values for the viscosity as they would have expected from other established methods, but with the added benefit of using the new, reversed-order limits. Specifically, they found that for a gas of particles with a small mass compared to the temperature (represented by the ratio $\xi = m/T$ ), the shear viscosity $\eta$ and bulk viscosity $\zeta$ follow specific patterns. For example, the ratio of shear viscosity to enthalpy density ( $\eta/h_0$ ) comes out to be $\tau_R (1/5 - 1/60\xi^2 + 1/96\xi^4) + O(\xi^5)$ , and the bulk viscosity ratio $\zeta/h_0$ is $5/432\tau_R\xi^4 + O(\xi^5)$ . These numbers matched up with previous studies, confirming that the new formulas are robust.

The authors also explored a subtle detail about "relaxing modes," which are like the different ways the soup can wiggle or settle down after being disturbed. They noticed a relationship between the speed of sound, the relaxation time, and the number of these modes, suggesting that the quantity $n_R = v_s^2/(\tau_R R_R)$ might represent the total number of relaxing modes. However, they are careful to note that while their model suggests there might be six such modes, they don't know for sure if there are exactly six or if this is a universal rule for all systems. They also point out that while their model works well, it relies on a specific approximation (the Anderson-Witting model) and doesn't yet account for particles whose mass changes depending on where they are in space and time, which is a much harder problem to solve.

Ultimately, this paper serves as a crucial "stress test" for the new Kubo formulas. By showing that a realistic kinetic theory model with mass satisfies these formulas, the authors provide strong evidence that the new mathematical approach is valid. They emphasize that these formulas rely on the deep structural consistency of how energy and momentum correlate in a fluid, rather than just being a static snapshot. While this makes the formulas powerful for theoretical physics, the authors also note that this complexity means they might be very difficult to implement directly in computer simulations used on supercomputers (lattice QCD). Nevertheless, this work bridges the gap between simple particle models and complex field theories, giving physicists a new, reliable way to calculate how the universe's primordial soup flows.

Technical Summary: Order-reversed Kubo formulas in relativistic kinetic theory

Problem Statement
Recent work (Ref. [1]) introduced eight new Kubo formulas for shear ( $\eta$ ) and bulk ( $\zeta$ ) viscosities. Unlike standard transport coefficient calculations in quantum field theory, which typically take the zero frequency ( $\omega \to 0$ ) and zero wavenumber ( $k \to 0$ ) limits simultaneously (or in an order that avoids hydrodynamic poles), these new formulas require a specific ordering of limits: taking $\omega \to 0$ first, followed by $k \to 0$ . This order-reversed limit is necessary to access specific analytic structures of stress-energy tensor correlation functions that exhibit hydrodynamic poles.

A critical challenge arises because standard perturbative methods in quantum field theory (such as the Chapman-Enskog expansion) often commute these limits without encountering poles, rendering them insufficient for deriving these specific formulas. Furthermore, while linearized kinetic theory in the relaxation time approximation (RTA) has been shown to reproduce ladder diagram resummations in field theory, it remained to be verified whether the correlation functions derived from kinetic theory with non-zero particle mass satisfy the rigorous analytic conditions required by the new Kubo formulas.

Methodology
The authors utilize the Anderson-Witting model of relativistic kinetic theory with a non-zero particle mass ( $m \neq 0$ ) in the high-temperature limit ( $m/T \ll 1$ ). The methodology proceeds as follows:

Linearized Kinetic Equation: The Boltzmann equation is linearized under a small metric perturbation $h_{\mu\nu}$ in a flat background. The authors employ Maxwell-Boltzmann statistics and the Landau matching condition to define local equilibrium.
Stress-Energy Tensor Correlators: The deviation of the distribution function ( $\delta f$ ) is solved in Fourier space. The stress-energy tensor deviation ( $\delta T^{\mu\nu}$ ) is calculated by integrating $\delta f$ over momentum space.
Response Functions: The authors derive explicit expressions for the stress-energy tensor correlation functions. These are decomposed into five independent response functions ( $G_L, G_{LT}, G_T, G_1, G_2$ ) which depend on frequency $\omega$ and wavenumber $k$ .
Analytic Verification: The authors analyze the small $\omega$ and small $k$ limits of these derived response functions. They verify that these functions satisfy the specific analytic constraints (pole structures and limit behaviors) required by the eight new Kubo formulas introduced in Ref. [1].
Limit Order Comparison: The study explicitly compares the results of taking the limits in the standard order ( $k \to 0$ then $\omega \to 0$ ) versus the order-reversed sequence ( $\omega \to 0$ then $k \to 0$ ) to demonstrate how different Kubo formulas emerge from different limit orderings.

Key Contributions and Results

Verification of Analytic Structures: The paper demonstrates that the stress-energy tensor response functions derived from the massive Anderson-Witting model satisfy the required analytic conditions. Specifically, they exhibit the necessary hydrodynamic poles in the small $\omega$ and $k$ limits, confirming their suitability for evaluating the new Kubo formulas.
Consistency of Kubo Formulas: By applying the order-reversed limits to the kinetic theory results, the authors verify that all eight new Kubo formulas yield consistent values for the transport coefficients.
- For the massless case ( $m=0$ ), the results recover the known vanishing bulk viscosity ( $\zeta = 0$ ), the speed of sound $v_s^2 = 1/3$ , and the specific shear viscosity $\eta/h_0 = \tau_R/5$ .
- For the massive case ( $m \neq 0$ ) in the limit $m/T \ll 1$ , the authors derive expansions for the specific viscosities in terms of $\xi = m/T$ :
  $\frac{\eta}{h_0} = \tau_R \left( \frac{1}{5} - \frac{1}{60}\xi^2 + \frac{1}{96}\xi^4 \right) + O(\xi^5)$
  $\frac{\zeta}{h_0} = \frac{5}{432\tau_R}\xi^4 + O(\xi^5)$
- These results agree with previous calculations in Ref. [20] and numerical results from Ref. [22].
Limit Order Dependence: The study clarifies how different Kubo formulas correspond to different limit orderings. For instance, taking $\omega \to 0$ first yields the eighth Kubo formula for $(\zeta + 4\eta/3)$ , whereas taking $k \to 0$ first yields the standard shear viscosity formula.
Relaxation Modes: The authors observe that the kinetic theory correlators imply a relationship $v_s^2 R_R = 6\tau_R$ (where $R_R$ is the sum of relaxation times). They speculate that the quantity $n_R = v_s^2 / (\tau_R R_R)$ may represent the total number of relaxing modes, noting that massive second-order hydrodynamics has 2 modes and massless third-order has 3, while the kinetic theory result suggests 6.

Significance and Claims
The authors state that the primary significance of this work is twofold:

Validation of New Formulas: It provides a concrete verification of the eight new Kubo formulas within a simple, non-trivial setting (kinetic theory with non-zero mass), demonstrating their consistency and veracity.
Bridge to Field Theory: It confirms that the analytic structures of stress-energy correlation functions derived in Ref. [1] are realized in linearized kinetic theory. Since ladder diagram resummations in quantum field theory can be re-organized to recover linearized kinetic theory, this work serves as an essential intermediate consistency check for future full field theory calculations.

The paper concludes that the new Kubo formulas are robust tools for extracting transport coefficients, provided the correlation functions possess the correct analytic structure (hydrodynamic poles) and the limits are taken in the correct order. The authors note that direct lattice implementation of these formulas would be difficult due to the frequency derivatives and specific limit ordering required. Future work is suggested to include space-time dependent quasi-particle masses and to extend these calculations to simple field theories via effective Boltzmann equations.

More like this