Next Article in Journal
State of Health Prediction of Electric Vehicles’ Retired Batteries Based on First-Life Historical Degradation Data Using Predictive Time-Series Algorithms
Next Article in Special Issue
Functional Solutions of Stochastic Differential Equations
Previous Article in Journal
HEDV-Greedy: An Advanced Algorithm for Influence Maximization in Hypergraphs
Previous Article in Special Issue
On the Bias of the Unbiased Expectation Theory
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Controlled Reflected McKean–Vlasov SDEs and Neumann Problem for Backward SPDEs

Department of Mathematics and Statistics, Hainan Normal University, Haikou 571158, China
*
Author to whom correspondence should be addressed.
Mathematics 2024, 12(7), 1050; https://doi.org/10.3390/math12071050
Submission received: 12 March 2024 / Revised: 26 March 2024 / Accepted: 29 March 2024 / Published: 31 March 2024
(This article belongs to the Special Issue First SDE: New Advances in Stochastic Differential Equations)

Abstract

:
This paper is concerned with the stochastic optimal control problem of a 1-dimensional McKean–Vlasov stochastic differential equation (SDE) with reflection, of which the drift coefficient and diffusion coefficient can be both dependent on the state of the solution process along with its law and control. One backward stochastic partial differential equation (BSPDE) with the Neumann boundary condition can represent the value function of this control problem. Existence and uniqueness of the solution to the above equation are obtained. Finally, the optimal feedback control can be constructed by the BSPDE.

1. Introduction

Let p , q R with p < q , and T > 0 . Consider the following controlled reflected McKean–Vlasov stochastic differential equation:
d X t = b t X t , L ( X t ) , α t d t + σ t X t , L ( X t ) , α t d B t + d L t d U t , t [ 0 , T ] ; X 0 = η ; L 0 = U 0 = 0 ; p X t q , a . s . ; 0 T ( X t p ) d L t = 0 T ( q X s ) d U s = 0 a . s . ,
where B is a 1-dimensional Brownian motion on a complete probability space ( Ω , F , P ) , F = ( F t ) t 0 is the natural filtration generated by B, augmented with an independent σ -algebra G F . The drift coefficient b and the diffusion coefficient σ are random. Let P ( R ) be the set of probability measures on R and L ( X t ) be the distribution of X t under P . Let η be a random variable satisfying E | η | 2 < . The control α valued in a Polish space is said to be an admissible control if it is a F -progressive process such that there is a unique solution to Equation (1). Denote by A the set of all admissible controls. We denote the Borel σ -algebra of A as B ( A ) . L and U are two nondecreasing processes. Through the paper, we aim to minimize a cost function of the control problem over A , whose form is
J t ( X t ; α ) = E [ t T f s X s , L ( X s ) , α s d s t T g s X s , L ( X s ) d L s + t T g s X s , L ( X s ) d U s + G X T , L ( X T ) | F t ] .
Define the value function of the reflected McKean–Vlasov control problem as
V t ( η ) : = ess inf α A J t ( X t ; α ) X t = η .
The optimal stochastic control problems for McKean–Vlasov SDEs have been researched by many authors. Pham and Wei [1,2] studied the stochastic control problem of McKean–Vlasov SDE, whose coefficients might be dependent on the joint distribution of the state of the solution process and its control, then they obtained a dynamic programming Bellman equation for a distribution-dependent optimal control problem. At last, they showed the uniqueness and viscosity property of the value function for the mean-field stochastic control problem by relying on a concept of lifted viscosity solutions to the Bellman equation. Bayraktar and Cosso [3] investigated an optimal stochastic control problem, of which the state process satisfied a McKean–Vlasov stochastic differential equation and the diffusion coefficient might be degenerate. They showed that the value function of the control problem had a nonlinear Feynman–Kac representation, then they strictly proved the dynamic programming principle by using this probabilistic representation. Hafayed and Meherrem [4] proved the sufficient and necessary conditions for optimal stochastic control of the systems driven by controlled McKean–Vlasov SDEs, where the coefficients are dependent on the state of the solution process along with its distribution and control.
The reflected optimal control problems were also investigated by many authors. Ferrari [5] studied the infinite time-horizon singular stochastic optimal control problems reflected at zero for a 1-dimensional diffusion and showed that there are qualitatively different kinds of optimal strategies dependent on the characteristics of the state-dependent instant reward. Baurdoux and Yamazaki [6] investigated the two-sided singular optimal control problems, offered a sufficient condition for an optimal double barrier scheme, and proved that the scheme is optimal when the running cost function is convex. Bayraktar and Qiu [7] studied the stochastic control problem for a 1-dimensional stochastic differential equation with reflection, characterized the value function of the control problem by using a BSPDE, and constructed the optimal control. Elliott et al. [8] investigated a Neumann problem for a BSPDE with the singular terminal condition; they showed that the equation can represent the value function of a constrained stochastic optimal control problem, and then they derived the optimal feedback control.
The stochastic optimal problem of reflected McKean–Vlasov equations has not been widely studied in this field. Some important models in finance, biology, cybernetics, and many other fields involve this problem, so the reflected McKean–Vlasov optimal control problem plays an essential role in practice. In this context, we consider that the drift coefficient and diffusion coefficient of the reflected McKean–Vlasov SDEs can be both dependent on the state of the solution process and its law and control. BSPDEs in the Wasserstein space of probability measures can be written in the Hilbert space by relying on the concept of lifted function. The functions G and g in the cost function depend on the state and its distribution. In addition, the result in [4] for the 1-dimensional case can be covered by our boundary reflected McKean–Vlasov optimal control problem. For future work, we will study the reflected McKean–Vlasov optimal control problem with jump. It is worth noting that [9] obtained the well-posedness of multidimensional reflected McKean–Vlasov SDEs.
The remaining organization of the paper is as follows. In Section 2, we provide some concepts and present the main assumptions and results. In Section 3, we study the existence and uniqueness of the strong solution to general nonlinear BSPDE for Neumann problems, on the basis of which we firstly obtain the priori estimates of strong solution to linear BSPDEs and then prove the existence and uniqueness of the strong solution to nonlinear BSPDEs by using a continuity method. Finally, in Section 4, we prove the main results.

2. Preliminary Knowledge and Main Result

2.1. Notations and Differentiability in Wasserstein Space

The set of probability measures ν on R is denoted by P 2 ( R ) , where ν 2 2 : = R | x | 2 ν ( d x ) < , that is, they are square integrable. For any ν P 2 ( R ) , the set of measurable functions φ : R R is denoted by L ν 2 ( R ) , where functions φ are square integrable with respect to ν . For any ν P 2 ( R ) , the set of measurable functions ψ : R × R R is denoted by L ν ν 2 ( R 2 ) , where functions ψ are square integrable with respect to the product measure ν ν . Define
ν ( φ ) : = R φ ( x ) ν ( d x ) , ν ν ( ψ ) : = R × R ψ ( x , x ) ν ( d x ) ν ( d x ) .
L ν ( R ) (respectively, L ν ν ( R × R ) ) is defined as the subset of elements φ L ν 2 ( R ) (respectively, L ν ν 2 ( R × R ) ), where elements φ are bounded ν -a.e. (respectively, ν ν -a.e.). We denote by · the essential supremum. The set of all square integrable random variables valued in R on ( Ω , G , P ) is denoted by L 2 ( G ; R ) . For any random variable X, its probability law is denoted by L ( X ) under P . Without any loss of generality, we assume that G is sufficiently rich to take R -valued random variables with square integrable law, that is, P 2 ( R ) = { L ( ξ ) : ξ L 2 ( G ; R ) } . We equip the measurable space P 2 ( R ) with the 2-Wasserstein distance
W 2 ( ν , ν ) : = inf R × R | x y | 2 π ( d x , d y ) 1 2 : π P 2 ( R × R ) with marginals ν a n d ν = inf E | ξ ξ | 2 1 2 : ξ , ξ L 2 ( G ; R ) with L ( ξ ) = ν , L ( ξ ) = ν ,
and endow the space P 2 ( R ) with the corresponding Borel σ -field B ( P 2 ( R ) ) .
We shall depend on the concept of derivative about probability measure, which is presented by P.L. Lions in [10]. The concept is established by lifting the functions u : P 2 ( R ) R into function u ˜ , which is defined on L 2 ( G ; R ) and satisfies u ˜ ( ξ ) = u ξ , L ( ξ ) . On the contrary, when function u ˜ is defined on L 2 ( G ; R ) , we call function u the inverse-lifted function of u ˜ , which is defined on P 2 ( R ) and satisfies u ( ξ , ν ) = u ˜ ( ξ ) for ν = L ( ξ ) . When the lifted function u ˜ is Fréchet differentiable (respectively, Fréchet differentiable with respect to continuous derivatives) on L 2 ( G ; R ) , then the function u is differentiable (respectively, C 1 ) on P 2 ( R ) . In terms of Riesz’ theorem: [ D u ˜ ] ( ξ ) ( Z ) = E [ D u ˜ ( ξ ) . Z ] , we use
D u ˜ ( ξ ) = ν u L ( ξ ) ( ξ ) ,
to represent the Fréchet derivative [ D u ˜ ] ( ξ ) in this case, which is regarded as an element D u ˜ ( ξ ) of L 2 ( G ; R ) , where derivative ν u L ( ξ ) : R R is called as the derivative of function u at ν = L ( ξ ) . Moreover, ν u ( ν ) L ν 2 ( R ) for ν P 2 ( R ) = { L ( ξ ) : ξ L 2 ( G ; R ) } . As mentioned in [11], if the function u is C 1 and for all ν P 2 ( R ) , a version of the mapping ξ R ν u ( ν ) ( ξ ) is continuous such that the mapping ( ν , ξ ) P 2 ( R ) × R ν u ( ν ) ( ξ ) is continuous at any point ( ν , ξ ) and ξ Supp ( ν ) , and if the mapping ξ R ν u ( ν ) ( ξ ) is differentiable, whose derivative is jointly continuous at any point ( ν , ξ ) and ξ Supp ( ν ) , then the function u is partially C 2 . We denote by x ν u ( ν ) ( ξ ) the gradient of ν u ( ν ) . If the function u is partially C 2 , x ν u ( ν ) L ν ( R ) , and for all compact set K of P 2 ( R ) ,
sup ν K R | ν u ( ν ) ( x ) | 2 ν ( d x ) + x ν u ( ν ) < ,
then we say that the function u C b 2 P 2 ( R ) . As mentioned in [11], when the lifted function u ˜ is twice continuously Fréchet differentiable on L 2 ( G ; R ) and has Lipschitz continuous Fréchet derivative, then we have u C b 2 P 2 ( R ) . In terms of Riesz’ theorem, the second Fréchet derivative D 2 u ˜ ( ξ ) in this case is viewed as a bilinear form or a symmetric operator (thus it is bounded) on L 2 ( G ; R ) , and we also obtain the equality ([12] Appendix A.2) that for all ξ L 2 ( G ; R ) , Z L 2 ( G ; R ) ,
E D 2 u ˜ ( ξ ) ( Z N ) . Z N = E x ν u L ( ξ ) ( ξ ) Z 2 ,
where N L 2 ( G ; R ) is independent of ( ξ , Z ) with zero mean and unit variance.
In order to guarantee the controlled reflected McKean–Vlasov, SDE (1) is well-posed, we state the following assumptions.
( H 1 )
(i) For all ( t , x , ν , α ) [ 0 , T ] × R × P 2 ( R ) × A , there is a constant κ such that | σ t ( x , ν , α ) | 2 κ > 0 .
(ii) The drift coefficient b and diffusion coefficient σ are measurable functions on [ 0 , T ] × R × P 2 ( R ) × A .
(iii) There is a positive constant C such that for any t [ 0 , T ] , x , x R , ν , ν P 2 ( R ) and α A ,
| b t ( x , ν , α ) b t ( x , ν , α ) | + | σ t ( x , ν , α ) σ t ( x , ν , α ) | C | x x | + W 2 ( ν , ν ) ,
| b t ( 0 , δ 0 , α ) | + | σ t ( 0 , δ 0 , α ) | C ,
where
δ 0 ( A ) = 1 , 0 A , 0 , else .
(iv) For all ( t , x , ν ) [ 0 , T ] × R × P 2 ( R ) , the functions α b t ( x , ν , α ) , σ t ( x , ν , α ) are continuous on A .
Under ( H 1 ), we see that b and σ meet the conditions in Theorem 3.2 of [13], and domain [ p , q ] satisfies the condition on domain D in Assumption 2.5 of [13], therefore Equation (1) has a unique solution, which is denoted by { X t α , 0 t T } for α A .
The items that appear in cost functional (2) are supposed to fulfill the following conditions.
( H 2 )
(i) The functions f , g , and G are measurable, respectively, on [ 0 , T ] × R × P 2 ( R ) × A , [ 0 , T ] × R × P 2 ( R ) and R × P 2 ( R ) .
(ii) There is a positive constant C such that for any ( t , x , ν , α ) [ 0 , T ] × R × P 2 ( R ) × A ,
| f t ( x , ν , α ) | + | g t ( x , ν ) | + | G ( x , ν ) | C 1 + | x | 2 + ν 2 2 .
(iii) The coefficients f , g , and G, respectively, on [ 0 , T ] × R × P 2 ( R ) × A , [ 0 , T ] × R × P 2 ( R ) , and R × P 2 ( R ) are continuous, and are also local Lipschitz continuous, uniformly with respect to A : there is a positive constant C such that for any t [ 0 , T ] , x , x R , ν , ν P 2 ( R ) and α A ,
| f t ( x , ν , α ) f t ( x , ν , α ) | + | g t ( x , ν ) g t ( x , ν ) | + | G ( x , ν ) G ( x , ν ) | C ( 1 + | x | + | x | + ν 2 + ν 2 ) | x x | + W 2 ( ν , ν ) .
Under ( H 2 ), for any ( t , ξ , α ) [ 0 , T ] × L 2 ( G ; R ) × A , the cost functional (2) is well-posed and finite.
Considering Peng’s pioneering research [14] about non-Markovian stochastic control, the dynamic programming principle related to the reflected McKean–Vlasov stochastic control indicates that the value function V in (3) is the first element of the pair ( v , ψ ) , which is the solution to the following BSPDE with Neumann boundary condition:
d v t L ( ξ ) = ess inf α A E [ ν v t L ( ξ ) ( ξ ) b t ξ , L ( ξ ) , α + 1 2 x ν v t L ( ξ ) ( ξ ) | σ t ξ , L ( ξ ) , α | 2 + f ˘ t L ( ξ ) , α ] d t ψ ˘ t L ( ξ ) d B t , t , L ( ξ ) [ 0 , T ] × P 2 ( [ p , q ] ) ; ν v t L ( ξ ) ( p ) = g ˘ t ( δ p ) , ν v t L ( ξ ) ( q ) = g ˘ t ( δ q ) ; v T L ( ξ ) = G ˘ L ( ξ ) , L ( ξ ) P 2 ( [ p , q ] ) ,
where d v t is the derivative of v with respect to t, and f ˘ , g ˘ , G ˘ , ψ ˘ is inverse-lifted function, respectively, of f ( · , L ( · ) ) , g ( · , L ( · ) ) , G ( · , L ( · ) ) , ψ .
Therefore, by the lifting identification, we regard the function v as a function on [ 0 , T ] × L 2 ( G ; R ) and reserve the same notation v t ( ξ ) = v t L ( ξ ) (notice that v is dependent on ξ just by its law), then we notice from the relation (4) and (5) between derivatives in the Wasserstein space P 2 ( R ) and in the Hilbert space L 2 ( G ; R ) that the BSPDE (6) is also written in [ 0 , T ] × L 2 ( G ; [ p , q ] ) as
d v t ( ξ ) = H t ξ , D v t ( ξ ) , D 2 v t ( ξ ) d t ψ t ( ξ ) d B t , ( t , ξ ) [ 0 , T ] × L 2 ( G ; [ p , q ] ) ; D v t ( p ) = g t ( p , δ p ) , D v t ( q ) = g t ( q , δ q ) ; v T ( ξ ) = G ξ , L ( ξ ) , ξ L 2 ( G ; [ p , q ] ) ,
with Hamiltonian function
H t ξ , D v t ( ξ ) , D 2 v t ( ξ ) = ess inf α A E { f t ξ , L ( ξ ) , α + D v t ( ξ ) b t ξ , L ( ξ ) , α + 1 2 D 2 v t ( ξ ) | σ t ξ , L ( ξ ) , α N | 2 } ,
for ( t , ξ ) [ 0 , T ] × L 2 ( G ; [ p , q ] ) ; here, N L 2 ( G ; [ p , q ] ) is independent of ξ .

2.2. Definition of Solutions to BSPDEs and Main Result

For a Banach space V , the set of G -measurable and square-integrable V -valued random variables is denoted by the space L 2 ( G ; V ) , and the set of P -measurable càdlàg V -valued process ( X t ) t [ 0 , T ] such that
X S G n ( 0 , T ; V ) n = E sup t [ 0 , T ] X t V n <
is denoted by S G n ( 0 , T ; V ) , where P is the predictable σ -algebra corresponding to { G t } t 0 on Ω × [ 0 , ) and n [ 1 , ) . The set of P -measurable V -valued processes ( v t ) t [ 0 , T ] satisfying
v L G n ( 0 , T ; V ) n = E 0 T v t V n d t < , n [ 1 , ) ;
v L G ( 0 , T ; V ) = sup ( ω , t ) Ω × [ 0 , T ] v t V < , n =
is denoted by L G n ( 0 , T ; V ) . For simplicity, we neglect the subscript for the space S G 2 ( 0 , T ; V ) and space L G 2 ( 0 , T ; V ) , particularly if there is no confusion about the adaptedness and filtration.
H m , n ( [ p , q ] ) is denoted as the Sobolev space of real-valued functions φ and its up to m-th order derivatives are in L n ( [ p , q ] ) , which is endowed with the Sobolev norm φ H m , n ( [ p , q ] ) , m N + and n [ 1 , ) . The space of trace-zero functions in H m , n ( [ p , q ] ) is denoted by H 0 m , n ( [ p , q ] ) . For m = 0 , H 0 , n ( [ p , q ] ) : = L n ( [ p , q ] ) . For simplicity, v 1 , . . . , v l H m , n ( [ p , q ] ) and v H m , n ( [ p , q ] ) n = j = 1 l v j H m , n ( [ p , q ] ) n can be expressed as v = ( v 1 , . . . , v l ) H m , n ( [ p , q ] ) . The norm and the inner product in the usual Hilbert space L 2 ( [ p , q ] ) are denoted, respectively, by · and · , · , and the duality between Hilbert space H m , 2 ( [ p , q ] ) and their dual spaces is denoted by · , · when there is no confusion. For m = 1 , 2 , we set
H = S G 2 0 , T ; L 2 ( [ p , q ] ) L G 2 0 , T ; H 1 , 2 ( [ p , q ] ) × L G 2 0 , T ; L 2 ( [ p , q ] ) ,
H m = S G 2 0 , T ; H m , 2 ( [ p , q ] ) L G 2 0 , T ; H m + 1 , 2 ( [ p , q ] ) × L G 2 0 , T ; H m , 2 ( [ p , q ] ) ,
and the two spaces are complete and equipped, respectively, with the norms
( v , ψ ) H 2 = v S G 2 ( 0 , T ; L 2 ( [ p , q ] ) ) 2 + v L G 2 ( 0 , T ; H 1 , 2 ( [ p , q ] ) ) 2 + ψ L G 2 ( 0 , T ; L 2 ( [ p , q ] ) ) 2 , ( v , ψ ) H ,
and
( v , ψ ) H m 2 = v S G 2 ( 0 , T ; H m , 2 ( [ p , q ] ) ) 2 + v L G 2 ( 0 , T ; H m + 1 , 2 ( [ p , q ] ) ) 2 + ψ L G 2 ( 0 , T ; H m , 2 ( [ p , q ] ) ) 2 , ( v , ψ ) H m .
Now, the notion of solutions to BSPDEs with general nonlinear coefficients is introduced as follows.
Definition 1. 
For any x , x 1 , x 2 R and any z , z 1 R , let R be a random function such that
R . ( · , x , x 1 , x 2 , z , z 1 ) : Ω × [ 0 , T ] × L 2 ( G ; [ p , q ] ) R
is P B ( [ p , q ] ) -measurable, and let G L 2 ( G ; L 2 ( [ p , q ] ) ) . We say that there is a weak solution ( v , ψ ) to the BSPDE:
d v t ( ξ ) = R t ( ξ , v , D v , D 2 v , ψ , D ψ ) d t ψ t ( ξ ) d B t , ( t , ξ ) [ 0 , T ] × L 2 ( G ; [ p , q ] ) ; D v t ( p ) = g t ( p , δ p ) , D v t ( q ) = g t ( q , δ q ) ; v T ( ξ ) = G ξ , L ( ξ ) , ξ L 2 ( G ; [ p , q ] ) ,
if ( v , ψ ) H and satisfies Equation (9) in the weak sense, that is, for all ϕ C c ( ( p , q ) ) ,
ϕ , R . ( · , v , D v , D 2 v , ψ , D ψ ) L G 1 ( 0 , T ; R )
and
ϕ , v t = ϕ , G + t T ϕ , R s ( · , v , D v , D 2 v , ψ , D ψ ) d s t T ϕ , ψ s d B s a . s . , t [ 0 , T ] .
If the regularity ( v , ψ ) H 1 has been proved, we call the above ( v , ψ ) a strong solution.
Particularly, we have a case of nonlinear term R with
R t ( ξ , v , D v , D 2 v , ψ , D ψ ) = H t ξ , D v t ( ξ ) , D 2 v t ( ξ ) ,
which is corresponding to BSPDE (7).
In order to obtain that Equation (7) is well-posed, we additionally take the assumptions below.
( H 3 )
(i) Let
G · , L ( · ) L 2 G ; H 2 , 2 ( [ p , q ] ) , D G · , L ( · ) g T · , L ( · ) L 2 G ; H 0 1 , 2 ( [ p , q ] ) ,
and along with another function Ψ , the pair ( g , Ψ ) belongs to H 1 and satisfies BSPDE
d g t ξ , L ( ξ ) = Q t ξ , L ( ξ ) d t Ψ t ξ , L ( ξ ) d B t , ξ L 2 ( G ; [ p , q ] ) ,
in the weak sense with Q L G 2 0 , T ; L 2 ( [ p , q ] ) .
(ii) For any u S G 2 0 , T ; H 1 , 2 ( [ p , q ] ) L G 2 0 , T ; H 2 , 2 ( [ p , q ] ) , we assume that H . ( · , u , D u ) , ( D H ) . ( · , u , D u ) L G 2 0 , T ; L 2 ( [ p , q ] ) , and there are nonnegative constants K 0 and K 1 such that for any u 1 , u 2 S G 2 0 , T ; H 1 , 2 ( [ p , q ] ) L G 2 0 , T ; H 2 , 2 ( [ p , q ] ) , ( t , ξ ) [ 0 , T ] × L 2 ( G , [ p , q ] ) ,
| H t ( ξ , u 1 , D u 1 ) H t ( ξ , u 2 , D u 2 ) | K 0 | u 1 u 2 | + K 1 | D ( u 1 u 2 ) | .
(iii) There is an A -valued P B ( R ) B ( R ) -measurable function Π satisfying
H t ( ξ , y , z ) = E [ f t ξ , L ( ξ ) , Π t ( ξ , y , z ) + b t ξ , L ( ξ ) , Π t ( ξ , y , z ) y + 1 2 | σ t ξ , L ( ξ ) , Π t ( ξ , y , z ) N | 2 z ] ,
that is,
Π t ( ξ , y , z ) arg ess inf α A E f t ξ , L ( ξ ) , α + b t ξ , L ( ξ ) , α y + 1 2 | σ t ξ , L ( ξ ) , α N | 2 z ,
and for any u S G 2 0 , T ; H 1 , 2 ( [ p , q ] ) L G 2 0 , T ; H 2 , 2 ( [ p , q ] ) , there exists a unique solution to the controlled reflected McKean–Vlasov SDE (1) associated with the drift coefficient b t X t , L ( X t ) , Π t ( X t , u t , D u t ) .
Finally, we summarize the main theorem. In the next section, we make some preparations to prove it.
Theorem 1. 
Suppose that ( H 1 ), ( H 2 ), and ( H 3 ) hold. There is a unique strong solution ( v , ψ ) to BSPDE (7). For the strong solution, we also obtain that ( v , ψ ) H 2 . In addition, for any η L 2 ( G ; [ p , q ] ) , E [ v t ( η ) ] = V t ( η ) , which is the value function (3). We obtain the optimal control α * and the associated state process X * , respectively, expressed as α * = Π t X t * , D v t ( X t * ) , D 2 v t ( X t * ) and
d X t * = b t X t * , L ( X t * ) , Π t X t * , D v t ( X t * ) , D 2 v t ( X t * ) d t + σ t X t * , L ( X t * ) , Π t X t * , D v t ( X t * ) , D 2 v t ( X t * ) d B t + d L t d U t ; X 0 * = ξ ; L 0 = U 0 = 0 ; p X t * q a . s . ; 0 T ( X t * p ) d L t = 0 T ( q X s * ) d U s = 0 a . s .
We take into consideration the conditions ( H 1 ) and ( H 2 ) as standing assumptions in this context, and they are standard to ensure that the BSPDE (7) is adapted and the controlled McKean–Vlasov SDE with reflection is well-posed.
By assumption ( H 3 ) (i), in order to obtain ( v , ψ ) H 2 , it is standard for the requests of G (by L p -theory of BSPDE of [15]); considering the Skorohod conditions of McKean–Vlasov SDE (1) with reflection, one has
0 T g s X s , L ( X s ) d L s = 0 T g s ( p , δ p ) d L s , 0 T g s X s , L ( X s ) d U s = 0 T g s ( q , δ q ) d U s ,
so the reflected control problem only involves g s ( p , δ p ) and g s ( q , δ q ) .
In assumption ( H 3 ) (ii), we suppose that the Hamilton function H t ( ξ , u , D u ) is Lipschitz continuous with respect to u and D u , which implies that u H . ( · , u , D u ) , D u H . ( · , u , D u ) L ( Ω × [ 0 , T ] × L 2 ( G , [ p , q ] ) for any u R .

3. Existence and Uniqueness of a Strong Solution to General Nonlinear BSPDEs

In this section, the existence and uniqueness of strong solution to general nonlinear BSPDEs with Neumann boundary condition is established. In this context, we only take into account the 1-dimensional case, although it is not essentially difficult with the multidimensional case.
Take into consideration the following Neumann problem for BSPDE:
d v t ( ξ ) = 1 2 | σ t ξ , L ( ξ ) , α | 2 D 2 v t ( ξ ) + Γ t ( ξ , v , D v , D 2 v , ψ , D ψ ) d t ψ t ( ξ ) d B t , D v t ( p ) = 0 , D v t ( q ) = 0 ; v T ( ξ ) = G ξ , L ( ξ ) , ξ L 2 ( G , [ p , q ] ) .
In this section, we take the following hypothesis.
( H 4 )
For any ( v , ψ ) H 2 , 2 ( [ p , q ] ) × H 1 , 2 ( [ p , q ] ) , Γ . ( · , v , D v , D 2 v , ψ , D ψ ) L 2 ( 0 , T ; L 2 ( [ p , q ] ) ) , and there are nonnegative constants L 1 and L 2 such that for any ( v i , ψ i ) H 2 , 2 ( [ p , q ] ) × H 1 , 2 ( [ p , q ] ) , i = 1 , 2 , for any t [ 0 , T ] , there holds
Γ t ( · , v 1 , D v 1 , D 2 v 1 , ψ 1 , D ψ 1 ) Γ t ( · , v 2 , D v 2 , D 2 v 2 , ψ 2 , D ψ 2 ) L 1 D 2 ( v 1 v 2 ) + D ( ψ 1 ψ 2 ) + L 2 v 1 v 2 H 1 , 2 ( [ p , q ] ) + ψ 1 ψ 2 L 2 ( [ p , q ] ) a . s . .
Theorem 2. 
Let G L 2 ( G ; H 1 , 2 ( [ p , q ] ) ) and assume that ( H 2 ) and ( H 4 ) hold. There is a positive constant L 0 depending on L 2 and T, such that when 0 L 1 < L 0 , there is a unique strong solution ( v , ψ ) to BSPDE (12) which satisfies
( v , ψ ) H 1 C G L 2 ( G ; H 1 , 2 ( [ p , q ] ) ) + Γ 0 L 2 ( 0 , T ; L 2 ( [ p , q ] ) ) ,
in which Γ 0 : = Γ . ( · , 0 , 0 , 0 , 0 , 0 ) and C is dependent on L 1 , L 2 , κ, and T.
Firstly, we shall establish some priori estimates for linear equations in Section 3.1 and then exploit the method of continuity to finish the proof of Theorem 2 in Section 3.2.

3.1. The Priori Estimates

For every γ [ 0 , 1 ] , we take into account the following linear BSPDE:
d v t ( ξ ) = γ 2 | σ t ξ , L ( ξ ) , α | 2 D 2 v t ( ξ ) + 1 γ 2 D 2 v t ( ξ ) + F t ( ξ ) d t ψ t ( ξ ) d B t , ( t , ξ ) [ 0 , T ] × L 2 ( G , [ p , q ] ) ; D v t ( p ) = 0 , D v t ( q ) = 0 ; v T ( ξ ) = G ξ , L ( ξ ) , ξ L 2 ( G , [ p , q ] ) .
We place the proof of the following Proposition 1 in Appendix A.
Proposition 1. 
Assume that ( H 2 ) holds and
F L 2 0 , T ; L 2 ( [ p , q ] ) , G L 2 G ; H 1 , 2 ( [ p , q ] ) .
Assume that ( v , ψ ) is a strong solution to BSPDE for a Neumann problem (14). Then, this strong solution is unique, satisfying
( v , ψ ) H 2 C 1 G L 2 ( G ; L 2 ( [ p , q ] ) ) 2 + E 0 T | F s , v s | d s ,
and
( v , ψ ) H 1 2 C 2 G L 2 ( G ; H 1 , 2 ( [ p , q ] ) ) 2 + E 0 T | F s , v s | + | F s , D 2 v s | d s C 3 G L 2 ( G ; H 1 , 2 ( [ p , q ] ) ) 2 + F L 2 ( 0 , T ; L 2 ( [ p , q ] ) ) 2 ,
where the constants C 1 , C 2 , and C 3 are only dependent on κ and T, and independent of γ [ 0 , 1 ] .

3.2. Existence and Uniqueness of Strong Solution to (12) (Proof of Theorem 2)

Firstly, we take into account the BSPDE for Neumann problem with Laplacian operator as follows:
d v t ( ξ ) = D 2 v t ( ξ ) + F t ( ξ ) d t ψ t ( ξ ) d B t , ( t , ξ ) [ 0 , T ] × L 2 ( G , [ p , q ] ) ; D v t ( p ) = 0 , D v t ( q ) = 0 ; v T ( ξ ) = G ( ξ , L ( ξ ) ) , ξ L 2 ( G , [ p , q ] ) .
Proposition 2. 
Set
F L G 2 0 , T ; L 2 ( [ p , q ] ) , G L 2 G ; H 1 , 2 ( [ p , q ] ) .
Then there is a unique strong solution ( v , ψ ) to BSPDE (15).
Proof. 
In view of Proposition 1, we can directly obtain the uniqueness of the strong solution. Hence, it just requires us to prove the existence of the strong solution.
Step 1. Further assume that F L G 2 0 , T ; H 1 , 2 ( [ p , q ] ) and D G L 2 G ; H 0 1 , 2 ( [ p , q ] ) . In view of the theory about the deterministic parabolic PDEs for the Neumann problem of Theorem 7.22 in [16], there is a unique strong solution v ^ to PDE:
t v ^ t ( ξ ) = D 2 v ^ t ( ξ ) + F t ( ξ ) , ( t , ξ ) [ 0 , T ] × L 2 ( G , [ p , q ] ) ; D v ^ t ( p ) = 0 , D v ^ t ( q ) = 0 ; v ^ T ( ξ ) = G ( ξ , L ( ξ ) ) , ξ L 2 ( G , [ p , q ] ) ,
such that v ^ , D v ^ , D 2 v ^ , t v ^ L 2 G ; L 2 ( [ 0 , T ] × [ p , q ] ) . We first take conditional expectations in Hilbert spaces ([17]) and then set
v t = E [ v ^ t | G t ] a . s . , for each t [ 0 , T ] ,
which admits a version in S G 2 0 , T ; L 2 ( [ p , q ] ) L G 2 0 , T ; H 2 , 2 ( [ p , q ] ) . Along with ψ L G 2 0 , T ; L 2 ( [ p , q ] ) , v satisfies L 2 ( [ p , q ] ) -valued BSDE:
d v t ( ξ ) = D 2 v t ( ξ ) + F t ( ξ ) d t ψ t ( ξ ) d B t ; v T ( ξ ) = G ( ξ , L ( ξ ) ) , ξ L 2 ( G , [ p , q ] ) .
Considering that the definition of v, v satisfies the zero-Neumann boundary conditions, in view of Definition 1 and Proposition 1, we may easily obtain that ( v , ψ ) is a weak solution to Equation (15).
Step 2. Now let us show that the weak solution ( v , ψ ) constructed above is indeed a unique strong solution to Equation (15). Using a similar method for Step 1, we can easily check that D v ^ is the strong solution of the Dirichlet problem:
t u ^ t ( ξ ) = D 2 u ^ t ( ξ ) + D F t ( ξ ) , ( t , ξ ) [ 0 , T ] × L 2 ( G , [ p , q ] ) ; u ^ t ( p ) = 0 , u ^ t ( q ) = 0 ; u ^ T ( ξ ) = D G ξ , L ( ξ ) , ξ L 2 ( G , [ p , q ] ) ,
and ( D v , D ψ ) satisfies L 2 ( [ p , q ] ) -valued BSDE (16) corresponding to the coefficients ( D F , D G ) . Particularly, we obtain
( D v , D ψ ) S G 2 0 , T ; L 2 ( [ p , q ] ) × L G 2 0 , T ; L 2 ( [ p , q ] ) ,
and, hence, ( v , ψ ) is a strong solution to Equation (15). For F L G 2 0 , T ; L 2 ( [ p , q ] ) and G L 2 G ; H 1 , 2 ( [ p , q ] ) , we can choose a sequence
{ ( F n , G n ) } n N + L G 2 0 , T ; H 1 , 2 ( [ p , q ] ) × L 2 G ; H 1 , 2 ( [ p , q ] )
with { D G n } n N + L 2 G ; H 0 1 , 2 ( [ p , q ] ) such that ( F n , G n ) converges to ( F , G ) in L G 2 ( 0 , T ;
L 2 ( [ p , q ] ) ) × L 2 G ; H 1 , 2 ( [ p , q ] ) . For each ( F n , G n ) , we obtain the strong solution ( v n , ψ n ) to BSPDE (15) associated with the coefficients ( F n , G n ) . Therefore, we can obtain the convergence of ( v n , ψ n ) and the existence of strong solution from the estimates in Proposition 1. □
Finally, we can exploit the method of continuity to complete the proof of Theorem 2.
Proof. 
Step 1. For each γ [ 0 , 1 ] and F L G 2 0 , T ; L 2 ( [ p , q ] ) , consider the following BSPDE:
d v t ( ξ ) = [ γ 2 | σ t ξ , L ( ξ ) , α | 2 D 2 v t ( ξ ) + γ Γ t ( ξ , v , D v , D 2 v , ψ , D ψ ) + ( 1 γ ) D 2 v t ( ξ ) + F t ( ξ ) ] d t ψ t ( ξ ) d B t , ( t , ξ ) [ 0 , T ] × L 2 ( G , [ p , q ] ) ; D v t ( p ) = 0 , D v t ( q ) = 0 ; v T ( ξ ) = G ξ , L ( ξ ) , ξ L 2 ( G , [ p , q ] ) .
When γ = 1 and F = 0 , we observe that the BSPDE (12) is a special case of the above equation. Then we shall extend priori estimates inferred from the linear case to the nonlinear BSPDE (17). Assume that ( v , ψ ) is a strong solution to Equation (17). For any t [ 0 , T ] , we apply Proposition 1 (see also estimates (A4) and (A7)) and then we obtain, by ( H 4 ),
E sup s [ t , T ] v s 2 + E t T ψ s 2 + D v s 2 d s C E G 2 + t T | F s + γ Γ s ( · , v , D v , D 2 v , ψ , D ψ ) , v s | d s C E G 2 + ε t T F s 2 + Γ s ( · , v , D v , D 2 v , ψ , D ψ ) 2 d s + 2 ε t T v s 2 d s C 1 E [ G 2 + ε t T F s 2 + Γ s 0 2 + v s 2 + D v s 2 + ψ s 2 d s + 2 ε t T v s 2 d s + ε L 1 2 t T D ψ s 2 + D 2 v s 2 d s ] ,
and
E sup s [ t , T ] v s H 1 , 2 ( [ p , q ] ) 2 + E t T v s H 2 , 2 ( [ p , q ] ) 2 + ψ s H 1 , 2 ( [ p , q ] ) 2 d s C E [ G H 1 , 2 ( [ p , q ] ) 2 + t T | F s + γ Γ s ( · , v , D v , D 2 v , ψ , D ψ ) , v s | d s + t T | F s + γ Γ s ( · , v , D v , D 2 v , ψ , D ψ ) , D 2 v s | d s ] C E [ G H 1 , 2 ( [ p , q ] ) 2 + 1 + 1 ε t T F s 2 + Γ s 0 2 + v s 2 + D v s 2 + ψ s 2 d s + t T L 1 2 ε + ε D 2 v s 2 + L 1 2 ε D ψ s 2 d s ] .
Letting 0 < ε < 1 C 1 , by (18) we have
E sup s [ t , T ] v s 2 + E t T ψ s 2 + D v s 2 d s C 2 E G 2 + t T F s 2 + Γ s 0 2 + 2 ε v s 2 + ε L 1 2 D ψ s 2 + D 2 v s 2 d s ,
with C 2 independent of ( ε , L 1 ) . From (20) and (19), there is a constant L 0 depending on L 2 and T such that as L 1 < L 0 , setting ε sufficiently small and applying Grownwall inequality, we have
( v , ψ ) H 1 C G L 2 ( G ; H 1 , 2 ( [ p , q ] ) ) + Γ 0 L ( 0 , T ; L 2 ( [ p , q ] ) ) + F L ( 0 , T ; L 2 ( [ p , q ] ) ) ,
with the constant C depending on κ , L 1 , L 2 , and T.
Step 2. Suppose that ( v 1 , ψ 1 ) and ( v 2 , ψ 2 ) are two strong solutions to Equation (17). Then the pair ( δ v , δ ψ ) = ( v 1 v 2 , ψ 1 ψ 2 ) satisfies the following BSPDE:
d δ v t ( ξ ) = { γ [ 1 2 | σ t ξ , L ( ξ ) , α | 2 D 2 δ v t ( ξ ) + Γ t ( ξ , v 1 , D v 1 , D 2 v 1 , ψ 1 , D ψ 1 ) Γ t ( ξ , v 2 , D v 2 , D 2 v 2 , ψ 2 , D ψ 2 ) ] + ( 1 γ ) D 2 δ v t ( ξ ) } d t δ ψ t ( ξ ) d B t ; D δ v t ( p ) = 0 , D δ v t ( q ) = 0 ; δ v T ( ξ ) = 0 .
Recalling ( H 4 ), we obtain
Γ t ( · , v 1 , D v 1 , D 2 v 1 , ψ 1 , D ψ 1 ) Γ t ( · , v 2 , D v 2 , D 2 v 2 , ψ 2 , D ψ 2 ) L 1 D 2 ( v 1 v 2 ) + D ( ψ 1 ψ 2 ) + L 2 v 1 v 2 H 1 , 2 ( [ p , q ] ) + ψ 1 ψ 2 L 2 ( [ p , q ] ) ,
a.s. for each t [ 0 , T ] . Using a similar method for Step 1 and Itô’s formula to the square norm of ( δ v , δ ψ ) , we can obtain the estimates (18)–(20), in addition to (21) where ( G , Γ 0 , F ) is substituted with zero values. This suggests that the strong solution, respectively, to Equation (17) and (12) is unique.
Step 3. Firstly, we notice that the priori estimate (21) holds in which the constant C is independent of γ [ 0 , 1 ] . Proposition 2 implies that Equation (17) has a unique strong solution ( v , ψ ) when γ = 0 . Noticing that BSPDE (12) corresponds to the special case of BSPDE (17) when γ = 1 , then we want to generalize the well-posedness of BSPDE through the interval [ 0 , 1 ] starting from γ = 0 .
Suppose that for some γ = γ 0 , BSPDE (17) which satisfies assumptions ( H 1 ), ( H 2 ), and ( H 4 ) has a unique strong solution ( v , ψ ) , which holds when γ 0 = 0 (by Proposition 2 as above). Then, for each ( v ˇ , ψ ˇ ) H 1 , the following BSPDE:
d v t ( ξ ) = { γ 0 1 2 | σ t ξ , L ( ξ ) , α | 2 D 2 v t ( ξ ) + Γ t ( ξ , v , D v , D 2 v , ψ , D ψ ) + ( γ γ 0 ) 1 2 | σ t ξ , L ( ξ ) , α | 2 D 2 v ˇ t ( ξ ) + Γ t ( ξ , v ˇ , D v ˇ , D 2 v ˇ , ψ ˇ , D ψ ˇ ) D 2 v ˇ t ( ξ ) + ( 1 γ 0 ) D 2 v t ( ξ ) } d t ψ t ( ξ ) d B t , ( t , ξ ) [ 0 , T ] × L 2 ( G ; [ p , q ] ) ; D v t ( p ) = 0 , D v t ( q ) = 0 ; v T ( ξ ) = G ξ , L ( ξ ) , ξ L 2 ( G ; [ p , q ] ) ,
corresponds to a special case of BSPDE (17) when γ = γ 0 and
F t ( ξ ) = ( γ γ 0 ) 1 2 | σ t ξ , L ( ξ ) , α | 2 D 2 v ˇ t ( ξ ) + Γ t ( ξ , v ˇ , D v ˇ , D 2 v ˇ , ψ ˇ , D ψ ˇ ) D 2 v ˇ t ( ξ ) ,
and it admits a unique strong solution ( v , ψ ) ; we can further define the following solution map:
M γ 0 : H 1 H 1 , ( v ˇ , ψ ˇ ) ( v , ψ ) .
For each ( v i , ψ i ) H 1 , i = 1 , 2 , using a similar method for Step 2, we obtain
( v 1 v 2 , ψ 1 ψ 2 ) H 1 C | γ γ 0 | 1 2 | σ t ξ , L ( ξ ) , α | 2 D 2 ( v ˇ 1 v ˇ 2 ) + Γ . ( · , v ˇ 1 , D v ˇ 1 , D 2 v ˇ 1 , ψ ˇ 1 , D ψ ˇ 1 ) Γ . ( · , v ˇ 2 , D v ˇ 2 , D 2 v ˇ 2 , ψ ˇ 2 , D ψ ˇ 2 ) D 2 ( v ˇ 1 v ˇ 2 ) L G 2 ( 0 , T ; L 2 ( [ p , q ] ) ) C ˜ | γ γ 0 | ( v ˇ 1 v ˇ 2 , ψ ˇ 1 ψ ˇ 2 ) H 1 ,
with the constant C ˜ being independent of ( γ , γ 0 ) . When | γ γ 0 | < 1 C ˜ , M γ 0 is a contraction mapping admitting a unique fixed point ( v , ψ ) H 1 , and it is a strong solution to Equation (17). With this method, when Equation (17) admits a strong solution for γ 0 , then it also satisfies | γ γ 0 | < 1 / C ˜ for any γ . We can arrive at γ = 1 in finite steps starting from γ = 0 . Applying the estimates (21) and the uniqueness obtained in Step 2, we can complete the proof. □

4. Proof of Theorem 1

Lemma 1 
([7] Itô-Kunita–Wentzell formula). Let
X t = x + 0 t η r d r + 0 t d A r + 0 t ρ r d B r , 0 t T ,
where ( η , ρ ) L G 2 ( 0 , T ; R ) × L G 2 ( 0 , T ; R ) and A is a continuous G -adapted bounded variation process with A 0 = 0 . Assume that almost surely p X t q and
v t ( x ) = v 0 ( x ) + 0 t μ r ( x ) d r + 0 t ψ r ( x ) d W r , f o r a l l ( t , x ) [ 0 , T ] × [ p , q ] ,
holds in the weak sense with ( v , μ , ψ ) belonging to
S 2 ( 0 , T ; H 2 , 2 ( [ p , q ] ) ) L 2 ( 0 , T ; H 3 , 2 ( [ p , q ] ) ) × L 2 ( 0 , T ; H 1 , 2 ( [ p , q ] ) ) × L 2 ( 0 , T ; H 2 , 2 ( [ p , q ] ) ) .
Thus, for any x [ p , q ] and t [ 0 , T ] ,
v t ( X t ) v 0 ( x ) = 0 t μ r ( X r ) + η r D v r ( X r ) + 1 2 | ρ | 2 D 2 v r ( X r ) d r + 0 t D v r ( X r ) d A r + 0 t ψ r ( X r ) + D v r ( X r ) ρ r d B r . a . s .
Applying Sobolev’s embedding theorem, H m , 2 ( R ) is embedded into continuous function space C m 1 continuously. Therefore, for any x [ p , q ] , Equation (23) holds. With this method, Lemma 1 is similar to the first expression of Kunita [18], pp. 118–119, when the whole line R is replaced by the bounded domain [ p , q ] . For brevity, we just take into account the 1-dimensional case in Lemma 1, and it is not essentially difficult to generalize it to the multidimensional case.
Now, we present a consequence of the Dirichlet problem of BSPDE as below, whose proof is the same as that of [19], Theorem 2, under our assumptions.
Lemma 2. 
Regard the Dirichlet problem as follows:
d v t ( ξ ) = 1 2 | σ t ξ , L ( ξ ) , α | 2 D 2 v t ( ξ ) + F t ( ξ ) d t ψ t ( ξ ) d B t , ( t , ξ ) [ 0 , T ] × L 2 ( G ; [ p , q ] ) ; v t ( p ) = 0 , v t ( q ) = 0 ; v T ( ξ ) = G ξ , L ( ξ ) , ξ L 2 ( G ; [ p , q ] ) .
with G L 2 ( G ; H 0 1 , 2 ( [ p , q ] ) ) and F L 2 ( 0 , T ; L 2 ( [ p , q ] ) ) . Under assumption ( H 2 ), there is a unique weak solution ( v , ψ ) to BSPDE (24) and it is unique with
( v , ψ ) H 1 C G L 2 ( G ; H 1 , 2 ( [ p , q ] ) ) + F L 2 ( 0 , T ; L 2 ( [ p , q ] ) ) ,
where the constant C depends on κ and T.
Finally, we begin to accomplish the proof of Theorem 1.
Proof. 
Firstly, the Neumann problem (7) can be reduced to the instance of zero Neumann boundary condition. By assumption ( H 3 ) (i) as well as Definition 1, we set
( g ^ t , Q ^ t , Ψ ^ t ) ( ξ ) = 0 ξ ( g t , Q t , Ψ t ) y , L ( y ) d y , for each ( t , ξ ) [ 0 , T ] × L 2 ( G ; [ p , q ] ) ,
and then we obtain ( g ^ , Ψ ^ ) as the strong solution to the following BSPDE:
d g ^ t ( ξ ) = 1 2 | σ t ξ , L ( ξ ) , α | 2 D 2 g ^ t ( ξ ) + f ^ t ( ξ ) d t Ψ ^ t ( ξ ) d B t , ( t , ξ ) [ 0 , T ] × L 2 ( G ; [ p , q ] ) ; D g ^ t ( p ) = g t ( p , δ p ) , D g ^ t ( q ) = g t ( q , δ q ) ; g ^ T ( ξ ) = 0 ξ g T y , L ( y ) d y , ξ L 2 ( G ; [ p , q ] ) ,
with f ^ L 2 ( 0 , T ; H 1 , 2 ( [ p , q ] ) ) being defined as f ^ t ( ξ ) = 1 / 2 | σ t ξ , L ( ξ ) , α | 2 D 2 g ^ t ( ξ ) + Q ^ t ( ξ ) . Hence, the existence and uniqueness of the strong solution ( v , ψ ) to Equation (7) is equal to that of strong solution ( v ˜ , ψ ˜ ) = ( v g ^ , ψ Ψ ^ ) satisfying the following BSPDE:
d v ˜ t ( ξ ) = [ 1 2 | σ t ξ , L ( ξ ) , α | 2 D 2 v ˜ t ( ξ ) + H t ξ , D v ˜ t ( ξ ) + D g ^ t ( ξ ) , D 2 v ˜ t ( ξ ) + D 2 g ^ t ( ξ ) f ^ t ( ξ ) ] d t ψ ˜ t ( ξ ) d B t , ( t , ξ ) [ 0 , T ] × L 2 ( G ; [ p , q ] ) ; D v ˜ t ( p ) = 0 , D v ˜ t ( q ) = 0 ; v ˜ T ( ξ ) = G ξ , L ( ξ ) g ^ T ( ξ ) , ξ L 2 ( G ; [ p , q ] ) .
In view of Theorem 2, there is a unique strong solution ( v ˜ , ψ ˜ ) to BSPDE (27). Taking derivatives of the above equation, we can easily check that
( u , ζ ) : = ( D v ˜ , D ψ ˜ ) = ( D v D g ^ , D ψ D Ψ ^ ) = ( D v g , D ψ Ψ )
is a weak solution to the Dirichlet problem as follows:
d u t ( ξ ) = 1 2 | σ t ξ , L ( ξ ) , α | 2 D 2 u t ( ξ ) + R t ( ξ , α ) d t ζ t ( ξ ) d B t , ( t , ξ ) [ 0 , T ] × L 2 ( G ; [ p , q ] ) ; u t ( p ) = 0 , u t ( q ) = 0 ; u T ( ξ ) = D G ξ , L ( ξ ) g T ξ , L ( ξ ) , ξ L 2 ( G ; [ p , q ] ) ,
with
R t ( ξ , α ) = D f ^ t ( ξ ) + ( ξ H ) t ξ , D v ˜ t ( ξ ) + g t ξ , L ( ξ ) , D 2 v ˜ t ( ξ ) + D g t ξ , L ( ξ ) + ( u 1 H ) t ξ , D v ˜ t ( ξ ) + g t ξ , L ( ξ ) , D 2 v ˜ t ( ξ ) + D g t ξ , L ( ξ ) D 2 v ˜ t ( ξ ) + D g t ξ , L ( ξ ) + ( u 2 H ) t ξ , D v ˜ t ( ξ ) + g t ξ , L ( ξ ) , D 2 v ˜ t ( ξ ) + D g t ξ , L ( ξ ) D 3 v ˜ t ( ξ ) + D 2 g t ξ , L ( ξ ) = D f ^ t ( ξ ) + ess inf α A E { D f t ξ , L ( ξ ) , α + D v ˜ t ( ξ ) + g t ξ , L ( ξ ) D b t ξ , L ( ξ ) , α + D 2 v ˜ t ( ξ ) + D g t ξ , L ( ξ ) b t ξ , L ( ξ ) , α + 1 2 D 3 v ˜ t ξ , L ( ξ ) + D 2 g t ξ , L ( ξ ) · | σ t ξ , L ( ξ ) , α N | 2 + D 2 v ˜ t ( ξ ) + D g t ξ , L ( ξ ) σ t ξ , L ( ξ ) , α N 2 D σ t ξ , L ( ξ ) , α } .
As ( v ˜ , ψ ˜ ) is a unique strong solution to BSPDE (27), we arrive at
v ˜ S G 2 0 , T ; H 1 , 2 ( [ p , q ] ) L G 2 0 , T ; H 2 , 2 ( [ p , q ] ) L 2 ( 0 , T ; L 2 ( [ p , q ] ) ) .
In view of assumptions ( H 1 ) and ( H 2 ), we have D f , b , D b , σ , D σ L 2 ( 0 , T ; L 2 ( [ p , q ] ) ) . By assumption (i) of ( H 3 ), we obtain g , D g , D 2 g L 2 ( 0 , T ; L 2 ( [ p , q ] ) ) , and D f ^ L 2 ( 0 , T ; L 2 ( [ p , q ] ) ) ; therefore, function R L 2 ( 0 , T ; L 2 ( [ p , q ] ) ) . Then, we may derive from Lemma 2 that ( u , ζ ) is a strong solution. Hence ( D v ˜ , D ψ ˜ ) = ( u , ζ ) H 1 ; furthermore, ( D v , D ψ ) = ( D v ˜ + g , D ψ ˜ + Ψ ) H 1 . Therefore, ( v , ψ ) H 2 . The regularity and assumption (iii) of ( H 3 ) imply that α t * = Π t ( X t * , D v t ( X t * ) , D 2 v t ( X t * ) ) is admissible. For each admissible control α , exploiting the general Itô-Kunita–Wentzell formula to v t ( X t α ) implies that for any η L 2 ( G ; [ p , q ] ) , there almost surely holds
E v t ( X t α ) v T ( X T α ) = E t T H r X r α , D v r ( X r α ) , D 2 v r ( X r α ) + D v r ( X r α ) b r X r α , L ( X r α ) , α r d r 1 2 E t T D 2 v r ( X r α ) | σ r X r α , L ( X r α ) , α r | 2 d r E t T D v r ( X r α ) d ( L r U r ) = E t T ess inf α ˜ A { D v r ( X r α ) b r X r α , L ( X r α ) , α ˜ r + 1 2 D 2 v r ( X r α ) | σ r X r α , L ( X r α ) , α ˜ r N | 2 + f r X r α , L ( X r α ) , α ˜ r } d r E t T [ D v r ( X r α ) b r X r α , L ( X r α ) , α r + 1 2 D 2 v r ( X r α ) · | σ r X r α , L ( X r α ) , α r | 2 ] d r E t T D v r ( X r α ) d ( L r U r ) E t T f r X r α , L ( X r α ) , α r d r + E t T D v r ( X r α ) d ( U r L r ) , t [ 0 , T ] ,
that is, for any t [ 0 , T ] ,
E v t ( X t α ) E G X T α , L ( X T α ) + t T f r X r α , L ( X r α ) , α r d r + t T g X r α , L ( X r α ) d ( U r L r ) .
Hence, for each admissible control α , it holds almost surely:
E v t ( X t α ) J t ( X t α ; α ) , t [ 0 , T ] .
Moreover, using a similar method for inequality (29), for all η L 2 ( G ; [ p , q ] ) and t [ 0 , T ] ,
E v t ( X t * ) = E G X T * , L ( X T * ) + t T f r X r * , L ( X r * ) , α r * d r + t T g X r * , L ( X r * ) d ( U r L r ) = J t ( X t * ; α * )
holds with probability 1. Therefore, by (30), for any η L 2 ( G ; [ p , q ] ) , E [ v t ( η ) ] = V t ( η ) and optimal stochastic control is given by α t * = Π t ( X t * , D v t ( X t * ) , D 2 v t ( X t * ) ) with the associated state process X t * , which satisfies BSPDE (11). The proof is completed. □

5. Conclusions

In this paper, we investigated the stochastic optimal control problem of a 1-dimensional McKean–Vlasov stochastic differential equation reflected at [ p , q ] , whose drift and diffusion coefficients are both dependent on the state of the solution process along with its law and control; thus, the results of [4,7] were generalized. We depend on the concept of lifted function presented by P.L. Lions in [10], then the BSPDEs in the Wasserstain space of probability measures can be written in the Hilbert space. Thus, the value function of our control problem can be represented by a BSPDE with Neumann boundary condition in the Hilbert space. The existence and uniqueness of the solution to the above BSPDE were obtained. Finally, we constructed the optimal feedback control of the stochastic control problem by the Itô-Kunita–Wentzell formula. For the future work, we will investigate the multidimensional reflected McKean–Vlasov optimal control problem with jump and singular terminal conditions, and the corresponding applications to population biology, finance, and physics will also be studied.

Author Contributions

Formal analysis, L.M., F.S. and X.H.; investigation, L.M., F.S. and X.H.; methodology, L.M., F.S. and X.H.; project administration, X.H.; software, X.H.; Validation, L.M., F.S. and X.H.; Writing—original draft, L.M. and F.S.; and writing—review and editing, L.M., F.S. and X.H. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by National Natural Science Foundation of China grant number 11861029, Hainan Provincial Natural Science Foundation of China grants number 122MS056 and 124MS056, and the special research fund of The Innovation Platform for Academicians of Hainan Province grant number YSPTZX202215.

Data Availability Statement

Data are contained within the article.

Conflicts of Interest

The authors declare no conflicts of interest.

Appendix A. Poof of Proposition 1

Proof. 
Step 1. Apply the Itô’s formula of Lemma A.1 in [7] to the square norm, which arises
v t 2 + t T ψ s 2 d s G 2 = t T v s , γ | σ s ξ , L ( ξ ) , α | 2 D 2 v s + ( 1 γ ) D 2 v s + 2 F s d s 2 t T v s , ψ s d B s a . s .
for any t [ 0 , T ] . Considering Neumann boundary conditions, we obtain
t T v s , ( 1 γ ) D 2 v s d s = ( 1 γ ) t T D v s 2 d s ,
and
t T v s , | σ s ξ , L ( ξ ) , α | 2 D 2 v s d s = t T D v s , | σ s ξ , L ( ξ ) , α | 2 D v s d s t T v s D | σ s ξ , L ( ξ ) , α | 2 , D v s d s t T D v s , | σ s ξ , L ( ξ ) , α | 2 D v s d s + 2 C κ t T v s 2 d s + κ 2 t T D v s 2 d s κ 2 t T D v s 2 d s + 2 C κ t T v s 2 d s .
Hence,
v t 2 + t T ψ s 2 d s G 2 = t T v s , γ | σ s ξ , L ( ξ ) , α | 2 D 2 v s d s ( 1 γ ) t T D v s 2 d s + t T v s , 2 F s d s 2 t T v s , ψ s d B s a . s .
for all t [ 0 , T ] . Notice that γ [ 0 , 1 ] and
κ γ 2 + ( 1 γ ) 1 κ 2 > 0 .
Next, we claim
E sup s [ t , T ] v s 2 + E t T ψ s 2 + D v s 2 d s C E G 2 + t T | F s , v s | d s .
In order to prove (A4), we estimate
δ E sup s [ t , T ] v s 2 + ( 1 δ ) E v t 2 + E t T ψ s 2 + D v s 2 d s
for δ = 0 and δ = 1 separately.
Firstly, for the case δ = 0 , by (A1) and (A2), we have
E v t 2 + t T ψ s 2 d s = E G 2 ( 1 γ ) t T D v s 2 d s + t T v s , γ | σ s ξ , L ( ξ ) , α | 2 D 2 v s + 2 F s d s E [ G 2 ( 1 γ ) t T D v s 2 d s γ κ 2 t T D v s 2 d s + 2 γ C κ t T v s 2 d s + 2 t T | v s , F s | d s ] ,
thus, by (A3), we have
E v t 2 + E t T ψ s 2 d s + 1 κ 2 E t T D v s 2 d s E v t 2 + E t T ψ s 2 d s + κ γ 2 + 1 γ E t T D v s 2 d s E G 2 + 2 t T | v s , F s | d s + 2 γ C κ t T v s 2 d s .
Applying the Gronwall inequality, we obtain
E v t 2 + E t T ψ s 2 d s + E t T D v s 2 d s C E G 2 + t T | v s , F s | d s .
Secondly, when δ = 1 , by (A1) and (A2), we can obtain that
v t 2 G 2 + t T v s , γ | σ s ξ , L ( ξ ) , α | 2 D 2 v s d s + 2 t T | v s , F s | d s 2 t T v s , ψ s d B s G 2 κ 2 t T D v s 2 d s + 2 C κ t T v s 2 d s + 2 t T | v s , F s | d s 2 t T v s , ψ s d B s .
Notice, in view of Burkholder–Davis–Gundy inequality ([20]), that
2 E sup τ [ t , T ] | τ T v s , ψ s d B s | 2 E | t T v s , ψ s d B s | + 2 E sup τ [ t , T ] | t τ v s , ψ s d B s | 2 E sup τ [ t , T ] | t τ v s , ψ s d B s | + 2 E sup τ [ t , T ] | t τ v s , ψ s d B s | 4 E sup τ [ t , T ] | t τ v s , ψ s d B s | 16 2 E t T v s 2 ψ s 2 d s 1 / 2 .
Hence,
E sup s [ t , T ] v s 2 E [ G 2 + 2 t T | v s , F s | d s + 2 C κ t T v s 2 d s + 16 2 E t T v s 2 ψ s 2 d s 1 / 2 E G 2 + 2 t T | v s , F s | d s + 16 2 t T ψ s 2 d s + 2 C κ t T v s 2 d s + 1 2 E sup s [ t , T ] v s 2 ,
hence,
1 2 E sup s [ t , T ] v s 2 E G 2 + 2 t T | v s , F s | d s + 16 2 t T ψ s 2 d s + 2 C κ t T v s 2 d s E G 2 + 1 + 2 C κ t T sup r [ s , T ] v r 2 d s + t T F s 2 d s + 16 2 t T ψ s 2 d s .
Applying the Gronwall inequality successively, we obtain
E sup s [ t , T ] v s 2 2 E G 2 + t T F s 2 d s + 16 2 t T ψ s 2 d s e ( 1 + 2 C / κ ) T .
Thus, by combining (A5), we can obtain that (A4) holds.
Step 2. We first take the partial derivatives with respect to ξ on both sides of BSPDE (14), and then we may obtain without difficulty that ( u , Φ ) : = ( D v , D ψ ) is a weak solution to the BSPDE for the Dirichlet problem as follows:
d u t ( ξ ) = γ 2 | σ t ξ , L ( ξ ) , α | 2 D 2 u t ( ξ ) + 1 γ 2 D 2 u t ( ξ ) + D F t ( ξ ) d t Φ t ( ξ ) d B t ; u t ( p ) = 0 , u t ( q ) = 0 ; u T ( ξ ) = D G ( ξ , L ( ξ ) ) .
We apply the Itô’s formula of Lemma A.1 in [7] to the square norm again, which arises as
u t 2 + t T Φ s 2 d s D G 2 = t T u s , γ | σ t ξ , L ( ξ ) , α | 2 D 2 u t ( ξ ) + ( 1 γ ) D 2 u s + D F s d s 2 t T u s , Φ s d B s a . s . t [ 0 , T ] .
Considering the zero-Dirichlet condition, we have
t T u s , ( 1 γ ) D 2 u s d s = ( 1 γ ) t T D u s 2 d s , t T u s , D F s d s = t T D u s , F s d s .
Using a method similar for Step 1, we have
E sup s [ t , T ] u s 2 + E t T D u s 2 + Φ s 2 d s C E D G 2 + t T | F s , D u s | d s ,
that is,
E sup s [ t , T ] D v s 2 + E t T D 2 v s 2 + D ψ s 2 d s C E D G 2 + t T | F s , D 2 v s | d s ,
together with (A4), which suggests that
E sup s [ t , T ] v s H 1 , 2 ( [ p , q ] ) 2 + E t T v s H 2 , 2 ( [ p , q ] ) 2 + ψ s H 1 , 2 ( [ p , q ] ) 2 d s C E G H 1 , 2 ( [ p , q ] ) 2 + t T | F s , v s | + | F s , D 2 v s | d s C E G H 1 , 2 ( [ p , q ] ) 2 + 8 C t T F s 2 d s + 1 2 C t T v s H 2 , 2 ( [ p , q ] ) 2 d s ,
with C being dependent on κ and T.
E sup s [ t , T ] v s H 1 , 2 ( [ p , q ] ) 2 + 1 2 E t T v s H 2 , 2 ( [ p , q ] ) 2 d s + E t T ψ s H 1 , 2 ( [ p , q ] ) 2 d s C E G H 1 , 2 ( [ p , q ] ) 2 + 8 C 2 E t T F s 2 d s ,
thus,
E sup s [ t , T ] v s H 1 , 2 ( [ p , q ] ) 2 + E t T v s H 2 , 2 ( [ p , q ] ) 2 + ψ s H 1 , 2 ( [ p , q ] ) 2 d s C E G H 1 , 2 ( [ p , q ] ) 2 + t T F s 2 d s .
By the estimates and the linearity of the associated BSPDE, we obtain the uniqueness of the strong solution. □

References

  1. Pham, H.; Wei, X. Bellman equation and viscosity solutions for mean-field stochastic control problem. Esaim Control Optim. Calc. Var. 2018, 1, 24. [Google Scholar] [CrossRef]
  2. Pham, H.; Wei, X. Dynamic programming for optimal control of stochastic McKean-Vlasov dynamics. SIAM J. Control Optim. 2016, 55, 1069–1101. [Google Scholar] [CrossRef]
  3. Bayraktar, E.; Cosso, A.; Pham, H. Randomized dynamic programming principle and Feynman-Kac representation for optimal control of McKean-Vlasov dynamics. Trans. Am. Math. Soc. 2018, 370, 2115–2160. [Google Scholar] [CrossRef]
  4. Hafayed, M.; Meherrem, S.; Eren, Ş.; Guçoglu, D.H. On optimal singular control problem for general Mckean-Vlasov differential equations: Necessary and sufficient optimality conditions. Optim. Control Appl. Methods 2018, 3, 39. [Google Scholar] [CrossRef]
  5. Ferrari, G. On a Class of Singular Stochastic Control Problems for Reflected Diffusions. J. Math. Anal. Appl. 2019, 473, 952–979. [Google Scholar] [CrossRef]
  6. Baurdoux, E.J.; Yamazaki, K. Optimality of doubly reflected Levy processes in singular control. Stoch. Process. Their Appl. 2015, 125, 2727–2751. [Google Scholar] [CrossRef]
  7. Bayraktar, E.; Qiu, J. Controlled reflected SDEs and Neumann problem for backward SPDEs. Ann. Appl. Probab. 2019, 29, 2819–2848. [Google Scholar] [CrossRef]
  8. Elliott, R.; Qiu, J.; Wei, W. Neumann problem for backward SPDEs with singular terminal conditions and application in constrained stochastic control under target zone. Stoch. Process. Their Appl. 2022, 148, 68–97. [Google Scholar] [CrossRef]
  9. Wang, F. Distribution Dependent Reflecting Stochastic Differential Equations. Sci. China Math. 2023, 66, 2411–2456. [Google Scholar] [CrossRef]
  10. Cardaliaguet, P. Notes on Mean Field Games. 2012. Available online: https://www.ceremade.dauphine.fr/~cardaliaguet/MFG20130420.pdf (accessed on 28 March 2024).
  11. Chassagneux, J.F.; Crisan, D.; Delarue, F. A Probabilistic approach to classical solutions of the master equation for large population equilibria. Mem. Am. Math. Soc. 2022, 280, 1–98. [Google Scholar] [CrossRef]
  12. Carmona, R.; Delarue, F. The Master Equation for Large Population Equilibriums. In Springer Proceedings in Mathematics & Statistics; Springer: Cham, Switzerland, 2014; Volume 100, pp. 77–128. [Google Scholar]
  13. Adams, D.; Reis, G.D.; Ravaille, R.; Salkeld, W.; Tugaut, J. Large Deviations and Exit-times for reflected McKean-Vlasov equations with self-stabilizing terms and superlinear drifts. Stoch. Process. Their Appl. 2022, 146, 264–310. [Google Scholar] [CrossRef]
  14. Peng, S. Stochastic Hamilton-Jacobi-Bellman Equations. Siam J. Control Optim. 1992, 30, 284–304. [Google Scholar] [CrossRef]
  15. Du, K.; Qiu, J.; Tang, S. Lp Theory for Super-Parabolic Backward Stochastic Partial Differential Equations in the Whole Space. Appl. Math. Optim. 2012, 65, 175–219. [Google Scholar] [CrossRef]
  16. Lieberman, G.M. Second Order Parabolic Differential Equations; World Scientific: Singapore, 1996. [Google Scholar]
  17. Prato, G.D.; Zabczyk, J. Stochastic equations in infinite dimensions. In Encyclopedia of Mathematics and Its Applications; Cambridge University Press: Cambridge, UK, 1992. [Google Scholar]
  18. Kunita, H. Some extensions of Itô’s formula. Semin. Probab. XV 1981, 850, 118–141. [Google Scholar]
  19. Du, K.; Tang, S. Strong solution of backward stochastic partial differential equation in C2 domains. Probab. Theory Relat. Fields 2012, 154, 255–285. [Google Scholar] [CrossRef]
  20. Mao, X. Stochastic Differential Equations and Applications, 2nd ed.; Horwood Publishing: Chichester, UK, 2007. [Google Scholar]
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Ma, L.; Sun, F.; Han, X. Controlled Reflected McKean–Vlasov SDEs and Neumann Problem for Backward SPDEs. Mathematics 2024, 12, 1050. https://doi.org/10.3390/math12071050

AMA Style

Ma L, Sun F, Han X. Controlled Reflected McKean–Vlasov SDEs and Neumann Problem for Backward SPDEs. Mathematics. 2024; 12(7):1050. https://doi.org/10.3390/math12071050

Chicago/Turabian Style

Ma, Li, Fangfang Sun, and Xinfang Han. 2024. "Controlled Reflected McKean–Vlasov SDEs and Neumann Problem for Backward SPDEs" Mathematics 12, no. 7: 1050. https://doi.org/10.3390/math12071050

APA Style

Ma, L., Sun, F., & Han, X. (2024). Controlled Reflected McKean–Vlasov SDEs and Neumann Problem for Backward SPDEs. Mathematics, 12(7), 1050. https://doi.org/10.3390/math12071050

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop