⚛️ phenomenology

Comparing critical regions in Polyakov-quark-meson models: Effects of flavor (two versus 2+1) under on-shell renormalization, and on-shell versus curvature-mass parameter fixing for the two-flavor case

This paper compares the critical regions around the chiral critical end point in two-flavor and 2+1-flavor Polyakov-quark-meson models, analyzing how different treatments of vacuum fluctuations (on-shell renormalization versus curvature-mass parametrization) and flavor content influence the extent and location of enhanced quark number susceptibility.

Original authors: Pooja Kumari, Suraj Kumar Rai, Vivek Kumar Tiwari

Published 2026-09-24

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Pooja Kumari, Suraj Kumar Rai, Vivek Kumar Tiwari

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

Deep within the heart of every atom, protons and neutrons are held together by a force so powerful it defies our everyday intuition. This force, known as the strong interaction, is carried by particles called gluons and acts upon even smaller building blocks called quarks. Under normal conditions, quarks are forever trapped inside their parent particles, unable to roam free. However, scientists believe that if you heat matter up to temperatures trillions of degrees hotter than the center of the sun, or squeeze it to densities far beyond what any star can withstand, this confinement breaks down. The protons and neutrons melt away, dissolving into a primordial soup of free-moving quarks and gluons known as the quark-gluon plasma. This state of matter is thought to have filled the entire universe just microseconds after the Big Bang, and it is recreated for fleeting moments in massive particle accelerators where scientists smash heavy ions together at nearly the speed of light.

Understanding exactly how matter transforms from solid particles into this fluid plasma is one of the great challenges of modern physics. A key question is whether this change happens smoothly, like ice melting into water, or abruptly, like water boiling into steam. Theoretical models suggest that as you change the temperature and the density of the matter, there is a specific point where the nature of the transition shifts. At this special location, called the critical end point, the smooth change turns into a sharp, explosive one. Finding this point is crucial because it would tell us how the early universe evolved and help us understand the dense cores of neutron stars. However, because the conditions inside these stars and the early universe are so extreme, we cannot observe the critical end point directly. Instead, physicists must rely on complex mathematical models to predict where it lies and what the surrounding landscape looks like.

In a recent study, researchers set out to refine these predictions by testing how different ways of handling the invisible quantum fluctuations of particles affect the map of this critical region. They focused on a specific type of model that describes how quarks interact with other particles called mesons. To make their calculations more realistic, they included the effects of the "vacuum," a seething sea of virtual particles that pop in and out of existence even in empty space. The team compared two different mathematical approaches to fixing the parameters of their models: one that treats particle masses as they are measured in experiments, and another that defines them based on how the energy of the system curves at its lowest point. They also tested how adding a third type of quark, the strange quark, changes the picture, and how including a mechanism that mimics the confinement of quarks alters the results.

The researchers discovered that the way they handled these vacuum fluctuations made a dramatic difference in the shape and size of the critical region. When they used the approach based on measured particle masses, the area where the critical fluctuations occur was smaller and shifted to higher temperatures and lower densities. In contrast, the method based on the curvature of the energy landscape produced a much larger, more spread-out region that sat at lower temperatures and higher densities. This difference is significant because it changes where scientists should look for the critical end point in their experiments. The study also revealed that the presence of the strange quark acts as a dampener. When the third flavor of quark was included, the critical region shrank considerably, becoming narrower and less extended. This suggests that the extra complexity of the real world, with its three types of light quarks, makes the critical fluctuations harder to detect than in simpler, two-quark scenarios.

A particularly interesting finding involved the relationship between the critical end point and another theoretical feature called the tricritical point. In the simpler two-quark models, the tricritical point sits comfortably inside the critical region, meaning it strongly influences the behavior of the system near the critical end point. However, when the third quark was added, this relationship changed. The tricritical point moved outside the critical region in some models, effectively cutting it off from influencing the critical fluctuations. This separation implies that the physics governing the critical end point in the real world might be cleaner and less complicated by other nearby phase transitions than previously thought in simpler models.

The team also calculated how the sensitivity of the system to changes in density behaves as it approaches the critical point. They found that the rate at which this sensitivity grows follows a specific mathematical pattern, known as a power law. By measuring the slope of this growth, they determined a critical exponent, a number that describes the universality of the transition. Their results showed that the exponent values were consistent with theoretical expectations for a system undergoing a phase transition, though they varied slightly depending on the direction from which the critical point was approached. These calculations provide a more precise set of numbers for experimentalists to compare against their data from heavy-ion collision experiments.

Ultimately, this work highlights that the details of how we account for the quantum vacuum are not just mathematical technicalities; they fundamentally reshape our understanding of where the critical end point might be hiding. The study suggests that if we want to find this elusive point in the real world, we must look in a region of temperature and density that is higher and less dense than some previous models predicted. Furthermore, the inclusion of the strange quark and the specific way quarks are confined suggests that the critical region is more compact than earlier simulations indicated. These refined maps give the global scientific community a clearer, more targeted guide for the ongoing search for the critical end point, bringing us one step closer to understanding the fundamental nature of matter under the most extreme conditions in the universe.

Technical Summary: Comparing Critical Regions in Polyakov-Quark-Meson Models

Problem Statement
The paper addresses the theoretical uncertainties surrounding the Critical End Point (CEP) and Tricritical Point (TCP) in the Quantum Chromodynamics (QCD) phase diagram. Specifically, it investigates how different treatments of quark one-loop vacuum fluctuations and the inclusion of the strange quark flavor (2 vs. 2+1 flavors) affect the morphology, size, and location of critical fluctuation regions. Previous studies often neglected vacuum fluctuations or used inconsistent parameter-fixing schemes (e.g., matching curvature masses derived from the effective potential while the potential itself contains vacuum corrections). Furthermore, the interplay between the axial $U_A(1)$ anomaly, the strange quark mass, and the choice of Polyakov-loop potential (with or without quark back-reaction) on the critical region's topography requires systematic clarification to aid experimental searches at facilities like RHIC-BES, FAIR, and NICA.

Methodology
The authors employ a comparative framework using four distinct two-flavor models and their extensions to 2+1 flavors:

Renormalized Quark-Meson (RQM) and Renormalized Polyakov-Quark-Meson (RPQM) Models: These utilize an on-shell (OS) renormalization scheme where counter-terms are matched between the OS and Modified Minimal Subtraction (MS) schemes. Parameters are fixed using physical pole masses ( $m_\pi, m_\sigma$ ) and the pion decay constant ( $f_\pi$ ).
Quark-Meson with Vacuum Term (QMVT) and Log-PQMVT Models: These use a curvature-mass-based parameter fixing scheme, where model parameters are determined by the second derivatives of the effective potential at its minimum.
Polyakov-Loop Potentials: The study compares two forms of the Polyakov-loop potential:
- Logarithmic (Log): Used in Log-RPQM and Log-PQMVT, representing confinement-deconfinement without quark back-reaction.
- PolyLog-glue: Used in PolyLog-glue-RPQM, incorporating quark back-reaction effects on the gluon sector.

The authors compute phase diagrams for both physical pion mass ( $m_\pi = 138$ MeV) and the chiral limit ( $m_\pi = 0$ ). They map critical regions by calculating contours of the enhanced normalized quark number susceptibility ( $R_q = \chi_q / \chi_q^{free}$ ) and scalar susceptibility ( $R_s = \chi_s / \chi_s^{vac}$ ). The study also computes critical exponents for the power-law divergence of $\chi_q$ near the CEP and TCP.

Key Contributions and Results

Impact of Renormalization Scheme (On-Shell vs. Curvature Mass):
- In the RQM model, the $\sigma$ meson curvature mass ( $m_{\sigma,c} \approx 417$ MeV) differs from its pole mass ($500$ MeV) due to on-shell renormalization. This results in a CEP located at higher temperature ( $T \approx 35.7$ MeV) and lower chemical potential ( $\mu \approx 276.7$ MeV) compared to the QMVT model.
- The QMVT model, where curvature and pole masses coincide, exhibits a stronger smoothing of the chiral transition, pushing the CEP to lower $T$ ($28.4$ MeV) and higher $\mu$ ($298.9$ MeV).
- Consequently, the critical region contours ( $R_q=2$ ) in the QMVT model are broader and more extended in the $T$ - $\mu$ plane than in the RQM model. The QMVT contours exhibit a "neck-like" structure below the CEP due to the sharp curvature of the phase boundary at high $\mu$ , whereas the RQM contours are more compact and lack this feature.
Effect of Polyakov-Loop Potentials:
- Integrating the Logarithmic potential shifts CEPs to significantly higher temperatures.
- The Log-RPQM model (no back-reaction) produces a CEP at $(270.6, 92.5)$ MeV with compressed, pinched critical contours.
- The Log-PQMVT model (curvature mass + Log potential) yields a CEP at $(296.1, 74.5)$ MeV. The critical region here is highly stretched along the $T$ -axis but sharply compressed in width, forming a pronounced neck-like structure.
- The PolyLog-glue-RPQM model (with back-reaction) moderates these effects. The back-reaction smooths the phase boundary, resulting in broader contours compared to the Log-RPQM case, though they remain narrower than the QMVT counterparts.
Flavor Dependence (2 vs. 2+1 Flavors):
- Introducing the strange quark (2+1 flavors) universally reduces the size and width of the critical regions across all models.
- In the RQM model, the 2+1 flavor contours shrink proportionally in both $T$ and $\mu$ directions, maintaining a broad geometry.
- In the Log-RPQM and PolyLog-glue-RPQM models, the addition of the third flavor causes a drastic compression in the $T$ -direction (pinching), while the $\mu$ -direction reduction is more moderate.
- TCP Influence: In all two-flavor scenarios, the Tricritical Point (TCP) lies well inside the $R_q=2$ contour of the CEP, indicating that the TCP influences critical fluctuations near the CEP. In contrast, for the 2+1 flavor cases, the TCP moves outside the $R_q=2$ contour (or sits on its boundary), effectively neutralizing its influence on the CEP fluctuations.
Critical Exponents:
- The critical exponent for the divergence of $\chi_q$ at the TCP in the RQM model is calculated as $\gamma \approx 0.54$ , consistent with the mean-field prediction and previous PTRG results.
- At the CEP, the exponents for the divergence of $\chi_q$ $χ_{q}$ as $\mu \to \mu_{CEP}$ $μ \to μ_{C E P}$ vary by direction and model:
  - RQM: $\epsilon \approx 0.706$ (from below) and $0.656$ (from above).
  - Log-RPQM: $\epsilon \approx 0.718$ (from below) and $0.655$ (from above).
  - PolyLog-glue-RPQM: $\epsilon \approx 0.756$ (from below) and $0.689$ (from above).
- These values bracket the mean-field field-like exponent ( $2/3 \approx 0.67$ ), with deviations attributed to the mixing of scaling axes and the specific vacuum fluctuation treatments.

Significance
The paper claims that the treatment of quark one-loop vacuum fluctuations and the choice of renormalization scheme significantly alter the predicted location and morphology of the QCD critical region. The on-shell renormalized models (RQM/RPQM) predict CEPs at higher temperatures and lower chemical potentials with more compact critical regions compared to curvature-mass-based models (QMVT/PQMVT). Furthermore, the inclusion of the strange flavor and the $U_A(1)$ anomaly systematically suppresses critical fluctuations and shifts the TCP away from the CEP, potentially reducing the TCP's influence on experimental observables. The study provides essential theoretical boundaries for interpreting experimental signatures of the CEP across different collision energies, highlighting that the "shape" of the critical region is highly sensitive to the underlying effective theory assumptions regarding vacuum corrections and flavor content.

More like this