🔢 mathematics

A Complete Proof of Mueller--Ho Conjecture

This paper provides a complete proof of the Mueller–Ho conjecture from 2002 by establishing the existence of three critical thresholds that determine the global minimizers of a combined theta function and shifted theta function, thereby characterizing the optimal vortex shapes and relative positions in two-component Bose-Einstein condensates.

Original authors: Senping Luo, Juncheng Wei

Published 2026-10-09

📖 4 min read🧠 Deep dive

CC BY 4.0

Original authors: Senping Luo, Juncheng Wei

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

In the quiet, invisible world of quantum fluids, atoms can be cooled to temperatures so low that they lose their individual identities and merge into a single, collective state known as a Bose-Einstein condensate. Within these super-cold clouds, swirling currents of fluid called vortices can form. When the condensate consists of two different types of atoms mixed together, these vortices do not just float randomly; they arrange themselves into precise, repeating geometric patterns, much like soldiers standing in formation. For decades, physicists have wondered exactly what shape this formation takes. Does it form a perfect triangle, a square, or something in between? The answer depends on how strongly the two types of atoms push or pull against each other. This question is not merely about abstract shapes; understanding these patterns helps scientists grasp the fundamental rules governing how matter organizes itself at the smallest scales, with implications for everything from superconductors to the exotic states of matter found in the universe.

For twenty years, a specific prediction about these patterns, proposed by physicists Matthew Mueller and Tin-Lun Ho in 2002, remained unproven. They suggested that as the interaction strength between the two atomic components changes, the lattice of vortices would not just shift slightly, but would undergo a series of dramatic, distinct transformations. They predicted that the structure would start as a hexagonal triangle, then morph into a diamond-like rhombus, flatten into a perfect square, and finally stretch into a rectangle. However, while computer simulations and experiments had hinted at this behavior, no one had been able to provide a rigorous mathematical proof that these transitions were inevitable and that no other hidden shapes could exist. The question was whether nature strictly followed this specific script or if there were other possibilities lurking in the mathematics.

In a new study, Senping Luo and Juncheng Wei have finally settled this debate with a complete mathematical proof. They have demonstrated that the Mueller-Ho conjecture is exactly correct. The researchers showed that there are three precise tipping points where the vortex lattice changes its shape. When the interaction between the two atomic components is very weak, the vortices settle into a triangular pattern, which is the most efficient way to pack circles on a flat surface. As the interaction strength increases and crosses the first threshold, the lattice suddenly snaps into a new configuration: a rhombus, or a slanted diamond shape, where the atoms shift their positions relative to one another. This is not a gradual slide but a sharp transition where the old triangular order is replaced by a new, centered arrangement.

As the interaction grows stronger, the lattice continues to evolve. At a second, higher threshold, the rhombus straightens out completely to form a perfect square. This square shape remains stable for a range of interaction strengths, holding its form rigidly. However, once the interaction crosses a third and final threshold, the square begins to stretch. It elongates into a rectangle, with the atoms pulling further apart in one direction while staying closer in the other. The researchers calculated the exact values for these three transition points, identifying them as specific numbers that dictate when the shape must change. They proved that between these points, the lattice cannot be anything other than the predicted shape, and that the transition happens precisely when the energy of the system is minimized.

The proof required the authors to navigate a complex landscape of possibilities, ruling out every other potential shape the lattice could theoretically take. They showed that for any given interaction strength, the system will always choose the configuration that uses the least amount of energy. By carefully analyzing the mathematical functions that describe the energy of these vortex arrangements, they demonstrated that the triangular, rhombic, square, and rectangular phases are the only stable options. They also determined exactly how the atoms shift their positions within these shapes, confirming that the relative placement of the two atomic components changes in a specific, predictable way as the lattice transforms.

This work is significant because it moves the field from observation and simulation to absolute certainty. While earlier studies could show that these shapes appeared in computer models, they could not guarantee that a more exotic, unexpected shape wouldn't appear under the right conditions. Luo and Wei's proof eliminates that doubt. They have mapped the entire journey of the vortex lattice from its weakest to its strongest interactions, confirming that nature follows a strict, elegant path. The result is a complete classification of how two-component quantum fluids organize themselves, providing a solid foundation for future experiments and a deeper understanding of the geometric principles that govern the quantum world.

Technical Summary: A Complete Proof of the Mueller–Ho Conjecture

Problem Statement
The paper addresses the minimization of the energy functional associated with two-component Bose-Einstein condensates (BECs) in the presence of a large number of vortices. The energy is modeled by the sum of a classical theta function (representing the lattice energy of a single component) and a shifted theta function (representing the interaction between two components), weighted by an interaction strength parameter $\alpha \in [0, 1]$ .

Specifically, the authors investigate the global minimization of the functional:
$E_\alpha(z; a, b) = \theta(1; z) + \alpha J(z; a, b)$
where:

$z = x + iy \in \mathbb{H}$ (the upper half-plane) parameterizes the shape of the vortex lattice.
$(a, b) \in \mathbb{R}^2$ represents the relative displacement (shift) of the two components.
$\theta(1; z)$ is the classical theta function associated with the lattice $\Lambda_z$ .
$J(z; a, b)$ is the shifted theta function (Gaussian lattice sum) defined by $\sum_{m,n \in \mathbb{Z}} e^{-\pi(m^2y + (mx-n)^2/y)} \cos(2\pi(ma+nb))$ .

The problem, originally conjectured by Mueller and Ho (2002), seeks to determine the optimal lattice shape ( $z$ ) and relative displacement ( $a, b$ ) that minimize this energy for varying interaction strengths $\alpha$ . Previous work had established the minimizers for $\alpha=0$ (triangular lattice, Montgomery's Theorem) and for the centered case ( $a=b=1/2$ ), but the full phase diagram involving simultaneous optimization of shape and displacement remained unproven.

Methodology
The proof employs a rigorous combination of analytic number theory, variational calculus, and modular form theory. The core strategy involves reducing the infinite-dimensional minimization problem to a finite set of candidate configurations through the following steps:

Modular Reduction and Compactness:
The authors utilize the modular invariance of the energy functional under the group generated by $z \mapsto -1/z$ , $z \mapsto z+1$ , and $z \mapsto -\bar{z}$ . This allows them to restrict the search for minimizers to a fundamental domain $D_G$ where $0 \le x \le 1/2$ , $x^2 + y^2 \ge 1$ , and $y < 2$ . Compactness arguments (Lemma 2.3) ensure that a minimum is attained within this reduced domain.
Displacement Localization and Fourier Estimates:
A critical component is the "Elimination of the displacement" (Lemma 2.4). By analyzing the Fourier series expansion of the energy, the authors derive lower bounds that depend on the lattice shape $(x, y)$ but are minimized at specific displacements. They prove that for certain ranges of $\alpha$ , the global minimizer must be "centered" (i.e., $a=b=1/2$ ) or located near the "deep holes" of the triangular lattice (i.e., $a=b=1/3$ ).
Hessian Analysis and Convexity:
To establish uniqueness and local minimality, the authors compute the Hessian of the energy functional with respect to the variables $(x, y, a, b)$ .
- Triangular Neighborhood: In the vicinity of the triangular lattice ( $z = e^{i\pi/3}$ ), they prove the joint Hessian is positive definite for specific $\alpha$ ranges (Lemma 4.2). This confirms that the triangular deep-hole configuration is the unique local minimizer in this region.
- Centered Neighborhood: For the centered configurations, they utilize the "centered identity" (Lemma 2.2), which expresses the energy as a linear combination of unshifted theta functions. This reduces the problem to known classifications of theta function minimizers.
Phase Transition Analysis:
The authors define three critical thresholds $\alpha_a, \alpha_b, \alpha_c$ by comparing the energy of the triangular deep-hole configuration against the minimum of the centered configurations.
- They prove the function $\alpha \mapsto \min_{\tau} E_\alpha(\tau; 1/2, 1/2) - E_\alpha(e^{i\pi/3}; 1/3, 1/3)$ is concave.
- Using certified numerical bounds (e.g., at $\alpha = 1723/10000$ and $\alpha = 7/40$ ), they locate the unique crossing point $\alpha_a$ where the global minimizer switches from the triangular deep-hole to the centered rhombic configuration.

Key Contributions and Results
The paper provides a complete classification of the minimizers for all $\alpha \in [0, 1]$ , confirming the Mueller–Ho conjecture. The results are summarized as follows:

Phase 1 ( $0 \le \alpha < \alpha_a$ ): The global minimizer is the triangular lattice ( $z = e^{i\pi/3}$ ) with relative displacement $(a, b) = (1/3, 1/3)$ . This corresponds to the "deep hole" of the triangular lattice.
Phase 2 ( $\alpha = \alpha_a$ ): A phase transition occurs. Both the triangular deep-hole configuration and a centered rhombic configuration (with $z = e^{i\phi_{\alpha_a}}$ and $a=b=1/2$ ) are global minimizers. The transition involves a discontinuous change in both lattice shape and displacement.
Phase 3 ( $\alpha_a < \alpha < \alpha_b$ ): The minimizer is a centered rhombic lattice ( $z = e^{i\phi_\alpha}$ with $\phi_\alpha \in (\phi_{\alpha_a}, \pi/2)$ ) with fixed displacement $(a, b) = (1/2, 1/2)$ . The lattice shape changes continuously with $\alpha$ .
Phase 4 ( $\alpha_b \le \alpha \le \alpha_c$ ): The minimizer is a centered square lattice ( $z = i$ ) with displacement $(a, b) = (1/2, 1/2)$ .
Phase 5 ( $\alpha_c < \alpha \le 1$ ): The minimizer is a centered rectangular lattice ( $z = iy_\alpha$ with $y_\alpha \in (1, \sqrt{3}]$ ) with displacement $(a, b) = (1/2, 1/2)$ . At $\alpha=1$ , the aspect ratio is $y_1 = \sqrt{3}$ .

The numerical values for the thresholds are provided as:
$\alpha_a \approx 0.172356$ , $\alpha_b \approx 0.373215$ , $\alpha_c \approx 0.925649$ .

Significance
The paper claims to provide the first complete mathematical proof of the Mueller–Ho conjecture regarding the ground state structure of two-component Bose-Einstein condensates.

Mathematical Rigor: It resolves a long-standing open problem in the minimization of theta functions on lattices with shifts, extending Montgomery's theorem (1988) and Theorem B (B´etermin et al., 2022) to the coupled two-component case.
Physical Relevance: The results rigorously explain the experimental observations of vortex lattice shape transitions (from hexagonal to square/rectangular) in rotating two-component BECs as a function of inter-component interaction strength.
Methodological Advance: The work introduces a robust framework combining modular reduction, sharp Fourier estimates, and Hessian analysis to handle minimization problems involving both lattice geometry and relative displacements, a class of problems previously difficult to treat analytically.

The authors emphasize that their proof relies on exact identities, convexity arguments, and rigorous error bounds, avoiding purely numerical heuristics for the final classification, though numerical verification is used to locate the specific transition points.

Technical Summary: A Complete Proof of the Mueller–Ho Conjecture

More like this