⚛️ nuclear theory

Quark and hybrid stars with renormalization group improvement of NNLO perturbative QCD

This paper extends renormalization-group-optimized perturbative QCD to $\beta$ -equilibrated matter to derive a compact equation of state for massive quarks, demonstrating that specific renormalization scale parameters can yield stable pure and hybrid quark stars compatible with recent astrophysical observations like PSR J0740+6620 and GW190814.

Original authors: Loïc Fernandez, Jean-Loïc Kneur, Marcus Benghi Pinto, Constança Providência, Claudia Ratti, Tulio E. Restrepo

Published 2026-06-30

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Loïc Fernandez, Jean-Loïc Kneur, Marcus Benghi Pinto, Constança Providência, Claudia Ratti, Tulio E. Restrepo

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 filled with incredibly dense objects called neutron stars. These are the leftover cores of massive stars that have collapsed, packing more mass than our Sun into a city-sized sphere. Inside these stars, the pressure is so immense that the usual building blocks of matter (protons and neutrons) might get crushed so hard that they dissolve into a "soup" of their even smaller components: quarks.

The paper you are asking about is a team of physicists trying to figure out exactly what happens inside this "quark soup" and how it affects the size and weight of these stars. They are using a specific mathematical tool called Quantum Chromodynamics (QCD), which is the rulebook for how quarks interact.

Here is a simple breakdown of their work:

1. The Problem: A Wobbly Rulebook

The physicists are trying to calculate the "Equation of State" (EoS). Think of the EoS as a recipe that tells you how much a star will squish under its own weight.

The Issue: When they use the standard version of the QCD rulebook (called "perturbative QCD"), the recipe changes depending on how you hold the measuring cup. In physics terms, the results change too much based on an arbitrary number called the "renormalization scale." It's like trying to bake a cake, but the recipe says the amount of flour depends on whether you measure it in the morning or the afternoon. This makes the predictions unreliable.
The Goal: They want a recipe that stays the same no matter how you measure it.

2. The Solution: A Better Measuring Tool (RGOPT)

To fix the wobbly measurements, the authors used a new method called Renormalization Group Optimized Perturbation Theory (RGOPT).

The Analogy: Imagine you are trying to tune a radio to get a clear signal, but there is a lot of static (noise). The old method (standard QCD) leaves you with a lot of static that changes depending on how you turn the dial. The new method (RGOPT) is like a smart antenna that automatically adjusts itself to cancel out the static, giving you a much clearer, more stable signal.
The Result: By using RGOPT, they created a "cleaner" version of the quark soup recipe. This new recipe is much less sensitive to the arbitrary measuring scales, making the predictions for how dense matter behaves much more trustworthy.

3. The "Pocket Formula"

The math behind this new recipe is incredibly complex, involving thousands of pages of equations. To make it useful for other scientists, the authors created a "pocket formula."

The Analogy: Think of the full math as a massive, heavy encyclopedia. The "pocket formula" is a handy cheat sheet or a smartphone app that gives you the exact same answer without needing to carry the whole encyclopedia. This allows other researchers to easily plug in numbers and calculate how these stars behave.

4. Testing the Stars: Pure Quark Stars vs. Hybrid Stars

The team used their new recipe to build two types of theoretical stars and compared them against real observations from telescopes and gravitational wave detectors.

A. Pure Quark Stars
These are stars made entirely of the "quark soup."

The Finding: To make these stars heavy enough to match what we see in the sky (some are twice as heavy as our Sun), the "clean" RGOPT recipe required a specific setting (a parameter called $X$ ) to be between 3.08 and 3.58.
Comparison: The old, wobbly recipe (standard QCD) needed a much tighter and lower range of settings to match the same stars. This suggests the old recipe was less stable and perhaps less accurate.

B. Hybrid Stars
These are stars with a normal "crust" of protons and neutrons, but a core of "quark soup."

The Experiment: They tested three different models for the "crust" (some soft and squishy, some hard and stiff) and combined them with their new quark recipe.
The Finding: They found that hybrid stars with a small quark core (just a tiny bit of soup in the middle) work well with the observations. However, if they tried to make the quark core very large (like a giant ball of soup in the center), the star would become too big or unstable to match what we see.
The GW190814 Mystery: There was a recent observation of a very light, compact object (possibly a neutron star) colliding with a black hole. If that object was a neutron star, the authors found that the "pure quark star" model would need a very specific, high setting ( $X \approx 4.1$ ) to fit.

5. The Bottom Line

The paper claims that by using this new, "cleaner" mathematical method (RGOPT), they have created a more reliable way to describe the extreme matter inside neutron stars.

Their new "pocket formula" allows scientists to easily calculate how these stars behave.
They found that pure quark stars are possible but require specific conditions to match real-world data.
Hybrid stars (with a small quark core) are very compatible with current observations, especially for massive pulsars like PSR J0740+6620.
The method significantly reduces the "guesswork" (scale dependence) that plagued previous calculations, giving astronomers a more solid foundation for understanding the densest objects in the universe.

In short, they didn't just find a new star; they built a better ruler to measure the stars we already know exist.

Problem Statement
The theoretical description of the Equation of State (EoS) for cold, dense quark matter is essential for modeling the interiors of neutron stars (NS) and hybrid stars. While perturbative Quantum Chromodynamics (pQCD) provides an ab-initio framework for high baryon densities, its application at the intermediate densities relevant to neutron star cores ( $\sim 5-10 \rho_0$ ) is hindered by a significant sensitivity to the arbitrary renormalization scale ( $\Lambda$ ). Specifically, the Next-to-Next-to-Leading Order (NNLO) pQCD pressure exhibits a large residual scale dependence, primarily originating from the "hard" sector contributions. This uncertainty complicates the determination of whether quark matter exists in the cores of observed massive pulsars (e.g., PSR J0740+6620) and constrains the nature of the hadron-quark phase transition. Furthermore, existing lattice QCD calculations are limited to low chemical potentials due to the sign problem, and effective models often lack asymptotic freedom.

Methodology
The authors extend the Renormalization Group Optimized Perturbation Theory (RGOPT) framework, previously applied to symmetric matter, to the case of non-symmetric, $\beta$ -equilibrated, and charge-neutral quark matter with $N_f = 2+1$ massive quarks.

RGOPT Framework: Instead of a standard weak-coupling expansion, RGOPT introduces a modified perturbative expansion around massive quasiparticle states. This involves a "RG-dressed" screening mass ( $\eta$ ) that incorporates all-order Renormalization Group (RG) properties. The method utilizes a specific substitution of parameters (Eq. 2.3) and fixes the arbitrary mass parameter $\eta$ by solving the RG equation (Eq. 2.6) alongside a stationarity condition, ensuring the solution respects asymptotic freedom and remains real-valued through a controlled renormalization scheme change (RSC).
Thermodynamic Consistency: To apply the theory to compact stars, the authors enforce $\beta$ -equilibrium and charge neutrality. They address the issue of thermodynamic inconsistency (where $dP/d\mu \neq \rho$ ) in the raw RGOPT pressure by integrating the difference between the partial and total derivatives of the pressure (Eq. 3.4).
Compact Formulas: Recognizing the complexity of the full NNLO RGOPT expression, the authors derive "pocket" fitting formulas for the thermodynamically consistent pressure and energy density. These formulas are valid for two distinct ranges of the renormalization scale parameter $X$ (where $\Lambda = X \mu_B/3$ ): $2 \le X \le 3$ (relevant for hybrid stars) and $3 \le X \le 6$ (relevant for pure quark stars).
Hybrid Star Construction: To model hybrid stars, the authors employ a Maxwell construction to describe a first-order phase transition between hadronic and quark phases. They combine the RGOPT quark EoS with three representative hadronic models based on Relativistic Mean Field (RMF) theory: DD2 (stiff), SFHo (soft), and DIDY (soft, including hyperons).

Key Contributions

Extension of RGOPT to $\beta$ -equilibrium: The paper provides the first NNLO RGOPT EoS for cold, dense matter with $N_f=2+1$ massive quarks in $\beta$ -equilibrium and charge neutrality, a necessary step for realistic stellar modeling.
Reduction of Scale Uncertainty: The study demonstrates that the RGOPT pressure is significantly less sensitive to variations in the renormalization scale parameter $X$ compared to standard NNLO pQCD.
Analytical Tools: The authors provide explicit, compact fitting formulas (Eqs. 3.9 and 3.13) for the RGOPT pressure and baryon density, enabling immediate application in astrophysical simulations without recalculating the complex perturbative integrals.
Systematic Comparison: A comprehensive comparison is made between pure quark stars and hybrid stars using both RGOPT and standard pQCD, constrained by recent astrophysical observations.

Results

Pure Quark Stars (QS):
- To reproduce maximum masses of $2.0, 2.3,$ and $2.6 M_\odot$ , the RGOPT requires scale parameters in the range $X \approx 3.08 - 3.58$ . In contrast, standard pQCD requires a narrower and lower range ( $X \approx 2.95 - 3.26$ ) and exhibits a stronger dependence of maximum mass on $X$ .
- If the low-mass object in the GW190814 event is interpreted as a neutron star, RGOPT requires $X \approx 4.10$ , while pQCD requires $X \approx 3.56$ .
- RGOPT predictions for pure quark stars are more stable against scale variations, with the maximum mass changing by $\sim 0.59 M_\odot$ per unit of $X$ , compared to $\sim 0.98 M_\odot$ for pQCD.
- Both approaches satisfy the Bodmer-Witten hypothesis for stable strange matter.
Hybrid Stars:
- Hybrid stars were constructed using the Maxwell construction with the three hadronic models.
- For hybrid stars compatible with the massive pulsar PSR J0740+6620, the RGOPT scale parameter $X$ ranges from approximately $1.99$ to $2.98$ depending on the hadronic model.
- The largest $X$ values yield hybrid stars with large quark cores (radii of $5$ to $8$ km). The smallest $X$ values yield stars with very small quark cores.
- The transition density and the strength of the phase transition depend on the hadronic model. For the stiff DD2 model, the transition occurs at lower densities ( $\sim 2 \rho_0$ ) compared to the soft SFHo and DIDY models ( $\sim 4 \rho_0$ ).
- Hybrid stars with large quark cores (compatible with PSR J0740+6620) generally have hadronic shells of $3-5$ km and are inconsistent with the low-mass compact object HESS J1731-347. Conversely, hybrid stars with small quark cores (small $X$ ) are compatible with HESS J1731-347.
- The speed of sound in the hybrid EoSs exceeds the conformal limit ( $\sqrt{1/3}$ ) in the hadronic branch but drops below it in the quark phase after the transition.

Significance and Claims
The paper claims that the RGOPT framework offers a more robust theoretical tool for describing the quark sector of compact stars than standard pQCD at the same perturbative order. By incorporating RG properties directly into the perturbative expansion, RGOPT mitigates the large residual scale dependence that plagues standard pQCD predictions at intermediate densities. This reduction in theoretical uncertainty allows for tighter constraints on the renormalization scale when matching to astrophysical data. The authors assert that their derived compact formulas provide a practical and accurate means to incorporate these improved QCD results into stellar structure calculations, facilitating the investigation of whether quark matter cores exist in observed neutron stars. The work highlights that while both RGOPT and pQCD can accommodate current observational constraints, RGOPT provides a more stable theoretical foundation for exploring the high-density equation of state.

1. The Problem: A Wobbly Rulebook

2. The Solution: A Better Measuring Tool (RGOPT)

3. The "Pocket Formula"

4. Testing the Stars: Pure Quark Stars vs. Hybrid Stars

5. The Bottom Line

More like this