Necessary and sufficient conditions of nonlinear causality in viscous anisotropic hydrodynamics

Technical Summary: Necessary and Sufficient Conditions of Nonlinear Causality in Viscous Anisotropic Hydrodynamics

Problem Statement
Relativistic hydrodynamics is the standard framework for describing the space-time evolution of matter produced in relativistic heavy-ion collisions, such as the quark-gluon plasma (QGP). However, the validity of this description at the earliest stages of evolution remains an open question, particularly when the system is far from equilibrium. In these regimes, large pressure anisotropies arise due to rapid longitudinal expansion. Viscous Anisotropic Hydrodynamics (VAH) has been proposed as an extension of conventional viscous hydrodynamics (VH) to treat these anisotropies non-perturbatively. While VAH has shown phenomenological success in describing flow observables and capturing non-equilibrium attractors better than traditional frameworks, its theoretical consistency in the far-from-equilibrium regime is not guaranteed. Specifically, the nonlinear causality of the VAH equations of motion has not been rigorously established. Without such guarantees, the theory may permit superluminal signal propagation, rendering it unphysical in early-stage collision dynamics.

Methodology
The authors investigate the nonlinear causality of the VAH framework by analyzing the characteristic velocities of the equations of motion. The study operates under a specific approximation: it neglects next-to-leading-order (NLO) dissipative corrections ($\delta \tilde{f}$), focusing solely on the leading-order anisotropic distribution function ($f_a$) which incorporates the relaxation of pressure anisotropies.

The methodology proceeds as follows:

1. Formulation: The energy-momentum tensor is decomposed to separate longitudinal ($P_L$) and transverse ($P_\perp$) pressures, along with associated dissipative currents. The resulting equations of motion (Eqs. 5–9) form a first-order quasilinear system, expressed in matrix form as $A^\alpha(\Psi)\partial_\alpha \Psi = F(\Psi)$.
2. Characteristic Analysis: To determine signal propagation speeds, the authors derive the characteristic equation $\det(A^\alpha \xi_\alpha) = 0$, where $\xi_\alpha$ is the normal vector to the characteristic surface. By decomposing $\xi_\alpha$ in the local rest frame, they define the characteristic velocity $v_c$.
3. Derivation of Conditions: The characteristic equation reduces to a quartic equation in $v_c$, which becomes a quadratic equation in $v_c^2$. The condition for nonlinear causality is that all real solutions for $v_c^2$ must lie within the interval $[0, 1]$ (ensuring subluminal propagation) for any propagation direction.
4. Mathematical Constraints: The authors translate the requirement that the quadratic function $f(v_c^2)$ has roots within $[0, 1]$ into a set of four algebraic conditions involving the function's value at boundaries, its symmetry axis, and its discriminant. By systematically solving these inequalities for arbitrary propagation directions, they derive the necessary and sufficient constraints on the transport coefficients.

Key Contributions and Results
The primary contribution of this work is the derivation of a closed set of necessary and sufficient conditions that ensure nonlinear causality in the VAH framework. These conditions are expressed as inequalities involving the transport coefficients ($\bar{\zeta}$) and the thermodynamic variables (energy density $E$ and pressures $P_L, P_\perp$):

1. Diagonal Constraints:
   $$-1 \leq \frac{\bar{\zeta}^L_z}{E + P_L} \leq 0 \quad \text{and} \quad -1 \leq \frac{\bar{\zeta}^\perp_\perp}{E + P_\perp} \leq 0$$
   These inequalities ensure that the characteristic velocities for purely longitudinal and purely transverse propagation remain causal.

2. Coupling Constraints:
   $$0 \leq \frac{\bar{\zeta}^L_\perp}{E + P_\perp} \frac{\bar{\zeta}^\perp_z}{E + P_L} \leq \frac{\bar{\zeta}^L_z}{E + P_L} \frac{\bar{\zeta}^\perp_\perp}{E + P_\perp} \leq 1$$
   These conditions govern the mixing between longitudinal and transverse modes. They ensure that when coupling terms are non-zero, the mixed eigenmodes remain real and causal.

The authors demonstrate that these conditions admit a clear physical interpretation: the breaking of isotropic symmetry splits the system into two distinct eigenmodes. The derived constraints ensure that for all propagation directions, the squared characteristic velocities of these mixed modes are real and non-negative, thereby maintaining both thermodynamic stability and causality.

Significance and Claims
The paper claims that these results establish the regime of validity for VAH, providing a theoretical foundation for its application to the early-time dynamics of relativistic heavy-ion collisions. The authors highlight that, unlike the complex conditions previously derived for conventional viscous hydrodynamics, the conditions for VAH take a "remarkably simple form." This simplicity is attributed to the neglect of NLO corrections in the current analysis.

The study concludes that while VAH offers a theoretically consistent way to treat anisotropic expansion non-perturbatively, its application is strictly limited to regions of parameter space where the derived inequalities are satisfied. The authors note that future work is required to incorporate the NLO corrections ($\delta \tilde{f}$) to establish causality conditions for the complete VAH framework and to perform numerical simulations comparing VAH with conventional VH to further validate the hydrodynamic description in early-stage collisions.