Ground state energy of the Bose--Hubbard model with large coordination number with a polaron-type quantum de Finetti theorem

This paper establishes the convergence of the ground state energy of the Bose-Hubbard model on graphs with large coordination number to a strong-coupling mean-field functional by developing and applying a novel "polaron-type" quantum de Finetti theorem tailored for tensor products involving bosonic Fock spaces.

The Gist

Imagine a giant, bustling city made of tiny apartments (the lattice sites). Inside each apartment, there are identical, hyper-social particles called bosons. These particles have two main behaviors:

1. Hopping: They love to run from their current apartment to a neighbor's.
2. Bumping: If too many of them try to squeeze into the same apartment, they get annoyed and push each other away (this is the "interaction").

This setup is called the Bose-Hubbard model. Physicists use it to understand how materials switch between being a solid (where particles are stuck in their spots, like a Mott Insulator) and a superfluid (where they flow freely like a liquid with zero friction, like a Superfluid).

The big question is: How do we predict what this city will do?

The Problem: Too Many Neighbors

In a real city (or a 3D crystal), every apartment has a few neighbors (maybe 6). But in this paper, the authors imagine a "super-city" where every apartment is connected to a huge number of other apartments (imagine thousands of neighbors).

When you have that many neighbors, calculating the exact behavior of every single particle is impossible. It's like trying to predict the exact path of every person in a stadium by tracking every handshake. You need a shortcut.

The Shortcut: The "Mean-Field" Guess

The standard shortcut in physics is called Mean-Field Theory. It's like saying: "Don't worry about who your specific neighbors are. Just assume you are surrounded by an 'average' crowd."

In most physics problems, this "average" works well when the interactions are weak. But in this specific model, the particles are strongly interacting (they really hate being crowded). Usually, the "average crowd" guess fails miserably here.

However, the authors prove something amazing: Even with strong interactions, if the number of neighbors is huge, the "average crowd" guess actually works perfectly!

The Magic Tool: The "Polaron" Theorem

To prove this, the authors invented a new mathematical tool called the Polaron-type Quantum de Finetti Theorem. Let's break down that scary name with an analogy:

The "Polaron" Analogy:
Imagine a famous celebrity (the "Core" particle) walking through a massive, dense crowd of fans (the "Shell" particles).

• The fans are all identical and interchangeable.
• The celebrity is unique.
• The fans are all jostling around the celebrity, but they are also jostling each other.

The de Finetti Theorem is a famous math rule that says: "If you have a huge group of identical people, they will eventually act like a perfect, random mix of independent individuals."

The New Twist:
The old rule didn't work for our celebrity because the celebrity changes the crowd's behavior. The authors created a new version of the rule specifically for this "Celebrity + Crowd" situation. They proved that even though the celebrity is special, the crowd of thousands of neighbors still behaves in a predictable, "averaged" way that can be described by a single, simple formula.

What Did They Actually Do?

1. The Setup: They took the complex math of the "super-city" with thousands of neighbors.
2. The Reduction: They showed that instead of tracking the whole city, you only need to look at one single apartment and its relationship to the "average" of all its neighbors.
3. The Proof: They used their new "Polaron Theorem" to prove that as the number of neighbors goes to infinity, the energy of the whole system converges exactly to the energy of that single, simplified apartment.

Why Does This Matter?

• It Validates a Popular Tool: Physicists have been using this "Mean-Field" shortcut for decades to design new materials and understand superconductors. This paper provides the rigorous mathematical "seal of approval" saying, "Yes, this shortcut is mathematically sound for these specific conditions."
• It Solves a Hard Problem: It proves that even when particles are strongly interacting (which usually breaks simple models), the sheer number of connections saves the day, making the system predictable.
• New Math: The "Polaron-type theorem" they invented is a new tool that other scientists can now use to solve different problems involving a "special" particle interacting with a "bath" of many others (like an electron moving through a crystal lattice).

In a Nutshell

The authors proved that in a world where every particle is connected to a massive number of neighbors, the chaos of the whole system simplifies into a single, elegant rule. They did this by inventing a new mathematical lens (the Polaron Theorem) that lets us see the forest (the whole system) clearly, even when we are standing right next to a single, complex tree.

Technical Summary

1. Problem Statement

The paper addresses the rigorous derivation of the mean-field limit for the Bose–Hubbard model on a graph with a large and homogeneous coordination number $z$ (the number of nearest neighbors per site).

• The Model: The Bose–Hubbard Hamiltonian describes bosons on a lattice with on-site interaction ($U$) and hopping between nearest neighbors ($J$).
• The Limit: The authors consider the limit $z \to \infty$. Crucially, the hopping amplitude is scaled as $J/z$ to ensure the kinetic energy remains comparable to the on-site terms in the limit.
• The Challenge: Standard mean-field theories for bosonic systems (like the Hartree limit) typically rely on the quantum de Finetti theorem, which requires the system to be symmetric under the exchange of particles. However, the Bose–Hubbard model is defined on a lattice where sites are distinct; the Hamiltonian is generally not symmetric under the exchange of lattice sites.
• Goal: To prove that the ground state energy per site of the Bose–Hubbard model converges to the minimizer of a specific mean-field energy functional as $z \to \infty$, and to establish the validity of the resulting phase diagram (Mott insulator vs. Superfluid).

2. Methodology

The authors employ a novel combination of translation invariance arguments and a new mathematical tool: the Polaron-type Quantum de Finetti Theorem.

A. Reduction to a Local "Polaron" System

Since the full lattice Hamiltonian is not symmetric under site exchange, the authors exploit the translation invariance of the graph.

1. They define a "core" site (index 0) and its $z$ neighboring sites (indices $1, \dots, z$).
2. They construct a reduced local Hamiltonian $H_{1,z}$ acting on the tensor product of the core Hilbert space and the symmetric tensor product of the $z$ neighbors.
3. This reduced Hamiltonian is symmetric with respect to the exchange of the $z$ neighbor sites.
4. The problem is thus reduced to analyzing a system consisting of one distinct particle (the core) coupled to a symmetric bath of $z$ indistinguishable particles (the neighbors). This structure is analogous to a polaron (an impurity coupled to a bosonic bath).

B. The Polaron-Type Quantum de Finetti Theorem

The core technical novelty is the development of a new version of the quantum de Finetti theorem (Theorem 7) tailored for this specific structure:

• Standard de Finetti: Applies to symmetric states of $N$ particles, stating that reduced density matrices converge to a convex combination of product states.
• Polaron-type de Finetti: Applies to a state in $\mathcal{H}_0 \otimes \mathcal{F}_+(\mathcal{H}_s)$, where $\mathcal{H}_0$ is the "core" space and $\mathcal{F}_+(\mathcal{H}_s)$ is the bosonic Fock space of the "shell" (neighbors).
• Result: As the number of shell particles $N \to \infty$, the reduced density matrix of the core and $k$ shell particles converges to an integral over a probability measure of product states:
  $$\gamma^{(1,k)} = \int_{S\mathcal{H}_s} \zeta(u) \otimes |u\rangle\langle u|^{\otimes k} \, dP(u)$$
  Here, $\zeta(u)$ is a density matrix on the core space dependent on the parameter $u$ (the state of the shell), and $|u\rangle$ represents the coherent state of the shell.
• Technical Scope: The theorem handles infinite-dimensional Hilbert spaces (unlike previous finite-dimensional versions) and explicitly constructs the de Finetti measure.

C. Energy Bounds

• Upper Bound: Proven trivially using a Gutzwiller product state (a tensor product of identical single-site states) as a trial state. This yields the mean-field energy functional directly.
• Lower Bound: This is the non-trivial part.
  1. They establish bounds on the moments of the number operators ($M_2$) to ensure the minimizing sequences are compact in the appropriate topology.
  2. They use the Polaron-type de Finetti theorem to represent the limiting reduced density matrix.
  3. By substituting this representation into the energy functional, they show that the energy is bounded below by the infimum of the mean-field functional.

3. Key Contributions

1. Rigorous Convergence of Ground State Energy: The paper provides the first rigorous proof that the ground state energy of the Bose–Hubbard model with large coordination number converges to the mean-field limit.
2. Polaron-Type de Finetti Theorem: The authors introduce a new theorem (Theorem 7) that generalizes the quantum de Finetti theorem to systems with a "core-shell" structure (one distinct subsystem coupled to a symmetric bath). This is expected to be a broadly useful tool for other many-body problems involving impurities or tracer particles.
3. Strong Coupling Regime: Unlike weak-coupling mean-field limits (where interactions are scaled down), this work treats the interaction $U$ as unscaled. This places the mean-field theory in the strong coupling regime, which is essential for capturing the Mott insulator phase transition.
4. Handling Non-Symmetric Lattices: The work successfully bridges the gap between lattice models (which lack particle exchange symmetry) and de Finetti techniques by utilizing translation invariance to create a local symmetric subsystem.

4. Main Results

• Theorem 2 (Convergence): For $U > 0$ and $J \geq 0$, the ground state energy per site of the Bose–Hubbard model converges to the minimum of the mean-field energy functional:
  $$\lim_{z \to \infty} \frac{1}{|V_z|} \inf_{\psi} \langle \psi, H_{V_z} \psi \rangle = \inf_{\phi \in \ell^2(\mathbb{C}), \|\phi\|=1} \langle \phi, h_\phi \phi \rangle$$
  where $h_\phi$ is the nonlinear mean-field operator derived from averaging the hopping term.
• Phase Diagram Validity: The results justify the use of mean-field theory (specifically the Gutzwiller approximation) to predict the phase diagram of the Bose–Hubbard model (Mott Insulator vs. Superfluid) for lattices with large coordination numbers (e.g., high-dimensional cubic lattices). The mean-field prediction is shown to be qualitatively correct even for moderate coordination numbers.

5. Significance

• Mathematical Physics: This work resolves a long-standing gap in the rigorous justification of mean-field limits for lattice bosons with strong interactions. Previous rigorous results were largely restricted to continuum systems or weak coupling.
• Connection to DMFT: The approximation derived is mathematically equivalent to the Dynamical Mean-Field Theory (DMFT) used widely in condensed matter physics. This paper provides the rigorous mathematical foundation for DMFT in the context of the Bose–Hubbard model.
• New Analytical Tool: The "Polaron-type de Finetti theorem" opens new avenues for analyzing composite quantum systems where a specific degree of freedom (like an impurity or a specific lattice site) interacts with a large, symmetric environment. This is applicable to polaron models, tracer particles in fluids, and other impurity-bath systems.
• Phase Transitions: By validating the mean-field limit in the strong coupling regime, the paper supports the existence and nature of the Mott insulator-to-superfluid transition in high-dimensional lattices, a phenomenon difficult to analyze with standard perturbation theory.

In summary, the paper successfully combines advanced functional analysis (de Finetti theorems) with lattice model physics to rigorously justify the mean-field description of strongly interacting bosons on high-coordination lattices, introducing a powerful new mathematical tool in the process.