⚛️ lattice

Tensor renormalization group study of cold and dense QCD in the strong coupling limit

Using the tensor renormalization group method, this study investigates the phase structure of (3+1)-dimensional cold and dense QCD in the strong coupling limit, determining consistent critical quark masses where first-order chiral and nuclear transitions terminate and confirming these findings against dual formulation Monte Carlo simulations and mean-field analysis.

Original authors: Yuto Sugimoto, Shinichiro Akiyama, Yoshinobu Kuramashi

Published 2026-06-19

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Yuto Sugimoto, Shinichiro Akiyama, Yoshinobu Kuramashi

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

Imagine the universe is made of tiny, invisible Lego bricks called quarks. Usually, these bricks are glued together so tightly by a force called the "strong force" that they can never be separated; they are stuck inside larger structures called protons and neutrons. This is how matter behaves in our everyday world.

However, physicists want to know what happens if you squeeze these bricks together with infinite pressure (like in the core of a neutron star) or heat them up to infinite temperatures (like right after the Big Bang). This is the realm of "cold and dense" matter. The problem is, we can't build a machine to squeeze matter that hard, and our best computer simulations get stuck because the math becomes too messy to solve (a problem known as the "complex action problem").

This paper is like a team of master architects who found a new, clever way to build a model of this extreme universe without getting stuck in the math mess.

The New Tool: The "Tensor Renormalization Group" (TRG)

Think of the standard way of simulating this universe as trying to guess the weather by flipping a coin millions of times. Sometimes the coin flips give you a confusing, contradictory result (the "complex action problem").

The authors used a different tool called the Tensor Renormalization Group (TRG).

The Analogy: Imagine you have a giant, intricate tapestry representing the universe. Instead of trying to look at every single thread at once (which is impossible), the TRG method is like a smart camera that takes a photo, zooms out, and blends the threads together into a simpler, smaller picture that still keeps the essential pattern. It repeats this process, zooming out further and further, until the whole complex tapestry is reduced to a manageable size that a computer can easily understand.
The Breakthrough: This method works perfectly even when the math gets "weird" or "imaginary," which is exactly what happens in dense matter.

What They Discovered: The "Phase Transitions"

The team wanted to see how the Lego bricks change their behavior as you add more and more of them (increasing density). They were looking for two specific "switches" or phase transitions:

The Nuclear Switch: Imagine the bricks are packed so tightly that they stop behaving like individual marbles and start flowing like a liquid. This is the "nuclear transition."
The Chiral Switch: Imagine the bricks suddenly lose their "handedness" (a property called chirality) and become symmetric. This is the "chiral transition."

In the real world, these switches might happen at different pressures. The big question was: Do they happen at the exact same time, or do they happen separately?

The Results: A Perfect Match

The authors ran their "smart camera" simulation on a massive grid (a digital lattice) to see what happens.

Finding the "Critical Point": They looked for a specific weight of the bricks (quark mass) where the "jump" between phases disappears. Below this weight, the transition is a sudden, violent snap (like ice cracking). Above this weight, the transition is a smooth slide (like melting butter).
The Discovery: They found that for the "Nuclear Switch" and the "Chiral Switch," this critical weight is almost exactly the same.
- The Metaphor: It's as if you have two different light switches in a room. Usually, you might flip one, then the other. But in this extreme universe, they found that both switches are wired to the exact same button. When you press it, both lights turn on at the exact same moment.

They also checked this on a giant grid (simulating an infinite universe) and confirmed that the "sudden snap" (first-order phase transition) really does happen, just as their simpler theories predicted.

Why This Matters (According to the Paper)

Solving the Puzzle: They proved that their new "smart camera" (TRG) method works for the 3D universe with 3 types of quarks (the real world), not just simplified 2D versions.
Consistency: Their results sit right in the middle of two other theories: one that uses a different math trick (Dual Formulation) and one that uses a rough approximation (Mean-Field). Their precise calculation suggests the truth is somewhere between those two guesses.
No Future Applications (Yet): The paper is purely about understanding the fundamental rules of the universe. They do not claim this will lead to new engines, medical treatments, or energy sources. It is a map of a territory we cannot physically visit, drawn using a new kind of compass.

In short: The authors built a new mathematical microscope to look at the densest matter in the universe. They found that two major changes in how matter behaves happen simultaneously at a specific density, confirming that our theoretical tools can finally see clearly into this hidden corner of physics.

Technical Summary: Tensor Renormalization Group Study of Cold and Dense QCD in the Strong Coupling Limit

Problem Statement
The phase structure of cold and dense Quantum Chromodynamics (QCD) is theoretically expected to be rich, featuring nuclear and chiral phase transitions, as well as potential color-superconducting phases at high densities. However, this regime is inaccessible experimentally, and standard lattice QCD simulations based on the Monte Carlo approach are obstructed by the complex action problem at finite baryon chemical potential. While the strong coupling limit offers a tractable starting point where the theory can be reformulated using dual variables to bypass the sign problem, previous studies have been limited to relatively small lattice sizes ( $O(10^4)$ ) or specific formulations (e.g., two-color QCD). There remains a need to investigate the (3+1)-dimensional QCD phase structure with dynamical fermions at larger scales to approach the thermodynamic limit and resolve the nature of phase transitions (first-order vs. crossover) and critical endpoints.

Methodology
This study employs the Tensor Renormalization Group (TRG) method to investigate (3+1)-dimensional cold and dense QCD with Kogut–Susskind quarks in the strong coupling limit. The approach involves the following technical steps:

Formulation: The partition function is derived by integrating out the $SU(N_c)$ link variables and the quark fields, resulting in a Grassmann tensor network representation. The authors utilize an anisotropic TRG (ATRG) algorithm, accelerated by bond-swapping techniques and multi-GPU parallelization.
Network Contraction: For finite-temperature calculations ( $N_\tau = 8$ ), the algorithm performs four-dimensional coarse-graining along the temporal direction (imaginary-time evolution) followed by the contraction of the resulting three-dimensional cubic tensor network. For zero-temperature calculations, the spatial directions are coarse-grained to reduce the staggered unit cell to a single tensor, followed by iterative contraction.
Observables: The study calculates the thermodynamic potential ( $f$ ), the chiral condensate ( $\langle \bar{\chi}\chi \rangle$ ), and the quark number density ( $\langle n \rangle$ ) as functions of the chemical potential ( $\mu$ ) and quark mass ( $m$ ).
Critical Endpoint Determination: The critical quark masses for the chiral ( $m_c^\chi$ ) and nuclear ( $m_c^n$ ) transitions are identified as the points where the discontinuities in the respective order parameters vanish, signaling the termination of the first-order phase transition.

Key Results

Critical Endpoints at $N_\tau = 8$ : The authors determined the critical quark masses where the first-order transitions terminate.
- For the nuclear transition (quark number density), the critical mass is found to be $m_c^n(N_\tau=8) = 2.075(23)$ .
- For the chiral transition (chiral condensate), the critical mass is $m_c^\chi(N_\tau=8) = 2.0545(34)$ .
- The results indicate that both transitions exhibit discontinuities at the same chemical potential for a given quark mass, and the derived critical masses are consistent with each other within fitting errors.
Comparison with Other Methods: The obtained critical mass ( $m_c \approx 2.06$ ) lies between the value estimated from Monte Carlo simulations in the dual formulation ( $m_c \approx 1.7$ ) and the mean-field (MF) prediction ( $m_c \approx 2.4$ ).
Zero-Temperature Limit: To confirm the nature of the phase transition in the thermodynamic limit, simulations were performed on a $1024^4$ lattice at zero temperature with a quark mass of $m=2.07$ . The results show a clear, sharp first-order phase transition in both the chiral condensate and quark number density at $\mu \approx 1.591$ , consistent with mean-field predictions. The transition point determined from the chiral condensate ( $1.592(1)$ ) agrees with that from the quark number density ( $1.5915(5)$ ).

Significance and Claims
The paper claims to successfully extend the TRG method from two-color QCD to (3+1)-dimensional QCD with three colors, effectively handling the increased degrees of freedom without encountering the sign or complex action problems that plague Monte Carlo methods. The study demonstrates the capability of the TRG approach to explore the cold and dense region of QCD by:

Providing precise determinations of critical endpoints for chiral and nuclear transitions in the strong coupling limit.
Confirming the existence of a first-order phase transition at finite quark mass in the thermodynamic limit, a regime previously difficult to access with high precision.
Establishing a baseline for future investigations, noting that the method is theoretically capable of incorporating dynamical non-Abelian gauge fields at finite couplings and could be extended to Wilson quarks to search for color-flavor locked phases.

The authors emphasize that while the current work is restricted to the strong coupling limit, it serves as a preparatory and validation step for applying tensor network methods to more realistic QCD scenarios where the complex action problem is a barrier.

The New Tool: The "Tensor Renormalization Group" (TRG)

What They Discovered: The "Phase Transitions"

The Results: A Perfect Match

Why This Matters (According to the Paper)

More like this