Nonlinear diffusion in relativistic kinetic theory

Technical Summary: Nonlinear Diffusion in Relativistic Kinetic Theory

Problem Statement
The paper addresses a fundamental inconsistency in relativistic kinetic theory regarding diffusion processes. While linear kinetic diffusion equations (such as the relativistic Kramers equation) have been proposed to model particle interactions, existing Lorentz-invariant formulations fail to preserve the conservation laws of energy and momentum. This violation renders such equations incompatible with General Relativity, as the conservation of the stress-energy tensor ($\nabla_\beta T^{\alpha\beta} = 0$) is a necessary constraint for the existence of solutions to the Einstein equations. Previous attempts to resolve this, such as introducing a cosmological scalar field to transfer energy-momentum, modify the standard Einstein equations. The author seeks a formulation that is Lorentz invariant, conserves particle number, energy, and momentum, and remains compatible with the standard Einstein equations without introducing additional matter fields.

Methodology
The paper employs a constructive approach, first establishing a consistent framework in Newtonian mechanics before generalizing it to Special and General Relativity.

1. Newtonian Framework (Section 2): The author begins with the linear kinetic diffusion equation (Fokker-Planck/Kramers equation) and identifies its failure to conserve energy and momentum. A nonlinear modification is introduced where the diffusion scalar field ($\Phi_f$) and the drift vector field ($W_f$) depend functionally on the macroscopic moments of the distribution function (particle density $n_f$, current $j_f$, and energy density $e_f$). Specifically, the drift is set proportional to the difference between the particle velocity and the local mean velocity, and the diffusion coefficient is adjusted to ensure the conservation laws hold. This nonlinear model is shown to be Galilean invariant and to converge to the standard Maxwellian distribution in equilibrium.

2. Special Relativistic Generalization (Section 3): The methodology is extended to Minkowski space. The author postulates a Lorentz-invariant equation for the kinetic density $f(x, u)$ defined on the mass shell. To ensure Lorentz invariance, the diffusion and drift terms are constructed using four-vectors and tensors that transform correctly under Lorentz transformations. The key step involves determining the specific forms of the diffusion scalar field $\Phi$ and the drift vector field $K$ such that the local conservation laws for particle number ($\nabla_\alpha N^\alpha = 0$) and energy-momentum ($\nabla_\beta T^{\alpha\beta} = 0$) are satisfied. The author derives explicit expressions for these fields in terms of the macroscopic moments ($N^\alpha$ and $T^{\alpha\beta}$) and a drag coefficient $\mu$.

3. General Relativistic Coupling (Section 4): The derived equation is generalized to a curved Lorentzian manifold. The partial derivatives are replaced by covariant derivatives involving the Levi-Civita connection and Ricci rotation coefficients. The author demonstrates that the resulting equation preserves the contracted Bianchi identities, thereby allowing it to be coupled directly to the Einstein field equations without modifying the gravitational sector.

Key Contributions and Results

• Nonlinear Lorentz Invariant Equation: The paper introduces a specific nonlinear kinetic diffusion equation (Eq. 3.8 and 3.18) that is consistent with the conservation of particle number, energy, and momentum. Unlike previous linear relativistic models, this equation does not violate conservation laws.
• Equilibrium Distribution: The equilibrium solution ($G_{eq}$) of the new nonlinear equation is derived. A significant result is that this equilibrium distribution is not the standard Jüttner distribution commonly employed in relativistic kinetic theory (which typically describes the equilibrium of linear relativistic kinetic models). Instead, for a particle distribution with zero average momentum, the equilibrium distribution is given by:
  $$G_{eq}(u) = \frac{c^{-3}Z(\lambda)}{(1 + c^{-2}|u|^2)^\lambda}$$
  where $\lambda$ is a dimensionless constant and $Z(\lambda)$ is a normalization factor. The paper proves that this distribution converges to the standard Maxwellian distribution in the Newtonian limit ($c \to \infty$).
• Compatibility with General Relativity: The primary theoretical contribution is the demonstration that this nonlinear diffusion equation satisfies $\nabla_\beta T^{\alpha\beta} = 0$. Consequently, it can be coupled to the Einstein equations ($R_{\alpha\beta} - \frac{1}{2}g_{\alpha\beta}R = \frac{8\pi G}{c^2}T_{\alpha\beta}$) without the need for auxiliary scalar fields or modifications to the Einstein tensor.
• Newtonian Limit: The relativistic model is shown to reduce to the nonlinear Newtonian model derived in Section 2 when $c \to \infty$.
• Spatially Homogeneous Solutions: The paper analyzes spatially homogeneous solutions, showing that while the relativistic equation remains nonlinear (unlike its Newtonian counterpart), it possesses time-independent solutions that converge to the non-central Maxwellian in the Newtonian limit.

Significance and Claims
The paper claims to resolve a long-standing incompatibility between kinetic diffusion models and General Relativity. By constructing a nonlinear diffusion equation that inherently respects the conservation of energy and momentum, the author provides a framework where kinetic theory and gravity can be coupled consistently within the standard Einstein theory. The work suggests that the standard Jüttner distribution may not be the correct equilibrium state for systems governed by such diffusion processes (specifically in the zero-momentum frame), proposing a new equilibrium distribution that depends on the specific diffusion parameters. The author notes that applications to cosmology will be discussed in future publications, but the current work focuses on the mathematical consistency and derivation of the equations.