---
title: Bulk Viscosity in Hyperonic Matter
url: https://www.emergentmind.com/topics/bulk-viscosity-in-hyperonic-matter
type: topic
---

# Bulk Viscosity in Hyperonic Matter

Searching arXiv for recent and foundational papers on hyperonic bulk viscosity, merger matter, and r-mode damping.
Bulk viscosity in hyperonic matter is the dissipative response generated when compression and rarefaction drive strange baryon populations out of chemical equilibrium and weak strangeness-changing reactions restore equilibrium with a finite lag. In neutron-star applications, the relevant medium is typically a baryon-lepton mixture containing nucleons together with hyperons such as \(\Lambda\), \(\Sigma^{-}\), and, in some models, \(\Xi^{-}\). The mechanism is important because it converts mechanical oscillation energy into heat, thereby damping stellar pulsations, \(r\)-modes, and related hydrodynamic motions. The modern literature treats two rather different regimes: hot merger matter, where full finite-temperature phase-space integrals and relativistic kinematics are required, and colder hyperonic cores of rotating neutron stars, where the same nonequilibrium chemistry strongly affects \(r\)-mode instability windows [2009.05181], [1911.08407].

## 1. Physical setting and equilibrium structure

Hyperonic bulk viscosity is studied in dense matter subject to baryon-number conservation, charge neutrality, and weak-interaction equilibrium constraints. In the merger-oriented treatment based on the hyperonic RMF model PK1+H, the composition contains \(n,p,\Lambda,\Sigma^{-}\), electrons, muons, and the mesons \(\sigma,\omega,\rho\), with comparison to GM1'B used as a robustness check [2009.05181]. In the neutron-star-core literature, related compositions include \(n,p,\Lambda,\Sigma^{-},e,\mu\) [0806.4914] and, in modern \(\Lambda\Xi^{-}\) scenarios, \(npe\mu\Lambda\Xi^{-}\) matter [1911.08407].

For the \(n,p,\Lambda,\Sigma^{-},e,\mu\) system, the equilibrium conditions can be written as
\[
n_B = n_n + n_p + n_{\Sigma^-} + n_\Lambda,
\]
\[
n_p = n_e + n_\mu + n_{\Sigma^-},
\]
\[
\mu_p = \mu_n - \mu_e,\qquad \mu_e = \mu_\mu,
\]
\[
\mu_{\Sigma^-} = \mu_n + \mu_e,\qquad \mu_\Lambda = \mu_n.
\]
In the RMF description used for merger matter, baryons satisfy the relativistic dispersion relation
\[
E_i=\sqrt{p_i^2+\left(M_i^*\right)^2} +g_{\omega i}\langle\omega_0\rangle +g_{\rho i}I_{i3}\langle\rho_{03}\rangle,
\qquad
M_i^* = M_i - g_{\sigma i}\langle \sigma \rangle .
\]
These relations determine the susceptibilities that enter bulk-viscosity formulae and the threshold densities for hyperon appearance [2009.05181].

The astrophysical environments differ markedly. Neutron-star mergers are characterized by densities of several times nuclear saturation density, temperatures of order MeV to tens of MeV, and density oscillation frequencies around \(1\ \mathrm{kHz}\) [2009.05181]. By contrast, the classical \(r\)-mode problem concerns colder rotating neutron stars, often in the \(10^8\)–\(10^9\) K range, where hyperonic reactions can dominate over nucleonic Urca channels in the damping budget [2208.14436], [1911.08407].

## 2. Nonequilibrium reactions and the origin of dissipation

The canonical non-leptonic \(|\Delta S|=1\) reactions in \(n,p,\Lambda,\Sigma^{-}\) matter are
\[
n+n \leftrightarrow p+\Sigma^-,
\]
\[
n+p \leftrightarrow p+\Lambda,
\]
\[
n+n \leftrightarrow n+\Lambda,
\]
\[
n+\Lambda \leftrightarrow \Lambda+\Lambda.
\]
A fast strong reaction,
\[
n+\Lambda \leftrightarrow p+\Sigma^-,
\]
is usually assumed to remain in equilibrium, reducing the number of independent chemical imbalances [2009.05181], [0806.4914].

Under that constraint, the relevant weak imbalance becomes
\[
\delta\mu = 2\mu_n-\mu_p-\mu_{\Sigma^-} = \mu_n-\mu_\Lambda.
\]
In the subthermal regime,
\[
\delta\mu \ll T,
\]
the net rate difference for each channel is linearized as
\[
\Delta\Gamma_i = \lambda_i\,\delta\mu_i,
\]
and, because the strong-equilibrium condition ties the channels together, all four reactions relax the same \(\delta\mu\), so that the total relaxation strength is
\[
\lambda = \lambda_{\text{I}}+\lambda_{\text{II}}+\lambda_{\text{III}}+\lambda_{\text{IV}} .
\]
This is the microscopic origin of bulk viscosity: density oscillations perturb the composition, weak reactions restore equilibrium with a phase lag, entropy is produced, and oscillation energy is dissipated [2009.05181].

In the \(\Lambda\Xi^{-}\) composition used in modern \(r\)-mode studies, the reaction network changes. If only \(\Lambda\) is present, the included processes are
\[
n + p \leftrightarrow \Lambda + p,\qquad
n + n \leftrightarrow \Lambda + n,\qquad
n + \Lambda \leftrightarrow \Lambda + \Lambda.
\]
Once \(\Xi^{-}\) appears, additional channels are
\[
n + \Xi^- \leftrightarrow \Lambda + \Xi^-,
\qquad
\Lambda + n \leftrightarrow \Xi^- + p.
\]
Strong-reaction equilibrium again reduces the dynamics to a common imbalance,
\[
\Delta\mu = \mu_n - \mu_\Lambda ,
\]
and the net source is governed by the summed coefficient \(\lambda=\sum_\alpha \lambda_\alpha\) [1911.08407].

A recurrent misconception is that only the two traditional channels \(n+n\leftrightarrow p+\Sigma^{-}\) and \(n+p\leftrightarrow p+\Lambda\) matter. Modern treatments explicitly stress that \(n+n\leftrightarrow n+\Lambda\) and \(n+\Lambda\leftrightarrow\Lambda+\Lambda\) can be important, especially when one-meson exchange is included [2009.05181], and that the entire reaction network changes when the first hyperons are \(\Lambda\) and \(\Xi^{-}\) rather than \(\Sigma^{-}\) [1911.08407].

## 3. Linear-response theory and standard forms of the viscosity

For non-superfluid hyperonic matter, the bulk viscosity is defined through the pressure lag relative to the equilibrium pressure,
\[
P-P_{\rm eq} \equiv -\xi\,{\rm div}\,{\bf u},
\]
or, in the relaxation-time notation common in \(r\)-mode studies,
\[
p-\tilde p = -\zeta\, \nabla\cdot \vec v .
\]
The linear-response result derived for hot merger matter is
\[
\Re \zeta = \frac{\lambda\, \frac{\partial P}{\partial x_H}\, \frac{\partial \delta\mu}{\partial n_B}}
{\omega^2+\lambda^2\left(\frac{\partial \delta\mu}{\partial x_H}\right)^2/n_B^2},
\]
which can be rewritten as
\[
\Re \zeta = n_B\left(\frac{\partial\delta\mu}{\partial x_H}\right)^{-1}
\frac{\partial\delta\mu}{\partial n_B}
\frac{\partial P}{\partial x_H}
\frac{\gamma}{\omega^2+\gamma^2},
\qquad
\gamma = B\lambda,\quad
B \equiv \frac{1}{n_B\,\frac{\partial\delta\mu}{\partial x_H}} .
\]
Here \(x_H=n_H/n_B\) with \(n_H=n_{\Sigma^-}+n_\Lambda\), \(\omega\) is the oscillation frequency, \(\lambda\) is the microscopic weak equilibration strength, \(B\) is a susceptibility, and the viscosity peaks when
\[
\gamma \approx \omega .
\]
This resonant structure is central throughout the subject [2009.05181].

Gusakov and Kantor derived the corrected non-superfluid formula for nucleon-hyperon matter as
\[
{\rm Re}\,\xi = -\frac{n_b^2}{\lambda}\,
\frac{\left(\frac{\partial P}{\partial x_H}\right)\left(\frac{\partial \delta\mu}{\partial n_b}\right)}
{\left(\frac{\partial \delta\mu}{\partial x_H}\right)^2}
\,
\frac{1}{1+\omega^2\tau^2},
\qquad
\tau \equiv \frac{n_b}{\lambda}
\left(\frac{\partial \delta\mu}{\partial x_H}\right)^{-1}.
\]
They emphasized that earlier work using the neutron fraction \(x_n\) as an independent variable neglected the source term associated with the fast reaction \(n+\Lambda\leftrightarrow p+\Sigma^{-}\), which is conceptually incorrect even if the numerical difference is typically only \(\sim 10\%-30\%\) [0806.4914].

In the \(\Lambda\Xi^{-}\) formulation used for instability-window calculations, the same resonant structure appears in the form
\[
\mathrm{Re}\,\zeta = \frac{z_{\max}\, 2(\lambda/\lambda_{\max})}{1+(\lambda/\lambda_{\max})^2},
\]
with
\[
z_{\max} \equiv \frac{n_b}{2\omega}
\left(\frac{\partial P}{\partial y_s}\right)
\left(\frac{\partial \Delta\mu}{\partial n_b}\right)
\left(\frac{\partial \Delta\mu}{\partial y_s}\right)^{-1},
\qquad
\lambda_{\max} \equiv n_b \omega \left(\frac{\partial \Delta\mu}{\partial y_s}\right)^{-1}.
\]
The meaning is standard: if reactions are too slow, composition is effectively frozen and \(\zeta\propto\lambda\); if they are too fast, matter remains near equilibrium and \(\zeta\propto 1/\lambda\); the maximum occurs when the equilibration rate matches the oscillation frequency [1911.08407].

## 4. Microscopic rate calculations: contact interaction, one-meson exchange, and finite-temperature phase space

The main modern advance in the merger-regime calculation is the replacement of the Fermi-surface approximation by a direct numerical evaluation of the full finite-temperature phase-space integral,
\[
\Gamma_{12\to 34} = \frac{1}{S}\int \frac{d^3p_1}{(2\pi)^3} \frac{d^3p_2}{(2\pi)^3}
\frac{d^3p_3}{(2\pi)^3} \frac{d^3p_4}{(2\pi)^3}
\frac{\sum_s |M_{1234}|^2}{2^4 E_1^*E_2^*E_3^*E_4^*}
(2\pi)^4 \delta(E_1+E_2-E_3-E_4) \delta^3(\mathbf p_1+\mathbf p_2-\mathbf p_3-\mathbf p_4)
\]
\[
\times f_1(E_1,\mu_1)f_2(E_2,\mu_2)\,[1-f_3(E_3,\mu_3)][1-f_4(E_4,\mu_4)],
\]
with \(S=2\) for reactions containing two identical baryons on one side and
\[
f_i(E_i,\mu_i)=\frac{1}{1+\exp[(E_i-\mu_i)/T]}.
\]
After analytic simplification, the remaining multidimensional integral is evaluated with the CUBA integration library. This is required because thermal smearing is large at merger temperatures, relativistic kinematics matter, the Fermi-surface approximation is invalid, and hyperons may be thermally populated below their \(T=0\) onset density [2009.05181].

The matrix element is treated in two ways. The contact interaction corresponds to integrated-out \(W\)-exchange between baryons. The one-meson-exchange (OME) treatment instead resolves a weak flavor-changing vertex in one baryon and a strong meson exchange with the second baryon, with
\[
M_{1234} = \left[ \bar u_3 F^S_{23}u_2\,\bar u_4F^W_{14}u_1\,D_\varphi(k_1^2)
- \bar u_3 F^S_{13}u_1\,\bar u_4F^W_{24}u_2\,D_\varphi(k_2^2) \right],
\]
\[
F_{ij}^W = G_F m_\pi^2(A_{ij}+B_{ij}\gamma_5), \qquad
F_{ij}^S = g_{ij}\gamma_5, \qquad
D_\varphi(k)=\frac{1}{k_0^2-k^2-m_\varphi^2}.
\]
In both the hot-merger and cold-star literatures, OME is found to be more effective than contact \(W\)-exchange. The merger study concludes that OME is typically an order of magnitude faster than the contact channel and that restricting the calculation to contact interactions while omitting processes III and IV can underestimate the total rate by a factor of about \(10^3\) [2009.05181]. In the \(\Lambda\Xi^{-}\) context, the OME result for \(np\leftrightarrow\Lambda p\) is about \(7\)–\(15\) times larger than the contact estimate, and some reactions have no \(W\)-exchange contribution at all but do proceed through meson exchange [1911.08407].

The Fermi-surface approximation remains useful at low \(T\). For momentum-independent matrix elements it yields
\[
\Gamma_{12\to 34} = \frac{T^3|M_{1234}|^2 (2\pi)^5 2^3 S}{} I(\xi)\,Q^{(4)},
\qquad \xi \equiv \delta\mu/T,
\]
with
\[
I(\xi)=\frac{e^\xi(e^\xi-1)(4\pi^2\xi+\xi^3)}{6},
\qquad
\lim_{\xi\to 0}I(\xi)=\frac{2\pi^2}{3}.
\]
However, the merger analysis explicitly states that this approximation can overestimate the full rate by up to an order of magnitude and completely fails below onset, where thermal hyperons dominate [2009.05181].

## 5. Superfluid hyperonic matter and the multi-coefficient viscosity tensor

In superfluid nucleon-hyperon matter, bulk viscosity is no longer described by a single scalar transport coefficient. Gusakov and Kantor extended dissipative relativistic hydrodynamics to mixtures in which neutrons, protons, \(\Lambda\), and \(\Sigma^{-}\) hyperons are all superfluid [0806.4914]. The particle currents are
\[
j_{(i)}^\mu = n_i u^\mu + Y_{ik} w_{(k)}^\mu,
\qquad
j_{(l)}^\mu = n_l u^\mu,\quad l=e,\mu,
\]
and the energy-momentum tensor contains both normal and superfluid contributions,
\[
T^{\mu\nu} = (P+\varepsilon)u^\mu u^\nu + P\eta^{\mu\nu}
+ Y_{ik}\!\left[ w_{(i)}^\mu w_{(k)}^\nu +\mu_i w_{(k)}^\mu u^\nu +\mu_k w_{(i)}^\nu u^\mu \right] +\tau^{\mu\nu}.
\]

The superfluid four-vectors satisfy the potentiality condition
\[
w_{(i)}^\mu = \partial^\mu \phi_i - q_i A^\mu -(\mu_i+\varkappa_i)u^\mu,
\]
or equivalently
\[
\partial^\nu\!\left[ w_{(i)}^\mu + q_iA^\mu + (\mu_i+\varkappa_i)u^\mu \right]
=
\partial^\mu\!\left[ w_{(i)}^\nu + q_iA^\nu + (\mu_i+\varkappa_i)u^\nu \right].
\]

Before imposing constraints, the dissipative expansion contains 25 bulk-viscosity coefficients. After using quasineutrality and the symmetry between protons and \(\Sigma^{-}\), the number reduces to 16. Onsager relations,
\[
\xi_{3\Lambda}=\xi_{5n},\qquad
\xi_{3\Sigma p}=\xi_{6n},\qquad
\xi_{4n}=\xi_{1n},\qquad
\xi_{5\Sigma p}=\xi_{6\Lambda},\qquad
\xi_{4\Lambda}=\xi_{1\Lambda},\qquad
\xi_{4\Sigma p}=\xi_{1\Sigma p},
\]
reduce the count to 10 independent coefficients. The fast-reaction equilibrium condition and entropy-production constraints then imply that only three are independent; a convenient choice is
\[
\xi_2,\qquad \xi_{1n},\qquad \xi_{1\Lambda}.
\]
This result is one of the key structural differences between normal and superfluid hyperonic matter [0806.4914].

Entropy production takes the compact form
\[
T\,\partial_\mu S^\mu = \delta\mu\,\Delta\Gamma = \lambda\,\delta\mu^2,
\]
which guarantees positivity and leads to the quadratic determinant constraints among the coefficients. A plausible implication is that any realistic calculation of oscillation damping in superfluid hyperonic cores must treat bulk viscosity as a coupled matrix in the superfluid fluxes rather than as a single scalar parameter.

## 6. Astrophysical regimes: mergers, \(r\)-modes, and instability windows

The astrophysical role of hyperonic bulk viscosity is regime dependent. For neutron-star mergers, the full finite-temperature OME-dominated calculation finds that, at \(T\gtrsim 1\)–\(2\) MeV and densities above saturation, the effective relaxation rate satisfies \(\gamma \gg \omega\) for oscillation frequencies around \(1\ \mathrm{kHz}\). Equilibration is therefore too fast to resonate with the macroscopic oscillation, and the bulk viscosity is small. The calculated dissipation times are seconds and longer in the MeV range. When the temperature is reduced to the keV range, the resonance condition
\[
\gamma \sim \omega
\]
can be met; the study finds a resonance around
\[
T \approx 4\ \mathrm{keV},
\]
where dissipation times fall to milliseconds. The same work concludes that hyperon bulk viscosity for temperatures in the MeV regime can probably be neglected in neutron-star mergers, but becomes highly relevant for keV-range temperatures [2009.05181].

In rotating neutron stars, the situation is different because the temperatures are much lower and the mode frequencies can lie close to the chemical relaxation rates. The \(r\)-mode stability problem is commonly formulated through
\[
\frac{1}{\tau_r(\Omega,T)}=
\frac{1}{\tau_{GR}(\Omega)}
+\frac{1}{\tau_H(\Omega,T)}
+\frac{1}{\tau_U(\Omega,T)}
+\frac{1}{\tau_\eta(\Omega,T)},
\]
with the critical angular velocity determined by
\[
\frac{1}{\tau_r(\Omega_C,T)}=0.
\]
In the modified chiral effective RMF model of rotating hyperonic stars, hyperonic bulk viscosity caused by non-leptonic weak interactions completely suppresses \(r\)-modes between \(\sim 10^8\) K and \(\sim 10^9\) K, producing a stable region between a low-temperature and a high-temperature instability window [2208.14436]. In the earlier effective chiral model, the instability remains significant around \(10^8\) K and the minimum of the critical curve is reported near \(T\approx 5\times 10^{10}\) K with \(\Omega_C\approx 0.04\,\Omega_K\) [1002.4253]. These contrasting outcomes reflect differences in the equation of state, composition, and rate model rather than disagreement about the underlying mechanism.

More recent relativistic work emphasizes two amplifying effects: hyperons open fast non-leptonic channels whose rates scale roughly as \((T/T_F)^2\), and relativistic nonbarotropic \(r\)-modes can enhance the dissipation by up to \(\sim 10^2\)–\(10^3\) compared with Newtonian estimates in the relevant region [2508.04226]. For the FSU2H EOS, stabilization of observed low-mass X-ray binaries is reported to be possible at roughly
\[
M \gtrsim 1.6\,M_\odot,
\]
while for TM1C it generally requires around
\[
M \gtrsim 1.8\,M_\odot
\]
once pairing is included [2508.04226]. In the OME-based \(\Lambda\Xi^{-}\) analysis, the bulk-viscosity peak occurs at temperatures around
\[
T_{\rm opt}\sim (0.5-1)\times 10^8\,{\rm K}
\]
in the unpaired case, which is close to the inferred internal temperatures of many low-mass X-ray binaries [1911.08407].

## 7. Limitations, controversies, and open directions

Several limitations recur across the literature. The merger calculation neglects superfluidity, justified at MeV temperatures but not at keV temperatures; neglects neutrino trapping, which may become important above roughly \(5\) MeV; treats only strangeness-changing non-leptonic reactions; and is restricted to the linear, subthermal regime [2009.05181]. The \(\Lambda\Xi^{-}\) \(r\)-mode study ignores baryon pairing in the microphysical rates, includes only the lightest exchanged mesons, and does not fully incorporate all hyperon species or superfluid hydrodynamic effects [1911.08407]. The superfluid hydrodynamic theory identifies another missing ingredient for detailed applications: the entrainment matrix \(Y_{ik}\) for superfluid nucleon-hyperon matter had not yet been calculated in that framework [0806.4914].

There are also model-dependent controversies about which hyperons appear first and which reactions dominate. Older calculations often assumed a \(\Sigma^{-}\Lambda\) composition and contact \(W\)-exchange rates; modern RMF equations of state frequently predict that \(\Lambda\) appears first, followed by \(\Xi^{-}\), while \(\Sigma^{-}\) appears much later or not at all [1911.08407]. This changes the reaction network and, correspondingly, the temperature of maximal damping. Likewise, the contact-interaction picture is now widely regarded as incomplete where OME channels are available, because OME can be substantially faster and can activate reactions absent in the contact description [2009.05181], [1911.08407].

A final source of variation is superfluid suppression. In the relativistic instability-window analysis, neutron superfluidity strongly suppresses weak reaction rates exponentially at low \(T\), while proton superconductivity has a smaller effect because it blocks only a subset of the relevant hyperonic reactions [2508.04226]. This suggests that quantitative stability boundaries remain sensitive to pairing models even when the qualitative role of hyperons is robust.

Taken together, the current literature supports a coherent picture. Hyperonic bulk viscosity is a resonance phenomenon controlled by the competition between compression timescales and strangeness-changing weak equilibration. In hot post-merger matter, equilibration is generally too rapid for strong dissipation at kilohertz frequencies [2009.05181]. In colder hyperonic neutron-star cores, the same reactions can dominate the damping of \(r\)-modes and substantially reshape instability windows, especially when meson-exchange microphysics, realistic hyperon compositions, relativistic mode structure, and superfluid suppression are treated consistently [2508.04226], [2208.14436], [1911.08407].

Source: https://www.emergentmind.com/topics/bulk-viscosity-in-hyperonic-matter