Observational Probes of the Neutron Star Equation of State with Hyperons, Bosonic Dark Matter, and Quark Matter

This study demonstrates that a hybrid neutron star model incorporating hyperons, deconfined quark matter, and bosonic sexaquark dark matter with a mass of approximately 1900 MeV successfully satisfies all current observational constraints on mass, radius, and tidal deformability.

The Gist

Imagine the universe is filled with invisible, ghostly stuff called Dark Matter. We know it's there because of how it pulls on galaxies, but we've never actually seen a single particle of it. It's like trying to figure out what's in a sealed, heavy box just by weighing the box and shaking it.

For decades, scientists have been trying to guess what's inside that box. One of the most exciting new ideas is a particle called the Sexaquark. Think of it as a "super-tight knot" made of six quarks (the tiny building blocks of protons and neutrons) stuck together so tightly that it acts like a single, heavy, invisible marble.

This paper is a detective story about whether these "Sexaquark marbles" could be hiding inside Neutron Stars.

The Setting: The Ultimate Pressure Cooker

Neutron stars are the dead, collapsed cores of massive stars. They are so dense that a teaspoon of their material would weigh a billion tons. It's the most extreme pressure cooker in the universe.

Inside these stars, scientists have long suspected that matter gets so squished that it changes form:

1. The Outer Layer: Normal stuff (protons and neutrons).
2. The Middle: Exotic stuff (hyperons, which are heavier cousins of neutrons).
3. The Core: A soup of free-floating quarks (Deconfined Quark Matter).

The problem is that when scientists try to build a mathematical model of this "pressure cooker," the stars they predict are often too big and too stiff. But when we look at real neutron stars through telescopes (like the NICER telescope), they seem smaller and squishier than our models predict. It's like our recipe for a cake predicts a fluffy sponge, but the real cake is a dense brownie.

The New Ingredient: The Sexaquark

The authors of this paper asked: What if we add Sexaquarks to the recipe?

They imagined that inside the crushing pressure of a neutron star, these Sexaquark marbles form. Because they are bosons (a type of particle that likes to crowd together), they act like a "softener" for the star's interior.

The Analogy:
Imagine a crowd of people in a room (the neutrons). If you tell them to stand perfectly still and rigid, the room feels very stiff and hard to squeeze.
Now, imagine you introduce a group of people who are very flexible and like to huddle together in a tight, soft ball (the Sexaquarks). Suddenly, the whole room becomes easier to compress. The "stiffness" of the room drops.

In the paper, the authors found that adding these Sexaquarks makes the neutron star model softer, which matches the real observations of how small these stars actually are.

The Investigation: Testing the Recipe

The team didn't just guess; they ran a massive simulation using a supercomputer. They treated the neutron star like a giant puzzle with three pieces:

1. Hyperons (the heavy cousins).
2. Quark Soup (the free-floating core).
3. Sexaquarks (the new dark matter ingredient).

They tested different "weights" for the Sexaquark.

• Too Light? The star would collapse or behave strangely.
• Too Heavy? The star would be too stiff, and the model would fail to match real telescope data.

They also looked at two specific "extreme" neutron stars to test their recipe:

• The "Tiny" Star (HESS J1731-347): A very small, light neutron star.
• The "Giant" Star (PSR J0952-0607): The heaviest neutron star we know of.

The Verdict: The Sweet Spot

The results were exciting. They found that if the Sexaquark has a mass of about 1,900 MeV (roughly twice the mass of a proton), the model works perfectly.

• The "Goldilocks" Zone: The paper concludes that Sexaquarks with a mass between 1,885 and 1,935 MeV are the "just right" candidates.
• Why it matters: Without these Sexaquarks, the models for normal neutron stars (around 1.4 times the mass of our Sun) were too big. With the Sexaquarks, the models shrink down to match exactly what the NICER telescope sees.

The Big Picture

This paper suggests that the "ghostly" dark matter we are looking for might not just be floating around in the space between stars. It might be hiding inside the densest objects in the universe, acting as a secret ingredient that changes how these stars behave.

In summary:
The universe is like a giant kitchen. For years, our recipes for neutron stars didn't taste right—they were too stiff. This paper suggests that if we add a pinch of "Sexaquark spice" (a specific type of dark matter particle), the recipe finally works, explaining why real neutron stars are the size and shape we observe. It's a strong clue that dark matter might be made of these six-quark knots, hiding in plain sight inside the hearts of dead stars.

Technical Summary

1. Problem Statement

Neutron stars (NSs) serve as extreme laboratories for testing the Equation of State (EOS) of dense matter. Current theoretical models face significant challenges in reconciling several observational constraints:

• The "Hyperon Puzzle": The inclusion of hyperons (strange baryons) in the EOS typically softens the pressure, often preventing the formation of neutron stars with masses $\ge 2 M_\odot$, which contradicts observations of heavy pulsars like PSR J0740+6620.
• Tidal Deformability: Gravitational wave (GW) observations from the GW170817 merger constrain the tidal deformability ($\Lambda_{1.4}$) of $1.4 M_\odot$ stars, favoring softer EOSs with smaller radii.
• Radius Constraints: Recent NICER X-ray observations of pulsars (e.g., PSR J0437–4715, PSR J0614–3329) suggest radii that are smaller than predicted by many standard hadronic models, particularly for canonical $1.4 M_\odot$ stars.
• Dark Matter (DM) Candidates: There is growing interest in whether exotic particles, specifically the sexaquark (S)—a hypothetical bosonic hexaquark ($uuddss$)—could exist within NS cores. If present, S could act as a DM component, potentially altering the EOS to satisfy these conflicting constraints.

The central question is whether a hybrid model incorporating hyperons, bosonic sexaquark dark matter, and deconfined quark matter can simultaneously satisfy all current multi-messenger observational constraints, and what mass range for the sexaquark is permitted.

2. Methodology

The authors constructed a unified, multi-phase EOS and performed a rigorous Bayesian analysis against observational data.

A. Theoretical Framework

• Hadronic Phase: Utilized the DD2Y-T model, a relativistic density functional that includes all $J=1/2$ octet baryons (nucleons and hyperons: $\Lambda, \Sigma, \Xi$).
• Dark Matter Interaction: The sexaquark (S) is treated as a bosonic degree of freedom. Instead of introducing complex meson couplings, the interaction with the baryonic medium is modeled via a density-dependent effective mass shift:
  $$m^*_S = m_S + \Delta m_S = m_S + m_S x_S \frac{n_B}{n_0}$$
  where $x_S$ is a slope parameter fixed at a minimal value ($0.03$) to represent weak coupling, and $n_B$ is the baryon density. This leads to Bose-Einstein condensation of S at high densities.
• Quark Matter Phase: Modeled using a nonlocal Nambu–Jona-Lasinio (nlNJL) model, which accounts for color superconductivity. Two parameter sets were used, derived from previous Bayesian analyses, to ensure consistency with high-mass constraints.
• Phase Transition: The transition from hadronic matter to deconfined quark matter is modeled as a smooth crossover using the Replacement Interpolation Construction (RIC) method. This avoids the abrupt first-order phase transition (Maxwell construction) and allows for a continuous change in composition, consistent with lattice QCD expectations at high temperatures.

B. Observational Constraints

The study incorporated a comprehensive set of multi-messenger data:

1. Mass Constraints:
   • PSR J0348+0432 ($2.01 M_\odot$) and PSR J0740+6620 ($2.08 M_\odot$) for the maximum mass limit.
   • PSR J0952–0607 ($2.35 M_\odot$), the most massive known pulsar (Black Widow).
2. Radius & Mass Constraints (NICER):
   • PSR J0030+0451, PSR J0437–4715, and the newly published PSR J0614–3329 ($1.44 M_\odot$, $R \approx 10.3$ km).
   • HESS J1731–347, an ultracompact, low-mass object ($M \approx 0.77 M_\odot$).
3. Tidal Deformability: The constraint from the GW170817 event ($\Lambda_{1.4} = 190^{+390}_{-120}$).

C. Statistical Analysis

A Bayesian Analysis (BA) was performed to determine the posterior probability distribution of the sexaquark mass ($m_S$). The likelihood function combined independent constraints from mass, radius, and tidal deformability data.

3. Key Contributions

• Unified Hybrid Model: The paper presents a self-consistent EOS that integrates hyperons, bosonic sexaquark DM, and deconfined quark matter within a single fluid framework.
• RIC Crossover Implementation: It demonstrates the application of the RIC method to model a smooth hadron-quark crossover in the presence of a bosonic DM component, avoiding unphysical discontinuities.
• Systematic Mass Scan: The authors scanned the sexaquark mass range from 1890 MeV to 2054 MeV (the theoretical upper bound for stability) to identify the viable parameter space.
• Inclusion of Latest Data: The study uniquely incorporates the very recent NICER data for PSR J0614–3329 and the extreme constraints from HESS J1731–347 and PSR J0952–0607.

4. Results

• EOS Softening: The inclusion of the sexaquark softens the EOS in the hadronic phase due to the onset of Bose-Einstein condensation. This softening is crucial for reducing the radii of canonical $1.4 M_\odot$ stars, bringing them into agreement with NICER measurements and the GW170817 tidal deformability limit.
• Mass-Radius Relations:
  • Models without DM (pure DD2Y-T or hybrid without S) fail to satisfy the radius constraints for PSR J0614–3329 and the tidal deformability limit.
  • Hybrid models with $m_S \lesssim 1935$ MeV successfully pass through the 68% credible regions of PSR J0437–4715 and PSR J0614–3329.
  • The presence of S delays the onset of deconfinement to higher densities, allowing the star to support high masses ($>2 M_\odot$) while maintaining small radii at lower masses.
• Tidal Deformability: The hybrid EOSs with sexaquarks satisfy the $\Lambda_{1.4}$ constraint for masses up to $\approx 2040$ MeV. Pure hadronic models and those with heavier S particles exceed the upper limit.
• Bayesian Posterior Distribution:
  • The analysis yields a strong preference for lower sexaquark masses.
  • The 68.3% credible region is $1885 \text{ MeV} \lesssim m_S \lesssim 1935 \text{ MeV}$.
  • Masses above 2000 MeV are strongly disfavored (falling outside the 95.5% credibility region).
  • The inclusion of the ultracompact HESS J1731–347 further tightens the preference toward the lower end of the mass range ($m_S \lesssim 1900$ MeV).

5. Significance

• Resolution of Tensions: The study demonstrates that the presence of a bosonic sexaquark DM component provides a viable mechanism to resolve the tension between the requirement for stiff EOSs (to support $2 M_\odot$ stars) and soft EOSs (to satisfy tidal deformability and radius constraints).
• New DM Constraints: It places the first astrophysical constraints on the sexaquark mass specifically derived from a hybrid star model, narrowing the viable mass window to the $\sim 1.9$ GeV scale. This distinguishes S from other DM candidates like axions (µeV-MeV) or WIMPs (GeV-TeV).
• Implications for NS Structure: The results suggest that if sexaquarks exist, they likely form a significant fraction (12–15%) of the core density in neutron stars, coexisting with hyperons and quark matter.
• Methodological Advance: The work highlights the necessity of using smooth crossover transitions and incorporating the latest multi-messenger data (especially from NICER and GW events) to accurately probe the internal composition of neutron stars.

In conclusion, the paper argues that a hybrid neutron star model containing hyperons, deconfined quark matter, and a bosonic sexaquark dark matter component with a mass in the range of 1885–1935 MeV is the most consistent with all current observational data, offering a compelling solution to the equation of state puzzle.