Linear causality and stability constraints on relativistic second-order magnetohydrodynamics

Abstract: In this work, we construct a theoretical framework for relativistic second-order magnetohydrodynamics based on entropy current analysis. The formalism consistently incorporates the relaxation dynamics of dissipative fluxes, ensuring the hyperbolic nature of the evolution equations. Utilizing linear mode analysis, we investigate the constraints imposed by causality and stability on this anisotropic system. By linearizing the theory around a homogeneous equilibrium state, we demonstrate that the excitation spectrum decomposes into magnetosonic, Alfvén, and charge-diffusion sectors. For each sector, we derive asymptotic dispersion relations in both the long-wavelength (small-$k$) and short-wavelength (large-$k$) regimes, validating them against exact numerical roots. Our numerical analysis confirms the accuracy of these asymptotic solutions and uncovers a nontrivial angular dependence, especially near special propagation directions where the ordinary momentum expansion becomes less reliable. By evaluating the large-$k$ behavior of the propagating branches alongside the damping properties of non-hydrodynamic modes, we delineate the corresponding causality constraints. We find that the admissible causal domain is governed by the interplay between anisotropic transport coefficients and relaxation times, with the resulting bounds being intrinsically mode-dependent. These findings provide a systematic theoretical foundation for developing stable and causal relativistic magnetohydrodynamics beyond the first-order approximation.

• The paper derives an entropy-consistent Israel–Stewart RMHD framework with seven dynamical dissipative fluxes, anisotropic transport channels, and nonnegative relaxation times constrained by entropy production.
• The paper decomposes the linearized 19-variable system into magneto-sonic, Alfvén, and charge-diffusion sectors, showing that causality depends on correlated transport and relaxation parameters rather than independent coefficient bounds.
• The paper identifies the charge-sector condition κ sin²θ ≤ n̄τν, confirms stability numerically for representative parameters, and shows that generic-angle analytic stability, nonlinear well-posedness, and microscopic second-order couplings remain open problems.

This paper develops a relativistic second-order magnetohydrodynamic (RMHD) framework derived from entropy-current analysis and characterizes its linear causality and stability properties in the presence of a uniform background magnetic field (2606.00649). The work addresses a gap in the literature: while second-order viscous and resistive RMHD formulations exist from kinetic theory, a systematic, mode-resolved determination of how anisotropic transport coefficients and relaxation times jointly constrain the linear causal domain had remained incomplete.

Entropy-current construction of ideal and first-order RMHD

The theory is built on the conservation of energy-momentum, particle number, and magnetic flux ($\partial_\mu \tilde F^{\mu\nu}=0$), with the inhomogeneous Maxwell equation treated as a constitutive relation rather than an independent conservation law—consistent with electric screening in a conducting plasma. Working in the Landau–Lifshitz frame for a parity-even plasma, the authors use the entropy-current principle to derive the ideal constraints, including the generalized thermodynamic identity $\epsilon + p_\perp = Ts + \mu n + B^\mu H_\mu$ and the pressure-anisotropy relation $p_\perp - p_\parallel = B^\mu H_\mu$, which identifies $1/\mu_m = (p_\perp - p_\parallel)/B^2$ as the inverse magnetic permeability.

At first order, Curie's principle splits the dissipative fluxes into parity-even $$\{\Pi_\perp, \Pi_\parallel, f^\mu, \pi_\perp^{\mu\nu}\}$$ and parity-odd $\{\nu_\perp^\mu, \ell^\mu, m_\perp^{\mu\nu}\}$ sectors, yielding anisotropic bulk viscosities with cross-couplings $\zeta_\times$, $\zeta_\times'$, coupled charge-diffusion/induction coefficients $\kappa$, $\kappa_\times$, $\epsilon + p_\perp = Ts + \mu n + B^\mu H_\mu$0, and transverse shear $\epsilon + p_\perp = Ts + \mu n + B^\mu H_\mu$1. Positivity of entropy production imposes positive-semi-definiteness conditions such as $\epsilon + p_\perp = Ts + \mu n + B^\mu H_\mu$2 (in the charge-conjugate symmetric limit) and $\epsilon + p_\perp = Ts + \mu n + B^\mu H_\mu$3.

Israel–Stewart-type second-order completion

The second-order extension promotes all seven dissipative fluxes to independent dynamical variables via a quadratic correction $\epsilon + p_\perp = Ts + \mu n + B^\mu H_\mu$4 to the entropy current, containing diagonal terms with coefficients $\epsilon + p_\perp = Ts + \mu n + B^\mu H_\mu$5 plus cross-coupling structures $\epsilon + p_\perp = Ts + \mu n + B^\mu H_\mu$6 and $\epsilon + p_\perp = Ts + \mu n + B^\mu H_\mu$7. The stability requirement $\epsilon + p_\perp = Ts + \mu n + B^\mu H_\mu$8 fixes the signs of the diagonal coefficients and hence non-negative relaxation times $\epsilon + p_\perp = Ts + \mu n + B^\mu H_\mu$9, $p_\perp - p_\parallel = B^\mu H_\mu$0, etc., while the off-diagonal coefficients remain unconstrained by the second law alone—a limitation the authors acknowledge, noting that their values require microscopic input from kinetic theory. Enforcing $p_\perp - p_\parallel = B^\mu H_\mu$1 yields full Israel–Stewart-type relaxation equations with additional anisotropic channels absent in unmagnetized systems.

Linear mode structure

Linearizing around homogeneous equilibrium with $p_\perp - p_\parallel = B^\mu H_\mu$2 along $p_\perp - p_\parallel = B^\mu H_\mu$3 reduces the system to a $p_\perp - p_\parallel = B^\mu H_\mu$4 secular problem that block-decomposes into magneto-sonic ($p_\perp - p_\parallel = B^\mu H_\mu$5), Alfvén ($p_\perp - p_\parallel = B^\mu H_\mu$6), and charge-diffusion ($p_\perp - p_\parallel = B^\mu H_\mu$7) sectors, with residual couplings controlled by $p_\perp - p_\parallel = B^\mu H_\mu$8, $p_\perp - p_\parallel = B^\mu H_\mu$9, and $1/\mu_m = (p_\perp - p_\parallel)/B^2$0. Under a conformal equation of state ($1/\mu_m = (p_\perp - p_\parallel)/B^2$1) and charge-conjugate-symmetric background ($1/\mu_m = (p_\perp - p_\parallel)/B^2$2, $1/\mu_m = (p_\perp - p_\parallel)/B^2$3), the sectors fully decouple.

The small-$1/\mu_m = (p_\perp - p_\parallel)/B^2$4 expansion shows that hydrodynamic branches emerge from zero modes at $1/\mu_m = (p_\perp - p_\parallel)/B^2$5: fast/slow magneto-sonic pairs with angle-dependent velocities, an Alfvén pair $1/\mu_m = (p_\perp - p_\parallel)/B^2$6, and a diffusive charge mode $1/\mu_m = (p_\perp - p_\parallel)/B^2$7. All remaining poles are gapped relaxation modes at $1/\mu_m = (p_\perp - p_\parallel)/B^2$8 receiving perturbative gradient corrections.

Large-$1/\mu_m = (p_\perp - p_\parallel)/B^2$9 asymptotics and causality conditions

The large-$$\{\Pi_\perp, \Pi_\parallel, f^\mu, \pi_\perp^{\mu\nu}\}$$0 expansion reveals the hyperbolic completion. In the charge sector, diffusion telegraphizes into a propagating pair with front velocity $$\{\Pi_\perp, \Pi_\parallel, f^\mu, \pi_\perp^{\mu\nu}\}$$1, yielding the transparent necessary causality bound $$\{\Pi_\perp, \Pi_\parallel, f^\mu, \pi_\perp^{\mu\nu}\}$$2. The magneto-sonic propagating speeds are roots of a cubic whose coefficients are lengthy functions of the bulk, shear, and induction transport data; the Alfvén sector similarly yields four propagating branches plus two finite relaxation poles with explicit angular-dependent damping rates. Symbolic verification of the Routh–Hurwitz-type positivity conditions succeeds only at the special angles $$\{\Pi_\perp, \Pi_\parallel, f^\mu, \pi_\perp^{\mu\nu}\}$$3 and $$\{\Pi_\perp, \Pi_\parallel, f^\mu, \pi_\perp^{\mu\nu}\}$$4; the generic-angle case is left analytically unverified due to coefficient complexity—an explicitly conceded limitation.

A notable angular finding is that the small-$$\{\Pi_\perp, \Pi_\parallel, f^\mu, \pi_\perp^{\mu\nu}\}$$5 expansion loses accuracy near $$\{\Pi_\perp, \Pi_\parallel, f^\mu, \pi_\perp^{\mu\nu}\}$$6 for the magneto-sonic and Alfvén sectors (consistent with the non-analytic angular structure identified in prior work on spin MHD), whereas the large-$$\{\Pi_\perp, \Pi_\parallel, f^\mu, \pi_\perp^{\mu\nu}\}$$7 expansion fails near $$\{\Pi_\perp, \Pi_\parallel, f^\mu, \pi_\perp^{\mu\nu}\}$$8 in the charge sector because $$\{\Pi_\perp, \Pi_\parallel, f^\mu, \pi_\perp^{\mu\nu}\}$$9 appears at order $\{\nu_\perp^\mu, \ell^\mu, m_\perp^{\mu\nu}\}$0; the limits $\{\nu_\perp^\mu, \ell^\mu, m_\perp^{\mu\nu}\}$1 and $\{\nu_\perp^\mu, \ell^\mu, m_\perp^{\mu\nu}\}$2 do not commute there.

Numerical causality analysis

Using a representative parameter set at $\{\nu_\perp^\mu, \ell^\mu, m_\perp^{\mu\nu}\}$3, exact numerical roots confirm both asymptotic expansions in their respective regimes, with a crossover window around $\{\nu_\perp^\mu, \ell^\mu, m_\perp^{\mu\nu}\}$4–$\{\nu_\perp^\mu, \ell^\mu, m_\perp^{\mu\nu}\}$5 where neither truncated expansion is uniformly accurate. Contour maps of the maximal asymptotic group velocity $\{\nu_\perp^\mu, \ell^\mu, m_\perp^{\mu\nu}\}$6 show that:

• Magneto-sonic sector: the luminal boundary is deformed nontrivially by the mixed bulk coupling $\{\nu_\perp^\mu, \ell^\mu, m_\perp^{\mu\nu}\}$7; causality cannot be inferred by scanning bulk coefficients independently. The constraint is a genuinely multi-parameter restriction among $\{\nu_\perp^\mu, \ell^\mu, m_\perp^{\mu\nu}\}$8, $\{\nu_\perp^\mu, \ell^\mu, m_\perp^{\mu\nu}\}$9, $\zeta_\times$0, $\zeta_\times$1, $\zeta_\times$2, and $\zeta_\times$3, often taking the form of correlated finite windows rather than independent upper bounds.
• Alfvén sector: pairwise scans among $\zeta_\times$4, $\zeta_\times$5, $\zeta_\times$6, $\zeta_\times$7 reveal smoother but still clearly correlated trade-off boundaries.
• Charge-diffusion sector: a single monotonic threshold in $\zeta_\times$8 reproduces the analytic condition exactly—the cleanest manifestation of Maxwell–Cattaneo telegraphization.

Throughout the scanned parameter space, all non-hydrodynamic poles remain in the lower-half complex $\zeta_\times$9-plane and phase velocities stay bounded, indicating no linear instabilities for the chosen parameters.

Limitations and open questions

Several restrictions qualify these results. The causal domain is established only as a necessary linear condition about a homogeneous equilibrium, not as well-posedness of the nonlinear theory. The diagonalized form of $\zeta_\times'$0 excludes inter-sectorial cross-couplings (e.g., $\zeta_\times'$1) for tractability, and the second-order coefficients $\zeta_\times'$2, $\zeta_\times'$3 lack microscopic determination. Analytic stability proofs are restricted to $\zeta_\times'$4; generic angles rest on numerical evidence. Anomalous/chiral transport is excluded by the parity-even assumption, and the residual gapless non-propagating root $\zeta_\times'$5 in the magneto-sonic sector is not analyzed further.

Conclusion

The paper provides a self-consistent entropy-current derivation of second-order RMHD together with a sector-resolved map of its linear causal structure. Its central result is that causality in magnetized relativistic fluids is governed by the correlated interplay of anisotropic transport coefficients and relaxation times within each propagation sector, rather than by any single coefficient—a conclusion quantified through matched analytical asymptotics and exact numerical spectra. Open problems include microscopic fixation of the second-order couplings, generalizations beyond static backgrounds, and extensions to chiral RMHD relevant to the chiral magnetic effect.