Optimality Conditions in Control Problems with Random State Constraints in Probabilistic or Almost Sure Form II: Explicit Optimality Conditions for Probabilistic Elliptic PDE State ConstraintsResearch Paper
Motivation
Optimal control of partial differential equations under uncertainty arises whenever a physical system is steered by a control while its sources, coefficients or loads are only known in distribution. A typical requirement is a state constraint: the state (a temperature, a concentration, a displacement) must stay below a threshold everywhere in the domain. When the source term is random, demanding this for every realization is often impossible or too conservative, and the natural relaxation is a probabilistic (chance) constraint: the state must stay below the threshold uniformly on the domain with probability at least .
Probabilistic constraints are classical in operations research (Prékopa, Stochastic Programming, 1995), but in a PDE setting the constraint is an infinite family of random inequalities, one per point of the domain, and the probability function is not differentiable in general. C. Geiersbach and R. Henrion (arXiv:2306.03965v2, Math. Oper. Res. 2025) derive optimality conditions for such problems in terms of the Clarke subdifferential, building on the spherical-radial decomposition of Gaussian vectors developed by Hantoute, Henrion and Pérez-Aros (Math. Program. 2019) and van Ackooij and Henrion. This mission formalizes the PDE half of §2 of that paper: explicit formulas for the subdifferential of the probability function and the resulting optimality conditions.
Setting
Let , , be a nonempty, open, bounded set of class S: there are and with for all and . For the Dirichlet problem in , on has a unique weak solution in ; the solution operator is linear and continuous.
Fix and a threshold . The control-to-state operator is
The random vector on a probability space is Gaussian with positive definite, and is the random state. With and the probability function , the problem is
with convex and Fréchet differentiable and .
The analysis uses the spherical-radial decomposition. Let be the uniform distribution on the unit sphere , the chi distribution with degrees of freedom and density , the square root of . The radial probability function is and the radius function is the supremum of that set. For , is defined by , and
is the Clarke subdifferential, with identified with , and integrals of set-valued maps are Aumann integrals.
Formalization targets
Goal: Theorem 2.16
Assume a control with (41), and let be feasible with expected state satisfying on (42). If solves (24), there is with
and these conditions are sufficient when is a singleton for -a.e. with , or the set is bounded. The goal has no hard-coded constants; it asserts the form of the optimality system.
Milestones
In attack order: well-posedness of the state equation with the bound (Lemma 1.2); measurability of the random state (Lemma 1.4); convexity and global Lipschitz continuity of (Lemma 2.9); moderate growth (Lemma 2.10); the Lipschitz and subdifferential conclusions of the general theory under (26) (Corollary 2.11); the optimality conditions (21), (22) for (24) (Corollary 2.12); when (Lemma 2.13); the explicit formula when (Lemma 2.14); and the upper estimate (40) of with its equality cases (Theorem 2.15).
Significance
The result turns an abstract statement — lies in a multiple of the Clarke subdifferential of a probability function that has no closed form — into a condition written in the problem data: Green-type functions , the expected state, and the radius function along each direction of the sphere. This is what makes the conditions usable for checking candidate controls and for designing algorithms, and Example 2.18 of the paper computes a solution this way. The sufficiency half identifies when the nonconvex chance-constrained problem behaves like a convex one.
The paper's proofs are on paper only; to our knowledge none of these results, nor the spherical-radial decomposition, Clarke calculus in Hilbert spaces, or Aumann integrals in , has a machine-checked proof. Formalizing Theorem 2.16 requires connecting elliptic regularity, Gaussian measure theory and nonsmooth analysis, and the milestones are reusable on their own: Lemma 1.2 is a regularity statement for the Poisson equation, and Lemmas 2.13 and 2.14 are max-function subdifferential formulas in an infinite-dimensional space.
Difficulty
The obvious route — differentiate and use a Lagrange multiplier — fails because is generally not differentiable even for smooth data (Example 2.17 of the paper), and because is only Lipschitz, being a maximum over . The general subdifferential estimate (Theorem 2.2, from Hantoute–Henrion–Pérez-Aros) needs a growth condition and only gives an inclusion, since the radial probability function is not Clarke regular; the radius function is an infimum over the open set where , which is neither compact nor independent of , so standard max-function formulas do not apply directly. Directions with require a separate limiting argument through the Mordukhovich subdifferential (Lemma 2.13). The regularity input, continuity of solutions on for class-S domains (Kinderlehrer–Stampacchia 1980), is not in Mathlib.
Formalization scope
All declarations are in the namespace ChanceStateCtrl.PDE. and are EuclideanSpace ℝ (Fin d) and EuclideanSpace ℝ (Fin m), indices are Fin m, is Lp ℝ 2 (volume.restrict D). Weak derivatives are the published HunterPDE.Shared.HasWeakDeriv; a weak solution is a function continuous on , zero on , in with weak first derivatives in , satisfying the weak form against .
Standing hypotheses on every theorem: open, bounded and nonempty, , class S, and CompactSpace (closure D) (a consequence of boundedness). Nonemptiness of is the one added hypothesis: Lean's supremum over an empty domain is , not , and Lemma 2.9's identification of events fails for . is a continuous linear map to together with the hypothesis that it returns the weak solution; it is never an arbitrary operator. , as the uniform distribution on requires. is the positive semidefinite root; is its operator norm. The Clarke directional derivative is (the page's denominator is a misprint). takes values in . clco is the norm closure of the convex hull. The integral in (40)/(43) is over the restriction of to , not a renormalization. A solution is a global minimizer over ; the page's "solution to (9)" refers to this instance of the abstract problem (9).
Trivializing formalizations are ruled out: every Clarke subdifferential is taken at a point where the function is locally Lipschitz (proved in Corollary 2.11 and Theorem 2.15), the state operator is tied to the PDE, the multiplier carries its sign and complementarity, Aumann selections must be scalarly integrable, and optimality is global.
Needed infrastructure: elliptic regularity in for class-S domains; the spherical-radial decomposition of ; Clarke calculus (chain rule, max functions) on Hilbert spaces; Clarke's theorem on subdifferentials of integral functionals. Contributions of any of these as stand-alone theorems are welcome.
Selected references
- C. Geiersbach, R. Henrion, Optimality conditions in control problems with random state constraints in probabilistic or almost-sure form, arXiv:2306.03965v2, 2024; Math. Oper. Res., 2025. https://arxiv.org/abs/2306.03965
- A. Hantoute, R. Henrion, P. Pérez-Aros, Subdifferential characterization of probability functions under Gaussian distribution, Math. Program. 174 (2019), 167–194. https://doi.org/10.1007/s10107-018-1237-9
- F. H. Clarke, Optimization and Nonsmooth Analysis, Wiley, 1983; SIAM reprint 1990. https://doi.org/10.1137/1.9781611971309
- D. Kinderlehrer, G. Stampacchia, An Introduction to Variational Inequalities and Their Applications, Academic Press, 1980; SIAM reprint 2000. https://doi.org/10.1137/1.9780898719451
- A. Prékopa, Stochastic Programming, Kluwer, 1995. https://doi.org/10.1007/978-94-017-3087-7