Hyperon Coupling Freedom Narrows the Neutron Star Radius Shift Seen in SU(6) Models

What earlier work established about hyperons in neutron stars

The baseline picture is well established: as density rises in a neutron star core, hyperons become energetically favorable once their Fermi energy exceeds their rest mass, softening the equation of state and reducing the maximum mass [2][5]. This creates the hyperon puzzle—the tension between the theoretical expectation of hyperons and the observational requirement that neutron stars reach at least ~2 solar masses [2][5]. Systematic comparisons within relativistic mean-field models had suggested that to compensate for hyperon softening, the underlying nucleonic equation of state must be stiffer, which in turn pushes radii upward for intermediate-mass stars [1]. The SU(6) symmetry relations, which fix hyperon-vector-meson couplings to nucleonic ones via valence quark counting, became the standard tool for this [3][4].

The SU(6) relations are attractive because they remove arbitrariness: they uniquely determine the vector couplings from the nucleonic ones [3]. But they are also a strong assumption. The strange and non-strange sectors need not respect naive quark-model symmetry at the baryonic level, and the scalar coupling still has to be fixed separately from hypernuclear data [4]. The new paper [1] directly tests what happens when that assumption is relaxed.

How the new metamodel separates coupling freedom from nucleonic uncertainty

The anchor study [1] uses a generalized relativistic density functional (GDFM) as a generator of a large ensemble of equations of state, rather than calibrating a single interaction. It performs a Bayesian analysis with constraints from chiral effective field theory for pure neutron matter, the AME2016 nuclear mass table, a Λ optical potential of −30 ± 5 MeV, the maximum TOV mass from pulsar measurements, and tidal deformability from GW170817 [1]. Two hyperonic coupling schemes are compared: the standard SU(6) relations and a 'Ratio' scheme in which the hyperonic couplings are proportional to the nucleonic ones but the ratios RσΛ, RωΛ, RφΛ, and Rσ*Λ are varied within ranges [1].

Crucially, the Ratio scheme does not correspond to fully free variation—the density dependence of the hyperonic couplings is still tied to the nucleonic ones—but it does decouple the strange and non-strange sectors at the level of coupling strength [1]. This is the key methodological move: it isolates the effect of SU(6) restriction from the effect of nucleonic equation-of-state uncertainty.

The radius shift collapses when SU(6) is relaxed

The central result is stark. In the SU(6) case, the most probable equation of state has a stiffer nucleonic part and higher proton fraction at low densities, and the predicted radii for intermediate-mass neutron stars are relatively high—qualitatively agreeing with previous studies [1]. But in the Ratio case, no apparent softening is observed in the representative equation of state despite similar hyperon fractions and TOV maximum masses [1]. The mass-radius contours for the two settings overlap almost completely for stars below ~1.4 solar masses, and the differences that appear at higher masses are modest [1].

The interpretation is that the radius shift attributed to hyperons in SU(6) models is not a generic consequence of hyperon appearance. It arises because the SU(6) restriction forces the nucleonic equation of state to be stiffer to satisfy the maximum mass constraint, and that stiffness—not the hyperons themselves—drives the larger radii [1]. When coupling freedom allows the hyperonic sector to adjust, the nucleonic part need not be as stiff, and the radius signature fades.

How this compares with other Bayesian hyperon coupling analyses

A competing Bayesian study by Sun et al. [4] used a set of stiff relativistic mean-field equations of state and inferred hyperon couplings from GW170817 tidal deformability and NICER mass-radius measurements of PSR J0030+0451 and PSR J0740+6620. They found that the laboratory hypernuclear constraint on the RσΛ–RωΛ correlation ensures a large enough scalar coupling to match the large vector coupling, and that the maximum mass of hyperon stars is at most 2.176 +0.085/−0.202 solar masses at 68% credibility [4]. Their analysis also relaxed SU(6) for the Λ sector but kept SU(6) for Σ and Ξ couplings [4].

The key difference is scope: Sun et al. [4] work within a fixed set of stiff nucleonic equations of state and focus on constraining the Λ coupling ratios, while the anchor paper [1] varies both nucleonic and hyperonic sectors simultaneously within a metamodel. This means the anchor paper can ask whether the radius shift survives when nucleonic stiffness is not pre-selected—and it finds that it does not [1]. The two studies agree that hyperons soften the equation of state and reduce the maximum mass, but they differ in how much of the radius signature is attributable to hyperons versus to the underlying nucleonic model.

Nucleonic functional choice matters more for composition than for radii

The anchor paper also compares two nucleonic functional forms, GDFM and TW, albeit only in the SU(6) case [1]. It finds that the equation-of-state distribution and global stellar properties, particularly the mass-radius relation and the radius of a 1.4 solar mass star, differ only modestly between the two functionals. The differences that do appear are already present at the nucleonic level rather than being introduced by hyperons [1]. Composition, by contrast—specifically the proton fraction—is markedly more sensitive to the functional form [1].

This pattern matches what was established for purely nucleonic models: different functionals can yield similar mass-radius sequences while differing substantially in their underlying proton fractions [1]. The implication is that the functional-form dependence in hyperonic results mainly affects the weakly constrained composition sector, leaving the equation of state and global observables comparatively robust. The paper also checks thermodynamic stability via a convexity analysis of the energy density and finds that instabilities only appear for Rσ*Λ ≳ 1, corresponding to an attractive ΛΛ coupling that overcomes the nucleonic one—a regime that does not seem realistic [1].

Where the conclusion stops and what remains open

The conclusion that the radius shift is strongly reduced by coupling freedom is limited to the metamodel framework and the density functionals used in the anchor paper [1]. The study does not perform a full Bayesian constraint against specific astronomical observations; it imposes the maximum mass and tidal deformability as constraints but does not derive posterior distributions for the hyperon couplings themselves [1]. The Ratio scheme is not fully free—the density dependence of the hyperonic couplings is still tied to the nucleonic ones—so a fully free exploration of the hyperonic parameter space remains beyond the scope of this work [1].

The hyperon puzzle itself is not resolved. The maximum mass still drops when hyperons appear, and the question of how to reconcile that drop with the observed ~2 solar mass pulsars remains open [1][2]. What the new work changes is the interpretation of the radius signature: it is not a generic prediction of hyperonic matter but a consequence of specific coupling assumptions. Future work with fully free hyperonic couplings and direct Bayesian comparison to NICER and gravitational-wave data will be needed to determine whether the SU(6) picture has any remaining observational support.

About These Sources

This research page is built on 5 studies (4 peer-reviewed, 1 preprint) — published from 2021 to 2025, 1 from 2024 or later, 3 in Q1 journals, collectively cited 323 times — selected as the most relevant from 13 studies that passed quality screening, drawn from 130 papers retrieved from a database of over 500 million.

Sources used in this answer