Homo-Energetic Solutions in Kinetic Theory

• Homo-energetic solutions are a class of velocity-shifted solutions to the Boltzmann–Rayleigh equation that incorporate spatial deformations through a prescribed shear field.
• They balance shear-induced transport with collision effects, leading to non-Maxwellian steady states under small shear and self-similar growth in high-shear regimes.
• Rigorous analysis using Markov semigroup theory, moment ODEs, and fixed-point arguments establishes well-posedness and the existence of stationary non-equilibrium measures.

Searching arXiv for the specified paper and closely related work on homoenergetic solutions and the Rayleigh–Boltzmann equation. Homoenergetic solutions are a particular class of solutions to the linear Boltzmann–Rayleigh equation in which the spatial dependence is absorbed into a velocity shift induced by a prescribed deformation field. In the formulation studied in “Homoenergetic solutions for the Rayleigh-Boltzmann equation: existence of a stationary non-equilibrium solution” (Miele et al., 24 Jan 2025), they are used to describe the dynamics of Boltzmann gases under mechanical deformations, with emphasis on simple shear. Their significance lies in the fact that their long-time behavior cannot be described by Maxwellian distributions and depends strongly on the homogeneity of the collision kernel. For cut-off collision kernels with homogeneity $\gamma \geq 0$, including Maxwell molecules ($\gamma=0$) and hard potentials with $0 \leq \gamma < 1$, the analysis establishes both well-posedness in the space of non-negative Radon measures and the existence of a stationary non-equilibrium solution under the non-equilibrium condition induced by shear (Miele et al., 24 Jan 2025).

1. Kinematic ansatz and reduction to a homogeneous problem

The starting point is the linear Boltzmann–Rayleigh equation for a one-particle distribution

$g=g(t,x,w)\ge 0$

on $$\mathbb{R}_+\times \mathbb{R}^3_x \times \mathbb{R}^3_w$$,

$$\partial_t g + w\cdot \nabla_x g = \mathcal{L}g(w).$$

A homoenergetic solution is introduced through the ansatz

$g(t,x,w)=f(t,v), \qquad v:=w-\Sigma(t,x),$

with deformation field

$\Sigma(t,x)=L(t)x,$

where the matrix $L(t)$ solves

$\dot L + L^2 = 0,\qquad L(0)=A\in M_3(\mathbb{R}).$

Substitution yields a spatially homogeneous equation for $\gamma=0$0,

$\gamma=0$1

For the classical choice of simple shear,

$\gamma=0$2

one has $\gamma=0$3 and therefore $\gamma=0$4 for all $\gamma=0$5. In this case the reduced equation becomes

$\gamma=0$6

This reduction isolates the competition between transport induced by shear and collisions encoded in $\gamma=0$7. The resulting equation is spatially homogeneous but not equilibrium-preserving in the Maxwellian sense. The paper identifies this as the mechanism behind the non-equilibrium stationary states studied later (Miele et al., 24 Jan 2025).

2. Collision operator and homogeneity of the kernel

The collision operator $\gamma=0$8 is obtained by linearizing the classical Boltzmann binary-collision operator around a Maxwellian background $\gamma=0$9. With post- and pre-collisional velocities denoted by $0 \leq \gamma < 1$0 and $0 \leq \gamma < 1$1, and with

$0 \leq \gamma < 1$2

the operator is written as

$0 \leq \gamma < 1$3

where the collision rule is

$0 \leq \gamma < 1$4

and the normalized Maxwellian is

$0 \leq \gamma < 1$5

Under Grad’s cut-off assumption,

$0 \leq \gamma < 1$6

The full homoenergetic equation under simple shear is then

$0 \leq \gamma < 1$7

The collision kernel is assumed homogeneous of degree $0 \leq \gamma < 1$8 in the relative speed,

$0 \leq \gamma < 1$9

so that for any $g=g(t,x,w)\ge 0$0,

$g=g(t,x,w)\ge 0$1

The summary distinguishes three regimes: $g=g(t,x,w)\ge 0$2 for Maxwell molecules, $g=g(t,x,w)\ge 0$3 for moderately hard potentials, and $g=g(t,x,w)\ge 0$4 for a hyperbolic-dominated regime. In the cases actually treated, namely $g=g(t,x,w)\ge 0$5, the qualitative balance between shear transport and collisions is central. The text states that $g=g(t,x,w)\ge 0$6 gives a perfect balance at large velocities, while $g=g(t,x,w)\ge 0$7 yields a collision-dominated regime (Miele et al., 24 Jan 2025).

3. Weak formulation and well-posedness in non-negative Radon measures

The well-posedness theory is developed for

$g=g(t,x,w)\ge 0$8

where $g=g(t,x,w)\ge 0$9 denotes the non-negative finite Radon measures. Weak solutions are defined by testing against $$\mathbb{R}_+\times \mathbb{R}^3_x \times \mathbb{R}^3_w$$0:

$$\mathbb{R}_+\times \mathbb{R}^3_x \times \mathbb{R}^3_w$$1

with $$\mathbb{R}_+\times \mathbb{R}^3_x \times \mathbb{R}^3_w$$2 the adjoint collision operator.

The main argument proceeds through the backward adjoint problem

$$\mathbb{R}_+\times \mathbb{R}^3_x \times \mathbb{R}^3_w$$3

which is then rewritten forward in time. The generator

$$\mathbb{R}_+\times \mathbb{R}^3_x \times \mathbb{R}^3_w$$4

is considered on the Banach space $$\mathbb{R}_+\times \mathbb{R}^3_x \times \mathbb{R}^3_w$$5 with a domain imposing sufficient decay or growth control. The paper shows that $$\mathbb{R}_+\times \mathbb{R}^3_x \times \mathbb{R}^3_w$$6 is a Markov generator, via the Hille–Yosida or Lumer–Phillips theorem, and therefore generates a positive, strongly continuous semigroup $$\mathbb{R}_+\times \mathbb{R}^3_x \times \mathbb{R}^3_w$$7 on $$\mathbb{R}_+\times \mathbb{R}^3_x \times \mathbb{R}^3_w$$8 (Miele et al., 24 Jan 2025).

Solutions are then defined by duality:

$$\mathbb{R}_+\times \mathbb{R}^3_x \times \mathbb{R}^3_w$$9

This construction automatically satisfies the weak formulation. Uniqueness follows by testing the difference of two solutions against the dual semigroup.

For $$\partial_t g + w\cdot \nabla_x g = \mathcal{L}g(w).$$0, an alternative approach is also described. One introduces the characteristic function

$$\partial_t g + w\cdot \nabla_x g = \mathcal{L}g(w).$$1

which satisfies a transport-plus-integral-operator equation. According to the summary, this equation is amenable to a Banach fixed-point argument in the space of characteristic functions. This suggests that the Maxwell case permits an analytically distinct treatment because the homogeneity of the kernel simplifies the Fourier-side structure (Miele et al., 24 Jan 2025).

4. Stationary non-equilibrium states under shear

The central stationary problem is the steady equation

$$\partial_t g + w\cdot \nabla_x g = \mathcal{L}g(w).$$2

In weak form, for every $$\partial_t g + w\cdot \nabla_x g = \mathcal{L}g(w).$$3,

$$\partial_t g + w\cdot \nabla_x g = \mathcal{L}g(w).$$4

The solution $$\partial_t g + w\cdot \nabla_x g = \mathcal{L}g(w).$$5 is therefore stationary in the co-moving shear frame but remains non-equilibrium because the shear deformation produces a persistent balance between transport and collision effects rather than relaxation to the background Maxwellian.

The construction differs in the two regimes considered in the paper.

Regime	Shear condition	Main compactness set
$$\partial_t g + w\cdot \nabla_x g = \mathcal{L}g(w).$$6	Small shear $$\partial_t g + w\cdot \nabla_x g = \mathcal{L}g(w).$$7	$$\partial_t g + w\cdot \nabla_x g = \mathcal{L}g(w).$$8
$$\partial_t g + w\cdot \nabla_x g = \mathcal{L}g(w).$$9	Arbitrary $g(t,x,w)=f(t,v), \qquad v:=w-\Sigma(t,x),$0	$g(t,x,w)=f(t,v), \qquad v:=w-\Sigma(t,x),$1

In the Maxwell case, the argument begins with closed ODEs for the second-order moment matrix

$g(t,x,w)=f(t,v), \qquad v:=w-\Sigma(t,x),$2

These satisfy

$g(t,x,w)=f(t,v), \qquad v:=w-\Sigma(t,x),$3

where $g(t,x,w)=f(t,v), \qquad v:=w-\Sigma(t,x),$4 is a linear operator and $g(t,x,w)=f(t,v), \qquad v:=w-\Sigma(t,x),$5 comes from loss-gain balance. For $g(t,x,w)=f(t,v), \qquad v:=w-\Sigma(t,x),$6, the spectrum of $g(t,x,w)=f(t,v), \qquad v:=w-\Sigma(t,x),$7 lies in the left half-plane, so there exists a unique stable stationary matrix

$g(t,x,w)=f(t,v), \qquad v:=w-\Sigma(t,x),$8

The set $g(t,x,w)=f(t,v), \qquad v:=w-\Sigma(t,x),$9 defined by mass one and fixed second moment is weak* compact in $\Sigma(t,x)=L(t)x,$0. The dual semigroup $\Sigma(t,x)=L(t)x,$1 leaves $\Sigma(t,x)=L(t)x,$2 invariant and is weak* continuous. Schauder’s fixed-point theorem then yields a nontrivial measure $\Sigma(t,x)=L(t)x,$3 with $\Sigma(t,x)=L(t)x,$4 for all $\Sigma(t,x)=L(t)x,$5, hence a steady solution (Miele et al., 24 Jan 2025).

For hard potentials with $\Sigma(t,x)=L(t)x,$6, the key input is a uniform-in-time bound on the second moment

$\Sigma(t,x)=L(t)x,$7

Using a Povzner-type estimate, the paper proves that for any $\Sigma(t,x)=L(t)x,$8,

$\Sigma(t,x)=L(t)x,$9

with $L(t)$0 for a suitable choice of constants. Therefore $L(t)$1 remains bounded. The weak* compact set $L(t)$2 is invariant under the semigroup, and Schauder’s theorem again yields a stationary solution.

In strong form, the stationary equation is

$L(t)$3

A plausible implication is that the existence mechanism depends less on explicit solvability than on the combination of invariant convex sets, moment control, and weak* compactness.

5. Maxwell molecules: small-shear stationarity and large-shear growth

A distinctive feature of the Maxwell case $L(t)$4 is the existence of two qualitatively different behaviors depending on the size of the shear parameter $L(t)$5. The paper states that there is a threshold $L(t)$6 separating small- and large-shear regimes (Miele et al., 24 Jan 2025).

For small shear, $L(t)$7, the operator governing the second-moment dynamics has strictly negative spectrum. As a consequence, all moments converge exponentially to those of the unique stationary matrix $L(t)$8. The summary further reports that numerically one can expand in powers of $L(t)$9 to obtain

$\dot L + L^2 = 0,\qquad L(0)=A\in M_3(\mathbb{R}).$0

so that $\dot L + L^2 = 0,\qquad L(0)=A\in M_3(\mathbb{R}).$1 is a small non-equilibrium perturbation of the Maxwellian. This does not mean the stationary state is Maxwellian; rather, it remains shear-driven and only approaches the Maxwellian configuration in a perturbative sense as $\dot L + L^2 = 0,\qquad L(0)=A\in M_3(\mathbb{R}).$2.

For large shear, $\dot L + L^2 = 0,\qquad L(0)=A\in M_3(\mathbb{R}).$3, the moment ODE matrix acquires a real eigenvalue

$\dot L + L^2 = 0,\qquad L(0)=A\in M_3(\mathbb{R}).$4

and the second moments grow like

$\dot L + L^2 = 0,\qquad L(0)=A\in M_3(\mathbb{R}).$5

The summary explicitly states that no bounded steady solution can exist in this regime. Instead, a self-similar asymptotic behavior is conjectured:

$\dot L + L^2 = 0,\qquad L(0)=A\in M_3(\mathbb{R}).$6

for suitable $\dot L + L^2 = 0,\qquad L(0)=A\in M_3(\mathbb{R}).$7, where $\dot L + L^2 = 0,\qquad L(0)=A\in M_3(\mathbb{R}).$8 solves

$\dot L + L^2 = 0,\qquad L(0)=A\in M_3(\mathbb{R}).$9

This dichotomy is one of the paper’s main conceptual outcomes. It indicates that, for Maxwell molecules, shear can be either weak enough to permit a stationary non-equilibrium balance or strong enough to drive unbounded moment growth. The summary describes this as a bifurcation at $\gamma=0$00 between steady-state relaxation and unbounded, self-similar growth. This suggests that the precise balance implied by $\gamma=0$01 is structurally marginal: neither unequivocally collision-dominated nor purely transport-dominated.

6. Conceptual significance and relation to non-equilibrium kinetic theory

The results place homoenergetic solutions in a framework where mechanical deformation is built directly into the ansatz, rather than treated as a perturbation of an equilibrium evolution. In this setting, the long-time behavior is governed by the interaction between the deformation field and the homogeneity of the collision kernel. The paper’s conclusions may be summarized as follows: well-posedness holds in the space of non-negative Radon measures; a stationary non-equilibrium state exists when collisions dominate sufficiently, namely for $\gamma=0$02 and arbitrary $\gamma=0$03, or for $\gamma=0$04 with small shear $\gamma=0$05; and the Maxwell case exhibits a transition to a non-stationary self-similar regime when $\gamma=0$06 (Miele et al., 24 Jan 2025).

A common misconception would be to identify stationary solutions of the reduced equation with equilibrium Maxwellians. The analysis does not support that interpretation. The stationary objects constructed here satisfy a balance law in the sheared frame and arise specifically under a non-equilibrium condition induced by the presence of the shear deformation. Even in the perturbative Maxwell case, where $\gamma=0$07 numerically, the stationary state is described as a non-equilibrium perturbation rather than equilibrium itself.

Another potential misconception is that the cut-off assumption or the linearization around a Maxwellian background trivializes the dynamics. The summary instead shows that the combination of transport, collision homogeneity, and moment propagation leads to distinct regimes requiring different tools: Markov semigroup theory, Fourier-transform methods, moment ODEs, Povzner-type estimates, and Schauder fixed-point arguments. This suggests that even in the linear Rayleigh–Boltzmann setting, homoenergetic solutions retain a rich non-equilibrium structure.

Within the scope described in the paper, homoenergetic solutions provide a mathematically controlled setting for studying non-Maxwellian steady states and deformation-driven asymptotics in kinetic theory. The analysis is limited to simple shear and cut-off collision kernels with $\gamma=0$08, particularly Maxwell molecules and hard potentials with $\gamma=0$09. A plausible implication is that extensions to other deformation fields or to different homogeneity classes would require new mechanisms for moment control and compactness beyond those established here.