\documentclass[10pt,letterpaper,oneside]{amsart} \usepackage{geometry} \geometry{letterpaper,textwidth=6.15in,textheight=8.65in,centering} \usepackage[T1]{fontenc} \usepackage[utf8]{inputenc} % Standard LaTeX text and mathematics fonts. \usepackage{microtype} \setlength{\emergencystretch}{2em} \usepackage[british]{babel} \usepackage{amssymb,amsmath,amsfonts,amsthm,mathrsfs,mathtools} \usepackage[dvipsnames]{xcolor} \usepackage{tikz} \usetikzlibrary{calc,arrows.meta} \usepackage{graphicx} \input{figures/colors.tex} % figure colours (sources: visual-final/tikz-ka) \renewcommand\figdata{figures/data/} \usepackage{enumitem} \usepackage{needspace} \definecolor{paperlink}{RGB}{0,0,114} \usepackage[colorlinks=true,linktoc=all,hypertexnames=false,linkcolor=paperlink,citecolor=paperlink,urlcolor=paperlink]{hyperref} \hypersetup{pdftitle={From D\"oblin to Aleksandrov},pdfauthor={Am\'elie Loher, Connor Mooney and Cl\'ement Mouhot}} \usepackage{parskip} \newcommand{\R}{\mathbb R} \newcommand{\1}{\mathbf 1} \newcommand{\dd}{\,\mathrm d} \newcommand{\abs}[1]{\lvert#1\rvert} \newcommand{\norm}[1]{\lVert#1\rVert} \newcommand{\Sym}{\operatorname{Sym}} \newcommand{\supp}{\operatorname{supp}} \newcommand{\sgn}{\operatorname{sgn}} \newcommand{\dist}{\operatorname{dist}} \newcommand{\TV}{\mathrm{TV}} \newcommand{\PA}{\mathcal L_A} % physical kinetic operator, matrix diffusion \newcommand{\Pa}{\mathcal L_a} % physical kinetic operator, scalar diffusion \newcommand{\KA}{\mathcal K_A} % reflected terminal operator, matrix \newcommand{\Ka}{\mathcal K_a} % reflected terminal operator, scalar \newcommand{\Lop}{\mathcal K} % parabolic operator with a transported variable \newcommand{\kin}{k} % K denotes the kinetic measure kernel; P denotes its velocity marginal. \newcommand{\Kf}[1]{\widehat K^{#1}} % Fourier kernel at frequency #1 \newtheorem{theorem}{Theorem}[section] \newtheorem{proposition}[theorem]{Proposition} \newtheorem{lemma}[theorem]{Lemma} \newtheorem{corollary}[theorem]{Corollary} \theoremstyle{remark} \newtheorem{remark}[theorem]{Remark} \numberwithin{equation}{section} \allowdisplaybreaks \setcounter{tocdepth}{1} \setlist{topsep=0.4em,itemsep=0.2em,parsep=0pt} \setlist[enumerate,1]{label=\textup{(\roman*)},ref=\textup{(\roman*)}} \title[From D\"oblin to Aleksandrov]{From D\"oblin to Aleksandrov} \author{Am\'elie Loher} \address{All Souls College, University of Oxford} \author{Connor Mooney} \address{Department of Mathematics, UC Irvine} \author{Cl\'ement Mouhot} \address{DPMMS, University of Cambridge} \date{2 October 2026} \begin{document} \begin{abstract} We prove Aleksandrov maximum principles for linear second-order kinetic equations in non-divergence form with measurable, uniformly elliptic diffusion coefficients. For coefficients $A(t,v)$ we obtain every source exponent $p>2d+1$. For autonomous coefficients $a(x,v)$ in one space dimension, we obtain every $p>p_*(\frac\Lambda\lambda)$, where $p_*(\frac\Lambda\lambda)<4$ depends only on the ellipticity ratio. The first step in both proofs is Döblin--Fourier cancellation: two positive parts of the kernel have nearly opposite phases, and a common lower bound for their velocity marginals allows us to cancel fixed portions of them. For the autonomous problem, this argument applies on velocity intervals away from zero. A preliminary Aleksandrov estimate gives H\"older continuity and then a Harnack comparison. A homogeneous Bellman barrier controls concentration in position and velocity; together with Harnack, it makes the local density estimates summable across zero velocity, including at $p=4$. We give the full proofs and the passage from the maximum principle to local H\"older continuity. \end{abstract} \maketitle \tableofcontents \section{Introduction} Diffusion changes velocity; transport turns those changes into motion in position. This interaction is the source of regularity for the kinetic operator $\partial_t+v\cdot\nabla_x-A:D_v^2$. Kolmogorov computed its fundamental solution for $A=I_d$~\cite{Kolmogorov1934}. For smooth coefficients, H\"ormander's theorem explains why diffusion in velocity also smooths position: commutators with transport supply the missing derivatives~\cite[Theorem~1.1]{Hormander1967}. Here the coefficient is only measurable. We ask which estimates follow from ellipticity alone, without differentiating the coefficient. For a uniformly elliptic operator $a^{ij}\partial_{ij}$ on a bounded domain $\Omega\subset\R^N$ the Aleksandrov--Bakelman--Pucci estimate~\cite{Aleksandrov1960,Bakelman1961,Pucci1966} answers this question in the form \refstepcounter{equation} \[ \sup_\Omega u\leq\sup_{\partial\Omega}u_++C\norm{f_+}_{L^N(\Omega)}\,, \tag{\theequation}\label{e:abp-elliptic} \] whenever $a^{ij}\partial_{ij}u\geq-f$; see~\cite[Chapter~9]{GilbargTrudinger}. The constant depends only on the dimension, ellipticity and domain. Krylov's parabolic estimate uses $L^{N+1}$ in space and time, with a factor $r^{\frac{N}{N+1}}$ on a cylinder of spatial radius $r$~\cite[Theorem~3]{Krylov1976}; see also~\cite{Tso1985,FerrettiSafonov}. These bounds control the effect of a source concentrated on a small set. That control is a starting point for the Krylov--Safonov Harnack theory~\cite{KrylovSafonov}. The classical proof does not carry over directly to the kinetic operator. It measures the source where the graph of $u$ is touched from above by an affine function, using the determinant of the Hessian to estimate the volume of the resulting slopes. The kinetic equation controls $D_v^2u$ but contains no second derivative in position, so this argument leaves the position directions uncontrolled. We instead use transport to turn diffusion in velocity into cancellation of spatial waves. This gives two estimates of the form~\eqref{e:abp-elliptic}, with the power of the radius fixed by kinetic scaling. In each coefficient class, uniform ellipticity and Borel measurability suffice: no continuity or derivative bound is assumed. The solution is continuous on the closed cylinder, of class $C^{1,1,2}$ inside it, with no integrability assumption on its derivatives at the boundary. \subsection{Main results} Let $d\geq1$, $P_0=(t_0,x_0,v_0)\in\R\times\R^d\times\R^d$ and $R>0$. The kinetic cylinder is \refstepcounter{equation} \[ Q_R^-(P_0) = \bigl\{(t,x,v):\ t_0-R^20$; on those two sets the estimates below prescribe nothing. Figure~\ref{fig:ka-boundary} shows these faces, and their reflection used in Section~\ref{s:maximum}. Let $A:\R\times\R^d\to\Sym_d$ be Borel measurable and let $0<\lambda\leq\Lambda<\infty$ with \refstepcounter{equation} \[ \lambda I_d\leq A(t,v)\leq\Lambda I_d \quad\text{for almost every }(t,v)\in\R\times\R^d\,. \tag{\theequation}\label{e:A-ellipticity} \] For $u\in C^{1,1,2}$ (see the notation fixed at the end of this introduction) set \refstepcounter{equation} \[ \PA u=\partial_tu+v\cdot\nabla_xu-A(t,v):D_v^2u\,. \tag{\theequation}\label{e:PA} \] \begin{theorem}[Coefficients depending on time and velocity] \label{t:A} Let $A$ satisfy~\eqref{e:A-ellipticity}, let $\PA$ be the operator~\eqref{e:PA}, let $2d+10\}}}_{L^p(Q_R^-(P_0))}$. \end{corollary} The powers of $R$ follow from kinetic scaling. The map $(t,x,v)\mapsto(R^2t,R^3x,Rv)$ carries $Q_1^-(0)$ onto $Q_R^-(0)$ and multiplies Lebesgue measure by $R^{4d+2}$, the homogeneous dimension being $2+3d+d$; a source rescaled from $Q_R^-$ to $Q_1^-$ acquires a factor $R^2$ and its $L^p$ norm a factor $R^{-\frac{4d+2}{p}}$. The resulting power of $R$ vanishes at $p=2d+1$, where the estimate would be scale invariant; Theorem~\ref{t:A} covers every larger exponent. The strong $L^{2d+1}$ endpoint fails even for constant coefficients, by the nested-bump construction in~\cite[Section~8.5]{DoeblinFourier}. For $d=1$ the scaling threshold is $p=3$. Theorem~\ref{t:a} reaches every $p>p_*(\frac\Lambda\lambda)$ with $p_*<4$, and in particular includes $p=4$. The improvement below four depends on the ellipticity ratio. The threshold $p_*$ is the one supplied by the Bellman argument of Section~\ref{s:below-four}; neither its optimality nor an estimate at $p=p_*$ is claimed. Appendix~\ref{s:counterexample} proves a stronger obstruction for autonomous coefficients. For every $p<4d$, a suitable fixed ellipticity ratio admits smooth examples whose height stays positive while their source norm tends to zero. The construction starts with a positive stationary profile $H$ of degree $0<\alpha<1$, singular only at the origin. Flattening $1-H$ there keeps its height and confines its source to a small kinetic annulus. A decaying velocity barrier removes the initial data without increasing the positive source. The source norm is bounded by $Cr^{\alpha-2+\frac{4d}{p}}$, which tends to zero when $\alpha$ is sufficiently close to one. The quantifiers matter: the coefficient bounds can depend on $p$, but stay fixed along each sequence. In dimension one, this rules out any threshold below four valid for all ellipticity ratios. It is consistent with Theorem~\ref{t:a}, whose improvement is for each fixed ratio and may shrink as the ratio grows. In higher dimensions the construction concerns $A(x,v)$, not the restricted class $A(t,v)$ of Theorem~\ref{t:A}. For smooth $A(t,v)$, the argument also gives the Lorentz-space endpoint with source in $L^{2d+1,1}$; Remark~\ref{r:lorentz} records the short additional argument. \subsection{Strategy} Both proofs start from the same Green-density estimate. We reverse time to use terminal-value kernels and reflect position to preserve the sign of transport. This changes neither the estimates nor their exponents; we return to the original variables at the end. Each operator then becomes an instance of \refstepcounter{equation} \[ \Lop=\partial_\sigma+B(\sigma,v):D_v^2+b(v)\cdot\nabla_z \qquad\text{on }\R_\sigma\times D_v\times\R^d_z\,, \tag{\theequation}\label{e:intro-Lop} \] Diffusion acts in $v$, transport acts in $z$, and the coefficients are independent of $z$. For $A=A(t,v)$, we take $z$ to be the reflected position, $b(v)=v$ and $D=\R^d$. For $a=a(x,v)$ in one dimension, we work on a velocity interval separated from zero. Position then moves in one direction and can serve as time. After normalisation, $D$ is a fixed bounded interval. The old time becomes the transported variable; autonomy is what makes the new diffusion coefficient independent of it. Fourier transformation in $z$ separates the waves. The transformed kernel $\Kf\xi$ is a complex measure in terminal velocity. Its variation is bounded by the positive velocity marginal $P$, obtained at frequency zero. With $v$ denoting the starting velocity, our key estimate is \refstepcounter{equation} \[ \norm{\Kf\xi_{\sigma,\tau}(v,\cdot)}_{\TV}\leq C\exp\bigl(-c\,(\tau-\sigma)\abs\xi^{\frac23}\bigr)\,, \tag{\theequation}\label{e:intro-decay} \] To obtain this decay, we choose two curves in velocity space whose accumulated displacements in $z$ give phases differing by $\pi$. The curves end at the same centre, and a barrier retains a fixed amount of the kernel in a small tube around each curve. The Lipschitz bound on $b$ keeps the phases close to those of the curves. The resulting terminal velocity distributions need not agree. We continue both in one fixed ball, where the parabolic Harnack inequality~\cite{KrylovSafonov,FerrettiSafonov} gives a common minorant: a positive measure bounded above by both velocity marginals. The portions supplied by this lower bound cancel up to the phase error. This removes a fixed fraction of total variation in a fixed time at frequency one. At frequency $\xi$, kinetic scaling gives the same contraction over a time proportional to $\abs\xi^{-\frac23}$: a velocity change of size $r=\abs\xi^{-\frac13}$ produces a displacement of size $r^3$ over time $r^2$. Composition repeats the contraction and gives~\eqref{e:intro-decay}. The comparison of two measures with a common positive part is the feature shared with D\"oblin's argument~\cite{Doeblin1940}; here the loss comes from their opposite Fourier phases (Lemma~\ref{l:cancel}). The positive marginal in $(\sigma,v)$ is controlled by the parabolic Aleksandrov estimate~\cite[Theorem~2.4]{FerrettiSafonov}: its occupation density belongs to $L^\gamma$ for $1\leq\gamma\leq\frac{d+1}{d}$. Fourier inversion combines this bound with~\eqref{e:intro-decay} to give a full Green density $G$ for unit initial mass. With $\tau$ denoting elapsed time, it satisfies \[ \int_T^{2T}\int_{D\times\R^d}G(\tau,w,z)^q\dd w\dd z\dd\tau \leq C T^{1-2d(q-1)}. \] Integrating this inequality against $\frac{\dd T}{T}$ is finite near zero precisely when $q<1+\frac1{2d}$. This gives the finite-horizon Green estimate in the whole space. On a bounded velocity interval, exponential loss at the endpoints also makes the integral finite at infinity; for $d=1$ the same small-time condition is $q<\frac32$. Integrating the source against this Green density gives a non-negative potential. Subtracting it from the reflected sub-solution leaves a function bounded by its exit-boundary values. H\"older's inequality then gives the Aleksandrov estimate with $p=\frac q{q-1}$. For $A(t,v)$, the available time is at most $R^2$, and the density bound gives exactly the factor $R^{2-\frac{4d+2}p}$. We first prove this for smooth coefficients. For measurable coefficients, we keep the classical sub-solution fixed and approximate the coefficient on compact inner cylinders, where its velocity Hessian is bounded. Exhaustion then gives the estimate up to the kinetic boundary. For $a=a(x,v)$, the difficulty is to combine the estimates on velocity intervals as they approach zero. A first count of repeated visits already gives the Aleksandrov estimate at $p=6$. This implies H\"older continuity, from which propagation along kinetic paths gives a Harnack comparison for the full-space autonomous evolution. We then use this regularity to improve the estimate from which it came. The improvement concerns the timing of returns. A convex function of velocity bounds the time-integrated mass in a small velocity interval. Harnack prevents that mass from concentrating in a short time interval: a value at elapsed time $t$ is bounded by an average over $[2t,3t]$. Together these give the bound $C\frac r{\sqrt t}$ for mass in an interval of size $r$, when the starting velocity has size $r$ and $t\geq Cr^2$. Consequently, the mass beginning new visits during successive intervals of length $r^2$ is bounded by $C(1+j)^{-\frac12}$, with time measured from the first arrival in the inner interval. Summing its $q$-th powers gives the density bound $C R^{\frac2q-1}r^{\frac4q-3}$. Its integral against $\frac{\dd r}{r}$ is finite for $q<\frac43$, giving $p>4$. This recovers the range proved in the short paper. To reach below four, we also record the position at which each visit begins. A homogeneous barrier bounds occupation of boxes of width $r^3$ in position and $r$ in velocity. Harnack converts this integrated bound into faster decay for the mass in each box. In the density sum, this contributes an additional power $r^{\varepsilon(q-1)}$, making the sum converge even at $q=\frac43$ and slightly beyond. The preliminary estimate therefore improves through its own regularity consequences, together with the new barrier; the final theorem is used only at the end. For general $A(t,x,v)$, neither choice of time leaves a transported variable absent from the coefficient. The Fourier waves then interact, and the argument above does not apply. In the autonomous case, changing the clock removes this difficulty only away from zero velocity; the estimates on returns are what let us cross that remaining region. \subsection{Literature} For kinetic equations in divergence form, Golse, Imbert, Mouhot and Vasseur~\cite{GIMV} proved H\"older continuity and a Harnack inequality with measurable uniformly elliptic diffusion and bounded drift and source. Their cylinders are those of~\eqref{e:cylinder}. Guerand and Mouhot~\cite{GuerandMouhot} gave short quantitative De Giorgi proofs of the weak Harnack and Harnack inequalities, and Guerand and Imbert~\cite{GuerandImbert} gave a constructive proof using a logarithmic transform. In~\cite{GIMV} and~\cite{GuerandImbert}, the source is controlled by its supremum. In non-divergence form, Abedin and Tralli~\cite{AbedinTralli} prove invariant Harnack inequalities under a Cordes--Landis condition on the eigenvalues of the coefficient or under a uniform modulus of continuity. Dong and Yastrzhembskiy~\cite{DongYastrzhembskiy} prove global $L^p$ estimates and solvability when the coefficient is measurable in time but of vanishing mean oscillation in position and velocity. An Aleksandrov maximum principle for general measurable $A(t,x,v)$ is announced in the forthcoming work~\cite{NiebelForthcoming}. Hao and Zhang~\cite{HaoZhang2026} prove such an estimate through a probabilistic occupation bound and derive Harnack and H\"older estimates for strong solutions. The present argument instead uses Fourier cancellation in the two restricted coefficient classes above. Zhu~\cite{Zhu2026} proves interior H\"older estimates for classical stationary solutions in one space dimension with a merely measurable coefficient $a(x,v)$. The autonomous coefficient in Theorem~\ref{t:a} is the more general, since our solutions may depend on time, and we obtain a sharper threshold on the Lebesgue exponent. The short paper~\cite{DoeblinFourier} presents the cancellation argument, its Aleksandrov estimates, and a mixing theorem for diffusion on the sphere. This companion gives the detailed proofs of the Aleksandrov estimates and the passage to H\"older continuity. It also improves the one-dimensional autonomous range from $p>4$ to $p>p_*(\frac\Lambda\lambda)$ with $p_*<4$. The spherical mixing theorem is proved in the short paper. We use its kernel notation and denote the space dimension by $d$ here; the relation between the two conventions for solution operators is explained below. \subsection{Outline} Sections~\ref{s:additive}--\ref{s:reconstruction} concern the operator~\eqref{e:intro-Lop} alone: Section~\ref{s:additive} fixes the setting and the kernels, defines the Green measure and states the density theorem; Section~\ref{s:fourier} proves the decay~\eqref{e:intro-decay}; Section~\ref{s:zero} the occupation estimate at frequency zero; Section~\ref{s:reconstruction} combines them. Section~\ref{s:maximum} reflects the variables and proves what the two applications share; Sections~\ref{s:time-velocity} and~\ref{s:autonomous} prove Theorems~\ref{t:A} and~\ref{t:a}. Section~\ref{s:time-velocity} also obtains local spatial smoothing and a direct H\"older estimate for $A(t,v)$. Section~\ref{s:holder} shows that the Aleksandrov principle, in the localised form of Corollary~\ref{c:localised}, implies local H\"older continuity in the kinetic quasi-distance: Theorem~\ref{t:holder} states this for an arbitrary uniformly elliptic Borel coefficient $A(t,x,v)$, Corollary~\ref{c:holder-a} uses this implication for the autonomous class, while Corollary~\ref{c:holder-A} follows from the direct local argument. Appendix~\ref{s:construction} constructs the solution operators, the Duhamel potentials and the Green and exit measures of a velocity strip. Appendix~\ref{s:bellman} proves the homogeneous adjoint degree gap and constructs the Bellman barrier used for the improvement below four. Appendix~\ref{s:counterexample} constructs the autonomous counterexamples. \subsection{Notation} Smooth means $C^\infty$; $\Sym_d$ denotes the real symmetric $d\times d$ matrices, acting by $A:D_v^2u=\sum_{i,j}A_{ij}\partial_{v_iv_j}u$. For an open set $O\subseteq\R^{1+2d}$ we write $C^{1,1,2}(O)$ for the set of functions on $O$ whose derivatives $\partial_tu$, $\nabla_xu$, $\nabla_vu$ and $D_v^2u$ exist and are continuous in $O$; no condition is imposed at the boundary of $O$. We write $g_+=\max\{g,0\}$, and every $L^p$ norm is taken in all the displayed variables with respect to Lebesgue measure. Constants $C$ and $c$ depend only on the parameters listed as subscripts or named in the statement in which they occur. \bigskip \noindent \textbf{Acknowledgements.} The first author acknowledges support of the Fondation Sciences Mathématiques de Paris. \noindent \textbf{AI disclosure.} The mathematics is a product of more than a year of work, discussions, and eventually assistance of AI. This manuscript here is written by Codex. It accompanies the shorter human-written file \cite{DoeblinFourier}. \section{A parabolic equation with a transported variable} \label{s:additive} Both coefficient classes reduce to an equation with diffusion in $v$ and transport in $z$, whose coefficients do not depend on $z$. This independence lets each spatial frequency evolve separately. We introduce its positive evolution kernel first, then its Fourier transform and its time integral, the Green measure. The density estimate for the last of these will give the Aleksandrov principle. \subsection{Setting} We fix $d\geq1$. Throughout Sections~\ref{s:additive}--\ref{s:reconstruction} we consider a domain $D\subseteq\R^d$ and a smooth map $b:\R^d\to\R^d$, in exactly one of the following two situations. \[ D=\R^d,\qquad b(v)=v,\qquad m=L_b=1\,; \tag{W}\label{case:W} \] \[ d=1,\qquad D=J\ \text{a bounded open interval},\qquad m\leq b'\leq L_b\ \text{on }\R,\qquad 00$, and put \refstepcounter{equation} \[ \sigma=\sigma_*+r^2\tau,\qquad v=v_*+rY,\qquad z=z_*+(\sigma-\sigma_*)b(v_*)+r^3Z\,. \tag{\theequation}\label{e:scaling} \] Let \[ D_r=r^{-1}(D-v_*), \quad B_r(\tau,Y)=B(\sigma_*+r^2\tau,v_*+rY), \quad b_r(Y)=r^{-1}\bigl(b(v_*+rY)-b(v_*)\bigr). \] Then $(D_r,b_r)$ is again of type \eqref{case:W}, respectively \eqref{case:I}, with the same constants $m,L_b$ (and $\abs{J_r}=r^{-1}\abs J$), $B_r$ satisfies~\eqref{e:B-ellipticity}, and if $u(\sigma,v,z)=U(\tau,Y,Z)$ then \refstepcounter{equation} \[ \Lop u=r^{-2}\bigl(\partial_\tau+B_r:D_Y^2+b_r\cdot\nabla_Z\bigr)U\,. \tag{\theequation}\label{e:scaled-operator} \] Moreover, if $\Kf{\xi}$ and $\Kf{\xi}_{(r)}$ denote the Fourier kernels of the two operators, then for all $\sigma\leq\sigma'$, $v\in D$ and $\xi\in\R^d$, \refstepcounter{equation} \[ \norm{\Kf\xi_{\sigma,\sigma'}(v,\cdot)}_{\TV} = \norm{\Kf{r^3\xi}_{(r),\tau,\tau'}(Y,\cdot)}_{\TV}\,, \qquad \sigma=\sigma_*+r^2\tau,\ \sigma'=\sigma_*+r^2\tau',\ v=v_*+rY\,. \tag{\theequation}\label{e:scaled-tv} \] \end{lemma} \begin{proof} The chain rule gives~\eqref{e:scaled-operator}; the shift in $z$ subtracts the constant transport $b(v_*)$. Also $Db_r(Y)=Db(v_*+rY)$, so the structural bounds are unchanged. Uniqueness of the terminal problem identifies the original kernel with the affine image of the rescaled kernel. Since \[ z'-z=(\sigma'-\sigma)b(v_*)+r^3(Z'-Z), \] their Fourier projections satisfy \[ \Kf\xi_{\sigma,\sigma'}(v,E) =e^{-i(\sigma'-\sigma)\xi\cdot b(v_*)}\, \Kf{r^3\xi}_{(r),\tau,\tau'}\bigl(Y,r^{-1}(E-v_*)\bigr). \] The affine pushforward and the constant phase preserve total variation, proving~\eqref{e:scaled-tv}. \end{proof} \section{Decay of the non-zero Fourier modes} \label{s:fourier} We choose two curves in velocity space that carry a spatial wave into opposite phases. A fixed amount of the kernel stays near each curve. After a short further evolution, the two terminal velocity distributions have a common positive lower bound by parabolic Harnack. The corresponding portions of the kernel cancel up to a small phase error. The resulting contraction is uniform in the starting point, so it can be repeated. Higher frequencies require less time for the same phase separation, giving the following decay. \begin{proposition}[Decay of the Fourier kernels] \label{p:decay} There are constants $C,c>0$ depending only on the structural constants such that for all $\sigma<\tau$, $v\in D$ and $\xi\in\R^d$: \refstepcounter{equation} \[ \norm{\Kf\xi_{\sigma,\tau}(v,\cdot)}_{\TV} \leq C\exp\bigl(-c\,(\tau-\sigma)\abs\xi^{\frac23}\bigr) \text{in case \eqref{case:W}}, \tag{\theequation}\label{e:decay-W} \] \refstepcounter{equation} \[ \norm{\Kf\xi_{\sigma,\tau}(v,\cdot)}_{\TV} \leq C\exp\bigl(-c\,(\tau-\sigma)\bigl(1+\abs\xi^{\frac23}\bigr)\bigr) \text{in case \eqref{case:I}}. \tag{\theequation}\label{e:decay-I} \] \end{proposition} In case~\eqref{case:W}, the dependence on ellipticity can be written in the form used in the short paper: \[ \norm{\Kf\xi_{\sigma,\tau}(v,\cdot)}_{\TV} \leq C_{d,\kappa}\exp\left(-c_{d,\kappa}\lambda^{\frac13} (\tau-\sigma)|\xi|^{\frac23}\right), \qquad \kappa=\frac\Lambda\lambda. \] Indeed, the change $\widehat\sigma=\lambda^{\frac13}\sigma$, $\widehat v=\lambda^{-\frac13}v$, $\widehat z=z$ leaves transport in the form $\widehat v\cdot\nabla_{\widehat z}$ after division by $\lambda^{\frac13}$, and normalises the diffusion bounds to $[1,\kappa]$. Position and its frequency are unchanged, while pushing forward the velocity measure preserves total variation. A velocity change of size $r$ over time $r^2$ produces a relative displacement of size $r^3$. Taking $r=|\xi|^{-\frac13}$ therefore separates phases in time $|\xi|^{-\frac23}$. We prove a fixed contraction at unit frequency in Proposition~\ref{p:block}, using moving-ball barriers and the common minorant. Near the endpoints of a bounded interval, where the balls may not fit, boundary loss replaces cancellation. Section~\ref{ss:decay} then scales and iterates the contraction; boundary loss also supplies the decay at frequency zero. \subsection{Two barriers in a moving ball} \label{ss:barriers} The first barrier gives a lower bound on the mass remaining in $B_\rho(\gamma(r))$. The second bounds its displacement in $z$, and hence its Fourier phase. Both use the following comparison principle. The velocity domain may move along a piecewise $C^1$ curve. Since position ranges over all of $\R^d$, comparison also requires control at spatial infinity. \begin{lemma}[Maximum principle on a moving domain] \label{l:maxprinciple} Let $\Omega\subseteq\R^d$ be an open ball or, when $d=1$, a bounded open interval, let $\gamma:[\sigma,\tau]\to\R^d$ be continuous and $C^1$ on each interval of a subdivision $\sigma=r_00$ and the boundary values remain non-positive. We fix $\sigma'>\sigma$. If $\Psi_\varepsilon$ has a positive value on $\{r\geq\sigma'\}$, it attains a positive maximum there: in the first alternative omit $z$, and in the second use the uniform decay at infinity to restrict $z$ to a compact set. At a maximiser $(r_*,v_*,z_*)$, the terminal and velocity boundary conditions give $r_*<\tau$ and $v_*\in\Omega_{r_*}$. Hence \[ \nabla_v\Psi_\varepsilon=0,\qquad D_v^2\Psi_\varepsilon\leq0,\qquad \nabla_z\Psi_\varepsilon=0. \] The curve $r\mapsto(r,v_*+\gamma(r)-\gamma(r_*),z_*)$ stays in the moving domain for $r>r_*$ close to $r_*$. The right derivative at its maximum is non-positive, giving \[ \partial_r\Psi_\varepsilon+\gamma'(r_*)\cdot\nabla_v\Psi_\varepsilon =\partial_r\Psi_\varepsilon\leq0. \] This also covers a maximiser at $r_*=\sigma'$. Positive semidefiniteness of $B$ now gives $\Lop\Psi_\varepsilon\leq0$, a contradiction. We let $\sigma'\downarrow\sigma$ and then $\varepsilon\downarrow0$. For piecewise $C^1$ curves, we apply this argument successively on the subdivision intervals, starting with the last. Continuity supplies the terminal inequality for each preceding interval. \end{proof} \begin{lemma}[Retention in a moving ball] \label{l:retention} Let $\gamma:[\sigma,\tau]\to\R^d$ be continuous and piecewise $C^1$ with $\abs{\gamma'}\leq H$, and let $\rho>0$ be such that $B_\rho(\gamma(r))\subseteq D$ for $\sigma\leq r\leq\tau$. Let $K^{\gamma,\rho}$ be the kernel of $\Lop$ with zero boundary values on $\partial B_\rho(\gamma(r))\times\R^d$. Then for every $z\in\R^d$, \refstepcounter{equation} \[ K^{\gamma,\rho}_{\sigma,\tau}\bigl((\gamma(\sigma),z),B_\rho(\gamma(\tau))\times\R^d\bigr) \geq \exp\bigl(-C_{d,\lambda,\Lambda}\,(\rho^{-2}+H^2)(\tau-\sigma)\bigr)\,. \tag{\theequation}\label{e:retention} \] \end{lemma} \begin{proof} We extend $\gamma$ constantly outside $[\sigma,\tau]$. By Proposition~\ref{p:evolution}\ref{i:marginal}, the mass in~\eqref{e:retention} is the mass of the parabolic kernel $P^{\gamma,\rho}$ on the moving ball, independently of the starting position. We choose an integer $\mathsf N\geq3$ with $2(\mathsf N-1)\lambda>d\Lambda$, and set $c_{\mathsf N}=4\mathsf N(\mathsf N-1)\lambda-2\mathsf N d\Lambda>0$. For $00\}$. We choose a smooth increasing $\chi:\R\to[0,1]$, zero on $(-\infty,0]$ and equal to one above the positive minimum of $\varkappa_+$ on $\supp\phi$. Then $F=\chi(\varkappa_+)$ satisfies $\Lop F\leq0$ and dominates $\phi$ at time $\tau$. The terminal solution $u=\mathsf S^{\gamma,\rho}_{r,\tau}\phi$ is bounded between zero and one and vanishes on the velocity boundary. To compare it with $F$ on the unbounded position domain, put \[ h_\eta(r,z)=\eta(1+\abs{z-z_0}^2)e^{(M_b+1)(\tau-r)}, \qquad M_b=\sup_{\sigma\leq r\leq\tau,\;v'\in\overline{B_\rho(\gamma(r))}}\abs{b(v')}. \] The supremum is finite, and $2|z-z_0|\leq1+|z-z_0|^2$ gives $\Lop h_\eta\leq-h_\eta$. Thus $u-F-h_\eta$ has non-positive terminal and lateral values, tends uniformly to $-\infty$ as $|z|\to\infty$, and has non-negative $\Lop$. Lemma~\ref{l:maxprinciple} gives $u\leq F+h_\eta$. At $(\sigma,v,z_0)$, $F=0$, so \[ \int\phi\dd K^{\gamma,\rho}_{\sigma,\tau}((v,z_0),\cdot) \leq \eta e^{(M_b+1)(\tau-\sigma)}. \] We let $\eta\downarrow0$. Inner regularity shows that the kernel gives no mass to $\{\varkappa_+(\tau,z)>0\}$. Repeating the argument with $\varkappa_-$ excludes the other side and proves~\eqref{e:phase}. \end{proof} \subsection{A common minorant at the terminal time} \label{ss:minorant} The barriers retain a positive amount of mass in each tube. A lower bound on two total masses does not yet give a common positive measure below both distributions. We obtain this common minorant by letting both evolve in one fixed ball and applying parabolic Harnack. For each non-negative terminal function of velocity, its integral against the parabolic kernel is a non-negative solution of the parabolic equation in the starting time and velocity. We apply Harnack to this solution after reversing time. Because the comparison holds for every such terminal function, it is an inequality between measures. This is the lower bound used in \S~\ref{ss:cancel}; the Fourier kernels themselves are complex and play no part in this application of Harnack. \begin{theorem}[Parabolic Harnack inequality~\cite{KrylovSafonov}; see~{\cite[Theorem~1.5]{FerrettiSafonov}}] \label{t:harnack} Let $N\geq1$ and let $B:(0,1)\times B_1(0)\to\Sym_N$ be continuous with $\lambda I_N\leq B\leq\Lambda I_N$. There is $\mathfrak h=\mathfrak h(N,\lambda,\Lambda)>0$ such that every non-negative $\mathfrak u\in C^{1,2}((0,1)\times B_1(0))$ solving $\partial_\theta \mathfrak u=B(\theta,v):D_v^2\mathfrak u$ in $(0,1)\times B_1(0)$ satisfies \refstepcounter{equation} \[ \sup_{(\frac14,\frac12)\times B_{\frac12}(0)}\mathfrak u\ \leq\ \mathfrak h^{-1}\inf_{(\frac34,1)\times B_{\frac12}(0)}\mathfrak u\,. \tag{\theequation}\label{e:harnack} \] \end{theorem} The two cylinders in~\eqref{e:harnack} are concentric and separated in time, so one application compares two times only at points that are close together; a chain of applications along a segment removes that restriction and gives the two-point comparison used below. \begin{lemma}[A common terminal minorant] \label{l:minorant} Let $v_0\in D$ and $\rho>0$ with $D_0=B_{4\rho}(v_0)\subseteq D$, let $T\in\R$, and let $P^{D_0}_{\theta,T}(v,\cdot)$ be the kernel of $\partial_\sigma+B:D_v^2$ with zero boundary values on $\partial B_{4\rho}(v_0)$. There are $h,q_0>0$ depending only on $d,\lambda,\Lambda$ such that \refstepcounter{equation} \[ \begin{aligned} P^{D_0}_{T-4\rho^2,T}(v,\cdot)&\geq h\,P^{D_0}_{T-\rho^2,T}(v_0,\cdot) \quad\text{for every }v\in B_\rho(v_0),\\ P^{D_0}_{T-\rho^2,T}\bigl(v_0,B_{4\rho}(v_0)\bigr)&\geq q_0\,. \end{aligned} \tag{\theequation}\label{e:minorant} \] \end{lemma} \begin{proof} By parabolic scaling we may take $\rho=1$ and $v_0=0$, retaining $B$ for the rescaled coefficient with the same ellipticity bounds. For $0\leq\phi\in C_c^\infty(B_4)$, the function \[ \mathfrak u(\vartheta,v)=\int\phi\dd P^{B_4}_{T-\vartheta,T}(v,\cdot) \] is bounded, non-negative and smooth for $\vartheta>0$ by Proposition~\ref{p:evolution}, and solves $\partial_\vartheta\mathfrak u=B(T-\vartheta,v):D_v^2\mathfrak u$. We join $(1,0)$ to $(4,v)$, $|v|<1$, in six Harnack steps. We set $\vartheta_j=1+\frac j2$, $v_j=\frac j6v$ and $\bar v_j=\frac12(v_j+v_{j+1})$. The cylinders \[ (\vartheta_j-\frac38,\vartheta_j+\frac58)\times B_1(\bar v_j), \qquad j=0,\ldots,5, \] lie in $(0,\infty)\times B_4$. The successive points have time offsets $\frac38$ and $\frac78$ and lie within $\frac1{12}$ of $\bar v_j$. Theorem~\ref{t:harnack} therefore gives $\mathfrak u(1,0)\leq\mathfrak h^{-6}\mathfrak u(4,v)$. Testing with every such $\phi$ proves the first measure inequality in~\eqref{e:minorant}, with $h=\mathfrak h^6$. For the second inequality, Lemma~\ref{l:retention} on $B_3$, followed by domain monotonicity, bounds the mass from $(T-1,0)$ below by $q_0=e^{-\frac{C_{d,\lambda,\Lambda}}9}$. Scaling back proves both claims. \end{proof} \subsection{Cancellation of two phases} \label{ss:cancel} The lemma below applies to the portions selected using the common minorant. Their velocity marginal is $\Theta$, and their phases lie near two values differing by $\pi$. Replacing each phase by its central value gives opposite Fourier contributions. The error is small, so the decrease in total variation is proportional to $\Theta(D)$. \begin{lemma}[D\"oblin--Fourier cancellation] \label{l:cancel} Let $M$, $F_1$, $F_2$ be finite positive Borel measures on $D\times\R^d$ with $F_1+F_2\leq M$, and suppose that each $F_i$ has velocity marginal $\Theta$. Let $\xi\in\R^d$, $z_0\in\R^d$, $0\leq\varepsilon<1$ and $\Phi_1,\Phi_2\in\R$ with $e^{-i\Phi_1}+e^{-i\Phi_2}=0$, and suppose that $F_i$ is supported in $\{(w,z'):\abs{\xi\cdot(z'-z_0)-\Phi_i}\leq\varepsilon\}$ for $i=1,2$. Then the complex measure $E\mapsto\int_{E\times\R^d}e^{-i\xi\cdot(z'-z_0)}\dd M$ satisfies \refstepcounter{equation} \[ \Bigl\lVert\int_{\,\cdot\,\times\R^d}e^{-i\xi\cdot(z'-z_0)}\dd M\Bigr\rVert_{\TV} \leq M(D\times\R^d)-2(1-\varepsilon)\,\Theta(D)\,. \tag{\theequation}\label{e:cancel} \] \end{lemma} \begin{proof} For a finite positive measure $N$ on $D\times\R^d$, denote its Fourier projection on $D$ by $\widehat{N}$. The phase has modulus one, so \refstepcounter{equation} \[ \norm{\widehat{N}}_{\TV}\leq N(D\times\R^d). \tag{\theequation}\label{e:trivial-tv} \] If $h$ is a bounded complex Borel function on $D$ with $\abs h\leq1$, the common marginal and the phase support give \[ \left|\int_D h\,\dd\bigl(\widehat F_i-e^{-i\Phi_i}\Theta\bigr)\right| \leq\int_{D\times\R^d}\abs{e^{-i\xi\cdot(z'-z_0)}-e^{-i\Phi_i}}\dd F_i \leq\varepsilon\Theta(D). \] Taking the supremum over $h$ yields $\norm{\widehat F_i-e^{-i\Phi_i}\Theta}_{\TV}\leq\varepsilon\Theta(D)$. Since $e^{-i\Phi_1}+e^{-i\Phi_2}=0$, it follows that $\norm{\widehat F_1+\widehat F_2}_{\TV}\leq2\varepsilon\Theta(D)$. The remainder $M-F_1-F_2$ is positive and has mass $M(D\times\R^d)-2\Theta(D)$. Thus~\eqref{e:trivial-tv} and the triangle inequality give \[ \norm{\widehat{M}}_{\TV} \leq M(D\times\R^d)-2\Theta(D)+2\varepsilon\Theta(D), \] which is~\eqref{e:cancel}. \end{proof} \subsection{A contraction at unit frequency} \label{ss:block} We now choose the two paths and apply the preceding estimates. In the whole space there is room for the moving tubes at every starting velocity. In a bounded interval they may not fit near an endpoint, but there the zero boundary condition already removes a fixed amount of mass. The next lemma quantifies this loss by comparison with the one-dimensional heat equation on a half-line. \begin{lemma}[Killing near an endpoint] \label{l:killing} Let $d=1$, let $J'$ be a bounded open interval, and let $P$ be the kernel of $\partial_\sigma+B\partial_{vv}$ with zero boundary values at $\partial J'$. For $\ell>0$ and $v\in J'$, \refstepcounter{equation} \[ \begin{aligned} P_{\sigma,\sigma+\ell}(v,J') \leq \operatorname{erf}\Bigl(\frac{\dist(v,\partial J')}{2\sqrt{\lambda\ell}}\Bigr)\,, \qquad \operatorname{erf}(x)=\frac{2}{\sqrt\pi}\int_0^xe^{-t^2}\dd t\,. \end{aligned} \tag{\theequation}\label{e:killing} \] In particular $\sup_{v\in J'}P_{\sigma,\sigma+\abs{J'}^2}(v,J')\leq\operatorname{erf}\bigl(\tfrac{1}{4\sqrt\lambda}\bigr)<1$. \end{lemma} \begin{proof} We fix $\tau=\sigma+\ell$ and an endpoint $\mathfrak e$ nearest to $v$. The distance $\mathfrak d(v')=|v'-\mathfrak e|$ is affine on $J'$, with $|\mathfrak d'|=1$. The half-line heat solution \[ W(\vartheta,x)=\operatorname{erf}\left(\frac{x}{2\sqrt{\lambda\vartheta}}\right) \] satisfies $W_\vartheta=\lambda W_{xx}$ and $W_{xx}\leq0$ for $x\geq0$. Consequently \refstepcounter{equation} \[ W_\vartheta-BW_{xx}=(\lambda-B)W_{xx}\geq0. \tag{\theequation}\label{e:erf-supersolution} \] For $\varepsilon>0$, the function $\mathcal W(r,v')=W(\tau-r+\varepsilon,\mathfrak d(v'))$ is a supersolution for the terminal problem. It is non-negative at both endpoints and has terminal value $\varphi_\varepsilon(v')=W(\varepsilon,\mathfrak d(v'))$. Comparison with smooth terminal data $0\leq\phi\leq\varphi_\varepsilon$, compactly supported in $J'$, gives \[ \int\phi\dd P_{\sigma,\tau}(v,\cdot) \leq W(\ell+\varepsilon,\mathfrak d(v)). \] We first approximate $\varphi_\varepsilon$ from below by such data. We then let $\varepsilon\downarrow0$; since $\varphi_\varepsilon\uparrow1$ in $J'$, monotone convergence gives~\eqref{e:killing}. Finally, $\mathfrak d(v)\leq\frac12|J'|$ and $\ell=|J'|^2$ give the stated strict contraction. \end{proof} The four ingredients now combine, at the unit frequency and over a single elapsed time whose length is fixed by the structural constants alone. \begin{proposition}[A contraction at unit frequency] \label{p:block} We assume \eqref{case:W}, or \eqref{case:I} with $J$ replaced by an arbitrary bounded open interval $J'$ on which~\eqref{case:I} holds. There are $L>0$ and $\delta\in(0,1)$, depending only on $d,\lambda,\Lambda,m,L_b$, such that for every unit vector $\xi\in\R^d$, every $\sigma\in\R$ and every $v\in D$, \refstepcounter{equation} \[ \norm{\Kf\xi_{\sigma,\sigma+L}(v,\cdot)}_{\TV}\leq1-\delta\,. \tag{\theequation}\label{e:block} \] \end{proposition} The constants do not depend on the length or position of the interval. This is important because the change of scale in Lemma~\ref{l:scaling} changes that interval. \begin{proof} We denote by $C_*=C_{d,\lambda,\Lambda}$ the constant of Lemma~\ref{l:retention} and by $h,q_0$ the constants of Lemma~\ref{l:minorant}. We put \refstepcounter{equation} \[ L_0=2+\frac{2\pi}{m}, \qquad \rho=\min\Bigl\{1,\ (128L_b)^{-\frac13}\Bigr\}, \qquad \eta=\min\Bigl\{\frac \rho 4,\ \frac{1}{8L_bL_0}\Bigr\}\,, \tag{\theequation}\label{e:block-constants} \] and \refstepcounter{equation} \[ L=L_0+4\rho^2, \qquad R_*=1+4\rho, \qquad m_0=\exp\Bigl(-C_*\Bigl(\eta^{-2}+\tfrac14\Bigr)L_0\Bigr)\,. \tag{\theequation}\label{e:block-constants-two} \] All of these depend only on $d,\lambda,\Lambda,m,L_b$. We fix a unit vector $\xi\in\R^d$, a time $\sigma\in\R$, a point $z_0\in\R^d$ and a point $v_0\in D$, and assume until Step~5 that $B_{R_*}(v_0)\subseteq D$; in case~\eqref{case:W} this holds for every $v_0$. We must bound the total variation of the complex measure $\Kf\xi_{\sigma,\sigma+L}(v_0,\cdot)$, which by~\eqref{e:fourier-kernel} is the Fourier transform, at the frequency $\xi$ and with base point $z_0$, of the positive measure \[ M=K_{\sigma,\sigma+L}\bigl((v_0,z_0),\cdot\bigr), \qquad M\bigl(D\times\R^d\bigr)\leq1\,, \] the bound on the mass being part of Proposition~\ref{p:evolution}~\ref{i:kernel}. \emph{Step 1: two paths with opposite phases.} For $t\in[0,L_0]$ we set $\mathsf f(t)=\min\{t,1,L_0-t\}$. Since $L_0\geq2$, the function $\mathsf f$ is continuous and piecewise affine with $\mathsf f(0)=\mathsf f(L_0)=0$, $0\leq\mathsf f\leq1$, $\abs{\mathsf f'}\leq1$ and \[ \int_0^{L_0}\mathsf f(t)\dd t=L_0-1=1+\frac{2\pi}{m}\,. \] For $\alpha\in[0,\frac12]$ define curves on $[\sigma,\sigma+L_0]$ by \[ \gamma_1(r)=v_0, \qquad \gamma_2^\alpha(r)=v_0+\alpha\,\mathsf f(r-\sigma)\,\xi\,. \] Both curves begin and end at $v_0$, have speed at most $\frac12$, and stay within distance $\frac12$ of $v_0$. We extend them constantly by $v_0$ outside the interval. Since $\eta\leq\frac14$, their tubes fit inside $D$: \refstepcounter{equation} \[ B_\eta\bigl(\gamma_i(r)\bigr)\subseteq B_{\frac12+\eta}(v_0)\subseteq B_1(v_0)\subseteq B_{R_*}(v_0)\subseteq D, \qquad i=1,2. \tag{\theequation}\label{e:corridor-inclusion} \] We define \[ \Delta(\alpha)=\int_0^{L_0}\xi\cdot\bigl(b(v_0+\alpha\mathsf f(t)\xi)-b(v_0)\bigr)\dd t. \] The bound~\eqref{e:b-bounds} gives \[ \Delta'(\alpha)\geq m\int_0^{L_0}\mathsf f(t)\dd t=m+2\pi. \] Since $\Delta(0)=0$ and $\Delta(\frac12)>\pi$, we can choose $\alpha\in(0,\frac12)$ with $\Delta(\alpha)=\pi$. We fix this value, write $\gamma_2=\gamma_2^\alpha$, and set \[ a_i=\int_\sigma^{\sigma+L_0}\xi\cdot b(\gamma_i(r))\dd r, \qquad a_2-a_1=\pi. \] \emph{Step 2: restricting the kernel to the two moving tubes.} For $i=1,2$, we denote by $K^{\gamma_i,\eta}$ the kernel of $\Lop$ on $[\sigma,\sigma+L_0]$ with zero boundary values on $\partial B_\eta(\gamma_i(r))\times\R^d$. It is available by~\eqref{e:corridor-inclusion}. From $(v_0,z_0)$, its terminal measure is carried by $B_\eta(v_0)\times\R^d$, since $\gamma_i(\sigma+L_0)=v_0$, and has mass at most one by Proposition~\ref{p:evolution}~\ref{i:kernel}. Lemma~\ref{l:retention}, applied with the curve $\gamma_i$, $H=\frac12$ and radius $\eta$, gives the lower bound \refstepcounter{equation} \[ m_0\leq m_i:=K^{\gamma_i,\eta}_{\sigma,\sigma+L_0} \bigl((v_0,z_0),D\times\R^d\bigr)\leq1\,. \tag{\theequation}\label{e:corridor-mass} \] We retain mass $\frac{m_0}{2}$ from each tube by setting \[ R_i:=\frac{m_0}{2m_i}K^{\gamma_i,\eta}_{\sigma,\sigma+L_0} \bigl((v_0,z_0),\cdot\bigr), \qquad R_i(D\times\R^d)=\frac{m_0}{2}. \] Domain monotonicity bounds each tube kernel by the full kernel. Since $\frac{m_0}{2m_i}\leq\frac12$, each $R_i$ is at most half of that kernel. Hence \refstepcounter{equation} \[ R_1+R_2\leq K_{\sigma,\sigma+L_0}((v_0,z_0),\cdot), \tag{\theequation}\label{e:corridors-dominated} \] even if the two tubes overlap. \emph{Step 3: selecting portions using a common minorant.} We put $T=\sigma+L$ and $D_0=B_{4\rho}(v_0)\subset D$. We continue $R_i$ with the kernel killed at $\partial D_0$: \[ E_i(E)=\int K^{D_0}_{\sigma+L_0,T}((v',z'),E)R_i(\dd v'\dd z'). \] The starting velocities lie in $B_\eta(v_0)\subset D_0$. Domain monotonicity,~\eqref{e:corridors-dominated} and composition give \refstepcounter{equation} \[ E_1+E_2\leq K_{\sigma,T}((v_0,z_0),\cdot)=M. \tag{\theequation}\label{e:continued-dominated} \] The velocity marginals are \[ \nu_i(E)=\int P^{D_0}_{\sigma+L_0,T}(v',E)R_i(\dd v'\dd z'). \] Since $\eta<\rho$, Lemma~\ref{l:minorant} applies throughout this integral. It gives \[ \nu_i\geq\Theta:=\frac{m_0h}{2}P^{D_0}_{T-\rho^2,T}(v_0,\cdot), \qquad \Theta(D)\geq\frac{m_0hq_0}{2}>0. \] We set $g_i=\frac{\dd\Theta}{\dd\nu_i}$, so $0\leq g_i\leq1$, and define \[ F_i(\dd w\dd z')=g_i(w)E_i(\dd w\dd z'). \] At each terminal velocity $w$, this keeps the fraction $g_i(w)$ of $E_i$. This operation is sometimes called trimming; it is a weighted restriction, which need not be a restriction to a set. It preserves every retained phase, gives $F_1+F_2\leq M$, and gives each $F_i$ the velocity marginal $\Theta$. The original marginals $\nu_1$ and $\nu_2$ need not agree. \emph{Step 4: the phase error and the cancellation.} We set $\Phi_i=a_i+4\rho^2\xi\cdot b(v_0)$, so $\Phi_2-\Phi_1=\pi$. Lemma~\ref{l:phase} bounds the phase error during the tube by $L_b\eta L_0$. During the common continuation, the velocity stays within $4\rho$ of $v_0$ for time $4\rho^2$, adding at most $16L_b\rho^3$. Thus $E_i$, and hence $F_i$, is supported where \[ |\xi\cdot(z'-z_0)-\Phi_i| \leq L_b\eta L_0+16L_b\rho^3\leq\frac14. \] Lemma~\ref{l:cancel} now gives \refstepcounter{equation} \[ \norm{\Kf\xi_{\sigma,\sigma+L}(v_0,\cdot)}_{\TV} \leq1-\frac32\Theta(D)\leq1-\delta_1, \qquad \delta_1=\frac34m_0hq_0. \tag{\theequation}\label{e:block-interior} \] This proves the contraction whenever the two tubes fit, in particular at every starting velocity in the whole space. \emph{Step 5: points near an endpoint.} In case~\eqref{case:I}, a point not covered above satisfies $\dist(v,\partial J')\leq R_*$. Domination by the velocity marginal and Lemma~\ref{l:killing} give \refstepcounter{equation} \[ \norm{\Kf\xi_{\sigma,\sigma+L}(v,\cdot)}_{\TV} \leq P_{\sigma,\sigma+L}(v,J') \leq\operatorname{erf}\left(\frac{R_*}{2\sqrt{\lambda L}}\right) =:1-\delta_2<1. \tag{\theequation}\label{e:block-boundary} \] This bound is independent of the interval's length and location. We take $\delta=\min\{\delta_1,\frac12\}$ in case~\eqref{case:W} and $\delta=\min\{\delta_1,\delta_2,\frac12\}$ in case~\eqref{case:I}. The interior and boundary estimates cover every starting velocity. All constants depend only on $d,\lambda,\Lambda,m,L_b$, as required. \end{proof} \subsection{Proof of the decay} \label{ss:decay} At frequency $\xi\ne0$, the choice $r=|\xi|^{-\frac13}$ turns Proposition~\ref{p:block} into a fixed contraction over time $Lr^2$. Composition repeats this contraction. On a bounded interval, Lemma~\ref{l:killing} also gives decay at frequency zero; we combine the two bounds by splitting the time interval in half. \begin{proof}[Proof of Proposition~\ref{p:decay}] For $\xi\ne0$, Lemma~\ref{l:scaling} with $r=|\xi|^{-\frac13}$ reduces the frequency to $r^3\xi=\frac\xi{|\xi|}$. The ellipticity and transport bounds are unchanged, so Proposition~\ref{p:block} and~\eqref{e:scaled-tv} give \refstepcounter{equation} \[ \norm{\Kf\xi_{s,s+Lr^2}(v,\cdot)}_{\TV}\leq1-\delta \qquad\text{for all }s\in\R,\ v\in D\,. \tag{\theequation}\label{e:block-any-frequency} \] We compose $k=\lfloor\frac{\tau-\sigma}{Lr^2}\rfloor$ such blocks. The remaining factor has norm at most one by Lemma~\ref{l:fourier-kernels}~\ref{i:kato}. Thus~\eqref{e:tv-composition} gives a bound $(1-\delta)^k$, and, with $C_1=(1-\delta)^{-1}$ and $c_1=\frac{-\log(1-\delta)}L$, \refstepcounter{equation} \[ \norm{\Kf\xi_{\sigma,\tau}(v,\cdot)}_{\TV} \leq (1-\delta)^{\frac{\tau-\sigma}{Lr^2}-1} = C_1\exp\bigl(-c_1(\tau-\sigma)\abs\xi^{\frac23}\bigr)\,. \tag{\theequation}\label{e:decay-frequency} \] The mass bound covers $\xi=0$, proving~\eqref{e:decay-W}. In case~\eqref{case:I}, Lemma~\ref{l:killing} gives \[ \sup_{v\in J}P_{s,s+|J|^2}(v,J)\leq1-\delta_0, \qquad \delta_0=1-\operatorname{erf}\Bigl(\frac1{4\sqrt\lambda}\Bigr)>0. \] The same iteration, now in blocks of length $|J|^2$, and $|\Kf\xi|\leq P$ imply \refstepcounter{equation} \[ \begin{aligned} &\norm{\Kf\xi_{\sigma,\tau}(v,\cdot)}_{\TV} \leq P_{\sigma,\tau}(v,J)\leq C_2\exp\bigl(-c_2(\tau-\sigma)\bigr) \end{aligned} \tag{\theequation}\label{e:decay-killing} \] where $C_2=(1-\delta_0)^{-1}$ and $c_2=\frac{-\log(1-\delta_0)}{|J|^2}$. Finally, we split $[\sigma,\tau]$ at its midpoint and apply~\eqref{e:decay-frequency} to the first factor and~\eqref{e:decay-killing} to the second. This proves~\eqref{e:decay-I} with $C=C_1C_2$ and $c=\frac12\min\{c_1,c_2\}$. \end{proof} \section{The velocity marginal} \label{s:zero} Integrating the positive kernel in $z$ removes transport and leaves the parabolic kernel $P_{\sigma,\tau}(v,\cdot)$. This is also its Fourier transform at frequency zero. The parabolic Aleksandrov estimate, by duality, bounds the density of the time integral of $P$ in $L^{\frac{d+1}{d}}$. Unlike the mass bound alone, this controls concentration in time and velocity. We also check that the estimate is uniform over the whole velocity space, as needed in Section~\ref{s:reconstruction}. \begin{theorem}[Parabolic Aleksandrov estimate~\cite{Krylov1976}; see~{\cite[Theorem~2.4]{FerrettiSafonov}}] \label{t:krylov} Let $N\geq1$, $00$. There is a non-negative density $g_\rho^S$ on $(0,S)\times D$ characterised, for bounded Borel $\phi$, by \refstepcounter{equation} \[ \int_0^S\!\int_D\phi(\tau,w)g_\rho^S(\tau,w)\dd w\dd\tau=\int_D\rho(\dd v)\int_0^S\!\!\int_D\phi(\tau,w)\,P_{\sigma_0,\sigma_0+\tau}(v,\dd w)\dd\tau \tag{\theequation}\label{e:occupation} \] and, for $1\leq\gamma\leq\gamma_0\coloneqq \frac{d+1}{d}$, \refstepcounter{equation} \[ \norm{g_\rho^S}_{L^\gamma((0,S)\times D)}\leq C_{d,\lambda,\Lambda,\gamma}\,M\,S^{\beta_\gamma}, \qquad \beta_\gamma=\frac{d+2}{2\gamma}-\frac d2\,. \tag{\theequation}\label{e:occupation-bound} \] \end{proposition} Proposition~\ref{p:occupation} is the dual form of the parabolic Aleksandrov estimate; the short localisation below makes its constant independent of the spatial extent of the source. \begin{proof} We first take $\sigma_0=0$, $S=1$ and $D=\R^d$. For $0\leq F\in C_c^\infty((0,1)\times\R^d)$, let \[ W(r,v)=\int_r^1\!\int_{\R^d}F(t,w)\,P_{r,t}(v,\dd w)\dd t, \qquad H=\sup_{[0,1]\times\R^d}W. \] Then $W\geq0$, $H\leq\norm F_\infty$, $W(1,\cdot)=0$ and $\partial_rW+B:D_v^2W=-F$. We choose $R^2=\max\{1,4d\Lambda\}$. For any $0\leq\alpha<1$ and $v_*$, let $V$ solve the same source equation on $(\alpha,1)\times B_R(v_*)$, with zero terminal and lateral values. Parabolic solvability~\cite[Theorem~5.14]{Lieberman} gives a classical solution $V\in C^{1,2}$ in the open cylinder, continuous on its closure; no regularity up to the boundary is used. After reflection $s=1-r$ and scaling, put $S=1-\alpha$, $\rho=(1-2\delta/S)^{\frac12}$ and $T_\delta(s,v)=(\delta+\rho^2s,\ v_*+\rho(v-v_*))$. Then $T_\delta$ maps the closed cylinder into the open one, so $V_\delta=V\circ T_\delta$ is $C^{1,2}$ on a neighbourhood of the closed cylinder and solves the same type of equation, with ellipticity bounds unchanged and source $\rho^2F\circ T_\delta$, whose $L^{d+1}$-norm is at most $\norm F_{L^{d+1}}$. Let $\eta_\delta=\sup|V_\delta|$ over the parabolic boundary; $\eta_\delta\to0$ by the uniform continuity of $V$ and $V=0$ there. Theorem~\ref{t:krylov} applied to $V_\delta-\eta_\delta$ gives $V_\delta\leq\eta_\delta+CR^{\frac d{d+1}}\norm F_{L^{d+1}}$, and $\delta\downarrow0$ yields \[ 0\leq V\leq CR^{\frac d{d+1}}\norm F_{L^{d+1}}. \] The quadratic function \[ \psi(r,v)=\frac{\abs{v-v_*}^2+2d\Lambda(1-r)}{R^2} \] is non-negative, at least one on the lateral boundary, and satisfies $(\partial_r+B:D_v^2)\psi\leq0$. Thus comparison gives $W\leq V+H\psi$ on this cylinder. At $(\alpha,v_*)$, $\psi\leq\frac12$. Taking the supremum over $\alpha,v_*$ and absorbing yields \[ H\leq CR^{\frac d{d+1}}\norm F_{L^{d+1}}+\frac12H, \qquad H\leq C_{d,\lambda,\Lambda}\norm F_{L^{d+1}}. \] This proves the whole-space bound without a time iteration. The right side of~\eqref{e:occupation}, with $S=1$ and $\phi=F$, equals $\int W(0,v)\rho(\dd v)$. Duality therefore gives the density $g_\rho^1$ with $\norm{g_\rho^1}_{L^{\frac{d+1}{d}}}\leq CM$. Its $L^1$ norm is at most $M$ by contraction. Interpolation gives every $1\leq\gamma\leq\gamma_0$. For general $S$, we use $(\tau,v)=(S\hat\tau,\sqrt S\,\hat v)$ and push forward the initial measure, retaining its mass $M$. Uniqueness identifies the rescaled parabolic kernels, and the time integral contributes a factor $S$. Consequently \[ g_\rho^S(\tau,v)=S^{-\frac{d}{2}}\widetilde g\left(\frac\tau S,\frac v{\sqrt S}\right), \qquad \norm{g_\rho^S}_{L^\gamma} =S^{\frac{d+2}{2\gamma}-\frac d2}\norm{\widetilde g}_{L^\gamma}, \] which proves~\eqref{e:occupation-bound}. Finally, the killed interval kernel is dominated by the whole-line kernel with the same coefficient, so the same bound holds there. \end{proof} \section{Reconstruction of the Green density} \label{s:reconstruction} The two preceding sections control different parts of the same positive measure. The parabolic estimate bounds its integral in position; Fourier decay bounds its height in position. Their combination gives an $L^q$ bound in all variables. We first work away from the starting time, on $T<\tau<2T$, and then integrate the resulting power of $T$ down to zero. Throughout this section, $\mu$ has mass $M$ and is prescribed at time $\sigma_0$. On each interval $T<\tau<2T$, the Green measure $\Gamma_\mu^\infty$ has mass at most $MT$. We retain $\gamma_0=\frac{d+1}{d}$ and $\beta_\gamma=\frac{d+2}{2\gamma}-\frac d2$ from Proposition~\ref{p:occupation}. \begin{lemma}[Fourier marginals on one elapsed-time interval] \label{l:slab-fourier} Let $\mu$ be as in Theorem~\ref{t:green} and $T>0$. On $T<\tau<2T$, the $(\tau,w)$-marginal of $\Gamma_\mu^\infty$ has a density $g_{\mu,T}$ with \refstepcounter{equation} \[ \norm{g_{\mu,T}}_{L^\gamma((T,2T)\times D)}\leq C_\gamma\,M\,T^{\beta_\gamma}, \qquad 1\leq\gamma\leq\gamma_0\,. \tag{\theequation}\label{e:slab-zero} \] For $\xi\in\R^d$, the complex measure $\Phi\mapsto\int_{\{T<\tau<2T\}}e^{-i\xi\cdot z}\Phi(\tau,w)\dd\Gamma_\mu^\infty$ on $(T,2T)\times D$ has a density $k^\xi_{\mu,T}$, and \refstepcounter{equation} \[ \begin{aligned} \norm{k^\xi_{\mu,T}}_{L^\gamma((T,2T)\times D)} &\leq C_\gamma\,M\,T^{\beta_\gamma}\,\epsilon_T(\xi),\\ \epsilon_T(\xi)&= \begin{cases} \exp\bigl(-cT\abs\xi^{\frac23}\bigr)&\text{in case \eqref{case:W}},\\ \exp\bigl(-cT(1+\abs\xi^{\frac23})\bigr)&\text{in case \eqref{case:I}}, \end{cases} \end{aligned} \tag{\theequation}\label{e:slab-fourier} \] for $1\leq\gamma\leq\gamma_0$. Consequently \refstepcounter{equation} \[ \frac{1}{(2\pi)^d}\int_{\R^d}\norm{k^\xi_{\mu,T}}_{L^\gamma((T,2T)\times D)}\dd\xi \leq C_\gamma\,M\,T^{\beta_\gamma-\frac{3d}2}\times \begin{cases}1&\text{in case \eqref{case:W}},\\ e^{-cT}&\text{in case \eqref{case:I}}.\end{cases} \tag{\theequation}\label{e:slab-height} \] \end{lemma} \begin{proof} The $(\tau,w)$-marginal is the restriction to $(T,2T)$ of the parabolic occupation measure issued from $\mu_v(E)=\mu(E\times\R^d)$. Proposition~\ref{p:occupation}, with $S=2T$, gives~\eqref{e:slab-zero}. For the Fourier marginal, we evolve first to elapsed time $\frac T2$ and define the complex velocity measure \[ \eta_T^\xi(E)=\int_{D\times\R^d}e^{-i\xi\cdot z}\, \Kf\xi_{\sigma_0,\sigma_0+\frac T2}(v,E)\,\mu(\dd v\dd z). \] Proposition~\ref{p:decay} gives $\norm{\eta_T^\xi}_{\TV}\leq CM\epsilon_T(\xi)$, after decreasing the rate $c$. We continue its positive total variation $|\eta_T^\xi|$ by the parabolic kernel $P$, and let $g$ be the resulting occupation density over elapsed times $(0,\frac{3T}{2})$. Composition of the Fourier kernels and their domination by $P$ give, for bounded Borel $\Phi$, \[ \left|\int_{\{T<\tau<2T\}}e^{-i\xi\cdot z}\Phi(\tau,w)\dd\Gamma_\mu^\infty\right| \leq\int_T^{2T}\!\int_D|\Phi(\tau,w)|\, g\left(\tau-\frac T2,w\right)\dd w\dd\tau. \] Thus the Fourier marginal has a density bounded in modulus by this translate of $g$. Proposition~\ref{p:occupation} gives~\eqref{e:slab-fourier}. This uses the parabolic estimate only on a positive measure, not on a complex solution. Finally, $\zeta=T^{\frac{3}{2}}\xi$ gives \[ \int_{\R^d}e^{-cT|\xi|^{\frac{2}{3}}}\dd\xi=C_{d,c}T^{-\frac{3d}{2}}, \] which proves~\eqref{e:slab-height}. Measurability of the integrand follows by expressing the $L^\gamma$ norm as a supremum over a countable norming family of bounded test functions; each associated Fourier pairing is continuous in $\xi$. \end{proof} The next lemma combines the marginal and Fourier bounds to construct the full density $G$. The input densities $g$ and $k^\xi$ are then its position integral and Fourier transform: \[ g=\int_{\R^d}G(\cdot,z)\dd z, \qquad k^\xi=\int_{\R^d}e^{-i\xi\cdot z}G(\cdot,z)\dd z. \] \begin{lemma}[Reconstruction from the marginal and the Fourier bounds] \label{l:reconstruction} Let $Y$ be a Borel subset of a Euclidean space with a $\sigma$-finite Borel measure $m$, and let $\Gamma$ be a finite positive Borel measure on $Y\times\R^d$ whose $Y$-marginal is $g\,\dd m$. Let $1q>(d+1)(q-1)$, both $\gamma$ and $q_1$ lie in $[1,\gamma_0]$. We apply Lemma~\ref{l:reconstruction} to the marginal bounds~\eqref{e:slab-zero} and~\eqref{e:slab-height}. The resulting power of $T$ is \[ (q-1)\left(\beta_\gamma-\frac{3d}{2}\right)+\beta_{q_1} =\frac{d+2}{2}\left(\frac{q-1}{\gamma}+\frac1{q_1}\right)-2d(q-1)-\frac d2 =1-2d(q-1). \] In case~\eqref{case:I}, the exponential factor is raised to $q-1$; decreasing $c$ gives~\eqref{e:slab}. \end{proof} \begin{proof}[Proof of Theorem~\ref{t:green}] The local densities from Proposition~\ref{p:slab} agree on overlaps, giving a density $G$ for $\Gamma_\mu^\infty$ on positive elapsed times. We put $\delta_q=1-2d(q-1)>0$. We integrate the estimates over all interval lengths: for every non-negative function $F$, \[ \int_0^S F(\tau)\dd\tau \leq\frac1{\log2}\int_0^S\frac{\dd T}{T}\int_T^{2T}F(\tau)\dd\tau. \] Indeed, each $\tau\in(0,S)$ is counted for $\frac\tau20$, Chebyshev's inequality on the two time ranges therefore gives \[ |\{G>a\}|\leq\frac ha+C a^{-q}h^{1-2d(q-1)}. \] Choosing $h=a^{-\frac1{2d}}$ bounds this by $Ca^{-q_*}$. Linearity gives $\|G_\mu^\infty\|_{L^{q_*,\infty}}\leq CM$ for initial mass $M$. Lorentz-space H\"older duality then bounds the source potential by $C\|f\|_{L^{2d+1,1}}$. The comparison argument of the next section gives the corresponding endpoint maximum principle for smooth $A(t,v)$. \end{remark} \section{From Green measures to maximum principles} \label{s:maximum} A positive source produces a non-negative potential. Subtracting that potential from a sub-solution leaves a function controlled by its boundary values; Hölder's inequality then turns a Green-density bound into the Aleksandrov estimate. This section makes that comparison on the kinetic boundary, where data are prescribed only on the initial, velocity and incoming position faces. We first reflect time and position to use the terminal kernels constructed above. We set $s=-t$, $X=-x$ and $U(s,X,v)=u(-s,-X,v)$, leaving the velocity unreflected. Given $Z_0=(s_0,X_0,v_0)\in\R\times\R^d\times\R^d$ and $R>0$, the forward cylinder is \refstepcounter{equation} \[ \mathcal Q_R(Z_0)=\bigl\{(s,X,v):\ s_00$, then \refstepcounter{equation} \[ \sup_{\overline{\mathcal Q}}Z\leq\sup_{\partial^+\mathcal Q}Z_+\,. \tag{\theequation}\label{e:comparison} \] \end{lemma} \begin{proof} We set $Z_\varepsilon=Z+\varepsilon(s-s_0-R^2)$. Then $Z_\varepsilon\leq Z$ and $\KA Z_\varepsilon\geq\varepsilon$ wherever $Z_\varepsilon>0$. In relative coordinates the operator is $\partial_s+V\cdot\nabla_Y+A:D_V^2$. We suppose that a positive maximum of $Z_\varepsilon$ on $\overline{\mathcal Q_\delta}$ occurs outside its exit boundary. The velocity is interior, so $D_V^2Z_\varepsilon\leq0$ and $A:D_V^2Z_\varepsilon\leq0$. The time is either interior or at the initial face; in the latter case its right derivative is non-positive. Hence $\partial_sZ_\varepsilon\leq0$. Finally, either the position is interior, or it lies on $|Y|=\rho_\delta^3R^3$ with $V\cdot Y<0$. In the second case $Y+hV$ lies inside the ball for small $h>0$, so the one-sided derivative gives $V\cdot\nabla_YZ_\varepsilon\leq0$; in the first this term is zero. Thus $\KA Z_\varepsilon\leq0$, a contradiction. It follows that $\sup_{\overline{\mathcal Q_\delta}}Z_\varepsilon\leq\sup_{\partial^+\mathcal Q_\delta}(Z_\varepsilon)_+$. We use~\eqref{e:retraction} and continuity to send $\delta\downarrow0$, then use $Z\leq Z_\varepsilon+\varepsilon R^2$ and send $\varepsilon\downarrow0$. \end{proof} The comparison uses only pointwise positive semidefiniteness of $A$. The following proposition adds a source potential with a uniform density bound, as constructed in Sections~\ref{s:time-velocity} and~\ref{s:autonomous}. \begin{proposition}[From a Green density to the Aleksandrov estimate] \label{p:abp-abstract} Let $\mathcal Q=\mathcal Q_R(Z_0)$, let $A$ be smooth on a neighbourhood of $\overline{\mathcal Q}$ with $\lambda I_d\leq A\leq\Lambda I_d$, let $10$. Suppose that to every $P\in\mathcal Q$ is associated a positive Borel measure $\Gamma_P$ on $\mathcal Q$ such that: \begin{enumerate} \item\label{i:potential} for every $0\leq g\in C_c^\infty(\mathcal Q)$, the function $W_g(P)=\int_{\mathcal Q}g\dd\Gamma_P$ belongs to $C^{1,1,2}(\mathcal Q)$, is non-negative, and satisfies $\KA W_g=-g$ in $\mathcal Q$; \item\label{i:density} $\Gamma_P$ has a density $G_P$ with $\norm{G_P}_{L^q(\mathcal Q)}\leq\mathsf K$ for every $P\in\mathcal Q$. \end{enumerate} Let $U\in C(\overline{\mathcal Q})\cap C^{1,1,2}(\mathcal Q)$ and let $g_0\geq0$ be a Borel function with $g_0\in L^p(\mathcal Q)$ and $\KA U\geq-g_0$ almost everywhere in $\mathcal Q$. Then \refstepcounter{equation} \[ \sup_{\overline{\mathcal Q}}U\leq\sup_{\partial^+\mathcal Q}U_+ +\mathsf K\norm{g_0\1_{\{U>0\}}}_{L^p(\mathcal Q)}\,. \tag{\theequation}\label{e:abp-abstract} \] \end{proposition} \begin{proof} The defect $g=(-\KA U)_+$ is continuous, satisfies $g\leq g_0$ almost everywhere, and makes $\KA U+g\geq0$ pointwise. We fix $\varepsilon>0$ and choose $0\leq\theta\leq1$ smooth, with $\theta=0$ below $\frac\varepsilon2$ and $\theta=1$ above $\varepsilon$. Then \[ 0\leq\theta(U)g\leq g_0\1_{\{U>0\}}\quad\text{almost everywhere}. \] We choose $\delta_j\downarrow0$. On each inner cylinder $\mathcal Q_j=\mathcal Q_{\delta_j}$ choose $0\leq g_j\in C_c^\infty(\mathcal Q)$ majorising $\theta(U)g$ on $\overline{\mathcal Q_j}$, with \[ \norm{g_j}_{L^p}\leq\norm{g_0\1_{\{U>0\}}}_{L^p}+\frac1j. \] Such majorants are obtained using a cutoff between zero and one on a neighbourhood of $\overline{\mathcal Q_j}$, mollifying, and adding its vanishing uniform approximation error times a fixed cut-off equal to one there. We set $W_j(P)=\int g_j\dd\Gamma_P$. Wherever $U-W_j-\varepsilon>0$, non-negativity of $W_j$ gives $U>\varepsilon$, hence $g_j\geq g$ and $\KA(U-W_j-\varepsilon)=\KA U+g_j\geq0$. Lemma~\ref{l:comparison} and H\"older's inequality yield, for $P\in\mathcal Q_j$, \[ U(P)\leq\sup_{\partial^+\mathcal Q_j}U_+ +\mathsf K\left(\norm{g_0\1_{\{U>0\}}}_{L^p}+\frac1j\right)+\varepsilon. \] The retractions~\eqref{e:retraction} and continuity pass the boundary term to $\partial^+\mathcal Q$. We let $j\to\infty$ and then $\varepsilon\downarrow0$. \end{proof} \begin{remark}[Viscosity subsolutions for smooth coefficients] \label{r:viscosity} For smooth $A$ and continuous $g_0\geq0$, Proposition~\ref{p:abp-abstract} also holds for continuous viscosity subsolutions of $-\KA U\leq g_0$, without classical derivatives. Consequently Theorems~\ref{t:A} and~\ref{t:a}, including Corollary~\ref{c:localised}, hold for smooth coefficients, after reflection, for continuous viscosity subsolutions with continuous source and with the same constants. Only comparison needs checking. In the relative coordinates of Lemma~\ref{l:comparison}, put $s_1=s_0+R^2$ and $M=\sup_{\partial^+\mathcal Q}Z_+$. If $\sup Z>M$, choose $\varepsilon>0$ so that $Z_\varepsilon=Z+\varepsilon(s-s_1)$ retains a positive excess at some point, and maximise \[ Z_\varepsilon-M-\frac{\delta}{s-s_0} -\frac{\delta}{R^6-|Y|^2}. \] For small $\delta$ the maximum retains a fixed positive excess and lies in the interior. Since $Z\geq Z_\varepsilon>0$ there, the smooth upper test gives \[ \varepsilon\leq-\frac{\delta}{(s-s_0)^2} +\frac{2\delta\,V\cdot Y}{(R^6-|Y|^2)^2}. \] As $\delta\downarrow0$, a limit of these maximisers cannot lie on the exit boundary, by continuity and the fixed positive excess. If its position is interior, the positive part of the spatial term tends to zero; if it lies on a position face, that face is incoming, so $V\cdot Y<0$. Both cases contradict the display. We let $\varepsilon\downarrow0$ to obtain comparison. The preceding proof now applies with smooth majorants of $\theta(U)g_0$ in place of $\theta(U)(-\KA U)_+$. This argument does not establish a viscosity extension for Borel coefficients. The approximation below uses the locally bounded classical velocity Hessian and cannot simply be transferred to that solution class. \end{remark} Our Green measures are constructed for smooth coefficients. The following lemma passes their uniform Aleksandrov estimates to Borel coefficients by approximation on compact inner cylinders, keeping the solution fixed. \begin{lemma}[From smooth to Borel coefficients] \label{l:borel} Let $1\leq p<\infty$ and $\alpha\in\R$, and let $\mathcal C$ be one of the two coefficient classes \begin{enumerate} \item[(a)] Borel maps $A:\R\times\R^d\to\Sym_d$ with $\lambda I_d\leq A(t,v)\leq\Lambda I_d$ almost everywhere, regarded as coefficients of $\PA$; \item[(b)] $d=1$ and Borel maps $a:\R\times\R\to[\lambda,\Lambda]$, regarded as coefficients of $\Pa$. \end{enumerate} We write $\mathcal L$ for the operator $\PA$, respectively $\Pa$, with the coefficient under consideration. Suppose that for every smooth member of $\mathcal C$ satisfying its bounds everywhere, every $P_0\in\R^{1+2d}$ and $R>0$, and every $u\in C(\overline{Q_R^-(P_0)})\cap C^{1,1,2}(Q_R^-(P_0))$ and Borel $f$ with $\mathcal Lu\leq f$ almost everywhere and $f_+\in L^p$, one has \refstepcounter{equation} \[ \sup_{\overline{Q_R^-(P_0)}}u\leq\sup_{\partial_{\kin}Q_R^-(P_0)}u_+ +C_0R^{\alpha}\norm{f_+}_{L^p(Q_R^-(P_0))} \tag{\theequation}\label{e:smooth-estimate} \] with $C_0$ independent of the coefficient. Then~\eqref{e:smooth-estimate} holds, with the same $C_0$, for every member of $\mathcal C$ and every such $u$. The same implication holds when $\norm{f_+}_{L^p(Q_R^-(P_0))}$ is replaced by $\norm{f_+\1_{\{u>0\}}}_{L^p(Q_R^-(P_0))}$ in hypothesis and conclusion. \end{lemma} \begin{proof} In case (a), we redefine $A=\lambda I_d$ on the null set where ellipticity fails. We mollify the coefficient in $(t,v)$ in case (a), and in $(x,v)$ in case (b), obtaining smooth coefficients $C_j$ in the same class, with the same ellipticity bounds and $C_j\to C$ almost everywhere. Tonelli's theorem lifts this convergence to phase space. We keep $u$ fixed and put \[ e_j=|(C_j-C):D_v^2u|. \] Then $\mathcal L_{C_j}u\leq f+e_j$ almost everywhere and $e_j\leq2\sqrt d\,\Lambda|D_v^2u|$, with $d=1$ in case (b). We take $Q_\delta$ to be the reflection of~\eqref{e:inner-explicit}; it is a kinetic cylinder of radius $\rho_\delta R$ with closure inside $Q_R^-(P_0)$. For fixed $\delta$, the classical Hessian is bounded there, so dominated convergence gives $\norm{e_j}_{L^p(Q_\delta)}\to0$. We apply~\eqref{e:smooth-estimate} on $Q_\delta$ and first let $j\to\infty$, using $(f+e_j)_+\leq f_++e_j$. This gives \[ \sup_{\overline{Q_\delta}}u\leq\sup_{\partial_{\kin}Q_\delta}u_++C_0(\rho_\delta R)^\alpha\norm{f_+}_{L^p(Q_\delta)}. \] We now let $\delta\downarrow0$. The reflected retractions~\eqref{e:retraction} carry the original kinetic boundary onto $\partial_{\kin}Q_\delta$ and converge uniformly to the identity. Continuity of $u$ therefore passes the boundary supremum to the limit; exhaustion passes the left-hand side and source norm, and $\rho_\delta R\to R$. No uniform Hessian bound as $\delta\downarrow0$ is needed. For the localised form, use $(f+e_j)_+\1_{\{u>0\}}\leq f_+\1_{\{u>0\}}+e_j$ on each inner cylinder and take the same two limits. \end{proof} \section{Coefficients depending on time and velocity} \label{s:time-velocity} Here the coefficient is independent of position on the whole space. Reflection gives the operator~\eqref{case:W}, with $z=X$, $B(s,v)=A(-s,v)$ and $b(v)=v$. The Green-density estimate therefore applies directly: H\"older's inequality bounds the source potential, and comparison gives the maximum principle. The second part of this section uses the same Fourier decay to obtain local smoothness in position and then H\"older continuity. \begin{proof}[Proof of Theorem~\ref{t:A}] For smooth $A$, reflection~\eqref{e:reflection} gives $U(s,X,v)=u(-s,-X,v)$ with $\mathcal K_{\widetilde A}U\geq-g_0$, where $g_0(s,X,v)=f_+(-s,-X,v)$. It preserves the source norm and maps the kinetic boundary onto the exit boundary of $\mathcal Q=\mathcal Q_R(-t_0,-x_0,v_0)$. For $P\in\mathcal Q$, we denote by $\Gamma_P$ the whole-space Green measure from $P$, restricted to $\mathcal Q$ and written in absolute time. Lemma~\ref{l:duhamel} supplies the potential identity in Proposition~\ref{p:abp-abstract}. With $q=\frac p{p-1}<1+\frac1{2d}$ and remaining time $S_P0$ and $\alpha\in(0,1)$, depending only on $d,\lambda,\Lambda$, such that, for $P=(t,x,v)$ and $P'=(s,y,w)$ in $Q_{\frac{R}{2}}^-(P_0)$, \[ |u(P)-u(P')| \leq C\operatorname{osc}_{Q_R^-(P_0)}u \left(\frac{|t-s|^{\frac{1}{2}}+|v-w|+|x-y-(t-s)v_0|^{\frac{1}{3}}}R\right)^\alpha. \] Here $C_m$ depends only on $d,\lambda,\Lambda,m$. \end{proposition} \begin{proof} We first suppose that $A$ is smooth and work with the reflected operator $\KA=\partial_s+v\cdot\nabla_X+A(s,v):D_v^2$. For a starting point $(s,X,v)$ in the image of $Q_{\frac{3R}{4}}^-(P_0)$ under reflection, put $D=B_R(v_0)$ and $T=s+\frac18R^2$. This terminal time is still inside the reflected outer cylinder; the spatial containment is checked below. \emph{The local representation.} Proposition~\ref{p:evolution} gives the kernel $K^D$ killed at $\partial D$. For a source $F$, we write $W_F^D$ for the Duhamel potential obtained from this kernel. It also defines the solution with prescribed terminal and lateral data: for a smooth compactly supported extension $\phi$ of those data, set \[ \mathsf E\phi=\phi+W_{\KA\phi}^D. \] Here signed potentials are defined by subtracting the potentials of the positive and negative parts. We use the bounded-source distributional identity proved in Appendix~\ref{s:construction}; it gives $\KA \mathsf E\phi=0$, and Theorem~\ref{t:hormander} gives interior smoothness. The quadratic barrier gives \[ |W_F^D(r,X,v)|\leq \norm F_\infty \min\left\{T-r,\frac{R^2-|v-v_0|^2}{2d\lambda}\right\}; \] indeed $\KA(R^2-|v-v_0|^2)=-2\operatorname{tr}A\leq-2d\lambda$. Approximation extends this bound from compactly supported smooth sources to bounded sources. Thus $\mathsf E\phi$ has the prescribed traces. Bounded comparison, using the growth barrier of Appendix~\ref{s:construction}, shows that it depends only on those traces and is a positive contraction. Smooth ambient traces are uniformly dense in continuous compactly supported boundary data: extend such data continuously with compact support and mollify. Uniform limits of their solutions remain distributionally homogeneous, hence smooth inside. Riesz representation therefore gives a positive boundary measure $\omega^D_{s,v}$ on \[ \bigl(\{T\}\times\R^d\times\overline D\bigr) \cup\bigl([s,T]\times\R^d\times\partial D\bigr). \] We write its position variable as $y=X'-X-(r-s)v_0$; translation invariance removes the starting position $X$. For a unit vector $e$, we put $\varkappa_\pm=\pm e\cdot y-R(r-s)$. Then $\KA\varkappa_\pm\leq0$. If non-negative boundary data $g\leq1$ have compact support where $\varkappa_+>0$, choose a smooth increasing $\chi$ with $\chi(0)=0$ and $\chi(\varkappa_+)\geq g$ on the boundary. This is a supersolution, so bounded comparison gives $\mathsf Eg\leq\chi(\varkappa_+)$ and $\mathsf Eg(s,X,v)=0$. The growth barrier justifies comparison on the unbounded position domain exactly as in Lemma~\ref{l:phase}. Inner regularity, the same argument for $\varkappa_-$, and a countable dense set of directions give \[ |y|\leq R(r-s)\leq\frac18R^3 \quad\text{on the support of }\omega^D_{s,v}. \] It gives the same support for the killed kernels. We choose $\eta\in C_c^\infty(\R^d)$ equal to one on a neighbourhood of $\{|y|\leq R(T-s)\}$ and use $\phi(r,X',w)=\eta(X'-X-(r-s)v_0)$, with time and velocity cutoffs equal to one near the closed strip. Then $\KA\phi=0$ on the killed Green support, so \[ 1=\mathsf E\phi(s,X,v)=\int\phi\dd\omega^D_{s,v}=\int1\dd\omega^D_{s,v}. \] Thus the boundary measure has mass one. Bounded classical homogeneous solutions are represented by this measure: cut off their boundary values outside the cone and use comparison. There is no initial-velocity-face mass for an interior starting velocity, as also follows from the short-time estimate below. \emph{Decay before the boundary is reached.} For $v\in B_{\frac{3R}{4}}(v_0)$ and $04$. We then retain the position of each visit. A homogeneous Bellman barrier supplies the additional concentration bound needed for every $p>p_*(\frac\Lambda\lambda)$. Classical parabolic Harnack is used for the Fourier minorant after changing the clock; the kinetic comparison controls returns in the original variables. Throughout this section $d=1$ and $a\in C^\infty(\R^2)$ satisfies $\lambda\leq a\leq\Lambda$ everywhere; the passage to a Borel coefficient is made at the very end, by Lemma~\ref{l:borel}. By the reflection convention of Section~\ref{s:maximum}, the forward coefficient is $(X,v)\mapsto a(-X,v)$, which is again smooth, independent of $s$, and $[\lambda,\Lambda]$-valued; we write $a$ for it as well, so that the operator~\eqref{e:KA} reads \[ \Ka U=\partial_sU+v\partial_XU+a(X,v)\partial_{vv}U\,. \] Its coefficient depends on the transported variable $X$, so $\Ka$ is not an operator of the form~\eqref{e:Lop} and Theorem~\ref{t:green} does not apply to it directly. We work on velocity strips unbounded in position, so that no position boundary interrupts the change of time variable. The Green measure records the evolution inside a strip, and the exit measure records the data needed when it reaches the strip boundary. Composing these measures will account for repeated returns to smaller velocity intervals. For a bounded open interval $H$ and $-\infty0\}$ and \refstepcounter{equation} \[ (\chi_e)_\#\bigl(\abs v\,\Gamma_e\bigr)=\abs{\bar v}\,r^2\,\widetilde\Gamma_e\,. \tag{\theequation}\label{e:pushforward} \] Consequently $\Gamma_e$ has a density $G_e$ with respect to $\dd s\dd X\dd v$, and for $10\}$. We take $0\leq F\in C_c^\infty((0,\infty)\times\R\times J^*)$, vanishing for $\sigma\geq\sigma_1$. We write $h$ for its Duhamel potential for $\Lop_e$. By Lemma~\ref{l:duhamel}, $h$ is smooth, non-negative, solves $\Lop_eh=-F$, vanishes for $\sigma\geq\sigma_1$, and has zero velocity-boundary values. It is bounded uniformly over all starting points by Theorem~\ref{t:green}\ref{i:green-I} and H\"older's inequality, using any fixed exponent in $(1,\frac32)$. We put $\phi=h\circ\chi_e$. It belongs to the test class of Proposition~\ref{p:strip}, and \[ \Ka\phi=-\frac{|v|}{|\bar v|r^2}\,F\circ\chi_e. \] We apply the Green identity on the strip with terminal time $s_+\geq s_e+2r^2\sigma_1$, and denote its Green measure by $\Gamma_{e,s_+}$. Its velocity interval is $I^*$. The exit term is zero: $\phi$ vanishes on the velocity faces, and on the support of the terminal exit measure the cone gives $\sigma\geq\sigma_1$. Hence, with $y_e=\frac{v_e-\bar v}{r}$, \[ |\bar v|r^2\int F\dd\widetilde\Gamma_e =|\bar v|r^2h(0,0,y_e) =\int |v|\,(F\circ\chi_e)\dd\Gamma_{e,s_+}. \] We let $s_+\uparrow\infty$. Both limiting measures are finite by~\eqref{e:strip-mass} and Theorem~\ref{t:green}\ref{i:green-I}; equality against all such $F$ therefore proves~\eqref{e:pushforward}. If $\widetilde G_e$ is the density of $\widetilde\Gamma_e$, the pushforward identity and the Jacobian~\eqref{e:clock-jacobian} give \[ G_e=\frac{|\bar v|}{|v|}r^{-4}\widetilde G_e\circ\chi_e,\qquad \int |G_e|^q \leq 2^q r^{6-4q}\int |\widetilde G_e|^q. \] Theorem~\ref{t:green}\ref{i:green-I} now gives~\eqref{e:pole-density}. \end{proof} The estimate for one starting point does not distinguish visits beginning at different times. The next corollary keeps track of these times. The quadratic velocity barrier makes the mass remaining in $I^*$ decay exponentially on the time scale $r^2$. Thus visits beginning far apart in time contribute little to the same part of the Green density. \begin{corollary}[Starting at different times] \label{c:entrance} In the setting of Lemma~\ref{l:clock}, let $\nu$ be a finite positive Borel measure supported in $\R\times\R\times I^*$, and let $\Gamma_\nu=\int\Gamma_e\dd\nu(e)$, where $\Gamma_e$ is as in Proposition~\ref{p:pushforward}. We put $m_j=\nu((jr^2,(j+1)r^2]\times\R\times I^*)$. Then $\Gamma_\nu$ has a density $G_\nu$, and for $10$, and write $\mathcal Q=\mathcal Q_R(Z_0)$ for the forward cylinder~\eqref{e:forward-cylinder}, $J_Q=(v_0-R,v_0+R)$ for its velocity interval, and $\mathscr S=\mathscr S_{J_Q}$ for the strip with $s_-=s_0$ and $s_+=s_0+R^2$, so that $\mathcal Q\subset\mathscr S$. We retain the domain superscript in $\Gamma_P^{J_Q}$ and $\omega_P^{J_Q}$. We use $P\in\mathcal Q$ for a starting point of these two measures, and $e$ for a starting point of the all-time measures of Proposition~\ref{p:pushforward}. At every time the killed kernel has mass at most one. Its Green measure therefore satisfies $\Gamma_P^{J_Q}(J\times\R\times J_Q)\leq |J|$ for every time interval $J$. We will use this bound to count visits beginning during a short interval, as well as over the whole time interval. We also need to control how much time is spent in a narrow velocity interval. A convex function whose second derivative is concentrated on that interval gives the following bound. In particular, the Green measure cannot put mass on the slice $v=0$, where position cannot serve as time. \begin{lemma}[Velocity occupation of the strip measure] \label{l:capacity} For every $P\in\mathcal Q$ and every bounded interval $E\subset\R$, \refstepcounter{equation} \[ \Gamma_P^{J_Q}\bigl(\mathscr S\cap\{v\in E\}\bigr)\leq\frac{2R}{\lambda}\abs E\,. \tag{\theequation}\label{e:capacity} \] In particular $\Gamma_P^{J_Q}(\{v=v_1\})=0$ for every $v_1\in\R$. \end{lemma} \begin{proof} For $|E|>0$, we choose $0\leq\chi\in C_c^\infty(\R)$ with $\chi\geq\1_E$ and $\int\chi\leq2|E|$. The smooth function $\Phi(v)=\int|v-\zeta|\chi(\zeta)\dd\zeta$ satisfies $\Phi''=2\chi$ and $\|\Phi'\|_\infty\leq2|E|$. Both $\Phi$ and $\Ka\Phi=2a\chi$ are bounded on the strip, so Proposition~\ref{p:strip}\ref{i:green-identity} gives \[ 2\lambda\Gamma_P^{J_Q}(\mathscr S\cap\{v\in E\}) \leq\int 2a\chi\dd\Gamma_P^{J_Q} =\int\bigl(\Phi(v)-\Phi(v_P)\bigr)\dd\omega_P^{J_Q} \leq4R|E|, \] since the exit measure has mass one and its velocities lie in the interval $\overline{J_Q}$ of length $2R$. This proves~\eqref{e:capacity}; shrinking $E$ to a point gives the final assertion. \end{proof} \needspace{8\baselineskip} The starting velocity may be outside $I$, and mass may leave and return repeatedly. We alternate evolution inside $I^*$ with evolution outside $\overline I$. The measure below records where and when each visit to $\overline I$ begins, including the initial point if it already lies in $\overline I$. A new visit is counted only after the mass has left $I^*$ and returned to $\overline I$. The gap between the two boundaries prevents an exit from being counted immediately as a return. Figure~\ref{fig:ka-visits} illustrates the decomposition. \begin{figure}[t] \centering \includegraphics{figures/ka2_visits.pdf} \caption{\small{Counting visits (Lemmas~\ref{l:entrance} and~\ref{l:timed-entrance}), with $\bar v=2r$. \textbf{(a)} One velocity path, drawn with a slower diffusion for legibility. A visit begins when the velocity reaches $\overline I$ and ends when it leaves $I^*$; re-crossing $\partial I$ during a visit does not begin another. The time bins $J_j$, of length $r^2$, are counted from the first arrival in $\overline I$. \textbf{(b)} The mass $m_j$ of visits beginning in $J_j$, estimated from $40\,000$ paths started in $\overline I$, for one rough autonomous coefficient. The orange reference curve has $C$ chosen to bound the plotted Monte Carlo estimates. The masses count repeated visits, so their sum exceeds one. The bound $m_j\leq C(1+j)^{-\frac{1}{2}}$ uses the kinetic Harnack comparison of Lemma~\ref{l:return-time}, which rests on the preliminary $p=6$ estimate; this time-only count gives Proposition~\ref{p:cylinder-density}, hence $p>4$.}} \label{fig:ka-visits} \end{figure} \begin{lemma}[Counting visits to a velocity interval] \label{l:entrance} Let $P\in\mathcal Q$ and let $I_{\mathrm c}\subset I\subset I^*$ be as in Lemma~\ref{l:clock} with $I_{\mathrm c}\cap J_Q\neq\emptyset$. We put $H_1=I^*\cap J_Q$ and $H_2=J_Q\setminus\overline I$, and let \[ \Sigma_{\mathrm{in}}=(s_0,s_0+R^2)\times\R\times(\partial I\cap J_Q), \qquad \Sigma_{\mathrm{out}}=(s_0,s_0+R^2)\times\R\times(\partial I^*\cap J_Q)\,. \] For an empty strip, the integrated Green and exit measures below are understood to be zero. We define finite positive measures $\gamma_k,\beta_k$ ($k\geq1$) by \refstepcounter{equation} \[ \gamma_1= \begin{cases} \delta_P&\text{if }v_P\in\overline I,\\ \omega^{H_2}_{\delta_P}|_{\Sigma_{\mathrm{in}}}&\text{otherwise}, \end{cases} \qquad \beta_k=\omega^{H_1}_{\gamma_k}|_{\Sigma_{\mathrm{out}}}, \qquad \gamma_{k+1}=\omega^{H_2}_{\beta_k}|_{\Sigma_{\mathrm{in}}}\,. \tag{\theequation} \] Then $\gamma_I=\sum_{k\geq1}\gamma_k$, which counts all visits, is a finite measure on $\{v\in\overline I\}$ with \refstepcounter{equation} \[ \gamma_I(\mathscr S)\leq C_{\lambda,\Lambda}\Bigl(1+\frac Rr\Bigr)\,. \tag{\theequation}\label{e:entrance-mass} \] For every time interval $J$ of length $r^2$ it also satisfies \refstepcounter{equation} \[ \gamma_I\bigl(J\times\R\times J_Q\bigr)\leq C_{\lambda,\Lambda}\,. \tag{\theequation}\label{e:entrance-time} \] Moreover, \refstepcounter{equation} \[ \Gamma_P^{J_Q}|_{\{v\in I_{\mathrm c}\}}\leq\Gamma_{\gamma_I}|_{\{v\in I_{\mathrm c}\}}\,, \tag{\theequation}\label{e:entrance-domination} \] where $\Gamma_{\gamma_I}$ is the measure of Corollary~\ref{c:entrance}. Consequently $\Gamma_P^{J_Q}|_{\{v\in I_{\mathrm c}\}}$ has a density $G_{P,I}$ with \refstepcounter{equation} \[ \norm{G_{P,I}}_{L^q(\mathscr S)}\leq C_{\lambda,\Lambda,q}\,r^{\frac6q-4}\Bigl(1+\frac Rr\Bigr)^{\frac1q}, \qquad 14. \end{aligned} \] The preceding estimates already give the starting point, a preliminary Aleksandrov principle at $p=6$. Indeed, at $q=\frac65$, Lemma~\ref{l:entrance} bounds the density on a core of scale $r\leq CR$ by $CR^{\frac56}r^{\frac16}$. The continuous covering~\eqref{e:log-cover} below therefore gives \[ CR^{\frac56}\int_0^{CR}r^{\frac16}\frac{\dd r}{r}\leq CR. \] There is no mass at $v=0$ by Lemma~\ref{l:capacity}; cylinders with $|v_0|>8R$ are covered by a single interval in Proposition~\ref{p:pushforward}. Proposition~\ref{p:abp-abstract} thus gives the localised maximum principle at $p=6$. Applied to $\psi-u$ and $\psi+u$ for a homogeneous solution, it verifies~\eqref{e:amp}. We may consequently use the conditional H\"older theorem and propagation lemma of Section~\ref{s:holder}, with constants depending only on $\lambda,\Lambda$. This uses only the preliminary estimate, before either improvement. We write $\mathsf S_t$ for the positive full-space evolution with generator $v\partial_X+a(X,v)\partial_{vv}$. The terminal kernels on $(-M,M)$ from Proposition~\ref{p:strip} increase with $M$. For a starting point with velocity $v_e$, testing with $(v-v_e)^2$ bounds the mass lost before time $T$ by $\frac{2\Lambda T}{(M-|v_e|)^2}$. Their limit therefore has mass one. The bounded-limit and comparison argument of Appendix~\ref{s:construction} gives its classical solutions and composition. Autonomy then gives $\partial_t(\mathsf S_tF)=v\partial_X(\mathsf S_tF)+a\partial_{vv}(\mathsf S_tF)$ for bounded smooth $F$. Position reflection, $u(t,x,v)=\mathsf S_tF(-x,v)$, gives $\partial_tu+v\partial_xu-a(-x,v)\partial_{vv}u=0$, the convention of Section~\ref{s:holder}. Here we write functions in the physical coordinate order $(X,v)$. In terms of the full-space kernel constructed above, \[ \mathsf S_tF(X,v)=\int F(X',w)\, K_{0,t}\bigl((v,X),\dd w\dd X'\bigr). \] For a non-negative function $F$, this is its integral against the kernel at elapsed time $t$. If $F$ is the indicator of a velocity interval, it is the mass in that interval. The next lemma bounds it by a fixed multiple of its value at the same point at any later time in $[2t,3t]$. Averaging will then turn a time-integrated occupation estimate into a bound at time $t$. \begin{lemma}[A return-time comparison] \label{l:return-time} For every bounded non-negative Borel function $F$ on $\R^2$, \refstepcounter{equation} \[ \mathsf S_tF(X,v)\leq C_{\lambda,\Lambda}\mathsf S_sF(X,v), \qquad 2t\leq s\leq3t,\quad |v|\leq\sqrt t. \tag{\theequation}\label{e:return-time} \] \end{lemma} \begin{proof} The argument has three parts. Propagation carries a positive lower bound from a small neighbourhood to a later point. H\"older continuity supplies that neighbourhood from a positive value. An iteration removes the larger supremum that occurs in the H\"older estimate. \emph{A small positive neighbourhood reaches a later point.} We work in the reflected position coordinate, with a non-negative solution $u$ of $\Pa u=0$. We suppose $u\geq k$ on $Q_\ell^-(z)$, where $z=(\tau,x_c,v_c)$ lies in a fixed bounded region. After scaling by $\ell$ and following the transport of its centre, a cubic position path joins this patch to any point of \[ Q_{2\ell}^-(\tau+8\ell^2,x_c+8\ell^2v_c,v_c). \] The travel times lie between $4\ell^2$ and $8\ell^2$, and the normalised endpoints and accelerations are bounded. On the preceding interval of length $\frac14\ell^2$, we take the path to follow the centre of $Q_\ell^-(z)$, with tube widths $\frac12\ell^3$ and $\frac12\ell$. The initial tube is then contained in that cylinder. Lemma~\ref{l:propagation}, with $T_0=\frac14\ell^2$ and $T_1=8\ell^2$, gives a fixed factor $c_*>0$ throughout the later cylinder. We repeat with radii $\ell,2\ell,4\ell,\ldots$ until reaching a fixed small radius. The elapsed time is bounded by a constant times the square of that final radius, while the number of steps is at most $C+\log_2\frac1\ell$. One last path of fixed width reaches any prescribed point $z_*$ in a bounded future region separated by a fixed positive time. Choosing the final radius small enough leaves time for this last path. Thus, for some $\beta>0$ depending only on ellipticity and the fixed regions, \refstepcounter{equation} \[ u(z_*)\geq c\ell^\beta k. \tag{\theequation}\label{e:patch-power} \] All tubes lie in full space. Each doubling loses only a fixed factor in height, and the number of doublings grows like $\log\frac1\ell$, so the total loss is a power of $\ell$. This is why we enlarge the neighbourhood as we propagate it; keeping its width equal to $\ell$ for the whole time would give a much weaker bound. \emph{H\"older continuity supplies the neighbourhood.} We first take $F$ smooth, scale $t$ to one and translate $X$ to zero, without shifting velocity. For the reflected solution $u$, fix the later point $z_*=(s,0,v)$ with $2\leq s\leq3$ and $|v|\leq1$. We choose nested open boxes $U_R$, $\frac12\leq R\leq1$, around $(1,0,v)$, whose times lie in $(\frac{15}{16},\frac{17}{16})$ and whose position and velocity radii are small fixed multiples of $R$. They can be chosen so that $Q_{2h}^-(z)\subset U_b$ for $z\in U_a$, where $h=c(b-a)$ and $a0$. Taking $c$ small makes the required surrounding cylinders lie in $U_b$. Applying~\eqref{e:patch-power} to this neighbourhood and taking points with $u(z)$ approaching $M_a$ gives \[ u(z_*)\geq c h^\beta M_a^{1+\frac\beta\alpha}M_b^{-\frac\beta\alpha}. \] \emph{Removing the larger supremum.} The last inequality still contains $M_b$, because H\"older continuity was measured relative to that supremum. We rearrange it and use Young's inequality, with $\theta=\frac{\alpha}{\alpha+\beta}$, to obtain \[ M_a\leq Ch^{-\beta\theta}u(z_*)^\theta M_b^{1-\theta} \leq\varepsilon M_b+C_\varepsilon h^{-\beta}u(z_*). \] We set $a_j=1-2^{-j-1}$, so $a_0=\frac12$ and $h_j=c(a_{j+1}-a_j)=c2^{-j-2}$. Repeated substitution gives \[ M_{\frac12}\leq\varepsilon^N M_{a_N} +C_\varepsilon u(z_*)\sum_{j=0}^{N-1}(\varepsilon2^\beta)^j. \] We choose $\varepsilon2^\beta<1$ and let $N\to\infty$. Boundedness of $u$ removes the first term, and the sum stays bounded. This proves $u(1,0,v)\leq Cu(z_*)$. Scaling back proves~\eqref{e:return-time}. Approximation of the finite terminal measures extends the inequality from smooth non-negative $F$ to bounded non-negative Borel $F$. \end{proof} We use this comparison only for the full-space evolution, not for a kernel killed at the boundary of a one-sign velocity interval. The latter cannot in general return to the same position at a later time. \begin{lemma}[The improved density bound] \label{l:timed-entrance} In Lemma~\ref{l:entrance}, suppose in addition that $|\bar v|=2r$ and $r\leq6R$. Then \refstepcounter{equation} \[ \norm{G_{P,I}}_{L^q(\mathscr S)} \leq C_{\lambda,\Lambda,q}R^{\frac2q-1}r^{\frac4q-3}, \qquad 10$ and $\varsigma=\pm1$, we use~\eqref{e:intervals} with $\bar v=\varsigma\rho$ and $r=\frac\rho2$. Its core is $I^\varsigma_{\mathrm c}(\rho)=\{v:\frac{7\rho}8\leq\varsigma v\leq\frac{9\rho}8\}$. For $\varsigma v>0$, \refstepcounter{equation} \[ \int_0^\infty\1_{I^\varsigma_{\mathrm c}(\rho)}(v)\frac{\dd\rho}\rho =\int_{\frac{8\varsigma v}9}^{\frac{8\varsigma v}7}\frac{\dd\rho}\rho =\log\frac97\,, \tag{\theequation}\label{e:log-cover} \] When $|v|\leq9R$, the integral is unchanged if truncated at $L_R=\frac{72R}7$. For $\rho\leq L_R$, Lemma~\ref{l:timed-entrance} applies since $r\leq\frac{36}7R<6R$. For a non-negative Borel function $F$ on $\mathcal Q$, extended by zero to $\mathscr S$, Tonelli's theorem and H\"older's inequality give, with $q'=\frac q{q-1}$, \[ \int_{\mathcal Q\cap\{\varsigma v>0\}}F\dd\Gamma_P^{J_Q} \leq C_{\lambda,\Lambda,q}R^{\frac2q-1}\|F\|_{L^{q'}(\mathcal Q)} \int_0^{L_R}\rho^{\frac4q-3}\frac{\dd\rho}\rho \leq C_{\lambda,\Lambda,q}R^{\frac6q-4}\|F\|_{L^{q'}(\mathcal Q)}. \] The integral converges precisely for $q<\frac43$. Lemma~\ref{l:capacity} gives zero mass to $\{v=0\}$; summing the two signs and using duality proves~\eqref{e:cylinder-density}. If $|v_0|>8R$, we instead use one triple with $\bar v=v_0$ and $r=4R\leq\frac23|v_0|$. Its core contains $J_Q$ and $I$ contains $\overline{J_Q}$. Domain and time monotonicity bound $\Gamma_P^{J_Q}$ by the all-time measure on $I^*$, so Proposition~\ref{p:pushforward} gives the result for every $11$ and $\varepsilon=\gamma-1$; these are uniform upper bounds, not the rates of this example. \textbf{(c)} The overlapping cores of Proposition~\ref{p:below-four-density}. The core estimates~\eqref{e:below-four-band} sum for $q<\frac{4-\varepsilon}{3-\varepsilon}$, which gives $p>4-\varepsilon=3+\alpha$; the admissible $\alpha$ give every $p>p_*(\kappa)=1+\beta_*(\kappa)<4$. Neither the endpoint nor optimality is claimed.}} \label{fig:ka-belowfour} \end{figure} \begin{proposition}[A homogeneous Bellman barrier] \label{p:bellman-barrier} There is $\beta_*=\beta_*(\frac\Lambda\lambda)\in[2,3)$ such that, for every $\beta_*-2<\alpha<1$, there are a function $\Phi\in C^2(\R^2\setminus\{0\})$ and a constant $c>0$ satisfying \refstepcounter{equation} \[ \begin{gathered} \Phi(r^3X,rv)=r^\alpha\Phi(X,v),\qquad \rho(X,v)=(X^2+v^6)^{\frac16},\\ v\Phi_X+b\Phi_{vv}\geq c\rho^{\alpha-2} \qquad\text{for every }b\in[\lambda,\Lambda]. \end{gathered} \tag{\theequation}\label{e:bellman-barrier} \] The function and constants depend only on $\lambda,\Lambda,\alpha$. \end{proposition} We seek one barrier that works for every diffusion coefficient $b\in[\lambda,\Lambda]$. Taking the least favourable coefficient gives the stationary Bellman inequality \[ v\Phi_X+\min_{\lambda\leq b\leq\Lambda}b\Phi_{vv} \geq c\rho^{\alpha-2}. \] Appendix~\ref{s:bellman} proves existence by separation: if no such homogeneous barrier existed, there would be a positive stationary adjoint measure of degree $2+\alpha$. The definition of $\beta_*$ excludes this when $2+\alpha>\beta_*$, and the strict gap $\beta_*<3$ allows us to choose $\alpha<1$. Testing the evolution against this barrier bounds the time spent in a small box in position and velocity. Together with the Harnack comparison and the quadratic velocity test, this controls the mass of visits beginning in each time and position cell. Retaining where visits begin, as well as when, supplies the additional estimate needed below $p=4$. For the rest of the argument we fix one admissible $\alpha$ and set \refstepcounter{equation} \[ \gamma=2-\alpha\in(1,2),\qquad \varepsilon=1-\alpha=\gamma-1>0. \tag{\theequation}\label{e:bellman-exponents} \] Until the end of this section, constants may depend on $\lambda,\Lambda,\alpha$ and on the displayed integrability exponent. The barrier extends continuously to the origin by zero. Homogeneity bounds its derivatives, and integration along the coordinate lines gives \refstepcounter{equation} \[ |\Phi(X,v)-\Phi(X',v')| \leq C\bigl(|X-X'|^{\frac\alpha3}+|v-v'|^\alpha\bigr). \tag{\theequation}\label{e:barrier-modulus} \] Indeed, $|\Phi_X|\leq C|X|^{\frac\alpha3-1}$ and $|\Phi_v|\leq C|v|^{\alpha-1}$ away from the respective zero coordinates, and both singularities are integrable. \begin{lemma}[Occupation and concentration in a phase-space box] \label{l:box-occupation} For fixed $A_0,M>0$, we write \[ \mathcal B_r(Y)=\{(X,v):|X-Y|\leq A_0r^3,\ |v|\leq A_0r\}. \] Uniformly in $X,Y,r,T$, we have \refstepcounter{equation} \[ \int_0^T\mathsf S_s\1_{\mathcal B_r(Y)}(X,v)\dd s \leq Cr^\gamma T^{\frac\alpha2},\qquad |v|\leq M\sqrt T. \tag{\theequation}\label{e:box-occupation} \] Consequently, if $|v|\leq3r$, then \refstepcounter{equation} \[ \mathsf S_t\1_{\mathcal B_r(Y)}(X,v) \leq C\left(1+\frac t{r^2}\right)^{-\frac\gamma2}. \tag{\theequation}\label{e:box-concentration} \] Constants may also depend on $A_0,M,\lambda,\Lambda,\alpha$. \end{lemma} \begin{proof} We justify the use of the singular barrier first. We convolve only in position, with a non-negative smooth approximate identity $\eta_\delta$, and put $\Phi_\delta=\eta_\delta*_X\Phi$. Appendix~\ref{s:bellman} verifies that the derivatives of $\Phi$ have no distribution supported at the origin. The bound \[ |\Phi_{vv}(X,v)|\leq C|X|^{\frac{\alpha-2}{3}} \] is locally integrable in $X$. Thus $\Phi_\delta$ is smooth in position and twice continuously differentiable in velocity, including at $v=0$. Since~\eqref{e:bellman-barrier} holds for every scalar $b$, we may hold $b=a(X,v)$ fixed while convolving the shifted barrier. This gives \[ L[\Phi_\delta(X-Y,v)] \geq c(\eta_\delta*_X\rho^{-\gamma})(X-Y,v). \] Neither a derivative of $a$ nor a transport error appears. For clarity, the moment bounds used to apply the Green identity can be obtained from the smooth diffusion representing the same kernel: \[ \mathrm dX_s=V_s\dd s,\qquad \mathrm dV_s=\sqrt{2a(X_s,V_s)}\dd B_s,\qquad (X_0,V_0)=(X,v), \] where $B_s$ is one-dimensional Brownian motion. Smoothness gives local existence and uniqueness, and bounded diffusion prevents explosion. Its terminal solutions agree with those of Proposition~\ref{p:evolution} by bounded comparison. In particular, \[ \mathbb E\sup_{s\leq T}|V_s-v|^2\leq8\Lambda T, \qquad \sup_{s\leq T}|X_s-X|\leq |v|T+T\sup_{s\leq T}|V_s-v|. \] These estimates also justify stopping on large compact sets and then removing the stops in the Green identity. The modulus~\eqref{e:barrier-modulus}, unchanged by position convolution, makes the endpoint differences uniformly integrable. We next let $\delta\downarrow0$ and use Fatou's lemma to obtain \[ \begin{aligned} c\,\mathbb E\int_0^T\rho(X_s-Y,V_s)^{-\gamma}\dd s &\leq \mathbb E\Phi(X_T-Y,V_T)-\Phi(X-Y,v)\\ &\leq C_M T^{\frac\alpha2}. \end{aligned} \] At the origin the regularised integrand tends to $+\infty$; elsewhere it tends to $\rho^{-\gamma}$. On $\mathcal B_r(Y)$ we have $\rho\leq C_{A_0}r$, which proves~\eqref{e:box-occupation}. For $t\geq9r^2$, Lemma~\ref{l:return-time} and~\eqref{e:box-occupation} give \[ \begin{aligned} \mathsf S_t\1_{\mathcal B_r(Y)}(X,v) &\leq\frac Ct\int_{2t}^{3t}\mathsf S_s\1_{\mathcal B_r(Y)}(X,v)\dd s\\ &\leq Cr^\gamma t^{\frac\alpha2-1} =C\left(\frac r{\sqrt t}\right)^\gamma. \end{aligned} \] For $t<9r^2$, the mass bound gives~\eqref{e:box-concentration}. Harnack is applied to the positive solution $\mathsf S_tF$ in its starting variables, not to an adjoint density. \end{proof} \subsection{Counting visits by time and position} We divide time into intervals of length $r^2$ and position into cells of width $r^3$. The quantity $n_{j,k}$ below counts the mass of visits beginning in one such cell; its sum over $k$ is the earlier time count $m_j$. We will bound the largest cell using the Bellman barrier and the sum using the return-time comparison. We use the intervals~\eqref{e:intervals} with $\bar v=\pm2r$, so $I^*\subset\{|v|<3r\}$. We first start at $e=(0,X_e,v_e)$ with $v_e\in\overline I$, and observe up to $T\leq R^2$. As in Lemma~\ref{l:timed-entrance}, we work in a bounded interval $J$ containing $\overline{I^*}$ in its interior. All constants are independent of $J$. A visit begins when the velocity reaches $\overline I$ and ends when it leaves $I^*$. Another visit is counted only after a return to $\overline I$. The measure $\nu$ is the sum of the distributions of these starting times, positions and velocities, as in Lemma~\ref{l:entrance}. Its total mass counts visits and can exceed one. The existing quadratic test gives \[ \nu(\R^3)\leq C\left(1+\frac Rr\right). \] The sum of the occupations during active visits, and of their terminal measures at any fixed time, is dominated by the full-space evolution. This follows from the alternating Green and exit identities, or from the fact that at most one visit is active along each path. We divide time into $J_0=[0,r^2]$ and $J_j=(jr^2,(j+1)r^2]$ for $j\geq1$, truncating the last interval at $T$. We divide position into $B_k=[kr^3,(k+1)r^3)$, $k\in\mathbb Z$, and put \refstepcounter{equation} \[ n_{j,k}=\nu(J_j\times B_k\times\overline I), \qquad m_j=\sum_k n_{j,k}. \tag{\theequation}\label{e:position-visit-masses} \] Thus $m_j$ counts all starting positions in one time interval, whereas $n_{j,k}$ retains a single position interval. \begin{lemma}[Visits beginning in one time and position interval] \label{l:position-visits} The masses in~\eqref{e:position-visit-masses} satisfy \refstepcounter{equation} \[ m_j\leq C(1+j)^{-\frac12},\qquad \sup_k n_{j,k}\leq C(1+j)^{-\frac\gamma2}. \tag{\theequation}\label{e:position-visit-count} \] In particular, for $10$, which is the stated range. Its upper endpoint is less than $\frac32$, as required in the local estimates. For $|v_0|>8R$, the single-interval argument at the end of Proposition~\ref{p:cylinder-density} already gives the result for every $q<\frac32$. \end{proof} At $q=\frac43$, the extra power in~\eqref{e:below-four-band} is $\frac\varepsilon3>0$. Thus the density belongs to strong $L^{\frac43}$, and the maximum principle includes strong $L^4$ sources. Nearby values of $q$ give $p>4-\varepsilon=3+\alpha$. The barrier provides the additional position bound; the D\"oblin argument still supplies the local density on each one-sign interval. Neither bound alone gives the conclusion. \begin{proof}[Proof of Theorem~\ref{t:a}] We fix $p>p_*=1+\beta_*$ and choose \[ \beta_*-2<\alpha<\min\{1,p-3\} \] in Proposition~\ref{p:bellman-barrier}; this is possible since $\beta_*<3$ and $p>1+\beta_*$. We set $\varepsilon=1-\alpha\in(0,1)$. These exponents depend only on $\frac\Lambda\lambda$ and $p$. Since $p>3+\alpha=4-\varepsilon$, the conjugate $q=\frac p{p-1}$ is in the range of Proposition~\ref{p:below-four-density}. We first assume $a$ smooth and reflect time and position as in Section~\ref{s:maximum}. The function $U(s,X,v)=u(-s,-X,v)$ satisfies $\mathcal K_{\widetilde a}U\geq-g_0$, where $\widetilde a(X,v)=a(-X,v)$ and $g_0(s,X,v)=f_+(-s,-X,v)$. The coefficient remains autonomous with the same bounds. For each $P\in\mathcal Q$, we use $\Gamma_P^{J_Q}|_{\mathcal Q}$. Proposition~\ref{p:strip}\ref{i:strip-potential} gives the source potential, and Proposition~\ref{p:below-four-density} bounds its density. Proposition~\ref{p:abp-abstract} therefore gives \[ \sup_{\overline{\mathcal Q}}U \leq\sup_{\partial^+\mathcal Q}U_+ +C_{\lambda,\Lambda,p}R^{\frac6q-4} \norm{g_0\1_{\{U>0\}}}_{L^p(\mathcal Q)}. \] Since $\frac6q-4=2-\frac6p$, reflection gives~\eqref{e:a-estimate}, including the localised source. Finally, Lemma~\ref{l:borel} applies: we mollify $a$ only in position and velocity, keep $u$ fixed, pass to the limit on compact inner cylinders, and then exhaust the original cylinder. All constants are independent of coefficient derivatives. \end{proof} The threshold of Theorem~\ref{t:a}, supplied by the Bellman construction, is \refstepcounter{equation} \[ p_*(\kappa)=1+\beta_*(\kappa)\in[3,4),\qquad \kappa=\frac\Lambda\lambda. \tag{\theequation}\label{e:bellman-critical-p} \] This definition is nonquantitative: the proof gives neither a numerical formula for $p_*$ nor the estimate at $p=p_*$. Whether $p_*$ is the optimal Aleksandrov threshold remains open. Since $\beta_*<3$, every $p\geq4$ is admissible; with $\alpha=\frac{\beta_*-1}{2}$ one obtains in particular every $p>4-\varepsilon$ with $\varepsilon=\frac{3-\beta_*}{2}\in(0,\frac12]$. \begin{proof}[Proof of Corollary~\ref{c:localised}] For smooth coefficients, both theorem proofs use Proposition~\ref{p:abp-abstract}, whose source is already restricted to $\{U>0\}$. Reflection preserves this set and its source norm. The localised form of Lemma~\ref{l:borel} passes the same estimate to measurable coefficients. \end{proof} \section{From the Aleksandrov principle to H\"older continuity} \label{s:holder} The Aleksandrov estimate first turns positivity on most of a cylinder into a lower bound on a smaller region. A moving barrier carries that lower bound to later points. A covering argument then shows that positivity on any fixed fraction of a cylinder suffices. Applying this to a solution and its negative reduces the oscillation on a smaller cylinder, and iteration gives Hölder continuity. This is the growth argument of Landis and Krylov--Safonov~\cite[Chapter~3, Sections~6 and~8]{Landis}\cite{KrylovSafonov}, adapted to kinetic transport. In this section the localised Aleksandrov estimate is an assumption. The implication is valid for coefficients depending on all three variables; we have proved its hypothesis only for the two coefficient classes of Theorems~\ref{t:A} and~\ref{t:a}. \subsection{Setting} Throughout this section $A:\R^{1+2d}\to\Sym_d$ is Borel measurable and satisfies \refstepcounter{equation} \[ \lambda I_d\leq A(t,x,v)\leq\Lambda I_d \quad\text{for almost every }(t,x,v)\in\R^{1+2d}\,, \tag{\theequation}\label{e:A-full-ellipticity} \] and, extending~\eqref{e:PA}, we write \refstepcounter{equation} \[ \PA u=\partial_tu+v\cdot\nabla_xu-A(t,x,v):D_v^2u\,. \tag{\theequation}\label{e:PA-full} \] The cylinder $Q_R^-(P_0)$ and its kinetic boundary $\partial_{\kin}Q_R^-(P_0)$ are those of~\eqref{e:cylinder} and~\eqref{e:kinetic-boundary}. For $P_0=(t_0,x_0,v_0)$ and $R>0$ let \refstepcounter{equation} \[ \Phi_{P_0,R}(t,x,v) = \bigl(t_0+R^2t,\ x_0+R^3x+R^2tv_0,\ v_0+Rv\bigr)\,, \tag{\theequation}\label{e:kinetic-affine} \] the composition of the kinetic dilation of ratio $R$ with the Galilean translation by $P_0$. It satisfies $\Phi_{P_0,R}\bigl(Q_1^-(0)\bigr)=Q_R^-(P_0)$, $\Phi_{P_0,R}\circ\Phi_{0,\varrho}=\Phi_{P_0,R\varrho}$ and $\abs{\Phi_{P_0,R}(E)}=R^{4d+2}\abs E$ for every Borel set $E$; the exponent $4d+2$ is the homogeneous dimension already met in~\eqref{e:A-estimate}. The kinetic quasi-distance between $P=(t,x,v)$ and $P'=(t',x',v')$ is \refstepcounter{equation} \[ d_{\kin}(P,P') = \abs{t-t'}^{\frac12}+\abs{v-v'} +\max\bigl\{\abs{x'-x-(t'-t)v},\ \abs{x-x'-(t-t')v'}\bigr\}^{\frac13}\,. \tag{\theequation}\label{e:dkin} \] For $m\in\mathbb N$ the \emph{forward stack} of $Q_r^-(P_0)$ is \refstepcounter{equation} \[ Q^{m,+}_r(P_0) = \bigl\{(t,x,v):\ 0u$. Let $O\subseteq\R^{1+2d}$ be open, let $1\leq p<\infty$ and $C_{\mathrm A}<\infty$. A function $u\in C(O)$ is an \emph{admissible supersolution in $O$ with Aleksandrov data $(p,C_{\mathrm A})$} if, for every $P_0\in\R^{1+2d}$ and $R>0$ with $\overline{Q_R^-(P_0)}\subset O$ and every $\psi$ smooth on a neighbourhood of $\overline{Q_R^-(P_0)}$, \refstepcounter{equation} \[ \sup_{\overline{Q_R^-(P_0)}}(\psi-u) \leq \sup_{\partial_{\kin}Q_R^-(P_0)}(\psi-u)_+ +C_{\mathrm A}R^{2-\frac{4d+2}{p}} \bigl\lVert(\PA\psi)_+\1_{\{\psi>u\}}\bigr\rVert_{L^p(Q_R^-(P_0))}\,. \tag{\theequation}\label{e:amp} \] A function $u$ is an \emph{admissible solution in $O$} if both $u$ and $-u$ are admissible supersolutions in $O$. If $u$ is an admissible supersolution, so is $\theta u+c$ for every $\theta>0$ and $c\in\R$, with the same data; taking $\psi=0$ in~\eqref{e:amp} gives the minimum principle for $u$ on every cylinder. Condition~\eqref{e:amp} is precisely the conclusion of Theorem~\ref{t:A}, in the localised form of Corollary~\ref{c:localised}, applied to $\psi-u$ with source $\PA\psi$. \begin{theorem}[The Aleksandrov principle implies H\"older continuity] \label{t:holder} Let $A$ satisfy~\eqref{e:A-full-ellipticity}, let $1\leq p<\infty$ and $C_{\mathrm A}<\infty$. There are $\alpha\in(0,1]$ and $C<\infty$, depending only on $d$, $\lambda$, $\Lambda$, $p$ and $C_{\mathrm A}$, with the following property. Let $\Omega\subseteq\R^{1+2d}$ be open, let $u$ be a bounded admissible solution in $\Omega$ with Aleksandrov data $(p,C_{\mathrm A})$, let $\mathcal C\subset\Omega$ be compact, and let $R_0>0$ be such that $\overline{Q_{2R_0}^-(P)}\subset\Omega$ for every $P\in \mathcal C$. Then \refstepcounter{equation} \[ \abs{u(P)-u(P')} \leq C\Bigl(\operatorname{osc}_{\Omega_{R_0,\mathcal C}}u\Bigr) \Bigl(\frac{d_{\kin}(P,P')}{R_0}\Bigr)^{\alpha} \qquad\text{for all }P,P'\in \mathcal C\,, \tag{\theequation}\label{e:holder} \] where $\Omega_{R_0,\mathcal C}=\bigcup_{P\in \mathcal C}Q^-_{2R_0}(P)$. \end{theorem} The rest of the section proves Theorem~\ref{t:holder} and applies it to the autonomous class. For $A(t,v)$, the final corollary uses the direct local estimate of Proposition~\ref{p:local-A}. \subsection{Positivity from a nearly full superlevel set} We suppose that $u\geq1$ on most of a cylinder. We choose a smooth cutoff $\psi$ between zero and one, equal to one on a smaller interior region. The source in~\eqref{e:amp} is then counted only where $\psi>u$, a subset of the small set $\{u<1\}$. Its norm is small, so $u$ stays close to one wherever $\psi=1$. We use the following fixed interior region, called the \emph{reference cap}: \refstepcounter{equation} \[ \mathcal B = \bigl\{(t,x,v):\ -\tfrac140$, $\ell\geq0$ and $\eta\in(0,1)$, and let $u\geq0$ be an admissible supersolution with Aleksandrov data $(p,C_{\mathrm A})$ in an open set containing $\overline{Q_r^-(P_0)}$. If \refstepcounter{equation} \[ \bigl\lvert\{u\geq\ell\}\cap Q_r^-(P_0)\bigr\rvert \geq (1-\eta)\bigl\lvert Q_r^-(P_0)\bigr\rvert\,, \tag{\theequation}\label{e:near-full-density} \] then \refstepcounter{equation} \[ \inf_{\Phi_{P_0,r}(\mathcal B)}u \geq \bigl(1-C_1\eta^{\frac1p}\bigr)\ell\,. \tag{\theequation}\label{e:near-full} \] \end{lemma} \begin{proof} We choose $\psi_0\in C_c^\infty(Q_1^-(0))$, $0\leq\psi_0\leq1$, with $\psi_0=1$ on $\mathcal B$, and put $\psi=\ell\,\psi_0\circ\Phi_{P_0,r}^{-1}$, extended by zero. The kinetic chain rule gives $|\PA\psi|\leq C_0\ell r^{-2}$, where $C_0$ depends only on $d,\Lambda$ and the fixed cutoff. The boundary term in~\eqref{e:amp} vanishes, and $\{\psi>u\}\subseteq\{u<\ell\}$ has measure at most $\eta r^{4d+2}|Q_1^-(0)|$. Therefore \[ \sup_{\overline{Q_r^-(P_0)}}(\psi-u) \leq C_{\mathrm A}C_0|Q_1^-(0)|^{\frac1p}\eta^{\frac1p}\ell. \] Since $\psi=\ell$ on $\Phi_{P_0,r}(\mathcal B)$, this is~\eqref{e:near-full}. \end{proof} Only one consequence of Lemma~\ref{l:near-full} is used below, at one fixed density. We fix $\eta_*\in(0,1)$ with \refstepcounter{equation} \[ C_1\eta_*^{\frac1p}\leq\tfrac14\,, \tag{\theequation}\label{e:eta-star} \] depending only on $d,\Lambda,p,C_{\mathrm A}$. Lemma~\ref{l:near-full} then gives $\inf u\geq\frac34\ell$ on the cap at density $1-\eta_*$. \subsection{Propagation along kinetic paths} We carry a positive lower bound along a position path with bounded acceleration, using its derivative as the velocity path. A Gaussian barrier centred on this pair expands with velocity width proportional to $h^{\frac12}$ and position width proportional to $h^{\frac32}$. Its differential inequality holds for every admissible coefficient. On each short time interval, we reduce its height so that it starts below the preceding barrier; comparison follows from~\eqref{e:amp}, whose boundary and source terms then vanish. Figure~\ref{fig:ka-propagation} shows the construction. Each comparison cylinder reaches slightly before its time interval. We therefore need the initial lower bound on a tube of positive duration. The ratios in~\eqref{e:propagation-data} describe its widths, the travel time and the acceleration in kinetic units. The propagation constant depends on these ratios and ellipticity, but not on $p$ or $C_{\mathrm A}$. \begin{lemma}[Propagation along kinetic paths] \label{l:propagation} Let $H,\kappa_x,\kappa_v>0$ and $00$ propagates to $P$; velocity is the derivative of position. The cylinders $Q^{(k)}$ have radius $R=\mathsf Wh^{\frac{1}{2}}$, end at the next path point, and are sheared by the velocity there; only some are drawn. The cylinder proportions are schematic: they are drawn with $\mathsf W=2$, whereas~\eqref{e:hstar} gives $\mathsf W=2+\frac32L\lambda^{\frac{1}{2}}+2HT_1^{\frac{1}{2}}$. \textbf{(A$'$)} Each cylinder extends $\mathsf W^2h$ into the past, beyond its block of length $h$. \textbf{(B)} The exact positive region $\{\Psi_k>0\}$ for $\lambda=\Lambda=1$ and $H=0$, where $\Xi_hh=\frac{L^2}{4}$, in path-centred coordinates $\frac{\mathsf D_h^{-1}Z}{L}$, within the bounds~\eqref{e:positivity-set}. \textbf{(C)} Schematic heights as the next barrier starts, $s_k-\frac{h}{128}0\}$, $\Psi_k\leq\ell_ke^{\frac{L^2}{512}}=\ell_{k-1}c_{\mathrm b}\leq\Psi_{k-1}$; both barriers are non-positive outside their positive regions. Repetition gives $u(P)\geq c_{\mathrm{prop}}\ell$.}} \label{fig:ka-propagation} \end{figure} \begin{proof} We construct the Gaussian in coordinates centred on the path. The matrix $\mathsf J$ describes transport from velocity to position, and $\mathsf V$ inserts diffusion in velocity. Their covariance matrix, computed in Step~1, determines the Gaussian. At each step we start with the same rescaled profile and a smaller amplitude, using a time cutoff to make the comparison smooth. Throughout, $Z=(y,V)\in\R^{2d}$, \[ \mathsf J=\begin{pmatrix}0&I_d\\0&0\end{pmatrix}, \qquad \mathsf V=\begin{pmatrix}0\\I_d\end{pmatrix}, \qquad \mathsf D_h=\begin{pmatrix}h^{\frac32}I_d&0\\0&h^{\frac12}I_d\end{pmatrix} \quad(h>0)\,, \] $\varepsilon=\frac{\lambda}{64}$, and $\chi\in C^\infty(\R)$ is a fixed function with $0\leq\chi\leq1$, $\chi=0$ on $\bigl(-\infty,-\frac1{128}\bigr]$ and $\chi=1$ on $[0,\infty)$; no constant below depends on the choice of $\chi$. To obtain smooth barriers, we extend $\bar x$ by its tangent lines at the endpoints and mollify. The acceleration bound $H$ is preserved, and convergence in $C^1$ keeps the tube with halved radii inside~\eqref{e:corridor}. It retains the initial lower bound. Step~5 allows a fixed positive endpoint error. We choose one mollification small enough to preserve the tube inclusion and satisfy both endpoint bounds. We use the halved radii below. \emph{Step 1: the covariance matrix.} For $s\in[-T_0,T_P]$ we use the coordinates \refstepcounter{equation} \[ Z=\bigl(x-\bar x(s),\ v-\dot{\bar x}(s)\bigr)\,, \tag{\theequation}\label{e:moving-frame} \] and, for $\Psi(t,x,v)=\Upsilon(s,Z)$ with $s=t-t_-$, the chain rule gives \refstepcounter{equation} \[ \PA\Psi = \partial_s\Upsilon+(\mathsf JZ)\cdot\nabla_Z\Upsilon -\ddot{\bar x}(s)\cdot\nabla_V\Upsilon-A:D_V^2\Upsilon\,. \tag{\theequation}\label{e:moved-operator} \] We define $\mathsf G_\sigma$ by $\dot{\mathsf G}_\sigma=\mathsf J\mathsf G_\sigma+\mathsf G_\sigma\mathsf J^{\mathsf T}+\lambda\mathsf V\mathsf V^{\mathsf T}$ with $\mathsf G_0=\varepsilon\mathsf D_h^2$, that is \refstepcounter{equation} \[ \mathsf G_\sigma = e^{\sigma\mathsf J}\bigl(\varepsilon\mathsf D_h^2\bigr)e^{\sigma\mathsf J^{\mathsf T}} +\lambda\,\mathsf G_\star(\sigma), \qquad \mathsf G_\star(\theta)=\int_0^\theta e^{\varsigma\mathsf J}\mathsf V\mathsf V^{\mathsf T}e^{\varsigma\mathsf J^{\mathsf T}}\dd\varsigma =\begin{pmatrix}\frac{\theta^3}{3}I_d&\frac{\theta^2}{2}I_d\\\frac{\theta^2}{2}I_d&\theta I_d\end{pmatrix}\,, \tag{\theequation}\label{e:gramian} \] valid for every real $\sigma$. Since $e^{\sigma\mathsf J}=I_{2d}+\sigma\mathsf J$, rescaling gives \refstepcounter{equation} \[ \mathsf G_\sigma=\mathsf D_h\widehat{\mathsf G}_{\frac{\sigma}{h}}\mathsf D_h, \qquad \widehat{\mathsf G}_\theta =\varepsilon\,e^{\theta\mathsf J}e^{\theta\mathsf J^{\mathsf T}}+\lambda\,\mathsf G_\star(\theta)\,; \tag{\theequation}\label{e:gramian-scaling} \] every bound below is therefore a bound on the fixed function $\theta\mapsto\widehat{\mathsf G}_\theta$, uniform in $h$ and in the block. We write $\mathfrak q_\sigma(Z)=Z^{\mathsf T}\mathsf G_\sigma^{-1}Z$, $\mathsf p_\sigma(Z)=\mathsf V^{\mathsf T}\mathsf G_\sigma^{-1}Z$ and $\mathsf R_\sigma=\mathsf V^{\mathsf T}\mathsf G_\sigma^{-1}\mathsf V$. Differentiating $\mathsf G_\sigma^{-1}$ and using the equation for $\mathsf G_\sigma$, \refstepcounter{equation} \[ \bigl(\partial_\sigma+(\mathsf JZ)\cdot\nabla_Z\bigr)\mathfrak q_\sigma=-\lambda\abs{\mathsf p_\sigma}^2, \qquad \nabla_V\mathfrak q_\sigma=2\mathsf p_\sigma, \qquad D_V^2\mathfrak q_\sigma=2\mathsf R_\sigma\,, \tag{\theequation}\label{e:q-identities} \] where the transport terms cancel in the first identity. Four properties of $\widehat{\mathsf G}$ are used: positive definiteness for $\theta\geq-\frac1{128}$, and \refstepcounter{equation} \[ \begin{aligned} &\widehat{\mathsf G}_\theta\leq2\lambda I_{2d}\quad &(\abs\theta\leq1), \\ &\operatorname{tr}\bigl(\mathsf V^{\mathsf T}\widehat{\mathsf G}_\theta^{-1}\mathsf V\bigr)\leq C_{\mathsf G}\coloneqq \frac{128\,d}{\lambda}\quad& (0\leq\theta\leq1), \\ &\widehat{\mathsf G}_\theta\leq\tfrac{13}{50}\,\widehat{\mathsf G}_{\theta+1}\quad &\bigl(-\tfrac1{128}\leq\theta\leq0\bigr)\,. \end{aligned} \tag{\theequation}\label{e:G-bounds} \] These are scalar $2\times2$ calculations. For $|\theta|\leq1$, row sums give $\|\widehat{\mathsf G}_\theta\|\leq3\varepsilon+\frac{11}{6}\lambda<2\lambda$. For $\theta\geq0$, the non-negative covariance term gives $\mathsf V^{\mathsf T}\widehat{\mathsf G}_\theta^{-1}\mathsf V \leq\varepsilon^{-1}(1+\theta^2)I_d$. For $-\frac1{128}\leq\theta\leq0$, the diagonal entries are at least $\frac\varepsilon2$ and the off-diagonal entry has modulus at most $\frac\varepsilon{128}$, so the matrix is positive definite. Finally, $\mathsf G_\star(a)=\mathsf D_a\mathsf G_\star(1)\mathsf D_a$ and $\mathsf G_\star(1)^{-1}\leq16I_{2d}$ give, in this last range, \[ \widehat{\mathsf G}_{\theta+1}\geq\frac{125\lambda}{2048}I_{2d}, \qquad \widehat{\mathsf G}_\theta\leq\frac{65\lambda}{4096}I_{2d}. \] Their quotient is $\frac{13}{50}$. \emph{Step 2: the barrier and its differential inequality.} Let $\ell'\geq0$, $h>0$, $L>0$ and let $s_\bullet$ be the initial time of a block. Writing $\sigma=s-s_\bullet$ and $\mathcal G=e^{-\Xi_h\sigma-\mathfrak q_\sigma(Z)}$, with $Z$ as in~\eqref{e:moving-frame}, put \refstepcounter{equation} \[ \Psi(t,x,v) = \ell'\,\chi\Bigl(\frac{\sigma}{h}\Bigr)\bigl(\mathcal G-e^{-L^2}\bigr), \qquad \Xi_h=\frac{2\Lambda C_{\mathsf G}}{h}+\frac{2H^2}{3\lambda}\,. \tag{\theequation}\label{e:barrier} \] Positive definiteness holds on a neighbourhood of $[-\frac1{128},\infty)$, so the cutoff makes $\Psi$, extended by zero, smooth everywhere. For $0<\sigma\leq h$, the cutoff is one. Equations~\eqref{e:moved-operator} and~\eqref{e:q-identities} give \refstepcounter{equation} \[ \begin{aligned} \PA\Psi &= \ell'\,\mathcal G\Bigl(-\Xi_h+\lambda\abs{\mathsf p_\sigma}^2+2\,\ddot{\bar x}(s)\cdot\mathsf p_\sigma +2A:\mathsf R_\sigma-4\,\mathsf p_\sigma\cdot A\mathsf p_\sigma\Bigr) \\ &\leq \ell'\,\mathcal G\Bigl(-\Xi_h-3\lambda\abs{\mathsf p_\sigma}^2+2H\abs{\mathsf p_\sigma}+\frac{2\Lambda C_{\mathsf G}}{h}\Bigr) \\ &\leq \ell'\,\mathcal G\Bigl(-\Xi_h-\frac{3\lambda}{2}\abs{\mathsf p_\sigma}^2+\frac{2H^2}{3\lambda}+\frac{2\Lambda C_{\mathsf G}}{h}\Bigr) = -\frac{3\lambda}{2}\,\ell'\,\mathcal G\,\abs{\mathsf p_\sigma}^2 \leq0\,, \end{aligned} \tag{\theequation}\label{e:barrier-sign} \] almost everywhere on $\{0<\sigma\leq h\}$. We used $\ell'\mathcal G\geq0$, ellipticity, the trace bound~\eqref{e:G-bounds}, $|\ddot{\bar x}|\leq H$ and Young's inequality. \emph{Step 3: the constants.} We fix, in this order, \refstepcounter{equation} \[ \begin{aligned} &\varepsilon=\frac{\lambda}{64}, \quad C_{\mathsf G}=\frac{128\,d}{\lambda}, \quad \Xi_*=2\Lambda C_{\mathsf G}+\frac{2H^2T_1}{3\lambda}, \quad L=2\,\Xi_*^{\frac12}, \quad \mathsf W=2+\tfrac32L\lambda^{\frac12}+2HT_1^{\frac12}, \\ &c_1=\mathsf W+H\mathsf W^2T_1^{\frac12}, \quad c_2=\mathsf W^3+\tfrac12H\mathsf W^4T_1^{\frac12}, \quad h_*=\min\biggl\{\frac{T_0}{\mathsf W^2},\ \Bigl(\frac{\kappa_v}{2c_1}\Bigr)^{2},\ \Bigl(\frac{\kappa_x}{2c_2}\Bigr)^{\frac23}\biggr\}\,, \end{aligned} \tag{\theequation}\label{e:hstar} \] and put \refstepcounter{equation} \[ c_{\mathrm b}=e^{-\frac34L^2}-e^{-L^2}>0, \qquad \gamma=c_{\mathrm b}\,e^{-\frac{L^2}{512}}\in(0,1)\,. \tag{\theequation}\label{e:block-parameters} \] Thus $\mathsf W\geq2$ and $\mathsf W^2h_*\leq T_0$; and for $00\}$ one has $\chi>0$ and $\mathcal G>e^{-L^2}$, hence $\sigma>-\frac1{128}h$ and $\mathfrak q_\sigma(Z)0\}\cap\overline Q\,. \tag{\theequation}\label{e:positivity-set} \] The acceleration bound places every point of $\overline Q$ within $R+H\mathsf W^2h\leq c_1h^{\frac12}$ in velocity and $R^3+\frac12H\mathsf W^4h^2\leq c_2h^{\frac32}$ in position of the path at the same time. Since $\mathsf W^2h\leq T_0$, the whole cylinder lies in the tube. On $\{\Psi>0\}$, moreover, $|\sigma-h|\leq2h$. Thus~\eqref{e:positivity-set} gives \[ \abs{v-\dot{\bar x}(s_\bullet+h)} \leq\left(\frac32L\lambda^{\frac12}+2HT_1^{\frac12}\right)h^{\frac12}u\}\subseteq\{0<\sigma\leq h\}$, where~\eqref{e:barrier-sign} holds. Both terms on the right of~\eqref{e:amp} vanish, proving $\Psi\leq u$ on $\overline Q$. \emph{Step 5: iteration along the path.} We set $n=\lceil\frac{T_P}{h_*}\rceil$, $h=\frac{T_P}n$, $s_k=kh$, $\ell_0=e^{-\frac{L^2}{512}}\ell$ and $\ell_k=\gamma\ell_{k-1}$. We denote by $\Psi_k$ and $Q^{(k)}$ the barrier and cylinder starting at $s_k$ with amplitude $\ell_k$, for $0\leq k0$ and $\sigma=s-s_k\leq0$. There $-\frac h{128}<\sigma\leq0$, so $\mathfrak q_\sigma\geq0$ and $\Xi_hh\leq\frac14L^2$ give \refstepcounter{equation} \[ \Psi_k\leq\ell_k\,\mathcal G\leq\ell_k\,e^{\Xi_h\frac{h}{128}}\leq\ell_k\,e^{\frac{L^2}{512}}\,. \tag{\theequation}\label{e:cutoff-bound} \] For $k=0$, this set lies in the initial tube by~\eqref{e:positivity-set} and~\eqref{e:hstar}. There $u\geq\ell=\ell_0e^{\frac{L^2}{512}}$, so~\eqref{e:cutoff-bound} gives the comparison. For $k\geq1$ that set lies in $\overline{Q^{(k-1)}}$: its times $s\in[s_k-\frac1{128}h,s_k]$ lie in $[s_k-\mathsf W^2h,s_k]$, and $\abs{s-s_k}\leq h$ turns~\eqref{e:positivity-set} into \[ \abs{v-\dot{\bar x}(s_k)}\leq\left(\frac32L\lambda^{\frac12}+HT_1^{\frac12}\right)h^{\frac12}0$, depending only on $d$, with the following property. Let $m\in\mathbb N$ with $m\geq1$, let $\eta,r_0\in(0,1)$, and let $E$ and $F$ be bounded measurable sets with $E\subseteq F\cap Q_1^-(0)$. We assume that every kinetic cylinder $Q_r^-(P)\subseteq Q_1^-(0)$ with \[ \bigl\lvert E\cap Q_r^-(P)\bigr\rvert\geq(1-\eta)\bigl\lvert Q_r^-(P)\bigr\rvert \] satisfies $r0$ and $\ell\geq0$, and let $u\geq0$ be an admissible supersolution with Aleksandrov data $(p,C_{\mathrm A})$ in an open set containing $\overline{\Phi_{P_0,r}(\mathcal T_*)}$. If \[ \bigl\lvert\{u\geq\ell\}\cap Q_r^-(P_0)\bigr\rvert \geq (1-\eta_*)\bigl\lvert Q^-_r(P_0)\bigr\rvert\,, \] then $u\geq\delta\ell$ on $Q^{m,+}_r(P_0)$. \end{lemma} \begin{proof} We assume $\ell>0$. Applying Lemma~\ref{l:near-full} directly on $Q_r^-(P_0)$, with the choice~\eqref{e:eta-star}, gives \refstepcounter{equation} \[ u\geq\tfrac34\ell \qquad\text{on }\Phi_{P_0,r}(\mathcal B)\,. \tag{\theequation}\label{e:cap-bound} \] We start every path at the image $P_-=(t_-,x_-,v_-)$ of $(-\frac3{16},0,0)$ under $\Phi_{P_0,r}$: \refstepcounter{equation} \[ t_-=t_0-\tfrac3{16}r^2, \qquad v_-=v_0, \qquad x_-=x_0-(t_0-t_-)v_0\,, \tag{\theequation}\label{e:cap-point} \] with zero acceleration before $t_-$. We use \refstepcounter{equation} \[ T_0=\tfrac1{32}r^2, \qquad \kappa_x=\tfrac12\,8^{-3}r^3, \qquad \kappa_v=\tfrac1{16}r\,. \tag{\theequation}\label{e:cap-box} \] In the coordinates of the cap, this initial tube has times in $[-\frac7{32},-\frac3{16}]$, position radius $\frac12 8^{-3}$ and velocity radius $\frac1{16}$. It lies inside~\eqref{e:cap}, so~\eqref{e:cap-bound} supplies its lower bound. For $P\in Q^{m,+}_r(P_0)$, the travel time $T_P=t_P-t_-$ lies between $\frac18r^2$ and $T_1=(m+1)r^2$, and exceeds $T_0$. We use the cubic Hermite path~\eqref{e:hermite} from $(x_-,v_-)$ to $(x_P,v_P)$. Its position mismatch is $x_P-x_--T_Pv_-=x_P-x_0-(t_P-t_0)v_0$, bounded by $(m+2)r^3$ by~\eqref{e:stack}; the velocity mismatch is at most $r$. Hence its acceleration is bounded by \[ H=\bigl(384(m+2)+32\bigr)r^{-1}. \] Integrating this bound, including the zero-acceleration extension before $t_-$, gives a constant $c_m\geq m+3$ such that \refstepcounter{equation} \[ \abs{\dot{\bar x}(s)-v_0}\leq c_mr, \qquad \bigl\lvert\bar x(s)-x_0-(t_-+s-t_0)v_0\bigr\rvert\leq c_mr^3, \qquad -T_0\leq s\leq T_P\,, \tag{\theequation}\label{e:path-bounds} \] with $c_m$ depending only on $m$. We take $\mathcal T_*=\{(t,x,v):-2t_0-2r^2$, place every tube inside $\Phi_{P_0,r}(\mathcal T_*)$. Lemma~\ref{l:propagation} now gives $u(P)\geq\frac34c_{\mathrm{prop}}\ell$. Its four ratios~\eqref{e:propagation-data} depend only on $m$, uniformly in $P,P_0,r$. Thus $\delta=\frac34c_{\mathrm{prop}}$ has the asserted dependence. \end{proof} \subsection{Positivity from every positive density} We no longer assume that the positive set fills most of a cylinder. Starting from any fixed positive fraction, we lower the level from $1$ to $\delta$, then to $\delta^2$, and so on. At each level, a cylinder on which the set is nearly full gives positivity on its forward stack by Lemma~\ref{l:visibility}. The covering theorem then forces the measure of the next superlevel set to grow, up to an error controlled by a small radius threshold. After finitely many steps there must be a nearly full cylinder of radius above that threshold. Its cap can be joined to every point of the target cylinder by Lemma~\ref{l:propagation}. In the proof we choose the number of levels from the initial fraction $\beta$, then choose the radius threshold small enough to control the accumulated covering error. This order is important: propagation from a smaller cap gives a weaker lower bound, but the chosen cap still has a positive lower bound on its radius. The sampling and target cylinders are fixed inside the original cylinder, with the target sharing its top time, so the final estimate can be iterated on smaller cylinders. \begin{theorem}[Growth lemma] \label{t:all-density} There are $\vartheta\in(0,1)$ and a kinetic cylinder $\Sigma$ with $\overline\Sigma\subset Q_1^-(0)$, depending only on $d$, $\lambda$, $\Lambda$, $p$ and $C_{\mathrm A}$, and, for every $\beta\in(0,1]$, a constant $\kappa_\beta>0$ depending only on $\beta$ and the same parameters, with the following property. Let $P_0\in\R^{1+2d}$ and $R>0$, and let $u\geq0$ be an admissible supersolution with Aleksandrov data $(p,C_{\mathrm A})$ in an open set containing $\overline{Q_R^-(P_0)}$. If \refstepcounter{equation} \[ \bigl\lvert\{u\geq1\}\cap\Phi_{P_0,R}(\Sigma)\bigr\rvert \geq \beta\bigl\lvert\Phi_{P_0,R}(\Sigma)\bigr\rvert\,, \tag{\theequation}\label{e:all-density-hypothesis} \] then \refstepcounter{equation} \[ u\geq\kappa_\beta \quad\text{on }Q^-_{\vartheta R}(P_0)\,. \tag{\theequation}\label{e:all-density} \] \end{theorem} \begin{proof} \emph{Step 1: the reference geometry.} Relative time, velocity and transported position scale under $\Phi_{P_0,s}$ as $s^2,s,s^3$. Thus \refstepcounter{equation} \[ \Phi_{P_0,s}\bigl(Q_1^-(P')\bigr)=Q^-_s\bigl(\Phi_{P_0,s}(P')\bigr), \qquad \Phi_{P_0,s}\circ\Phi_{P',r}=\Phi_{\Phi_{P_0,s}(P'),\,sr}\,, \tag{\theequation}\label{e:phi-composition} \] and the same scaling carries forward stacks to forward stacks. We put $P^\sharp=(-(m+2),0,0)$ and $\Sigma^\sharp=Q_1^-(P^\sharp)$, whose times fill $\bigl(-(m+3),-(m+2)\bigr)$. A kinetic cylinder $Q^-_r(P)\subseteq\Sigma^\sharp$ has $r\leq1$ and top time at most $-(m+2)$, so its forward stack lies in $\{t\leq-2\}$, whereas $Q_1^-(0)$ lies in $\{-10$ small enough that $\Phi_{0,\varepsilon_0}(\mathcal G)\subseteq Q_1^-(0)$ and $\overline{\Phi_{0,\varepsilon_0}(\Sigma^\sharp)}\subset Q_1^-(0)$. We set $\Sigma=\Phi_{0,\varepsilon_0}(\Sigma^\sharp)$, $\vartheta=\varepsilon_0$ and $\Phi_\star=\Phi_{P_0,R\varepsilon_0}$. These choices depend only on the stated parameters, and~\eqref{e:phi-composition} gives \refstepcounter{equation} \[ \begin{aligned} &\Phi_\star(\mathcal G)\subseteq Q_R^-(P_0), \\ &\Phi_\star\bigl(\Sigma^\sharp\bigr)=\Phi_{P_0,R}(\Sigma)\eqqcolon \Sigma_R =Q^-_{R\varepsilon_0}\bigl(\Phi_\star(P^\sharp)\bigr), \\ &\Phi_\star\bigl(Q_1^-(0)\bigr)=Q^-_{\vartheta R}(P_0)\,. \end{aligned} \tag{\theequation}\label{e:transported-geometry} \] In particular $u\geq0$ is an admissible supersolution in an open set containing $\overline{\Phi_\star(\mathcal G)}$. \emph{Step 2: a dense cylinder whose radius is bounded below.} Let $r_0\in(0,1)$ and $j\geq0$, and suppose that some kinetic cylinder $Q^-_r(P)\subseteq\Sigma_R$ satisfies \refstepcounter{equation} \[ r\geq r_0R\varepsilon_0, \qquad \bigl\lvert\{u\geq\delta^j\}\cap Q^-_r(P)\bigr\rvert\geq(1-\eta_*)\bigl\lvert Q_r^-(P)\bigr\rvert\,. \tag{\theequation}\label{e:macroscopic} \] The cap construction in Lemma~\ref{l:visibility} gives $P_-\in\Sigma_R$ and $u\geq\frac34\delta^j$ on an initial tube with the data~\eqref{e:cap-box}. We join it to any $P'\in Q^-_{\vartheta R}(P_0)$ by the cubic Hermite path, continued with zero acceleration before its start. In the coordinates of $\Phi_\star$, the travel time lies in $[m+1,m+3]$ and the acceleration is at most $6(m+5)+8$. The radius becomes $\frac r{R\varepsilon_0}\in[r_0,1]$, so the initial duration and the two tube widths in~\eqref{e:cap-box} have positive lower bounds depending only on $r_0$. These tubes lie in $\Phi_\star(\mathcal G)$ by~\ref{i:G-tube}; the affine map preserves the path class and the ratios~\eqref{e:propagation-data}. Using the largest travel time and acceleration and the smallest duration and widths, Lemma~\ref{l:propagation} gives one constant $c_{\mathrm{mac}}>0$, depending only on $d,\lambda,\Lambda,m,r_0$, such that \refstepcounter{equation} \[ u\geq\tfrac34c_{\mathrm{mac}}\delta^j \qquad\text{on }Q^-_{\vartheta R}(P_0)\,. \tag{\theequation}\label{e:macroscopic-conclusion} \] \emph{Step 3: iteration over the levels.} Let $\beta\in(0,1]$ and let $u$ satisfy~\eqref{e:all-density-hypothesis}. We choose $N\in\mathbb N$ with $a^{-N}\beta>1-\eta_*$, which is possible because $a<1$ by~\eqref{e:m-star}, and then $r_0\in(0,1)$ so small that \refstepcounter{equation} \[ L_0\sum_{k=0}^{N-1}a^{-k}1-\eta_*\,. \] The cylinder $\Sigma_R$ itself, of radius $R\varepsilon_0\geq r_0R\varepsilon_0$, therefore satisfies~\eqref{e:macroscopic} at the level $N$, and~\eqref{e:macroscopic-conclusion} gives $u\geq\frac34c_{\mathrm{mac}}\delta^N$ on $Q^-_{\vartheta R}(P_0)$. Both alternatives thus yield~\eqref{e:all-density} with $\kappa_\beta=\frac34c_{\mathrm{mac}}\delta^N>0$, a constant depending only on $\beta$, $d$, $\lambda$, $\Lambda$, $p$ and $C_{\mathrm A}$. \end{proof} \subsection{Oscillation contraction} We normalise the solution between $-1$ and $1$. At least one of the sets $\{u\leq0\}$ and $\{u\geq0\}$ occupies half the sampling cylinder. In the first case the growth lemma applied to $1-u$ lowers the upper bound on the target; in the second, applied to $1+u$, it raises the lower bound. Either way the oscillation decreases by a fixed factor. Only affine changes of the solution are used, so no truncation or extension by zero is needed. \begin{proposition}[Oscillation contraction] \label{p:oscillation} Let $\vartheta$ be as in Theorem~\ref{t:all-density}. There is $q\in(0,1)$, depending only on $d$, $\lambda$, $\Lambda$, $p$ and $C_{\mathrm A}$, such that, for every $P_0\in\R^{1+2d}$ and $R>0$, every admissible solution $u$ with Aleksandrov data $(p,C_{\mathrm A})$ in an open set containing $\overline{Q^-_R(P_0)}$ satisfies \refstepcounter{equation} \[ \operatorname{osc}_{Q^-_{\vartheta R}(P_0)}u \leq q\operatorname{osc}_{Q^-_R(P_0)}u\,. \tag{\theequation}\label{e:oscillation} \] \end{proposition} \begin{proof} We may assume $\kappa_{\frac12}\leq1$ and normalise $u$ so that $-1\leq u\leq1$ and its oscillation on $Q_R^-(P_0)$ is two; the constant case is immediate. At least one of $\{u\leq0\}$ and $\{u\geq0\}$ fills half the sampling cylinder. By replacing $u$ with $-u$, we may assume it is $\{u\leq0\}$. If $O$ is the open set of admissibility, then for every $\epsilon>0$ the function $1-\frac u{1+\epsilon}$ is a positive admissible supersolution on $O\cap\{u<1+\epsilon\}$, which contains $\overline{Q_R^-(P_0)}$. It is at least one on $\{u\leq0\}$. Theorem~\ref{t:all-density} gives \[ u\leq(1+\epsilon)(1-\kappa_{\frac12}) \quad\text{on }Q^-_{\vartheta R}(P_0). \] Letting $\epsilon\downarrow0$ reduces the oscillation to at most $2-\kappa_{\frac12}$. Undoing the normalisation proves~\eqref{e:oscillation} with $q=1-\frac12\kappa_{\frac12}$. \end{proof} \begin{proof}[Proof of Theorem~\ref{t:holder}] We use $\vartheta,q$ from Proposition~\ref{p:oscillation} and put $\alpha=\min\{1,\frac{\log\frac1q}{\log\frac1\vartheta}\}$. Iterating at the same centre gives $\operatorname{osc}_{Q^-_{\vartheta^kR}(P_1)}u \leq q^k\operatorname{osc}_{Q_R^-(P_1)}u$. For $\vartheta^{k+1}R<\rho\leq\vartheta^kR$, inclusion of the cylinders and $q\leq\vartheta^\alpha$ therefore imply \refstepcounter{equation} \[ \operatorname{osc}_{Q^-_\rho(P_1)}u \leq\vartheta^{-\alpha}\left(\frac\rho R\right)^\alpha \operatorname{osc}_{Q_R^-(P_1)}u. \tag{\theequation}\label{e:oscillation-decay} \] For $P,P'\in \mathcal C$, we order the times so that $P$ is earlier and put $\rho_*=d_{\kin}(P,P')$. If $\rho_*>0$, both points lie in $\overline{Q^-_{\rho_*}(P')}$ by~\eqref{e:dkin}. The set \[ \overline{\Omega_{R_0,\mathcal C}} =\bigcup_{P''\in \mathcal C}\overline{Q^-_{2R_0}(P'')} \] is compact as the image of $\mathcal C\times\overline{Q_1^-(0)}$ under the kinetic affine map, and is contained in $\Omega$. For $\rho_*u\}$ on which the source is counted, so that every $p<\infty$ serves and Theorem~\ref{t:holder} needs no lower bound on $p$. A bound by the supremum of the source alone, having no gain from the measure of that set, does not start the argument. \end{remark} \subsection{The two coefficient classes} For the classical solution class used here, the maximum principles apply on every compact inner cylinder. \begin{corollary}[H\"older continuity for coefficients depending on time and velocity] \label{c:holder-A} Let $A$ satisfy~\eqref{e:A-ellipticity} and let $\Omega\subseteq\R^{1+2d}$ be open. Every bounded $u\in C^{1,1,2}(\Omega)$ solving $\PA u=0$ almost everywhere in $\Omega$ satisfies~\eqref{e:holder}, with $C$ and $\alpha$ depending only on $d$, $\lambda$ and $\Lambda$. \end{corollary} \begin{proof} We apply Proposition~\ref{p:local-A}. For $P,P'\in \mathcal C$ with $d_{\kin}(P,P')<\frac14R_0$, centre the cylinder at whichever point has the later time, say $P$. Both points lie in the closure of its inner half-cylinder. The spatial increment in Proposition~\ref{p:local-A} is one of those defining $d_{\kin}$, and continuity permits the point on the top face. Its outer cylinder lies in $\Omega_{R_0,\mathcal C}$, so the proposition gives~\eqref{e:holder}. For pairs at distance at least $\frac14R_0$, the oscillation bound gives the same conclusion after increasing $C$. \end{proof} \begin{corollary}[H\"older continuity for autonomous coefficients in one space dimension] \label{c:holder-a} Let $d=1$, let $a:\R\times\R\to[\lambda,\Lambda]$ be Borel measurable and let $\Omega\subseteq\R^3$ be open. Every bounded $u\in C^{1,1,2}(\Omega)$ solving $\Pa u=0$ almost everywhere in $\Omega$ satisfies~\eqref{e:holder}, with $C$ and $\alpha$ depending only on $\lambda$ and $\Lambda$. \end{corollary} \begin{proof} For every smooth test $\psi$, apply Theorem~\ref{t:a} and Corollary~\ref{c:localised}, with $p=6$, to $\psi-u$ and $\psi+u$ on each compactly contained kinetic cylinder. Their sources are $\Pa\psi$, which is bounded and Borel. Thus $u$ and $-u$ satisfy~\eqref{e:amp}, and Theorem~\ref{t:holder} applies. \end{proof} For a Borel coefficient these remain a priori estimates for solutions of the stated classical regularity; no existence assertion is made. \appendix \section{Construction of the solution operators} \label{s:construction} This appendix proves the three statements that were used but not proved in the body: the terminal solution operator and its kernels (Proposition~\ref{p:evolution}), the Duhamel potential (Lemma~\ref{l:duhamel}), and the Green and exit measures on a velocity strip (Proposition~\ref{p:strip}). \begin{theorem}[H\"ormander~{\cite[Theorem~1.1]{Hormander1967}}] \label{t:hormander} Let $X_0,X_1,\dots,X_k$ be smooth real vector fields on an open set $\mathcal O\subseteq\R^N$ such that the Lie algebra generated by $X_0,\dots,X_k$ spans $\R^N$ at every point of $\mathcal O$. Then every distribution $u$ on $\mathcal O$ with $(X_0+\sum_{i=1}^kX_i^2+c)u\in C^\infty(\mathcal O)$, $c\in C^\infty(\mathcal O)$, belongs to $C^\infty(\mathcal O)$; here $c$ is real-valued. \end{theorem} For our operator, the positive square root $\beta=B^{\frac12}$ is smooth. We put \[ X_i=\sum_j\beta_{ij}\partial_{v_j},\qquad X_0=\partial_\sigma+b(v)\cdot\nabla_z-\sum_k c_k\partial_{v_k}, \qquad c_k=\sum_{i,j}\beta_{ij}\partial_{v_j}\beta_{ik}. \] Then $\Lop=X_0+\sum_iX_i^2$. The $X_i$ span the velocity directions, and the $z$-components of $[X_i,X_0]$ are $Db(v)\beta_i$, where $\beta_i$ is the $i$th column of $\beta$. These span the $z$-directions by~\eqref{e:b-bounds}. Modulo these directions, $X_0$ gives $\partial_\sigma$. Thus Theorem~\ref{t:hormander} applies, also when $B$ depends on $z$. \begin{proof}[Proof of Proposition~\ref{p:evolution}] We fix a finite time interval $[\sigma_-,\tau]$ and $F\in C_c^\infty(\Omega_\tau\times\R^d)$. We write $\Lop_\varepsilon=\Lop+\varepsilon\Delta_z$, $0\leq\varepsilon\leq1$. We construct its terminal solution and then let $\varepsilon\downarrow0$. \emph{Comparison and boundary barriers.} Straightening the moving domain by $Y=v-\gamma(\sigma)$ adds the drift $-\gamma'(\sigma)\cdot\nabla_Y$. On a bounded cylinder the usual maximum argument, applied to $u+\delta\sigma$, gives comparison for $\Lop_\varepsilon$, including $\varepsilon=0$. We apply it successively across the finitely many smooth pieces of $\gamma$. Since $b$ has at most linear growth, for $C$ sufficiently large \[ \Phi=e^{C(\tau-\sigma)}(1+|Y|^2+|z|^2), \qquad \Lop_\varepsilon\Phi\leq-\Phi \quad(0\leq\varepsilon\leq1). \] Here $C$ depends only on bounds for $\gamma,\gamma'$, $b(0)$, $L_b,d,\Lambda$ on the fixed interval. Comparing $u-\delta\Phi$ on bounded truncations and then letting $\delta\downarrow0$ proves comparison for bounded functions on the full moving domain. We suppose $\Omega=B_{r_0}$ and put $\mathsf d=r_0-|v-\gamma(\sigma)|$. With $L_\gamma$ a Lipschitz constant for $\gamma$, we set \[ \kappa=\frac{1+L_\gamma}{\lambda}, \qquad w(t)=\frac{1-e^{-\kappa t}}{\kappa}. \] In a boundary collar, with $e=\frac{v-\gamma}{|v-\gamma|}$, \[ \Lop_\varepsilon w(\mathsf d) =w''e\cdot Be-\frac{w'}{|v-\gamma|} \bigl(\operatorname{tr}B-e\cdot Be\bigr)+w'\gamma'\cdot e \leq w'(L_\gamma-\lambda\kappa)\leq0. \] We choose $\mathsf d_*>0$ smaller than $\frac14r_0$ and half the distance of $\supp F$ from the terminal moving velocity boundary. The functions $\|F\|_\infty\frac{w(\mathsf d)}{w(\mathsf d_*)}$ and $\pm F+(\tau-\sigma)M_F$, where $M_F=\sup_{0\leq\varepsilon\leq1}\|\Lop_\varepsilon F\|_\infty$, are lateral and terminal barriers. Comparison gives \refstepcounter{equation} \[ |u_\varepsilon|\leq\|F\|_\infty \min\left\{1,\frac{w(\mathsf d)}{w(\mathsf d_*)}\right\}, \qquad |u_\varepsilon-F|\leq(\tau-\sigma)M_F \quad\text{near }\sigma=\tau. \tag{\theequation}\label{e:two-barriers} \] Compact support of $F$, its positive terminal-boundary distance and the common Lipschitz bound on the approximating curves give a terminal neighbourhood uniform in the exhaustion and in $\varepsilon$. For $\Omega=\R^d$ there is no lateral barrier and the terminal bound holds throughout the interval. \emph{Construction.} For $\varepsilon>0$, we approximate $\gamma$ uniformly by smooth curves $\gamma_k$ with the same Lipschitz bound. We choose $\alpha_k\downarrow0$ with $\alpha_{k+1}\leq\frac14\alpha_k$ and $\|\gamma_k-\gamma\|_\infty\leq\frac12\alpha_k$. The balls $B_{r_0-\alpha_k}(\gamma_k(\sigma))$ then increase to $\Omega_\sigma$. In the straightened variables $V=v-\gamma_k(\sigma)$, we use the ellipsoids \[ \frac{|V|^2}{(r_0-\alpha_k)^2}+\frac{|z|^2}{L^2}<1 \] and let $L\to\infty$. For $\Omega=\R^d$, we use balls in $\R^{2d}$. These are fixed smooth bounded spatial domains. We take zero lateral data and the transformed terminal datum $F$, choosing $k,L$ large enough that its compact support lies strictly inside the terminal domain. The datum and all its derivatives then vanish near the terminal--lateral corner. After time reflection, the coefficients are smooth and uniformly parabolic for each fixed $\varepsilon>0$, and the data satisfy the boundary compatibility conditions. Classical Dirichlet solvability~\cite[Theorem~5.14]{Lieberman} therefore applies. Its regularity constants may depend on the truncation and on $\varepsilon$; the uniform bounds used in the limit come from our comparison barriers. The barriers are uniform along the exhaustion. For fixed $\delta>0$, we choose $m$ with $\frac32\alpha_m<\delta$ and compare on $B_{r_0-\delta}(\gamma(\sigma))\times\{|z|0$, put $Y_j=\sqrt\varepsilon\,\partial_{z_j}$. Then $\Lop_\varepsilon=X_0+\sum_iX_i^2+\sum_jY_j^2$, and $X_0,X_1,\dots,X_d,Y_1,\dots,Y_d$ already span $\R^{1+2d}$ at every point, so Theorem~\ref{t:hormander} applies to $u_\varepsilon$ without brackets. We now take a weak-$*$ limit in $L^\infty$ as $\varepsilon\downarrow0$. For every compactly supported smooth test function, \refstepcounter{equation} \[ \int u_\varepsilon(\Lop^*\psi+\varepsilon\Delta_z\psi)=0, \qquad \Lop^*\psi=-\partial_\sigma\psi +\sum_{i,j}\partial_{v_iv_j}(B_{ij}\psi)-b(v)\cdot\nabla_z\psi. \tag{\theequation}\label{e:adjoint} \] The last term carrying $\varepsilon$ tends to zero, so the limit solves $\Lop u=0$ and is smooth by Theorem~\ref{t:hormander}. The bounds~\eqref{e:two-barriers} pass first almost everywhere, then everywhere by continuity. They give the required terminal and lateral traces, including the corners. Positivity and contraction pass from the approximations. Growth comparison gives uniqueness, independence of the subsequence, and independence of $\sigma_-$. This proves~\ref{i:classical}. \emph{Kernels and composition.} At a fixed starting point, $F\mapsto u(\sigma,v,z)$ is a positive functional on $C_c(\Omega_\tau\times\R^d)$ of norm at most one. Riesz representation gives $K_{\sigma,\tau}$; we use its integral to define $\mathsf S_{\sigma,\tau}$ on bounded Borel data. For $\sigma0$ on $(0,\infty)$. For $\delta>0$ put $\rho_\pm=\pm(X-X_e)-(V_H+\delta)(s-s_e)$. Then $\Ka\chi(\rho_\pm)\leq-\delta\chi'(\rho_\pm)$ and $\chi(\rho_\pm(e))=0$. The Green identity is a sum of two non-negative integrals, so both measures vanish on $\{\rho_\pm>0\}$. We let $\delta\downarrow0$ to obtain the cone. If the velocity has sign $\varsigma$, use also $\rho=-\varsigma(X-X_e)+(v_H-\delta)(s-s_e)$ to obtain the lower cone bound. There is no exit mass on the initial velocity face. To see this, let $m_H$ and $\ell$ be the midpoint and half-length of $H$, and $d_0=\ell-|v_e-m_H|>0$. The function \[ \Upsilon_{\theta,\kappa}(s,v) =e^{\Lambda\kappa^2(s_-+\theta-s)} \frac{\cosh(\kappa(v-m_H))}{\cosh(\kappa\ell)} \] satisfies $\Ka\Upsilon_{\theta,\kappa}\leq0$ and is at least one on $\{v\in\partial H,\ s\leq s_-+\theta\}$. The Green identity therefore gives \[ \omega_e^H\{v\in\partial H,\ s\leq s_-+\theta\} \leq2e^{\Lambda\kappa^2\theta-\kappa d_0} =2e^{-\frac{d_0^2}{4\Lambda\theta}}, \qquad \kappa=\frac{d_0}{2\Lambda\theta}. \] Letting $\theta\downarrow0$ shows that $\omega_e^H$ is carried by $\partial^\dagger\mathscr S$. Uniqueness of bounded homogeneous solutions with prescribed exit values now also proves its independence of $s_-$. For later use we extend~\eqref{e:strip-identity-c} to every bounded $\phi\in C^{1,1,2}(\mathscr S)\cap C(\mathscr S\cup\partial^\dagger\mathscr S)$ with bounded $\Ka\phi$, at an interior starting point. The function $u=\phi-W_{-\Ka\phi}$ is a bounded smooth homogeneous solution with exit value $\phi$. To approximate its boundary values, take a smooth increasing $\chi_1$, zero below $0$ and one above $1$, and put \[ \Theta_j=\sum_{\pm} \chi_1\bigl(\pm(X-X_e)-(V_H+1)(s-s_e)-j\bigr). \] Then $\Ka\Theta_j\leq0$, $\Theta_j(e)=0$, and $\Theta_j\geq1$ when $|X-X_e|\geq j+1+(V_H+1)(s_+-s_-)$. We choose a smooth compactly supported position cutoff $0\leq\rho_j\leq1$, equal to one when $|X-X_e|\leq j+1+(V_H+1)(s_+-s_-)$, and a smooth time cutoff $0\leq\eta_j\leq1$, zero for $s\leq s_-+\frac1{2j}$ and one for $s\geq s_-+\frac1j$. For large $j$, the boundary datum $h_j=\rho_j\eta_j\phi$, extended by zero at the initial velocity face, belongs to $C_c(\mathcal E)$ and converges to $\phi$ on $\partial^\dagger\mathscr S$. On $\partial^\dagger\mathscr S$, its discarded part is bounded by $\|\phi\|_\infty(\Theta_j+\Upsilon_{\frac{1}{j},\kappa_j})$, where $\kappa_j=\frac{jd_0}{2\Lambda}$: the first barrier controls the position cutoff, and the second is at least one on the early velocity face. Growth comparison yields \[ |u(e)-\mathsf Eh_j(e)| \leq2\|\phi\|_\infty e^{-\frac{jd_0^2}{4\Lambda}}\longrightarrow0. \] Dominated convergence against $\omega_e^H$ proves the Green identity for this class. The class in~\ref{i:green-identity} also permits a starting point with $s_e=s_-$: we first lower $s_-$ and use the independence just proved. \emph{Monotonicity and decomposition.} Formula~\eqref{e:strip-green} and domain comparison give monotonicity in $s_+$ and $H$. Integrating in the starting measure gives~\ref{i:monotone-strip}, and Lemma~\ref{l:duhamel} gives~\ref{i:strip-potential}. For $H_1\subseteq H_2$, we take $F\in C_c^\infty(\mathscr S_{H_2})$ and $h$ continuous and compactly supported on its closed exit boundary. The function \[ W(e')=\int F\dd\Gamma_{e'}^{H_2}+\int h\dd\omega_{e'}^{H_2} \] is bounded and smooth, solves $\Ka W=-F$, and has exit value $h$. Its restriction to the smaller strip belongs to the enlarged test class just proved, including across the internal boundary $\Sigma$. We apply the Green identity there, split the boundary integral along $\Sigma$, and integrate against $\nu$. This gives \[ \int F\dd\Gamma_\nu^{H_2}+\int h\dd\omega_\nu^{H_2} =\int F\dd(\Gamma_\nu^{H_1}+\Gamma_{\tau\nu}^{H_2}) +\int h\dd(\omega_\nu^{H_1}|_{\partial^\dagger\mathscr S_{H_2}} +\omega_{\tau\nu}^{H_2}). \] Taking $h=0$ and $F=0$ separately proves~\ref{i:nested} by Riesz uniqueness. For a finite union of intervals, we add the component identities. Finally, we let $s_+\uparrow\infty$. The Green measures increase, with total mass bounded by $\frac{\mathfrak q(v_e)}{2\lambda}$. Their limit inherits the cone and monotonicity. For tests vanishing after a finite time $T$, the Green integral is already independent of $s_+>T$. The identity then shows, by Riesz uniqueness, that the exit measures restricted to $\{sT$. These consistent restrictions define the infinite-horizon exit measure and give the claimed limiting identity. \end{proof} \section{Construction of the Bellman barrier} \label{s:bellman} We prove Proposition~\ref{p:bellman-barrier}. The first step excludes positive stationary adjoint measures with a certain homogeneity. The second turns that exclusion into a barrier by separation. We work with measures throughout; no smooth adjoint density or extremising coefficient is assumed. \subsection{The homogeneous adjoint exponent} We write $\delta_r(X,v)=(r^3X,rv)$ and $\rho(X,v)=(X^2+v^6)^{\frac16}$. We consider nonzero pairs of positive Radon measures $(\mu,\eta)$ on the punctured plane satisfying \refstepcounter{equation} \[ \begin{gathered} -\partial_X(v\mu)+\partial_{vv}\eta=0, \qquad \lambda\mu\leq\eta\leq\Lambda\mu,\\ \mu(\delta_r E)=r^{4-\beta}\mu(E),\qquad \eta(\delta_r E)=r^{4-\beta}\eta(E). \end{gathered} \tag{\theequation}\label{e:homogeneous-adjoint} \] The exponent $4$ is the homogeneous dimension of position and velocity, without time. We call $\beta$ the density degree, meaning the scaling in this display; a density need not exist. We define $\beta_*$ as the supremum of the degrees for which such a pair exists. The time integral of the constant-coefficient Gaussian kernel gives a positive stationary fundamental measure of degree two. Its source is at the origin, which is excluded here, so $\beta_*\geq2$. \begin{lemma}[A strict gap below degree three] \label{l:adjoint-degree-gap} We have $2\leq\beta_*<3$, and $\beta_*$ depends only on $\frac\Lambda\lambda$. \end{lemma} \begin{proof} We first exclude every $\beta\geq3$. On $X>0$ we write $X=s^3$ and $v=sy$. Homogeneity gives \[ \dd\mu=s^{3-\beta}\dd s\,F(\dd y),\qquad \dd\eta=s^{3-\beta}\dd s\,H(\dd y),\qquad \lambda F\leq H\leq\Lambda F. \] Testing~\eqref{e:homogeneous-adjoint} with separated functions of $s$ and $y$ yields \refstepcounter{equation} \[ \int\varphi''\dd H-\frac13\int y^2\varphi'\dd F +\frac{\beta-2}{3}\int y\varphi\dd F=0. \tag{\theequation}\label{e:angular-adjoint} \] Indeed, the transport derivative of $\zeta(s)\varphi(y)$ is $\frac{y}{3s}\zeta'\varphi-\frac{y^2}{3s^2}\zeta\varphi'$; integration by parts in $s$ gives the last term in~\eqref{e:angular-adjoint}. The distribution $J=H'+\frac13y^2F$ satisfies $J'=-\frac{\beta-2}{3}yF$, so it has a locally bounded-variation representative. Thus $H'$ is a Radon measure. A measure whose derivative is a measure on an interval has a bounded-variation density: subtracting a primitive of the derivative leaves a constant distribution. We obtain $H=h\dd y$, and comparison gives $F=f\dd y$ with $\lambda f\leq h\leq\Lambda f$. The identities then imply $h,J\in W^{1,1}_{\mathrm{loc}}$. We set \[ c_\beta=\frac{\beta-2}{3},\qquad d_\beta=\frac{\beta-3}{3},\qquad I(y)=\int_0^y zf(z)\dd z,\qquad M(y)=\int_0^y z^2f(z)\dd z. \] Two integrations give \refstepcounter{equation} \[ h(y)=h_0+j_0y-c_\beta yI(y)+d_\beta M(y), \qquad h_0=h(0),\quad j_0=J(0). \tag{\theequation}\label{e:angular-integrated} \] Local finiteness on rectangles meeting $X=0$ away from the origin gives the tail condition \refstepcounter{equation} \[ \int_{|y|>2}f(y)|y|^{\beta-4}\dd y<\infty. \tag{\theequation}\label{e:angular-tail} \] For example, on a rectangle with $1\leq v\leq2$, substituting $s=\frac vy$ in $s^{3-\beta}\dd s\,f(y)\dd y$ gives the factor $y^{\beta-4}$ after integration in $v$. Negative velocities give the other tail. For $y>0$, we have $I(y)\geq0$, $0\leq M(y)\leq yI(y)$ and $c_\beta-d_\beta=\frac13$. Non-negativity in~\eqref{e:angular-integrated} implies \[ 0\leq\frac{h(y)}y\leq\frac{h_0}y+j_0-\frac13 I(y). \] Hence $I(\infty)<\infty$. Dominated convergence gives $\frac{M(y)}y\to0$, and therefore \[ \lim_{y\to\infty}\frac{h(y)}y=j_0-c_\beta I(\infty)\geq0. \] A positive limit would force $f(y)\geq cy$ for large $y$, contradicting~\eqref{e:angular-tail}. Thus $j_0=c_\beta I(\infty)$, and \[ h(y)=h_0+c_\beta y\bigl(I(\infty)-I(y)\bigr)+d_\beta M(y)\geq h_0. \] Since $\beta\geq3$, the weight in~\eqref{e:angular-tail} has divergent integral at infinity. Comparison forces $h_0=0$. Near zero,~\eqref{e:angular-integrated} now gives $h(y)=j_0y+o(|y|)$. Non-negativity on both sides of zero forces $j_0=0$. Consequently $I(\infty)=0$, and $f=h=0$ on the positive half-line. For $y=-z<0$, we have $0\leq-M(-z)\leq zI(-z)$. With $h_0=j_0=0$,~\eqref{e:angular-integrated} gives \[ h(-z)\geq(c_\beta-d_\beta)zI(-z)=\frac z3I(-z). \] Any positive mass of $f$ on the negative half-line would make $I(-z)$ bounded below for large $z$, again contradicting~\eqref{e:angular-tail}. Both angular measures vanish. Reflection $(X,v)\mapsto(-X,-v)$ gives the same conclusion on $X<0$. The remaining pair, if any, is supported on $X=0$, $v\ne0$. Writing $\mu=\delta_0(\dd X)m(\dd v)$ and $\eta=\delta_0(\dd X)n(\dd v)$, the equation becomes \[ -\delta_0'(X)\otimes(vm)+\delta_0(X)\otimes n''=0. \] The independent normal distributions force $vm=0$. Since the origin is excluded, $m=0$, and comparison gives $n=0$. This contradicts nontriviality. Excluding degree three alone does not yet give a strict gap. We suppose, for a contradiction, that admissible degrees $\beta_j\uparrow3$ exist. We choose $0\leq\chi=\chi_0(\rho)\leq1$, equal to one on $1\leq\rho\leq2$ and supported in $\frac12<\rho<4$, and normalise by $\int\chi\dd\mu_j=1$. The normalising integral is positive because a dilation of this annulus meets every nonzero homogeneous measure. On $A_0=\{1\leq\rho\leq2\}$ we have $\mu_j(A_0)\leq1$, and homogeneity gives \[ \mu_j(\delta_{2^k}A_0)=2^{k(4-\beta_j)}\mu_j(A_0) \qquad(k\in\mathbb Z). \] Every compact subset of the punctured plane is covered by finitely many of these annuli. The degrees are bounded, so their masses are uniformly bounded; the same holds for $\eta_j\leq\Lambda\mu_j$. A subsequence therefore converges vaguely to a pair $(\mu,\eta)$. The normalisation passes to the limit, giving $\int\chi\dd\mu=1$. Testing with compactly supported smooth functions passes the adjoint equation and ellipticity inequalities to the limit. For a fixed $r>0$ and compactly supported continuous $F$, we pass to the limit in \[ \int F\circ\delta_r\dd\mu_j=r^{\beta_j-4}\int F\dd\mu_j, \] and likewise for $\eta_j$. Thus the limit is a nonzero pair of degree three, which we have just excluded. This proves $\beta_*<3$. Finally, the map $T(X,v)=(\lambda X,v)$ commutes with $\delta_r$. The transformed measures \[ \widetilde\mu=T_\#\mu,\qquad \widetilde\eta=\lambda^{-1}T_\#\eta \] satisfy the same stationary equation, have the same degree, and obey $\widetilde\mu\leq\widetilde\eta\leq\frac\Lambda\lambda\widetilde\mu$. This invertible correspondence proves the claimed dependence on the ratio alone. \end{proof} \subsection{From the degree gap to a universal barrier} We seek one homogeneous function $\Phi$ for which $v\Phi_X+b\Phi_{vv}>0$ for every $b\in[\lambda,\Lambda]$. If no such function existed, separation would supply a positive measure that annihilates all these expressions. Extending that measure by homogeneity produces a stationary adjoint pair of degree $2+\alpha$. Choosing $2+\alpha>\beta_*$ rules out this alternative and yields the barrier. \begin{proof}[Proof of Proposition~\ref{p:bellman-barrier}] We fix $\beta_*-2<\alpha<1$ and put $\Sigma=\{\rho=1\}$. In the space $C(\Sigma\times[\lambda,\Lambda])$, we consider the linear subspace of functions \[ (z,b)\longmapsto v\Phi_X(z)+b\Phi_{vv}(z),\qquad z=(X,v), \] where $\Phi$ ranges over the $C^2$, $\alpha$-homogeneous functions on the punctured plane. If this subspace contained no strictly positive function, neither would its closure: strict positivity on a compact set is open in the uniform norm. Hahn--Banach separation from that open cone then gives a nonzero continuous functional which vanishes on the subspace and is non-negative on the cone. By approximation with strictly positive functions it is positive on all non-negative functions. Riesz representation and normalisation give a probability measure $\Pi$ on $\Sigma\times[\lambda,\Lambda]$. We extend it radially to two positive measures: \[ \begin{aligned} \overline\mu(F)&=\int\int_0^\infty r^{1-\alpha}F(\delta_rz)\dd r\dd\Pi(z,b),\\ \overline\eta(F)&=\int\int_0^\infty br^{1-\alpha}F(\delta_rz)\dd r\dd\Pi(z,b). \end{aligned} \] They are nonzero and locally finite on the punctured plane, obey the ellipticity comparison, and have density degree $2+\alpha$. For $\zeta\in C_c^2(\R^2\setminus\{0\})$, the function \[ \Phi_\zeta(z)=\int_0^\infty r^{-1-\alpha}\zeta(\delta_rz)\dd r \] is $C^2$ and $\alpha$-homogeneous. The identity \[ (v\partial_X+b\partial_{vv})[\zeta(\delta_rz)] =r^2\bigl[(rv)\zeta_X(\delta_rz)+b\zeta_{vv}(\delta_rz)\bigr] \] shows that its annihilation by $\Pi$ is exactly \[ \int v\zeta_X\dd\overline\mu+\int\zeta_{vv}\dd\overline\eta=0. \] Thus $(\overline\mu,\overline\eta)$ is admissible in~\eqref{e:homogeneous-adjoint} with degree $2+\alpha>\beta_*$, a contradiction. The subspace therefore contains a strictly positive function. Its positive minimum on the compact set, followed by homogeneity, gives~\eqref{e:bellman-barrier}. \end{proof} We finally check the origin, as required in Lemma~\ref{l:box-occupation}. The derivatives $\Phi_X$, $\Phi_v$ and $\Phi_{vv}$ have homogeneous degrees $\alpha-3$, $\alpha-1$ and $\alpha-2$. All are greater than $-4$, so the derivatives are locally integrable. On $\partial\{\rho<\delta\}$, the normal surface elements in the $X$ and $v$ directions have sizes $O(\delta)$ and $O(\delta^3)$. Integration by parts gives an inner boundary term $O(\delta^{\alpha+1})$ for one position derivative, and terms $O(\delta^{\alpha+3})$ and $O(\delta^{\alpha+2})$ for two velocity derivatives. They tend to zero. No distribution supported at the origin is omitted, and~\eqref{e:bellman-barrier} holds distributionally on the whole plane. \section{The autonomous counterexample} \label{s:counterexample} A stationary singularity gives a stronger obstruction than the shrinking smooth bump. We construct a positive homogeneous function $H$ which solves the homogeneous equation away from the origin and vanishes there. Flattening $1-H$ near the origin preserves its height and confines its source to a small set. A decaying velocity barrier then removes the initial data without increasing the positive source. \begin{proposition}[No smaller threshold uniform over ellipticity ratios] \label{p:counterexample} For every $d\geq1$ and $1\leq p<4d$, there are fixed $0<\lambda\leq\Lambda$, a fixed kinetic cylinder $Q$, an autonomous measurable matrix $A(x,v)$ with $\lambda I_d\leq A\leq\Lambda I_d$, and smooth non-negative functions $U_j$ such that \[ U_j=0\ \text{on the initial and full lateral boundary of }Q, \qquad \sup_Q U_j\geq\frac14, \qquad \| (\mathcal L_AU_j)_+\|_{L^p(Q)}\longrightarrow0. \] The same conclusion holds with smooth autonomous matrices $A_j$ in place of $A$, with the same ellipticity bounds. Here $\mathcal L_A=\partial_t+v\cdot\nabla_x-A:D_v^2$. The bounds may depend on $d,p$, but do not change along the sequence. \end{proposition} \subsection{A positive homogeneous profile} We fix $0<\alpha<1$ and seek a coefficient $A$ and a function $H$ with \refstepcounter{equation} \[ H(r^3x,rv)=r^\alpha H(x,v),\qquad H\asymp\rho^\alpha,\qquad A:D_v^2H=v\cdot\nabla_xH, \qquad \rho=(|x|^2+|v|^6)^{\frac16}. \tag{\theequation}\label{e:counter-profile} \] The equation holds away from the origin. We construct $H$ separately for $d\geq2$ and $d=1$. \emph{Dimensions at least two.} For $x\ne0$, we put $e=\frac{x}{|x|}$, $y=\frac{v}{|x|^{\frac13}}$ and \[ H(x,v)=|x|^{\frac\alpha3}\phi(y,e). \] Then $D_v^2H=|x|^{\frac{\alpha-2}{3}}D_y^2\phi$, while $v\cdot\nabla_xH=|x|^{\frac{\alpha-2}{3}}T_\alpha\phi$, where \[ T_\alpha\phi=(y\cdot e)\left(\frac\alpha3\phi-\frac13y\cdot D_y\phi\right) +(y-(y\cdot e)e)\cdot D_e\phi. \] Here $D_e$ is the tangential gradient on the unit sphere. We arrange for $D_y^2\phi$ to have a positive eigenvalue everywhere, and for $T_\alpha\phi>0$ wherever that Hessian has no negative eigenvalue. A positive matrix can then weight its eigendirections to give the required trace. We choose $\sigma>0$ with $\frac{\sigma}{\sqrt{1-\alpha}}<\frac14$ and start with \[ \psi(y,e)=C_0+(\sigma^2+|y-e|^2)^{\frac\alpha2}. \] Writing $s=|y-e|$, its tangential and radial Hessian eigenvalues are \[ \alpha(\sigma^2+s^2)^{\frac\alpha2-1},\qquad \alpha(\sigma^2+s^2)^{\frac\alpha2-2}\bigl(\sigma^2-(1-\alpha)s^2\bigr). \] Thus the Hessian has both signs outside $s\leq\frac{\sigma}{\sqrt{1-\alpha}}$. On that small ball, $y\cdot e\geq\frac34$. Increasing $C_0$ adds $\frac\alpha3(y\cdot e)C_0$ to the transport term, so we choose it to make $T_\alpha\psi\geq1$ there. To extend through $x=0$, we remove the dependence on $e$ for large $y$. We choose a smooth radial cutoff $\chi_R$, equal to one on $|y|\leq R$ and zero on $|y|\geq2R$, with $|D^k\chi_R|\leq C_kR^{-k}$, and set \[ \phi(y,e)=|y|^\alpha+\chi_R(y)(\psi(y,e)-|y|^\alpha). \] Near zero this means $\phi=\psi$, so the apparent singular terms cancel. On the transition annulus, \[ \|D_y^2\phi-D_y^2|y|^\alpha\| \leq C\bigl(C_0R^{-2}+R^{\alpha-3}+\sigma^2R^{\alpha-4}\bigr). \] The eigenvalues of $D_y^2|y|^\alpha$ are $\alpha|y|^{\alpha-2}$ in the tangential directions and $-\alpha(1-\alpha)|y|^{\alpha-2}$ radially. Since $\alpha>0$, choosing $R$ large preserves both signs on the annulus. Hence $M=D_y^2\phi$ always has a positive eigenvalue, and $b=T_\alpha\phi>0$ wherever $M\geq0$. Outside $|y|\geq2R$, $b=0$. We write $m_+=\operatorname{tr}M_+$, $m_-=\operatorname{tr}M_-$ and $\Pi_+,\Pi_-,\Pi_0$ for the spectral projections of $M$. Thus $m_+>0$ and $m_-=0$ implies $b>0$. On the compact set $\{|y|\leq2R\}\times\mathbb S^{d-1}$, we can choose $c_-\geq1$ so large that $b+c_-m_->0$: on the closed set $\{b\leq0\}$, $m_-$ has a positive minimum. The matrix \[ A=\frac{b+c_-m_-}{m_+}\Pi_++c_-\Pi_-+\Pi_0 \] is measurable, positive definite and satisfies $A:M=b$. Its eigenvalues have fixed positive lower and finite upper bounds on that compact set. Outside it, $\frac{b+c_-m_-}{m_+}=\frac{c_-(1-\alpha)}{d-1}$, so the bounds are global. We use this matrix at $(y,e)$ when $x\ne0$. Near $(0,v)$ with $v\ne0$, the profile equals $|v|^\alpha$. It extends smoothly there. At $x=0$, $v\ne0$, we set \[ A(0,v)=c_-\widehat v\otimes\widehat v +\frac{c_-(1-\alpha)}{d-1}(I_d-\widehat v\otimes\widehat v), \qquad \widehat v=\frac v{|v|}, \] so the equation continues to hold. We set $A(0,0)=I_d$ and $H(0,0)=0$, including the value $1$ in the ellipticity bounds. Positivity and compactness on $\{\rho=1\}$ prove $H\asymp\rho^\alpha$; homogeneity also gives $|D_vH|\leq C\rho^{\alpha-1}$. \emph{Dimension one.} We put $\gamma=\frac\alpha3$ and choose \[ a_\Lambda(x,v)=\begin{cases}\Lambda,&xv<0,\\1,&xv\geq0.\end{cases} \] Here $\lambda=1$, and the ratio $\Lambda>1$ will be fixed below. The values on the axes do not affect the almost-everywhere equation. For $x>0$, the ansatz $H(x,v)=x^\gamma F(s)$, $s=-\frac{v}{x^{\frac13}}$, reduces the equation to \[ a(s)F''-\frac{s^2}{3}F'+\gamma sF=0, \qquad a(s)=\Lambda\ (s>0),\quad a(s)=1\ (s<0). \] With $z=\frac{s^3}{9a}$ this is Kummer's equation. We use its standard solutions $M(a,b,z)$ and $U(a,b,z)$~\cite[Sections~13.2, 13.4 and 13.7]{DLMF}, and set \[ \begin{aligned} F_+(s)&=U\left(-\gamma,\frac23,\frac{s^3}{9\Lambda}\right),\qquad s>0,\\ F_-(s)&=A_\gamma M\left(-\gamma,\frac23,\frac{s^3}{9}\right) +\frac{B_\gamma s}{(9\Lambda)^{\frac13}} M\left(\frac13-\gamma,\frac43,\frac{s^3}{9}\right),\qquad s<0,\\ A_\gamma&=\frac{\Gamma(\frac13)}{\Gamma(\frac13-\gamma)}, \qquad B_\gamma=\frac{\Gamma(-\frac13)}{\Gamma(-\gamma)}. \end{aligned} \] The connection formula gives the same value and first derivative at $s=0$; the equation gives $F''(0)=0$. Both constants $A_\gamma,B_\gamma$ are positive. We extend by $H(-x,-v)=H(x,v)$. The large-argument expansions give \[ \begin{aligned} F_+(s)&=(9\Lambda)^{-\gamma}s^\alpha(1+O(s^{-3})),\qquad s\to+\infty,\\ F_-(s)&=9^{-\gamma}(P_\gamma-Q_\gamma\Lambda^{-\frac13})|s|^\alpha +O(|s|^{\alpha-3}),\qquad s\to-\infty,\\ P_\gamma&=\cos(\pi\gamma)-\frac{\sin(\pi\gamma)}{\sqrt3}>0, \qquad Q_\gamma=\frac{2\sin(\pi\gamma)}{\sqrt3}>0. \end{aligned} \] Continuity across $x=0$ therefore requires \[ P_\gamma-Q_\gamma\Lambda^{-\frac13}=\Lambda^{-\gamma}. \] The left side minus the right side is negative at $\Lambda=1$ and tends to $P_\gamma>0$ as $\Lambda\to\infty$. We fix a finite root. The common trace is $H(0,v)=(9\Lambda)^{-\gamma}|v|^\alpha$. The transformation $U(-\gamma,\frac23,z)=z^{\frac13}U(\frac13-\gamma,\frac43,z)$ and the positive integral representation of the latter give $F_+>0$. The function $F_-$ is positive at zero and for large negative $s$. A negative interior minimum would satisfy $F'=0$, $F''\geq0$, whereas its equation gives $F''=-\gamma sF<0$. A zero minimum is excluded by uniqueness of the ordinary differential equation. Hence $H>0$ away from the origin. The differentiated expansions follow from $M'=\frac abM(a+1,b+1,\cdot)$ and $U'=-aU(a+1,b+1,\cdot)$~\cite[13.3.15 and 13.3.22]{DLMF}, applying the same expansions with shifted parameters and the chain rule for the arguments $\frac{s^3}{9\Lambda}$ and $\frac{s^3}{9}$. They give, with $C_\alpha=(9\Lambda)^{-\gamma}$, \[ \frac{\dd^j}{\dd s^j}\bigl(F_\pm(s)-C_\alpha|s|^\alpha\bigr) =O(|s|^{\alpha-3-j}),\qquad j=0,1,2, \] as $s\to\pm\infty$. In $H_x=x^{\gamma-1}(\gamma F-\frac s3F')$, the leading terms cancel. The remainder is bounded by $C|v|^{\alpha-3}$, locally bounded away from $v=0$. The same expansions give the common traces of $H_v,H_{vv}$ across $x=0$. Across $v=0$, $F,F',F''$ agree. Thus $H_x,H_v,H_{vv}$ are bounded weak derivatives on compact sets away from the origin, with no singular measures on the axes. The equation holds almost everywhere and the weak chain rule applies there. Positivity and homogeneity again give~\eqref{e:counter-profile} and $|H_v|\leq C\rho^{\alpha-1}$. \subsection{Flattening the singularity and removing the initial data} We use either profile and put $\Omega=\{H<1\}$. We choose a smooth convex $\Psi$ on $[0,\infty)$, zero on $[0,1]$, with $0\leq\Psi'\leq1$ and $\Psi(s)=s-c_\Psi$ for $s\geq2$, where $10$, we set \[ u_r=1-c_\Psi r^\alpha-r^\alpha\Psi\left(\frac{H}{r^\alpha}\right). \] Then $0\leq u_r\leq1$ on $\Omega$, $u_r=0$ on its boundary and $u_r(0,0)\to1$. The singularity has been flattened. The chain rule gives the non-negative stationary source \[ f_r=v\cdot\nabla_xu_r-A:D_v^2u_r =r^{-\alpha}\Psi''\left(\frac{H}{r^\alpha}\right)\langle AD_vH,D_vH\rangle. \] Its support lies in $\{r^\alpha\leq H\leq2r^\alpha\}$, a set of volume at most $Cr^{4d}$. There $|D_vH|\leq Cr^{\alpha-1}$, so \refstepcounter{equation} \[ \|f_r\|_{L^p(\Omega)}\leq Cr^{\alpha-2+\frac{4d}{p}}. \tag{\theequation}\label{e:counter-source} \] For the given $p<4d$, we choose $\max\{0,2-\frac{4d}{p}\}<\alpha<1$. All coefficient and geometric choices are now fixed; only $r$ tends to zero. To obtain zero initial data, we choose $R>\sup_{\overline\Omega}|v|$ and set \[ B(t,v)=e^{-\mu t}\left(2-\frac{|v|^2}{R^2}\right), \qquad \mu=\frac{d\lambda}{R^2},\qquad T=\frac{\log8}{\mu}. \] We have $B(0,\cdot)>1$, $0T$ and $\overline\Omega\subset B_{S^3}\times B_S$, and use the fixed cylinder \[ Q=Q_S^-(T,0,0)=(T-S^2,T)\times B_{S^3}\times B_S. \] We extend $U_r$ by zero for $t\leq0$ and outside $\Omega$. This extension has a fixed zero neighbourhood of the initial and lateral faces: $\inf_{\overline\Omega}B(0,\cdot)>1$, and $B\geq\frac18$ on $[0,T]\times\overline\Omega$, whereas $u_r=1-H$ near $\partial\Omega$. For $d\geq2$ the extension is smooth. For $d=1$, we first fix $r$ and convolve in $(x,v)$ at scale $h$. We also extend $f_r$ by zero outside $(0,T)\times\Omega$. The bounded weak derivatives give \[ \mathcal L_A(U_r*\eta_h)=(\mathcal L_AU_r)*\eta_h+e_{r,h}, \qquad \|e_{r,h}\|_{L^p}\longrightarrow0. \] Indeed, the transport error is bounded by $Ch\|\nabla_xU_r\|_{L^p}$. For diffusion, both $(A D_v^2U_r)*\eta_h$ and $A(D_v^2U_r*\eta_h)$ tend to $A D_v^2U_r$ in $L^p$. Since the mollifier is non-negative, \[ (\mathcal L_A(U_r*\eta_h))_+\leq f_r*\eta_h+|e_{r,h}|. \] We choose $h=h(r)$ small enough to keep the error below $r$, retain height at least $\frac14$ and preserve the zero boundary neighbourhood. This gives smooth functions with vanishing positive source. Finally, mollifying $A$ preserves ellipticity. For each fixed smooth function, the additional error tends to zero in $L^p$ because its velocity Hessian is bounded. Choosing this second mollification scale after $h(r)$ proves Proposition~\ref{p:counterexample} with smooth $A_j$ as well. This construction excludes every $p<4d$ for a suitable fixed ratio. It says nothing at $p=4d$, and does not exclude all $p<4d$ at each fixed ratio. In particular, it is consistent with the improvement below four in Theorem~\ref{t:a}. \begin{thebibliography}{99} \raggedright \bibitem{AbedinTralli} F.~Abedin and G.~Tralli, Harnack inequality for a class of Kolmogorov--Fokker--Planck equations in non-divergence form, \emph{Arch. Ration. Mech. Anal.} \textbf{233} (2019), no.~2, 867--900. \href{https://doi.org/10.1007/s00205-019-01370-z}{doi:10.1007/s00205-019-01370-z} \bibitem{Aleksandrov1960} A.~D. Aleksandrov, Certain estimates for the Dirichlet problem, \emph{Dokl. Akad. Nauk SSSR} \textbf{134} (1960), 1001--1004 (Russian); English transl.\ \emph{Soviet Math. Dokl.} \textbf{1} (1960), 1151--1154. \bibitem{Bakelman1961} I.~Ya. Bakelman, Theory of quasilinear elliptic equations, \emph{Sibirsk. Mat. Zh.} \textbf{2} (1961), 179--186 (Russian). \bibitem{DLMF} NIST Digital Library of Mathematical Functions, Chapter~13: Confluent hypergeometric functions. \href{https://dlmf.nist.gov/13}{https://dlmf.nist.gov/13}. \bibitem{Doeblin1940} W.~Doeblin, \emph{\'El\'ements d'une th\'eorie g\'en\'erale des cha\^ines simples constantes de Markoff}, Ann. Sci. \'Ecole Norm. Sup. (3) \textbf{57} (1940), 61--111. \bibitem{DongYastrzhembskiy} H.~Dong and T.~Yastrzhembskiy, Global $L_p$ estimates for kinetic Kolmogorov--Fokker--Planck equations in nondivergence form, \emph{Arch. Ration. Mech. Anal.} \textbf{245} (2022), no.~1, 501--564. \href{https://doi.org/10.1007/s00205-022-01786-0}{doi:10.1007/s00205-022-01786-0} \bibitem{FerrettiSafonov} E.~Ferretti and M.~V. Safonov, Growth theorems and Harnack inequality for second order parabolic equations, in: \emph{Harmonic Analysis and Boundary Value Problems} (Fayetteville, AR, 2000), Contemp. Math., vol.~277, American Mathematical Society, Providence, RI, 2001, pp.~87--112. \href{https://doi.org/10.1090/conm/277/04540}{doi:10.1090/conm/277/04540} \bibitem{GilbargTrudinger} D.~Gilbarg and N.~S. Trudinger, \emph{Elliptic Partial Differential Equations of Second Order}, 2nd ed., Grundlehren der mathematischen Wissenschaften, vol.~224, Springer, Berlin, 1983; reprinted in Classics in Mathematics, Springer, Berlin, 2001. \href{https://doi.org/10.1007/978-3-642-61798-0}{doi:10.1007/978-3-642-61798-0} \bibitem{GIMV} F.~Golse, C.~Imbert, C.~Mouhot and A.~F. Vasseur, Harnack inequality for kinetic Fokker--Planck equations with rough coefficients and application to the Landau equation, \emph{Ann. Sc. Norm. Super. Pisa Cl. Sci. (5)} \textbf{19} (2019), no.~1, 253--295. \href{https://doi.org/10.2422/2036-2145.201702_001}{doi:10.2422/2036-2145.201702\_001} \bibitem{GuerandImbert} J.~Guerand and C.~Imbert, Log-transform and the weak Harnack inequality for kinetic Fokker--Planck equations, \emph{J. Inst. Math. Jussieu} \textbf{22} (2023), no.~6, 2749--2774. \href{https://doi.org/10.1017/S1474748022000160}{doi:10.1017/S1474748022000160} \bibitem{GuerandMouhot} J.~Guerand and C.~Mouhot, Quantitative De Giorgi methods in kinetic theory, \emph{J. \'Ec. polytech. Math.} \textbf{9} (2022), 1159--1181. \href{https://doi.org/10.5802/jep.203}{doi:10.5802/jep.203} \bibitem{HaoZhang2026} Z.~Hao and X.~Zhang, Harnack estimates for nondivergence kinetic equations, preprint (2026), \href{https://arxiv.org/abs/2609.15534}{arXiv:2609.15534}. \bibitem{Hormander1967} L.~H\"ormander, Hypoelliptic second order differential equations, \emph{Acta Math.} \textbf{119} (1967), 147--171. \href{https://doi.org/10.1007/BF02392081}{doi:10.1007/BF02392081} \bibitem{ImbertSilvestre} C.~Imbert and L.~Silvestre, The weak Harnack inequality for the Boltzmann equation without cut-off, \emph{J. Eur. Math. Soc. (JEMS)} \textbf{22} (2020), no.~2, 507--592. \href{https://doi.org/10.4171/JEMS/928}{doi:10.4171/JEMS/928} \bibitem{Kolmogorov1934} A.~Kolmogoroff, Zuf\"allige Bewegungen (zur Theorie der Brownschen Bewegung), \emph{Ann. of Math. (2)} \textbf{35} (1934), no.~1, 116--117. \href{https://doi.org/10.2307/1968123}{doi:10.2307/1968123} \bibitem{Krylov1976} N.~V. Krylov, Sequences of convex functions and estimates of the maximum of the solution of a parabolic equation, \emph{Sibirsk. Mat. Zh.} \textbf{17} (1976), no.~2, 290--303 (Russian); English transl.\ \emph{Siberian Math. J.} \textbf{17} (1976), no.~2, 226--236. \href{https://doi.org/10.1007/BF00967569}{doi:10.1007/BF00967569} \bibitem{KrylovSafonov} N.~V. Krylov and M.~V. Safonov, A certain property of solutions of parabolic equations with measurable coefficients, \emph{Izv. Akad. Nauk SSSR Ser. Mat.} \textbf{44} (1980), no.~1, 161--175, 239 (Russian); English transl.\ \emph{Math. USSR-Izv.} \textbf{16} (1981), no.~1, 151--164. \href{https://doi.org/10.1070/IM1981v016n01ABEH001283}{doi:10.1070/IM1981v016n01ABEH001283} \bibitem{Landis} E.~M. Landis, \emph{Second order equations of elliptic and parabolic type}, Translations of Mathematical Monographs, vol.~171, American Mathematical Society, Providence, RI, 1998. \bibitem{Lieberman} G.~M. Lieberman, \emph{Second Order Parabolic Differential Equations}, World Scientific, River Edge, NJ, 1996. \href{https://doi.org/10.1142/3302}{doi:10.1142/3302} \bibitem{DoeblinFourier} A.~Loher, C.~Mooney and C.~Mouhot, D\"oblin--Fourier cancellation and kinetic Aleksandrov estimates, arXiv:2610.03083 (2026). \bibitem{NiebelForthcoming} A.~Loher, C.~Mooney, C.~Mouhot, L.~Niebel, L.~Silvestre and Y.~Zhu, Forthcoming work. \bibitem{Pucci1966} C.~Pucci, Limitazioni per soluzioni di equazioni ellittiche, \emph{Ann. Mat. Pura Appl. (4)} \textbf{74} (1966), 15--30. \href{https://doi.org/10.1007/BF02416445}{doi:10.1007/BF02416445} \bibitem{Tso1985} K.~Tso, On an Aleksandrov--Bakel'man type maximum principle for second-order parabolic equations, \emph{Comm. Partial Differential Equations} \textbf{10} (1985), no.~5, 543--553. \href{https://doi.org/10.1080/03605308508820388}{doi:10.1080/03605308508820388} \bibitem{Zhu2026} Y.~Zhu, Stationary kinetic Krylov--Safonov estimates in space dimension one, preprint (2026), \href{https://arxiv.org/abs/2608.31158}{arXiv:2608.31158}. \end{thebibliography} \end{document}