Chiral first order phase transition at finite baryon density and zero temperature from self-consistent pole masses in the linear sigma model with quarks

Using the two-flavor Linear Sigma Model with quarks at one-loop order, this study demonstrates that the chiral phase transition at zero temperature and finite baryon density is of first order, occurring when the chemical potential equals the vacuum quark mass, as evidenced by discontinuities in the chiral condensate, particle masses, couplings, and the speed of sound.

The Gist

Imagine the universe as a giant, cosmic kitchen. Inside this kitchen, there are tiny, fundamental ingredients called quarks and gluons. Under normal conditions (like inside a proton or neutron), these ingredients are locked together in a tight, invisible cage. They can't move freely; they are "confined." This is the state of matter we see every day.

However, physicists want to know what happens if you squeeze these ingredients incredibly hard or heat them up to extreme temperatures. Do they break free? Do they change their personality? This is the study of the QCD Phase Diagram.

This paper is like a sophisticated recipe book that tries to predict exactly what happens when you squeeze these ingredients at zero temperature (so cold it's absolute stillness) but with infinite pressure (high density).

Here is the story of their findings, broken down with simple analogies:

1. The Problem: The "Frozen" Kitchen

Usually, when scientists try to simulate this high-pressure environment using supercomputers (a method called Lattice QCD), they hit a wall. It's like trying to solve a puzzle where half the pieces are invisible or change color when you look at them. This is called the "sign problem," and it makes it impossible to see what happens at high density.

To get around this, the authors used a Linear Sigma Model. Think of this model as a simplified, cartoon version of the real kitchen. Instead of tracking every single sub-atomic particle, they group them into three main characters:

• The Quarks: The heavy workers.
• The Pions: Light, ghostly messengers.
• The Sigma: A heavy, invisible "glue" field that holds everything together.

2. The Method: A Self-Consistent Dance

In previous attempts, scientists used a "lazy" approximation. They assumed the particles were static statues. But in reality, these particles are dancing. If the Quark moves, it changes the Pion's dance steps, which changes the Sigma's rhythm, which changes the Quark's steps again.

The authors did something new: they set up a self-consistent loop.

• The Analogy: Imagine a group of friends trying to decide on a movie. Friend A says, "I'll watch a comedy if Friend B watches a drama." Friend B says, "I'll watch a drama if Friend C watches a horror." Friend C says, "I'll watch horror if Friend A watches a comedy."
• The Solution: Instead of guessing, the authors solved this loop simultaneously. They calculated how the "mass" (the weight or heaviness) of each particle changes based on the pressure (chemical potential), and how that change feeds back into the others. They did this mathematically at the "one-loop" level, which is like checking the first round of feedback in the conversation.

3. The Discovery: The "Snap" (First-Order Phase Transition)

The most exciting result is what happens when the pressure (chemical potential) reaches a specific critical point.

• The Scenario: As they increased the pressure, the "chiral condensate" (the glue holding the quarks together) stayed strong and steady. It was like a rubber band being stretched but not breaking.
• The Snap: Suddenly, at a specific pressure point (where the chemical potential equals the mass of a quark in a vacuum), the rubber band didn't just stretch; it snapped.
• The Result: The glue vanished instantly. The quarks, which were heavy and stuck together, suddenly became light and free. The "Pions" and "Sigmas" (the messengers and glue) swapped their behaviors and became identical twins.

This is called a First-Order Phase Transition.

• Analogy: Think of water turning into ice. Usually, it happens gradually. But a first-order transition is like a light switch. You flip it, and snap—it's either fully on or fully off. There is no "half-on" state. The paper found that the transition from "confined matter" to "free quark matter" happens exactly like flipping a switch, not like melting ice cream.

4. The Evidence: The Sound of the Switch

How do they know it was a snap and not a slow melt? They looked at the Speed of Sound.

• The Analogy: Imagine a crowd of people. If they are all holding hands tightly (confined matter), a rumor (sound) travels very slowly because everyone is stuck. If they let go and run freely (deconfined quark matter), the rumor travels much faster.
• The Finding: The authors calculated the "speed of sound" in this material.
  • Before the switch: The speed was zero (because the matter was rigid and static).
  • At the switch: The speed jumped instantly to a new value.
  • After the switch: The speed slowly climbed higher and higher, approaching a theoretical limit (the "conformal limit"), which represents a gas of free, massless particles moving at the speed of light.

This sudden jump in the speed of sound is the "smoking gun" that proves the transition was a sharp, first-order event.

Summary: Why Does This Matter?

This paper tells us that if you squeeze nuclear matter hard enough at zero temperature, it doesn't slowly turn into a soup of free quarks. Instead, it undergoes a violent, sudden transformation.

• For the Universe: This helps us understand the cores of Neutron Stars. These stars are so dense that they might be sitting right on the edge of this "snap." If the transition is sudden, it could cause the star to collapse or explode in a specific way.
• For the Future: The authors have built a new, more accurate tool (the self-consistent loop) that avoids the "lazy" approximations of the past. They plan to use this to map out the entire "map" of the universe's phase diagram, including what happens when you add heat and magnetic fields.

In short: The universe has a secret "light switch" for matter. Squeeze it hard enough, and snap—the rules of the game change instantly.

Technical Summary

1. Problem Statement

The paper addresses the nature of the chiral phase transition in Quantum Chromodynamics (QCD) at finite baryon chemical potential ($\mu_B$) and zero temperature ($T=0$).

• Context: While the transition at low density is known to be a smooth crossover, the behavior at high density is uncertain. It is hypothesized that a critical end point (CEP) exists where the crossover turns into a first-order phase transition.
• Limitations of Previous Approaches:
  • Lattice QCD: Cannot probe this region due to the "sign problem."
  • Effective Models (Standard): Previous applications of the Linear Sigma Model with quarks (LSMq) often relied on the ring-diagram approximation (static self-energies). This approximation is valid only at high temperatures and fails to accurately describe the system at low $T$ and arbitrary $\mu$.
  • Mass Instability: In standard tree-level treatments, as chiral symmetry is restored, the pion mass squared can become negative, signaling a breakdown of the classical description.
• Goal: To compute the particle pole masses self-consistently beyond the static approximation to determine the order of the phase transition and the behavior of thermodynamic quantities at $T=0$.

2. Methodology

The authors employ the Two-Flavor Linear Sigma Model with Quarks (LSMq) as an effective description of QCD. The calculation is performed at the one-loop order using a fully self-consistent approach.

• Lagrangian: The model includes $N_f=2$ quark flavors, $N_c=3$ colors, a scalar sigma field ($\sigma$), and a pion triplet ($\vec{\pi}$). The Lagrangian incorporates spontaneous chiral symmetry breaking and explicit breaking (via PCAC).
• Self-Consistent Pole Masses:
  • Unlike previous works that used fixed tree-level masses, the authors treat the masses ($M_\pi, M_\sigma, M_f$) as dynamical, $\mu$-dependent quantities.
  • They derive a system of three coupled integral equations for the pole masses. These equations arise from the Dyson series, where the self-energy of each particle depends on the masses of the others.
  • The calculation goes beyond the static ring-diagram approximation by including diagrams with external momentum dependence.
• Effective Potential:
  • The chiral order parameter (vacuum expectation value $v$) is determined by minimizing the one-loop effective potential ($V_{eff}$).
  • The potential is separated into vacuum (renormalized) and matter contributions. At $T=0$, only the fermion matter contribution survives for $\mu > M_f$.
• Numerical Implementation:
  • The system of coupled equations (masses + minimization condition) is solved numerically.
  • The authors use an iterative algorithm starting from high $\mu$ (where solutions are stable) and working backward, using random seeds to ensure global convergence.
  • Principal value integrals are handled carefully to manage singularities near the pole.

3. Key Contributions

1. Beyond Static Approximation: The paper implements a rigorous one-loop calculation for particle masses at $T=0$ that includes dynamic momentum dependence, overcoming the limitations of the high-temperature ring-diagram approximation.
2. Self-Consistent Framework: It establishes a framework where particle masses, couplings, and the chiral condensate are determined simultaneously, eliminating the unphysical issue of negative squared pion masses near the transition.
3. Thermodynamic Consistency: The method naturally yields thermodynamic quantities (pressure, density, speed of sound) directly from the effective potential without ad-hoc assumptions.

4. Key Results

The numerical solution reveals a first-order phase transition occurring at a specific critical chemical potential.

• Critical Chemical Potential: The transition occurs when the quark chemical potential reaches the vacuum quark mass:
  $$\mu_c = m_{f}^{vac} \approx 310 \text{ MeV}$$
  (Note: $\mu_B = 3\mu$).
• Discontinuous Behavior (First-Order Signatures):
  • Chiral Condensate ($v$): Drops discontinuously from its vacuum value ($\approx 92$ MeV) to a lower value at $\mu_c$, then decreases smoothly as $\mu$ increases further.
  • Particle Masses: The quark mass ($M_f$) and meson masses ($M_\pi, M_\sigma$) exhibit a discontinuous jump at $\mu_c$.
    • $M_f$ decreases monotonically after the transition, approaching zero (chiral restoration).
    • $M_\pi$ and $M_\sigma$ become degenerate at high densities, signaling the restoration of chiral symmetry.
  • Couplings and Parameters: The Yukawa coupling ($g$) and quartic coupling ($\lambda$) show discontinuities. Crucially, the mass parameter $a^2$ (which controls the curvature of the potential) changes sign from positive to negative, indicating the transition from the broken phase (Mexican hat potential) to the symmetric phase (parabolic potential).
• Thermodynamics:
  • Speed of Sound ($c_s^2$): At low density ($\mu < \mu_c$), $c_s^2 = 0$ (hadronic matter). At the transition, it jumps discontinuously to a finite value. As $\mu$ increases, $c_s^2$ rises smoothly and approaches the conformal limit ($1/3$) from below, consistent with the emergence of a gas of non-interacting relativistic massless quarks.

5. Significance

• Validation of First-Order Transition: The study provides strong theoretical evidence within the LSMq framework that the chiral phase transition at zero temperature is of the first order, occurring precisely when the chemical potential equals the constituent quark mass.
• Methodological Advancement: By solving the self-consistent pole mass equations, the authors demonstrate a robust method to handle the interplay between medium effects and chiral symmetry restoration, avoiding the pathologies of tree-level approximations.
• Phenomenological Implications: The results offer specific predictions for the equation of state (EoS) of dense matter. The discontinuity in the speed of sound and the specific value of $\mu_c$ are relevant for interpreting data from heavy-ion collision experiments (like STAR BES-I/II and the future NICA/MPD) and for modeling the interiors of neutron stars.
• Future Outlook: The framework is designed to be extended to finite temperatures ($T > 0$) to map the full QCD phase diagram and locate the Critical End Point (CEP), as well as to include isospin chemical potentials.

In summary, this work refines the understanding of dense QCD matter by applying a self-consistent, one-loop effective field theory approach, confirming a first-order chiral phase transition at zero temperature and providing a consistent thermodynamic description of the transition region.