A Hessian-dependent functional with free boundaries and applications to mean-field games
Abstract
We study a Hessian-dependent functional driven by a fully nonlinear operator. The associated Euler-Lagrange equation is a fully nonlinear mean-field game with free boundaries. Our findings include the existence of solutions to the mean-field game, together with Hölder continuity of the value function and improved integrability of the density. In addition, we prove the reduced free boundary is a set of finite perimeter. To conclude our analysis, we prove a -convergence result for the functional.
Keywords: Hessian-dependent functionals, fully nonlinear mean-field games with free boundaries, regularity theory.
MSC2020: 35B65; 35J35; 35R35; 35A01; 35Q89.
1 Introduction
We examine Hessian-dependent functionals of the form
(1) |
where is a uniformly elliptic operator, is a fixed constant, , and stands for the space of symmetric matrices of order . Our results include the existence of minimizers for (1), amounting to the existence of solutions to a fully nonlinear mean-field game with free boundaries. We prove Hölder-continuity of minimizers and improved integrability of the density. We also prove the free boundary has finite perimeter. Finally, we establish a result on the -convergence of and examine its consequences.
The functional in (1) is inspired by the usual one-phase Bernoulli problem, driven by the Dirichlet energy. To a limited extent, we understand as a Hessian-dependent counterpart of that problem. See [2]; see also [11].
The analysis of (1) relates closely with the system
(2) |
where denotes the derivative of with respect to the entry of . Here, the unknown is a pair solving the problem in a sense we make precise further.
The system in (2) amounts to the Euler-Lagrange equation associated with (1). We notice that (2) satisfies an adjoint structure. Its double-divergence equation is the formal adjoint, in the -sense, of the linearized fully nonlinear problem. Due to such a distinctive pattern, we refer to (2) as a fully nonlinear mean-field game with free boundary. Interpreting the first equation in (2) as a Hamilton-Jacobi equation, it becomes natural to ask about the underlying stochastic optimal control problem. We do not examine this matter in the present paper.
Fully nonlinear elliptic operators and equations in the double-divergence form have been studied by many authors. An attempt to put together a comprehensive list of references on those topics is unrealistic. For that reason, we mention solely the monographs [8, 12] and the references therein.
Our analysis sits at the intersection of Hessian-dependent functionals, free boundary problems, and mean-field games systems. Hence we proceed with some context on those classes of problems. Hessian-dependent functionals play an essencial role in various contexts. From a purely mathematical viewpoint, they are useful to produce examples of conformally invariant energies. In dimension , this is the case of
whose Euler-Lagrange equation is the biharmonic operator; see [15, 14].
When it comes to applications, we mention the realm of mechanics of solids. In particular, the analysis of energy-driven pattern formation and nonlinear elasticity. For example, Hessian-dependent models play a role in studying the occurrence of wrinkles in a twisted ribbon [24]. The energy modeling the system depends on the thickness of the ribbon, denoted with , and two symmetric tensors and . It has the form
Although depends on its arguments only through lower-order terms, the tensor depends on . Another instance where Hessian-dependent functionals appear is the analysis of blister patterns in thin films on compliant substrates [6]. Here the phenomena are modeled in terms of lower-order quantities, a small Hessian-dependent perturbation, and a parameter . An important question concerning this class of problems is the limiting behavior ; in fact, one expects that the lower and upper bounds of the functional scale similarly. We refer the reader to [18, 19, 29]. In this context, (1) amounts to an energy penalized by the measure of the positive phase. See also [23, 5, 4].
As noticed before, one can state the Euler-Lagrange equation associated with (1) in terms of the fully nonlinear mean-field game system with free boundaries (2). Mean-field games comprise a set of methods and techniques to model strategic interactions involving many players [25, 26, 27]; see also [28]. At the intersection of partial differential equations (PDE), stochastic analysis and numerical methods, this class of problems has attracted the attention of several authors, who have developed the theory in various directions.
Under assumptions on the stochastic dynamics governing the problem (e.g., independence of the Brownian motions among the population of players), it is possible to write a mean-field game in terms of a coupling. It comprises a Hamilton-Jacobi equation, accounting for the value of the game, and a Fokker-Planck equation describing the evolution of the population. For further references on the topics, we mention the monographs [13, 7, 22, 1]. Fully nonlinear mean-field games are the subject of [3, 16].
The interesting aspect in (2) concerns the appearance of a free boundary. At least heuristically, the game is played only in the regions where the value function is strictly positive. Combined with the free boundary condition, (2) models a game in which players optimize in the region where the value function is positive and might face extinction according to a flux condition endogenously determined.
Our first contribution is to prove the existence of solutions for the mean-field game system in (2). To that end, we impose natural assumptions on the data of the problem. Namely, we require the operator driving the problem to be uniformly elliptic (Assumption A1) and to satisfy a norm-like growth condition (Assumption A3). We also require the operator to be convex (Assumption (A2) and the boundary data to be a function in (Assumption A4). We report our findings in the following theorem.
Theorem 1 (Existence and regularity of solutions).
If, in addition, is strictly convex and , we can prove that is not only integrable but is indeed an -function, with estimates; c.f. [21]. To establish Theorem 1, we start with the direct method in the calculus of variations and the existence of minimizers for (1). Then we turn our attention to (2). First, we resort to the theory of weak solutions available for equations in the double-divergence form. Finally, elements in the -viscosity theory lead to the existence of solutions to the system. To complete the proof, we resort to a delicate application of Sobolev inequalities.
Once we have established the existence of solutions for (2) and produced a regularity result, we examine the free boundary. Regularity results for the solutions build upon ingredients of geometric measure theory to ensure the reduced free boundary is a set of finite perimeter. We summarize our findings in this direction in the following result.
Theorem 2 (Free boundary condition and finite perimeter).
The remainder of this paper is organized as follows. Section 2.1 details our main assumptions, whereas Section 2.2 gathers preliminary material and results. Section 3 presents the proof of Theorem 1. In Section 4, we examine the free boundary and put forward the proof of Theorem 2. A final section closes the paper with a -convergence result and some consequences.
2 Preliminaries
This section presents the main assumptions under which we work and collects some preliminary notions and results.
2.1 Main assumptions
We proceed with a condition on the uniform ellipticity of the operator .
A 1 (Uniform ellipticity).
We suppose the operator is -elliptic for some . That is,
(3) |
for every , with . We also suppose . Finally, we require , for every , where stands for the derivative of with respect to the entry of .
Remark 1.
Next, we impose a convexity condition on the operator .
A 2 (Convexity of the operator ).
We suppose the operator to be convex with respect to .
Part of our arguments requires to satisfy a coercivity condition in the entire . To that end, we strength A1 as follows.
A 3 (Growth condition).
We suppose there exists such that the operator satisfies
for every .
The typical example of an operator , satisfying the former assumptions, depends on through its norm. For instance, let be a constant and consider
For a more general operator, including explicit dependence on the space-variable , we consider and define as
For further examples related to the previous one see [23, 6].
Notice that A3 implies ; this fact is important when studying the weak lower semicontinuity of our functional. We conclude this section with an assumption on the boundary data.
A 4 (Boundary data).
We suppose the function is non-negative and non-trivial.
2.2 Preliminary notions and results
For the sake of completeness, we recall definitions and former results we use throughout the manuscript. We continue with a definition
Definition 1 (Affine Sobolev spaces).
Let and . We say that
if and .
From a PDE perspective, having is tantamount to prescribe on in the Sobolev sense. This interpretation will be helpful when relating (1) and (2).
As usual in the literature on mean-field games [25, 26, 27], a solution to (2) relies on two distinct definitions – namely, the notions of viscosity and weak (distributional) solutions. We proceed by recalling the definition of -viscosity solution of a fully nonlinear elliptic equation; see [10, Definition 2.1].
Definition 2 (-viscosity solutions).
Let be a fully nonlinear operator satisfying A1 and , for . A function is an -viscosity sub-solution (resp. super-solution) of
if, for all , whenever , is open, and
then cannot have a local maximum (resp. minimum) in . Moreover, if is both an -viscosity sub-solution and an -viscosity super-solution, is said to be an -viscosity solution.
The definition of -viscosity solution is necessary since -functions might not be defined at the points where the usual conditions must be tested. For a comprehensive account of this notion, we refer the reader to [10]. We continue with the definition of weak solutions for double-divergence equations.
Definition 3 (Weak solution).
Let and denote . Suppose
We say is a weak solution to
if, for every we have
Definition 4 (Solution for the MFG system).
The pair is a weak solution to (2) if the following hold:
-
1.
We have and , with ;
-
2.
The function is an -viscosity solution to
-
3.
The function is a weak solution to
Next, we recall the Poincaré’s inequality for functions lacking compact support. In particular, we are interested in .
Lemma 1 (Poincaré’s inequality).
Let and be the Poincaré’s constant associated with and the dimension . Then for every , there exists and such that
For the detailed proof of this fact, we refer the reader to [20, Lemma 2.7, p. 22]. It follows from and the usual Poincaré’s inequality.
In Section 3, we deal with the existence of minimizers for in . Our reasoning uses the weak lower-semicontinuity of the functional
this is the content of the following lemma.
3 Existence of solutions
In this section, we present the proof of Theorem 1; we start by establishing the existence of minimizers for (1).
Proposition 1 (Existence of minimizers).
Proof.
Under Assumptions A2 and A3, the existence of minimizers follows from the direct method in the calculus of variations. We split the argument into three steps.
Step 1 - We first examine
In view of the Remark 1, . Furthermore, since ,
Hence, . Let be a minimizing sequence; there exists such that
for every . Therefore, for all ,
In the next step the upper bound for builds upon properties of the functional.
Step 2 - As a consequence of the former inequality, we infer that is uniformly bounded in . Since , the embedding is compact. Furthermore, we conclude that is uniformly bounded in ; it follows from Lemma 1 combined with general facts [17]. Hence, there exists such that
(4) |
and
(5) |
The result follows at once if we ensure that
(6) |
and
(7) |
hold. Notice that Lemma 2 combines the convergence mode in (4) to yield (6). In the sequel, we establish (7).
Step 3 - Because of the strong convergence (5), there exists a subsequence, also denoted with , and a negligible subset , such that for every in . As a consequence, if , there exists such that for every . If , then . Therefore,
(8) |
for almost every . Hence,
which completes the proof. ∎
Remark 2.
We notice the minimizing sequence is uniformly bounded in . As a consequence, it is also uniformly bounded in some Hölder space. Therefore, we could have used uniform convergence in (4).
We close this section with the proof of Theorem 1.
Proof of Theorem 1.
We split the proof into four steps.
Step 1 - Let be the minimizer for (1) whose existence follows from Proposition 1. There exists such that is well-defined for every , with . This fact, combined with A3, implies that for almost every . Therefore, satisfies in the -viscosity sense; see [10, Lemma 2.6].
Step 2 - By considering a variation of compactly supported in , we obtain
(9) |
for every . Set ; we infer that is well-defined and satisfies for almost every . In addition,
that is, . Finally, we notice the integral in (9) is well-defined and leads to
for every .
Step 3 - It remains to check that is an -viscosity solution to the first equation in (2). The definition of implies that satisfies
for almost every . As before, an application of [10, Lemma 2.6] ends the proof.
Step 4 - We prove that for every . We start by recalling the Gagliardo-Nirenberg inequality for bounded domains. Being a minimizer for (1), there exists such that
(10) |
provided
(11) |
for some and . We notice the -norm of appears in (11) because
because of Assumption A3 and the fact that is a competitor for .
Remark 3 (Improved integrability for ).
Let be a weak solution to the fully nonlinear MFG (2). In case is strictly convex and , we claim that . In fact, is defined almost everywhere in as . Under the strict convexity of and , solutions to the Euler-Lagrange equation are minimizers for the functional (1). Hence, A3 transmits the integrability of to , and the claim follows. Compare with [21]; see also [8]. Re-writing the exponent above as we quantify the improved integrability of in face of the -regime.
Remark 4 (Improved regularity for the value function).
The value function is -Hölder-continuous, for every . Hence, the regularity established in the former argument amounts to an improvement of the usual Krylov-Safonov regularity theory implied by uniform ellipticity.
4 Information on the free boundary
In the sequel, we examine local properties of the free boundary and present the proof of Theorem 2. The following corollary connects the regularity of minimizers with information on the free boundary. We refer to it when proving the first part of Theorem 2.
Corollary 1.
Proof.
We split the argument into four steps and begin by proving that, for given (fixed and to be chosen later), one gets
(14) |
for some universal constant
Step 1 - We begin by fixing a function such that
and
for some universal constant . For a fixed , consider the functions
and
which, by the minimality condition (12), give
(15) |
Step 2 - Now we calculate in ; to do so, first notice that can be writen as
Therefore,
where . Combining the former equality with Assumption A3 one gets
(16) |
In the next step we detail an involved chain of inequalities used in the argument.
Step 3 - We combine (16) and (15), set
and resort to Assumption A3 to compute
Therefore
(17) |
where
and
The argument in the proof of Theorem 1 ensures that
Therefore, there exists a universal constant such that
(18) |
Recall that . Because of the bounds imposed on and the estimates available for , we conclude
Hence, by requiring , the Hölder inequality and the Theorem 1 yield
where and are positive, universal constants. Hence, there exists a universal constant such that
Thus one finds a universal constant for which
(19) |
By combining (17), (18) and (19) we have that
(20) |
Finally, set . Theorem 1 ensures that for every one obtains
where the constant is now universal.
Step 4 - Now, Assumption A3 yields
To estimate we recall the Galiardo-Nirenberg inequality for bounded domains. If , there exists universal constants such that
(21) |
provided
for some and , where . Given and it is always possible to find and satisfying (21), which implies that there exists such that
Also,
By combining the former inequalities, we get
where is an universal constant. A straightforward application of the area formula yields
and finishes the proof. ∎
4.1 Proof of Theorem 2
In what follows, we organize the previous results and present the proof of Theorem 2. The Sobolev regularity of minimizers and its corollary leads to the finite perimeter of the reduced free boundary.
5 Perturbation analysis via -convergence
This section specializes the operator to be the norm and considers small values of the parameter in (1). We regard the functional
(22) |
as a free boundary perturbation of
(23) |
Denote with a minimizer for (22) and with the minimizer for (23). We are interested in the behavior of , as . In particular, we search for the topologies where the convergence is available. Our starting point is a -convergence result. Namely, we first prove that as .
Although interesting on its own merits, the -convergence problem is motivated by its potential consequences on the regularity theory of minimizers to (22). Indeed, we use properties of -convergence to prove an approximation result. It states that minimizers are close, in a suitable topology, to a minimizer of the -limit (see Proposition 2). This type of approximation result is central to perturbative methods in regularity theory; see, for instance, [9, 12]. We believe the -convergence analysis can be used as an ingredient in the study of improved regularity for Hessian-dependent functionals through approximation methods. We proceed with some auxiliary lemmas.
Lemma 3 (Equicoerciveness).
Let be fixed and be a sequence such that , as . Define the functional as
if , and in case . Let be such that
(24) |
for every and some . Then , uniformly in , for some .
Proof.
Before continuing, we introduce the functional , given by
if , and if . The next lemma relates and .
Lemma 4.
Let be fixed and be a sequence of real numbers so that , as . For each there exists a sequence converging strongly to in , such that
(25) |
Proof.
Let be given and , for every . If , we get
and (25) is immediately satisfied. Conversely, suppose . In that case, we have
∎
Lemma 5.
Let be fixed and be a sequence of real numbers so that , as . Given and , with strongly in , we have
(26) |
Proof.
To deduce (26) from the strong convergence, suppose first . Then
and (26) follows. Otherwise, suppose .
Through a subsequence, if necessary, we can suppose the in (26) is in fact a limit. If such a limit is not finite, then (26) trivially holds. Suppose otherwise; if this limit is finite, there exists such that
for every , large enough (and therefore for every ). As a consequence, is uniformly bounded; evoking once again standard inequalities for Sobolev functions, one infers the existence of a constant such that
The weakly lower semi-continuity of the -norm yields
and completes the proof. ∎
Theorem 3 (Gamma Convergence).
Let be fixed and be a sequence of real numbers so that , as . Then .
In the sequel, we explore a consequence of the -convergence result. It consists of an approximation result by -regular functions.
5.1 Regular approximations
We have proved that minimizers for (1) are Hölder-continuous. However, the use of -convergence allows us to arbitrarily approximate minimizers by -regular functions. This is the content of the following proposition
Proposition 2 (-approximation).
Let be fixed. Given , there exists such that, if and be a minimizer for (1), one can find satisfying
Proof.
We use a contradiction argument. Suppose the statement of the proposition is false. In this case, there exist a real number and sequences and such that
as ,
for every and every , but
(27) |
for every , and every .
However,
for some . Hence, there exists such that converges , weakly in and strongly in . That is tantamount to say that is an accumulation point for the sequence .
Acknowledgements The authors are grateful to Giovanni Bellettini for his comments on the material in this paper. This work was partially supported by the Centre for Mathematics of the University of Coimbra - UIDB/00324/2020, funded by the Portuguese Government through FCT/MCTES. JC is funded by FAPERJ-Brazil (# E26/202.075/2020). EP is partly funded by FAPERJ-Brazil (Grant # E26/200.002/2018), ICTP-Trieste and Instituto Serrapilheira (Grant # 1811-25904). This study was financed in part by the Coordenação de Aperfeiçoamento de Pessoal de Nível Superior - Brazil (CAPES) - Finance Code 001.
References
- [1] Y. Achdou, P. Cardaliaguet, F. Delarue, A. Porretta, and F. Santambrogio. Mean field games, volume 2281 of Lecture Notes in Mathematics. Springer, Cham; Centro Internazionale Matematico Estivo (C.I.M.E.), Florence, [2020] ©2020. Edited by Pierre Cardaliaguet and Alessio Porretta, Fondazione CIME/CIME Foundation Subseries.
- [2] H. W. Alt and L. A. Caffarelli. Existence and regularity for a minimum problem with free boundary. J. Reine Angew. Math., 325:105–144, 1981.
- [3] P. Andrade and E. A. Pimentel. Stationary fully nonlinear mean-field games. J. Anal. Math., 145(1):335–356, 2021.
- [4] P. Aviles and Y. Giga. A mathematical problem related to the physical theory of liquid crystal configurations. In Miniconference on geometry and partial differential equations, 2 (Canberra, 1986), volume 12 of Proc. Centre Math. Anal. Austral. Nat. Univ., pages 1–16. Austral. Nat. Univ., Canberra, 1987.
- [5] P. Aviles and Y. Giga. The distance function and defect energy. Proc. Roy. Soc. Edinburgh Sect. A, 126(5):923–938, 1996.
- [6] J. Bedrossian and R. V. Kohn. Blister patterns and energy minimization in compressed thin films on compliant substrates. Comm. Pure Appl. Math., 68(3):472–510, 2015.
- [7] A. Bensoussan, J. Frehse, and P. Yam. Mean field games and mean field type control theory. Springer Briefs in Mathematics. Springer, New York, 2013.
- [8] V. I. Bogachev, N. V. Krylov, M. Röckner, and S. V. Shaposhnikov. Fokker-Planck-Kolmogorov equations, volume 207 of Mathematical Surveys and Monographs. American Mathematical Society, Providence, RI, 2015.
- [9] L. Caffarelli. Interior a priori estimates for solutions of fully nonlinear equations. Ann. of Math. (2), 130(1):189–213, 1989.
-
[10]
L. Caffarelli, M. G. Crandall, M. Kocan, and A. Świ
‘ e - [11] L. Caffarelli and S. Salsa. A geometric approach to free boundary problems, volume 68 of Graduate Studies in Mathematics. American Mathematical Society, Providence, RI, 2005.
- [12] L. A. Caffarelli and X. Cabré. Fully nonlinear elliptic equations, volume 43 of American Mathematical Society Colloquium Publications. American Mathematical Society, Providence, RI, 1995.
- [13] P. Cardaliaguet. Notes on mean-field games. 2013.
- [14] S.-Y. A. Chang, M. J. Gursky, and P. C. Yang. Regularity of a fourth order nonlinear PDE with critical exponent. Amer. J. Math., 121(2):215–257, 1999.
- [15] S.-Y. A. Chang, L. Wang, and P. C. Yang. A regularity theory of biharmonic maps. Comm. Pure Appl. Math., 52(9):1113–1137, 1999.
- [16] I. Chowdhury, E. R. Jakobsen, and M. Krupski. On fully nonlinear parabolic mean field games with examples of nonlocal and local diffusions, 2021.
- [17] A. Cianchi and V. Maz’ya. Sobolev inequalities in arbitrary domains. Adv. Math., 293:644–696, 2016.
- [18] S. Conti and F. Maggi. Confining thin elastic sheets and folding paper. Arch. Ration. Mech. Anal., 187(1):1–48, 2008.
- [19] S. Conti, F. Maggi, and S. Müller. Rigorous derivation of Föppl’s theory for clamped elastic membranes leads to relaxation. SIAM J. Math. Anal., 38(2):657–680, 2006.
- [20] G. Dal Maso. An introduction to -convergence, volume 8 of Progress in Nonlinear Differential Equations and their Applications. Birkhäuser Boston, Inc., Boston, MA, 1993.
- [21] E. B. Fabes and D. W. Stroock. The -integrability of Green’s functions and fundamental solutions for elliptic and parabolic equations. Duke Math. J., 51(4):997–1016, 1984.
- [22] D. Gomes, E. Pimentel, and V. Voskanyan. Regularity theory for mean-field game systems. SpringerBriefs in Mathematics. Springer, [Cham], 2016.
- [23] R. V. Kohn. Energy-driven pattern formation. In International Congress of Mathematicians. Vol. I, pages 359–383. Eur. Math. Soc., Zürich, 2007.
- [24] R. V. Kohn and E. O’Brien. The wrinkling of a twisted ribbon. J. Nonlinear Sci., 28(4):1221–1249, 2018.
- [25] J.-M. Lasry and P.-L. Lions. Jeux à champ moyen. I. Le cas stationnaire. C. R. Math. Acad. Sci. Paris, 343(9):619–625, 2006.
- [26] J.-M. Lasry and P.-L. Lions. Jeux à champ moyen. II. Horizon fini et contrôle optimal. C. R. Math. Acad. Sci. Paris, 343(10):679–684, 2006.
- [27] J.-M. Lasry and P.-L. Lions. Mean field games. Jpn. J. Math., 2(1):229–260, 2007.
- [28] P.-L. Lions. Cours au collège de france. www.college-de-france.fr.
- [29] S. Venkataramani. Lower bounds for the energy in a crumpled elastic sheet—a minimal ridge. Nonlinearity, 17(1):301–312, 2004.
Julio C. Correa-Hoyos
Instituto de Matemática e Estatística
Universidade do Estado do Rio de Janeiro
20550-013, Maracanã, Rio de Janeiro - RJ, Brazil.
julio.correa@ime.uerj.br.
Edgard A. Pimentel (Corresponding Author)
University of Coimbra
CMUC, Department of Mathematics
3001-501 Coimbra, Portugal
edgard.pimentel@mat.uc.pt