An enhanced term in the Szegő-type asymptotics for the free massless Dirac operator

This paper establishes enhanced Szegő-type asymptotic expansions for the regularized Fermi projection of the massless Dirac operator on $d$-dimensional cubes, demonstrating that while analytic test functions yield a $d$-term expansion with a logarithmic error bound, polynomial test functions of degree three or less admit a $(d+1)$-term expansion with a constant-order error, featuring an additional logarithmic term whose coefficient is independent of regularization.

The Gist

Imagine you are a physicist trying to count the number of electrons in a giant, invisible box. But this isn't a normal box; it's a mathematical "universe" where the rules of quantum mechanics apply. Specifically, we are looking at massless Dirac particles—think of them as ultra-fast, weightless messengers (like photons, but with a bit more spin) that zip around at the speed of light.

The paper by Leon Bollmann is about figuring out exactly how many of these particles are "packed" into a specific region of space as that region gets bigger and bigger.

Here is the story of the paper, broken down into simple concepts:

1. The Setup: The "Fermi Sea" and the "Box"

In quantum physics, particles fill up energy levels like water filling a glass. The "Fermi energy" is the waterline. Below this line, the particles are packed tight; above it, it's empty.

• The Problem: The author is looking at a special case where the waterline is exactly at zero energy. For massless particles, this is a tricky spot because their energy graph looks like two cones touching at a single point (the origin).
• The Box: To count the particles, the author puts them inside a giant $d$-dimensional cube (a square in 2D, a cube in 3D, or a hyper-cube in higher dimensions) and scales this box up to infinity.

2. The Smooth vs. The Rough (The "Discontinuity")

Usually, when you count things in a big box, the answer is simple:

• Volume Term: The total number is mostly just the size of the box (Volume).
• Surface Term: There is a smaller correction based on the surface area of the box.

However, in this specific case, the "symbol" (the mathematical rulebook describing the particles) has a glitch or a discontinuity right at the center (the origin).

• Analogy: Imagine trying to paint a wall. If the wall is smooth, you need a predictable amount of paint. But if there is a sudden, jagged crack right in the middle of the wall, the amount of paint you need changes in a weird way.
• In math terms, this "crack" usually creates a logarithmic correction (a term involving $\log L$, where $L$ is the size of the box). It's a "bonus" term that grows very slowly, but it's there.

3. The Main Discovery: Peeling the Onion

The author wanted to know: "If we have this glitch, can we predict the number of particles even more precisely? Can we see the next layer of the onion?"

Most previous studies could only see the first two layers (Volume and Surface). This paper goes deeper.

• The Result: The author proves that for a box shaped like a cube, you can actually predict $d+1$ terms in the expansion.
  • Terms 1 to $d$: These are the standard volume and surface corrections (like the smooth wall).
  • The Special Term: The $(d+1)$-th term is the logarithmic enhancement. It's the "bonus" caused by the glitch at the center.

4. The "Corner" Connection

Why does the shape of the box matter?

• The Metaphor: Imagine the glitch at the center of the universe is a tiny, sharp point. The box has corners. The author discovered that the "bonus" term depends on how the corners of the box interact with that central glitch.
• It's like shouting in a canyon. If the canyon walls are smooth, the echo is one thing. If the canyon has sharp, jagged corners, the echo changes. The author calculated exactly how the "echo" of the central glitch bounces off the corners of the cube to create that specific logarithmic term.

5. The "Magic" of Polynomials

The author had a limitation: they could only prove this detailed formula if the "test function" (the way they weigh the particles) was a simple polynomial (like $x$, $x^2$, or $x^3$).

• Why? The math gets incredibly messy (like trying to untangle a knot of headphones) when the function gets too complex.
• The Good News: Even with this restriction, they proved that the "bonus" logarithmic term is universal. It doesn't depend on how we smoothed out the glitch (the regularization parameter). It's a fundamental property of the system, just like the speed of light.

6. The Big Picture: Why Should We Care?

This isn't just abstract math; it connects to Entanglement Entropy.

• The Real World: In quantum computing and condensed matter physics, scientists want to know how "entangled" two parts of a system are.
• The Law: Usually, entanglement scales with the surface area of the region (the "Area Law").
• The Twist: This paper shows that for massless particles, there is a logarithmic enhancement to this area law. It's a subtle but crucial correction that tells us more about the quantum nature of the vacuum.

Summary in One Sentence

Leon Bollmann figured out that if you count massless quantum particles in a giant cube, the "glitch" at the center of their energy spectrum creates a tiny, special "echo" (a logarithmic term) that depends on the cube's corners, and he proved exactly how to calculate this echo for simple cases.

The "Takeaway" Analogy:
If counting particles in a box is like counting the number of tiles on a floor, most people just count the area. This paper says, "Wait, there's a weird stain in the middle of the floor that makes the tiling pattern slightly different near the corners. I've figured out the exact formula for how that stain changes the count, and it turns out the stain's effect is independent of how we cleaned it up!"

Technical Summary

1. Problem Statement and Context

The paper investigates the asymptotic behavior of the trace of functions of Wiener–Hopf operators associated with the free massless Dirac operator in $d$-dimensional space ($d \geq 2$). Specifically, it focuses on the regularized Fermi projection at Fermi energy zero ($E_F = 0$).

• The Operator: The Hamiltonian is the massless Dirac operator $\mathfrak{D} = -i \sum \alpha_k \partial_k$. Its symbol is $D(\xi) = \sum \alpha_k \xi_k$.
• The Singularity: Unlike the massive case or the Schrödinger operator, the massless Dirac operator has a dispersion relation consisting of two cones meeting at the origin. At Fermi energy zero, the symbol of the corresponding Wiener–Hopf operator features a zero-dimensional discontinuity (a point singularity) at the origin in momentum space.
• The Challenge: Standard Szegő-type asymptotics for smooth symbols yield an expansion in powers of the scaling parameter $L$ (volume term $L^d$, surface term $L^{d-1}$, etc.). For symbols with $(d-1)$-dimensional discontinuities (e.g., Fermi surfaces), the Widom–Sobolev formula predicts a logarithmically enhanced area term ($L^{d-1} \log L$).
• The Specific Question: How does a zero-dimensional discontinuity affect the asymptotic expansion? Previous work established that in 1D, this leads to an $O(\log L)$ term. In higher dimensions, it was unclear if such a term exists and how it interacts with the geometry of the spatial domain (specifically, the corners of a cube).

The goal is to determine the asymptotic expansion of:
$$\text{tr} \left[ g \left( 1_{\Lambda_L} \text{Op}(\psi^{(b)} D) 1_{\Lambda_L} \right) \right]$$
as $L \to \infty$, where $\Lambda_L = [0, 2L]^d$ is a scaled cube, $g$ is an analytic test function with $g(0)=0$, and $\psi^{(b)}$ is a smooth ultraviolet cutoff.

2. Methodology

The author employs a combination of techniques from spectral theory, pseudo-differential operator theory, and harmonic analysis, adapting methods from both continuous and discontinuous symbol cases.

A. Kernel Decay Estimates

The core difficulty is that the zero-dimensional discontinuity causes the integral kernel of the operator to decay only as $|x-y|^{-d}$, which is the borderline case for trace-class properties in $d$ dimensions.

• Hilbert-Schmidt Estimates: The author derives precise Hilbert-Schmidt norm estimates for operators involving products of projections and the Dirac operator.
• Geometric Decomposition: The cube $\Lambda_L$ is decomposed into regions associated with its faces and vertices. The author constructs specific "thickened" neighborhoods of the boundaries of these regions using power functions $f(x) = x^q$ ($0 < q < 1$) to control the decay of the kernel across different scales. This allows for the separation of contributions from the bulk, faces, and corners.

B. Commutation Strategy (Momentum Space)

To isolate the singular contribution, the author performs a commutation procedure to separate the smooth cutoff $\psi$ from the singular Dirac symbol $D$.

• Goal: Transform the operator $\text{Op}(\psi D)$ into a form where the smooth function $\psi$ is commuted to the right, leaving an operator involving only the singular symbol $D$ and the smooth function $g(\psi)$.
• Technique: This involves iteratively commuting $\text{Op}(\psi)$ past projections $1_{H_k}$ (half-spaces). The error terms generated by these commutators are shown to be of constant order $O(1)$ using the specific decay properties of the kernel and the structure of the Dirac matrices (vanishing trace, anti-commutation relations).

C. Local Asymptotic Formula

For the remaining singular term, the author derives a local asymptotic formula.

• Scaling: The problem is reduced to a local analysis near the vertices of the cube using unitary dilation.
• Kernel Properties: For polynomial test functions (specifically degree $\leq 3$), the author explicitly analyzes the integral kernel of the operator $\text{Op}(D)^k$. They prove that for $g(z)=z^2$ and $z^3$, the kernel is continuous on the positive quadrant, homogeneous of degree $-d$, and has a non-vanishing trace.
• Generalization of Widom's Proof: The author adapts the one-dimensional proof by Widom (which uses the Hilbert transform and Mellin transform) to higher dimensions. Instead of the Hilbert transform, the analysis relies on Riesz transforms. Due to the complexity of the higher-dimensional structure, a full generalization for arbitrary analytic functions is not achieved; instead, a logarithmic upper bound is proven for the remainder.

3. Key Contributions and Results

The paper establishes two main theorems regarding the asymptotic expansion as $L \to \infty$:

Theorem 2.1: General Analytic Test Functions

For any entire function $g$ with $g(0)=0$, the trace admits a $d$-term expansion plus a logarithmic error bound:
$$\text{tr}[\dots] = \sum_{m=0}^{d-1} (2L)^{d-m} A_{m,g,b} + O(\log L)$$

• Coefficients $A_{m,g,b}$: These depend on the regularization parameter $b$ and correspond to the geometric structure of the cube (faces of dimension $d-m$).
• Error Term: The remainder is bounded by $O(\log L)$. Crucially, the constant in this bound is independent of the regularization parameter $b$.

Theorem 2.3: Polynomial Test Functions (Degree $\leq 3$)

If $g$ is a polynomial of degree $\leq 3$, the expansion is sharpened to a $(d+1)$-term expansion with a constant error term:
$$\text{tr}[\dots] = \sum_{m=0}^{d-1} (2L)^{d-m} A_{m,g,b} + 2^d \log L \int_{S^{d-1}_+} \text{tr}_{\mathbb{C}^n}[K_g(y,y)] \, dy + O(1)$$

• The Enhanced Term: The coefficient of the $\log L$ term is given by an integral over the positive part of the unit sphere $S^{d-1}_+$.
• Independence: This logarithmic coefficient is independent of the regularization parameter $b$.
• Specific Ratios: For monomials $g_2(z)=z^2$ and $g_3(z)=z^3$, the coefficients satisfy a specific ratio ($3/2$), a consequence of the Dirac matrix algebra (vanishing trace and anti-commutation).

4. Significance and Implications

1. Resolution of the Zero-Dimensional Discontinuity: The paper confirms that a zero-dimensional discontinuity in the symbol (characteristic of the massless Dirac operator at $E_F=0$) generates a logarithmically enhanced term in the asymptotic expansion, analogous to the $(d-1)$-dimensional discontinuities in the Widom–Sobolev formula, but interacting with the vertices (0-dimensional corners) of the spatial domain.
2. Entanglement Entropy: The results are directly relevant to the study of bipartite entanglement entropy in relativistic fermionic systems. The logarithmic enhancement suggests a specific scaling law for the entanglement entropy of massless Dirac fermions in higher dimensions, distinct from the area law of gapped systems or the standard area law of free Fermi gases.
3. Regularization Independence: A major physical insight is that the coefficient of the logarithmic term is universal; it does not depend on the specific ultraviolet regularization ($b$) used to define the Fermi projection. This suggests the term is a robust physical observable.
4. Technical Advancement: The paper develops new technical tools for handling integral kernels with critical decay ($|x|^{-d}$) in multi-dimensional settings, specifically the construction of power-function neighborhoods to manage Hilbert-Schmidt estimates. It also extends the scope of Szegő-type asymptotics to operators with matrix-valued symbols and zero-dimensional singularities.

In summary, Bollmann provides a rigorous mathematical foundation for the "enhanced area law" in relativistic quantum systems, demonstrating that the geometry of the spatial domain (specifically its corners) and the spectral properties of the massless Dirac operator conspire to produce a universal logarithmic correction to the standard asymptotic expansion.