In this note, we prove that minimizers of convex functionals with a convexity constraint and a general class of Lagrangians can be approximated by solutions to fourth-order equations of Abreu type. Our result generalizes that of Le (Twisted Harnack inequality and approximation of variational problems with a convexity constraint by singular Abreu equations. Adv. Math. 434 (2023)) where the case of quadratically growing Lagrangians was treated.
On approximation of convex functionals with a convexity constraint and general Lagrangians Paper ID : 2504.07783Title : On approximation of convex functionals with a convexity constraint and general LagrangiansAuthor : Young Ho KimClassification : math.AP (Analysis of PDEs)Publication Date : October 13, 2025 (arXiv version)Paper Link : https://arxiv.org/abs/2504.07783v2 This paper proves that minimization problems for convex functionals with convexity constraints under the class of general Lagrangian functions can be approximated by solutions of Abreu-type fourth-order equations. This result generalizes the work of Le (2023), which only addressed the case of quadratic growth Lagrangian functions.
The research addresses variational problems with convexity constraints:
inf u ∈ S [ ϕ , Ω 0 ] ∫ Ω 0 F ( x , u ( x ) , D u ( x ) ) d x \inf_{u \in S[\phi,\Omega_0]} \int_{\Omega_0} F(x, u(x), Du(x)) dx inf u ∈ S [ ϕ , Ω 0 ] ∫ Ω 0 F ( x , u ( x ) , D u ( x )) d x
where the constraint set is:
S [ ϕ , Ω 0 ] = { u : Ω → R convex , u = ϕ on Ω ∖ Ω 0 } S[\phi,\Omega_0] = \{u : \Omega \to \mathbb{R} \text{ convex}, u = \phi \text{ on } \Omega \setminus \Omega_0\} S [ ϕ , Ω 0 ] = { u : Ω → R convex , u = ϕ on Ω ∖ Ω 0 }
Economic Applications : Such variational problems arise in economics in monopolist problems (Rochet-Choné model)Numerical Computation Difficulties : The presence of convexity constraints makes it difficult to derive tractable Euler-Lagrange equationsNumerical Scheme Challenges : Such variational problems are difficult to handle numericallyDimensional Restrictions : Previous results were either limited to two-dimensional casesGrowth Condition Constraints : Required Lagrangian functions to satisfy quadratic growth conditionsNarrow Applicability : Could not handle general classes of Lagrangian functionsGeneralize the results to general Lagrangian functions in n ≥ 2 dimensional spaces by modifying the quadratic term in the approximation scheme.
Extended Approximation Theory : Generalized results from Le et al. from quadratic growth Lagrangian functions to general Lagrangian function classesExpanded Dimensional Range : Established approximation results in n ≥ 2 dimensional spacesImproved Approximation Scheme : Introduced a new penalty term G(u-φ̃ε) replacing the original quadratic termEstablished Existence Theory : Proved existence and convergence of solutions to the modified Abreu equationConsider the variational problem:
min u ∈ S [ ϕ , Ω 0 ] J ( u ) = min u ∈ S [ ϕ , Ω 0 ] ∫ Ω 0 F ( x , u ( x ) , D u ( x ) ) d x \min_{u \in S[\phi,\Omega_0]} J(u) = \min_{u \in S[\phi,\Omega_0]} \int_{\Omega_0} F(x, u(x), Du(x)) dx min u ∈ S [ ϕ , Ω 0 ] J ( u ) = min u ∈ S [ ϕ , Ω 0 ] ∫ Ω 0 F ( x , u ( x ) , D u ( x )) d x
where F satisfies general growth conditions rather than quadratic growth conditions.
The paper assumes F satisfies:
(F1) F is smooth and convex in variables z and p(F2) The derivatives of F satisfy growth estimates:
∣ ∂ F ∂ z ∣ + ∣ ∂ F ∂ p i ∣ ≤ f 0 ( ∣ z ∣ ) g 0 ( ∣ p ∣ ) \left|\frac{\partial F}{\partial z}\right| + \left|\frac{\partial F}{\partial p_i}\right| \leq f_0(|z|)g_0(|p|) ∂ z ∂ F + ∂ p i ∂ F ≤ f 0 ( ∣ z ∣ ) g 0 ( ∣ p ∣ ) Original approximation functional:
J ( v ) + 1 2 ε ∫ Ω ∖ Ω 0 ( v − ϕ ) 2 d x − ε ∫ Ω log det D 2 v d x J(v) + \frac{1}{2\varepsilon}\int_{\Omega\setminus\Omega_0}(v-\phi)^2 dx - \varepsilon\int_\Omega \log \det D^2v dx J ( v ) + 2 ε 1 ∫ Ω ∖ Ω 0 ( v − ϕ ) 2 d x − ε ∫ Ω log det D 2 v d x
Modified approximation functional:
J ε ( u ) = ∫ Ω 0 F ( x , u , D u ) d x + 1 ε ∫ Ω ∖ Ω 0 G ( u − ϕ ~ ε ) d x − ε ∫ Ω log det D 2 u d x J_\varepsilon(u) = \int_{\Omega_0} F(x,u,Du) dx + \frac{1}{\varepsilon}\int_{\Omega\setminus\Omega_0} G(u-\tilde{\phi}_\varepsilon) dx - \varepsilon\int_\Omega \log \det D^2u dx J ε ( u ) = ∫ Ω 0 F ( x , u , D u ) d x + ε 1 ∫ Ω ∖ Ω 0 G ( u − ϕ ~ ε ) d x − ε ∫ Ω log det D 2 u d x
Boundary Function Correction : ϕ ~ ε ( x ) = ϕ ( x ) + ε 1 / ( 3 n 2 ) ( e ρ ( x ) − 1 ) \tilde{\phi}_\varepsilon(x) = \phi(x) + \varepsilon^{1/(3n^2)}(e^{\rho(x)} - 1) ϕ ~ ε ( x ) = ϕ ( x ) + ε 1/ ( 3 n 2 ) ( e ρ ( x ) − 1 ) Penalty Function : G ( x ) = ∫ 0 x 2 H ( t ) d t G(x) = \int_0^{x^2} H(t) dt G ( x ) = ∫ 0 x 2 H ( t ) d t , where H ( x ) = x ( 1 + f 0 ( x ) g 0 ( x ) + f 2 ( x ) g 2 ( x ) + x f 3 ( x ) g 3 ( x ) ) H(x) = x(1 + f_0(x)g_0(x) + f_2(x)g_2(x) + xf_3(x)g_3(x)) H ( x ) = x ( 1 + f 0 ( x ) g 0 ( x ) + f 2 ( x ) g 2 ( x ) + x f 3 ( x ) g 3 ( x )) Non-quadratic Penalty Term : Replacing the quadratic term with G(u-φ̃ε) enables control of Lagrangian functions with general growthBoundary Correction Technique : Using φ̃ε to make the function "sufficiently" uniformly convex, handling non-uniformly convex Lagrangian functionsGrowth Function Matching : The construction of G precisely matches the growth properties of FTheorem 1.1 : Let Ω₀ and Ω be smooth convex domains in ℝⁿ, with Ω uniformly convex and Ω₀ ⋐ Ω. Let φ ∈ C⁵(Ω), ψ ∈ C³(Ω), φ convex, and min∂Ω ψ > 0. Let F satisfy (F1)-(F2). If 0 < ε < ε₀ < 1, then:
(i) The second boundary value problem (1.14) has a uniformly convex W⁴'ˢ(Ω) solution uε for all s ∈ (n,∞)
(ii) The solution sequence (uεₖ) converges uniformly to the minimizer of the variational problem on compact subsets of Ω
L∞ Estimates : Establish uniform boundedness of uεGradient Estimates : Utilize convexity to obtain gradient boundsHessian Estimates : Obtain upper and lower bounds for det D²uε via the ABP maximum principleRegularity Theory : Combine Hölder estimates and W⁴'ˢ estimatesDegree Theory : Use Leray-Schauder degree theory to establish existenceCompactness : Obtain subsequence convergence using uniform estimatesBoundary Conditions : Prove the limit function satisfies boundary conditionsMinimality : Demonstrate the limit function is a minimizer via variational inequalitiesChallenge : Under general growth conditions, the integral ∫ Ω 0 F p i x i ( u ε − ϕ ~ ε ) d x \int_{\Omega_0} F_{p_ix_i}(u_\varepsilon - \tilde{\phi}_\varepsilon) dx ∫ Ω 0 F p i x i ( u ε − ϕ ~ ε ) d x cannot be controlled by quadratic terms.
Solution : Construct a growth-matching penalty function G such that:
G ′ ( u ε − ϕ ~ ε ) ( u ε − ϕ ~ ε ) ≥ H ( C 8 ( 1 + ∣ u ε ∣ ) ) G'(u_\varepsilon - \tilde{\phi}_\varepsilon)(u_\varepsilon - \tilde{\phi}_\varepsilon) \geq H(C_8(1 + |u_\varepsilon|)) G ′ ( u ε − ϕ ~ ε ) ( u ε − ϕ ~ ε ) ≥ H ( C 8 ( 1 + ∣ u ε ∣ ))
Challenge : Non-uniform convexity of the Lagrangian function causes ellipticity to degenerate.
Solution : Correct the boundary function φ̃ε to introduce sufficient uniform convexity.
Challenge : The term ∂ ∂ x i ( ∂ F ∂ p i ) \frac{\partial}{\partial x_i}\left(\frac{\partial F}{\partial p_i}\right) ∂ x i ∂ ( ∂ p i ∂ F ) in equation (1.14) depends on D²uε, which is only a measure under the convexity assumption.
Solution : Use transformation techniques and twisted Harnack inequalities to handle singularities.
The paper is primarily theoretical work, verified through rigorous mathematical proofs of the following key results:
Uniform Estimates : Proved L∞, W¹'∞, W⁴'ˢ estimates for the solution sequenceConvergence : Established uniform convergence on compact subsetsMinimality : Proved the limit function is indeed a minimizer of the original variational problemA specific example mentioned in the paper is the Lagrangian function:
F ( x , z , p ) = e ∣ p ∣ 2 F(x,z,p) = e^{|p|^2} F ( x , z , p ) = e ∣ p ∣ 2
This function does not satisfy quadratic growth conditions but satisfies the general growth conditions (F1)-(F2) of the paper.
Carlier-Radice (2019) : First introduced Abreu equation approximation, but only for F = F(x,z)Le (2020-2023) : Extended to F(x,z,p) = F⁰(x,z) + F¹(x,p), but required quadratic growthLe-Zhou (2021) : General Lagrangian functions in two dimensionsThis Paper : General Lagrangian functions in n ≥ 2 dimensionsMonge-Ampère Equation Theory : Utilized boundary regularity theory from Savin, Trudinger-Wang, and othersLinearized Monge-Ampère Equation : Based on boundary regularity results of Le-SavinTwisted Harnack Inequality : Used recent results from LeSuccessfully generalized convex functional approximation theory to general Lagrangian function classes Established complete existence and convergence theory in n ≥ 2 dimensional spaces Provided effective methods for handling non-quadratic growth Lagrangian functions Smoothness Requirements : Requires high regularity φ ∈ C⁵, ψ ∈ C³Geometric Restrictions : Requires domain Ω to be uniformly convexGrowth Conditions : While generalizing quadratic growth, still requires specific growth estimatesNumerical Implementation : Convert theoretical results into practically computable numerical schemesReduced Regularity : Investigate results under lower regularity conditionsNon-convex Domains : Extend to non-convex domain casesTheoretical Breakthrough : Significantly generalizes existing theory, removing quadratic growth restrictionsTechnical Innovation : The construction of penalty function G is clever and precisely matches the growth of the Lagrangian functionRigorous Proofs : Employs advanced techniques from modern partial differential equationsClear Exposition : Mathematically precise and logically coherentLimited Application Verification : Lacks concrete numerical experiments validating theoretical resultsStrong Conditions : High requirements on function regularity and domain geometryComplex Construction : The construction of correction function φ̃ε and G is quite technicalTheoretical Contribution : Provides important generalization to variational problem approximation theoryMethodological Value : Techniques employed may apply to similar problemsApplication Prospects : Provides theoretical foundation for economics, optimization, and related fieldsTheoretical Research : Partial differential equations, calculus of variations, convex analysisNumerical Computation : Provides theoretical guidance for developing new numerical methodsApplied Mathematics : Optimal transport, monopoly theory, and other economic problemsThe paper cites 16 important references, primarily including:
Abreu's pioneering work on Kähler metrics on toric manifolds Carlier-Radice's foundational work on Abreu equation approximation Le's important series of results on singular Abreu equations Regularity theory for Monge-Ampère equations by Savin, Trudinger-Wang, and others Overall Assessment : This is a high-quality theoretical mathematics paper that achieves significant progress in variational problem approximation theory. While primarily theoretical, it provides a solid mathematical foundation for related application fields. The paper has high technical content, rigorous proofs, and represents an important contribution to the field.