Many-body effects on dense matter with hyperons at finite temperature

This paper presents the first extension of the Many-Body Forces (MBF) Model to finite temperature, introducing new hyperon coupling schemes to analyze the thermodynamic properties of beta-equilibrated nuclear matter and the mass-radius relations of compact stars, thereby establishing a new framework for describing proto-neutron stars.

The Gist

Imagine the universe is filled with a cosmic "soup" of matter so dense that a single teaspoon would weigh as much as a mountain. This is the stuff inside neutron stars, the collapsed cores of dead, massive stars. For a long time, scientists have tried to figure out exactly how this soup behaves, but it's incredibly difficult to study because we can't recreate such extreme conditions in a lab.

This paper is like a new, upgraded recipe book for that cosmic soup. Specifically, the authors are updating a theoretical model called the Many-Body Forces (MBF) Model to include two things that were previously missing or handled roughly: heat and strange particles.

Here is a breakdown of what they did, using simple analogies:

1. The Problem: The "Hard" Math of the Universe

To understand how matter behaves at these extreme densities, physicists usually rely on a fundamental theory called Quantum Chromodynamics (QCD). However, using QCD to describe a neutron star is like trying to solve a puzzle where every piece is constantly changing shape and talking to every other piece at the same time. It's mathematically impossible to solve directly.

So, scientists use "effective theories." Think of these as simplified maps. Instead of drawing every single tree and rock (quarks and gluons), the map just shows the roads and cities (protons, neutrons, and other particles). The authors use a specific map called the MBF Model.

2. The Upgrade: Adding Heat and "Strange" Guests

The authors took their existing map and added two major features:

• Finite Temperature (Heat): Most previous models assumed the star was "cold" (frozen in time). But when a star is born (a "proto-neutron star"), it is incredibly hot—like a furnace. The authors updated their model to simulate this heat.
  • Analogy: Imagine a crowded dance floor. In a "cold" model, everyone is standing still in a rigid formation. In this new "hot" model, everyone is dancing wildly, bumping into each other, and moving around. This changes how the crowd pushes against the walls (pressure).
• Hyperons (The Strange Guests): In normal matter, you have protons and neutrons. But in the deep, dense core of a star, it becomes energetically favorable to create heavier, "strange" particles called hyperons.
  • Analogy: Imagine a party where the room gets so crowded that the host decides to let in some larger, heavier guests (hyperons). These new guests take up space and change the dynamic of the room. The paper explores how different "rules" for how these guests interact with the original party-goers change the outcome.

3. The Experiment: Testing Different "Rules"

The authors didn't just run one simulation; they tested different scenarios to see which one makes the most sense:

• The "Stiffness" Knob: They adjusted a parameter (called $\zeta$) that controls how "stiff" or "squishy" the matter is.
  • Stiff Matter: Like a solid steel block. It resists being squeezed.
  • Soft Matter: Like a sponge. It squishes easily.
  • They tested a "stiff" setting and a "soft" setting to see how the star reacts.
• The Interaction Schemes: They tried three different ways the "strange guests" (hyperons) interact with the "regular guests" (protons/neutrons).
  • Universal: Everyone interacts the same way.
  • Moszkowski: A specific rule based on particle composition.
  • SU(6): A complex rule based on symmetry and flavor.

4. The Results: What Happens to the Star?

By running these simulations, they calculated how the star's pressure, speed of sound, and size change.

• The "Hyperon Puzzle": A big mystery in physics is that hyperons usually make matter "soft" (squishy). If matter is too squishy, the star collapses under its own gravity, and the model predicts a maximum mass that is too small (less than 2 times the mass of our Sun). But we know neutron stars exist that are heavier than that.
• The Solution: The authors found that if they use the "stiff" setting ($\zeta = 0.040$) in their model, the matter stays strong enough to support heavy stars, even with the strange guests present.
• The "Soft" Setting Failure: If they used the "soft" setting ($\zeta = 0.129$), the star collapsed too easily, and the model failed to match the heavy stars we actually observe in the sky.
• Heat Helps: Interestingly, the heat in the early stages of a star's life (the proto-neutron star phase) acts like a temporary support beam. It keeps the star slightly larger and prevents it from collapsing as quickly as a cold star would.

5. The Conclusion: A Better Map for the Cosmos

The paper concludes that their updated model is a powerful tool. It successfully describes how dense matter behaves when it is both hot and filled with strange particles.

• The "Stiff" version of their model matches real-world observations of heavy neutron stars perfectly.
• The "Soft" version does not.

Essentially, they have provided a more accurate "recipe" for the densest matter in the universe. This helps astronomers understand how neutron stars are born, how they evolve as they cool down, and why some of them are massive enough to survive without collapsing into black holes.

In short: They updated the math to include heat and strange particles, tested different interaction rules, and found that a specific "stiff" version of their model is the only one that explains the heavy neutron stars we see in the universe today.

Technical Summary

Technical Summary: Many-body effects on dense matter with hyperons at finite temperature

Problem Statement
The equation of state (EoS) of dense matter at high densities and finite temperatures remains a significant challenge in nuclear physics and astrophysics. While Quantum Chromodynamics (QCD) is the fundamental theory of strong interactions, it is currently intractable for direct calculation in the high-density, finite-temperature regime relevant to neutron star (NS) interiors due to the "sign problem" in lattice QCD. Consequently, effective field theories, such as Relativistic Mean-Field (RMF) models, are employed. However, standard RMF models often struggle to simultaneously reproduce nuclear matter saturation properties, the existence of massive neutron stars ($>2 M_\odot$), and the presence of hyperons (the "hyperon puzzle"). Furthermore, previous applications of the Many-Body Forces (MBF) Model, which accounts for many-body interactions through field-dependent coupling constants, were restricted to zero temperature. This work addresses the need for a consistent description of dense matter that includes hyperons, many-body forces, and finite temperature effects, which are critical for modeling proto-neutron stars (PNS) and supernova dynamics.

Methodology
The authors extend the Many-Body Forces (MBF) Model to finite temperature within a relativistic quantum hadrodynamics formalism. The core of the model involves a parameterizable derivative coupling where meson-baryon coupling constants depend on the scalar meson fields, thereby introducing an implicit density dependence that mimics many-body correlations.

• Formalism: The Lagrangian density includes the baryon octet (nucleons and hyperons $\Lambda, \Sigma, \Xi$), leptons ($e, \mu$), and meson fields ($\sigma, \omega, \rho, \delta, \sigma^*, \phi$). The model utilizes the mean-field approximation to solve the equations of motion.
• Thermodynamics: The system is treated as charge-neutral, beta-equilibrated matter. The authors calculate thermodynamic quantities (energy density, pressure, entropy) using Fermi-Dirac distributions for both particles and antiparticles, ensuring the inclusion of thermal effects.
• Scenarios: Three specific evolutionary snapshots of a compact star are analyzed based on fixed entropy per baryon ($S/A$):
  1. $S/A = 1$ with a fixed lepton fraction ($Y_l = 0.4$), representing a hot PNS with trapped neutrinos.
  2. $S/A = 2$ with zero neutrino chemical potential, representing a later stage with deleptonization and "Joule heating."
  3. $S/A = 0$ (or $T=0$), representing a cold, fully evolved neutron star.
• Parameterization: The study focuses on the scalar (S) version of the MBF Model, controlled by a single adjustable parameter $\zeta$. Two extreme values are compared: $\zeta = 0.040$ (stiffer EoS) and $\zeta = 0.129$ (softer EoS).
• Hyperon Couplings: Three distinct coupling schemes for hyperons are investigated: Universal (U), Moszkowski (M), and Spin-Flavor Symmetry (SU(6)). Additionally, the role of the strange scalar meson $\sigma^*$ is tested by varying its coupling $g_{\sigma^*Y}$ between 0 and 5.
• Astrophysical Constraints: The resulting EoSs are used to solve the Tolman-Oppenheimer-Volkoff (TOV) equations to determine mass-radius relations. These are compared against multimessenger observational data, including NICER X-ray timing and LIGO/Virgo gravitational wave measurements.

Key Contributions

1. First Finite-Temperature Extension of MBF: This work presents the first application of the MBF Model to finite-temperature dense matter, providing a framework to study the thermodynamic evolution of proto-neutron stars.
2. New Hyperon Coupling Schemes: The authors introduce and implement three different hyperon coupling schemes (U, M, SU(6)) within the MBF formalism for the first time, allowing for a systematic study of how many-body forces interact with different hyperon interaction assumptions.
3. Unified Thermodynamic Analysis: The paper provides a comprehensive calculation of key physical quantities—speed of sound, compressibility, adiabatic index, and particle populations—across the full range of densities and temperatures relevant to PNS evolution.
4. Resolution of the Hyperon Puzzle (Conditional): The study demonstrates that the MBF Model, specifically with the stiffer parameterization ($\zeta = 0.040$), can support neutron stars with masses exceeding $2 M_\odot$ even in the presence of hyperons, effectively addressing the hyperon puzzle without requiring exotic repulsive mechanisms beyond the model's standard many-body terms.

Results

• Equation of State (EoS): The inclusion of hyperons generally softens the EoS. However, the degree of softening is highly dependent on the coupling scheme and the parameter $\zeta$. The Universal coupling (U) yields the stiffest hyperonic EoSs, while SU(6) yields the softest. The parameter $\zeta = 0.040$ consistently produces stiffer EoSs than $\zeta = 0.129$.
• Particle Populations: At finite temperature and entropy, hyperons appear at lower densities compared to the $T=0$ case due to thermal smearing. Neutrino trapping ($S/A=1$) delays the onset of negatively charged hyperons by increasing the electron chemical potential. The SU(6) scheme predicts the earliest and most abundant hyperon population, while the Universal scheme suppresses them the most.
• Thermodynamic Properties:
  • Speed of Sound ($c_s^2$): The model strictly respects causality ($c_s^2 \leq 1$) in all cases. The speed of sound exhibits non-monotonic behavior (peaks and drops) corresponding to the onset of new hyperon species, particularly at $T=0$. These features are smoothed out at finite temperatures.
  • Compressibility ($K$) and Adiabatic Index ($\Gamma$): Both quantities show sensitivity to hyperon thresholds. The adiabatic index remains above the stability limit ($\Gamma > 4/3$) for all configurations. Thermal effects generally increase compressibility, though this is modulated by the hyperon fraction.
• Mass-Radius Relations:
  • Cold Stars ($T=0$): The $\zeta = 0.040$ parameterization supports maximum masses $> 2 M_\odot$ for all coupling schemes, consistent with observations of heavy pulsars (e.g., PSR J0740+6620). In contrast, the $\zeta = 0.129$ parameterization fails to support masses above $\sim 1.8 M_\odot$ when hyperons are included, conflicting with current data.
  • Proto-Neutron Stars: Finite temperature and neutrino trapping lead to larger stellar radii compared to cold stars. The "thermal puffing" effect is significant, with radii increasing by $\sim 25\%$ in hot scenarios. Neutrino trapping delays hyperonization, temporarily stiffening the EoS and allowing hot stars to sustain masses comparable to or slightly higher than their cold counterparts.
• Strange Scalar Meson ($\sigma^*$): The inclusion of the $\sigma^*$ meson has a subtle effect on global properties, slightly enhancing attraction and reducing maximum masses, but does not qualitatively alter the ordering of results or the model's ability to satisfy observational constraints when $\zeta = 0.040$.

Significance
The paper claims that the MBF Model, extended to finite temperature, offers a robust and versatile tool for describing dense baryonic matter. Its primary significance lies in demonstrating that many-body forces, parameterized via field-dependent couplings, can naturally reconcile the presence of hyperons with the existence of massive neutron stars ($>2 M_\odot$) without invoking ad hoc repulsive terms. The work highlights that the choice of the adjustable parameter $\zeta$ and the hyperon coupling scheme are critical determinants of the EoS stiffness. The results suggest that the "hyperon puzzle" can be resolved within standard effective field theories if the many-body interactions are sufficiently strong (as in the $\zeta = 0.040$ case). Furthermore, the study provides a self-consistent thermal treatment essential for modeling the early evolution of proto-neutron stars, showing that thermal effects and neutrino trapping significantly alter the internal composition and macroscopic structure of these objects. The authors conclude that the MBF Model is capable of reproducing a wide range of EoSs consistent with current multimessenger astrophysical constraints, provided the model parameters are appropriately selected.