From $χ$EFT to Multi-Region Modeling: Neutron star structure with a polytropic extension of $χ$EFT and MUSES Calculation Engine multi-layer modeling

This paper presents a comparative study of neutron star structure by contrasting a Chiral Effective Field Theory ($\chi$EFT) approach with the MUSES Calculation Engine's multi-region polytropic modeling, analyzing their respective mass-radius relations, advantages, and limitations.

The Gist

The Cosmic Puzzle: What's Inside a Neutron Star?

Imagine a neutron star. It's the collapsed core of a massive star that exploded. It's so dense that if you took a teaspoon of its material, it would weigh about a billion tons on Earth. It's essentially a giant atomic nucleus the size of a city.

The big mystery physicists are trying to solve is: What happens to matter when you squeeze it that hard?

We know how atoms behave in a lab, but inside a neutron star, the pressure is so extreme that the rules change. The "Equation of State" (EoS) is just a fancy name for a rulebook that tells us how this super-dense matter reacts to pressure. If we have the right rulebook, we can predict how big a neutron star can get before it collapses into a black hole.

The Two Approaches: The "Strict Architect" vs. The "Modular Builder"

This paper compares two different ways of building this rulebook. The author, Federico Nola, is essentially testing two different construction crews to see who builds the most realistic neutron star.

1. The "Strict Architect" (The $\chi$EFT + Polytrope Method)

The Analogy: Imagine you are building a tower. You start with a very precise, scientific blueprint for the bottom 10 floors (the lower density layers). This blueprint is based on Chiral Effective Field Theory ($\chi$EFT), which is a highly respected, math-heavy theory derived from the fundamental laws of particle physics.

However, the architect knows that the blueprint stops working once you get too high up because the physics gets too weird to calculate perfectly. So, for the top part of the tower, they switch to a "rule of thumb." They say, "Okay, for the next 5 floors, let's just assume the material gets stiffer in a smooth, predictable curve (a polytrope)." Then, for the very top, they add a safety cap to make sure the tower doesn't violate the laws of physics (like moving faster than light).

• Pros: It's mathematically clean and stays very close to the known laws of physics at the bottom.
• Cons: The top part is a bit of a guess. It's a smooth curve, but it might not capture the complex, chaotic changes happening deep inside the star.

The Result: This method built a neutron star that is smaller and heavier (about 2.17 times the mass of our Sun, but only 10 km wide). It's a very compact, dense object.

2. The "Modular Builder" (The MUSES Calculation Engine)

The Analogy: Imagine a different construction crew that uses a high-tech, modular toolkit called MUSES. Instead of using one rule for the whole tower, they use three different, specialized teams for three different zones:

• Zone 1 (The Crust): They use a model designed for the "crusty" outer layer.
• Zone 2 (The Core): They use the same strict scientific blueprint ($\chi$EFT) as the first crew for the middle section.
• Zone 3 (The Deep Core): Instead of a simple "rule of thumb," they switch to a completely different, advanced physics engine called Chiral Mean Field (CMF). This engine is designed to handle the extreme chaos of the deepest core, simulating how particles might behave in a soup of quarks and gluons.

They stitch these three different models together seamlessly.

• Pros: It's more realistic because it acknowledges that the physics changes drastically as you go deeper. It doesn't just guess; it uses a different, sophisticated theory for the deep core.
• Cons: It's much more complex to run and requires more computing power.

The Result: This method built a neutron star that is slightly lighter and wider (about 2.04 solar masses, but 11.3 km wide). It's a bit "fluffier" than the first one.

The Showdown: Who Won?

The author compared the two results against real-world data from telescopes and gravitational wave detectors (like LIGO).

• Both methods work: Both teams built stars that are consistent with what we see in the universe.
• The difference: The "Strict Architect" made a tighter, denser star. The "Modular Builder" made a slightly larger, less dense star.

The paper concludes that neither method is perfect, and neither is strictly "better." They are complementary tools:

• The Polytrope method is great for testing "what if" scenarios. It helps us see how much the assumptions about the deep core change the size of the star.
• The MUSES method is better for deep physical realism. If we want to know if there are exotic particles (like axions or quarks) hiding in the core, MUSES is the better tool because it can swap out the physics models easily.

The Big Picture

Think of this like trying to predict the weather.

• Method 1 says, "We know the weather perfectly for the last 10 miles. After that, let's just assume it gets a little hotter and stick to a simple formula."
• Method 2 says, "We know the weather for the last 10 miles. After that, we'll switch to a super-computer model that simulates ocean currents, jet streams, and humidity separately."

Both give you a forecast, but Method 2 gives you a more detailed picture of why the weather is changing, while Method 1 is faster and easier to tweak to see how sensitive the forecast is to your assumptions.

The Takeaway: To truly understand the universe's most extreme objects, we need both approaches. We need the simple, controlled models to test our limits, and the complex, modular models to explore the unknown depths of the neutron star core.

Technical Summary

1. Problem Statement

Neutron stars (NS) represent extreme environments where matter reaches densities comparable to atomic nuclei ($\bar{\rho} \sim 10^{15}$ g cm$^{-3}$). The internal structure of these stars is governed by the Equation of State (EoS) of cold, dense matter, which remains unknown at supranuclear densities.

• The Challenge: While Chiral Effective Field Theory ($\chi$EFT) provides a rigorous, model-independent description of nuclear interactions at low densities, it breaks down at high densities (typically $n \gtrsim 0.3\text{--}0.4$ fm$^{-3}$) due to the convergence of the chiral expansion and the potential emergence of non-nucleonic degrees of freedom (e.g., hyperons, quarks).
• The Gap: There is a need to compare different strategies for extending the EoS from the reliable $\chi$EFT regime to the high-density core to understand how high-density modeling uncertainties affect macroscopic observables like the Mass-Radius ($M-R$) relation.

2. Methodology

The authors construct and compare two distinct approaches to modeling the neutron star EoS, both starting from the same microscopic $\chi$EFT input for low densities.

Approach A: $\beta$-stabilized $\chi$EFT with Polytropic Extension ($P_{\chi\text{EFT-npe}\mu}$)

This method uses a purely nuclear approach extended phenomenologically:

1. Low Density ($n \le n_1$): Uses a $\beta$-stable EoS derived from $\chi$EFT (up to N$^3$LO for two-body and N$^2$LO for three-body forces) combined with Goldstone perturbation theory. It includes neutrons, protons, electrons, and muons ($npe\mu$).
2. Intermediate Density ($n_1 \le n \le n_2$): Matches the $\chi$EFT results to a single stiff polytropic segment ($P \propto n^{\Gamma_1}$).
3. High Density ($n \ge n_2$): Extends the EoS using a Speed-of-Sound (SoS) parametrization constrained by causality ($v_s^2 \le 1$). The SoS profile transitions smoothly from the polytrope to a constant high-density plateau.
4. Crust: A standard crust EoS (Negele-Vautherin/Baym-Pethick-Sutherland) is attached at low pressures.

Approach B: MUSES Calculation Engine Multi-Layer Modeling ($MUSES\text{-}npe\mu$)

This method utilizes the MUSES Calculation Engine, a modular framework that constructs the EoS by stitching together different physical models valid in specific density regimes:

1. Outer Crust: Uses Crust Density Functional Theory (Crust-DFT) based on extended Skyrme functionals constrained by $\chi$EFT and nuclear masses.
2. Inner Crust: Employs a Compressible Liquid Drop Model (CLDM) to handle the "pasta" phases and the transition to uniform matter.
3. Core (Saturation to High Density):
   • Near saturation: Uses $\chi$EFT.
   • Beyond saturation: Uses the Chiral Mean Field (CMF) model (specifically the CMF++ implementation). This is a relativistic, non-linear $\sigma$-model based framework that allows for the study of phase transitions and high-density matter while maintaining causality.
4. Synthesis: The modules are unified via thermodynamic matching (using hyperbolic tangent functions) to ensure continuity in energy density and pressure. The lepton module ensures charge neutrality and $\beta$-equilibrium throughout.

Stellar Structure Calculation

Both EoSs are used to solve the Tolman-Oppenheimer-Volkoff (TOV) equations to generate Mass-Radius ($M-R$) sequences. The study restricts the core to nucleonic degrees of freedom ($n, p$) to isolate the impact of the high-density modeling strategy rather than composition changes (like hyperons).

3. Key Contributions

• Comparative Framework: The paper provides a direct comparison between a parametric extrapolation (polytrope + SoS) and a microphysics-based multi-layer framework (MUSES) starting from identical low-density $\chi$EFT inputs.
• MUSES Implementation: It demonstrates the application of the MUSES Calculation Engine for constructing a complete, high-precision neutron star model that explicitly accounts for the different physical nature of the crust, inner crust, and core.
• Parameter Sensitivity Analysis: The authors perform a sensitivity analysis (Appendix A) on the high-density extension parameters (adiabatic index $\Gamma_1$, matching densities, and SoS cap), showing how they influence the stiffness of the EoS and the resulting $M-R$ curves.

4. Results

The study yields the following quantitative and qualitative results:

• Mass-Radius Relations:
  • $P_{\chi\text{EFT-npe}\mu}$ (Polytropic): Predicts a more compact star with a maximum mass $M_{\text{max}} \approx 2.17 M_\odot$ and a radius at maximum mass $R \approx 10.13$ km.
  • $MUSES\text{-}npe\mu$ (Multi-layer): Predicts a slightly less massive but larger star with $M_{\text{max}} \approx 2.04 M_\odot$ and $R \approx 11.3$ km.
• Consistency: Both models are consistent with observational constraints (e.g., from NICER and LIGO/Virgo) and the 50%/90% confidence regions derived from Gaussian Process analyses of pulsar data.
• Divergence at High Density: The models agree well at low densities (where $\chi$EFT is reliable) but diverge significantly at high densities. The polytropic extension produces a stiffer EoS at the highest densities compared to the CMF-based MUSES approach.
• Uncertainty Source: The results highlight that for a fixed composition (nucleonic only), the dominant source of uncertainty in global stellar properties arises from the treatment of the high-density regime.

5. Significance and Conclusion

• Complementarity: The paper concludes that neither approach is universally superior; they are complementary tools.
  • The Polytropic Extension is an economical, flexible tool for bracketing the impact of high-density stiffness and testing causality constraints.
  • The MUSES Engine offers a more physically grounded route to the high-density regime, embedding the core in a framework that can naturally incorporate phase transitions and non-nucleonic degrees of freedom (e.g., quarks) in future studies.
• Future Applications: The MUSES framework is highlighted as a versatile platform for exploring hybrid stars and testing particle physics theories (e.g., axion-nucleon coupling limits from cooling) where the EoS is a critical input.
• Scientific Impact: This work clarifies the respective roles of purely parametric extrapolations versus microphysics-based frameworks, providing a consistent baseline for future multi-messenger astrophysics studies.