 Solution

Write $K(u,f)=\Vert u{\Vert }_{C^2(B_1)}+\Vert f{\Vert }_{C^{0,\alpha }(B_1)}$; the paper's $|u{|}_{2;B_1}$ denotes a $C^2$ norm, not an $L^2$ norm. We use a blow-up compactness proof of an interior Schauder estimate. If the estimate failed for fixed $\delta >0$, normalization would produceChoose $x_j,y_j\in B_{1/2}$ with $d_j=|y_j-x_j|$ and $|D^2{u}_j(y_j)-D^2{u}_j(x_j)|\ge \frac{1}{2}{d}_j^{\alpha }$. Since $\Vert D^2{u}_j{\Vert }_{\infty }\to 0$, necessarily $d_j\to 0$. Subtract the quadratic Taylor polynomial $P_j$ at $x_j$ and setThe domains expand to ${\mathbb{R}}^n$. At zero, $w_j,D{w}_j,D^2{w}_j$ vanish; the Hölder seminorm of $D^2{w}_j$ is at most ${\delta }^{-1}$. Integration along segments gives uniform $C^{2,\alpha }$ bounds on each fixed ball. Alsolocally uniformly, since $[f_j{]}_{\alpha }\to 0$. By compact embedding of Hölder spaces and a diagonal subsequence, $w_j$ converges locally in $C^{2,\beta }$ to an entire harmonic function $w$, with $[D^{2}w{]}_{\alpha }\le {\delta }^{-1}$ and $D^{2}w(0)=0$. The unit vectors $(y_j-x_j)/d_j$ have a convergent subsequence; its limit $e$ satisfies $|D^{2}w(e)|\ge 1/2$.

Each second derivative of $w$ is an entire harmonic function with finite global Hölder seminorm. The supplied polynomial-growth Liouville theorem for harmonic functions, at growth exponent $\alpha <1$, makes it constant, and its value at zero makes it zero. This contradicts the nonzero Hessian at $e$, proving the estimate with the small $\delta$ term.

To remove that term, rescale the estimate to all contained balls and apply the Simon absorption lemma. Explicitly, with $d(x)=1-|x|$ and $d_{xy}=\min (d(x),d(y))$, defineFor $|x-y|<d_{xy}/8$, apply the rescaled estimate on $B_{d(x)/2}(x)$; its larger-ball seminorm is controlled by $M$ because all points there have boundary distance at least $d(x)/2$. For more separated pairs use $2\Vert D^{2}u{\Vert }_{\infty }$. The result is $M\le C_0\delta M+C_{\delta }K(u,f)$. Choose $C_0\delta <1/2$ and use $d_{xy}\ge 1/2$ on $B_{1/2}$ to obtainThe forcing's Hölder seminorm cannot be discarded. For integers $k\to \infty$, takeTheir $C^2$ and forcing supremum norms stay bounded, but the $11$ entry of the Hessian matrix differs by two at $x_1=\pm \pi /(2k)$, giving $[D^2{u}_k{]}_{\alpha ;B_{1/2}}\ge 2(k/\pi {)}^{\alpha }\to \infty$.