Next Article in Journal
Modeling and Numerical Simulation for Covering the Fractional COVID-19 Model Using Spectral Collocation-Optimization Algorithm
Next Article in Special Issue
Fast Compact Difference Scheme for Solving the Two-Dimensional Time-Fractional Cattaneo Equation
Previous Article in Journal
Coefficient Estimates for the Functions with Respect to Symmetric Conjugate Points Connected with the Combination Binomial Series and Babalola Operator and Lucas Polynomials
Previous Article in Special Issue
Application of the Explicit Euler Method for Numerical Analysis of a Nonlinear Fractional Oscillation Equation
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Strong Solution for Fractional Mean Field Games with Non-Separable Hamiltonians

College of Mathematics and Statistics, Shenzhen University, Shenzhen 518060, China
*
Author to whom correspondence should be addressed.
Submission received: 20 May 2022 / Revised: 18 June 2022 / Accepted: 21 June 2022 / Published: 29 June 2022

Abstract

:
In this paper, we establish the existence and uniqueness of a strong solution to a fractional mean field games system with non-separable Hamiltonians, where the fractional exponent σ ( 1 2 , 1 ) . Our result is new for fractional mean field games with non-separable Hamiltonians, which generalizes the work of D.M. Ambrose for the integral case. The important step is to choose the new appropriate fractional order function spaces and use the Banach fixed-point theorem under stronger assumptions for the Hamiltonians.

1. Introduction

In this paper, we consider the following time-dependent fractional mean field games systems with non-separable Hamiltonians:
u t ( Δ ) σ u + H t , x , D u , m = 0 t . x 0 , T × T n ,
m t + ( Δ ) σ m + div m H p t , x , D u , m = 0 t . x 0 , T × T n ,
where T n is the tours R n / Z n , n 2 , u is a value function, m is a probability distribution, H = H ( t , x , D u , m ) is a Hamiltonian, H p denotes p H t , x , p , m , and σ ( 1 2 , 1 ) . Here, we define the fractional Laplacian ( Δ ) σ by the Fourier decomposition. For any u T n , if
u ( x ) = k Z n u ^ ( k ) e i k x with u ^ ( k ) = 1 ( 2 π ) n T n u ( x ) e i k x d x ,
where i 2 = 1 , then its fractional Laplacian is defined by
( Δ ) σ u = k Z n k 2 σ u ^ ( k ) e i k x .
We especially consider the initial-terminal problem of Equations (1) and (2) (i.e., the initial value of m and the terminal value of u are prescribed functions):
m 0 , x = m 0 x , u T , x = u T ( x ) .
The mean field games (MFG) system describes systems with very large numbers of identical agents in noncooperative differential games. Assume each agent wants to control his or her own trajectory in the same state space, which is affected by a stochastic differential equation:
d X t = α t d t + d W t ,
where W t is a stochastic noise. Meanwhile, each agent is rational and aims to minimize the following cost functional:
J ( t , x , α ) = E t T L X s , α s , m ( s ) d s + G X T , m ( T ) ,
where T > 0 is the finite horizon of the problem and L and G are given continuous maps. Originally, there were two different approaches to solving such a problem from two different points of view, which were proposed independently by Lasry and Lions [1,2,3] in the mathematics community and Huang, Malhamé, and Caines [4] in the engineering field almost at the same time. The main idea of MFG is to implement strategies based on the distribution of the other agents. In recent years, MFG theory has attracted more and more interest and been used in more fields widely. Lachapelle, Salomon, and Turinici [5] presented a model for the choice of insulation technology in households using an MFG model. They obtained an existing result for the associated optimization problem and gave a monotonic algorithm to find the mean field equilibria. Based on the SIR model, Lee et al. [6] introduced an effective mean-field game model for controlling the propagation of epidemics and provided fast numerical algorithms based on proximal primal-dual methods. There are more and more applications in engineering, finance, AI optimization problems, pandemic and vaccine control, etc. We refer the reader to [7,8,9,10] for a fairly large description of the current literature on the models and their applications.
If the stochastic noise W ( t ) in Equation (4) is a Brownian motion, we define the value function u as
u ( t , x ) = inf α J ( t , x , α ) .
By the dynamic programming principle, we formally obtain the integral order equation MFG system:
u t + Δ u + H t , x , D u , m = 0 ,
m t Δ m + div m H p t , x , D u , m = 0 ,
where the Hamiltonian H : R + × R n × R n × P ( R n ) R is defined by
H ( t , x , p , m ) : = sup a A L ( x , a , m ) p · b ( x , a , m ) .
The MFG system of partial differential equations takes the form of a backward Hamilton–Jacobi equation coupled with a forward Fokker–Planck equation. We say the system is separable if the Hamiltonian H can be written in the following form:
H t , x , p , m = H t , x , p + F t , x , m .
Otherwise, we say the system is non-separable. In the case of a separable Hamiltonian, MFG systems are well analyzed. The existence and uniqueness of smooth solutions for the second-order MFG system with non-local terms is given in [2,3]. The existence and uniqueness of weak solutions with local terms and different types of boundary conditions were proven in [11,12,13]. In particular, Porretta [11,12] obtained the existence and uniqueness of weak solutions to one kind of planning problem, where the equation is only prescribed with the initial and terminal conditions for the density m (i.e., m ( 0 ) = m 0 , m ( T ) = m 1 , where m 0 , m 1 > 0 are smooth functions). In [13], Porretta also developed a complete weak framework for the well-posed nature of weak solutions. The key work of his was to find new results for Fokker–Planck equations under minimal assumptions on the drift. We also refer to the MFG system with standard diffusion terms [9,14], degenerated diffusion terms [15], and first-order systems [16,17].
Recently, fractional MFG systems have also interested many researchers. In such cases, the stochastic noise W ( t ) in Equation (4) is modeled by a symmetric σ stable noise, which is a jump process or anomalous diffusion. In [18], the subcritical order of the fractional case σ ( 1 2 , 1 ) was studied in the case of local and non-local coupling between equations. In [19], the authors established the existence and uniqueness of solutions to the evolution fractional MFG system σ ( 0 , 1 ) with regularized coupling by the vanishing viscosity method. In [20], the authors studied the fractional and non-local parabolic MFG systems driven by jump Lévy processes in the whole space. They obtained the existence and uniqueness of classical solutions of MFG systems with local and non-local couplings for separable Hamiltonians.
In the case of non-separable Hamiltonians, many works are concerned with the standard diffusion term. The existence theorem of classical solutions by the Schauder fixed-point theorem is shown in [21]. The existence theorem of weak solutions of stationary problems was proven in [22,23]. In [24,25], the author proved the existence and uniqueness theorem for strong solutions to time-dependent problems, and in [26], the author also studied the existence and uniqueness theory for non-separable mean field games in Sobolev spaces. To the best of our knowledge, there are few theories involved fractional mean field games with non-separable Hamiltonians.
In this paper, we study the existence and uniqueness results for strong solutions of the fractional MFG system in the case of a non-separable Hamiltonian when σ ( 1 2 , 1 ) . Our main ideas stem from [24,25]. We adopt the similar function space based on Wiener algebra and the Banach fixed-point theorem for a short time of existence. While different from [24,25], we find that when the value function u and the probability distribution m are in the same function space, it becomes tricky to extend the result to the fractional Laplacian. Therefore, we consider enhancing the regularity of the value function u to a higher case than the probability distribution m as in [26]. In addition, as long as the non-separable Hamiltonian satisfies the assumption A 1 in Section 3, we can draw a conclusion: if the initial measure m 0 is close enough to its uniform measure on T n , and the terminal value u T is small enough, then the solutions to the problems in Equations (1)–(3) exists and are unique in a ball about the origin.
This paper is organized as follows. In Section 2, we show some preliminaries, including the property of the norm in the function space and the estimation of operators to the case of a fractional Laplacian. In Section 3, we give the main result and use the contraction principle to prove our first main theorem (Theorem 1). Finally, in Section 4, for the payoff boundary conditions, we make the assumption that the payoff function G satisfies the condition A 2 and obtain the second theorem of this paper.

2. Preliminaries

In this section, we state some necessary lemmas about function spaces as in [24,25] that will be used below. Indeed, our main work is to extend the integer order space to a fractional order function space.
Let T > 0 and α 0 , T 2 be given. As in [24,25], the function β : [ 0 , T ] [ 0 , α ] is defined by
β ( s ) = 2 α s / T , s [ 0 , T / 2 ] , 2 α 2 α s / T , s [ T / 2 , T ] .
Let s be a positive real number. The space B s consists of continuous functions f from T n to R such that the norm f B s is finite (i.e., the following is true):
f B s = k Z n 1 + | k | s | f ^ ( k ) | < .
In addition, similar to the function space defined in [24,25], an extended space-time version B α s consists of all functions in C ( [ 0 , T ] ; B s ) such that f B a s is finite (i.e., we have the following):
f B a s = k Z n sup t [ 0 , T ] 1 + | k | s e β ( t ) | k | | f ^ ( t , k ) | < .
Before we prove some properties of B α j , we need the following algebra property:
Lemma 1.
For any j N , there exists c j > 0 such that if f B α j and g B α j , then f g B α j with the estimate
f g B a j c j f B a j g B a j .
The proof of Lemma 1 follows along the same lines as that in [25], except for minor modifications which are omitted here.
Since the MFG system is a couple of backward and forward evolution equations, we introduce the operators I σ + and I σ for the fractional nonhomogeneous heat equation, which are defined as follows:
( I σ + h ) ( t , · ) = 0 t e Δ σ t s h ( s , · ) d s ,
( I σ h ) ( t , · ) = t T e Δ σ s t h ( s , · ) d s .
For any j N and σ ( 1 2 , 1 ) , we show that I σ + and I σ are bounded linear operators from B α j to B α j + 2 σ :
Lemma 2.
For any σ ( 1 2 , 1 ) , if h B α j with T d h ( t , x ) d x = 0 , then we have ( I σ + h ) ( t , · ) B α j + 2 σ and
( I σ + h ) ( t , · ) B α j + 2 σ C h ( t , · ) B α j .
Proof. 
Through Equation (9) and a convolution formula, we have
( I σ + h ^ ) ( t , · ) = 0 t e ( t s ) | k | 2 σ h ^ ( s , · ) d s ,
and
I σ + h B α j + 2 σ = k Z n 1 + | k | j + 2 σ sup t [ 0 , T ] e β ( t ) | k | 0 t e t s | k | 2 σ h ^ s , · d s .
Since T d h ( t , x ) d x = 0 , h ^ s , 0 = 0 , then we have
I σ + h B a j + 2 σ 2 k Z n \ { 0 } | k | j + 2 σ sup t [ 0 , T ] e β ( t ) | k | 0 t e | k | 2 σ ( t s ) e β ( s ) | k | e β ( s ) | k | h ^ ( s , k ) d s 2 k Z n \ { 0 } | k | j + 2 σ sup t [ 0 , T ] e β ( t ) | k | 0 t e | k | 2 σ ( t s ) e β ( s ) | k | sup τ [ 0 , T ] e β ( τ ) | k | | h ^ ( τ , k ) | d s 2 k Z n \ { 0 } | k | j sup τ [ 0 , T ] e β ( τ ) | k | | h ^ ( τ , k ) | [ sup t 0 , T ] | k | 2 σ e β ( t ) | k | 0 t e | k | 2 σ ( t s ) β ( s ) | k | d s ] 2 sup t [ 0 , T ] sup k Z n \ { 0 } | k | 2 σ e β ( t ) | k | 0 t e | k | 2 σ ( t s ) β ( s ) | k | d s h B α j .
Therefore, to verify that I σ + is a bounded linear operator between B α j and B α j + 2 σ , we have to prove that
sup t [ 0 , T ] sup k Z n \ { 0 } | k | 2 σ e β ( t ) | k | t | k | 2 σ 0 t e | k | 2 σ s β ( s ) | k | d s < .
First, when t 0 , T / 2 , we have β s = 2 α s / T by (7). Then, we obtain that
0 t e | k | 2 σ s β ( s ) | k | d s = 0 t e | k | 2 σ s 2 α s | k | / T d s = exp | k | 2 σ t 2 α t | k | / T 1 | k | 2 σ 2 α | k | / T .
Since α ( 0 , T / 2 ) and σ ( 1 2 , 1 ) , by Equation (12), we have
sup t [ 0 , T / 2 ] sup k Z n \ { 0 } | k | 2 σ e β ( t ) | k | t | k | 2 σ 0 t e | k | 2 σ s β ( s ) | k | d s = sup t [ 0 , T / 2 ] sup k Z n \ { 0 } | k | 2 σ exp 2 α t | k | / T t | k | 2 σ · exp | k | 2 σ t 2 α t | k | / T 1 | k | 2 σ 2 α | k | / T = sup t [ 0 , T / 2 ] sup k Z n \ { 0 } 1 exp 2 α t | k | / T t | k | 2 σ 1 2 α ( | k | 1 2 σ ) / T 1 1 2 α / T .
Next, we consider the case of t T / 2 , T . From Equation (7), we know β s = 2 α 2 α s / T . Therefore, we have
0 t e | k | 2 σ s β ( s ) | k | d s = 0 T / 2 e | k | 2 σ s β ( s ) | k | d s + T / 2 t e | k | 2 σ s β ( s ) | k | d s .
The second integral on the right side of (14) can be written as
T / 2 t e | k | 2 σ s β ( s ) | k | d s = T / 2 t e | k | 2 σ s 2 α | k | + 2 α s | k | / T d s = e 2 α | k | exp | k | 2 σ t + 2 α t | k | / T exp | k | 2 σ T / 2 + α | k | | k | 2 σ + 2 α | k | / T .
By combining Equation (12) with Equation (15), we have found that for t [ T / 2 , T ] , the following equality holds:
0 t e | k | 2 s β ( s ) | k | d s = exp | k | 2 σ t 2 α t | k | / T 1 | k | 2 σ 2 α | k | / T + e 2 α | k | exp | k | 2 σ t + 2 α t | k | / T exp | k | 2 σ T / 2 + α | k | | k | 2 σ + 2 α | k | / T .
When t T / 2 , T , we set
sup t T / 2 , T ] sup k Z n \ { 0 } | k | 2 σ e β ( t ) | k | t | k | 2 σ 0 t e | k | 2 σ s β ( s ) | k | d s I + I I .
where I and I I are given by
I = sup t [ T / 2 , T ] sup k Z n \ { 0 } | k | 2 σ exp ( 2 α | k | 2 α | k | t / T t | k | 2 σ ) exp ( | k | 2 σ T / 2 α | k | ) 1 | k | 2 σ 2 α | k | / T , I I = sup t [ T / 2 , T ] sup k Z n \ { 0 } | k | 2 σ exp ( 2 α | k | t / T t | k | 2 σ ) · exp ( | k | 2 σ t + 2 α t | k | / T ) exp ( | k | 2 σ T / 2 + α | k | ) | k | 2 σ + 2 α | k | / T .
Finally, we will prove that I and I I are bounded. It is obvious that
I = sup t [ 0 , T / 2 ] sup k Z n \ { 0 } exp ( 2 α | k | 2 α | k | t / T t | k | 2 σ ) 1 2 α | k | 1 2 σ / T exp ( | k | 2 σ T / 2 α | k | ) 1 = sup t [ 0 , T / 2 ] sup k Z n \ { 0 } exp ( α | k | 2 α | k | t / T ( t T 2 ) | k | 2 σ ) exp ( 2 α | k | 2 α | k | t / T t | k | 2 σ ) 1 2 α | k | 1 2 σ / T .
Using the conditions t T / 2 , σ ( 1 2 , 1 ) , α 0 , T / 2 , and k Z n \ { 0 } , we have
exp | k | 2 σ T / 2 α | k | 1 > 0 , α | k | 2 α | k | t / T ( t T 2 ) | k | 2 σ < 0 .
Then, we have
I sup k Z n \ { 0 } 1 1 2 α | k | 1 2 σ / T 1 1 2 α / T = T T 2 α .
Now, let us consider I I . As in the above discussion, we have
1 exp | k | 2 σ T 2 t + α | k | 1 2 t T > 0 .
and
I I = sup t [ T / 2 , T ] sup k Z n 1 1 + 2 α | k | 1 2 σ / T 1 exp | k | 2 σ T 2 t + α | k | 1 2 t T sup k Z n \ { 0 } 1 1 + 2 α | k | 1 2 σ / T = 1 .
Based on Equations (13), (16) and (17), we proved that I σ + is a bounded linear operator and the operator norm satisfies
I σ + B α j B α j + 2 σ 2 T T 2 α + 2 .
Regarding I σ , we have a similar result, and we omitted its proof.
Lemma 3.
For any σ ( 1 2 , 1 ) , if h B α j with T d h ( t , x ) d x = 0 , then we have ( I σ h ) ( t , · ) B α j + 2 σ and
( I σ h ) ( t , · ) B α j + 2 σ C h ( t , · ) B α j .
At last, we have the following lemma:
Corollary 1.
For any σ ( 1 2 , 1 ) , if f B 2 σ , then
sup t [ 0 , T ] e Δ σ t f B α 2 σ + sup t [ 0 , T ] e Δ σ T t f B α 2 σ C f B 2 σ .
Proof. 
Using the definition B 2 σ , we have
sup t [ 0 , T ] e Δ σ t f B α 2 σ = sup t [ 0 , T ] k Z n 1 + | k | 2 σ e β ( t ) | k | e t | k | 2 σ f ^ k Z n 1 + | k | 2 σ f ^ sup t [ 0 , T ] e β ( t ) | k | t | k | 2 σ C | f | B 2 σ .
Using similar methods, we can find e Δ σ T t f B α 2 σ . □

3. Main Result and Its Proof

In this section, we first state the definition of the strong solution for the problems in Equations (1)–(3). Then, we show the existence and uniqueness theorem and give its proof.

3.1. Strong Solution Formulation

Let P be the projection operator, which removes the mean value of the periodic function (i.e., the following is true):
P f = f 1 vol T n T n f ( x ) d x .
We define w = P u and μ = P m = m m ¯ , where m ¯ = 1 / vol T n . Since m is a probability measure at each time, its integral in the spatial domain will always be equal to one. Then, from Equation (1), we can find the evolution equation of w:
w t Δ σ w + P H ( t , x , D w , μ ) = 0 ,
Additionally, from Equation (2), we have the following evolution equation for μ :
μ t + Δ σ μ + div μ H p ( t , x , D w , μ ) + m ¯ div H p ( t , x , D w , μ ) = 0 .
Of course, the initial-terminal condition is taken as w T , · = w T : = P u T , μ ( 0 , · ) = μ 0 : = m 0 m ¯ .
Before writing the Duhamel formula for w and μ . We assume that P H can be expressed as
P H = P b ( t , x ) μ + P Υ t , x , D w , μ .
The conditions satisfied by b and Υ will be given below.
We then have the following Duhamel formula for μ :
μ t , · = e Δ σ t μ 0 + I σ + div μ H p · , · , D w , μ t , · + m ¯ ( I σ + ( div ( H p ( · , · , D w , μ ) ) ) ) ( t , · ) ,
In addition, by integrating backward in time, we find
w t , · = e Δ σ T t w T I σ P Υ · , · , D w , μ I σ P b μ t , · .
Indeed, by substituting Equation (19) into Equation (20), we obtain the following equation for w:
w t , · = e Δ σ T t w T I σ ( P Υ ( · , · , D w , μ ) ) I σ ( P ( b e ( Δ ) σ · μ 0 ) ) ( t ) I σ P ( b I σ + div ( μ H p ( · , · , D w , μ ) ) ( · ) ) ( t ) m ¯ I σ P ( b I + div ( H p ( · , · , D w , μ ) ) ( · ) ) ( t ) .
We say ( u , m ) is a strong solution of the problems in Equations (1)–(3) if ( w , μ ) B α 2 σ + 1 × B α 2 σ such that the integral equalities in Equations (21) and (19) hold.

3.2. Main Theorem

In order to obtain the existence of strong solutions for the problems in Equations (1)–(3), we assume that the Hamiltonian H satisfies the following assumption:
(A1) Let σ ( 1 2 , 1 ) , Υ ( · , · , 0 , 0 ) = 0 and H p ( · , · , 0 , 0 ) = 0 . There exists a continuous function Φ 1 : B α 2 σ + 1 × B α 2 σ R such that as w 1 , w 2 , μ 1 , μ 2 0 , we have Φ 1 w 1 , w 2 , μ 1 , μ 2 0 and
H p ( · , · , D w 1 , μ 1 ) H p ( · , · , D w 2 , μ 2 ) B α 1 n Φ 1 ( w 1 , w 2 , μ 1 , μ 2 ) D w 1 D w 2 B α 2 σ n + μ 1 μ 2 B α 2 σ .
There exists a continuous function Φ 2 : B α 2 σ + 1 × B α 2 σ R so that as w 1 , w 2 , μ 1 , μ 2 0 , we have Φ 2 w 1 , w 2 , μ 1 , μ 2 0 and
P Υ ( · , · , D w 1 , μ 1 ) P Υ ( · , · , D w 2 , μ 2 ) B α 1 Φ 2 ( w 1 , w 2 , μ 1 , μ 2 ) D w 1 D w 2 ( B α 2 σ ) n + μ 1 μ 2 B α 2 σ .
From the preparations above, we now give the first main theorem. Our proof depends on the Banach fixed-point theorem as in [25]:
Theorem 1.
Let T > 0 and α ( 0 , T / 2 ) be given. Let assumption ( A 1 ) be satisfied, and assume b B α 1 . There exists δ > 0 such that if u T and the probability measure m 0 are such that w T = P u T and μ 0 = m 0 m ¯ satisfy w T B 2 σ + 1 + μ 0 B 2 σ < δ , then the system in Equations (1)–(3) has a strong, locally unique solution ( u , m ) B α 2 σ + 1 × B α 2 σ .
Proof. 
Let ( a 0 , b 0 ) be
( a 0 , b 0 ) = e ( Δ ) σ ( T t ) w T I σ ( P ( b e ( Δ ) σ · μ 0 ) ) t , e ( Δ ) σ t μ 0 .
We set X to be the closed ball in B α 2 σ + 1 × B α 2 σ centered at a point ( a 0 , b 0 ) with a radius r * , where r * will be determined later:
X = ( w , μ ) B α 2 σ + 1 × B α 2 σ : w a 0 B 2 σ + 1 + μ b 0 B 2 σ < r * .
We define a mapping T ( w , μ ) = ( T 1 ( w , μ ) , T 2 ( w , μ ) ) on X based on Equations (21) and (20) (i.e., we have the following):
T 1 ( w , μ ) = e ( Δ ) σ ( T t ) w T I σ ( P Υ ( · , · , D w , μ ) ) I σ ( P ( b e ( Δ ) σ · μ 0 ) ) ( t ) I σ ( P ( b I σ + div ( μ H p ( · , · , D w , μ ) ) ( · ) ) ) ( t ) m ¯ I σ ( P ( b I σ + div ( H p ( · , · , D w , μ ) ) ( · ) ) ) ( t ) ,
T 2 ( w , μ ) = e ( Δ ) σ t μ 0 + I σ + ( div ( μ H p ( · , · , D w , μ ) ) ) ( t , · ) + m ¯ ( I σ + ( div ( H p ( · , · , D w , μ ) ) ) ) ( t , · ) .
First, we prove that T maps X to X. Because of μ 0 B 2 σ , it is easy to know that b 0 = e Δ σ t μ 0 B α 2 σ from the inference in Corollary 1. By the definition of a 0 in Equation (24) and the property of I σ given earlier in the previous section, we have
a 0 B α 2 σ + 1 = e ( Δ ) σ ( T t ) w T I σ ( P ( b e ( Δ ) σ · μ 0 ) ) ( t ) B α 2 σ + 1 e ( Δ ) σ ( T t ) w T B α 2 σ + 1 + I σ ( P ( b e ( Δ ) σ · μ 0 ) ) ( t ) B α 2 σ + 1 e ( Δ ) σ ( T t ) w T B α 2 σ + 1 + C b B α 1 e ( Δ ) σ · μ 0 B α 2 σ .
Using b B α 1 and w T B 2 σ + 1 , it is easy to know that a 0 B α 2 σ + 1 .
To demonstrate that T maps X to X, we will only need to show that
T 1 ( w , μ ) a 0 B α 2 σ + 1 r * / 2 and T 2 ( w , μ ) b 0 B α 2 σ r * / 2 .
We begin with T 1 . Recalling the definition in Equation (25), we only need to prove the following inequality:
I σ ( P Υ ( · , · , D w , μ ) ) B α 2 σ + 1 r * / 6 ,
I σ P ( b I + div ( μ H p ( · , · , D w , μ ) ) ( · ) ) ( t ) B α 2 σ + 1 r * / 6 ,
m ¯ I σ P ( b I + div ( H p ( · , · , D w , μ ) ) ( · ) ) ( t ) B α 2 σ + 1 r * / 6 .
By assuming Equation (22) in ( A 1 ) and the property of I σ , we have
I σ ( P Υ ( · , · , D w , μ ) ) B α 2 σ + 1 C P Υ ( · , · , D w , μ ) B α 1 C Φ 2 ( w , 0 , μ , 0 ) ( D w B α 2 σ + μ B α 2 σ ) C Φ 2 ( w , 0 , μ , 0 ) ( w B α 2 σ + 1 + μ B α 2 σ ) .
Since w X and μ X , we have
w B α 2 σ + 1 + μ B α 2 σ a 0 B α 2 σ + 1 + b 0 B α 2 σ + r * .
Let μ 0 and w T be small enough that
a 0 B α 2 σ + 1 + b 0 B α 2 σ r * .
By combining Equations(30), (38) and (31), we have
I σ ( P Υ ( · , · , D w , μ ) ) B α 2 σ + 1 2 C r * Φ 2 ( w , 0 , μ , 0 ) .
Since Φ 2 is continuous and Φ 2 ( 0 , 0 , 0 , 0 ) = 0 , we may have μ 0 , w T and r * at small enough values such that
max ( w , μ ) X Φ 2 ( w , 0 , μ , 0 ) 1 12 C r * .
Then, Equation (32) implies Equation (27).
Next, we demonstrate Equation (28). By the properties of I σ and I σ + , and the fact that the divergence is a first-order operator, combined with the algebra porperty of B α j in Lemma 1, we have
I σ ( P ( b I + div ( μ H p ( · , · , D w , μ ) ) ( · ) ) ) ( t ) B α 2 σ + 1 C ( P ( b I + div ( μ H p ( · , · , D w , μ ) ) ( · ) ) ) ( t ) B α 1 C b B α 1 ( I + div ( μ H p ( · , · , D w , μ ) ) ( · ) ) ( t ) B α 1 C b B α 1 ( I + div ( μ H p ( · , · , D w , μ ) ) ( · ) ) ( t ) B α 2 σ C b B α 1 ( div ( μ H p ( · , · , D w , μ ) ) ( · ) ) ( t ) B α 0 C b B α 1 μ B α 1 H p ( · , · , D w , μ ) ( · ) ( t ) B α 1 .
Since 1 2 < σ < 1 , w X , and μ X , we have
μ B α 1 μ B α 2 σ a 0 B α 2 σ + 1 + b 0 B α 2 σ + r * .
From assumption (23) in A 1 , we have
H p · , · , D w , μ ( · ) ( t ) B α 1 Φ 1 w , 0 , μ , 0 D w B α 2 σ + μ B α 2 σ Φ 1 w , 0 , μ , 0 w B α 2 σ + 1 + μ B α 2 σ Φ 1 w , 0 , μ , 0 a 0 B α 2 σ + 1 + b 0 B α 2 σ + r * .
We again require that μ 0 and w T are small enough that
a 0 B α 2 σ + 1 + b 0 B α 2 σ r * .
By combining Equations (33) and (34) with Equation (35), we have
I σ P b I + div μ H p · , · , D w , μ ( · ) ( t ) B α 2 σ + 1 4 C ( r * ) 2 max w , μ X Φ 1 w , 0 , μ , 0 b B α 1 .
Thus, if r * is chosen to be a value small enough that
4 C ( r * ) 2 max w , μ X Φ 1 w , 0 , μ , 0 b B α 1 r * 6 ,
then Equation (36) implies Equation (28).
Finally, we will demonstrate Equation (29). Similar to the proof of Equation (28), we have
m ¯ I σ ( P ( b I + div ( H p ( · , · , D w , μ ) ) ( · ) ) ) ( t ) B α 2 σ + 1 C m ¯ ( P ( b I + div ( H p ( · , · , D w , μ ) ) ( · ) ) ) ( t ) B α 1 C m ¯ b B α 1 ( I + div ( H p ( · , · , D w , μ ) ) ( · ) ) ( t ) B α 1 C m ¯ b B α 1 H p ( · , · , D w , μ ) ( · ) ( t ) B α 1 C m ¯ b B α 1 max ( w , μ ) X Φ 1 ( w , 0 , μ , 0 ) ( a 0 B α 2 σ + 1 + b 0 B α 2 σ + r * ) .
If μ 0 and w T are chosen to be values small enough that
a 0 B α 2 σ + 1 + b 0 B α 2 σ r * and max w , μ X Φ 1 w , 0 , μ , 0 1 12 C m ¯ b B α 1 ,
then Equation (37) implies Equation (29).
Next, we consider T 2 . Recalling the definition in Equation (26), we only need to prove the following inequalities:
I σ + div μ H p ( · , · , D w , μ ) ( t , · ) B α 2 σ r * 4 ,
m ¯ ( I σ + ( div ( H p ( · , · , D w , μ ) ) ) ) ( t , · ) B α 2 σ r * 4 .
We begin by establishing Equation (38). It is obvious that
I σ + div μ H p ( · , · , D w , μ ) ( t , · ) B α 2 σ C ( div ( μ H p ( · , · , D w , μ ) ) ) ( t , · ) B α 0 C μ B α 1 H p ( · , · , D w , μ ) ( t , · ) B α 1 .
Adopting a method similar to that in Equation (28), we know if r * is chosen to be small enough that
4 C ( r * ) 2 max ( w , μ ) X Φ 1 ( w , 0 , μ , 0 ) r * 4 ,
then Equation (40) implies Equation (38).
Then, we demonstrate Equation (39). It is immediate that
m ¯ ( I σ + ( div ( H p ( · , · , D w , μ ) ) ) ) ( t , · ) B α 2 σ C m ¯ H p ( · , · , D w , μ ) ( t , · ) B α 1 C m ¯ max ( w , μ ) X Φ 1 ( w , 0 , μ , 0 ) ( a 0 B α 2 σ + 1 + b 0 B α 2 σ + r * ) .
If r * , μ 0 and w T are chosen to be small enough that
max ( w , μ ) X Φ 1 ( w , 0 , μ , 0 ) min 1 12 C m ¯ b B α 1 , 1 8 C m ¯ ,
then Equation (41) implies Equation (39). Therefore, we proved that when choosing r * , μ 0 , and w T appropriately, T maps X to X.
Now, we are ready to demonstrate the contraction estimate (i.e., we will demonstrate that if r * , μ 0 , and w T are sufficiently small, then there exists λ ( 0 , 1 ) such that for all ( w 1 , μ 1 ) X and ( w 2 , μ 2 ) X , we have the following):
T w 1 , μ 1 T w 2 , μ 2 B α 2 σ + 1 × B α 2 σ λ ( w 1 w 2 B α 2 σ + 1 + μ 1 μ 2 B α 2 σ ) .
Considering the definitions in Equations (25) and (26), we will demonstrate that
T 1 w 1 , μ 1 T 1 w 2 , μ 2 B α 2 σ + 1 λ 1 w 1 w 2 B α 2 σ + 1 + μ 1 μ 2 B α 2 σ , T 2 w 1 , μ 1 T 2 w 2 , μ 2 B α 2 σ λ 2 w 1 w 2 B α 2 σ + 1 + μ 1 μ 2 B α 2 σ .
It can be seen that repeated use of the triangle inequality implies that it is sufficient to establish the following bounds:
I σ P Υ · , · , D w 1 , μ 1 Υ · , · D w 2 , μ 2 B α 2 σ + 1 1 6 w 1 w 2 B α 2 σ + 1 + μ 1 μ 2 B α 2 σ , I σ P b I + div μ 1 H p · , · D w 1 , μ 1 μ 2 H p · , · , D w 2 , μ 2 B α 2 σ + 1
1 6 w 1 w 2 B α 2 σ + 1 + μ 1 μ 2 B α 2 σ , m ¯ I σ P b I σ + div H p · , · D w 1 , μ 1 H p · , · , D w 2 , μ 2 B 2 σ + 1
1 6 w 1 w 2 B α 2 σ + 1 + μ 1 μ 2 B α 2 σ , I σ + div μ 1 H p · , · , D w 1 , μ 1 μ 2 H p · , · , D w 2 , μ 2 B α 2 σ
1 6 w 1 w 2 B α 2 σ + 1 + μ 1 μ 2 B α 2 σ , m ¯ I σ + div H p · , · D w 1 , μ 1 H p · , · D w 2 , μ 2 B α 2 σ
1 6 w 1 w 2 B α 2 σ + 1 + μ 1 μ 2 B α 2 σ .
We begin by establishing Equation (43). By following the assumptions and using the tools from before, such as triangle inequalities, the properties of B α j , and the mapping properties of I σ + and I σ , we have
I σ ( P ( Υ ( · , · , D w 1 , μ 1 ) Υ ( · , · D w 2 , μ 2 ) ) ) B α 2 σ + 1 C P ( Υ ( · , · , D w 1 , μ 1 ) Υ ( · , · D w 2 , μ 2 ) ) B α 1 C Φ 2 ( w 1 , w 2 , μ 1 , μ 2 ) ( D w 1 D w 2 B α 2 σ + μ 1 μ 2 B α 2 σ ) C Φ 2 ( w 1 , w 2 , μ 1 , μ 2 ) ( w 1 w 2 B α 2 σ + 1 + μ 1 μ 2 B α 2 σ ) .
Since Φ 2 is continuous and Φ 2 ( 0 , 0 , 0 , 0 ) = 0 , if r * , μ 0 , and w T are small enough such that
max ( w 1 , μ 1 ) , ( w 2 , μ 2 ) X Φ 2 ( w 1 , w 2 , μ 1 , μ 2 ) 1 6 C ,
then Equation (48) implies Equation (43).
Now, we will demonstrate Equation (44). Similar to the above proof, we have the following facts:
I σ P b I + ( div ( μ 1 H p ( · , · D w 1 , μ 1 ) μ 2 H p ( · , · , D w 2 , μ 2 ) ) ) B α 2 σ + 1 C b B α 1 ( μ 1 H p ( · , · D w 1 , μ 1 ) μ 2 H p ( · , · , D w 2 , μ 2 ) ) B α 1 C b B α 1 ( μ 1 H p ( · , · D w 1 , μ 1 ) μ 1 H p ( · , · , D w 2 , μ 2 ) ) B α 1 + C ( μ 1 H p ( · , · D w 2 , μ 2 ) μ 2 H p ( · , · , D w 2 , μ 2 ) ) B α 1 C b B α 1 ( μ 1 B α 1 H p ( · , · D w 1 , μ 1 ) H p ( · , · , D w 2 , μ 2 ) B α 1 + μ 1 μ 2 B α 1 H p ( · , · , D w 2 , μ 2 ) B α 1 ) C r * b B α 1 ( H p ( · , · D w 1 , μ 1 ) H p ( · , · , D w 2 , μ 2 ) B α 1 + C r * max ( w , μ ) X Φ 1 ( w , 0 , μ , 0 ) μ 1 μ 2 B α 2 σ ) C r * b B α 1 max ( w 1 , μ 1 ) , ( w 2 , μ 2 ) X Φ 1 w 1 , w 2 , μ 1 , μ 2 w 1 w 2 B α 2 σ + 1 + C r * b B α 1 max ( w 1 , μ 1 ) , ( w 2 , μ 2 ) X Φ 1 ( w 1 , w 2 , μ 1 , μ 2 ) + max w , μ X Φ 1 w , 0 , μ , 0 μ 1 μ 2 B α 2 σ .
Since Φ 1 is continuous and Φ 1 ( 0 , 0 , 0 , 0 ) = 0 , if r * , μ 0 , and w T are small enough, we can obtain
max ( w 1 , μ 1 ) , ( w 2 , μ 2 ) X Φ 1 ( w 1 , w 2 , μ 1 , μ 2 ) + max ( w , μ ) X Φ 1 ( w , 0 , μ , 0 ) 1 6 C r * b B α 1 .
Then, Equation (49) implies Equation (44).
Next, we will demonstrate Equation (45). Similar to the above proof, we have the following facts:
m ¯ I σ ( P ( b I σ + ( div ( H p ( · , · D w 1 , μ 1 ) H p ( · , · , D w 2 , μ 2 ) ) ) ) ) B α 2 σ + 1 C m ¯ b B α 1 ( H p ( · , · D w 1 , μ 1 ) H p ( · , · , D w 2 , μ 2 ) ) B α 1 C m ¯ b B α 1 max w 1 , μ 1 , w 2 , μ 2 X Φ 1 w 1 , w 2 , μ 1 , μ 2 w 1 w 2 B α 2 σ + 1 + μ 1 μ 2 B α 2 σ .
If r * , μ 0 , and w T are small enough such that
max ( w 1 , μ 1 ) , ( w 2 , μ 2 ) X Φ 1 ( w 1 , w 2 , μ 1 , μ 2 ) min 1 6 C m ¯ b B α 1 , 1 12 C r * b B α 1 ,
then Equation (50) implies Equation (45).
Our task now moves to proving Equation (46). It is immediate that
I σ + div ( μ 1 H p ( · , · , D w 1 , μ 1 ) μ 2 H p ( · , · , D w 2 , μ 2 ) ) B α 2 σ C μ 1 H p · , · , D w 1 , μ 1 μ 2 H p · , · , D w 2 , μ 2 B α 1 C μ 1 B α 1 H p · , · , D w 1 , μ 1 H p · , · , D w 2 , μ 2 B α 1 + C H p · , · , D w 2 , μ 2 B α 1 μ 1 μ 2 B α 1 C r * max ( w 1 , μ 1 ) , ( w 2 , μ 2 ) X Φ 1 ( w 1 , w 2 , μ 1 , μ 2 ) w 1 w 2 B 2 σ + 1 + C r * max ( w 1 , μ 1 ) , ( w 2 , μ 2 ) X Φ 1 ( w 1 , w 2 , μ 1 , μ 2 ) + max ( w , μ ) X Φ 1 ( w , 0 , μ , 0 ) μ 1 μ 2 B α 2 σ .
If r * , μ 0 , and w T are small enough, we will have
max ( w 1 , μ 1 ) , ( w 2 , μ 2 ) X Φ 1 ( w 1 , w 2 , μ 1 , μ 2 ) min 1 6 C m ¯ b B α 1 , 1 12 C r * b B α 1 , 1 12 C r * .
Then, Equation (51) implies Equation (46).
Finally we have to show Equation (47). It is easy to find that
m ¯ I σ + ( div ( H p ( · , · D w 1 , μ 1 ) H p ( · , · D w 2 , μ 2 ) ) ) B α 2 σ m ¯ H p · , · D w 1 , μ 1 H p · , · D w 2 , μ 2 B α 1 C m ¯ max ( w 1 , μ 1 ) , ( w 2 , μ 2 ) X Φ 1 ( w 1 , w 2 , μ 1 , μ 2 ) ( w 1 w 2 B α 2 σ + 1 + μ 1 μ 2 B α 2 σ ) .
If r * , μ 0 and w T are small enough, we have
max ( w 1 , μ 1 ) , ( w 2 , μ 2 ) X Φ 1 ( w 1 , w 2 , μ 1 , μ 2 ) min 1 6 C m ¯ b B α 1 , 1 12 C r * b B α 1 , 1 12 C r * , 1 6 C m ¯ .
Then, Equation (52) implies Equation (47).
We established the estimates in Equations (43)–(47) and successfully proved the constant λ = 5 6 in Equaiton (42). This completes the proof of Theorem 1. □

4. The Payoff Problem

The above formulation and existence theorem can be readily adapted to the payoff problem, in which Equation (3) is replaced by
m 0 , x = m 0 x , u T , x = G m T , · .
The needed modification in the formulation is that Equation (21) is replaced with
w t , · = e Δ σ T t P G μ ( T , · ) I σ P Υ · , · , D w , μ I σ P b e Δ σ · μ 0 t I σ P b I σ + div μ H p · , · , D w , μ · t m ¯ I σ P b I σ + div H p · , · , D w , μ · t .
In addition, we need to make assumptions about the payoff function P G :
( A 2 ) P G ( 0 ) = 0 , and P G is in the neighborhood of the origin in B 2 σ + 1 . Specifically, we assume that there exists c > 0 and ϵ > 0 such that for all a 1 , a 2 satisfies a 1 B 2 σ < ϵ and a 2 B 2 σ < ϵ :
P G a 1 P G a 2 B 2 σ + 1 c a 1 a 2 B 2 σ .
For example, G ( a ) = Δ 1 2 a and thus certainly satisfies this assumption:
Theorem 2.
Let T > 0 and α 0 , T / 2 be given. Let assumptions ( A 1 ) and ( A 2 ) be satisfied, and assume b B α 1 . There exists δ > 0 such that if u T and the probability measure m 0 are such that w T = P u T and μ 0 = m 0 m ¯ satisfy | μ 0 | B 2 σ < δ , then the system in Equations (1), (2) and (53) has a strong, locally unique solution u , m B α 2 σ + 1 × B α 2 σ .
The proof of this result is quite similar to that given earlier for Theorem 1, and thus it was omitted.

Author Contributions

Conceptualization, Q.L. and H.Y.; writing—original draft preparation, W.Z.; writing—review and editing, Q.L. and W.Z.; supervision, H.Y. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by NNSFC of Guang Dong, numbers 2020A1515010554 and 2022A1515010348.

Data Availability Statement

Not applicable.

Acknowledgments

The authors would like to express their sincere thanks to the referees and the editor for their helpful comments on the original version of this paper, which greatly improved the quality of this paper.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Lasry, J.-M.; Lions, P.-L. Jeux à champ moyen. I-le cas stationnaire. C. R. Math. Acad. Sci. Paris 2006, 343, 619–625. [Google Scholar] [CrossRef]
  2. Lasry, J.-M.; Lions, P.-L. Jeux à champ moyen. II-Horizon fini et contrôle optimal. C. R. Math. Acad. Sci. Paris 2006, 343, 679–684. [Google Scholar] [CrossRef]
  3. Lasry, J.-M.; Lions, P.-L. Mean field games. Jpn. J. Math. 2007, 2, 229–260. [Google Scholar] [CrossRef] [Green Version]
  4. Huang, M.; Caines, P.E.; Malham, R.P. Large-population cost-coupled LQG problems with nonuniform agents: Individual-mass behavior and decentralized ε-Nash equilibria. IEEE Trans. Automat. Contr. 2007, 52, 1560–1571. [Google Scholar] [CrossRef]
  5. Lachapelle, A.; Salomon, J.; Turinici, G. Computation of mean field equilibria in economics. Math. Mod. Meth. Appl. Sci. 2010, 20, 567–588. [Google Scholar] [CrossRef] [Green Version]
  6. Lee, W.; Liu, S.; Tembine, H.; Li, W.; Osher, S. Controlling Propagation of Epidemics via Mean-Field Control. SIAM J. Appl. Math. 2021, 81, 190–207. [Google Scholar] [CrossRef]
  7. Banez, R.A.; Li, L.; Yang, C.; Han, Z. Mean Field Game and Its Applications in Wireless Networks; Springer: New York, NY, USA, 2021. [Google Scholar]
  8. Barreiro-Gomez, J.; Tembine, H. Mean-Field-Type Games for Engineers; Taylor& Franis Group: Oxfordshire, UK, 2021. [Google Scholar]
  9. Guéant, O.; Lasry, J.M.; Lions, P.L. Mean field games and applications. In Paris-Princeton Lectures on Mathematical Finance; Springer: Berlin/Heidelberg, Germany, 2011; pp. 205–266. [Google Scholar]
  10. Huang, M.; Malham, R.P.; Caines, P.E. Large population stochastic dynamic games: Closed-loop McKean-Vlasov systems and the Nash certainty equivalence principle. Commun. Inf. Syst. 2006, 6, 221–252. [Google Scholar]
  11. Porretta, A. On the planning problem for a class of mean field games. C. R. Math. Acad. Sci. Paris 2013, 351, 457–462. [Google Scholar] [CrossRef]
  12. Porretta, A. On the planning problem for the mean field games system. Dyn. Games Appl. 2014, 4, 231–256. [Google Scholar] [CrossRef]
  13. Porretta, A. Weak solutions to Fokker-Planck equations and mean field games. Arch. Rat. Mech. Anal. 2015, 216, 1–62. [Google Scholar] [CrossRef]
  14. Gomes, D.A.; Pimentel, E.A.; Voskanyan, V. Regularity Theory for Mean-Field Game Systems; Springer: New York, NY, USA, 2016. [Google Scholar]
  15. Cardaliaguet, P.; Graber, P.J.; Porretta, A.; Tonon, D. Second order mean field games with degenerate diffusion and local coupling. Nonlinear Differ. Equ. Appl. 2015, 22, 1287–1317. [Google Scholar] [CrossRef] [Green Version]
  16. Cardaliaguet, P.; Graber, P.J. Mean field games systems of first order. ESAIM Control Optim. Calc. Var. 2015, 21, 690–722. [Google Scholar] [CrossRef] [Green Version]
  17. Cardaliaguet, P.; Porretta, A.; Tonon, D. Sobolev regularity for the first order Hamilton-Jacobi equation. Calc. Var. Partial. Differ. Equ. 2015, 54, 3037–3065. [Google Scholar] [CrossRef] [Green Version]
  18. Cesaroni, A.; Cirant, M.; Dipierro, S.; Novaga, M.; Valdinoci, E. On stationary fractional mean field games. J. Math. Pures Appl. 2019, 122, 1–22. [Google Scholar] [CrossRef] [Green Version]
  19. Cirant, M.; Goffi, A. On the existence and uniqueness of solutions to time-dependent fractional MFG. SIAM J. Math. Anal. 2019, 51, 913–954. [Google Scholar] [CrossRef]
  20. Ersland, O.; Jakobsen, E.R. On fractional and nonlocal parabolic mean field games in the whole space. J. Differ. Equ. 2021, 301, 428–470. [Google Scholar] [CrossRef]
  21. Cardaliaguet, P.; Porretta, A. An Introduction to Mean Field Game Theory. In Mean Field Games; Springer: Berlin, Germany, 2020; pp. 1–158. [Google Scholar]
  22. Ferreira, R.; Gomes, D. Existence of weak solutions to stationary mean-field games through variational inequalities. SIAM J. Math. Anal. 2018, 50, 5969–6006. [Google Scholar] [CrossRef]
  23. Gomes, D.A.; Patrizi, S.; Voskanyan, V. On the existence of classical solutions for stationary extended mean field games. Nonlinear Anal. 2014, 99, 49–79. [Google Scholar] [CrossRef] [Green Version]
  24. Ambrose, D.M. Small strong solutions for time-dependent mean field games with local coupling. C. R. Math. Acad. Sci. Paris 2016, 354, 589–594. [Google Scholar] [CrossRef]
  25. Ambrose, D.M. Strong solutions for time-dependent mean field games with non-separable Hamiltonians. J. Math. Pures Appl. 2018, 113, 141–154. [Google Scholar] [CrossRef] [Green Version]
  26. Ambrose, D.M. Existence theory for non-separable mean field games in Sobolev spaces. Indiana U. Math. J. 2022, 71, 611–647. [Google Scholar] [CrossRef]
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Ye, H.; Zou, W.; Liu, Q. Strong Solution for Fractional Mean Field Games with Non-Separable Hamiltonians. Fractal Fract. 2022, 6, 362. https://0-doi-org.brum.beds.ac.uk/10.3390/fractalfract6070362

AMA Style

Ye H, Zou W, Liu Q. Strong Solution for Fractional Mean Field Games with Non-Separable Hamiltonians. Fractal and Fractional. 2022; 6(7):362. https://0-doi-org.brum.beds.ac.uk/10.3390/fractalfract6070362

Chicago/Turabian Style

Ye, Hailong, Wenzhong Zou, and Qiang Liu. 2022. "Strong Solution for Fractional Mean Field Games with Non-Separable Hamiltonians" Fractal and Fractional 6, no. 7: 362. https://0-doi-org.brum.beds.ac.uk/10.3390/fractalfract6070362

Article Metrics

Back to TopTop