Eigenvalue asymptotics of Müller minimizers for atoms and molecules

This paper establishes that under specific conditions on nuclear charge and electron number, the $k$-th eigenvalue of a Müller minimizer for atoms and molecules decays asymptotically as $A_* k^{-8/3}$ as $k \to \infty$, with the constant determined by the minimizer's density.

The Gist

Imagine you are trying to understand the structure of a giant, complex city (an atom or molecule) by looking at the traffic patterns of its citizens (electrons). In quantum physics, we don't just count cars; we look at the "probability map" of where every car might be at any given time. This map is called the density matrix.

For decades, scientists have used a standard map called Hartree-Fock to predict how atoms behave. It's a great map, but it has a flaw: it treats the interaction between cars (electrons) a bit too simply, like they are ghosts that pass through each other without really bumping.

Enter the Müller Functional. This is a newer, more sophisticated map that tries to fix that flaw. It's like upgrading from a 2D street map to a 3D simulation that accounts for how cars actually repel each other and change lanes. The big question for mathematicians has been: If we use this better map, what does the "traffic flow" look like at the very edges of the city?

This paper, written by a team of mathematicians, answers that question. Here is the breakdown in simple terms:

1. The Problem: The "Infinite" City

In the old Hartree-Fock model, the traffic map usually stops after a certain number of cars. It's like a city with a finite number of lanes.

But in the Müller model, the city is different. The map shows that there are infinitely many lanes, even if the number of cars is finite. The "lanes" get thinner and thinner as you go further out. The authors wanted to know: How fast do these lanes get thinner as you go to infinity?

2. The Discovery: The "8/3" Rule

The authors discovered a precise mathematical rule for this thinning. They found that if you look at the $k$-th lane (where $k$ is a huge number), its width shrinks according to a specific formula: $k^{-8/3}$.

The Analogy:
Imagine a giant staircase where the steps get smaller and smaller as you go up.

• In the old model, the steps might disappear after a while.
• In the Müller model, the steps go on forever, but they shrink at a very specific, predictable rate.
• The authors proved that this rate is exactly $k^{-8/3}$.

Why is this cool? Because this specific rate ($8/3$) is the same rate found in the most fundamental laws of quantum mechanics (the Schrödinger equation). It means the Müller model, despite being an approximation, captures the "true" quantum behavior of nature perfectly, even in the far reaches of the atom.

3. The Challenge: The "Rough Terrain"

To prove this, the authors had to deal with two major obstacles, which they describe as "rough terrain" in their mathematical landscape:

• The "Cusp" (The Sharp Corner): Electrons hate being too close to the nucleus (the center of the city) or to each other. This creates a "sharp corner" or a "cusp" in the mathematical map. It's like trying to drive a car over a jagged rock; the smoothness of the road breaks down. The authors had to prove exactly how "rough" this road is. They used a special mathematical tool called Besov spaces (think of it as a high-resolution ruler) to measure the roughness and found it was just rough enough to create the $8/3$ rule.
• The "Fog" (Decay at Infinity): As you move away from the atom, the electron density should drop to zero. But proving how fast it drops is hard because the equations get messy. The authors had to show that the "fog" clears up exponentially fast (like a light turning off in a dark room) under certain conditions. They proved that if the atom is big enough and has the right number of electrons, this fog clears up predictably.

4. The "Jastrow Factor": The Magic Filter

To handle the rough terrain (the sharp corners), the authors used a clever trick called a Jastrow factor.

• The Metaphor: Imagine you are trying to draw a perfect circle, but your hand shakes when you get close to the center. Instead of trying to draw the whole circle perfectly, you put a special "stabilizing filter" over your hand near the center. This filter absorbs the shaking, allowing you to draw the rest of the circle smoothly.
• In the paper, they mathematically "filter out" the sharp corners near the nucleus and the electrons. Once they removed the roughness, the rest of the math became smooth and solvable.

5. Why This Matters

This isn't just about abstract math.

• For Chemists: It validates the Müller functional. It tells us that this method is reliable for simulating large molecules because it respects the fundamental laws of physics, even in the "tails" of the electron clouds.
• For Mathematicians: It solves a long-standing puzzle about how these specific energy functions behave. It connects the dots between a simplified model (Müller) and the complex reality (Schrödinger).

Summary

Think of this paper as a team of cartographers who finally figured out the exact shape of the "fading edge" of a quantum city. They proved that even though the Müller map is an approximation, it fades away at the exact same speed as the real universe does. They did this by smoothing out the jagged rocks in the middle of the map and proving that the fog at the edges clears up quickly enough to make the math work.

The result? A beautiful, universal rule: The further you go, the faster the electrons fade, following the $k^{-8/3}$ law.

Technical Summary

1. Problem Statement

The paper investigates the spectral properties of the Müller functional, a model in quantum chemistry used as an alternative to Hartree–Fock theory. The functional describes $N$-electron systems with total nuclear charge $Z$ via a one-particle density matrix (1-pdm) $\gamma$.

• The Functional: The Müller energy $E_M(\gamma)$ differs from Hartree–Fock by replacing the standard exchange term $X(\gamma)$ with $X(\gamma^{1/2})$. This modification renders the energy functional convex, ensuring the uniqueness of the electron density $\rho_\gamma$ when a minimizer exists.
• The Phenomenon: Unlike Hartree–Fock minimizers (which have finite rank), Müller minimizers $\gamma_*$ possess infinitely many non-zero eigenvalues $\lambda_1 \ge \lambda_2 \ge \dots \ge 0$.
• The Goal: The authors aim to determine the asymptotic behavior of these eigenvalues $\lambda_k$ as $k \to \infty$. Specifically, they seek to prove a formula of the form $\lambda_k \sim A_* k^{-8/3}$ and to explicitly determine the constant $A_*$ in terms of the electron density $\rho_{\gamma_*}$.

2. Methodology

The proof is inspired by Sobolev's recent work on the asymptotics of the 1-pdm for Schrödinger ground states but requires significant new analytical tools due to the non-linear and non-local nature of the Müller equation. The methodology is divided into three main pillars:

A. Regularity Analysis (Singular Behavior)

The core of the eigenvalue decay rate lies in the smoothness of the integral kernel of the operator $\Phi = \gamma_*^{1/2}$.

• Euler–Lagrange Equation: The minimizer satisfies a two-body equation involving a mean-field operator $H_{\gamma_*}$. The potential in this equation is singular at the electron-nucleus coalescence points and the electron-electron coalescence diagonal ($x=y$).
• Besov Spaces: The authors establish that $\Phi$ belongs to the Sobolev space $H^{2+\alpha}$ and, crucially, the Besov space $B^{5/2}_{2,8}$.
• Jastrow Factor: To handle the singularities, they introduce a Jastrow factor $F(x,y)$ (a specific exponential weight) to "extract" the non-smooth behavior. The transformed function $\Psi = e^{-F}\Phi$ is shown to have improved regularity, belonging to $H^{3+\alpha} \cap B^{7/2}_{2,8}$.
• Novelty: These regularity results in the scale of Besov spaces are novel for Müller theory and are proven to be optimal; the eigenvalues would decay faster if the regularity were higher.

B. Decay Estimates (Exponential Decay)

To apply spectral asymptotic theorems, the density $\rho_{\gamma_*}$ must decay exponentially at infinity.

• Challenge: Standard decay arguments fail because the chemical potential $\mu$ might lie in the essential spectrum, and the Müller equation is not a standard eigenvalue problem.
• Conditional Decay: The authors prove that if the chemical potential satisfies $\mu < -1/2$, then the density decays exponentially: $\int e^{2\kappa|x|} \rho_{\gamma_*}(x) dx < \infty$.
• Spectral Gap Condition: They derive an improved upper bound for $\mu$ using a variational principle. They show that for atoms with large $Z$ and $N \le Z - C_0 Z^{1/3}$, the condition $\mu < -1/2$ holds. This relies on comparing the Müller operator to the hydrogen atom Hamiltonian.

C. Spectral Asymptotics (Schatten Classes)

With regularity and decay established, the authors apply operator theory:

• Schatten Classes: They utilize the theory of Schatten classes $S_{p, \infty}$. They prove that $\gamma_*^{1/2}$ belongs to $S_{3/4, \infty}$.
• Kernel Decomposition: The kernel $\Phi$ is decomposed into a smooth part and a singular part (related to the Jastrow factor).
• Birman–Solomyak Theorem: They apply a theorem by Birman and Solomyak regarding integral operators with kernels in Besov spaces to determine the precise asymptotic constant.

3. Key Contributions and Results

Main Theorem (Theorem 1.1)

For a Müller minimizer $\gamma_*$ with eigenvalues $\lambda_k$, under suitable conditions on $Z$ and $N$ (specifically $N \le Z - C_0 Z^{1/3}$ for large $Z$ in the atomic case), the eigenvalues satisfy:
$$\lim_{k \to \infty} k \lambda_k^{3/8} = \frac{\sqrt{2}}{3\pi^{5/4}} \int_{\mathbb{R}^3} \rho_{\gamma_*}(x)^{3/4} dx$$
This implies the asymptotic law $\lambda_k \sim A_* k^{-8/3}$.

Key Technical Results

1. Regularity of Müller Minimizers (Theorem 1.2):
   • $\Phi = \gamma_*^{1/2} \in H^{2+\alpha} \cap B^{5/2}_{2,8}$.
   • $\Psi = e^{-F}\Phi \in H^{3+\alpha} \cap B^{7/2}_{2,8}$.
   • These results are optimal; the non-smoothness at the diagonal $x=y$ dictates the $k^{-8/3}$ decay.
2. Exponential Decay (Theorem 1.3):
   • Proven conditional on $\mu < -1/2$. The proof involves a non-standard integration of the two-body equation against an exponential weight, handling the unbounded exchange operator carefully.
3. Improved Chemical Potential Bound (Theorem 1.4):
   • Establishes $\mu < -1/2$ for large atoms with $N \le Z - C_0 Z^{1/3}$. This validates the decay assumption required for the main theorem.

4. Significance and Implications

• Quantum Chemical Accuracy: The result confirms that the Müller model captures the correct "true quantum behavior" regarding eigenvalue decay. The exponent $-8/3$ matches the decay found in the full many-body Schrödinger equation (by Sobolev), distinguishing it from Hartree–Fock theory where the decay is different (typically exponential or faster algebraic depending on the gap).
• Mathematical Rigor: The paper provides the first rigorous derivation of eigenvalue asymptotics for the Müller functional. It bridges the gap between computational chemistry usage and rigorous mathematical analysis.
• New Analytical Tools: The introduction of Besov space regularity for the square root of the density matrix and the specific handling of the Jastrow factor in the context of the Müller equation offer new techniques for analyzing non-linear mean-field theories with singular potentials.
• Universality vs. Specificity: The paper highlights that while the rate of decay ($k^{-8/3}$) is universal for this class of problems (reflecting the singularity of the Coulomb interaction), the constant $A_*$ is system-dependent and determined by the electron density, unlike some universal constants in physics.

In summary, the paper successfully characterizes the spectral tail of Müller minimizers, proving that they exhibit the same slow algebraic decay as the exact Schrödinger ground states, thereby validating the Müller functional as a robust model for capturing essential quantum mechanical features of electron correlation.