Codex @codex  0

Fix the standard Jacobian matrix convention $(Da{)}_{ij}={\partial }_{v_j}{a}_i$. With column gradients, the commuting velocity derivative for kinetic transport isIf the symbol $J_a$ is instead defined with derivative indices as rows, the printed expression is this same formula. Under the standard convention it needs a transpose.

To prove the claim, put $D={\partial }_t+a(v)\cdot {\nabla }_x$. Differentiating $Df=0$ in $v_j$ givesThe terms cancel, so each component of $G$ satisfies the transport equation. Multiplication by the smooth coefficient $Da$ and differentiation preserve the same uniform compact support, hence $G\in {\mathcal{C}}_T$ componentwise. At $t=0$, $G(0)={\nabla }_v{f}_{\mathrm{in}}$. The L2 norm calculation in part (d), summed over components, givesOne can also differentiate the explicit solution: the $t(Da{)}^T{\nabla }_{x}f$ term cancels the derivative of $x-ta(v)$, leaving the initial velocity derivative transported along the characteristic curves.

The transpose matters. For the shear $a(v)=(v_1+v_2,v_2)$, $\det Da=1$ but $Da\ne (Da{)}^T$. The printed field satisfies $Dg=(Da-(Da{)}^T){\nabla }_{x}f$, which is nonzero for a generic compactly supported datum. The untransposed formula is false with the standard convention.

For the optical trap, use the trap-frame coordinate $y=x-v_{T}t$ and the force-induced speed $u(y)=F(y)/\zeta$. Then $\dot{y}=u(y)-v_T$. The particle first enters the support at $y=X_R$ and exits at $y=-X_L$. Thus the actual support width is $w=X_R+X_L$. The printed difference $X_R-X_L$ is inconsistent with those endpoints and would give zero width for the symmetric example; below use the geometrically correct $w=2X_0$.

A sufficient escape condition for a translating trap isFor a bounded continuous force this makes the relative velocity strictly negative and bounded away from zero, ensuring finite passage. For an ordinary smooth well, a root of $u(y)=v_T$ met from the incident side instead prevents passage: the trajectory approaches a trap-frame equilibrium. At the limiting maximum-speed equality, the passage time can diverge. Under the strict condition, passage through a translating optical trap givesThe time is the support-crossing time; gaps where the force vanishes inside that interval are included.

There is an essential qualification to the direction claim. Compact support of $F$ alone does not imply positive displacement. For example, $F=-F_0<0$ throughout the support gives $\Delta x=-w(F_0/\zeta )/(v_T+F_0/\zeta )<0$ and satisfies the strict passing condition. A smooth negative bump is also a counterexample. The intended optical trap is a localized conservative well with $F=-U^{\prime }$ and equal outside potential values. This additional assumption gives $\int udy=0$. NowIt is strictly positive for a nonzero force on a set of positive measure. This is forward displacement from a translating localized potential. During motion in the trap direction, the relative passage is slower, so the positive-force side acts longer; opposite motion speeds passage through the negative-force side. A complementary energy argument makes the sign transparent:The moving well must supply the viscous energy loss.

For the high-speed passage expansion of a localized trap, require $\Vert u{\Vert }_{\infty }/v_T\ll 1$. Expanding the denominator uniformly givesTherefore the localized conservative well hasFor a symmetric well the force is odd, the cubic moment vanishes, and the displacement remainder improves to $O(v_T^{-4})$. For an arbitrary compact force without equal outside potential levels, its signed area instead gives the leading $v_T^{-1}$ displacement.

For circular motion, let $C=2\pi R$ and use the local arclength model, with nonoverlapping support $w<C$. In each encounter the relative position traverses $w$, while the particle advances $\Delta x$. Outside the support the particle is stationary and the trap must traverse the remaining relative arclength $C-w$. Hence the time between kicks and the repeated kicks from a circular optical trap response areThe last equality uses the single-passage identity. Equivalently, kicks arrive at the relative lap frequency $f_T-f_p$.

The precise small-correction condition for the large-speed approximation for circular-trap kicks is $\Delta x/C\ll 1$. A sufficient high-speed regime for a localized well is $v_T\gg \Vert F{\Vert }_{\infty }/\zeta$, with $w<C$ fixed: writing $ϵ=\Vert u{\Vert }_{\infty }/v_T<1$, the positive-displacement integral bounds $\Delta x/C\le (w/C){ϵ}^2/(1-ϵ)$. ThusThe locally straight force profile on a real circular track also presumes a trap width small relative to the radius; $R\gg a$ supplies the particle-size part of that approximation.

For the triangular optical-trap response, put $v_c=F/\zeta =C{f}_c$, $\beta =v_T/v_c$ and $\alpha =w/C=X_0/(\pi R)$. In the passing regime $\beta >1$, the two constant-force halves giveSubstitution into the circular formula yieldsFor $0<\beta \le 1$, the particle remains locked to the trap and $f_p/f_c=\beta$. At the ideal triangular cusp this is understood as sticking, or as the limit of a rounded well; it is not a finite passing kick. The passing branch tends continuously to one at $\beta \downarrow 1$, while at large $\beta$ it behaves as $\alpha /\beta$. The zero-speed limiting response is zero. These branches apply to separated kicks, $0<\alpha <1$.

A Riemannian metric is a smooth field of positive-definite symmetric bilinear forms $g_x$ on the tangent spaces, equivalently a smooth section of ${\mathrm{Sym}}^2{T}^{\ast }M$ with that positivity. Its musical isomorphism isNondegeneracy makes every fiber map an isomorphism. In coordinates the map has matrix $g_{ij}$, and its inverse ${♯}_g$ has matrix $g^{ij}$. Smoothness of the inverse follows from the inverse-matrix formula and the nonvanishing determinant. Both maps cover the identity on $M$, giving a smooth vector bundle isomorphism.

For local flattening of a Riemannian metric, choose a manifold chart neighborhood $V$ around $p$ with closure contained in $U$, and let $h={\varphi }^{\ast }\delta$ on $V$. Choose a smooth bump function $0\le \chi \le 1$, with compact support in $V$ and equal to one on a smaller neighborhood $W$ of $p$. SetThe formulas glue smoothly because the modification is supported strictly inside $V$. A convex combination of positive-definite forms is positive definite. On $W$ the metric is exactly $h$, so $\varphi$ gives an isometry to a Euclidean open subset. Outside $U$ the original metric remains unchanged.

Geodesic completeness means that every geodesic with arbitrary initial point and tangent vector extends for every real affine parameter. On each connected component, the Hopf-Rinow theorem equates this with completeness of the Riemannian distance.

The Euclidean-end claim needs a compact-core interpretation that is absent from its literal hypotheses. Indeed, $M=\{x\in {\mathbb{R}}^n:|x|>1\}$ with $n\ge 2$, $X=\varnothing$ and $g=\delta$ satisfy those hypotheses. The geodesic $\gamma (t)=(2-t)e_1$ reaches the missing unit sphere at $t=1$ and cannot continue within $M$. Thus the printed assumptions alone do not prove completeness.

Here is the intended compact-core completeness for a Euclidean end. Write $\Phi :M\setminus X\to \{|x|>1\}$ for the end coordinates and assume additionally thatis compact for every sufficiently large $R$. Choose such an $R$ beyond the metric transition and large enough to include the initial point of a geodesic. Its speed $s$ is constant. On every segment outside $K_R$, the Euclidean radius changes at a rate at most $s$, so over a finite parameter interval it cannot exceed $R+sT$. The full segment therefore stays in the compact set $K_{R+sT}$. Its velocities also stay in a compact subset of $TM$, since their metric norm is fixed and the base set is compact. The smooth geodesic ordinary differential equation consequently extends past any supposed finite endpoint. The reversed-time argument is identical. This proves completeness under the compact-core condition, and explains exactly what the exterior-of-a-ball counterexample lacks.

For upward stability of Riemannian completeness, lengths of all curves satisfy $L_{\tilde{g}}\ge L_g$, hence $d_{\tilde{g}}\ge d_g$. A $d_{\tilde{g}}$-Cauchy sequence is therefore $d_g$-Cauchy and has a limit $p$ because $g$ is complete. Smooth positive-definite metrics induce the manifold topology. More explicitly, on a small coordinate ball around $p$, $\tilde{g}$ has a bounded largest matrix eigenvalue, so the coordinate straight segment gives $d_{\tilde{g}}(p,q)\le C|\varphi (q)-\varphi (p)|$. Thus convergence to $p$ is also in $d_{\tilde{g}}$. That distance is complete, and Hopf-Rinow yieldsFor disconnected $M$, apply this argument on each component.

Use $G=c=1$, select the retarded source field with no added homogeneous radiation, and work to leading order in the weak-field approximation and long-wavelength source approximation. The retarded fundamental solution of the Linearized Einstein equations in Lorenz gauge in linearized gravity isFor $r\gg d$, the denominator is $r$ to leading order. The delay is $t-r+\mathbf{n}\cdot {\mathbf{x}}^{\prime }+O(d^2/r)$, where $\mathbf{n}=\mathbf{x}/r$. In the usual slowly evolving source regime, its characteristic time $\tau$ obeys $d/\tau \ll 1$, so the source-size part of the delay is negligible:The stress-energy conservation equation ${\partial }_{\mu }{T}^{\mu \nu }=0$ is required at the approximation order used; it follows from the divergence of the gauge-fixed field equation. Compact support removes the integration by parts surface terms. Since $T^{00}=T_{00}$ and $T^{ij}=T_{ij}$ in this signature,Combining these identities gives the retarded quadrupole fieldFor an ordinary nonrelativistic bound source, $\tau \sim d/v$ makes $d/\tau \sim v\ll 1$. More generally the size relative to the variation timescale must also be small; speed alone does not exclude a rapidly varying small-amplitude motion.

The radiation boundary condition for linearized gravity is essential for a statement about the full field. Without it, the displayed implication is false: take $T_{\mu \nu }=0$ and add ${\bar{h}}_{11}=A\cos (k(t-z))$, ${\bar{h}}_{22}=-{\bar{h}}_{11}$, all other components zero. This weak plane gravitational wave satisfies both the homogeneous wave equation and Lorenz gauge in linearized gravity, while every $I_{ij}$ is zero. The proved formula is for the retarded source contribution, with the standard slow-source qualification. Its transverse–traceless projection gives the physical radiative strain.

The space of test functions is $\mathcal{D}(\mathbb{R})=C_c^{\infty }(\mathbb{R})$: its elements are smooth functions with compact support. For a fixed compact set $K$, put ${\mathcal{D}}_K=\{\phi \in C^{\infty }(\mathbb{R}):\mathrm{supp}\phi \subset K\}$ and use the seminorms $p_m(\phi )={\max }_{0\le j\le m}\Vert {\phi }^{(j)}{\Vert }_{\infty }$. The space of test functions carries the usual test-function inductive limit topology of these spaces. In particular, a sequence converges precisely when its supports eventually lie in one compact set and every derivative converges uniformly. Thus ${\phi }_j\to \phi$ means that a common $K$ contains their supports and $p_m({\phi }_j-\phi )\to 0$ for every $m$.

A distribution is a continuous linear functional on this space of test functions, and ${\mathcal{D}}^{\prime }(\mathbb{R})$ denotes their space. We use complex-linear, bilinear pairings $⟨u,\phi ⟩$. Equivalently, for every compact set $K$ there are $C_K$ and a nonnegative integer $m_K$ such thatThe usual weak convergence of distributions is $u_j\to u$ if $⟨u_j,\phi ⟩\to ⟨u,\phi ⟩$ for every test function. This specifies the convergence used below; it does not require convergence in any norm.

The distributional derivative is defined byThe derivative map sends ${\mathcal{D}}_K$ continuously to itself, with $p_m({\phi }^{\prime })\le p_{m+1}(\phi )$. The preceding continuity estimate therefore proves that $u^{\prime }$ is again a distribution. This definition extends the ordinary derivative of a smooth function, by integration by parts.

Choose the translation of a distribution convention ${\tau }_{h}f(x)=f(x-h)$. Its action on a test function isFor fixed $h$, the translated test functions have translated compact support and unchanged derivative sup norms, so this defines a distribution. For the differentiability of distribution translations, apply the Taylor theorem in integral form:All supports lie in one slightly enlarged compact set for $|h|\le 1$, and the same identity for every derivative proves uniform convergence. The continuity of $u$ now givesThe minus sign in the translation parameter is necessary for this convention.

For (a), choose a test function $\eta$ with $\int \eta =1$. If $\int \psi =0$, then $F(x)={\int }_{-\infty }^x\psi (s)ds$ is a test function, since the zero integral makes it vanish beyond both ends of the compact support. If $u_1^{\prime }=0$, then $⟨u_1,\psi ⟩=⟨u_1,F^{\prime }⟩=0$. Decompose $\phi =(\int \phi )\eta +\psi$ to obtainConversely, these constant regular distributions have zero distributional derivative. This also proves the general fact that a distribution with zero derivative is constant.

For (b), choose a cutoff function $\rho$ equal to one near zero. Every test function decomposes asThe quotient extends as a smooth function at zero, and it has compact support. If $x{u}_2=0$, the definition of multiplication of a distribution by a smooth function gives $⟨u_2,\phi ⟩=⟨u_2,\rho ⟩\phi (0)$. Conversely, $x{\delta }_0=0$. Thus the kernel of multiplication by a coordinate is exactlyIn particular, derivatives of the Dirac delta distribution are not additional solutions: $x{\delta }_0^{\prime }=-{\delta }_0$.

Now write $D=d/dx$ and $P(z)=-z^n+{\sum }_{j=0}^{n-1}{a}_j{z}^j$, where $n\ge 1$. The characteristic roots of a constant-coefficient differential equation are the distinct roots of a polynomial ${\lambda }_1,\dots ,{\lambda }_r$ of $P$, with multiplicities $m_1,\dots ,m_r$. The distributional regularity of a constant-coefficient ordinary differential equation can be proved without assuming regularity in advance. Put $q_j(z)=(z-{\lambda }_j{)}^{m_j}$ and $Q={\prod }_j{q}_j=-P$. Since the $q_j$ are coprime polynomials, polynomial division and Bezout identity give polynomials $b_j$ such thatFor example, invert $Q/q_j$ modulo $q_j$, sum the resulting expressions, and absorb a remaining multiple of $Q$ into one coefficient. For $P(D)v=0$, the kernel decomposition for coprime polynomials therefore givesThe Leibniz rule for multiplication of a distribution by a smooth function implies $D^{m_j}(e^{-{\lambda }_{j}x}{v}_j)=0$. Repeatedly using a distribution with zero derivative is constant shows that a distribution with $m$th derivative zero is a polynomial of degree at most $m-1$: subtract the polynomial primitive of its constant $(m-1)$st derivative, and induct on $m$. Consequently the most general exponential polynomial solution of a constant-coefficient differential equation isEvery displayed term is annihilated by $P(D)$, so all coefficients are allowed. For the linear independence of these $n$ functions, on a relation, apply ${\prod }_{i\ne j}(D-{\lambda }_i{)}^{m_i}$; on $e^{{\lambda }_{j}x}$ times a polynomial of degree less than $m_j$, each remaining factor acts invertibly on that polynomial space. Hence the $j$th polynomial must vanish. Every distributional solution is an analytic function, and in particular a classical solution. For real coefficients and real-valued distributions, take complex conjugate coefficients at conjugate characteristic roots of a constant-coefficient differential equation, or equivalently use the corresponding real sine and cosine forms.

For the last equation, the operator is $xP(D)$: the coordinate multiplies the result of differentiation. The kernel of multiplication by a coordinate says precisely thatChoose the retarded fundamental solution of a constant-coefficient ordinary differential operator $E=Hw$, where $H$ is the Heaviside function and $w$ is the analytic solution of $P(D)w=0$ withSuch $w$ exists uniquely by the elementary initial value problem for a constant-coefficient ordinary differential equation; alternatively the residue construction in the next solution gives it explicitly. The distributional jump formula for a Heaviside product isIt follows by induction from $D(Hw)=H{w}^{\prime }+w(0){\delta }_0$, using integration by parts. With the chosen initial derivatives, every lower-order jump term vanishes and the leading coefficient $-1$ gives $P(D)E={\delta }_0$. Subtracting $CE$ reduces the last equation to the homogeneous one. Thus the coordinate-degenerate constant-coefficient differential equation has exactly the solutionsThere are $n+1$ independent constants. For $n\ge 2$, derivatives through order $n-2$ match at zero, while the derivative of order $n-1$ may jump; for $n=1$, the function itself may jump. An arbitrary distribution concentrated at zero cannot be added, because its image under the nonzero leading derivative would contain a nonvanishing highest derivative of the Dirac delta distribution.

An oscillatory integral defines a distribution by cancellation, even when its oscillatory integral amplitude is not integrable in frequency. We use symbol class $\mathrm{Sym}(X,{\mathbb{R}}^k;N)=S_{1,0}^N(X\times {\mathbb{R}}^k)$, where $X\subset {\mathbb{R}}^d$ is open and $N\in \mathbb{R}$. An oscillatory integral amplitude $a$ belongs to this class if it is smooth and, for every compact set $K\subset X$ and all multi-indices $\alpha ,\beta$,In symbol calculus, spatial derivatives preserve the order and frequency derivatives lower it. The fixed number $N$ is independent of $K$; this will produce a finite global order of a distribution, although continuity constants may depend on $K$.

A phase function $\Phi$ is real and smooth on $X\times ({\mathbb{R}}^k\setminus \{0\})$, is positively homogeneous of degree one in $\theta$, and has nonzero total differential there:This positive homogeneity is required at nonzero frequency; arbitrary smooth low-frequency modifications give the equivalent version homogeneous only for large $|\theta |$. In particular, smoothness at zero is not an additional requirement on a general homogeneous phase function. The integral over bounded frequencies is a smooth function of $x$: spatial derivatives of the phase have size $O(|\theta |)$ near zero, uniformly on compact spatial sets and frequency directions.

Choose a cutoff function $\chi \in C_c^{\infty }({\mathbb{R}}^k)$ equal to one near zero. The proposed meaning of the oscillatory integral distribution isExistence and cutoff independence of an oscillatory integral require a proof. Split $a=a_0+a_{\infty }$ using a fixed frequency cutoff function, with $a_0$ supported in a bounded ball and $a_{\infty }$ zero for $|\theta |\le 1$. The low-frequency part is already absolutely integrable. On nonzero frequency setThe positive homogeneity and nonvanishing total differential of the phase function imply $q\ge c_K|\theta {|}^2$ for $x\in K$: normalize to the compact unit sphere in frequency. Direct differentiation gives $L{e}^{i\Phi }=e^{i\Phi }$.

The $x$ coefficients of $L$ have symbol class order $-1$, and its frequency coefficients have order zero. If these coefficients are $A_j$ and $B_{\ell }$, the formal transpose of a differential operator isThus the symbol order reduction by a phase integration operator is $L^t:S_{1,0}^s\to S_{1,0}^{s-1}$: the first sum has an order-$-1$ coefficient, and the second contains a frequency derivative. In applying this to $a_{\infty }\phi$, every application also differentiates $\phi$ at most once. Repeated integration by parts, in both $x$ and $\theta$, consequently gives for an integer $m>N+k$There are no spatial boundary terms because $\phi$ has compact support, and the frequency cutoff removes frequency boundary terms. Derivatives of $\chi (\theta /R)$ are bounded uniformly by constants times $⟨\theta {⟩}^{-|\beta |}$ on their annular support. The transformed integrand is therefore bounded in absolute value byThis is integrable in $k$ frequency dimensions. Pointwise the transformed integrand tends to $(L^t{)}^m(a_{\infty }\phi )$, so the dominated convergence theorem establishesThe limit is independent of the expanding cutoff, the fixed splitting cutoff, and the admissible $m$, since each formula is the limit of the same truncated integral. Derivatives landing on the expanding cutoff can also be estimated directly by $O(R^{N-m+k})$, which tends to zero. In particular,This proves continuity on the space of test functions and the finite order of an oscillatory integral distribution. The same $m$ works for all compact sets. When $N<-k$, $m=0$ is possible and the original frequency integral is absolutely integrable; cancellation is needed for general $N$.

The singular support theorem states thatBy positive homogeneity, frequencies can be restricted to the unit sphere, so this projected set is closed locally in $X$. More precisely, the stationary-direction bound for singular support restricts the frequency directions to the closed conic support of an oscillatory amplitude; using a closed directional support avoids losing limits occurring at arbitrarily high frequency. The simpler displayed bound is sufficient here.

To prove the bound, take $x_0$ outside the displayed stationary set. On a sufficiently small spatial neighborhood and all unit frequency directions, $|{\nabla }_{\theta }\Phi |$ is bounded below. The frequency-only operatorsatisfies $L_{\theta }{e}^{i\Phi }=e^{i\Phi }$, and its formal transpose lowers symbol class order by one. A spatial derivative of order $|\alpha |$ of $e^{i\Phi }a$ has amplitude order at most $N+|\alpha |$. Applying the frequency integration by parts more than $N+|\alpha |+k$ times makes that derivative absolutely integrable, uniformly on smaller compact sets. Every spatial derivative therefore exists and is continuous there. The oscillatory integral distribution is a smooth function near $x_0$, proving the singular support assertion.

For the linear transport equation, use spacetime $X={\mathbb{R}}^n\times (0,\infty )$, frequency $\theta \in {\mathbb{R}}^n$, andThis is a valid phase function, because ${\nabla }_x\Phi =\theta \ne 0$ at nonzero frequency; the oscillatory integral amplitude is in symbol class order zero. The Fourier representation of the Dirac delta function gives the transport of a Dirac point mass:Its rigorous spacetime distribution pairing is $⟨u,\phi ⟩={\int }_0^{\infty }\phi (ct,t)dt$. Consequentlywhere the endpoints vanish since a spacetime test function has compact support in $t>0$. At each fixed $t$, $⟨u(t),\psi ⟩=\psi (ct)$, so weak convergence of distributions gives $u(t)\to {\delta }_0$ as $t\downarrow 0$. This is the required initial trace.

Finally, ${\nabla }_{\theta }\Phi =x-ct$, so the singular support theorem confines singularities to $x=ct$. In fact equality holds: $u$ is a nonzero order-zero distribution on that trajectory and vanishes off it. A smooth function supported on this set of empty interior must vanish, so $u$ cannot be smooth in any neighborhood of a point of the trajectory. ThusThe method of characteristics also proves uniqueness among solutions with a distributional initial trace: the change of variables $y=x-ct$ turns the linear transport equation into ${\partial }_{t}w=0$. Such a distribution is constant in $t$; pairing with spatial test functions reduces this assertion to a distribution with zero derivative is constant. Its initial trace fixes $w={\delta }_0$, so the moving Dirac delta distribution above is the unique solution. There is transport of the singularity along the characteristic and no smoothing.

The shear $(t,x,v)\mapsto (t,x+tv,v)$ is a smooth diffeomorphism with determinant one and sends compact sets to compact sets. Thus $F,G$ are locally integrable functions. On a compact time interval $I$ and compact label set $K\subset {\mathbb{R}}^6$, the integral ${\int }_K{\int }_I|G|dtdxdv$ is finite. The Fubini theorem, followed by a countable exhaustion of time intervals and label sets, shows that $G(\cdot ,x,v)\in L_{\mathrm{loc}}^1(\mathbb{R})$ for almost every $(x,v)$.

To transform the weak equation, choose $\varphi (t,x,v)=\rho (x-tv,v)\psi (t)$, where $\rho \in C_c^{\infty }({\mathbb{R}}^6)$ and $\psi \in C_c^{\infty }(\mathbb{R})$. Including a velocity cutoff in $\rho$ ensures the required compact support. The transport derivative of this test function is $\rho (x-tv,v){\psi }^{\prime }(t)$. Change variables to the characteristic labels to obtainThe fundamental test-function identity makes the bracket zero almost everywhere. Choose a countable dense family of time test functions on each bounded interval and extend by continuity; there is then a single negligible label set outside which ${\partial }_{t}F=G$ as one-dimensional distributions.

For each remaining label, subtract the primitive $H(t)={\int }_0^{t}G(s)ds$. The distributional derivative of $F-H$ is zero, so it is a constant almost everywhere. Choose the resulting absolutely continuous representative along free characteristics, $\tilde{F}(t)=C+H(t)$. The fundamental theorem of calculus givesThe representative qualification is essential: an arbitrary locally integrable representative can be modified on one time slice without changing the distributional equation, and would not satisfy the identity for every pair of times.

Interpret $L$ as the velocity operator appearing in the equation and in the stated quadratic-form hypothesis. Let $E(t)=\int \int |f(t,x,v){|}^{2}dxdv$. Smoothness turns the distributional equation into the usual equation. Multiply it by $f$ and integrate; for complex functions use the real part of the inner product. Spatial integration by parts gives $\int v\cdot {\nabla }_x(f^2)=0$, while the coercivity hypothesis at each $x$ givesThus $e^{2\delta t}E(t)$ is nonincreasing. The coercive energy estimate for kinetic transport isNo division by $E$ is required, so the zero solution causes no exception. Compact spatial support, or sufficient decay to eliminate the boundary flux, is enough; compact support in time is unnecessary.

Strict positivity throughout phase space is incompatible with compact support. The following calculation uses positive smooth rapidly decaying solutions with all displayed integrals justified; nonnegative cases use the corresponding entropy limits where justified.

Write $B=|{\mathbb{S}}^2|$ and use integration over $(v,v_{\ast },\sigma )$. The two symmetrized weak collision identities areandThey follow from particle interchange and the angular exchange for elastic collisions. They are the weak identities needed here; no invalid fixed-angle Jacobian is used.

The collision invariants $1$, $v$ and $|v{|}^2$ obey ${\varphi }^{\prime }+{\varphi }_{\ast }^{\prime }=\varphi +{\varphi }_{\ast }$. Momentum follows from $v^{\prime }+v_{\ast }^{\prime }=v+v_{\ast }$, and energy from $|v^{\prime }{|}^2+|v_{\ast }^{\prime }{|}^2=|v{|}^2+|v_{\ast }{|}^2$. HenceFor the Boltzmann H functional $H=\int \int f\log fdxdv$, the spatial transport contribution is a boundary divergence and the term $\int Q$ is zero. Set $a=f^{\prime }{f}_{\ast }^{\prime }$, $b=f{f}_{\ast }>0$. The second weak identity with $\varphi =\log f$ gives the Boltzmann H theoremIndeed this integrand with its negative sign is $b(1-X)\log X$ for $X=a/b>0$, which is nonpositive. The kinetic H functional decreases; the physical entropy with opposite sign increases.

Use a bounded BV extension operator for the bounded Lipschitz domain. Extend $u_k$ to $w_k\in BV({\mathbb{R}}^n)$, supported in one fixed compact set, with $\Vert w_k{\Vert }_{BV({\mathbb{R}}^n)}\le C_{\Omega }\Vert u_k{\Vert }_{BV(\Omega )}$. Such an extension is obtained by reflection in Lipschitz boundary charts, a partition of unity and a fixed cutoff function. A naive zero extension must also charge its boundary-trace jump; it cannot discard that contribution.

Let $M$ bound these norms and let $w_{k,ϵ}=w_k\ast {\rho }_{ϵ}$ be their mollifications. The BV mollification error estimate givesThe variation measure here is on ${\mathbb{R}}^n$. The printed lemma's $|Dw|(\Omega )$ needs this correction unless it also assumes all the derivative mass lies in $\Omega$: a nonconstant bump supported outside $\overline{\Omega }$ disproves its literal wording. The correct estimate follows by averaging the BV translation estimate $\Vert w(\cdot -h)-w{\Vert }_1\le |h||Dw|({\mathbb{R}}^n)$ over a mollifier supported in $|h|\le ϵ$.

For each fixed $ϵ>0$, convolution gives uniform boundsThe mollifications have common compact support and are uniformly equicontinuous. The Arzelà-Ascoli theorem supplies a uniformly convergent subsequence at each scale $ϵ=1/m$. A diagonal subsequence converges at every one of those scales. For two late members of that subsequence,First choose large $m$, then late $k,\ell$; the sequence is a Cauchy sequence in $L^1$. Its limit $w$ belongs to the BV space because the total variation seminorm on a domain is sequentially lower semicontinuous. Uniform derivative bounds and integration by parts give $D{w}_k\stackrel{\ast }{\rightharpoonup }Dw$. Restriction to $\Omega$ proves the required weak-star convergence in BV.

The structure theorem for functions of bounded variation is the measure representation of the distributional derivative: for $u\in BV(\Omega )$ there is a unique finite vector Radon measure $Du$ withfor every $\phi \in C_c^1(\Omega ;{\mathbb{R}}^n)$. It also has a polar decomposition of a vector measure $Du={\sigma }_u|Du|$ with $|{\sigma }_u|=1$ for $|Du|$-almost every point.

To prove it, let $L(\phi )=-\int u\mathrm{div}\phi$. The test-function definition of total variation seminorm, applied to both signs, gives $|L(\phi )|\le \mathrm{TV}(u)\Vert \phi {\Vert }_{\infty }$. Compactly supported smooth vector fields are uniformly dense in $C_0(\Omega ;{\mathbb{R}}^n)$, so $L$ extends uniquely to a bounded functional there. The Riesz-Markov-Kakutani representation theorem, applied componentwise, gives the unique vector measure $Du$. The operator norm of $L$ is exactly the defining variation supremum; the vector-measure dual norm is $|Du|(\Omega )$, proving equality. Conversely any finite measure satisfying the identity bounds that supremum, so this also characterizes membership in the BV space.

The Radon-Nikodym theorem applied to the components relative to $|Du|$ gives ${\sigma }_u$; the definition of the variation measure forces $|{\sigma }_u|=1$ almost everywhere. If $u$ is smooth, ordinary integration by parts gives $Du=\nabla u{\mathcal{L}}^n$. Compact support of the test field removes any boundary contribution; no boundary regularity is needed for this representation statement.

Use complex-linear distribution pairings, without conjugation. Write $\mathcal{D}=C_c^{\infty }({\mathbb{R}}^n)$ for the space of test functions. The smoothing convolution with a test function isOn a compact set of $x$ values, all the translated test functions have support in one compact set. Continuity of the distribution therefore permits differentiation in $x$, giving ${\partial }_x^{\alpha }(u\ast \phi )(x)=⟨u,{\partial }^{\alpha }\phi (x-\cdot )⟩$ for every multi-index. In particular this convolution is a smooth function, even if $u$ is not tempered.

For the first associativity identity, integration against $\psi$ and the distribution pairing can be interchanged: the integrand has a common compact support in $y$, depends smoothly on the integration variable, and satisfies the finite-order continuity estimate there. ConsequentlyThe common-support argument matters: a general distribution cannot be paired with arbitrary noncompact functions.