Derivation of Gibbs measure from Gibbs state with the fractional Bessel interaction in Two Dimensions

This paper establishes the derivation of the classical Gibbs measure on $\mathbb{T}^2$ with a fractional Bessel interaction from a renormalized grand-canonical quantum Bose gas for the entire range $\frac{3}{2}<\beta\leq2$, overcoming the non-summability of the potential and divergent self-energy through zero-mode renormalization and detailed high-frequency analysis to prove the convergence of both relative free energies and reduced density matrices.

The Gist

The Big Picture: From Quantum Chaos to Classical Order

Imagine you are trying to understand the weather.

• The Quantum View: You are looking at every single air molecule, tracking its exact position, speed, and quantum spin. It's a chaotic, jittery mess of trillions of particles interacting in complex ways. This is the Quantum Gibbs State.
• The Classical View: You look at the weather map. You see smooth clouds, wind patterns, and temperature gradients. You don't care about individual molecules; you care about the "field" of air. This is the Classical Gibbs Measure.

The Goal of the Paper: The authors want to prove that if you take the chaotic quantum system and slowly "turn down the volume" on quantum weirdness (a process called the limit $\lambda \downarrow 0$), the system naturally settles into the smooth, classical weather pattern we expect. They want to show the math that bridges the gap between the jittery micro-world and the smooth macro-world.

The Problem: The "Infinite Noise" Wall

Usually, this bridge is easy to build if the particles interact gently. But in this paper, the particles interact via something called a Fractional Bessel Potential.

The Analogy: Imagine the particles are people at a party.

• Normal Interaction: People talk to their neighbors. If you move away, the conversation fades.
• This Interaction: These people have a super-power. They can hear everyone in the room, no matter how far away they are, and the signal gets stronger the more people there are. In math terms, the interaction is "non-summable."

The Crisis: When the authors tried to do the math, they hit a wall. The "noise" (or self-energy) from these super-connections was infinite. It was like trying to calculate the total volume of a party where everyone is shouting at once; the number explodes to infinity, and the equations break.

The Solution: The "Renormalization" Trick

To fix the infinite noise, the authors had to be very clever. They used a technique called Renormalization.

The Analogy: The "Centered" Party Host
Imagine the party has a host (the "zero mode") who is the loudest person in the room.

1. The Problem: If you try to measure the total noise, the host's voice drowns out everything else and breaks the microphone (the math).
2. The Fix: The authors decided to "center" the host. They subtracted the host's average noise level from the equation. Now, instead of measuring the raw shouting, they measure the fluctuations—how much the host is louder or quieter than their usual self.
3. The Result: The infinite noise disappears, and the math becomes manageable again.

The Strategy: The "Onion Peeling" Method

The paper is a massive engineering project to prove that the quantum system becomes the classical one. They did this by peeling the problem apart like an onion, dealing with different layers of frequency (how fast the particles are jiggling) one by one.

1. The Low-Frequency Core (The Smooth Part)

• What it is: The slow, big movements of the particles.
• The Strategy: They isolated this part and showed it behaves exactly like the classical "weather map" (the Gibbs measure). They used a tool called the Quantum de Finetti Theorem, which is like saying, "If you look at enough particles, their collective behavior looks like a smooth, random field."

2. The High-Frequency Shell (The Jittery Middle)

• What it is: The fast, jittery movements.
• The Problem: This is where the "infinite noise" usually hides.
• The Strategy: They split this layer into two:
  • The Shell: The middle layer of jitters. They proved that even though these particles are jittering, their fluctuations cancel each other out in a way that makes their net effect vanish as the system gets larger.
  • The Tail: The very fastest, highest-energy jitters. They showed that these are so rare and so weak (due to the specific math of the Bessel potential) that they don't matter. They essentially "fade into the background."

The "Double Commutator" Secret Weapon

One of the most technical parts of the paper involves something called a Double Commutator Estimate.

The Analogy: Imagine trying to balance a stack of cards (the particles) on a wobbly table (the quantum fluctuations).

• If you just push the stack, it falls.
• But the authors found a specific way to push and pull (a "commutator") that creates a stabilizing force. They proved that the "wobble" of the quantum system actually creates a hidden rigidity that keeps the stack from collapsing. This allowed them to prove that the high-frequency errors are small enough to ignore.

The Conclusion: The Bridge is Built

After all this hard work—fixing the infinite noise, peeling the onion, and proving the high-frequency jitters vanish—the authors reached their destination.

The Result:
They proved that as you turn down the quantum effects:

1. The Energy Matches: The total energy of the quantum system converges perfectly to the energy of the classical system.
2. The Patterns Match: The probability of finding particles in certain arrangements (the "density matrices") converges to the smooth, classical probability distribution.

In Simple Terms:
They took a system that looked mathematically impossible to solve because of infinite noise, cleaned it up by removing the "average" noise, and proved that underneath all that quantum chaos, there is a beautiful, smooth, classical order waiting to be found. They successfully built the bridge from the quantum world to the classical world for this specific, difficult type of interaction.

Technical Summary

1. Problem Statement and Motivation

The paper addresses the rigorous derivation of a classical nonlinear Gibbs measure from a renormalized grand-canonical quantum Bose gas in two dimensions ($T^2$). Specifically, it focuses on the fractional Bessel interaction potential defined by its Fourier coefficients:
$$\hat{v}_\beta(k) = \langle k \rangle^{-\beta} = (1 + |k|^2)^{-\beta/2}, \quad k \in \mathbb{Z}^2$$
with the exponent range $3/2 < \beta \le 2$.

Key Challenges:

• Non-Summability: In this range, the potential is not summable ($\sum_k \hat{v}_\beta(k) = \infty$). This prevents the standard "density-square" rewriting of the interaction (i.e., $W \sim \int :\rho^2:$) because it leads to a divergent self-energy term proportional to the particle number $N$.
• Ultraviolet Singularity: Unlike previous works dealing with bounded or regular potentials, the singularity here is intrinsic to the quantum Hamiltonian itself. The interaction cannot be regularized simply by removing a constant; the zero mode must be handled differently.
• Limit Regime: The goal is to show that as the temperature parameter $\lambda \downarrow 0$ (high temperature/mean-field limit), the quantum Gibbs state converges to the classical Gibbs measure associated with the renormalized interaction $D(u) = \frac{1}{2} \int :|u|^2: v_\beta :|u|^2:$.

2. Methodology and Strategy

The authors employ a variational approach combined with low-frequency localization and high-frequency remainder analysis. The proof is structured around controlling the relative free energy and the reduced density matrices.

A. Renormalization Scheme

Since the standard density-square form diverges, the authors define the renormalized interaction $W^{re}$ directly in terms of the genuine two-body operator, isolating and centering only the zero mode:
$$W^{re} = W^{(2)}_\beta + c_\lambda N + C_\lambda$$
where $c_\lambda$ and $C_\lambda$ are chosen to cancel the divergent self-energy of the zero mode, specifically using a centered number-fluctuation term $(N - N_0)^2$. The non-zero modes are kept in their original quartic form.

B. Decomposition of Errors

The core of the proof involves decomposing the interaction error into three distinct parts based on frequency cutoffs ($\Lambda = \lambda^{-\alpha}$ and $\Lambda_1 = \lambda^{-\alpha_1}$):

1. Zero-Mode Remainder ($HF_0$): High-frequency fluctuations of the zero mode.
2. Shell Contribution ($E_{shell}$): Interactions where at least one particle is in the "shell" (intermediate frequencies between $\Lambda$ and $\Lambda_1$).
3. Tail Contribution ($E_{tail}$): Interactions involving the "far tail" (frequencies $> \Lambda_1$).

C. Key Technical Tools

• Second-Order Correlation Inequality: Used extensively (Theorem 4.1) to bound fluctuations of observables. This relates the variance of an observable to its double commutator with the Hamiltonian.
• Asymmetric Quartic Hilbert-Schmidt Bounds: Since the interaction is not summable, the authors cannot use standard $L^2$ bounds. They develop a novel technique to estimate quartic operators by splitting them into blocks based on whether momentum legs cross the frequency cutoffs. This allows them to exploit the presence of at least one low-frequency leg to gain decay.
• Quantitative de Finetti Theorem: Used to relate the low-frequency quantum reduced density matrices to the classical moments of the finite-dimensional Gibbs measure.
• Block Decomposition: The double commutators are analyzed by classifying quartic monomials into "blocks" (e.g., $P \to P$, $P \to Q$, $Q \to Q$). Crucially, they show that the purely high-high block ($Q \to Q$) vanishes in the commutator structure, leaving only terms with at least one low-frequency leg that can be controlled.

3. Key Contributions

1. Extension to Non-Summable Potentials: This is the first rigorous derivation for the fractional Bessel interaction in the regime where $\sum \hat{v}_\beta(k) = \infty$. Previous results (e.g., [LNR15, LNR21]) were limited to summable interactions.
2. Novel Renormalization: The paper introduces a specific renormalization that centers the zero mode via a number-fluctuation term while keeping non-zero modes in the genuine quartic form, avoiding the need for a formal density-square representation that would diverge.
3. Multi-Scale High-Frequency Analysis: The authors develop a sophisticated multi-scale analysis to handle the non-positive and non-negligible "shell" terms. They separate the problem into:
   • Zero-mode fluctuations: Controlled via double commutators and kinetic energy bounds.
   • Shell terms: Controlled by polarizing fluctuation observables and using mode-by-mode estimates.
   • Tail terms: Controlled by exploiting the positivity of the full Bessel operator and bounding the mass of the far-tail number operator.
4. Convergence of Reduced Density Matrices: Beyond free energy, the paper proves the Hilbert-Schmidt convergence of the rescaled reduced density matrices to the classical correlation operators, a stronger result than scalar moment convergence.

4. Main Results

Theorem 2.1 (Main Result):
As $\lambda \downarrow 0$ with $3/2 < \beta \le 2$:

1. Free Energy Convergence: The relative free energy converges to the negative log of the classical partition function:
   $$-\log \frac{Z^{re}_\lambda}{Z_0} \longrightarrow -\log z_{cl}$$
2. Density Matrix Convergence: For every $k \ge 1$, the rescaled reduced density matrices converge in the Hilbert-Schmidt norm ($S_2$) to the classical correlation operators:
   $$k! \lambda^k \Gamma^{(k)}_\lambda \longrightarrow \int |u^{\otimes k}\rangle\langle u^{\otimes k}| \, d\nu(u)$$
   where $\nu$ is the classical Gibbs measure.

Technical Estimates:

• Theorem 5.1: Quantitative vanishing of the zero-mode high-frequency remainder.
• Theorem 6.1: Vanishing of the shell contribution.
• Theorem 7.1: Quantitative lower bound showing the tail contribution is negligible (specifically $E_{tail} \ge -o(1)$).

5. Significance

• Bridging Quantum and Classical: The paper solidifies the bridge between bosonic quantum statistical mechanics and classical random fields (specifically the Gaussian Free Field with singular interactions) in a regime previously inaccessible due to ultraviolet divergences.
• Methodological Advancement: The techniques developed for handling non-summable kernels, particularly the block decomposition of double commutators and the asymmetric Hilbert-Schmidt estimates, provide a new toolkit for analyzing singular quantum many-body systems.
• Completeness: By establishing both free energy convergence and the convergence of reduced density matrices (including the operator topology), the paper provides a complete picture of the quantum-to-classical limit for this class of interactions.
• Future Directions: The authors note that their methods can likely be extended to inhomogeneous settings (with external potentials) and potentially to the Coulomb potential ($\beta=2$), though the current paper focuses on the fractional Bessel case to isolate the core difficulties of non-summability.

In summary, this work resolves a long-standing technical hurdle in the derivation of Gibbs measures for singular interactions in 2D, providing a rigorous framework that handles the divergence of the self-energy through a careful renormalization and a multi-scale analysis of high-frequency remainders.