Phase structure of (3+1)-dimensional dense two-color QCD at $T=0$ in the strong coupling limit with the tensor renormalization group

Using the tensor renormalization group method, this study investigates the phase structure of (3+1)-dimensional strong coupling two-color QCD at zero temperature and finite chemical potential, revealing that the critical exponents for the diquark condensate align with mean-field theory predictions.

The Gist

The Big Picture: Solving the "Impossible" Puzzle

Imagine you are trying to understand how a crowd of people behaves when they are packed incredibly tightly together in a room. In the world of subatomic particles, this is like studying Quarks (the building blocks of protons and neutrons) when they are squeezed into a very dense state, like inside a neutron star.

For decades, scientists have been stuck on this problem. The standard way to simulate these particles on a computer is like trying to roll dice to see what happens. But when you add "density" (squeezing them together), the math turns into a nightmare of "negative probabilities" and imaginary numbers. It's like trying to roll a die that sometimes lands on "minus three." Computers can't handle this; the simulation crashes. This is known as the Sign Problem.

The Solution:
The authors of this paper decided to stop rolling dice and start building with Legos. Instead of simulating the whole room at once, they broke the problem down into tiny, manageable blocks (called Tensors) and connected them together like a giant 3D puzzle. This method is called the Tensor Renormalization Group (TRG).

The Experiment: Two-Color QCD

To test their new Lego-building method, they didn't try to solve the full, messy real-world problem immediately. Instead, they built a "training version" of the universe called Two-Color QCD.

• Real World (3-Color QCD): Quarks come in three "colors" (Red, Green, Blue). It's like a complex traffic jam with three lanes of cars.
• The Training Version (2-Color QCD): They removed one color. Now there are only two lanes. This simplifies the math enough to avoid the "negative probability" crash, but it still keeps the essential physics of dense matter.

They also turned the temperature down to absolute zero ($T=0$). Think of this as freezing the crowd in place so they aren't jiggling around, making it easier to see how they arrange themselves when just the pressure (density) is applied.

The Lego Method: How TRG Works

Imagine you have a massive, 4-dimensional grid of Legos representing space and time.

1. The Grid: Each tiny block in the grid represents a tiny piece of the universe.
2. The Problem: If you try to look at the whole grid at once, it's too big for any computer.
3. The Trick (Renormalization): The authors use a "zoom-out" technique. They take a small cluster of Legos, figure out how they interact, and then replace that whole cluster with a single, slightly bigger "super-Lego" that acts the same way.
4. Repeating: They do this over and over, zooming out until the entire 4D universe is represented by a single, manageable number. This allows them to calculate the state of the system without the computer crashing.

What They Found: The Phase Structure

By zooming out on their giant Lego grid, they watched how the particles behaved as they increased the pressure (chemical potential). They discovered three distinct "states of matter," similar to how water changes from ice to liquid to steam:

1. The Empty State (Low Pressure): When the pressure is low, the particles are lazy. They don't pair up, and the "chiral condensate" (a measure of how much the particles are interacting with the vacuum) stays high. It's like a quiet library where everyone is sitting still.
2. The Super-Condensate State (Medium Pressure): As they squeeze harder, something magical happens. The particles suddenly decide to pair up (forming "diquarks"). This is like the library patrons suddenly grabbing hands and dancing in pairs. This is a superfluid state. The authors found the exact pressure point where this happens.
3. The Saturated State (High Pressure): If they squeeze even harder, the particles get so packed that they can't move anymore. The density hits a maximum limit (like a crowded elevator where no one can fit another person).

The Results: Did the Math Work?

The authors compared their Lego results with old, theoretical predictions (called Mean Field Theory).

• The Good News: Their results matched the old theories very well. The "critical exponents" (mathematical numbers that describe how the system changes at the tipping point) were almost identical to the predictions.
• The Significance: This proves that their "Lego method" (TRG) works perfectly for dense matter. It's a dry run.

Why Does This Matter?

This paper is a dress rehearsal.

The authors successfully used this method on the "training version" (2 colors). Now that they know the method works and can handle the massive computer power required (they used a grid size of $1024^4$, which is astronomically large), they are ready to tackle the real thing: Three-Color QCD (our actual universe).

If they can apply this same Lego method to the real, 3-color universe, they might finally be able to simulate what happens inside neutron stars or the early moments of the Big Bang, solving a mystery that has stumped physicists for 40 years.

Summary in One Sentence

The authors built a giant, 4D digital Lego puzzle to simulate dense matter, proving that this new method can accurately predict how particles behave under extreme pressure, paving the way to finally solve the mysteries of neutron stars and the early universe.

Technical Summary

1. Problem Statement

The primary challenge in studying Quantum Chromodynamics (QCD) at finite baryon density and zero temperature ($T=0$) is the complex action problem (or sign problem). In standard Monte Carlo simulations, the presence of a finite chemical potential ($\mu$) makes the action complex, rendering the probability measure non-positive and preventing importance sampling.

While Two-Color QCD (QC2D) is free from the sign problem (allowing standard Monte Carlo simulations), investigating the $T=0$ limit requires extremely large lattice sizes to approach the thermodynamic limit, which is computationally prohibitive for standard methods. Previous analytical studies in the strong coupling limit (using Mean Field theory and $1/d_s$ expansions) provided predictions, but a rigorous non-perturbative numerical verification in $(3+1)$ dimensions at $T=0$ was lacking.

2. Methodology

The authors employ the Tensor Renormalization Group (TRG) method, specifically the Grassmann-anisotropic TRG (ATRG) algorithm, to overcome the sign problem and handle large lattice volumes.

• Model Formulation:
  • Theory: $(3+1)$-dimensional QC2D with $N_c=2$ colors.
  • Action: Strong coupling limit ($g \to \infty$), where the gluonic action vanishes, leaving only the fermionic part.
  • Fermions: Kogut-Susskind (staggered) quarks with finite chemical potential $\mu$ and mass $m$.
  • Symmetry Breaking: An explicit $U(1)_V$ symmetry breaking term ($\lambda$) is added to the action to detect diquark condensation.
• Tensor Network Representation:
  • The partition function is rewritten as a product of local Grassmann tensors.
  • Auxiliary Grassmann variables ($\zeta, \xi$) are introduced to handle the hopping terms.
  • Link variables ($U_\nu$) are integrated out analytically using Weingarten calculus, resulting in a local tensor $T$ that depends on fermion occupation numbers and color indices.
  • The initial tensor size is large ($16^8$ components) due to the color degree of freedom, necessitating efficient handling.
• Numerical Algorithm:
  • Anisotropic TRG (ATRG): Used to coarse-grain the tensor network. Unlike the Bond-Weighted TRG (BTRG) used for $(1+1)$D, ATRG is suitable for $(3+1)$D.
  • Optimization: The initial tensor is constructed as a sparse matrix by exploiting Kronecker delta constraints from the integration, avoiding naive loops over all indices.
  • Truncation: Singular Value Decomposition (SVD) is used to truncate the bond dimension to $D=55$.
  • Hardware: Multi-GPU parallelization is utilized to handle the computational load on a lattice volume of $1024^4$.

3. Key Contributions

1. First TRG Application to $(3+1)$D Strong Coupling QC2D: The paper successfully extends the TRG method from lower dimensions and effective models to full $(3+1)$D QC2D, overcoming the computational difficulty of the enlarged initial tensor size (increased by a factor of $2^{16}$ compared to NJL models).
2. Thermodynamic Limit at $T=0$: By utilizing the TRG method, the authors achieved simulations on a $1024^4$ lattice, effectively reaching the thermodynamic limit at zero temperature, a regime inaccessible to standard Monte Carlo methods.
3. Critical Exponent Determination: The study provides a non-perturbative calculation of critical exponents associated with the diquark condensate, testing the validity of Mean Field (MF) theory predictions in the strong coupling limit.

4. Key Results

The study evaluated the chiral condensate ($\langle \bar{\chi}\chi \rangle$), diquark condensate ($\langle \chi\chi \rangle$), and quark number density ($\langle n \rangle$) as functions of the chemical potential $\mu$ at fixed mass $m=1.0$.

• Phase Structure:

  • Silver Blaze Region: For $\mu < \mu_{c}^{low}$, the system remains in the vacuum state (no baryon density).
  • Diquark Condensation Phase: For $\mu_{c}^{low} \leq \mu \leq \mu_{c}^{up}$, a diquark condensate forms ($\langle \chi\chi \rangle \neq 0$). The chiral condensate decreases monotonically in this region.
  • Saturated Phase: For $\mu \geq \mu_{c}^{up}$, the diquark condensate vanishes, and the quark number density saturates at $\langle n \rangle/2 = 1$ (consistent with the Pauli exclusion principle for two-color QCD).
  • Agreement with MF: The qualitative behavior and the boundaries of the phases ($\mu_{c}^{low} \approx 1.095$, $\mu_{c}^{up} \approx 1.19$) show strong consistency with previous Mean Field analytical predictions.

• Critical Exponents:

  • $\beta_m$ (Order Parameter Exponent): Fitting $\langle \chi\chi \rangle \propto (\mu - \mu_c)^{\beta_m}$ yields $\beta_m = 0.514(27)$. This is consistent with the Mean Field prediction of $\beta_m = 0.5$.
  • $\delta$ (Critical Isotherm Exponent): Fitting the free energy dependence on $\lambda$ yields $\delta = 2.44(6)$ (with systematic errors). This is slightly lower than the MF prediction of $\delta = 3$. The authors attribute the discrepancy to the difficulty in precisely determining the critical chemical potential $\mu_c$. However, measurements slightly away from the critical point yield $\delta \approx 2.91$, closer to the MF value.

5. Significance and Outlook

• Validation of TRG: The results demonstrate that the Tensor Renormalization Group is a robust tool for investigating finite-density QCD-like theories, successfully bypassing the sign problem and accessing the $T=0$ thermodynamic limit.
• Benchmark for Analytical Methods: The numerical results validate the Mean Field theory predictions for the phase structure and critical exponents in the strong coupling limit of QC2D.
• Future Directions: This work serves as a preparatory step for tackling full $(3+1)$D QCD (with $N_c=3$), which suffers from the sign problem. The authors suggest future work will incorporate finite gauge coupling and extend the method to Wilson quarks.

In summary, this paper establishes a new benchmark for non-perturbative lattice QCD studies at finite density, proving that TRG can handle the complexities of $(3+1)$D gauge theories with fermions and providing precise numerical data on the phase diagram of dense two-color matter.