Sharp upper bounds for the density of relativistic atoms: Noninteracting case

This paper establishes an optimal upper bound for the electron density of noninteracting relativistic infinite Bohr atoms described by Chandrasekhar and Dirac operators, including a separate analysis for each angular momentum channel.

The Gist

Imagine the atom not as a tiny solar system, but as a bustling, chaotic city. In this city, the nucleus is a massive, heavy skyscraper in the center, and the electrons are millions of tiny, hyper-fast commuters zooming around it.

For a long time, physicists have tried to predict exactly how crowded the streets are at any given distance from the skyscraper. This "crowd density" is what the paper calls the electron density.

Here is the problem: When these electrons get very close to the skyscraper (the nucleus), they don't just slow down; they start moving at speeds close to the speed of light. This is where the rules of "Relativity" kick in, making the math incredibly messy.

This paper, written by Rupert Frank and Konstantin Merz, is like a master architect finally drawing up the perfect, strict building code for how crowded this city can get. They prove the absolute tightest possible limit on how many electrons can squeeze into a specific spot, especially near the center.

Here is a breakdown of their discovery using simple analogies:

1. The Two Types of "Traffic Laws"

The authors looked at two different models of how these electrons behave:

• The Chandrasekhar Model: Think of this as a "simplified" version of the city where the electrons are like cars that can't spin or turn sharply. They just zoom in straight lines.
• The Dirac Model: This is the "realistic" version. Here, the electrons are like acrobats. They have a property called "spin" (like a spinning top) and they move in a more complex, twisting way.

The authors wanted to know: No matter how many electrons we have, is there a hard limit to how dense the traffic can get near the skyscraper?

2. The "Singularity" (The Black Hole of Traffic)

In the old, non-relativistic world (where electrons are slow), the traffic density near the center is manageable. But in the relativistic world, as you get closer to the nucleus, the density tries to shoot up to infinity.

Imagine a funnel. As you pour water (electrons) into the top, it gets tighter and tighter at the bottom.

• The Old Math: Previous studies said, "Okay, the funnel gets tight, but we can't be sure exactly how tight it gets right at the very tip." They had a "fuzzy" estimate.
• The New Math: Frank and Merz have calculated the exact shape of the funnel tip. They proved that the density grows at a specific rate (like $1/r^{2\eta}$), and they found the precise number $\eta$ that controls this growth.

The Analogy: It's like finally measuring the exact width of a drainpipe. Before, we knew it was narrow, but now we know it is exactly 2.3 millimeters wide, and we know that no amount of water can force it to be narrower without breaking the laws of physics.

3. The "Heat" Analogy (How they solved it)

How did they prove this? They used a technique called Heat Kernel Analysis.

Imagine you drop a single drop of hot ink into a cold pool of water.

• The ink spreads out over time.
• The speed and shape of that spreading depend on the properties of the water and the container.

The authors treated the electrons like that spreading ink. By studying how this "quantum ink" spreads out from the nucleus over a tiny fraction of a second, they could reverse-engineer the rules of the city. They showed that the "ink" cannot spread faster or slower than a specific limit, which directly translates to a limit on how crowded the electrons can be.

4. Why Does This Matter?

You might ask, "Why do we care about the exact density of electrons in a theoretical atom?"

• The "Scott Correction": In the 1950s, a physicist named Scott predicted that if you have a huge atom (with a massive nucleus), the energy of the system has a specific "correction" term. This paper provides the missing piece of the puzzle to prove that prediction is rock-solid for relativistic atoms.
• Stability of Matter: If the density gets too high, the atom might collapse. By proving these sharp upper bounds, the authors are essentially proving that relativistic atoms are stable. They have shown that even with the extreme speeds of light, the electrons won't crash into the nucleus and destroy the atom.

Summary

Think of this paper as the ultimate traffic report for the most extreme city in the universe.

• Before: We knew traffic got bad near the center, but we didn't know the exact limit.
• Now: We have a precise, mathematical "speed limit" sign that says, "No matter how fast you go, the crowd density cannot exceed this specific number."

It's a victory for precision. It takes a messy, complex problem involving the speed of light and quantum mechanics and gives us a clean, sharp answer that holds up under the most extreme conditions.

Technical Summary

1. Problem Statement

The paper addresses the mathematical problem of establishing optimal (sharp) upper bounds for the electron density of relativistic atoms in the non-interacting limit (infinite Bohr atoms). Specifically, the authors investigate the behavior of the spectral density function associated with two fundamental relativistic operators:

1. The Chandrasekhar Operator ($C_\kappa$): A pseudo-relativistic operator describing a hydrogen-like atom with effective nuclear charge $\kappa$.
2. The Dirac–Coulomb Operator ($D_\nu$): The standard relativistic operator for a single electron in a Coulomb field.

The central object of study is the density of states at negative energies (for Chandrasekhar) or energies in $[0, 1)$ (for Dirac), defined as the sum of the squares of the normalized eigenfunctions:
$$\rho(x) = \sum_{E_n < 0} |\psi_n(x)|^2$$
The authors aim to determine the precise asymptotic behavior of this density $\rho(x)$ as $x \to 0$ (near the nucleus) and $x \to \infty$ (far from the nucleus). Previous work had established bounds, but they were either suboptimal near the origin or did not cover the critical coupling constant cases ($\kappa = 2/\pi$ for Chandrasekhar and $\nu = 1$ for Dirac).

2. Methodology

The authors employ a combination of heat kernel analysis, spectral theory, and angular momentum decomposition.

• Heat Kernel Bounds: A primary technique involves bounding the spectral density using the heat kernel (semigroup) $e^{-tH}$. By utilizing the spectral theorem, the density is bounded by the diagonal of the heat kernel at a specific time $t$. The authors derive sharp pointwise bounds for the heat kernels of the homogeneous operators $L_\kappa = (-\Delta)^{\alpha/2} - \kappa|x|^{-\alpha}$ and extend these to the perturbed operators using Trotter's product formula and subordination arguments.
• Hardy–Kato–Herbst Inequalities: The proofs rely heavily on the sharp Hardy–Kato–Herbst inequality to ensure the operators are well-defined and bounded from below for specific ranges of the coupling constants.
• Angular Momentum Decomposition:
  • For the Chandrasekhar case, the problem is reduced to radial operators in each angular momentum channel $\ell$. The authors analyze the singularity strength $\eta$ as a function of $\ell$, showing that higher angular momentum weakens the singularity.
  • For the Dirac case, the authors utilize the partial wave decomposition of the Dirac–Coulomb operator. They derive explicit expressions for eigenvalues and eigenfunctions involving Laguerre polynomials and confluent hypergeometric functions. They then apply Szegő's inequality for Laguerre polynomials to bound the sum of eigenfunctions.
• Summation Techniques: To recover the full density from the angular momentum channels, the authors perform careful summations over $\ell$ (for Chandrasekhar) and the quantum number $k$ (for Dirac), splitting the sums into regimes based on the distance $|x|$ to optimize the bounds.

3. Key Contributions and Results

A. The Chandrasekhar Operator (Theorem 1.1)

The authors prove a sharp upper bound for the density $1_{(-\infty, 0)}(C_\kappa)(x, x)$ for $0 < \kappa \le 2/\pi$.

• Result: There exists a constant $A_\kappa$ such that:
  $$\rho(x) \le A_\kappa \left( |x|^{-2\eta_\kappa} \mathbb{1}_{|x| \le 1} + |x|^{-3/2} \mathbb{1}_{|x| > 1} \right)$$
  where $\eta_\kappa \in (0, 1]$ is the unique solution to $(1-\eta_\kappa)\tan(\frac{\pi \eta_\kappa}{2}) = \kappa$.
• Novelty:
  • Near the origin ($|x| \le 1$): Previous results provided bounds of $|x|^{-3/2}$ for small $\kappa$ and $|x|^{-2\eta_\kappa - \epsilon}$ for larger $\kappa$. This paper establishes the optimal exponent $-2\eta_\kappa$ without the $\epsilon$ loss, valid for all $\kappa \le 2/\pi$, including the critical case $\kappa = 2/\pi$.
  • Far field ($|x| > 1$): The bound $|x|^{-3/2}$ is confirmed for all $\kappa$, including the critical case.
• Optimality: The bound is shown to be best possible near the origin because the square of the ground state eigenfunction behaves exactly as $|x|^{-2\eta_\kappa}$.

B. The Dirac–Coulomb Operator (Theorem 1.3)

The authors prove the analogue for the Dirac operator $D_\nu$ with $0 < \nu \le 1$.

• Result: There exists a constant $A_\nu$ such that:
  $$\text{Tr}_{\mathbb{C}^4} 1_{[0, 1)}(D_\nu)(x, x) \le A_\nu \left( |x|^{-2\Sigma_\nu} \mathbb{1}_{|x| \le 1} + |x|^{-3/2} \mathbb{1}_{|x| > 1} \right)$$
  where $\Sigma_\nu = 1 - \sqrt{1-\nu^2}$.
• Novelty: Similar to the Chandrasekhar case, this provides the sharp bound $|x|^{-2\Sigma_\nu}$ near the origin for all $\nu \in (0, 1]$, resolving previous suboptimal estimates that required an $\epsilon$-loss or failed at $\nu=1$.
• Methodology: The proof involves explicit summation of Laguerre polynomial expansions of the eigenfunctions, utilizing Szegő's inequality to control the growth of the series.

C. Generalization (Theorem 2.1)

The paper establishes a general bound for a class of operators $C_\kappa = (1-\Delta)^{\alpha/2} - 1 - \kappa|x|^{-\alpha}$ in dimension $d$. This generalizes the specific physical cases ($d=3, \alpha=1$) and provides a robust framework for analyzing singular potentials in fractional Laplacian settings.

4. Significance

• Resolution of the Strong Scott Conjecture: These sharp bounds are a critical ingredient in the proof of the Strong Scott Conjecture for relativistic atoms. The conjecture describes the asymptotic behavior of the ground state energy and density of heavy atoms as the nuclear charge $Z \to \infty$. The pointwise density bounds derived here were previously the "missing link" required to rigorously prove the convergence of the rescaled electron density to the relativistic Thomas-Fermi limit.
• Critical Coupling Constants: The results are significant because they hold up to the critical coupling constants ($\kappa = 2/\pi$ and $\nu = 1$), where the operators are on the verge of instability (collapse). Previous methods often broke down or lost precision at these limits.
• Optimality: By proving that the upper bounds match the lower bounds (derived from the ground state behavior), the authors confirm that the singularities of the relativistic electron density are exactly characterized by the exponents $\eta_\kappa$ and $\Sigma_\nu$.
• Future Applications: The authors conjecture that these bounds will facilitate proving the existence of limits for the rescaled density at $x \to 0$ and $x \to \infty$, further refining the understanding of relativistic atomic structure.

In summary, this paper provides the definitive mathematical characterization of the electron density singularity in relativistic non-interacting atoms, utilizing advanced heat kernel techniques and explicit spectral analysis to close gaps left by previous literature.