On a Class of Global Solutions to 3D Free-Boundary Relativistic Euler Equations with a Physical Vacuum Boundary

Technical Summary: On a Class of Global Solutions to 3D Free-Boundary Relativistic Euler Equations with a Physical Vacuum Boundary

Problem Statement
This work addresses the global-in-time dynamics of the relativistic Euler equations in Minkowski spacetime ($M^{1+3}$) for a gas with a physical vacuum boundary. The system models a relativistic fluid where the density $\rho$ vanishes at the free boundary $\partial \Omega_t$, distinguishing it from liquid models where density remains strictly positive. The study focuses on isentropic, barotropic, and polytropic gases governed by the equation of state $p = \rho^{1+\kappa}$ with $\kappa \in (0, 2/3]$.

The central challenge lies in establishing the existence of future-global solutions that expand asymptotically linearly in time. While local well-posedness for such systems has been recently established (e.g., Disconzi–Ifrim–Tataru [12]), the long-time behavior remains largely unexplored. Specifically, the authors investigate whether the expansion of the fluid bulk can suppress the formation of singularities (such as shocks or implosions) that typically plague compressible fluid dynamics, even in the presence of a vacuum boundary.

Methodology and Framework
The authors employ a rigorous energy method within a weighted functional framework, adapting techniques developed for classical compressible Euler equations (e.g., Jang–Masmoudi [27], Guo–Hadžić–Jang [18]) to the relativistic setting. The methodology involves several critical steps:

1. Mass-Critical Scaling Ansatz:
   Unlike the classical compressible Euler equations, the relativistic system lacks a scaling symmetry in the variables $(\rho, u^\mu)$ due to the nonlinear constraint $m_{\mu\nu}u^\mu u^\nu = -1$ and the specific scaling of the pressure term. To overcome this, the authors enforce a "mass-critical" scaling ansatz. They introduce a time-dependent scaling factor $\lambda(t)$ satisfying an ODE $\partial_t^2 \lambda = \delta \lambda^{-3\kappa-1}$, which models asymptotically linear expansion ($\lambda(t) \sim t$). The unknowns are rescaled as $\rho(t,r) = \lambda(s)^{-3}\tilde{\rho}(s,z)$ and $v(t,r) = \lambda(s)^{-3/2\kappa}\tilde{v}(s,z)$, where $z = r/\lambda(t)$.

2. Lagrangian Formulation and Modified Velocity:
   The authors introduce a modified velocity $V$ to capture perturbations around the background expansion. They transition to Lagrangian coordinates defined by the flow map $\eta$ generated by $V$. A key innovation is the treatment of the mass equation. In the classical case, the quantity $fF$ (related to the Jacobian) is conserved. In the relativistic case, this quantity evolves non-trivially. The authors introduce a "relativistic corrector" $G$ (and a related variable $\bar{G}$) to reformulate the mass equation. This allows them to express the Lagrangian density $f$ almost explicitly in terms of fundamental quantities, effectively decoupling the mass evolution from the top-order derivatives of the momentum equation.

3. Perturbation Analysis and Energy Estimates:
   The problem is reduced to analyzing a perturbation variable $\Theta = \eta - \zeta$ around the expanding background. The authors derive a high-order system of equations for $\Theta$ involving a degenerate elliptic operator $L_0$ and various error terms.

   • Decomposition of Errors: The forcing terms are decomposed into "Type I" errors (which have favorable decay but may lose derivatives) and "Type II" errors (which are low-frequency, large-size errors without decay).
   • Cancellation Structures: A significant portion of the technical work involves identifying precise cancellation structures in the highest-order terms, particularly those arising from the relativistic correction $G$ and the velocity component $u^0$. These cancellations are essential to avoid derivative loss when integrating by parts.
   • Weighted Norms: The analysis utilizes a hierarchy of weighted Sobolev norms $E_N$ and damping functionals $D_N$, incorporating admissible weight functions $w$ that vanish at the vacuum boundary to handle the degeneracy.

Key Contributions and Results
The primary result is the construction of an open class of initial data that leads to future-global solutions.

• Theorem 1.1 (Main Theorem): For $\kappa \in (0, 2/3]$ and initial expanding rates $\lambda_1 \in (0, 1)$, there exists an open set of spherically symmetric initial data (supported in a ball $B(0, \lambda_0)$ with small initial mass $\delta^{1/\kappa}$) that generates global solutions.
• Asymptotic Behavior: These solutions expand linearly in time with an asymptotic rate $\bar{\lambda} \sim \lambda_1$. Crucially, the expansion rate can be arbitrarily close to the speed of light (since $\lambda_1 < 1$ is the only constraint imposed by causality).
• Stability: The solutions are shown to be nonlinearly stable under small perturbations. The authors prove that the expansion mechanism is sufficient to suppress singularity formation, even with the presence of a physical vacuum boundary and relativistic effects.

Significance and Claims
The paper claims to provide the first rigorous construction of an open class of initial data leading to future-global solutions for the relativistic Euler equations with a physical vacuum boundary.

• Beyond Perturbation Theory: The authors emphasize that their result is not merely a perturbation of existing classical results. The lack of scaling symmetry in the relativistic equations and the specific structure of the vacuum boundary necessitate a novel approach involving the relativistic corrector $G$ and specific cancellation mechanisms that do not exist in the classical isentropic Euler case.
• Causality Saturation: The result is significant because it covers the full regime of causality allowed by the theory ($\lambda_1 < 1$), demonstrating that global existence holds even for expansion rates arbitrarily close to the speed of light.
• Physical Context: While the specific expanding solutions do not model static stars, the authors note that understanding long-time dynamics and stability mechanisms in relativistic fluids is a prerequisite for broader mathematical studies of stellar stability, black hole formation, and supernovae. The solutions satisfy standard physical requirements, including non-negative energy density and subluminal sound speed.

In summary, the work bridges a gap in the mathematical theory of relativistic fluids by demonstrating that linear expansion acts as a robust stabilizing mechanism against singularity formation, even in the challenging regime of physical vacuum boundaries and near-light-speed expansion.