On electric fields in hot QCD: infrared regularization dependence

This paper resolves a historical discrepancy in the definition of electric susceptibility in hot QCD by utilizing an exact fermion propagator and improved perturbation theory to demonstrate that the disagreement stems from the choice of thermodynamic ensemble and infrared regularization, while also constructing the susceptibility within a hadron resonance gas model for the low-temperature regime.

The Gist

Imagine you are holding a giant, invisible sponge. This sponge is made of a hot, bubbling soup of tiny, charged particles (like electrons or quarks). Now, imagine you turn on a powerful electric fan and blow a steady stream of air (an electric field) through this sponge.

What happens? The sponge reacts. But here's the twist: scientists have been arguing for years about exactly how to measure that reaction.

This paper, titled "On electric fields in hot QCD," is the story of how two teams of physicists finally solved a long-standing mystery about why their measurements didn't match. They discovered that the disagreement wasn't because one team was wrong, but because they were asking slightly different questions about the nature of the "sponge."

Here is the breakdown in simple terms:

1. The Setup: The Hot Soup and the Fan

In the early moments of the universe (or in particle smashers like the Large Hadron Collider), matter exists as a super-hot plasma.

• The Players: Charged particles floating in a hot soup.
• The Disturbance: A background electric field (like the fan).
• The Goal: Measure the Electric Susceptibility. Think of this as a score that tells us how easily the soup gets "polarized" or how much it resists the fan's wind.

2. The Problem: Two Different Scores

For a long time, physicists used two different mathematical recipes to calculate this score.

• Recipe A (The "Schwinger" method): This method looks at the exact, complex path every single particle takes in the electric field.
• Recipe B (The "Weldon" method): This method uses a simpler, averaged approach, looking at how the particles interact with the field in a "cloud."

The Mystery: When they applied these recipes to hot plasma, they got different numbers.

• Recipe A said: "The soup reacts this way."
• Recipe B said: "No, it reacts that way."
  It was like two weather forecasters predicting the temperature, but one said 70°F and the other said 80°F, and neither could figure out why.

3. The Culprit: The "Infinite Room" vs. The "Oscillating Wind"

The authors of this paper realized the problem wasn't the math; it was the setup. They identified two hidden traps in how we imagine the experiment:

Trap 1: The Infinite Room (Infrared Divergence)

Imagine trying to fill an infinitely large room with a perfectly uniform wind blowing from one side to the other.

• In a real, finite room, the wind pushes the air, but the walls stop it, and the air piles up at the back.
• In an infinite room with a uniform wind, the air would have to pile up forever. The density of particles would become infinite at the edges. This creates a mathematical "singularity" (a breakdown).
• The Fix: You can't actually have an infinite room with a uniform wind. You have to pretend the room has a size (a finite volume) or pretend the wind isn't perfectly uniform.

Trap 2: The Order of Operations (The "Non-Commutative" Trap)

This is the most important part. The authors showed that the order in which you do things matters.

Imagine you are trying to measure the average wind speed in a room.

• Scenario A: You measure the wind speed at every single point in the room first, and then you average them up.
• Scenario B: You pretend the wind is a gentle, wavy breeze (oscillating field), measure the average, and then make the waves infinitely long (turning it into a steady wind).

The paper proves that Scenario A and Scenario B give different results.

• If you average the room first and then make the wind uniform, you get Result X (Recipe A).
• If you make the wind uniform first and then average, you get Result Y (Recipe B).

In the real world, these correspond to two different physical situations:

1. The "Grand Canonical" Ensemble: The system can swap particles with the outside world (like a room with an open door). This leads to Result X.
2. The "Canonical" Ensemble: The system is closed; the number of particles is fixed (like a sealed room). This leads to Result Y.

4. The Analogy: The Crowd in a Hallway

Let's use a crowd of people in a hallway to visualize this.

• The Electric Field: A strong wind blowing down the hallway.
• The Particles: People walking.

The "Schwinger" View (Recipe A):
You look at the hallway as a giant, open space. The wind pushes people, and they pile up at the far end. You calculate the pressure based on this massive pile-up. This is what happens if you let people enter and leave the hallway freely.

The "Weldon" View (Recipe B):
You look at a specific, sealed section of the hallway. The wind pushes people, but they can't leave the section. Instead of piling up infinitely, they just shift slightly within the walls. You calculate the pressure based on this shift.

The Discovery:
The authors showed that the "pile-up" (the macroscopic charge rearrangement) is a real physical effect, but it's not the "quantum response" (the microscopic jiggling of particles) that we are usually interested in.

• If you want to know how the material itself reacts to the wind (ignoring the pile-up), you must use the Weldon method (Recipe B).
• If you want to know the total effect including the pile-up, you use the Schwinger method (Recipe A).

The disagreement happened because scientists were mixing these two definitions without realizing it.

5. The Solution: A New Model

To prove their point, the authors built a simplified model using "Hadron Resonance Gas" (imagine the soup is made of simple, heavy balls like pions instead of complex quarks).

• They calculated the susceptibility using their new, clarified rules.
• They compared their math to real-world data from Lattice QCD (supercomputer simulations of the strong force).
• The Result: Their new, clarified math matched the supercomputer data perfectly!

The Takeaway

The paper teaches us that in physics, how you define your experiment is just as important as the math you use.

• The "Infinite Volume" limit (pretending the universe is endless) and the "Homogeneous Field" limit (pretending the wind is perfectly steady) don't always play nice together.
• Depending on whether you fix the number of particles or let them flow, and whether you average the field first or the space first, you get different answers.
• Both answers are "correct," but they describe different physical realities.

In short: The universe isn't confusing; our definitions were just a bit fuzzy. By tightening the definitions, the authors resolved the conflict and showed us exactly how hot, charged matter behaves when you blow an electric fan on it.

Technical Summary

1. Problem Statement

The paper addresses a long-standing discrepancy in thermal field theory regarding the calculation of the electric susceptibility ($\xi$) of a hot plasma of charged particles in the presence of a background electric field. The susceptibility, defined as the second derivative of the thermodynamic potential $\Omega$ with respect to the electric field ($E$) at $E=0$, characterizes the medium's linear response.

Historically, two distinct formalisms have yielded different results for $\xi$ at non-zero temperature ($T$):

1. Schwinger's Approach ($\xi_S$): Uses the exact fermion propagator in a homogeneous electric background (often via the imaginary time formalism or real-time methods).
2. Weldon's Approach ($\xi_W$): Uses a gradient expansion of the effective action or photon vacuum polarization diagrams at $E=0$.

The results differ by a constant term in the high-temperature limit:
$$\xi_S - \xi_W = \frac{1}{6\pi^2} + \mathcal{O}(m^2/T^2)$$
Previous literature conjectured this mismatch arose from infrared (IR) divergences caused by the singular equilibrium charge profile in an infinite volume under a homogeneous electric field. However, the precise origin of the disagreement—specifically regarding the ordering of limits and the choice of thermodynamic ensemble—remained unresolved.

2. Methodology

The authors employ a rigorous perturbative approach at the one-loop level to resolve the discrepancy. Their methodology involves:

• Finite Volume and IR Regularization: Instead of working immediately in the infinite volume limit, the authors introduce IR regularization by considering a system of finite length $L_3$ in the direction of the electric field. They also consider an alternative regularization using an oscillatory (sinusoidal) electric field with wavenumber $k$.
• Exact Propagator Derivation: A key technical achievement is the derivation of the exact fermion propagator in a finite volume ($L_3$) at non-zero temperature ($T$) and in the presence of a homogeneous (imaginary) electric field. This involves twisted boundary conditions to handle the gauge potential.
• Thermodynamic Ensembles: The authors distinguish between two ensembles to isolate the "quantum response" from the "macroscopic charge rearrangement":
  • Local Grand Canonical Ensemble: Allows particle exchange; the effective chemical potential $\bar{\mu}$ is fixed to a constant (zero for a neutral system), effectively canceling the macroscopic charge redistribution.
  • Local Canonical Ensemble: Fixes the charge density profile; the chemical potential is adjusted locally to maintain a fixed density.
• Limit Ordering Analysis: The core of the analysis involves examining the non-commutativity of limits:
  1. Infinite volume limit ($L_3 \to \infty$) vs. Weak field limit ($E \to 0$).
  2. Infinite wavelength limit ($k \to 0$) vs. Spatial averaging ($\int dx_3$).

3. Key Contributions

A. Resolution of the $\xi_S \neq \xi_W$ Discrepancy

The authors demonstrate that the mismatch is not due to a fundamental error in either method but arises from different physical definitions of the observable based on how IR limits are taken:

• Homogeneous Fields (Finite Volume): When taking the infinite volume limit ($L_3 \to \infty$) before or after spatial averaging in a homogeneous field, the result is consistently $\xi_S$ in both canonical and grand canonical ensembles.
• Oscillatory Fields: When using an oscillatory field ($A_0 \sim \sin(kx)$) and taking the limit $k \to 0$:
  • If spatial averaging is performed before the $k \to 0$ limit, the result is $\xi_W$.
  • If the $k \to 0$ limit is taken before spatial averaging, the result is $\xi_S$.
• Ensemble Dependence: In the canonical ensemble with oscillatory fields, the limits $k \to 0$ and spatial averaging do not commute. This non-commutativity is the root cause of the difference between $\xi_S$ and $\xi_W$. The grand canonical ensemble yields $\xi_S$ regardless of the ordering in the oscillatory case, while the canonical ensemble yields $\xi_W$ if averaging precedes the limit.

B. Physical Interpretation of IR Divergences

The paper clarifies that the "singular equilibrium charge profile" in an infinite homogeneous field creates an IR divergence. The two methods effectively measure different things:

• $\xi_S$ corresponds to the response of a system where the macroscopic charge redistribution is suppressed (or the system is treated as locally neutral with fixed density).
• $\xi_W$ corresponds to a setup where the spatial averaging of the oscillatory field includes the contribution of the inhomogeneous charge profile before the field becomes homogeneous.

C. Hadron Resonance Gas (HRG) Model

The authors construct an electric susceptibility model for the low-temperature QCD regime using a Hadron Resonance Gas consisting of charged pions (scalar fields).

• They calculate $\xi$ using the Weldon-type definition (oscillatory fields, canonical ensemble, averaging before $k \to 0$).
• This choice is motivated by the fact that Lattice QCD simulations typically measure susceptibility via current-current correlators at fixed mass, which corresponds to the $\xi_W$ definition.
• The model predicts a negative electric susceptibility ($\xi < 0$), indicating that thermal fluctuations generate a polarization opposing the external field.

4. Results

• Table 1 Summary: The paper provides a comprehensive table (Table 1) mapping the resulting susceptibility values based on the ensemble (Canonical vs. Grand Canonical) and the order of limits (Volume averaging vs. IR regularization removal).
  • Grand Canonical: Yields $\xi_S$ in all limit orderings.
  • Canonical: Yields $\xi_S$ for finite volume limits but yields $\xi_W$ if the oscillatory field limit ($k \to 0$) is taken after spatial averaging.
• Analytic Expressions:
  • $\xi_S = -\frac{1}{6\pi^2} \int_m^\infty d\omega \frac{1}{\sqrt{\omega^2-m^2}} [2n_F(\omega) + \omega n'_F(\omega)]$
  • $\xi_W = -\frac{1}{6\pi^2} \int_m^\infty d\omega \frac{1}{\sqrt{\omega^2-m^2}} [2n_F(\omega) - \omega n'_F(\omega)]$
• Lattice Comparison: The HRG model results for both electric ($\xi$) and magnetic ($\chi$) susceptibilities show excellent agreement with continuum-extrapolated Lattice QCD data up to temperatures of $T \approx 120$ MeV.
  • The magnetic susceptibility is negative (diamagnetic), consistent with the spinless nature of pions.
  • The electric susceptibility decreases (becomes more negative) as temperature increases.

5. Significance

1. Theoretical Consistency: The paper resolves a decades-old ambiguity in thermal field theory, proving that $\xi_S$ and $\xi_W$ are not contradictory but correspond to different measurement protocols (different orderings of limits and ensemble choices).
2. Experimental Relevance: The findings are crucial for interpreting:
   • Heavy-Ion Collisions: Understanding the early-stage QCD plasma response to strong electromagnetic fields generated by colliding nuclei.
   • Laser Physics: Interpreting the response of QED plasmas to intense laser pulses.
3. Lattice QCD Connection: By identifying that Lattice QCD simulations effectively probe the $\xi_W$ definition (via current correlators in a fixed volume), the authors provide a rigorous theoretical framework to compare perturbative calculations with non-perturbative lattice results.
4. Methodological Rigor: The derivation of the exact finite-volume propagator and the detailed analysis of IR regularization limits set a new standard for handling singularities in background field theories.

In conclusion, the paper establishes that the "disagreement" in electric susceptibility is a manifestation of the non-commutativity of infrared limits in the canonical ensemble, and it provides a unified framework to reconcile perturbative QCD/QED calculations with lattice simulations and experimental phenomenology.