Conformal Data for the $O(2)$ Wilson-Fisher CFT in $(2+1)$-Dimensional Spacetime from Exact Diagonalization and Matrix Product States on the Fuzzy Sphere

This paper utilizes exact diagonalization and matrix product state techniques on a fuzzy sphere to extract conformal data for the $(2+1)$-dimensional $O(2)$ Wilson-Fisher CFT, identifying 32 primary operators and verifying their scaling dimensions against conformal bootstrap predictions and large charge expansion results.

The Gist

Imagine you are trying to understand the rules of a massive, chaotic dance party. The dancers are atoms, and the music is the laws of physics. Usually, when you try to study this party, you have to look at it through a foggy window or from a very far distance, which makes it hard to see the individual moves.

This paper is like a team of physicists who built a perfect, crystal-clear virtual room to watch this dance party up close, specifically to understand a very special moment called a "phase transition."

Here is the breakdown of their adventure, using simple analogies:

1. The Problem: The "Foggy Window"

In the real world, studying these atomic dances is hard because our computers usually simulate them on a grid (like a chessboard). But atoms don't actually live on a grid; they live in a smooth, round world. When you force them onto a square grid, you break the natural symmetry of the universe, making it hard to see the true "dance moves" (mathematical rules) that govern them.

2. The Solution: The "Fuzzy Sphere"

To fix this, the authors used a clever trick called the Fuzzy Sphere.

• The Analogy: Imagine trying to paint a perfect sphere using only a finite number of square tiles. You can't do it perfectly; you'll have gaps or overlaps.
• The Trick: Instead of tiles, they used "fuzzy" tiles. Think of it like a low-resolution digital image of a sphere. It's not perfectly smooth, but it preserves the spirit of the sphere's roundness. This allows them to keep the most important rule of the dance: Rotational Symmetry. No matter how you spin the room, the rules stay the same.

3. The Dance: The "O(2) Wilson-Fisher" Party

The specific dance they are watching is called the O(2) Wilson-Fisher Conformal Field Theory (CFT).

• What is it? It's the universal rulebook for how certain materials change states. For example, think of liquid helium turning into a superfluid (a liquid that flows with zero friction) or a magnet losing its magnetism as it gets hot.
• The Critical Point: There is a specific temperature where the material is neither solid nor liquid, neither magnetized nor unmagnetized. It's a "tipping point." At this exact moment, the system becomes scale-invariant.
• The Analogy: Imagine a crowd of people. If they are calm, they stand still. If they are angry, they run in a chaotic mob. At the "critical point," they are doing a synchronized, rhythmic wave that looks the same whether you zoom in on a few people or zoom out to see the whole stadium. This paper maps out the exact steps of that wave.

4. The Tools: "Exact Diagonalization" and "Matrix Product States"

To solve the math of this dance, they used two powerful computer techniques:

• Exact Diagonalization (ED): This is like trying to solve a giant puzzle by looking at every single piece at once. It's incredibly accurate but very heavy. They could only do this for a "small room" (about 13 dancers).
• Matrix Product States (MPS): This is a smart shortcut. Instead of looking at every piece, it looks at how the dancers are connected to their neighbors. It's like understanding a conversation by listening to pairs of people talking rather than the whole crowd shouting. This allowed them to simulate a much larger room (up to 28 dancers).

5. The Discovery: The "State-Operator Correspondence"

This is the magic key of the paper. In this special "Fuzzy Sphere" room, there is a direct link between Energy and Mathematical Rules.

• The Analogy: Usually, energy is just how fast the dancers are moving. But in this theory, every specific energy level corresponds to a specific "move" or "operator" in the rulebook.
• By measuring the energy of the system, they could directly read off the "scaling dimensions" (how the dance moves scale up or down) and "OPE coefficients" (how likely two dancers are to bump into each other and create a new move).

6. The Results: Finding the "Hidden Moves"

The team successfully identified 32 primary dance moves (operators) and their variations.

• They found the "Order Parameter" (the move that tells you if the material is magnetized or not).
• They found the "Stress-Energy Tensor" (the move that describes how the dance floor itself is stretching and squeezing).
• They found "Currents" (the flow of the dance).

The Big Win: Their numbers matched incredibly well with predictions made by a different, very famous method called the Conformal Bootstrap (which is like trying to solve the puzzle by only looking at the edges and guessing the middle). This proves their "Fuzzy Sphere" method works perfectly.

7. The "Large Charge" Connection

Finally, they looked at dancers carrying a lot of "charge" (like a group of dancers all holding a specific color balloon).

• The Theory: There is a theory that says if you have a huge number of these balloons, the dance becomes predictable and follows a simple formula (like a sound wave or a "phonon").
• The Finding: They checked this theory and found it holds up even for small numbers of balloons! They saw that the "phonon" moves (the sound waves of the dance) appear exactly where the theory predicted, connecting the chaotic critical point to the orderly superfluid phase.

Summary

In short, this paper built a perfectly round, digital playground to watch atoms dance at the moment they change states. By using smart math tricks, they mapped out the exact rules of the dance, proving that their method is as accurate as the best theories we have. It's a major step forward in understanding the fundamental laws that govern how matter behaves at its most critical moments.

Technical Summary

1. Problem Statement

The paper addresses the challenge of extracting universal conformal data (scaling dimensions and Operator Product Expansion coefficients) for the O(2) Wilson-Fisher Conformal Field Theory (CFT) in (2+1) dimensions. This CFT describes the critical point of various physical systems, including the superfluid-Mott insulator transition in the Bose-Hubbard model and the $\lambda$-transition in liquid helium-4.

Key Challenges:

• Symmetry Breaking: Traditional numerical methods (like Exact Diagonalization on a torus) break continuous rotational symmetry ($SO(3)$), making it difficult to utilize the state-operator correspondence, which maps local CFT operators to energy eigenstates on a sphere.
• Finite-Size Effects: Standard lattice discretizations introduce significant finite-size corrections that obscure the extraction of precise CFT data.
• Need for Precision: While the Conformal Bootstrap (CB) and Monte Carlo (MC) methods have provided high-precision data, independent numerical verification using different systematic uncertainties is crucial, especially given existing tensions between experimental measurements (e.g., specific heat of helium-4) and theoretical predictions.

2. Methodology

The authors employ a fuzzy-sphere regularization combined with Exact Diagonalization (ED) and Density Matrix Renormalization Group (DMRG) techniques.

A. The Fuzzy Sphere Model

• Concept: Instead of a lattice, the system is defined on a fuzzy sphere where the algebra of functions is replaced by a finite-dimensional matrix algebra. This is realized by placing charged fermions in the Lowest Landau Level (LLL) on a sphere pierced by a magnetic monopole flux.
• Symmetry Preservation: This geometry preserves the full spatial $SO(3)$ rotational symmetry and the internal $O(2)$ symmetry, allowing for a direct application of the state-operator correspondence.
• Hamiltonian: The authors construct a microscopic $S=1$ quantum spin-1 XY model on the fuzzy sphere. The Hamiltonian includes:
  • Density-density interactions ($V_0$) to enforce single occupancy per orbital.
  • Ferromagnetic XY exchange interactions ($V_1$).
  • Single-ion anisotropy ($D$) to tune the system.
• Filling: The system is set at a filling fraction of $\nu = 1/3$ (one-third of the orbitals occupied), which maps the fermionic system to the desired spin-1 model.

B. Numerical Techniques

• Exact Diagonalization (ED): Used for smaller system sizes ($N = 6$ to $13$ electrons/orbitals).
• Matrix Product States (DMRG): Used for larger system sizes ($N = 16$ to $28$), utilizing the ITensor library with bond dimensions up to $\chi = 3072$.
• Conformal Perturbation Theory (CPT): To determine the critical point ($D_c$) and the speed of light ($c$), the authors use CPT. They assume the energy spectrum follows the scaling ansatz:
  $$\delta E_o(R, D) = \frac{c(R, D)}{R} \Delta_o(R) + 4\pi g_\epsilon(R, D) f_{o\epsilon o}(R)$$
  By tuning $D$ such that the coupling to the relevant perturbation ($g_\epsilon$) vanishes, they isolate the scaling dimensions $\Delta_o$.

C. Data Extraction

• Scaling Dimensions: Extracted from the energy gaps at the critical point $D_c$.
• OPE Coefficients: Extracted from the linear response of energy levels to the perturbation $g_\epsilon$ (specifically the derivative $\partial \delta E_o / \partial g_\epsilon$).
• Large-Charge Expansion: The authors compare their results for high-charge sectors ($S_z$) against the theoretical large-charge expansion, which predicts the scaling of operators with large global $U(1)$ charge.

3. Key Contributions

1. First Comprehensive Fuzzy Sphere Study of O(2): This work provides the first extensive numerical extraction of the operator spectrum for the O(2) Wilson-Fisher CFT using the fuzzy sphere, identifying 32 primary operators and their descendants.
2. High-Precision Scaling Dimensions: The authors report scaling dimensions for primary operators (e.g., $\sigma, \epsilon, t, \chi$) that show excellent agreement with Conformal Bootstrap and Monte Carlo results.
   • Example: The scaling dimension of the order parameter $\sigma$ is fixed to the CB value ($0.519088$) for calibration, while the leading irrelevant operator $\epsilon'$ is found to have $\Delta \approx 3.79$, consistent with CB.
3. OPE Coefficients: The paper extracts diagonal OPE coefficients ($f_{o\epsilon o}$) for various sectors, including the stress-energy tensor and conserved currents, verifying consistency with bootstrap predictions.
4. Large-Charge Verification: The study numerically verifies the large-charge expansion predictions. It confirms that the scaling dimensions of operators with large $S_z$ follow the semiclassical expansion $\Delta_Q \sim Q^{3/2} + \beta Q^{1/2} + \dots$.
5. Phonon Primaries: The authors identify "phonon primaries" at the critical point, which correspond to the Goldstone modes of the ordered superfluid phase. They observe that these phonon states exist for angular momentum $L \leq S_z$, providing a numerical check on theoretical predictions regarding the crossover from phonon to vortex-antivortex excitations.

4. Key Results

• Critical Point Determination: The critical anisotropy $D_c$ was determined for various system sizes $N$, showing convergence as $N \to \infty$.
• Operator Spectrum:
  • Identified the stress-energy tensor $T_{\mu\nu}$ with $\Delta = 3$ (exact) and the conserved current $j_\mu$ with $\Delta = 2$ (exact), confirming the preservation of symmetries.
  • Found new primary operators not previously discussed in CB/MC literature, such as the Lorentz scalar/pseudoscalar $j$ ($\Delta \approx 5.23$) and higher-rank tensors.
• Agreement with Theory:
  • The extracted scaling dimensions for low-lying operators match Conformal Bootstrap results within error margins.
  • The large-charge expansion coefficients $\alpha$ and $\beta$ were fitted to be $0.33094$ and $0.27856$, respectively, matching previous Monte Carlo estimates.
• Finite-Size Scaling: The data exhibits clear convergence patterns when plotted against $R^{-\omega}$ (where $\omega$ is the dimension of the leading irrelevant operator), validating the CPT framework.

5. Significance

• Independent Verification: This work provides a crucial, independent check on the Conformal Bootstrap results using a completely different numerical approach (Hamiltonian diagonalization vs. correlation function bounds). This is vital for resolving the long-standing tension between theoretical predictions and experimental measurements of the O(2) critical exponents (specifically the correlation length exponent $\nu$).
• Methodological Advancement: It demonstrates the power of the fuzzy sphere as a tool for studying CFTs in (2+1)D, preserving continuous symmetries that are usually broken in lattice simulations.
• Connecting Phases: The study successfully bridges the gap between the ordered superfluid phase (Goldstone modes) and the critical point (phonon primaries), offering a concrete numerical realization of how gapless excitations in a broken-symmetry phase map to the CFT spectrum.
• Future Directions: The authors suggest that combining fuzzy sphere methods with Quantum Monte Carlo (which can access larger $N$) and constructing explicit infrared conformal generators could further improve precision, potentially resolving the experimental-theoretical discrepancy in the O(2) universality class.

In summary, this paper establishes the fuzzy sphere as a robust platform for high-precision CFT data extraction, successfully mapping the spectrum of the O(2) Wilson-Fisher fixed point and validating theoretical predictions regarding large-charge operators and phonon modes.