--- title: Sigma-2 Hessian Equation url: https://www.emergentmind.com/topics/sigma-2-hessian-equation type: topic --- # Sigma-2 Hessian Equation Searching arXiv for recent and foundational papers on the sigma-2 Hessian equation to ground the article in the cited literature. The sigma-2 Hessian equation is the fully nonlinear elliptic partial differential equation governed by the second elementary symmetric function of the eigenvalues of the Hessian. For a twice-differentiable function \(u\), with Hessian eigenvalues \(\lambda_1,\dots,\lambda_n\), the operator is \[ \sigma_2(D^2u)=\sum_{1\le i0,\ \sigma_2(\lambda)>0\}, \] so the natural admissibility condition is \(2\)-convexity, namely \(\lambda(D^2u(x))\in\Gamma_2\) at every point. The equation occupies an intermediate position between the Poisson equation (\(k=1\)) and the Monge–Ampère equation (\(k=n\)), and its modern theory spans boundary blow-up, interior \(C^2\) estimates, viscosity regularity, global rigidity, very weak nonuniqueness, and numerical approximation [1508.04653]. ## 1. Operator, admissibility, and elliptic structure The sigma-2 operator is the second elementary symmetric function of the Hessian eigenvalues. In diagonal coordinates it is \[ \sigma_2(D^2u)=\sum_{i0,\ \sigma_2(\lambda)>0\}, \] equivalently the class of \(2\)-convex functions. On this cone, the linearization is positive definite. For the constant-right-hand-side equation, the linearized operator is \[ F^{ij}(D^2u)=\frac{\partial\sigma_2}{\partial u_{ij}} =\sigma_1(D^2u)\delta_{ij}-u_{ij}, \] which is elliptic precisely on \(\Gamma_2\) [1906.10588]. This admissibility restriction is not a peripheral technicality; it is the mechanism that distinguishes the equation from uniformly elliptic models. The operator is elliptic only on a restricted set of functions, and the PDE theory therefore intertwines second-order estimates with convexity-type constraints, viscosity admissibility, and geometric properties of level sets [1502.04969]. In dimension three, one equivalent description of \(\Gamma_2\) is that the Laplacian restricted to every coordinate plane is nonnegative together with \(\sigma_2>0\), a characterization that is useful both analytically and numerically [1502.04969]. A recurrent theme in the literature is the distinction between \(2\)-convexity and stricter forms of positivity. One line of work studies classical or viscosity solutions on the admissible branch; another studies semiconvex or dynamically semiconvex regimes, where lower bounds on the minimal eigenvalue relative to \(\Delta u\) or \(\sigma_2\) restore enough effective ellipticity for higher-dimensional interior estimates [2305.12587]. ## 2. Boundary value problems and Keller–Osserman blow-up theory A fundamental nonlinear boundary phenomenon is the boundary blow-up problem \[ \sigma_2(D^2u)=g(u)\quad\text{in }\Omega,\qquad \lim_{x\to\partial\Omega}u(x)=+\infty, \] posed on a smooth bounded domain \(\Omega\). For the sigma-2 case, a necessary and sufficient existence criterion is the Keller–Osserman integral condition \[ \int_{\beta}^{\infty}\Bigl(\int_0^t g(s)\,ds\Bigr)^{-1/3}\,dt<\infty \] under the hypotheses that \(g:\mathbb R\to[0,\infty)\) is \(C^{2,\alpha}\) on \(\{g>0\}\), convex, nondecreasing, \(g(s)=0\) for \(s\le0\), \(g(s)>0\) for \(s>0\), and \(\partial\Omega\) is \(C^{4+\alpha}\) and strictly \(1\)-convex [1508.04653]. The proof uses a subsolution–supersolution framework. A strict subsolution \(\phi\) is obtained from a Dirichlet problem with constant boundary value \(c\) and right-hand side \(g(c)+1\), exploiting monotonicity of \(g\). A supersolution \(\psi\) is built from the large-boundary-data linear problem \[ \Delta \psi = n\,g(\psi),\qquad \psi=+\infty\ \text{on }\partial\Omega, \] together with Maclaurin’s inequality \(\sigma_2^{1/2}\le (\Delta u)/n\). Solving finite Dirichlet problems \[ \sigma_2(D^2u_N)=g(u_N),\qquad u_N=N\ \text{on }\partial\Omega, \] and passing to the monotone limit yields a blow-up solution between the two barriers [1508.04653]. Two model nonlinearities exhibit the sharpness of the criterion. For \(g(u)=u^p\), the Keller–Osserman integral converges exactly when \(p>2\), so existence of a positive boundary-blow-up solution is equivalent to \(p>2\). For \(g(u)=e^u\), the integral always converges, hence the equation \[ \sigma_2(D^2u)=e^u \] admits at least one positive blow-up solution [1508.04653]. This theory shows that the sigma-2 equation shares the classical Keller–Osserman dichotomy but with the exponent \(-1/3\) in the integral condition, reflecting the algebraic structure of the second Hessian operator rather than the Laplacian or Monge–Ampère cases. A plausible implication is that boundary singularity formation for \(\sigma_2\) is already rigidly encoded at the level of one-dimensional integral growth conditions. ## 3. Interior Hessian estimates and regularity theory Interior \(C^2\) estimates are central because once \(|D^2u|\) is controlled, the equation becomes uniformly elliptic and concave in the admissible cone, so Evans–Krylov and Schauder theory yield higher regularity. In dimension three, a priori interior \(C^2\) bounds for \[ \sigma_2(D^2u)=f(x,u,Du)>0 \] were established using a combination of differentiated identities, a mean-value inequality on a graph submanifold, and a maximum-principle argument for a refined scalar test function [1711.00948]. For semiconvex solutions of the constant equation \[ \sigma_2(D^2u)=1, \] there is an interior Hessian estimate of exponential type: \[ \sup_{x\in B_{R/2}(0)}|D^2u(x)| \le C(n,K)\exp\bigl(C(n,K)\,\|u\|_{L^\infty(B_R)}/R^2\bigr), \] assuming \(u\in C^4(B_R)\) and \(D^2u\ge -K I\) [1911.04246]. An equivalent version uses \(\|Du\|_{L^\infty}/R\). The proof rests on three ingredients: a Jacobi-type differential inequality for \(b=\ln\lambda_1\), the Legendre–Löwy transform \[ y=Du(x)+Kx,\qquad w(y)=x\cdot y-(u(x)+K|x|^2/2), \] which converts the nonuniformly elliptic linearization into a uniformly elliptic equation for \(w\), and a mean-value inequality for the transformed operator [1911.04246]. The novelty is that semiconvexity \(D^2u\ge -KI\) with arbitrary \(K\) suffices, improving earlier “almost convex” hypotheses. Dimension four marks a threshold in the known a priori theory. For admissible \(u\in C^4(B_r(x_0))\) solving \[ \sigma_2(D^2u)=f(x,u,Du) \quad\text{in }B_r(x_0)\subset\mathbb R^4, \] with \(f>0\) and \(f\in C^{1,1}\), there is an interior estimate \[ \|D^2u(x_0)\|\le C\Bigl(\|u\|_{C^1(B_r)},\ \|f\|_{C^{1,1}},\ \|1/f\|_{L^\infty},\ r\Bigr) \] [2509.03217]. The proof combines an almost-Jacobi inequality for \(b=\log\Delta u\), a Pogorelov-type doubling argument based on a test function of the form \[ P(x)=2\log(9-|x|^2)+\alpha(x\cdot Du-u)+\tfrac\beta2|Du|^2+\log(b-\inf b), \] and a small-scale control obtained from almost-everywhere twice differentiability of Lipschitz viscosity solutions together with a generalized Savin small-perturbation theorem [2509.03217]. A closely related four-dimensional estimate for \(f\equiv1\) and then \(f(x)\) was established earlier by a doubling-plus-blow-up method [2305.12587]. In dimensions \(n\ge5\), the general interior regularity question remains open, and Pogorelov-type examples obstruct dimension-free estimates for all admissible solutions. Existing higher-dimensional \(C^2\) bounds require extra dynamic semi-convexity conditions, such as a lower bound on \(\lambda_{\min}(D^2u)\) in terms of \(\Delta u\) or \(\sigma_2(D^2u)\) [2305.12587]. This suggests that the borderline between regularity and singularity in higher dimensions is not merely technical but structural. Recent work adds an a priori interior \(C^2\) estimate in terms of the \(W^{2,p}\) norm, for any \(p>2\): \[ \|D^2u\|_{L^\infty(B_{1/2})} \le C\bigl(n,p,\|D^2u\|_{L^p(B_1)}\bigr) \] for smooth solutions of \(\sigma_2(D^2u)=1\) in \(B_1\) [2505.14586]. The proof rewrites the linearized equation in divergence form as a Laplace–Beltrami equation, uses Krylov–Safonov and ABP on the sublevel sets of \(\Delta u\), and then exploits the algebraic relation \(\sigma_2(D^2u)=1\) to bound all eigenvalues once \(\Delta u\) is controlled [2505.14586]. ## 4. Convexity, strict \(2\)-convexity, and geometric interpretation The sigma-2 equation has a strong geometric side. For convex viscosity solutions of \[ \sigma_2(D^2u)\ge1, \] strict \(2\)-convexity means that every supporting linear function has a contact set of dimension at most \(n-2\). This property was proved by constructing paraboloid barriers in thin cylinders and showing that any \((n-1)\)-dimensional contact set would contradict the viscosity supersolution condition [2006.05940]. The same work derives short proofs of smoothness and interior \(C^2\) estimates for convex viscosity solutions of \[ \sigma_2(D^2u)=1, \] yielding \[ |D^2u(0)|\le C\bigl(n,\|u\|_{L^\infty(B_1)}\bigr) \] and then full interior \(C^\infty\) regularity by Evans–Krylov and Schauder bootstrap [2006.05940]. The argument passes from qualitative strict \(2\)-convexity to a quantitative control on the geometry of sublevel sets and then invokes classical Dirichlet solvability and Pogorelov-type interior estimates. A more geometric strengthening states that viscosity solutions of \[ \sigma_2(D^2u)=1\quad\text{in }B_1 \] cannot touch a harmonic function on a smooth embedded piece of a minimal hypersurface from below [2505.14586]. The barrier used there has the form \[ v_K(x)=K\,d(x)+h(\pi(x))+\varepsilon f(d(x)), \] where \(d(x)\) is signed distance to the minimal surface, \(\pi(x)\) is nearest-point projection, and \(f(s)=s|\log s|\). The non-contact theorem can be viewed as a form of strict \(2\)-convexity and, in the convex setting, forces the contact set with any supporting hyperplane to have dimension at most \(n-2\) [2505.14586]. The equation also connects to special Lagrangian geometry. In dimension three, the critical-phase special Lagrangian equation \[ \sum_{i=1}^3 \arctan\lambda_i=\frac{\pi}{2} \] is equivalent to \(\sigma_2(\lambda)=1\), so rigidity for sigma-2 solutions implies a Bernstein-type theorem for global special Lagrangian graphs in \(\mathbb R^3\times\mathbb R^3\) under quadratic growth [1809.02902]. More generally, the reformulation of \(\sigma_2(D^2u)=f(x)^2\) as a twisted special Lagrangian equation in a weighted Euclidean metric yields interior \(C^{2,\alpha}\) regularity for continuous viscosity solutions in dimension three when \(f\) is positive and Lipschitz [2311.14260]. These results show that “strict \(2\)-convexity” is not merely an eigenvalue condition. It manifests through contact geometry, sublevel-set thickness, non-contact with minimal-surface harmonic data, and the geometry of gradient graphs. ## 5. Entire solutions, Liouville theorems, and rigidity A major global theme is rigidity of entire admissible solutions. If \(u\in C^4(\mathbb R^n)\), \(n\ge3\), is an entire \(2\)-convex solution of \[ \sigma_2(D^2u)=1 \] with quadratic growth from below, \[ u(x)\ge c|x|^2-b\quad\text{for }|x|\gg1, \] then \(u\) must be a quadratic polynomial [1906.10588]. The proof rescales \(u\) on sublevel sets \(\Omega_R=\{y:u(Ry)\le R^2\}\), uses a global gradient estimate on mean-convex domains together with a Pogorelov estimate, and obtains a uniform bound on \(\Delta u\). Evans–Krylov theory then implies higher regularity, and the limit argument forces \(D^2u\) to be constant [1906.10588]. A related rigidity theorem holds in all dimensions for entire \(2\)-convex solutions of \(\sigma_2(D^2u)=1\) that satisfy quadratic growth and a one-sided bound \[ \sigma_3(D^2u)\ge -A. \] Under these assumptions, the solution is again a quadratic polynomial; in dimension three the \(\sigma_3\)-assumption is redundant [1809.02902]. The proof proceeds through a local Pogorelov-type estimate in bounded domains, \[ \max_{\Omega}(-u)^\alpha |D^2u|\le C, \] followed by scaling and then Evans–Krylov and Schauder theory [1809.02902]. The literature frames these results as an extension of Jörgens–Calabi–Pogorelov-type rigidity from the Monge–Ampère equation to the intermediate Hessian setting [1906.10588]. A plausible implication is that, for admissible entire sigma-2 solutions, the main obstruction to nonquadratic behavior is not lack of global convexity alone but the combination of growth, admissibility, and interior \(C^2\) control. Higher-dimensional regularity results also feed into Liouville theory. In four dimensions, interior Hessian estimates imply that \(C^0\)-viscosity solutions are \(C^\infty\) in the interior, and entire solutions with quadratic growth under the dynamic semi-convexity condition are necessarily quadratic polynomials [2305.12587]. ## 6. Weak solution theories, nonuniqueness, and computational methods The sigma-2 equation supports several inequivalent weak frameworks. The viscosity theory is stable under uniform limits and coincides with classical theory on smooth \(2\)-convex solutions [2505.14586]. By contrast, “very weak solutions” may be defined through a double-divergence identity: \[ -\sum_{i,j=1}^n \int_\Omega \sigma_2^{ij}(D^2\phi)\,\partial_i v\,\partial_j v = \int_\Omega f\,\phi \] for all \(\phi\in C_c^\infty(\Omega)\), where \(v\in W^{1,2}_{\mathrm{loc}}(\Omega)\) [2403.09356]. Using convex integration and a cut-off technique, it was shown that for \(\alpha<1/(1+n+n^2)\) there exist infinitely many \(C^{1,\alpha}\) very weak solutions of the Dirichlet problem with prescribed boundary value, under \(f>0\), \(\Omega\in C^{2,\kappa}\), and \(f\in C^{0,\kappa}\) [2403.09356]. This establishes low-regularity flexibility and nonuniqueness in sharp contrast with the classical viscosity theory, where smooth or admissible solutions are unique under the usual structural assumptions [2403.09356]. The contrast is substantive rather than terminological: the convex integration solutions are very weak solutions, not viscosity solutions. This addresses a common misconception that all reasonable weak formulations of \(\sigma_2(D^2u)=f\) should share the same uniqueness and regularity properties. The computational literature reflects the same admissibility issues. A wide-stencil monotone finite-difference scheme for the three-dimensional Dirichlet problem is provably convergent to the viscosity solution via the Barles–Souganidis framework, while a standard centered-difference scheme is more accurate but not monotone and thus lacks an unconditional convergence proof [1502.04969]. Both are typically solved with Newton’s method. For smooth solutions, the accurate scheme exhibits second-order behavior in experiments, whereas the monotone scheme is usually first order unless very wide stencils are used [1502.04969]. A nine-point finite-difference discretization for \(k\)-Hessian equations yields a locally unique discrete solution with \(O(h^2)\) error for smooth nondegenerate solutions, and Newton’s method converges quadratically once iterates enter the local contraction ball [1406.5366]. That work also introduced parameter-free subharmonicity-preserving and nonlinear Gauss–Seidel-type iterations for the \(2\)-Hessian case, designed to preserve discrete \(2\)-convexity and to handle non-smooth solutions in practice [1406.5366]. ## 7. Dimension dependence, extensions, and open directions The theory of the sigma-2 Hessian equation is sharply dimension-dependent. In dimensions \(2,3,4\), the available interior Hessian estimates are strong enough to force smoothness of viscosity solutions in broad regimes; in higher dimensions the general interior regularity problem remains open [2505.14586]. Dimension four is a particularly important threshold: the constant-right-hand-side equation and positive \(C^{1,1}\) inhomogeneous equations admit interior Hessian estimates without extra structural assumptions there, while in \(n\ge5\) one must impose dynamic semi-convexity or related lower bounds on the smallest eigenvalue [2509.03217]. Several extensions move beyond the pure equation \(\sigma_2(D^2u)=1\). In dimension three, there are interior \(C^{2,\alpha}\) estimates for \(C^0\)-viscosity solutions of \[ \sigma_2(D^2u)=f(x)^2 \] with positive Lipschitz \(f\), obtained through the twisted special Lagrangian formulation [2311.14260]. More recent work establishes interior Hessian estimates for \[ \sigma_2(D^2u)=\psi(x,u,\nabla u) \] in dimension three for \(2\)-convex solutions with \(C^1\) bounds, and in higher dimensions under semiconvexity assumptions [2602.14064]. In dimension four, positive \(C^{1,1}\) right-hand sides \(f(x,u,Du)\) can also be treated [2509.03217]. Another active direction concerns Pogorelov-type \(C^2\) estimates. For admissible solutions of the Dirichlet problem \[ \sigma_2(D^2u)=f(x,u,Du),\qquad u=0\ \text{on }\partial\Omega, \] with \(00\), one has an interior estimate \[ \sup_{\Omega'}|D^2u| \le C\bigl(\|u\|_{C^1(\overline\Omega)},\|f\|_{C^2},c_0,\operatorname{dist}(\Omega',\partial\Omega)\bigr) \] and a Pogorelov bound \[ (-u)^\beta \Delta u\le C \] [2504.06711]. The dependence on the lower bound \(c_0\) reflects the intermediate-Hessian nature of the equation; the paper notes that this hypothesis can be removed when \(k=n\), but not for \(k=2\) [2504.06711]. Across these developments, a consistent pattern emerges. The sigma-2 Hessian equation is governed by a tension between ellipticity restricted to \(\Gamma_2\), geometric strengthening such as strict \(2\)-convexity, and the search for dimension-sensitive mechanisms that restore effective uniform ellipticity. This suggests that future progress will likely continue to depend on hybrid tools: Jacobi inequalities, convexity transforms, doubling arguments, small-perturbation theory, and geometric barrier constructions rather than a single universal regularity principle. Source: https://www.emergentmind.com/topics/sigma-2-hessian-equation