⚛️ phenomenology

Effects of Axion Interactions on Quark Stars in 4D Einstein-Gauss-Bonnet Gravity

This paper demonstrates that combining axion-induced interactions within the Polyakov-Nambu-Jona-Lasinio model and four-dimensional Einstein-Gauss-Bonnet gravity stiffens the quark matter equation of state, resulting in quark stars with higher maximum masses, larger radii, and reduced compactness that better align with current observational constraints than standard general relativity models.

Original authors: Rui Zhou, Ming-zheng-xuan Wu, Xin-ran Yang, Chong-long Xie, Min Zhou, Zhi-yang Liu, Shi-jun Mao, Guo-yun Shao

Published 2026-07-28

📖 6 min read🧠 Deep dive

CC BY 4.0

Original authors: Rui Zhou, Ming-zheng-xuan Wu, Xin-ran Yang, Chong-long Xie, Min Zhou, Zhi-yang Liu, Shi-jun Mao, Guo-yun Shao

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 as a giant cosmic playground where the rules of physics get stretched to their absolute limits. Deep inside the hearts of dead stars, called neutron stars or quark stars, matter is squeezed so tightly that atoms can't even hold their shape. It's like trying to pack an entire city into a single sugar cube. For decades, scientists have used a rulebook called General Relativity to predict how these stars behave, but that rulebook might be missing a few pages. There are also whispers of invisible particles called "axions" that might be hiding in the dark matter of the universe, and tiny, ghostly forces that could change how these stars hold themselves together. The big question is: if we tweak the rules of gravity and add these hidden particles to the mix, do these stars get heavier, bigger, or do they collapse?

This paper takes a wild ride through that question by building a new model of a "quark star." The authors combine two big ideas: a theory of gravity that adds a little extra "curvature" to the fabric of space (called 4D Einstein-Gauss-Bonnet gravity) and a microscopic model where axions interact with the soup of quarks inside the star. They don't just guess; they run complex computer simulations to see how these changes affect the star's size and weight. Their results suggest a fascinating possibility: when you mix these specific gravity tweaks with axion interactions, quark stars might be able to support much more mass than we thought possible, potentially solving a puzzle about why we see some incredibly heavy stars in the sky that shouldn't exist according to the old rulebook.

The Cosmic Stress Test

To understand what the authors did, let's imagine a quark star as a super-dense, self-contained ball of cosmic dough. In the standard version of our universe (General Relativity), gravity is like a relentless, invisible hand squeezing this dough. The dough pushes back with its own internal pressure. If the gravity gets too strong, the dough collapses into a black hole. If the dough is stiff enough, it holds its shape as a star.

The authors decided to test two new ingredients in this cosmic recipe. First, they added a "gravity modifier" known as 4D Einstein-Gauss-Bonnet (EGB) gravity. You can think of this as adding a little bit of "elasticity" to the universe's fabric. In this modified gravity, the invisible hand of gravity doesn't squeeze quite as hard as it does in the standard model. It's like the star is wearing a slightly more comfortable, stretchy suit that lets it expand a bit more before the pressure becomes too much.

Second, they introduced "axions." Imagine axions as tiny, invisible ghosts that float around inside the star's dough. These ghosts interact with the quarks (the tiny particles making up the dough) and change how the dough behaves. Specifically, the authors found that these axion interactions make the dough "stiffer." In everyday terms, it's like switching from soft, squishy jelly to a firm, rubbery gel. A stiffer material is much harder to compress, meaning it can resist the crushing weight of gravity much better.

The Results: Bigger, Heavier, and Less Squished

When the authors ran their simulations, the results were quite exciting. They found that both of these new ingredients work together to make the stars more robust.

First, the "stiffening" effect of the axions means the star can hold more weight. In their model, as they increased the strength of the axion interaction (represented by a parameter called $\theta$ ), the maximum weight the star could hold went up. For example, with no axion effects, the star could hold about 1.67 times the mass of our Sun ( $1.67 M_\odot$ ). But with the strongest axion effects they tested, the star could hold up to 1.94 times the mass of our Sun ( $1.94 M_\odot$ ). That's a significant jump!

Second, the modified gravity (EGB) acted like a safety net. Because it weakens the gravitational squeeze, the star doesn't need to be as tightly packed to stay stable. This allowed the stars to grow larger in size (radius) while still holding massive amounts of weight. The authors noted that this combination leads to stars that are "less compact," meaning they are bigger and fluffier than the super-dense, tiny stars predicted by standard gravity.

Checking Against the Real Universe

The real test for any theory is whether it matches what we actually see in the sky. The authors compared their new, super-stiff, gravity-modified stars against real observations of massive pulsars (spinning neutron stars) like PSR J0740+6620 and PSR J1614-2230. These real stars are incredibly heavy, and standard models sometimes struggle to explain how they don't collapse.

The paper suggests that when you combine the axion effects with the EGB gravity correction, the resulting quark stars fit the observational data much better. Specifically, for certain values of the axion parameter (like $\theta \ge \pi/3$ ) and the gravity coupling (like $\alpha = 6.0 \text{ km}^2$ ), the simulated stars can easily reach the heavy masses observed in PSR J0740+6620 without breaking the rules of physics. They also found that these stars would have a smaller "surface gravitational redshift," which is a fancy way of saying that light escaping from their surface wouldn't be stretched out as much as it would be from a standard, super-dense star.

The Takeaway

In short, this paper suggests that the universe might be playing by slightly different rules than we thought. By mixing a modified theory of gravity with the potential influence of axion particles, the authors propose a new way to explain how massive quark stars can exist. They aren't claiming to have proven that these stars definitely work this way, but their simulations suggest that this combination of factors provides a very plausible explanation for the heavy, stable stars we observe. It's a reminder that in the extreme depths of the cosmos, the interplay between invisible particles and the curvature of space might just be the key to unlocking the secrets of the heaviest objects in the universe.

Technical Summary: Effects of Axion Interactions on Quark Stars in 4D Einstein-Gauss-Bonnet Gravity

Problem Statement
The internal structure and global properties of compact stars, particularly quark stars, are governed by the equation of state (EOS) of dense matter and the underlying theory of gravity. Two significant open questions remain: the nature of microscopic interactions in dense quark matter (specifically the potential role of axions as dark matter candidates or solutions to the strong CP problem) and the validity of General Relativity (GR) in the strong-field regime. While axion interactions are known to influence the QCD phase structure, and 4D Einstein-Gauss-Bonnet (EGB) gravity offers a higher-curvature extension to GR, the combined phenomenological impact of axion-induced $\theta$ effects and 4D EGB corrections on the structure of quark stars has not been previously explored.

Methodology
The authors construct a theoretical framework combining a microscopic model of quark matter with a macroscopic modified gravity theory:

Microscopic Model (Axion-Extended PNJL): The study employs a $(2+1)$ -flavor Polyakov–Nambu–Jona-Lasinio (PNJL) model extended to include axion interactions. The Lagrangian incorporates a $\theta$ -dependent 't Hooft term, where $\theta = a/f_a$ represents the effective axion-induced CP-violating phase. The model includes scalar, vector, and instanton-induced interactions. The thermodynamic potential is minimized to determine quark condensates and the Polyakov loop, yielding the pressure ( $P$ ) and energy density ( $\epsilon$ ) for cold, $\beta$ -equilibrated, charge-neutral quark matter.
Macroscopic Framework (4D EGB Gravity): The stellar structure is modeled using the regularized 4D EGB gravity, formulated as an effective scalar-tensor theory of Horndeski type. The authors utilize the modified Tolman–Oppenheimer–Volkoff (TOV) equations derived from this theory, which include a Gauss-Bonnet coupling constant $\alpha$ .
Numerical Setup: The study solves the modified TOV equations for various values of the axion parameter $\theta$ ($0$ to $\pi$ ) and the Gauss-Bonnet coupling $\alpha$ (up to $6.0 \text{ km}^2$ , consistent with observational bounds). Key observables calculated include the mass-radius ( $M-R$ ) relations, maximum mass ( $M_{\text{max}}$ ), stellar radius ( $R$ ), compactness ( $C=M/R$ ), surface gravitational redshift ( $Z$ ), and central pressure/density profiles.

Key Contributions and Results

Effect of Axion Interactions on the EOS: The inclusion of axion interactions (increasing $\theta$ ) stiffens the EOS of quark matter. This is evidenced by an increase in pressure at fixed energy density and a rise in the squared sound speed ( $c_s^2$ ) and adiabatic index ( $\gamma$ ) in specific density regimes. Consequently, the maximum mass of quark stars increases, and the stellar radius expands. For instance, in standard GR ( $\alpha=0$ ), increasing $\theta$ from $0$ to $\pi$ raises $M_{\text{max}}$ from $1.67 M_\odot$ to $1.94 M_\odot$ .
Effect of 4D EGB Gravity on Stellar Structure: The Gauss-Bonnet correction ( $\alpha > 0$ ) effectively weakens gravitational compression. This allows quark stars to support larger masses before reaching the onset of instability. Unlike the axion effect which stiffens the matter, the EGB correction modifies the gravitational potential, leading to larger radii for a given mass. This results in a decrease in both compactness and surface gravitational redshift compared to standard GR predictions.
Combined Effects and Observational Compatibility: The interplay between axion-induced stiffening and 4D EGB gravitational corrections yields mass-radius sequences that align more closely with current astrophysical constraints than standard GR models with conventional EOS.
- The combined model allows quark stars to satisfy the mass constraints of massive pulsars PSR J1614-2230 and PSR J0740+6620 for specific parameter choices (e.g., $\alpha = 6.0 \text{ km}^2$ and $\theta \geq \pi/3$ ).
- These configurations remain compatible with radius constraints inferred from the gravitational wave event GW170817 and the X-ray source XMMU J173203.3-344518.
- The study identifies that the turning point of stability (onset of instability) shifts toward higher masses and, in the case of increasing $\theta$ , toward lower central energy densities.

Significance and Claims
The paper claims that the simultaneous consideration of axion dynamics and 4D EGB gravity provides a viable phenomenological framework for describing massive quark stars. The authors assert that while axion interactions and 4D EGB corrections operate through distinct mechanisms—stiffening the EOS versus weakening gravitational compression, respectively—their combined effect shifts theoretical predictions toward better agreement with observational data. Specifically, this framework offers a potential explanation for the existence of very massive compact stars without requiring exotic forms of matter beyond the standard quark description or violating current gravitational wave constraints. The authors conclude that future multimessenger observations could further constrain the model parameters $\alpha$ and $\theta$ .

The Cosmic Stress Test

The Results: Bigger, Heavier, and Less Squished

Checking Against the Real Universe

The Takeaway

More like this