alphaFeedback

WorksMathematics

A Liouville Theorem for the Real k-Hessian Equation with a One-Sided Growth Condition

AI contributions
Not recorded for this paper
View details

Abstract

We establish a scale-invariant interior C2C^2 estimate for smooth solutions of the real kk-Hessian equation. As an application, we prove that every smooth entire solution of σk(D2u)=1σ_k(D^2u)=1 satisfying u(x)≥b∣x∣2−Bu(x)\ge b|x|^2-B for some constants b,Bb,B with b>0b>0 is a quadratic polynomial. For smooth solutions, admissibility is forced by the growth condition. The results also hold for viscosity solutions and for weak solutions in Trudinger's sense. Admissibility is required in the latter case.

Introduction

We study the rigidity of entire solutions to the real kk-Hessian equation

σk(D2u)=1in Rn,2≤k≤n.(1)\sigma_k(D^2u)=1 \qquad\text{in }\mathbb{R}^n,\qquad2\le k\le n. \tag*{(1)}

Here σj\sigma_j denotes the jj-th elementary symmetric polynomial, evaluated at the eigenvalues of the real Hessian D2uD^2u. We use σ0=1\sigma_0=1 and σj=0\sigma_j=0 when jj exceeds the number of variables.

The operators σk\sigma_k connect the Laplacian and the real Monge–Ampère operator. Their elliptic branches are described by the cones

Γj(m)={λ∈Rm:σℓ(λ)>0 for 1≤ℓ≤j},1≤j≤m.(2)\Gamma_j^{(m)}=\{\lambda\in\mathbb{R}^m:\sigma_\ell(\lambda)>0\text{ for }1\le\ell\le j\},\qquad1\le j\le m. \tag*{(2)}

We write H∈ΓjH\in\Gamma_j when the eigenvalues of a symmetric matrix HH belong to this cone. For smooth solutions of (1.1), we use kk-admissible to mean D2u∈ΓkD^2u\in\Gamma_k at every point.

The weak admissibility conventions are specified in Sections 6.1 and 7.1.

A natural hypothesis is positive quadratic lower growth:

u(x)≥b∣x∣2−Bfor some constants b>0 and B≥0.(3)u(x)\ge b|x|^2-B \qquad\text{for some constants } b>0 \text{ and } B\ge0. \tag*{(3)}

This condition makes large sublevel sets bounded. It does not directly bound the oscillation of a rescaled solution on a fixed Euclidean ball, which is needed in the usual interior-estimate argument.

Earlier estimates and rigidity theorems

Liouville theorems for nonlinear partial differential equations have been studied extensively. For Hessian equations, a central question is quadratic rigidity under growth or convexity assumptions. Bao, Chen, Guan, and Ji proved quadratic rigidity for strictly convex entire solutions under (1.3) [1 Theorem 3.2]. For 2≤k≤n−12 \le k \le n - 1, Li, Ren, and Wang extended this result to (k+1)(k + 1)-admissible solutions using a Pogorelov estimate [11].

Chu and Dinew treated a general class of Hessian equations through weighted gradient and Hessian estimates. For kk-admissible solutions of (1.1), their result applies under (1.3) and the additional bound σk+1(D2u)≥−A\sigma_{k+1}(D^2u) \ge-A [6 Corollary 1.3].

Within the natural class of kk-admissible solutions, Wang and Bao proved rigidity under positive two-sided quadratic growth [13 Theorem 1.1]. Their proof combines interior estimates with rescaling and Hölder continuity of the Hessian.

Chen and Xiang proved that every 2-admissible solution u∈C4(R3)u \in C^4(\mathbb{R}^3) of σ2(D2u)=1\sigma_2(D^2u) = 1 satisfying (1.3) is a quadratic polynomial [4 Theorem 1.4].

Bin Wang proved a lower-growth rigidity theorem for kk-admissible solutions when 2k>n2k > n [12 Theorem 1.2]. His argument establishes a Pogorelov estimate independent of the gradient bound.

For the complex Monge–Ampère equation det⁡(vijˉ)=1\det(v_{i\bar{j}}) = 1, Chu and Liu proved quadratic rigidity for smooth entire plurisubharmonic solutions under positive two-sided quadratic growth [7 Theorem 1.2]. At the same time, using a different method, Hu and Sheng established interior C2C^2 estimates and proved quadratic rigidity for locally bounded entire plurisubharmonic Bedford–Taylor solutions under weak quadratic growth conditions [8 Theorems 1.1 and 1.3]. In particular, the Liouville theorem in that paper permits a positive quadratic lower bound together with an upper bound of order (1+∣z∣2)log⁡(e+∣z∣2)(1 + |z|^2)\log(e + |z|^2). We expect that this upper growth bound can be further weakened.

Motivated by these developments, particularly the recent work of Wang [12], we ask whether the upper growth bound can be removed entirely for the real kk-Hessian equation, so that (1.3) alone implies quadratic rigidity.

The main results

We first state the Liouville theorem, followed by the interior estimate used in its proof.

Theorem 1.1 (Liouville theorem). Let 2≤k≤n2 \le k \le n, and suppose u∈C∞(Rn)u \in C^\infty(\mathbb{R}^n) satisfies

σk(D2u)=1in Rn.(4)\sigma_k(D^2u) = 1 \qquad\text{in } \mathbb{R}^n. \tag*{(4)}

Suppose there exist constants b>0b > 0 and B≥0B \ge0 such that

u(x)≥b∣x∣2−B.(5)u(x) \ge b|x|^2 - B. \tag*{(5)}

Then

u(x)=12xTAx+ℓ⋅x+c,σk(A)=1,A≥2bI,(6)u(x) = \tfrac{1}{2}x^T A x + \ell\cdot x + c,\qquad\sigma_k(A) = 1,\qquad A \ge2bI, \tag*{(6)}

where AA is a constant symmetric matrix, ℓ∈Rn\ell\in\mathbb{R}^n, and c∈Rc \in\mathbb{R}.

The following estimate is essential to the proof of the Liouville theorem. Its constants depend only on the dimension, the order of the equation, and the coefficient of the quadratic lower bound.

Theorem 1.2 (Interior C2C^2 estimate). Let 2≤k≤n2 \le k \le n and b>0b > 0. There exist constants ρ=ρ(n,k,b)∈(0,1)\rho= \rho(n,k,b) \in(0,1) and C=C(n,k,b)>0C = C(n,k,b) > 0 with the following property.

For R>0R > 0, suppose u∈C∞(BR)u \in C^\infty(B_R) satisfies

σk(D2u)=1in BR,u(0)=0,(7)\sigma_k(D^2u) = 1 \quad\text{in } B_R,\qquad u(0) = 0, \tag*{(7)}

and, for some constant BB,

u(x)≥b∣x∣2−B(x∈BR),0≤B≤bR24.(8)u(x) \ge b|x|^2-B \quad(x\in B_R), \qquad0\le B\le\frac{bR^2}{4}. \tag*{(8)}

Then

sup⁡BρR(∣u∣R2+∣Du∣R+∣D2u∣)≤C.(9)\sup_{B_{\rho R}}\left(\frac{|u|}{R^2}+\frac{|Du|}{R}+|D^2u|\right)\le C. \tag*{(9)}

The Liouville theorem and the interior estimate extend to viscosity solutions, as stated in the following two theorems.

Theorem 1.3 (Liouville theorem for viscosity solutions). Let 2≤k≤n2\le k\le n, and let u∈C(Rn)u\in C(\mathbb{R}^n) be a viscosity solution, in the sense of § 6.1, of

σk(D2u)=1in Rn.(10)\sigma_k(D^2u)=1 \quad\text{in }\mathbb{R}^n. \tag*{(10)}

Suppose that, for some constants b>0b>0 and B≥0B\ge0,

u(x)≥b∣x∣2−B.(11)u(x)\ge b|x|^2-B. \tag*{(11)}

Then

u(x)=12xTAx+ℓ⋅x+c,σk(A)=1,A≥2bI,(12)u(x)=\frac{1}{2}x^T A x+\ell\cdot x+c,\qquad\sigma_k(A)=1,\qquad A\ge2bI, \tag*{(12)}

where AA is a constant symmetric matrix, ℓ∈Rn\ell\in\mathbb{R}^n, and c∈Rc\in\mathbb{R}.

Theorem 1.4 (Interior estimate for viscosity solutions). Let 2≤k≤n2\le k\le n and b>0b>0. There exist constants ρ=ρ(n,k,b)∈(0,1)\rho=\rho(n,k,b)\in(0,1) and C=C(n,k,b)>0C=C(n,k,b)>0 with the following property.

For R>0R>0, suppose u∈C(BR)u\in C(B_R) is a viscosity solution, in the sense of § 6.1, of

σk(D2u)=1in BR,u(0)=0.(13)\sigma_k(D^2u)=1 \quad\text{in }B_R,\qquad u(0)=0. \tag*{(13)}

Assume that, for some constant BB,

u(x)≥b∣x∣2−B(x∈BR),0≤B≤bR24.(14)u(x)\ge b|x|^2-B \quad(x\in B_R),\qquad0\le B\le\frac{bR^2}{4}. \tag*{(14)}

Then u∈C∞(BρR)u\in C^\infty(B_{\rho R}) and

sup⁡BρR(∣u∣R2+∣Du∣R+∣D2u∣)≤C.(15)\sup_{B_{\rho R}}\left(\frac{|u|}{R^2}+\frac{|Du|}{R}+|D^2u|\right)\le C. \tag*{(15)}

Remark 1.5 (Pogorelov’s example and the growth condition). Pogorelov’s example shows that viscosity solutions can have an unbounded Hessian when the growth hypothesis is absent; see [10 Section 4.4]. We use the construction in kk variables and extend it independently of the remaining n−kn-k variables. For 3≤k≤n3\le k\le n, write

x=(y,t,z)∈Rk−1×R×Rn−k,αk=2−2k.(16)x=(y,t,z)\in\mathbb{R}^{k-1}\times\mathbb{R}\times\mathbb{R}^{n-k},\qquad\alpha_k=2-\frac{2}{k}. \tag*{(16)}

Choose a positive smooth solution hh on (−T,T)(-T,T), for some T>0T>0, of

αkk−1hk−2((αk−1)hh′′−αk(h′)2)=1,h(0)=1,h′(0)=0.(17)\alpha_k^{k-1}h^{k-2}\left((\alpha_k-1)hh''-\alpha_k(h')^2\right)=1,\qquad h(0)=1,\qquad h'(0)=0. \tag*{(17)}

The function

P(y,t,z)=∣y∣αkh(t)(18)P(y,t,z)=|y|^{\alpha_k}h(t) \tag*{(18)}

is defined on the slab Rk−1×(−T,T)×Rn−k\mathbb{R}^{k-1}\times(-T,T)\times\mathbb{R}^{n-k}. For r=∣y∣>0r=|y|>0, differentiation gives

det⁡Dy,t2P=αkk−1hk−2((αk−1)hh′′−αk(h′)2)rkαk−2k+2=1.\det D^2_{y,t}P=\alpha_k^{k-1}h^{k-2}\left((\alpha_k-1)hh''-\alpha_k(h')^2\right)r^{k\alpha_k-2k+2}=1.

Here kαk−2k+2=0k\alpha_k-2k+2=0. The k−2k-2 angular eigenvalues are positive, and the radial–tt block has positive radial entry and determinant. Thus the active k×kk\times k Hessian block is positive definite.

The remaining n−kn-k eigenvalues vanish, so σk(D2P)=1\sigma_k(D^2P)=1 off {y=0}\{y=0\}. Since αk>1\alpha_k>1, the gradient extends continuously across this set; restriction to lines therefore shows that PP is convex.

We verify the viscosity equation on {y=0}\{y=0\}, where PP is not C2C^2.

At y=0y=0, no C2C^2 upper touching test exists because 1<αk<21<\alpha_k<2. A lower touching test has Hessian AA with A∣S≤0A|_S\leq0 on the flat subspace S={y=0}S=\{y=0\}. If A∈ΓkA\in\Gamma_k and σk(A)>0\sigma_k(A)>0, Maclaurin’s inequalities give A∈ΓkA\in\Gamma_k.

Restrict successively to k−1k-1 hyperplanes within the preceding subspace, each containing SS. Equation (62) gives A∣S∈Γ1(n−k+1)A|_S\in\Gamma_1^{(n-k+1)}, so tr⁡(A∣S)>0\operatorname{tr}(A|_S)>0, contradicting A∣S≤0A|_S\leq0.

Thus every admissible lower test has σk(A)=0\sigma_k(A)=0. This verifies the supersolution inequality on the singular set.

The growth condition in Theorem 1.4 excludes this example. Indeed, take 0<R<T0<R<T and x∗=(0,3R/4,0)∈BRx_*=(0,3R/4,0)\in B_R. Since P(0)=P(x∗)=0P(0)=P(x_*)=0, the assumed lower bound would imply

0=P(x∗)≥9bR216−B≥5bR216>0.(19)0=P(x_*)\geq\frac{9bR^2}{16}-B\geq\frac{5bR^2}{16}>0. \tag*{(19)}

Thus the quantitative lower bound prevents the flat singular axis of Pogorelov’s example. This example is defined on a slab and concerns local regularity.

The results also extend to admissible weak solutions in the sense of Trudinger [15]. We recall his definition in § 7.1.

Theorem 1.6 (Liouville theorem for admissible weak solutions). Let 2≤k≤n2\leq k\leq n, and let u∈C(Rn)u\in C(\mathbb{R}^n) be an admissible weak solution, in Trudinger’s sense stated in § 7.1, of

σk(D2u)=1in Rn.(20)\sigma_k(D^2u)=1 \qquad\text{in }\mathbb{R}^n. \tag*{(20)}

Suppose that, for some constants b>0b>0 and B≥0B\geq0,

u(x)≥b∣x∣2−B.(21)u(x)\geq b|x|^2-B. \tag*{(21)}

Then

u(x)=12xTAx+ℓ⋅x+c,σk(A)=1,A≥2bI,(22)u(x)=\frac{1}{2}x^TAx+\ell\cdot x+c,\qquad\sigma_k(A)=1,\qquad A\geq2bI, \tag*{(22)}

where AA is a constant symmetric matrix, ℓ∈Rn\ell\in\mathbb{R}^n, and c∈Rc\in\mathbb{R}.

Theorem 1.7 (Interior estimate for admissible weak solutions). Let 2≤k≤n2\leq k\leq n and b>0b>0. There exist constants ρ=ρ(n,k,b)∈(0,1)\rho=\rho(n,k,b)\in(0,1) and C=C(n,k,b)>0C=C(n,k,b)>0 with the following property. For R>0R>0, suppose u∈C(BR)u\in C(B_R) is an admissible weak solution, in Trudinger’s sense stated in § 7.1, of

σk(D2u)=1in BR,u(0)=0.(23)\sigma_k(D^2u)=1 \qquad\text{in }B_R,\qquad u(0)=0. \tag*{(23)}

Assume that, for some constant BB,

u(x)≥b∣x∣2−B(x∈BR),0≤B≤bR24.(24)u(x)\geq b|x|^2-B\quad(x\in B_R),\qquad0\leq B\leq\frac{bR^2}{4}. \tag*{(24)}

Then u∈C∞(BρR)u\in C^\infty(B_{\rho R}) and

sup⁡BρR(∣u∣R2+∣Du∣R+∣D2u∣)≤C.(25)\sup_{B_{\rho R}}\left(\frac{|u|}{R^2}+\frac{|Du|}{R}+|D^2u|\right)\leq C. \tag*{(25)}

Remark 1.8 (Distributional solutions without admissibility). The conclusions of Theorems 1.6 and 1.7 fail for unrestricted distributional solutions of the 2-Hessian equation. We give an example in every dimension n≥2n \ge2.

For w∈Wloc1,2w \in W_{\mathrm{loc}}^{1,2}, the distribution

S2[w]=12∑i,j∂ij(wiwj)−12Δ∣Dw∣2(26)S_2[w] = \frac{1}{2}\sum_{i,j}\partial_{ij}(w_iw_j) - \frac{1}{2}\Delta|Dw|^2 \tag*{(26)}

agrees with σ2(D2w)\sigma_2(D^2w) for smooth ww; see [9 Definition 1.1 and Section 2.1], with our normalization. Its pairing with a test function ϕ\phi is given by the gradient integrals on the right-hand side of (95).

Set

c=(n2)−1/2,q(x)=c2∣x∣2+cx1,w(x)=∣q(x)∣.(27)c = \binom{n}{2}^{-1/2}, \qquad q(x) = \frac{c}{2}|x|^2 + cx_1, \qquad w(x) = |q(x)|. \tag*{(27)}

Since Dw=sgn⁡(q)DqDw = \operatorname{sgn}(q)Dq almost everywhere, we obtain

Dw⊗Dw=Dq⊗Dq,S2[w]=S2[q]=(n2)c2=1.(28)Dw \otimes Dw = Dq \otimes Dq, \qquad S_2[w] = S_2[q] = \binom{n}{2}c^2 = 1. \tag*{(28)}

Moreover, completing the square gives

w(x)≥q(x)≥c4∣x∣2−c.(29)w(x) \ge q(x) \ge\frac{c}{4}|x|^2-c. \tag*{(29)}

Thus ww has positive quadratic lower growth but is not quadratic. On {q<0}\{q < 0\}, its Hessian is −cI-cI, so ww is not 2-admissible.

The conclusion of Theorem 1.7 also fails for this distributional class: w(0)=0w(0) = 0, and the choices b=c/4b = c/4, B=cB = c, and R=4R = 4 satisfy B=bR2/4B = bR^2/4. Nevertheless, w(te1)=∣ct+ct2/2∣w(te_1) = |ct+ct^2/2| is not differentiable at the origin.

The upper-test inequality itself implies admissibility in Theorem 1.3 and Theorem 1.4; see Lemma 6.1. In Theorem 1.6 and Theorem 1.7, admissibility is part of Trudinger’s weak-solution definition.

For the smooth solution in Theorem 1.2, the growth bound ensures that uu attains a minimum inside the ball. At this point its Hessian is positive semidefinite, and the equation gives admissibility. Maclaurin’s inequalities then give positive lower bounds for the lower symmetric polynomials. These bounds prevent the Hessian from reaching the boundary of Γk\Gamma_k; see § 4.2.

To prove Theorem 1.2, we choose a regular sublevel component whose diameter is bounded by a multiple of RR. A distance barrier shows that this component contains a ball of radius comparable to RR.

The same barrier separates every fixed negative sublevel from the boundary. The interior gradient estimate of Chou and Wang then provides the gradient bound needed in their weighted Hessian estimate [5 Theorems 3.2 and 1.5].

The construction has two geometric features. When k<nk < n, the principal curvatures of a (k−1)(k-1)-convex boundary may be negative. The algebraic inequality in Lemma A.2 controls their fractional transformation.

The distance function is also nonsmooth at the cut locus. Local semiconvexity of the distance power passes its almost-everywhere Hessian inequalities to smooth convolutions, which can be used in the comparison argument.

Organization of the paper

§ 2 constructs the admissible distance power, and § 3 proves the uniform boundary estimate. § 4 proves Theorem 1.2.

§ 5 applies this estimate and the Evans–Krylov theorem to prove Theorem 1.1. § 6 proves the viscosity results, Theorems 1.3 and 1.4. § 7 treats admissible weak solutions and proves Theorems 1.6 and 1.7. Appendix § 8 contains the two algebraic lemmas.

Admissible powers of the distance function

We construct a function that vanishes on the boundary and whose Hessian belongs to the closed admissible cone. At smooth points of the distance function, this follows from Lemma A.2.

For 2≤k≤n2 \le k \le n, fix

N=n−k+1,θ=2k(1−12N),α=1−θ,β=αk.(30)\begin{aligned} N &= n-k+1, \qquad\theta= \frac{2}{k}\left(1-\frac{1}{2N}\right), \qquad\alpha= 1-\theta, \qquad\beta= \frac{\alpha}{k}. \tag*{(30)} \end{aligned}

These constants satisfy 0<θ<10 < \theta< 1 and 0<β<α<10 < \beta< \alpha< 1.

Let Ω⊂Rn\Omega\subset\mathbb{R}^{n} be a nonempty bounded smooth domain whose outward principal curvatures belong to Γk−1(n−1)\Gamma_{k-1}^{(n-1)}. We take the outward principal curvatures to be positive on spheres.

Set

D=diam⁡Ω,d(x)=dist⁡(x,∂Ω),v(x)=−d(x)α.(31)D = \operatorname{diam}\Omega, \qquad d(x) = \operatorname{dist}(x,\partial\Omega), \qquad v(x) = -d(x)^{\alpha}. \tag*{(31)}

A function ff is locally semiconvex if, on a convex neighborhood of each point, f(x)+C2∣x∣2f(x) + \frac{C}{2}|x|^{2} is convex for some C≥0C \ge0. Local semiconcavity means that −f-f is locally semiconvex.

The spectral matrix cone Γj\Gamma_j is convex and is preserved under addition of positive-semidefinite matrices. The same statements hold for Γ‾j\overline{\Gamma}_j. These standard properties are summarized in [12 Section 2].

At smooth points of dd, the normal eigenvalue of D2vD^{2}v is positive. Lemma A.2 shows that it compensates for the possible negative tangential eigenvalues, with a choice of α\alpha depending only on nn and kk.

Lemma 2.1 (The distance power). The function v=−dαv=-d^{\alpha} is locally semiconvex. At almost every point,

D2v∈Γ‾k,Δv≥c0:=α(1−α)Dα−2.(32)D^{2}v \in\overline{\Gamma}_{k}, \qquad\Delta v \ge c_{0} := \alpha(1-\alpha)D^{\alpha-2}. \tag*{(32)}

Its convolutions with nonnegative mollifiers satisfy the same inequalities wherever the convolution is defined inside Ω\Omega.

Proof. The distance dd is Lipschitz and hence differentiable almost everywhere. If yy is a nearest boundary point to such an xx, then d(z)≤∣z−y∣d(z) \le|z-y| with equality at xx. Differentiation gives

Dd(x)=x−yd(x),y=x−d(x)Dd(x).Dd(x) = \frac{x-y}{d(x)}, \qquad y = x-d(x)Dd(x).

Thus the nearest boundary point is unique. The focal values of the smooth normal map have measure zero by Sard’s theorem, so it remains to compute at points with a unique, nonfocal nearest boundary point.

Let ν\nu be the outward unit normal. The minimizing segment is normal and points inward, so

x=y−d(x)ν(y).x = y-d(x)\nu(y).

Write Ψ(y,s)=y−sν(y)\Psi(y,s)=y-s\nu(y), and let eie_i be a principal direction at yy. With our sign convention for the curvatures,

DyΨ(y,s)ei=(1−sκi(y))ei.D_y\Psi(y,s)e_i=(1-s\kappa_i(y))e_i.

Choose a unit-speed boundary curve γi\gamma_i through yy with γi′(0)=ei\gamma_i'(0)=e_i. Since yy minimizes the distance to xx, our curvature convention gives

0≤d2dt2∣x−γi(t)∣22∣t=0=1−dκi(y).0\leq\left.\frac{d^2}{dt^2}\frac{|x-\gamma_i(t)|^2}{2}\right|_{t=0}=1-d\kappa_i(y).

Nonfocality makes these inequalities strict, so the normal map has a smooth local inverse.

Let xm→xx_m\to x and choose any nearest boundary points ymy_m. Every convergent subsequence of ymy_m has a nearest point to xx as its limit, because distance is continuous.

Uniqueness of the nearest point and compactness give ym→yy_m\to y. Also d(xm)→d(x)d(x_m)\to d(x), so the pairs (ym,d(xm))(y_m,d(x_m)) remain in the local inverse neighborhood for all sufficiently large mm.

In this neighborhood, the inverse normal coordinate ss equals dd. Thus dd is smooth there, with

Dd(Ψ(y,s))=−ν(y).Dd(\Psi(y,s))=-\nu(y).

Differentiating this local gradient identity in each principal direction gives the eigenvalues of D2dD^2d:

−κ1(y)1−dκ1(y),…,−κn−1(y)1−dκn−1(y), 0,1−dκi(y)>0.(33)-\frac{\kappa_1(y)}{1-d\kappa_1(y)},\ldots,-\frac{\kappa_{n-1}(y)}{1-d\kappa_{n-1}(y)},\ 0,\qquad1-d\kappa_i(y)>0. \tag*{(33)}

To compute the Hessian of vv, apply the chain rule:

D2v=−αdα−1D2d+α(1−α)dα−2Dd⊗Dd.D^2v=-\alpha d^{\alpha-1}D^2d+\alpha(1-\alpha)d^{\alpha-2}Dd\otimes Dd.

The second term acts only in the normal direction, since DdDd is a unit normal. Introduce ai=dκi(y)a_i=d\kappa_i(y) and μi=ai/(1−ai)\mu_i=a_i/(1-a_i). The eigenvalues of D2vD^2v are therefore

αdα−2(μ1,…,μn−1,1−α).(34)\alpha d^{\alpha-2}(\mu_1,\ldots,\mu_{n-1},1-\alpha). \tag*{(34)}

Here a∈Γk−1a\in\Gamma_{k-1} and ai<1a_i<1. Since 1−α=θ1-\alpha=\theta, Lemma A.2 gives D2v∈ΓkD^2v\in\Gamma_k.

The same lemma gives σ1(μ)>0\sigma_1(\mu)>0. Using d≤Dd\leq D and α−2<0\alpha-2<0, we obtain

Δv=αdα−2(σ1(μ)+1−α)≥α(1−α)Dα−2=c0.(35)\Delta v=\alpha d^{\alpha-2}\bigl(\sigma_1(\mu)+1-\alpha\bigr)\geq\alpha(1-\alpha)D^{\alpha-2}=c_0. \tag*{(35)}

The distance dd is locally semiconcave in Ω\Omega; see [3 Proposition 2.2.2]. We can check the needed semiconvexity of vv directly on any convex neighborhood V⋐ΩV\Subset\Omega.

Choose a>0a>0 such that dist⁡(V,∂Ω)≥a\operatorname{dist}(V,\partial\Omega)\geq a. For y∈∂Ωy\in\partial\Omega, write r=∣x−y∣r=|x-y| and e=(x−y)/re=(x-y)/r. Since r≥ar\geq a, differentiation gives

Dx2(−∣x−y∣α)=αrα−2((2−α)e⊗e−I)≥−αaα−2I.D_x^2(-|x-y|^\alpha)=\alpha r^{\alpha-2}\bigl((2-\alpha)e\otimes e-I\bigr)\geq-\alpha a^{\alpha-2}I.

Hence each function −∣x−y∣α+12αaα−2∣x∣2-|x-y|^\alpha+\frac{1}{2}\alpha a^{\alpha-2}|x|^2 is convex on VV. Their supremum is convex as well, and

v(x)=sup⁡y∈∂Ω(−∣x−y∣α).v(x)=\sup_{y\in\partial\Omega}(-|x-y|^\alpha).

This proves local semiconvexity, including at points where the nearest boundary point is not unique.

The distributional Hessian of the locally semiconvex function vv is a matrix-valued Radon measure; see [3 Theorem 2.3.1]. Its Lebesgue decomposition is

D2v=H dx+Ms,H∈Γ‾k a.e.,tr⁡H≥c0.(36)D^2v=H\,dx+M_s,\qquad H\in\overline{\Gamma}_k\ \text{a.e.},\qquad\operatorname{tr}H\geq c_0. \tag*{(36)}

The density HH agrees almost everywhere with the classical Hessian computed above. Locally, v+C2∣x∣2v+\frac{C}{2}|x|^2 is convex for some C≥0C\geq0, so D2v+CI dxD^2v+C I\,dx is a positive matrix-valued measure.

The added quadratic changes only the absolutely continuous part. Hence the singular part MsM_s is positive semidefinite.

Let φρ\varphi_\rho be a nonnegative mollifier of integral one. For vρ=v∗φρv_\rho= v * \varphi_\rho, distributional differentiation gives

D2vρ=H∗φρ+Ms∗φρ.(37)D^2v_\rho= H * \varphi_\rho+ M_s * \varphi_\rho. \tag*{(37)}

The first term is a nonnegative average of matrices in Γ‾k\overline{\Gamma}_k, so it belongs to Γ‾k\overline{\Gamma}_k by convexity. The second term is positive semidefinite and therefore preserves this cone under addition.

Taking the trace in (37) also gives

Δvρ=(tr⁡H)∗φρ+(tr⁡Ms)∗φρ≥c0∫φρ=c0.\Delta v_\rho= (\operatorname{tr} H) * \varphi_\rho+ (\operatorname{tr} M_s) * \varphi_\rho\ge c_0 \int\varphi_\rho= c_0.

These arguments apply whenever the support of the translated mollifier is contained in Ω\Omega.

A uniform boundary estimate

We perturb the distance power by a quadratic function. The perturbation makes the barrier strictly admissible, while its trace lower bound supplies the strict differential inequality used in comparison.

Proposition 3.1 (Uniform boundary estimate). Let Ω⊂Rn\Omega\subset\mathbb{R}^n be a nonempty bounded smooth domain whose outward principal curvatures belong to Γk−1(n−1)\Gamma_{k-1}^{(n-1)}. Suppose U∈C2(Ω)∩C(Ω‾)U \in C^2(\Omega) \cap C(\overline{\Omega}) is kk-admissible and

σk(D2U)=1in Ω,U=0on ∂Ω.(38)\sigma_k(D^2U) = 1 \quad\text{in } \Omega, \qquad U = 0 \quad\text{on } \partial\Omega. \tag*{(38)}

Set D=diam⁡ΩD = \operatorname{diam}\Omega and d(x)=dist⁡(x,∂Ω)d(x) = \operatorname{dist}(x,\partial\Omega). There is a constant C0=C0(n,k)>0C_0 = C_0(n,k) > 0 such that

0≤−U(x)≤C0D2(d(x)D)β(x∈Ω),(39)0 \le-U(x) \le C_0D^2\left(\frac{d(x)}{D}\right)^\beta\qquad(x \in\Omega), \tag*{(39)}

where β\beta is defined in (2.1).

The constant is independent of the boundary curvature and its smoothness norms. We prove the proposition in the following two steps.

Since D2U∈ΓkD^2U \in\Gamma_k, we have ΔU>0\Delta U > 0. The maximum principle and the boundary condition give U≤0U \le0.

Prepare the comparison function

Fix x0∈Ωx_0 \in\Omega and ε>0\varepsilon> 0. Introduce the barrier

w=Kv+ε(∣x−x0∣2−D2),K=Ln,kD2−αε1−k,(40)w = Kv + \varepsilon\left(|x-x_0|^2-D^2\right), \qquad K = L_{n,k}D^{2-\alpha}\varepsilon^{1-k}, \tag*{(40)}

where

Ln,k=22k−1(n−1k−1)α(1−α).(41)L_{n,k} = \frac{2}{2^{k-1}\binom{n-1}{k-1}\alpha(1-\alpha)}. \tag*{(41)}

On ∂Ω\partial\Omega, we have v=0v = 0 and ∣x−x0∣≤D|x-x_0| \le D. Consequently, w≤0=Uw \le0 = U there.

For a symmetric matrix A∈Γ‾kA \in\overline{\Gamma}_k, expansion of the symmetric polynomial gives

σk(KA+2εI)=∑j=0k(n−jk−j)(2ε)k−jKjσj(A).(42)\sigma_k(KA+2\varepsilon I) = \sum_{j=0}^{k} \binom{n-j}{k-j}(2\varepsilon)^{k-j}K^j\sigma_j(A). \tag*{(42)}

To see the coefficient, choose jj eigenvalue factors from KAKA in each product defining σk\sigma_k; the remaining k−jk-j factors contribute 2ε2\varepsilon and can be chosen in (n−jk−j)\binom{n-j}{k-j} ways.

Since A∈Γ‾kA \in\overline{\Gamma}_k, every σj(A)\sigma_j(A) in this expansion is nonnegative. Keeping the j=1j=1 term gives

σk(KA+2εI)≥(n−1k−1)(2ε)k−1KtrA.(43)\sigma_k(KA+2\varepsilon I)\ge\binom{n-1}{k-1}(2\varepsilon)^{k-1}K\mathop{\rm tr} A. \tag*{(43)}

Apply Lemma 2.1 to vρv_\rho, and set wρ=Kvρ+ε(∣x−x0∣2−D2)w_\rho=Kv_\rho+\varepsilon(|x-x_0|^2-D^2). Adding 2εI2\varepsilon I makes all the symmetric polynomials of orders 1,…,k1,\ldots,k strictly positive, by the same expansion.

The powers of DD and ε\varepsilon in KK cancel those in (2ε)k−1c0(2\varepsilon)^{k-1}c_0. Thus (2.3) and (3.4) give the fixed strict inequality

D2wρ∈Γk,σk(D2wρ)≥(n−1k−1)(2ε)k−1Kc0=2.(44)D^2w_\rho\in\Gamma_k,\qquad\sigma_k(D^2w_\rho)\ge\binom{n-1}{k-1}(2\varepsilon)^{k-1}Kc_0=2. \tag*{(44)}

Exclude the interior maximum

We prove w≤Uw\le U. Suppose M:=max⁡Ω‾(w−U)>0M:=\max_{\overline{\Omega}}(w-U)>0. Continuity gives δ>0\delta>0 such that w−U<M/4w-U<M/4 whenever d≤2δd\le2\delta.

The compact set Kδ={x∈Ω:d(x)≥δ}K_\delta=\{x\in\Omega:d(x)\ge\delta\} contains a point xMx_M where w−U=Mw-U=M. Every point of ∂Kδ\partial K_\delta satisfies d=δd=\delta, so the collar estimate applies there.

Choose ρ<δ/2\rho<\delta/2 so small that ∥wρ−w∥L∞(Kδ)<M/4\lVert w_\rho-w\rVert_{L^\infty(K_\delta)}<M/4.

At xMx_M we have wρ−U>3M/4w_\rho-U>3M/4, whereas on ∂Kδ\partial K_\delta we have wρ−U<M/2w_\rho-U<M/2. Therefore its maximum on KδK_\delta occurs at an interior point, where the mollified barrier is smooth.

At that point, D2wρ≤D2UD^2w_\rho\le D^2U. Monotonicity of σk\sigma_k on Γk\Gamma_k and (3.7) give

2≤σk(D2wρ)≤σk(D2U)=1.(45)2\le\sigma_k(D^2w_\rho)\le\sigma_k(D^2U)=1. \tag*{(45)}

This contradiction excludes a positive maximum and proves w≤Uw\le U.

Evaluating at x0x_0 yields

−U(x0)≤Ln,kD2−αε1−kd(x0)α+εD2.(46)-U(x_0)\le L_{n,k}D^{2-\alpha}\varepsilon^{1-k}d(x_0)^\alpha+\varepsilon D^2. \tag*{(46)}

We balance the two terms by choosing ε=Ln,k1/k(d(x0)/D)α/k\varepsilon=L_{n,k}^{1/k}(d(x_0)/D)^{\alpha/k}. Since β=α/k\beta=\alpha/k, this choice satisfies

Ln,kD2−αε1−kd(x0)α=εD2=Ln,k1/kD2(d(x0)D)β.L_{n,k}D^{2-\alpha}\varepsilon^{1-k}d(x_0)^\alpha=\varepsilon D^2=L_{n,k}^{1/k}D^2\left(\frac{d(x_0)}{D}\right)^\beta.

Substitution into (3.9) proves (3.2) with C0=2Ln,k1/kC_0=2L_{n,k}^{1/k}. The point x0x_0 was arbitrary.

Proof of the interior estimate

We first convert Proposition 3.1 into an interior estimate for the zero-boundary problem. We then apply that estimate to a sublevel component of the solution in Theorem 1.2.

Estimate the solution away from the boundary

Lemma 4.1 (Interior estimate on a zero-boundary domain). Let Ω\Omega satisfy the assumptions of Proposition 3.1, and let U∈C∞(Ω)∩C(Ω‾)U\in C^\infty(\Omega)\cap C(\overline{\Omega}) be an admissible solution of (3.1). For $D=\operatorname{diam}\Omega

and 0<δ<10 < \delta< 1,

sup⁡{d≥δD}(∣U∣D2+∣DU∣D+∣D2U∣)≤C(n,k,δ),(47)\sup_{\{d \ge\delta D\}}\left(\frac{|U|}{D^2}+\frac{|DU|}{D}+|D^2U|\right)\le C(n,k,\delta), \tag*{(47)}

whenever the set in the supremum is nonempty.

Proof. Translation and the rescaling U(Dy)/D2U(Dy)/D^2 reduce the proof to D=1D=1. The Hessian is unchanged by this scaling, while the gradient is divided by DD and the function by D2D^2.

Proposition 3.1 then gives

0≤−U≤C0dβ,∥U∥L∞(Ω)≤C0.(48)0\le-U\le C_0d^\beta,\qquad\|U\|_{L^\infty(\Omega)}\le C_0. \tag*{(48)}

Fix η>0\eta>0. If −U(x)≥η/4-U(x)\ge\eta/4, the boundary estimate forces C0d(x)β≥η/4C_0d(x)^\beta\ge\eta/4. Solving for d(x)d(x) gives

−U(x)≥η4⟹d(x)≥rη:=(η4C0)1/β.(49)-U(x)\ge\frac{\eta}{4}\quad\Longrightarrow\quad d(x)\ge r_\eta:=\left(\frac{\eta}{4C_0}\right)^{1/\beta}. \tag*{(49)}

For the constant right-hand side, all derivatives vanish, so condition (3.13) of [5] holds with h=0h=0. Apply their interior gradient estimate [CW01, Theorem 3.2]} on balls of radius rη/2r_\eta/2 centered in this sublevel.

Together with (4.2), it gives

sup⁡{−U≥η/4}∣DU∣≤C(n,k,η).(50)\sup_{\{-U\ge\eta/4\}}|DU|\le C(n,k,\eta). \tag*{(50)}

The balls lie inside Ω\Omega by (4.3). Their radius and the oscillation of UU depend only on n,k,ηn,k,\eta, which explains the uniform constant in (4.4).

If the sublevel is empty, the corresponding estimate is vacuous.

Choose s∈(η/4,η/2)s\in(\eta/4,\eta/2) such that −s-s is a regular value of UU, and put

Ωs={U<−s},V=U+s.(51)\Omega_s=\{U<-s\},\qquad V=U+s. \tag*{(51)}

The set Ωs\Omega_s lies in {−U≥η/4}\{-U\ge\eta/4\}, so (4.3) places it a positive distance from ∂Ω\partial\Omega. Regularity of the chosen value makes its boundary smooth.

Adding ss leaves the Hessian unchanged. Thus VV is admissible, solves the same equation, and vanishes on ∂Ωs\partial\Omega_s.

The regular level is (k−1)(k-1)-convex, by the tangential Hessian identity (62) below. Equations (4.2) and (4.4) bound the C1C^1 norm of VV.

Since Ωs⋐Ω\Omega_s\Subset\Omega, the function VV is smooth on a neighborhood of Ωs‾\overline{\Omega_s}. The right-hand side is identically one, so the remaining hypotheses of the estimate are satisfied.

Apply the zero-boundary estimate of Chou and Wang [5 Theorem 1.5] on each component of Ωs\Omega_s:

(−V)4∣D2V∣≤C(n,k,η).(52)(-V)^4|D^2V|\le C(n,k,\eta). \tag*{(52)}

The constant in their estimate is independent of the domain. For a point with −U≥η-U\ge\eta, the choice s<η/2s<\eta/2 gives

−V=−U−s≥η−s>η2.-V=-U-s\ge\eta-s>\frac{\eta}{2}.

Dividing (4.6) by this positive fourth power, and using D2V=D2UD^2V=D^2U, gives

sup⁡{−U≥η}∣D2U∣≤C(n,k,η).(53)\sup_{\{-U\ge\eta\}}|D^2U|\le C(n,k,\eta). \tag*{(53)}

To express this bound in terms of distance, fix x∈Ωx \in\Omega and 0<r<d(x)0 < r < d(x). On Br(x)B_r(x) introduce

q(y)=an,k(∣y−x∣2−r2),an,k=12(nk)1/k.(54)q(y) = a_{n,k}(|y-x|^2-r^2), \qquad a_{n,k} = \frac{1}{2\binom{n}{k}^{1/k}}. \tag*{(54)}

All eigenvalues of D2qD^2q equal 2an,k>02a_{n,k}>0, and hence

σk(D2q)=(nk)(2an,k)k=1.\sigma_k(D^2q)=\binom{n}{k}(2a_{n,k})^k=1.

On ∂Br(x)\partial B_r(x), we have U≤0=qU \le0=q. Both functions are admissible solutions of the same equation, so the comparison principle gives U≤qU \le q throughout the ball.

Evaluate at xx and let r↑d(x)r \uparrow d(x). This yields

−U(x)≥an,kd(x)2.(55)-U(x) \ge a_{n,k}d(x)^2. \tag*{(55)}

Taking η=an,kδ2\eta=a_{n,k}\delta^2 in (53) bounds the Hessian on {d≥δ}\{d\ge\delta\}.

The bound for UU follows from (48). The gradient bound follows from (50), because (55) places {d≥δ}\{d\ge\delta\} inside {−U≥η}\{-U\ge\eta\}.

Scaling back proves (47).

Select the admissible branch

Let uu satisfy the hypotheses of Theorem 1.2. Introduce the rescaled function uR(y)=R−2u(Ry)u_R(y)=R^{-2}u(Ry) on B1B_1. Its derivatives satisfy

DuR(y)=R−1Du(Ry),D2uR(y)=D2u(Ry).Du_R(y)=R^{-1}Du(Ry), \qquad D^2u_R(y)=D^2u(Ry).

Thus the equation is unchanged, and the lower bound becomes uR(y)≥b∣y∣2−B/R2u_R(y)\ge b|y|^2-B/R^2. The hypothesis on BB gives 0≤B/R2≤b/40\le B/R^2\le b/4.

It therefore suffices to prove (9) when R=1R=1. In this normalization, u(0)=0u(0)=0, u(x)≥b∣x∣2−Bu(x)\ge b|x|^2-B, and 0≤B≤b/40\le B\le b/4.

The growth bound gives

u(0)=0,lim inf⁡∣x∣→1u(x)≥b−B≥3b4>0.(56)u(0)=0, \qquad\liminf_{|x|\to1}u(x)\ge b-B\ge\frac{3b}{4}>0. \tag*{(56)}

Choose r0∈(1/2,1)r_0\in(1/2,1). The quadratic lower bound is positive on B1∖B‾r0B_1\setminus\overline{B}_{r_0}, whereas u(0)=0u(0)=0. A minimum on B‾r0\overline{B}_{r_0} is therefore also a minimum in B1B_1.

At such a point x∗∈B1x_*\in B_1, the second derivative test gives

D2u(x∗)≥0.(57)D^2u(x_*)\ge0. \tag*{(57)}

Since σk(D2u(x∗))=1\sigma_k(D^2u(x_*))=1, at least kk of its nonnegative eigenvalues are positive. Hence σj(D2u(x∗))>0\sigma_j(D^2u(x_*))>0 for 1≤j≤k1\le j\le k, so D2u(x∗)∈ΓkD^2u(x_*)\in\Gamma_k.

We now use a continuity argument. Fix x∈B1x\in B_1 and follow the Hessian along the segment from x∗x_* to xx:

H(t)=D2u((1−t)x∗+tx),0≤t≤1.(58)H(t)=D^2u((1-t)x_*+tx), \qquad0\le t\le1. \tag*{(58)}

Define

G={t∈[0,1]:H(t)∈Γk}.(59)G=\{t\in[0,1]:H(t)\in\Gamma_k\}. \tag*{(59)}

Since H(0)∈ΓkH(0)\in\Gamma_k, we have 0∈G0\in G. If t0∈Gt_0\in G, then every number σj(H(t0))\sigma_j(H(t_0)), 1≤j≤k1\le j\le k, is positive. Continuity preserves these finitely many strict inequalities for tt sufficiently close to t0t_0, so GG is relatively open. To prove closedness, we use Maclaurin’s inequalities: for λ∈Γk\lambda\in\Gamma_k,

(σs−1(λ)(ns−1))1/(s−1)≥(σs(λ)(ns))1/s,2≤s≤k.(60)\left(\frac{\sigma_{s-1}(\lambda)}{\binom{n}{s-1}}\right)^{1/(s-1)} \ge\left(\frac{\sigma_s(\lambda)}{\binom{n}{s}}\right)^{1/s}, \qquad2 \le s \le k. \tag*{(60)}

Iterating these inequalities and using σk(H(t))=1\sigma_k(H(t))=1, we obtain

σj(H(t))≥(nj)(nk)−j/k>0,t∈G,1≤j<k.(61)\sigma_j(H(t)) \ge\binom{n}{j}\binom{n}{k}^{-j/k}>0, \qquad t \in G,\quad1 \le j<k. \tag*{(61)}

The lower bounds in (61) are independent of tt. Let tm∈Gt_m \in G converge to t∈[0,1]t \in[0,1]. Continuity gives

σj(H(t))=lim⁡m→∞σj(H(tm))≥(nj)(nk)−j/k>0,1≤j<k.\sigma_j(H(t))=\lim_{m\to\infty}\sigma_j(H(t_m))\ge\binom{n}{j}\binom{n}{k}^{-j/k}>0,\qquad1\le j<k.

The equation also gives σk(H(t))=1\sigma_k(H(t))=1. Thus all defining inequalities of Γk\Gamma_k remain strict at H(t)H(t), which proves t∈Gt\in G and the relative closedness of GG.

Since [0,1][0,1] is connected and GG is nonempty, open, and closed in [0,1][0,1], we conclude that G=[0,1]G=[0,1]. In particular, D2u(x)=H(1)∈ΓkD^2u(x)=H(1)\in\Gamma_k.

As x∈B1x\in B_1 was arbitrary, uu is admissible throughout the ball. We can now apply admissible comparison and the interior estimates.

Describe the regular sublevel

We first explain the curvature of an admissible level set. For a symmetric matrix HH and 0≤j≤n−10\le j\le n-1, its jj-th Newton tensor is the symmetric matrix Tj(H)T_j(H) defined by

Dσj+1(H)[K]=∑p,q=1nTjpq(H)Kpq,K=KT.D\sigma_{j+1}(H)[K]=\sum_{p,q=1}^{n}T_j^{pq}(H)K_{pq},\qquad K=K^T.

If H∈ΓkH\in\Gamma_k and E=ν⊥E=\nu^\perp for a unit vector ν\nu, then

σj(H∣E)=⟨Tj(H)ν,ν⟩>0,1≤j≤k−1.(62)\sigma_j(H|_E)=\langle T_j(H)\nu,\nu\rangle>0,\qquad1\le j\le k-1. \tag*{(62)}

To see the equality, choose orthonormal coordinates with last vector ν\nu. Differentiating σj+1(H)\sigma_{j+1}(H) with respect to HnnH_{nn} leaves the sum of all jj-dimensional principal minors of H∣EH|_E. Since H∈Γk⊂Γj+1H\in\Gamma_k\subset\Gamma_{j+1}, ellipticity makes Tj(H)T_j(H) positive definite. Evaluating its quadratic form on ν\nu proves positivity in (62).

By Sard’s theorem, choose a regular value t∈(b/4,b/2)t\in(b/4,b/2). Since u(0)=0<tu(0)=0<t, the set {u<t}\{u<t\} has a component Ωt\Omega_t containing the origin.

For x∈Ωtx\in\Omega_t, the lower bound gives b∣x∣2−B≤u(x)<tb|x|^2-B\le u(x)<t. Using t<b/2t<b/2 and B≤b/4B\le b/4, we obtain

Ωt⊂B(t+B)/b⊂B3/2,Dt:=diam⁡Ωt≤3.(63)\Omega_t\subset B_{\sqrt{(t+B)/b}}\subset B_{\sqrt{3}/2},\qquad D_t:=\operatorname{diam}\Omega_t\le\sqrt{3}. \tag*{(63)}

Its boundary is smooth and compactly contained in B1B_1. The outward normal and second fundamental form are

ν=Du∣Du∣,II∂Ωt=D2u∣T∂Ωt∣Du∣.(64)\nu=\frac{Du}{|Du|},\qquad\mathrm{II}_{\partial\Omega_t}=\frac{D^2u|_{T\partial\Omega_t}}{|Du|}. \tag*{(64)}

The second identity follows by differentiating ν=Du/∣Du∣\nu=Du/|Du| in tangent directions. The term from differentiating ∣Du∣−1|Du|^{-1} vanishes when paired with a tangent vector.

Equation (62) places the tangential Hessian in Γk−1\Gamma_{k-1}. Division by the positive number ∣Du∣|Du| preserves this cone, so the principal curvatures belong to Γk−1\Gamma_{k-1}. The function U=u−tU=u-t is admissible, satisfies σk(D2U)=1\sigma_k(D^2U)=1, and vanishes on ∂Ωt\partial\Omega_t. It therefore satisfies the hypotheses of Proposition 3.1 and Lemma 4.1.

Apply the domain estimate

Set dt(x)=dist⁡(x,∂Ωt)d_t(x)=\operatorname{dist}(x,\partial\Omega_t). At the origin, the boundary estimate and (63) imply

b4<t=−U(0)≤C0Dt2−βdt(0)β≤C031−β/2dt(0)β.(65)\frac{b}{4}<t=-U(0)\le C_0D_t^{2-\beta}d_t(0)^\beta\le C_0 3^{1-\beta/2}d_t(0)^\beta. \tag*{(65)}

Define constants depending only on n,k,bn,k,b by

γ=min⁡{14,(b4C031−β/2)1/β},ρ=γ2.(66)\gamma=\min\left\{\frac{1}{4},\left(\frac{b}{4C_0 3^{1-\beta/2}}\right)^{1/\beta}\right\},\qquad\rho=\frac{\gamma}{2}. \tag*{(66)}

Then dt(0)≥γd_t(0)\ge\gamma, so Bγ⊂ΩtB_\gamma\subset\Omega_t. Indeed, a segment from the origin to a point outside Ωt\Omega_t must first meet ∂Ωt\partial\Omega_t, at distance at least dt(0)d_t(0).

Distance to a fixed closed set is 1-Lipschitz. Consequently, for x∈Bρx\in B_\rho we have

dt(x)≥dt(0)−∣x∣≥ρ≥ρ3Dt.(67)d_t(x)\ge d_t(0)-|x|\ge\rho\ge\frac{\rho}{\sqrt{3}}D_t. \tag*{(67)}

Apply Lemma 4.1 with δ=ρ/3\delta=\rho/\sqrt{3}. Since Dt≤3D_t\le\sqrt{3} and u=U+tu=U+t with t<b/2t<b/2, we obtain

sup⁡Bρ(∣u∣+∣Du∣+∣D2u∣)≤C(n,k,b).(68)\sup_{B_\rho}(|u|+|Du|+|D^2u|)\le C(n,k,b). \tag*{(68)}

Apply (68) to uRu_R and use the derivative identities from § 4.2. The three terms become R−2∣u∣R^{-2}|u|, R−1∣Du∣R^{-1}|Du|, and ∣D2u∣|D^2u| on BρRB_{\rho R}, which proves Theorem 1.2.

Proof of the Liouville theorem

We apply Theorem 1.2 on expanding balls. The resulting Hessian bound makes the equation uniformly elliptic. The Evans–Krylov estimate then forces the Hessian to be constant.

Obtain a global Hessian bound

Let uu satisfy the hypotheses of Theorem 1.1. Normalize its value at the origin by setting

v=u−u(0),B~=B+u(0)≥0.(69)v=u-u(0),\qquad\widetilde{B}=B+u(0)\ge0. \tag*{(69)}

The last inequality follows from (5) at x=0x=0. For every sufficiently large RR, we have

v(x)≥b∣x∣2−B~(x∈BR),0≤B~≤bR24.(70)v(x)\ge b|x|^2-\widetilde{B}\quad(x\in B_R),\qquad0\le\widetilde{B}\le\frac{bR^2}{4}. \tag*{(70)}

Theorem 1.2 therefore gives

sup⁡BρR∣D2u∣≤C(n,k,b),sup⁡Rn∣D2u∣≤M:=C(n,k,b).(71)\sup_{B_{\rho R}}|D^2u|\le C(n,k,b),\qquad\sup_{\mathbb{R}^n}|D^2u|\le M:=C(n,k,b). \tag*{(71)}

For each fixed x∈Rnx\in\mathbb{R}^n, choose RR large enough that B~≤bR2/4\widetilde{B}\le bR^2/4 and ∣x∣<ρR|x|<\rho R. The first estimate then bounds ∣D2u(x)∣|D^2u(x)| by the same constant MM.

The branch argument in § 4.2 applies on each of these balls. Since the balls exhaust Rn\mathbb{R}^n, it also gives D2u∈ΓkD^2u\in\Gamma_k everywhere.

Rescale the uniformly elliptic equation

Set F(H)=σk(H)1/kF(H)=\sigma_k(H)^{1/k} and introduce the set containing all Hessian values:

K={H∈Γ‾k: σk(H)=1, ∣H∣≤M}.(72)\mathcal{K}=\{H\in\overline{\Gamma}_k:\ \sigma_k(H)=1,\ |H|\leq M\}. \tag*{(72)}

The Maclaurin inequalities (4.14) extend by continuity to Γ‾k\overline{\Gamma}_k. Since σk(H)=1\sigma_k(H)=1 on K\mathcal{K}, they give σj(H)>0\sigma_j(H)>0 for every H∈KH\in\mathcal{K} and 1≤j<k1\leq j<k.

Thus K\mathcal{K} is a compact subset of Γk\Gamma_k. The derivative matrix is

Fij(H)=1kσk(H)1/k−1Tk−1ij(H).F^{ij}(H)=\frac{1}{k}\sigma_k(H)^{1/k-1}T_{k-1}^{ij}(H).

It is positive definite on Γk\Gamma_k and continuous on K\mathcal{K}. Compactness therefore gives constants 0<λ≤Λ0<\lambda\leq\Lambda, depending only on nn, kk, MM, such that λI≤(Fij(H))≤ΛI\lambda I\leq(F^{ij}(H))\leq\Lambda I.

For completeness, extend the operator to all symmetric matrices by

F~(H)=inf⁡A∈K{1+∑i,jFij(A)(Hij−Aij)}.(73)\widetilde{F}(H)=\inf_{A\in\mathcal{K}}\left\{1+\sum_{i,j}F^{ij}(A)(H_{ij}-A_{ij})\right\}. \tag*{(73)}

The matrices AA and (Fij(A))(F^{ij}(A)) are bounded on K\mathcal{K}, so this infimum is finite for every HH. An infimum of affine functions is concave.

For every positive-semidefinite matrix PP, each affine function in (5.5) increases by a number between λtr⁡P\lambda\operatorname{tr}P and Λtr⁡P\Lambda\operatorname{tr}P.

Taking infima preserves these bounds, so

λtr⁡P≤F~(H+P)−F~(H)≤Λtr⁡P,P≥0.\lambda\operatorname{tr}P\leq\widetilde{F}(H+P)-\widetilde{F}(H)\leq\Lambda\operatorname{tr}P,\qquad P\geq0.

Thus F~\widetilde{F} is uniformly elliptic on all symmetric matrices. For H,A∈KH,A\in\mathcal{K}, concavity of FF gives

1=F(H)≤F(A)+∑i,jFij(A)(Hij−Aij).1=F(H)\leq F(A)+\sum_{i,j}F^{ij}(A)(H_{ij}-A_{ij}).

Every expression in the infimum is therefore at least one. Choosing A=HA=H attains one, which proves F~(H)=1\widetilde{F}(H)=1 on K\mathcal{K}.

Subtract the affine Taylor polynomial at the origin and rescale:

wR(y)=u(Ry)−u(0)−RDu(0)⋅yR2,y∈B1.(74)w_R(y)=\frac{u(Ry)-u(0)-RDu(0)\cdot y}{R^2},\qquad y\in B_1. \tag*{(74)}

The normalization gives wR(0)=0w_R(0)=0 and DwR(0)=0Dw_R(0)=0. Also, D2wR(y)=D2u(Ry)D^2w_R(y)=D^2u(Ry), so (5.3) bounds its Hessian independently of RR.

Integrating along the segment from 00 to yy yields

∣DwR(y)∣≤C(n)M∣y∣,∣wR(y)∣≤C(n)M∣y∣2.|Dw_R(y)|\leq C(n)M|y|,\qquad|w_R(y)|\leq C(n)M|y|^2.

Hence ∥wR∥C1,1(B1)≤C(n)M\|w_R\|_{C^{1,1}(B_1)}\leq C(n)M. Its Hessian takes values in K\mathcal{K}, and therefore F~(D2wR)=1\widetilde{F}(D^2w_R)=1.

To use the convex formulation of [2 Theorem 1], consider −wR-w_R and the operator H↦1−F~(−H)H\mapsto1-\widetilde{F}(-H). This operator is convex and uniformly elliptic.

Its value on D2(−wR)D^2(-w_R) is zero. The cited Evans–Krylov estimate therefore gives constants ζ∈(0,1)\zeta\in(0,1) and C>0C>0, independent of RR, such that

[D2wR]Cζ(B1/2)≤C,[D2u]Cζ(BR/2)≤CR−ζ.(75)[D^2w_R]_{C^\zeta(B_{1/2})}\leq C,\qquad[D^2u]_{C^\zeta(B_{R/2})}\leq CR^{-\zeta}. \tag*{(75)}

For fixed x,y∈Rnx,y\in\mathbb{R}^n, both points lie in BR/2B_{R/2} for all sufficiently large RR. The second estimate then gives

∣D2u(x)−D2u(y)∣≤CR−ζ∣x−y∣ζ.|D^2u(x)-D^2u(y)|\leq CR^{-\zeta}|x-y|^\zeta.

Letting R→∞R \to\infty makes the right-hand side zero. Thus D2u(x)=D2u(y)D^2u(x) = D^2u(y) for every pair of points. Write this constant Hessian as AA. Integration gives u(x)=12xTAx+ℓ⋅x+cu(x) = \frac{1}{2}x^T A x + \ell\cdot x + c, and the equation gives σk(A)=1\sigma_k(A) = 1.

Finally, take x=rex = re in (1.5), where ∣e∣=1|e| = 1. Substitution of the quadratic expression and division by r2r^2 give

12eTAe+ℓ⋅er+cr2≥b−Br2.\frac{1}{2}e^T A e + \frac{\ell\cdot e}{r} + \frac{c}{r^2} \ge b - \frac{B}{r^2}.

Letting r→∞r \to\infty removes the affine and constant terms, yielding

12eTAe≥b.(76)\frac{1}{2}e^T A e \ge b. \tag*{(76)}

Hence A≥2bIA \ge2bI, proving Theorem 1.1.

Viscosity solutions

We specify the viscosity formulation and prove automatic admissibility. Approximation by smooth solutions then gives the interior estimate in Theorem 1.4 and the Liouville theorem in Theorem 1.3.

The viscosity formulation

Let Ω⊂Rn\Omega\subset\mathbb{R}^n be open. A proper upper semicontinuous function u:Ω→[−∞,∞)u : \Omega\to[-\infty,\infty) is kk-convex if every C2C^2 function touching uu from above has Hessian in Γ‾k\overline{\Gamma}_k at the contact point.

Here proper means that uu is not identically −∞-\infty on any component. We denote this class by Φk(Ω)\Phi_k(\Omega) and use its upper semicontinuous representatives throughout. For weak functions, the admissibility condition means membership in Φk(Ω)\Phi_k(\Omega).

A continuous function uu is a viscosity solution of σk(D2u)=1\sigma_k(D^2u) = 1 if every C2C^2 upper touching test ϕ\phi satisfies

σk(D2ϕ(x))≥1,(77)\sigma_k(D^2\phi(x)) \ge1, \tag*{(77)}

and every C2C^2 lower touching test ϕ\phi with D2ϕ(x)∈Γ‾kD^2\phi(x) \in\overline{\Gamma}_k satisfies

σk(D2ϕ(x))≤1.(78)\sigma_k(D^2\phi(x)) \le1. \tag*{(78)}

Touching from above means that u−ϕu - \phi has a local maximum zero at the contact point; touching from below means a local minimum zero. The lower tests are taken on the elliptic branch.

The following lemma shows that the upper-test inequality already forces admissibility, giving the usual formulation in [16 Section 2].

Lemma 6.1 (Automatic viscosity admissibility). Let u∈C(Ω)u \in C(\Omega) satisfy (6.1) for every C2C^2 upper touching test. Then every such test has Hessian in Γ‾k\overline{\Gamma}_k at its contact point. In particular, u∈Φk(Ω)u \in\Phi_k(\Omega).

Proof. Fix an upper touching test ϕ\phi at x0x_0, and put A=D2ϕ(x0)A = D^2\phi(x_0). Choose T>0T > 0 such that A+TI>0A + TI > 0, and define

G={t∈[0,T]:A+tI∈Γ‾k}.(79)G = \{t \in[0,T] : A + tI \in\overline{\Gamma}_k\}. \tag*{(79)}

The matrix A+TIA + TI is positive definite, so T∈GT \in G. For t0∈Gt_0 \in G, continuity keeps all the finitely many inequalities σj(A+tI)>0\sigma_j(A + tI) > 0 valid when tt is sufficiently close to t0t_0. Thus GG is relatively open in [0,T][0,T].

For every t∈[0,T]t \in[0,T], introduce

ϕt(x)=ϕ(x)+t2∣x−x0∣2.\phi_t(x) = \phi(x) + \frac{t}{2}|x - x_0|^2.

The added term is nonnegative and vanishes at x0x_0, so ϕt\phi_t still touches uu from above. Applying (77) to this test gives σk(A+tI)≥1\sigma_k(A+tI)\ge1 for the whole interval.

For t∈Gt\in G, Maclaurin’s inequality (60) yields

σj(A+tI)≥(nj)(nk)−j/k>0,1≤j<k.(80)\sigma_j(A+tI)\ge\binom{n}{j}\binom{n}{k}^{-j/k}>0,\qquad1\le j<k. \tag*{(80)}

To prove closedness, take tm∈Gt_m\in G with tm→t∈[0,T]t_m\to t\in[0,T]. Continuity of each symmetric polynomial preserves the positive lower bounds in (80).

Together with σk(A+tI)≥1\sigma_k(A+tI)\ge1, these bounds imply A+tI∈ΓkA+tI\in\Gamma_k. Hence t∈Gt\in G, and GG is relatively closed.

Connectedness gives G=[0,T]G=[0,T], so A∈ΓkA\in\Gamma_k. Since the upper touching test was arbitrary, u∈Φk(Ω)u\in\Phi_k(\Omega). □\square

This argument does not use the growth condition. For the classical theorems, admissibility follows instead from the growth condition and continuity of the Hessian, as in Section 4.2.

A margin in the smooth estimate

We need a small margin in the growth condition to accommodate the error in approximation. The proof of Theorem 1.2 gives the following variant with constants depending on the same parameters.

Lemma 6.2. Let u∈C∞(BR)u\in C^\infty(B_R) satisfy σk(D2u)=1\sigma_k(D^2u)=1 and u(0)=0u(0)=0. Suppose b>0b>0 and

u(x)≥b∣x∣2−B,0≤B≤bR22.(81)u(x)\ge b|x|^2-B,\qquad0\le B\le\frac{bR^2}{2}. \tag*{(81)}

There are constants η=η(n,k,b)∈(0,1)\eta=\eta(n,k,b)\in(0,1) and C=C(n,k,b)C=C(n,k,b) such that

sup⁡BηR(∣u∣R2+∣Du∣R+∣D2u∣)≤C.(82)\sup_{B_{\eta R}}\left(\frac{|u|}{R^2}+\frac{|Du|}{R}+|D^2u|\right)\le C. \tag*{(82)}

Proof. Rescale to R=1R=1 and choose r0∈(1/2,1)r_0\in(1/\sqrt{2},1). For r0≤∣x∣<1r_0\le|x|<1, the lower bound gives

u(x)≥br02−B≥b(r02−1/2)>0=u(0).u(x)\ge br_0^2-B\ge b(r_0^2-1/2)>0=u(0).

Hence a minimum on B‾r0\overline{B}_{r_0} is an interior minimum in B1B_1. The continuity argument in Section 4.2 then gives D2u∈ΓkD^2u\in\Gamma_k throughout B1B_1.

Choose a regular value t∈(b/8,b/4)t\in(b/8,b/4), and let Ωt\Omega_t be the component of {u<t}\{u<t\} containing the origin. Then

Ωt⊂B(t+B)/b⊂B3/2,diam⁡Ωt≤3.(83)\Omega_t\subset B_{\sqrt{(t+B)/b}}\subset B_{\sqrt{3}/2},\qquad\operatorname{diam}\Omega_t\le\sqrt{3}. \tag*{(83)}

As in Section 4, its boundary is (k−1)(k-1)-convex. For U=u−tU=u-t and dt(x)=dist⁡(x,∂Ωt)d_t(x)=\operatorname{dist}(x,\partial\Omega_t), Proposition 3.1 gives

b8<t=−U(0)≤C031−β/2dt(0)β.(84)\frac{b}{8}<t=-U(0)\le C_0 3^{1-\beta/2}d_t(0)^\beta. \tag*{(84)}

Choose

γ0=min⁡{14,(b8C031−β/2)1/β},η=γ02.\gamma_0=\min\left\{\frac14,\left(\frac{b}{8C_0 3^{1-\beta/2}}\right)^{1/\beta}\right\},\qquad\eta=\frac{\gamma_0}{2}.

Then (84) gives dt(0)≥γ0d_t(0)\ge\gamma_0. The segment argument following (66) gives Bγ0⊂ΩtB_{\gamma_0}\subset\Omega_t.

Thus Bη⊂ΩtB_\eta\subset\Omega_t. Since distance is 1-Lipschitz and η=γ0/2\eta=\gamma_0/2, every x∈Bηx\in B_\eta satisfies dt(x)≥ηd_t(x)\ge\eta. Writing Dt=diam⁡Ωt≤3D_t = \operatorname{diam}\Omega_t \le\sqrt{3}, we therefore have dt(x)≥(η/3)Dtd_t(x) \ge(\eta/\sqrt{3})D_t on BηB_\eta. Lemma 4.1 applies with δ=η/3\delta= \eta/\sqrt{3}.

Finally, u=U+tu = U + t and t<b/4t < b/4, so the resulting bounds for UU, DUDU, and D2UD^2U give (6.6) when R=1R = 1. The change of variables x=Ryx = Ry restores the stated powers of RR.

Approximation by solutions of the equation

Proof of Theorem 1.4. By Lemma 6.1, the continuous viscosity solution uu is admissible. Rescale to R=1R = 1 and set r=3/2r = \sqrt{3}/2. The smaller ball has closure inside B1B_1, so uu is continuous on Br‾\overline{B_r} and its boundary values can be approximated uniformly.

Choose smooth functions φm\varphi_m on ∂Br\partial B_r with

εm:=∥φm−u∥L∞(∂Br)⟶0.(85)\varepsilon_m := \|\varphi_m-u\|_{L^\infty(\partial B_r)} \longrightarrow0. \tag*{(85)}

The sphere ∂Br\partial B_r is connected, and its principal curvatures are all 1/r1/r. Consequently,

σk−1(κ)=(n−1k−1)r1−k>0.\sigma_{k-1}(\kappa) = \binom{n-1}{k-1}r^{1-k} > 0.

The smooth Dirichlet theorem of Caffarelli, Nirenberg, and Spruck [14 Theorem 3] therefore supplies admissible solutions, smooth up to the boundary, satisfying

σk(D2um)=1in Br,um=φmon ∂Br.(86)\sigma_k(D^2u_m) = 1 \quad\text{in } B_r,\qquad u_m = \varphi_m \quad\text{on } \partial B_r. \tag*{(86)}

Maclaurin’s inequality (60) bounds every lower symmetric polynomial of D2umD^2u_m away from zero. Continuity preserves these bounds on ∂Br\partial B_r.

For fixed mm, it follows that D2um(Br‾)⋐ΓkD^2u_m(\overline{B_r}) \Subset\Gamma_k. Thus there exists δm>0\delta_m > 0 such that perturbations by ±2δI\pm2\delta I preserve admissibility whenever 0<δ<δm0 < \delta< \delta_m.

We explain the comparison that transfers the boundary error into the ball. Define

qm+=um+εm−δ(∣x∣2−r2),qm−=um−εm+δ(∣x∣2−r2).q_m^+ = u_m + \varepsilon_m - \delta\left(|x|^2-r^2\right),\qquad q_m^- = u_m - \varepsilon_m + \delta\left(|x|^2-r^2\right).

Take 0<δ<δm0 < \delta< \delta_m. For any symmetric matrix AA,

ddtσk(A+tI)=(n−k+1)σk−1(A+tI)>0when A+tI∈Γk.\frac{d}{dt}\sigma_k(A+tI) = (n-k+1)\sigma_{k-1}(A+tI) > 0 \quad\text{when } A+tI \in\Gamma_k.

Apply this identity on the admissible segments from D2umD^2u_m to D2um±2δID^2u_m \pm2\delta I. Since σk(D2um)=1\sigma_k(D^2u_m) = 1, it gives σk(D2qm+)<1<σk(D2qm−)\sigma_k(D^2q_m^+) < 1 < \sigma_k(D^2q_m^-).

The boundary approximation gives qm−≤u≤qm+q_m^- \le u \le q_m^+ on ∂Br\partial B_r. If u−qm+u-q_m^+ had a positive interior maximum, a constant translate of qm+q_m^+ would violate the upper-test inequality. Similarly, a positive interior maximum of qm−−uq_m^- - u would provide an admissible lower touching test with σk>1\sigma_k > 1, contradicting (78). Hence qm−≤u≤qm+q_m^- \le u \le q_m^+ in the ball.

Letting δ↓0\delta\downarrow0 proves

∥um−u∥L∞(Br)≤εm.(87)\|u_m-u\|_{L^\infty(B_r)} \le\varepsilon_m. \tag*{(87)}

Normalize at the center by putting vm=um−um(0)v_m = u_m-u_m(0). The equation is unchanged by subtraction of this constant. Since u(0)=0u(0)=0, we have ∣um(0)∣≤εm|u_m(0)| \le\varepsilon_m, and hence ∥vm−u∥L∞(Br)≤2εm\|v_m-u\|_{L^\infty(B_r)} \le2\varepsilon_m.

Combining this estimate with the lower bound for uu, we obtain

vm(x)≥b∣x∣2−B−2εm,B+2εm≤br22=3b8.(88)v_m(x) \ge b|x|^2-B-2\varepsilon_m,\qquad B+2\varepsilon_m \le\frac{br^2}{2}=\frac{3b}{8}. \tag*{(88)}

for all sufficiently large mm: the assumption B≤b/4B \le b/4 leaves a margin 3b/8−b/4=b/83b/8-b/4=b/8, which absorbs 2εm2\varepsilon_m. Lemma 6.2, applied on BrB_r, now gives

sup⁡Bηr(∣vm∣+∣Dvm∣+∣D2vm∣)≤C(n,k,b).(89)\sup_{B_{\eta r}} \left(|v_m|+|Dv_m|+|D^2v_m|\right) \le C(n,k,b). \tag*{(89)}

The Hessians lie in a common compact subset of Γk\Gamma_k, by (60) and σk(D2vm)=1\sigma_k(D^2v_m)=1. The construction in § 5 gives a common concave, uniformly elliptic extension of σk1/k\sigma_k^{1/k} on this Hessian range.

The Evans–Krylov estimate gives a uniform C2,ζC^{2,\zeta} bound on Bηr/2B_{\eta r/2} for some ζ>0\zeta>0. Compactness gives a subsequence converging in C2C^2 on every compact subset of this ball.

Its limit is uu, because vm→uv_m\to u uniformly. Continuity of σk\sigma_k therefore gives σk(D2u)=1\sigma_k(D^2u)=1 classically on Bηr/2B_{\eta r/2}.

Pass to the limit at each point of this smaller ball in (89), and then take the supremum. This gives

sup⁡Bηr/2(∣u∣+∣Du∣+∣D2u∣)≤C(n,k,b).\sup_{B_{\eta r/2}} \left(|u|+|Du|+|D^2u|\right) \le C(n,k,b).

For smoothness, differentiate the equation of the smooth approximants. With F=σk1/kF=\sigma_k^{1/k}, each coordinate derivative satisfies

∑i,jFij(D2vm) ∂ij(∂ℓvm)=0,1≤ℓ≤n.\sum_{i,j}F^{ij}(D^2v_m)\,\partial_{ij}(\partial_\ell v_m)=0,\qquad1\le\ell\le n.

The coefficients have uniform positive ellipticity bounds and uniform C0,ζC^{0,\zeta} norms on smaller balls. Interior Schauder estimates therefore give uniform C3,ζC^{3,\zeta} bounds for vmv_m there.

Iterating on compactly contained balls gives bounds of every order, so u∈C∞(Bηr/2)u\in C^\infty(B_{\eta r/2}). Rescaling gives (1.15) with ρ=ηr/2\rho=\eta r/2.

The Liouville theorem

Proof of Theorem 1.3. Let uu be the given entire viscosity solution. Normalize by v=u−u(0)v=u-u(0) and put B~=B+u(0)≥0\widetilde B=B+u(0)\ge0. For all sufficiently large RR, we have

v(x)≥b∣x∣2−B~,0≤B~≤bR24.(90)v(x)\ge b|x|^2-\widetilde B,\qquad0\le\widetilde B\le\frac{bR^2}{4}. \tag*{(90)}

Theorem 1.4 gives u∈C∞(BρR)u\in C^\infty(B_{\rho R}), where ρ>0\rho>0 depends only on n,k,bn,k,b. These balls exhaust Rn\mathbb{R}^n, so u∈C∞(Rn)u\in C^\infty(\mathbb{R}^n).

The quadratic lower bound is unchanged. Theorem 1.1 therefore gives (1.12).

Admissible weak solutions

We use Trudinger’s definition of admissible weak solutions by approximation. His viscosity characterization gives the two weak-solution theorems. We then explain the integration-by-parts identities satisfied by these solutions.

Trudinger's definition

Following [15 p. 1252, (1.4)], a function u∈C(Ω)u \in C(\Omega) is an admissible weak solution of σk(D2u)=1\sigma_k(D^2u)=1 if there exist um∈C2(Ω)u_m \in C^2(\Omega) satisfying

D2um∈Γ‾k in Ω,um→ulocally uniformly,D^2u_m \in\overline{\Gamma}_k \text{ in } \Omega,\qquad u_m \to u \quad\text{locally uniformly},
σk(D2um)→1in Lloc1(Ω).(91)\sigma_k(D^2u_m) \to1 \quad\text{in } L^1_{\mathrm{loc}}(\Omega). \tag*{(91)}

The approximants are kk-admissible in the weak sense, including points where their Hessians lie on the boundary of Γk\Gamma_k. Their locally uniform limit belongs to Φk(Ω)\Phi_k(\Omega).

The Liouville theorem and the interior estimate

For a continuous right-hand side, Trudinger characterizes admissible weak solutions by the viscosity inequalities [15 p. 1255]. We apply this characterization with right-hand side equal to one.

Proof of Theorem 1.6. The function uu is continuous by definition. Trudinger's characterization makes it a viscosity solution of σk(D2u)=1\sigma_k(D^2u)=1 throughout Rn\mathbb{R}^n. The lower bound is unchanged, so all hypotheses of Theorem 1.3 hold. That theorem gives (1.22).

Proof of Theorem 1.7. Trudinger's characterization makes uu a continuous viscosity solution on BRB_R. The value u(0)=0u(0)=0 and the lower bound with 0≤B≤bR2/40 \le B \le bR^2/4 are unchanged. Apply Theorem 1.4 to obtain u∈C∞(BρR)u \in C^\infty(B_{\rho R}) and (1.25). Its constants depend only on n,k,bn,k,b.

Integration by parts and Hessian measures

For smooth uu, write Tk−1ij(D2u)T^{ij}_{k-1}(D^2u) for the entries of the Newton tensor Dσk(D2u)D\sigma_k(D^2u), as in § 4. Its divergence-free property and the homogeneity of σk\sigma_k give

∑i∂iTk−1ij(D2u)=0,∑i,jTk−1ij(D2u)uij=kσk(D2u).(92)\sum_i \partial_i T^{ij}_{k-1}(D^2u)=0,\qquad\sum_{i,j}T^{ij}_{k-1}(D^2u)u_{ij}=k\sigma_k(D^2u). \tag*{(92)}

The second identity follows by differentiating σk(tH)=tkσk(H)\sigma_k(tH)=t^k\sigma_k(H) at t=1t=1. The first identity removes the derivative of the Newton tensor in the product rule, giving

∑i,j∂i(Tk−1ij(D2u)uj)=kσk(D2u).\sum_{i,j}\partial_i\left(T^{ij}_{k-1}(D^2u)u_j\right)=k\sigma_k(D^2u).

Multiply by ϕ∈Cc∞(Ω)\phi\in C^\infty_c(\Omega) and integrate first in xix_i. Then integrate in xjx_j; symmetry of the Newton tensor gives the corresponding divergence identity in that index. Thus

∫Ωϕσk(D2u) dx=−1k∑i,j∫ΩϕiTk−1ij(D2u)uj dx=1k∑i,j∫ΩuTk−1ij(D2u)ϕij dx.(93)\begin{aligned} \int_\Omega\phi\sigma_k(D^2u)\,dx &=-\frac{1}{k}\sum_{i,j}\int_\Omega\phi_iT^{ij}_{k-1}(D^2u)u_j\,dx\\ &=\frac{1}{k}\sum_{i,j}\int_\Omega uT^{ij}_{k-1}(D^2u)\phi_{ij}\,dx. \tag*{(93)} \end{aligned}

Both boundary terms vanish because ϕ\phi has compact support. The divergence identities and their mixed versions are discussed in [TW02, Section 2].

The identity extends to C2C^2 functions by local smooth approximation. Apply it to each umu_m in (91). The convergence of the Hessian densities in Lloc1L^1_{\mathrm{loc}} yields

1klim⁡m→∞∑i,j∫ΩumTk−1ij(D2um)ϕij dx=∫Ωϕ dx,ϕ∈Cc∞(Ω).(94)\frac{1}{k}\lim_{m\to\infty}\sum_{i,j}\int_{\Omega}u_mT_{k-1}^{ij}(D^2u_m)\phi_{ij}\,dx=\int_{\Omega}\phi\,dx,\qquad\phi\in C^\infty_c(\Omega). \tag*{(94)}

Only the contracted integral is passed to the limit. No pointwise second derivatives of uu or separate limits of the products umTk−1ij(D2um)u_mT_{k-1}^{ij}(D^2u_m) are needed.

For k=2k=2, the Newton tensor has the explicit form T1ij(D2v)=Δv δij−vijT_1^{ij}(D^2v)=\Delta v\,\delta_{ij}-v_{ij}. For smooth vv, integrate the two terms separately:

∫vΔv Δϕ=−∫∣Dv∣2Δϕ−∫vDv⋅DΔϕ,\int v\Delta v\,\Delta\phi=-\int|Dv|^2\Delta\phi-\int vDv\cdot D\Delta\phi,
−∑i,j∫vvijϕij=∑i,j∫vivjϕij+∫vDv⋅DΔϕ.-\sum_{i,j}\int vv_{ij}\phi_{ij}=\sum_{i,j}\int v_iv_j\phi_{ij}+\int vDv\cdot D\Delta\phi.

The last terms cancel. Dividing the resulting identity by two gives

12∑i,j∫vT1ij(D2v)ϕij=12∑i,j∫vivjϕij−12∫∣Dv∣2Δϕ.(95)\frac{1}{2}\sum_{i,j}\int vT_1^{ij}(D^2v)\phi_{ij}=\frac{1}{2}\sum_{i,j}\int v_iv_j\phi_{ij}-\frac{1}{2}\int|Dv|^2\Delta\phi. \tag*{(95)}

If v∈Wloc1,2v\in W^{1,2}_{\mathrm{loc}}, ordinary mollifications converge strongly in this space. Their quadratic gradient products consequently converge in Lloc1L^1_{\mathrm{loc}}, by the Cauchy–Schwarz inequality.

Thus the right-hand side is the pairing ⟨S2[v],ϕ⟩\langle S_2[v],\phi\rangle used in Remark 1.8.

We also use the Hessian measure associated with u∈Φk(Ω)u\in\Phi_k(\Omega) by Trudinger and Wang [16 Theorem 1.1]. Its classical density and weak continuity are

μk[u]=σk(D2u) dxif u∈C2(Ω),μk[um]⇀μk[u]if um→u in Lloc1(Ω).(96)\begin{aligned} \mu_k[u]&=\sigma_k(D^2u)\,dx &&\text{if }u\in C^2(\Omega),\\ \mu_k[u_m]&\rightharpoonup\mu_k[u] &&\text{if }u_m\to u\text{ in }L^1_{\mathrm{loc}}(\Omega). \tag*{(96)} \end{aligned}

Here all functions in the convergence statement belong to Φk(Ω)\Phi_k(\Omega). Weak convergence is tested against continuous compactly supported functions.

For the approximants in (91), the densities also converge to one in Lloc1L^1_{\mathrm{loc}}. Hence every admissible weak solution in Trudinger’s sense satisfies

μk[u]=dx.\mu_k[u]=dx.

Why almost-everywhere equality is insufficient

The necessary condition μk[u]=dx\mu_k[u]=dx detects singular Hessian mass that is invisible in an almost-everywhere equation. For example, set

c=(nk)−1/k,w(x)=c2∣x∣2+∣x1∣.(97)c=\binom{n}{k}^{-1/k},\qquad w(x)=\frac{c}{2}|x|^2+|x_1|. \tag*{(97)}

This function is convex, has positive quadratic lower growth, and satisfies σk(D2w)=1\sigma_k(D^2w)=1 almost everywhere. It is not a quadratic polynomial.

Indeed, on either half-space {x1>0}\{x_1>0\} or {x1<0}\{x_1<0\}, the term ∣x1∣|x_1| is affine. Thus D2w=cID^2w=cI there and σk(D2w)=(nk)ck=1\sigma_k(D^2w)=\binom{n}{k}c^k=1.

Its Hessian measure contains a nonzero singular part:

μk[w]=dx+2(n−1k−1)ck−1Hn−1↾{x1=0}.(98)\mu_k[w]=dx+2\binom{n-1}{k-1}c^{k-1}\mathcal{H}^{n-1}\mathbin{\upharpoonright}\{x_1=0\}. \tag*{(98)}

To compute it, let χε\chi_\varepsilon be a one-dimensional mollification of ∣⋅∣|\cdot| and set wε=c∣x∣2/2+χε(x1)w_\varepsilon=c|x|^2/2+\chi_\varepsilon(x_1). Then

D2wε=cI+χε′′(x1)e1⊗e1,D^2w_\varepsilon=cI+\chi_\varepsilon''(x_1)e_1\otimes e_1,
σk(D2wε)=1+(n−1k−1)ck−1χε′′(x1).\sigma_k(D^2w_\varepsilon)=1+\binom{n-1}{k-1}c^{k-1}\chi_\varepsilon''(x_1).

Since (∣t∣)′′=2δ0(|t|)''=2\delta_0 in distributions, χε′′(t) dt⇀2δ0\chi_\varepsilon''(t)\,dt\rightharpoonup2\delta_0. Integration in the remaining coordinates gives the hypersurface measure in (98).

Weak continuity identifies this limit with μk[w]\mu_k[w]. Since the hypersurface term is nonzero, μk[w]≠dx\mu_k[w]\ne dx. Thus ww is not an admissible weak solution in Trudinger’s sense.

Two algebraic lemmas

We prove the two algebraic inequalities used to construct the distance barrier. The first controls a quadratic symmetric polynomial under a reciprocal constraint.

The second transfers this inequality to the fractional transformation of the principal curvatures. Its conclusion gives the admissibility of the distance power in § 2.

Lemma A.1 (A reciprocal inequality). Let N≥1N\ge1 and η∈RN\eta\in\mathbb{R}^N satisfy

ηi>−1(1≤i≤N),∑i=1Nηi1+ηi≥0.(99)\eta_i>-1\quad(1\le i\le N),\qquad\sum_{i=1}^{N}\frac{\eta_i}{1+\eta_i}\ge0. \tag*{(99)}

Then

σ2(η)+(1−12N)σ1(η)≥0.(100)\sigma_2(\eta)+\left(1-\frac{1}{2N}\right)\sigma_1(\eta)\ge0. \tag*{(100)}

Proof. When N=1N=1, the denominator 1+η11+\eta_1 is positive, so the assumption gives η1≥0\eta_1\ge0. Since σ2(η)=0\sigma_2(\eta)=0 in one variable, the conclusion follows.

Assume N≥2N\ge2. We compare the sum of squares with the square of the sum. Set S=∑iηiS=\sum_i\eta_i and E=∑iηi2/(1+ηi)E=\sum_i\eta_i^2/(1+\eta_i). The identities

ηi1+ηi=1−11+ηi,ηi21+ηi=ηi−ηi1+ηi,\frac{\eta_i}{1+\eta_i}=1-\frac{1}{1+\eta_i},\qquad\frac{\eta_i^2}{1+\eta_i}=\eta_i-\frac{\eta_i}{1+\eta_i},

convert the hypothesis into

∑i11+ηi≤N,0≤E=S−∑iηi1+ηi≤S.(101)\sum_i\frac{1}{1+\eta_i}\le N,\qquad0\le E=S-\sum_i\frac{\eta_i}{1+\eta_i}\le S. \tag*{(101)}

In particular, S≥0S\ge0. Fix an index ii and apply Cauchy–Schwarz to the N−1N-1 positive numbers 1+ηj1+\eta_j, with j≠ij\ne i:

(N−1)2≤(∑j≠i(1+ηj))(∑j≠i11+ηj).(N-1)^2\le\left(\sum_{j\ne i}(1+\eta_j)\right)\left(\sum_{j\ne i}\frac{1}{1+\eta_j}\right).

The second factor is at most NN by (A.3). Dividing by this positive factor gives

∑j≠i(1+ηj)≥(N−1)2∑j≠i(1+ηj)−1≥(N−1)2N.(102)\sum_{j\ne i}(1+\eta_j)\ge\frac{(N-1)^2}{\sum_{j\ne i}(1+\eta_j)^{-1}}\ge\frac{(N-1)^2}{N}. \tag*{(102)}

Subtracting from ∑j(1+ηj)=N+S\sum_j(1+\eta_j)=N+S, we find $1+\eta_i\le S+2-N^{-1}foreveryfor everyi. Multiply this bound by the nonnegative number ηi2/(1+ηi)\eta_i^2/(1+\eta_i) and sum over ii. Using E≤SE \le S and S+2−N−1>0S+2-N^{-1}>0, we obtain

∑iηi2≤(S+2−1N)E≤S2+(2−1N)S.(103)\sum_i \eta_i^2 \le\left(S+2-\frac{1}{N}\right)E \le S^2+\left(2-\frac{1}{N}\right)S. \tag*{(103)}

Finally, expand S2S^2 into its diagonal and off-diagonal terms:

2σ2(η)=S2−∑iηi2≥−(2−1N)S.2\sigma_2(\eta)=S^2-\sum_i\eta_i^2\ge-\left(2-\frac{1}{N}\right)S.

Dividing by two and using S=σ1(η)S=\sigma_1(\eta) proves (A.2).

Lemma A.2 (A fractional transformation of the symmetric cone). Let 2≤k≤n2\le k\le n, a∈Γ‾k−1(n−1)a\in\overline{\Gamma}^{(n-1)}_{k-1}, and ai<1a_i<1 for every ii. Define

μi=ai1−ai,N=n−k+1,θ=2k(1−12N).(104)\mu_i=\frac{a_i}{1-a_i},\qquad N=n-k+1,\qquad\theta=\frac{2}{k}\left(1-\frac{1}{2N}\right). \tag*{(104)}

Then 0<θ<10<\theta<1, μ∈Γ‾k−1(n−1)\mu\in\overline{\Gamma}^{(n-1)}_{k-1}, and

σk(μ)+θσk−1(μ)≥0.(105)\sigma_k(\mu)+\theta\sigma_{k-1}(\mu)\ge0. \tag*{(105)}

In particular, (μ1,…,μn−1,θ)∈Γ‾k(n)(\mu_1,\ldots,\mu_{n-1},\theta)\in\overline{\Gamma}^{(n)}_k.

Proof. We connect the transformed vector μ\mu to the original vector aa by introducing

μi(s)=μi1+sμi=ai1−(1−s)ai,0≤s≤1.(106)\mu_i(s)=\frac{\mu_i}{1+s\mu_i}=\frac{a_i}{1-(1-s)a_i},\qquad0\le s\le1. \tag*{(106)}

Here μ(0)=μ\mu(0)=\mu and μ(1)=a\mu(1)=a. The denominator in the second expression is (1−s)(1−ai)+s(1-s)(1-a_i)+s, which is positive for 0≤s≤10\le s\le1. Moreover,

μi(s)−ai=(1−s)ai21−(1−s)ai≥0.(107)\mu_i(s)-a_i=\frac{(1-s)a_i^2}{1-(1-s)a_i}\ge0. \tag*{(107)}

Thus μ(s)\mu(s) is obtained from a∈Γk−1a\in\Gamma_{k-1} by adding a vector with nonnegative entries. Monotonicity of the cone gives μ(s)∈Γk−1\mu(s)\in\Gamma_{k-1}, including μ=μ(0)\mu=\mu(0).

Now introduce the polynomials

P(s)=∏i=1n−1(1+sμi),Q(s)=P(k−2)(s).(108)P(s)=\prod_{i=1}^{n-1}(1+s\mu_i),\qquad Q(s)=P^{(k-2)}(s). \tag*{(108)}

In a nonzero term of P(j)P^{(j)}, the derivatives fall on jj distinct linear factors. Each index set II with jj elements occurs j!j! times; dividing by P(s)P(s) leaves ∏i∈Iμi/(1+sμi)\prod_{i\in I}\mu_i/(1+s\mu_i). Summing over all choices of jj indices therefore gives

P(j)(s)=j!P(s)σj(μ(s)),0≤j≤k.(109)P^{(j)}(s)=j!P(s)\sigma_j(\mu(s)),\qquad0\le j\le k. \tag*{(109)}

Since 1+μi=(1−ai)−1>01+\mu_i=(1-a_i)^{-1}>0, all factors of PP are positive on [0,1][0,1]. The cone membership just proved gives

Q=(k−2)!Pσk−2(μ(s))>0,Q′=(k−1)!Pσk−1(μ(s))>0.Q=(k-2)!P\sigma_{k-2}(\mu(s))>0,\qquad Q'=(k-1)!P\sigma_{k-1}(\mu(s))>0.

For k=2k=2, the first identity uses the convention σ0=1\sigma_0=1.

The polynomial PP has only real roots. Repeated application of Rolle’s theorem, counting multiplicities, shows that QQ also has only real roots.

Since Q′>0Q' > 0 on [0,1][0,1], the polynomial QQ is nonconstant. In particular, deg⁡P≥k−1\deg P \ge k-1, so differentiation reduces its degree by exactly k−2k-2. Hence

deg⁡Q=deg⁡P−(k−2)≤n−k+1=N.(110)\deg Q = \deg P - (k-2) \le n-k+1 = N. \tag*{(110)}

Let r1,…,rdeg⁡Qr_1,\ldots,r_{\deg Q} be the roots of QQ, counted with multiplicity. Since Q>0Q > 0 on [0,1][0,1], each root lies in (−∞,0)∪(1,∞)(-\infty,0) \cup(1,\infty).

For ηj=−1/rj\eta_j=-1/r_j, a negative root gives ηj>0\eta_j>0, while a root greater than one gives −1<ηj<0-1<\eta_j<0. Factoring QQ and normalizing its value at zero yields

Q(s)Q(0)=∏j=1N(1+sηj),ηj>−1.(111)\frac{Q(s)}{Q(0)}=\prod_{j=1}^{N}(1+s\eta_j),\qquad\eta_j>-1. \tag*{(111)}

If deg⁡Q<N\deg Q<N, we add zero entries to η\eta.

Logarithmic differentiation at s=1s=1 yields

∑j=1Nηj1+ηj=Q′(1)Q(1)>0.(112)\sum_{j=1}^{N}\frac{\eta_j}{1+\eta_j}=\frac{Q'(1)}{Q(1)}>0. \tag*{(112)}

Thus η\eta satisfies both hypotheses of Lemma A.1. To express its conclusion in terms of μ\mu, differentiate (A.13) at zero:

σ1(η)=Q′(0)Q(0),σ2(η)=Q′′(0)2Q(0).\sigma_1(\eta)=\frac{Q'(0)}{Q(0)},\qquad\sigma_2(\eta)=\frac{Q''(0)}{2Q(0)}.

Since P(0)=1P(0)=1, formula (A.11) evaluates these three derivatives. We obtain

σ1(η)=(k−1)σk−1(μ)σk−2(μ),σ2(η)=k(k−1)σk(μ)2σk−2(μ).(113)\begin{aligned} \sigma_1(\eta)&=\frac{(k-1)\sigma_{k-1}(\mu)}{\sigma_{k-2}(\mu)},\\ \sigma_2(\eta)&=\frac{k(k-1)\sigma_k(\mu)}{2\sigma_{k-2}(\mu)}. \tag*{(113)} \end{aligned}

The denominator σk−2(μ)\sigma_{k-2}(\mu) is positive. Substitute these identities into (A.2) and multiply by 2σk−2(μ)/(k(k−1))2\sigma_{k-2}(\mu)/(k(k-1)) to obtain

σk(μ)+2k(1−12N)σk−1(μ)≥0.\sigma_k(\mu)+\frac{2}{k}\left(1-\frac{1}{2N}\right)\sigma_{k-1}(\mu)\ge0.

This is (A.7). When k=nk=n, the second identity in (A.15) has both sides zero, since deg⁡Q≤1\deg Q\le1 and σn(μ)=0\sigma_n(\mu)=0.

Since μ∈Γk−1\mu\in\Gamma_{k-1} and θ>0\theta>0, for 1≤j<k1\le j<k we have

σj(μ,θ)=σj(μ)+θσj−1(μ)>0.(114)\sigma_j(\mu,\theta)=\sigma_j(\mu)+\theta\sigma_{j-1}(\mu)>0. \tag*{(114)}

For every ε>0\varepsilon>0, (A.7) also gives

σk(μ,θ+ε)≥εσk−1(μ)>0.(115)\sigma_k(\mu,\theta+\varepsilon)\ge\varepsilon\sigma_{k-1}(\mu)>0. \tag*{(115)}

Thus (μ,θ+ε)∈Γk(\mu,\theta+\varepsilon)\in\Gamma_k. Letting ε↓0\varepsilon\downarrow0 proves the closure statement. The formula for θ\theta gives 0<θ<10<\theta<1. □\square

Acknowledgements

The research is assisted by ChatGPT 6.

Jingchen Hu: [email protected], [email protected].

Heming Jiao: [email protected].

References

  1. [1]J. Bao, J. Chen, B. Guan, and M. Ji, Liouville property and regularity of a Hessian quotient equation, Amer. J. Math. 125 (2003), no. 2, 301–316.DOI
  2. [2]L. Caffarelli and L. Silvestre, On the Evans–Krylov theorem, Proc. Amer. Math. Soc. 138 (2010), no. 1, 263–265.arxiv.org/abs/0905.1336
  3. [3]P. Cannarsa and C. Sinestrari, Semiconcave Functions, Hamilton–Jacobi Equations, and Optimal Control, Progress in Nonlinear Differential Equations and Their Applications, vol. 58, Birkhäuser, Boston, 2004.DOI
  4. [4]L. Chen and N. Xiang, Rigidity theorems for the entire solutions of 2-Hessian equation, J. Differential Equations 267 (2019), no. 9, 5202–5219.DOI
  5. [5]K.-S. Chou and X.-J. Wang, A variational theory of the Hessian equation, Comm. Pure Appl. Math. 54 (2001), no. 9, 1029–1064.DOI
  6. [6]J. Chu and S. Dinew, Liouville theorem for a class of Hessian equations, preprint, 2023. arXiv:2306.13825.arxiv.org/abs/2306.13825
  7. [7]J. Chu and Y. Liu, Quadratic growth rigidity for the complex Monge–Ampère equation, preprint, 2026. arXiv:2610.26016.arxiv.org/abs/2610.26016
  8. [8]J. Hu and L. Sheng, Interior C² estimates and a Liouville theorem for the complex Monge–Ampère equation, preprint, 2026. arXiv:2610.28375.arxiv.org/abs/2610.28375
  9. [9]T. Li and G. Qiu, Very weak solutions of the Dirichlet problem for 2-Hessian equation, J. Differential Equations 444 (2025), Paper No. 113577.
  10. [10]C. Mooney, The Monge–Ampère equation, Rend. Semin. Mat. Univ. Politec. Torino 76 (2018), no. 1, 93–113.DOI
  11. [11]M. Li, C. Ren, and Z. Wang, An interior estimate for convex solutions and a rigidity theorem, J. Funct. Anal. 270 (2016), no. 7, 2691–2714.arxiv.org/abs/1504.01206
  12. [12]B. Wang, A Liouville type rigidity result for entire solutions of the k-Hessian equation, preprint, 2026. arXiv:2608.29392.arxiv.org/abs/2608.29392
  13. [13]C. Wang and J. Bao, Liouville property and existence of entire solutions of Hessian equations, Nonlinear Anal. 223 (2022), Paper No. 113020.arxiv.org/abs/2203.02646
  14. [14]L. Caffarelli, L. Nirenberg, and J. Spruck, The Dirichlet problem for nonlinear second order elliptic equations, III: Functions of the eigenvalues of the Hessian, Acta Math. 155 (1985), 261–301.DOI
  15. [15]N. S. Trudinger, Weak solutions of Hessian equations, Comm. Partial Differential Equations 22 (1997), nos. 7–8, 1251–1261.DOI
  16. [16]N. S. Trudinger and X.-J. Wang, Hessian measures II, Ann. of Math. (2) 150 (1999), no. 2, 579–604.

Paper details

vlnc:W_5tc5f7k8
v1 · latestHistory (1)

published · Fresh — added to the Hub within the last 30 days

added ·

Usage

views · 0

downloads · 0

A view counts once per browser per day, a download each time the PDF is requested.

CC BY 4.0

Original source

  • Mathematics

Not recorded for this paper · View statement →

Looking for related papers…
Contents