Skip to main content

Mean-field backward doubly stochastic differential equations and related SPDEs

Abstract

Existence and uniqueness result of the solutions to mean-field backward doubly stochastic differential equations (BDSDEs in short) with locally monotone coefficients as well as the comparison theorem for these equations are established. As a preliminary step, the existence and uniqueness result for the solutions of mean-field BDSDEs with globally monotone coefficients is also established. Furthermore, we give the probabilistic representation of the solutions for a class of stochastic partial differential equations by virtue of mean-field BDSDEs, which can be viewed as the stochastic Feynman-Kac formula for SPDEs of mean-field type.

1 Introduction

In this paper, we study a new kind of stochastic partial differential equations (SPDEs):

u ( t , x ) = E [ Φ ( X T 0 , x 0 , x ) ] + t T L u ( s , x ) d s + t T E [ f ( s , X s 0 , x 0 , x , u ( s , X s 0 , x 0 ) , u ( s , x ) , σ ˆ u ( s , X s 0 , x 0 ) , σ ˆ u ( s , x ) ) ] d s + t T E [ g ( s , X s 0 , x 0 , x , u ( s , X s 0 , x 0 ) , u ( s , x ) , σ ˆ u ( s , X s 0 , x 0 ) , σ ˆ u ( s , x ) ) ] d B s ,
(1.1)

where σ ˆ is the transpose of σ ˆ which is defined by σ ˆ :=E[σ(s, X s 0 , x 0 ,x)], and is a second-order differential operator given by ( L u ) i = ( L u i ) 1 i n with

L:= 1 2 i , j = 1 d a i j 2 x i x j + i = 1 d E [ b i ( t , X t 0 , x 0 , x ) ] x i ,

and

a:=( a i , j )= ( E [ σ ( t , X t 0 , x 0 , x ) ] E [ σ ( t , X t 0 , x 0 , x ) ] ) .

Here, the function u(t,x):[0,T]× R d R n is the unknown function, and { B t ,0tT} is an l-dimensional Brownian motion process defined on a given complete probability space (Ω,F,P). X t 0 , x 0 , a stochastic process starting at x 0 when t=0, is the solution of one class of stochastic differential equations (SDEs), and E denotes expectation with respect to the probability P. In this paper, we call this kind of equations (1.1) McKean-Vlasov SPDEs, because they are analogous to McKean-Vlasov PDEs except the stochastic term d B t .

McKean-Vlasov PDEs involving models of large stochastic particle systems with mean-field interaction have been studied by stochastic methods in recent years (see [14] and the references therein). Mean-field approaches have applications in many areas such as statistical mechanics and physics, quantum mechanics and quantum chemistry. Recently, Lasry and Lions introduced mean-field approaches for high-dimensional systems of evolution equations corresponding to a large number of ‘agents’ or ‘particles’. They extended the field of such mean-field approaches to problems in economics, finance and game theory (see [5] and the references therein).

As is well known, to give a probabilistic representation (Feynman-Kac formula) of quasilinear parabolic SPDEs, Pardoux and Peng [6] introduced a new class of backward stochastic differential equations (BSDEs) called backward doubly stochastic differential equations which have two different types of stochastic integrals: a standard (forward) stochastic integral d W t and a backward stochastic integral d B t . They proved the existence and uniqueness for solutions of BDSDEs under uniformly Lipschitz coefficients. When the coefficients are smooth enough, they also established the connection between BDSDEs and a certain kind of quasilinear SPDEs. BDSDEs have a practical background in finance. The extra noise B can be regarded as some extra inside information in a derivative security market. Since 1990s, BDSDEs have drawn more attention from many authors (cf. [713] and the references therein). Shi, Gu and Liu gave the comparison theorem of BDSDEs and investigated the existence of solutions for BDSDEs with continuous coefficients in [11]. To relax the Lipschitz conditions, Wu and Zhang studied two kinds of BDSDEs under globally (respectively, locally) monotone assumptions and obtained the uniqueness and existence results of the solutions (see [12]).

Mean-field BSDEs are deduced by Buckdahn, Djehiche, Li and Peng [14] when they studied a special mean-field problem with a purely stochastic method. Later, Buckdahn, Li and Peng [15] investigated the properties of these equations in a Markovian framework, obtained the uniqueness of the solutions of mean-field BSDEs as well as the comparison theorem and also gave the viscosity solutions of a class of McKean-Vlasov PDEs in terms of mean-field BSDEs.

In this paper, we study a new type of BDSDEs, that is, the so called mean-field BDSDEs, under the globally (respectively locally) monotone coefficients. We obtain the existence and uniqueness result of the solution by virtue of the technique proposed by Wu and Zhang [12] and the contraction mapping theorem under certain conditions. Also, the comparison principle for mean-field BDSDEs is discussed when the coefficients satisfy some stricter assumptions. A comparison theorem is a useful result in the theory of BSDEs. For instance, it can be used to study viscosity solutions of PDEs. Here, we point out that it is more delicate to prove the comparison theorem for mean-field BDSDEs because of the mean-field term.

We also present the connection between McKean-Vlasov SPDEs and mean-field BDSDEs. In detail, let { X s t , x ,tsT} be the solution of

{ d X s t , x = E [ b ( s , ( X s 0 , x 0 ) , X s t , x ) ] d s + E [ σ ( s , ( X s 0 , x 0 ) , X s t , x ) ] d W s , s [ t , T ] , X t t , x = x .

Assume that Eq. (1.1) has a classical solution. Then the couple ( Y s t , x , Z s t , x ), where Y s t , x =u(s, X s t , x ) and Z s t , x = E [ σ ( s , ( X s 0 , x 0 ) , X s t , x ) ] Du(s, X s t , x ), verifies the following mean-field BDSDE :

Y s t , x = E [ Φ ( ( X T 0 , x 0 ) , X T t , x ) ] + s T E [ f ( r , ( X r 0 , x 0 ) , X r t , x , ( Y r 0 , x 0 ) , Y r t , x , ( Z r 0 , x 0 ) , Z r t , x ) ] d r + s T E [ g ( r , ( X r 0 , x 0 ) , X r t , x , ( Y r 0 , x 0 ) , Y r t , x , ( Z r 0 , x 0 ) , Z r t , x ) ] d B r s T Z r t , x d W r .
(1.2)

In Eq. (1.2), the integral d W t is a forward Itô integral, and the integral d B t denotes a backward Itô integral. { W t ,0tT} and { B t ,0tT} are two mutually independent standard Brownian motion processes with values respectively in R d and in R l . This conclusion gives a probabilistic representation of McKean-Vlasov SPDEs (1.1), which can be regarded as a stochastic Feynman-Kac formula for Mckean-Vlasov SPDEs.

Our paper is organized as follows. In Section 2, we present the existence and uniqueness results about mean-field BDSDEs with globally monotone coefficients. We investigate the properties of mean-field BDSDEs with locally monotone assumptions in Section 3. We first prove the existence and uniqueness of the solutions of mean-field BDSDEs and then derive the comparison theorem when the mean-field BDSDEs are one-dimensional. In Section 4, we introduce the decoupled mean-field forward-backward doubly stochastic differential equation and study the regularity of its solution with respect to x, which is the initial condition of the McKean-Vlasov SDE. Finally, Section 5 is devoted to the formulation of McKean-Vlasov SPDEs and provides the relationship between the solutions of SPDEs and those of mean-field BDSDEs.

2 Mean-field BDSDEs with globally monotone coefficients

In this section, we study mean-field BDSDEs with globally monotone coefficients, which is helpful for the case of locally monotone coefficients. To this end, we firstly introduce some notations and recall some results on mean-field BSDEs obtained by Buckdahn, Li and Peng [15].

Let { W t ,0tT} and { B t ,0tT} be two mutually independent standard Brownian motion processes, with values respectively in R d and R l , defined over some complete probability space (Ω,F,P), where T is a fixed positive number throughout this paper. Moreover, let N denote the class of P-null sets of . For each 0sT, we define

F t F 0 , t W F t , T B

with F s , t W =σ{ W r ;srt}N and F t , s B =σ{ B r B t ;trs}N.

Note that { F t ,t[0,T]} is not an increasing family of σ-fields, so it is not a filtration.

We will also use the following spaces:

  • For any nN, let H F 2 (0,T; R n ) denote the set of (classes of dP×dt a.e. equal) n-dimensional jointly measurable random processes { ψ t ;t[0,T]} which satisfy: Evidently, H F 2 (0,T; R n ) is a Banach space endowed with the canonical norm ψ= { E 0 T | ψ s | 2 d s } 1 2 .

  1. (i)

    E 0 T | ψ t | 2 dt<+,

  2. (ii)

    ψ t is F t measurable, for a.e. 0tT.

  • We denote similarly by S F 2 ([0,T]; R n ) the set of continuous n-dimensional random processes { ψ t ;t[0,T]} which satisfy:

  1. (i)

    E( sup 0 t T | ψ t | 2 )<+,

  2. (ii)

    ψ t is F t measurable, for a.e. 0tT.

  • L 0 (Ω,F,P; R n ) denotes the space of all R n valued -measurable random variables.

  • For 1p<, L p (Ω,F,P; R n ) is the space of all R n valued -measurable random variables such that E[ | ξ | p ]<.

Let ( Ω ¯ , F ¯ , P ¯ )=(Ω×Ω,FF,PP) be the (non-completed) product of (Ω,F,P) with itself, and we define F ¯ ={ F ¯ t =F F t ,0tT} on this product space. A random variable ξ L 0 (Ω,F,P; R n ) originally defined on Ω is extended canonically to Ω ¯ : ξ ( ω ,ω)=ξ( ω ), ( ω ,ω) Ω ¯ =Ω×Ω. For any θ L 1 ( Ω ¯ , F ¯ , P ¯ ), the variable θ(,ω):Ω R n belongs to L 1 (Ω,F,P), P(dω)-a.s., whose expectation is denoted by

E [ θ ( , ω ) ] = Ω θ ( ω , ω ) P ( d ω ) .

Notice that E [θ]= E [θ(,ω)] L 1 (Ω,F,P) and

E ¯ [θ] ( = Ω ¯ θ d P ¯ = Ω E [ θ ( , ω ) ] P ( d ω ) ) =E [ E [ θ ] ] .

Moreover, for all ( y , z ,y,z), f=f( ω ,ω,t, y , z ,y,z): Ω ¯ ×[0,T]× R n × R n × d × R n × R n × d R n , g=g( ω ,ω,t, y , z ,y,z): Ω ¯ ×[0,T]× R n × R n × d × R n × R n × d R n × l are two F ¯ t -measurable functions which satisfy

Assumption 2.1 (A1) g(t,0,0,0,0) H F ¯ 2 (0,T; R n × l ), and there exist L>0 and 0<α< 1 2 such that

(A2) for any fixed ( ω ,ω,t), f( ω ,ω,t,,,,) is continuous;

(A3) there exist a process f ¯ t H F ¯ 2 (0,T;R) and a constant L>0 such that

| f ( t , y , z , y , z ) | f ¯ t +L ( | y | + | z | + | y | + | z | ) ;

(A4) there exist constants λ 1 , λ 2 R such that for all t[0,T], y i , y i R n , z i , z i R n × d (i=1,2),

( y 1 y 2 ) ( f ( t , y 1 , z , y 1 , z ) f ( t , y 2 , z , y 2 , z ) ) λ 1 ( y 1 y 2 ) ( y 1 y 2 ) + λ 2 | y 1 y 2 | 2 ;

(A5) there exists L>0 such that, P ¯ -a.s., for all t[0,T], y ,y R n , z 1 , z 2 , z 1 , z 2 R n × d ,

| f ( t , y , z 1 , y , z 1 ) f ( t , y , z 2 , y , z 2 ) | 2 L ( | z 1 z 2 | 2 + | z 1 z 2 | 2 ) .

We now consider the following mean-field BDSDEs with the form:

Y t = ξ + t T E [ f ( s , Y s , Z s , Y s , Z s ) ] d s + t T E [ g ( s , Y s , Z s , Y s , Z s ) ] d B s t T Z s d W s , 0 t T .
(2.1)

Remark 2.1 Due to our notation, the coefficients of (2.1) are interpreted as follows:

Remark 2.2 If coefficient f meets the following Lipschitz assumption: There exists a constant L>0 such that, P ¯ -a.s., for all t[0,T], y i , y i R n , z i , z i R n × d (i=1,2),

| f ( t , y 1 , z 1 , y 1 , z 1 ) f ( t , y 2 , z 2 , y 2 , z 2 ) | 2 L ( | y 1 y 2 | 2 + | z 1 z 2 | 2 + | y 1 y 2 | 2 + | z 1 z 2 | 2 ) ,

then it must satisfy conditions (A4) and (A5).

Definition 2.1 A pair of F t -measurable processes {( Y t , Z t );0tT} is called a solution of mean-field BDSDE (2.1) if (Y,Z) S F 2 ([0,T]; R n )× H F 2 (0,T; R n × d ) and it satisfies mean-field BDSDE (2.1).

The main result of this section is the following theorem.

Theorem 2.1 For any random variable ξ L 2 (Ω, F T ,P; R n ), under Assumption  2.1, mean-field BDSDE (2.1) admits a unique solution (Y,Z) S F 2 ([0,T]; R n )× H F 2 (0,T; R n × d ).

Proof Step 1: For any (y,z) H F 2 (0,T; R n × R n × d ), BDSDE

Y t = ξ + t T E [ f ( s , y s , z s , Y s , Z s ) ] d s + t T E [ g ( s , y s , z s , Y s , Z s ) ] d B s t T Z s d W s , 0 t T
(2.2)

has a unique solution. In order to get this conclusion, we define

φ ( y , z ) (s,μ,ν):= E [ φ ( s , y s , z s , μ , ν ) ] ,for φ=f,g.

Then (2.2) can be rewritten as

Y t =ξ+ t T f ( y , z ) ( Y s , Z s )ds+ t T g ( y , z ) ( Y s , Z s )d B s t T Z s d W s .

Due to Assumption 2.1, for all (μ,ν),( μ 1 , ν 1 ),( μ 2 , ν 2 ) R n × R n × d , g satisfies

| g ( y , z ) ( μ 1 , ν 1 ) g ( y , z ) ( μ 2 , ν 2 ) | 2 L | μ 1 μ 2 | 2 +α | ν 1 ν 2 | 2 ,

and f fulfills

According to Theorem 2.2 in [12], BDSDE (2.2) has a unique solution.

Step 2: Now, we introduce a norm on the space H F 2 (0,T; R n × R n × d ) which is equivalent to the canonical norm

( y , z ) β = { E 0 T e β s ( c ¯ | y s | 2 + | z s | 2 ) d s } 1 2 , c ¯ ,β>0.

The parameters c ¯ and β are specified later.

From Step 1, we can introduce the mapping ( Y , Z )=I[( y , z )]: H F 2 (0,T; R n × R n × d ) H F 2 (0,T; R n × R n × d ) through the equation

Y t = ξ + t T E [ f ( s , y s , z s , Y s , Z s ) ] d s + t T E [ g ( s , y s , z s , Y s , Z s ) ] d B s t T Z s d W s , 0 t T .

For any ( y 1 , z 1 ),( y 2 , z 2 ) H F 2 (0,T; R n × R n × d ), we set ( Y 1 , Z 1 )=I[( y 1 , z 1 )], ( Y 2 , Z 2 )=I[( y 2 , z 2 )], ( y ¯ , z ¯ )=( y 1 y 2 , z 1 z 2 ) and ( Y ¯ , Z ¯ )=( Y 1 Y 2 , Z 1 Z 2 ). Then applying Itô’s formula to e β s | Y ¯ s | 2 and by virtue of Y 1 , Y 2 S F 2 ([0,T]; R n ), we have

(2.3)

From condition (A4) and noting that E [ Y s ( ω )]=E[ Y s (ω)], for any M>0, we get

Then we have

If we set M= 2 L 1 2 α , c ¯ = 4 1 + 2 α ( | λ 1 | M +L), β=| λ 1 |M+2 λ 2 + +3M+L+ 1 2 c ¯ , then it yields

E [ t T e β s ( c ¯ | Y ¯ s | 2 + | Z ¯ s | 2 ) d s ] 1 + 2 α 2 E [ t T e β s ( c ¯ | y ¯ s | 2 + | z ¯ s | 2 ) d s ] .

Consequently, I is a strict contraction on H F 2 (0,T; R n × R n × d ) equipped with the norm β for 0<α< 1 2 . With the contraction mapping theorem, there admits a unique fixed point (Y,Z) H F 2 (0,T; R n × R n × d ) such that I(Y,Z)=(Y,Z). On the other hand, from Step 1, we know that if I(Y,Z)=(Y,Z), then (Y,Z) S F 2 ([0,T]; R n )× H F 2 (0,T; R n × d ), which is the unique solution of Eq. (2.1). □

Suppose that: For some f: Ω ¯ ×[0,T]× R n × R n × d × R n × R n × d R n satisfying (A2)-(A5), the generators f i , i=1,2 are of the form

f i ( s , ( Y s i ) , ( Z s i ) , Y s i , Z s i ) =f ( s , ( Y s i ) , ( Z s i ) , Y s i , Z s i ) + φ i (s),dsd P ¯ -a.e.,i=1,2,

where φ i H F ¯ 2 (0,T; R n ). Then we have the following corollary.

Corollary 2.1 Suppose that ( Y i , Z i ) is the solution of mean-field BDSDE (2.1) with data ( ξ i , f i ,g), i=1,2, where ξ 1 , ξ 2 L 2 (Ω, F T ,P; R n ) are two arbitrary terminal values. The difference of ( Y 1 , Z 1 ) and ( Y 2 , Z 2 ) satisfies the following estimate:

(2.4)

where β=2| λ 1 |+2 λ 2 + + 4 L 1 2 α +2+2L.

The proof of the above corollary is similar to that of Theorem 2.1 and is therefore omitted.

3 Mean-field BDSDEs with locally monotone coefficients

In this section, we investigate mean-field BDSDEs with locally monotone coefficients. The results can be regarded as an extension of the results in [12] to the mean-field type.

We assume

(A3′) there exist L>0 and 0γ<1 such that |f(t, y , z ,y,z)|L(1+ | y | γ + | z | γ + | y | γ + | z | γ );

(A4′) for any NN, there exist constants λ N , λ ¯ N R such that, y i , y i , z ,z satisfying | y i |,| y i |,| z |,|z|N (i=1,2), we have

( y 1 y 2 ) ( f ( t , y 1 , z , y 1 , z ) f ( t , y 2 , z , y 2 , z ) ) λ N ( y 1 y 2 ) ( y 1 y 2 ) + λ ¯ N | y 1 y 2 | 2 ;

(A5′) NN, there exists L N >0 such that, for any y , y, z i , z i satisfying | y |,|y|,| z i |,| z i |N (i=1,2), it holds

| f ( t , y , z 1 , y , z 1 ) f ( t , y , z 2 , y , z 2 ) | 2 L N ( | z 1 z 2 | 2 + | z 1 z 2 | 2 ) .

Remark 3.1 Since | x | γ 1+|x|, γ[0,1), (A3′) implies that

| f ( t , y , z , y , z ) | L ( 5 + | y | + | z | + | y | + | z | ) .

We need the following lemma, which plays an important role in the proof of the main result.

Lemma 3.1 Under (A2), (A3′)-(A5′) there exists a sequence of { f m } m = 1 such that

  1. (i)

    for fixed mN, ω , ω, t, f m ( ω ,ω,t,,,,) is continuous;

  2. (ii)

    m, | f m (t, y , z ,y,z)||f(t, y , z ,y,z)|L(1+ | y | γ + | z | γ + | y | γ + | z | γ );

  3. (iii)

    N, ρ N ( f m f)0 as m, where ρ m 2 (f):=E[ 0 T sup | y | , | z | , | y | , | z | m | f ( t , y , z , y , z ) | 2 dt];

  4. (iv)

    m, f m is globally monotone in y; moreover, for any m, N with mN, it holds that

    ( y 1 y 2 ) ( f m ( t , y 1 , z , y 1 , z ) f m ( t , y 2 , z , y 2 , z ) ) λ N ( y 1 y 2 ) ( y 1 y 2 ) + λ ¯ N | y 1 y 2 | 2 ,

for any t, y i , y i , z , z satisfying | y i |,| y i |,| z |,|z|N (i=1,2);

  1. (v)

    m, f m is globally Lipschitz in z , z; moreover, for any m, N with mN, it holds that

    | f m ( t , y , z 1 , y , z 1 ) f m ( t , y , z 2 , y , z 2 ) | 2 L N ( | z 1 z 2 | 2 + | z 1 z 2 | 2 ) ,

for any t, y , y, z i , z i satisfying | y |,|y|,| z i |,| z i |N (i=1,2).

Proof We define f m by

f m ( t , y , z , y , z ) =f ( t , y , z , y , z ) ϕ m ( y ) φ m ( z ) ψ m (y) η m (z),

where ϕ m : R n R + is a sequence of smooth functions such that 0 ϕ m 1, ϕ m (x)=1 for |x|m, and ϕ m (x)=0 for |x|m+1. Similarly, we define the sequences φ m : R n × d R + , ψ m : R n R + , η m : R n × d R + . It should be pointed out that ϕ m , φ m , ψ m and η m are continuously differentiable with bounded derivatives for each m. The conclusion of this lemma can be easily obtained by arguments similar to those of Lemma 3.3 in [12]. □

We now present the main result of this section.

Theorem 3.1 Let (A1), (A2), (A3′)-(A5′) hold. Assume, moreover,

1 + exp ( 2 L + 2 | λ N | + 2 λ ¯ N + + 2 L N θ 1 + 2 ) N 2 ( 1 γ ) 0, as N,
(3.1)

where θ is an arbitrarily fixed constant such that 0<θ<12α. Then mean-field BDSDE (2.1) has a unique solution (Y,Z) S F 2 ([0,T]; R n )× H F 2 (0,T; R n × d ).

Proof We now construct an approximate sequence. Let f m be associated to f by Lemma 3.1. Then for each m, f m is globally monotone in y and globally Lipschitz in z. By Theorem 2.1, the following mean-field BDSDE

Y t m = ξ + t T E [ f m ( s , ( Y s m ) , ( Z s m ) , Y s m , Z s m ) ] d s + t T E [ g ( s , ( Y s m ) , ( Z s m ) , Y s m , Z s m ) ] d B s t T Z s m d W s , 0 t T ,
(3.2)

admits a unique solution ( Y m , Z m ) S F 2 ([0,T]; R n )× H F 2 (0,T; R n × d ) for each mN. Applying Itô’s formula to | Y t m | 2 yields

where

and

Hence,

Then it follows from Gronwall’s inequality and the B-D-G inequality that

where C>0 only depends on T, α, L and is independent of m.

For any m,kN, set

A : = { ( ω , ω , s ) : | ( Y s m ) | + | ( Z s m ) | + | ( Y s k ) | + | ( Z s k ) | + | Y s m | + | Z s m | + | Y s k | + | Z s k | N } ,

and A ¯ :=ΩA.

Next, we will conclude that ( Y m , Z m ) is a Cauchy sequence in S F 2 ([0,T]; R n )× H F 2 (0,T; R n × d ). Actually, since mean-field BDSDE

Y t k = ξ + t T E [ f k ( s , ( Y s k ) , ( Z s k ) , Y s k , Z s k ) ] d s + t T E [ g ( s , ( Y s k ) , ( Z s k ) , Y s k , Z s k ) ] d B s t T Z s k d W s , 0 t T ,

admits a unique solution ( Y k , Z k ) S F 2 ([0,T]; R n )× H F 2 (0,T; R n × d ) for each kN. Applying Itô’s formula to | Y t m Y t k | 2 , we have

(3.3)

where

We next estimate I, II and III.

For the first term I, based on Hölder’s inequality and Chebyshev’s inequality, we have

I E t T | Y s m Y s k | 2 d s + E [ t T E [ | f m ( s , ( Y s m ) , ( Z s m ) , Y s m , Z s m ) f k ( s , ( Y s k ) , ( Z s k ) , Y s k , Z s k ) | 2 ] I A d s ] E t T | Y s m Y s k | 2 d s + C N 2 ( 1 γ ) ,
(3.4)

where C>0 depends on T, L, α and E[ 0 T | E [ g ( s , 0 , 0 , 0 , 0 ) ] | 2 ds].

For the second term II, due to the local monotonicity of f m in y and the local Lipschitz condition of f m in z, we obtain that for M>0, the following holds:

II = 2 E [ t T ( Y s m Y s k ) E [ f m ( s , ( Y s m ) , ( Z s m ) , Y s m , Z s m ) f m ( s , ( Y s k ) , ( Z s m ) , Y s k , Z s m ) ] I A ¯ d s ] + 2 E [ t T ( Y s m Y s k ) E [ f m ( s , ( Y s k ) , ( Z s m ) , Y s k , Z s m ) f m ( s , ( Y s k ) , ( Z s k ) , Y s k , Z s k ) ] I A ¯ d s ] 2 E t T [ λ N ( Y s m Y s k ) E ( Y s m Y s k ) + λ ¯ N | Y s m Y s k | 2 ] d s + M E t T | Y s m Y s k | 2 d s + L N M E t T [ E [ | ( Z s m ) ( Z s k ) | 2 ] + | Z s m Z s k | 2 ] d s ( 2 | λ N | + 2 λ ¯ N + + M ) E t T | Y s m Y s k | 2 d s + 2 L N M E t T | Z s m Z s k | 2 d s .
(3.5)

For the last term, we have

IIIE t T | Y s m Y s k | 2 ds+ ρ N 2 ( f m f k ).
(3.6)

Choose M>0 such that θ:= 2 L N M <12α. Then from (3.3)-(3.6), we obtain

Applying Gronwall’s inequality and the B-D-G inequality to the above inequality yields

where c>0 is independent of m, k. Now passing to the limit successively on m, k and N, we see that ( Y m , Z m ) is a Cauchy (hence convergent) sequence in S F 2 ([0,T]; R n )× H F 2 (0,T; R n × d ); denote the limit by (Y,Z), which satisfies

E [ sup 0 t T | Y t m Y t | 2 ] +E [ 0 T | Z s m Z s | 2 d s ] 0

as m.

Next, we show that (Y,Z) is the solution of mean-field BDSDE (2.1). To this end, we only need to prove that the following conclusion holds along a subsequence:

(3.7)

Set

A m := { ( ω , ω , s ) : | ( Y s m ) | + | ( Z s m ) | + | Y s | + | Z s | + | Y s m | + | Z s m | + | Y s | + | Z s | N } ,

and A ¯ m :=Ω A m .

Since

f m ( s , ( Y s m ) , ( Z s m ) , Y s m , Z s m ) f ( s , Y s , Z s , Y s , Z s ) = I 1 (m,s)+ I 2 (m,s)+ I 3 (m,s)+ I 4 (m,s),

where

then we have

(3.8)

where C>0 is independent of m. As

E [ sup 0 t T | Y t m Y t | 2 ] +E [ 0 T | Z s m Z s | 2 d s ] 0and sup m N E [ sup 0 t T | Y t m | 2 ] <,

there exists a subsequence of Y m , still denoted by Y m , such that Y t m Y t a.e., a.s. It then follows from the continuity of f in y and the dominated convergence theorem that

E [ t T E [ | f ( s , ( Y s m ) , Z s , Y s m , Z s ) f ( s , Y s , Z s , Y s , Z s ) | 2 ] d s ] 0

as m.

Now, passing to the limit as m and N in (3.8) successively, it follows that (3.7) holds. Then letting m in (3.2) yields

Y t = ξ + t T E [ f ( s , Y s , Z s , Y s , Z s ) ] d s + t T E [ g ( s , Y s , Z s , Y s , Z s ) ] d B s t T Z s d W s , 0 t T .

Therefore, we come to the conclusion of this theorem. □

Now, we discuss the comparison theorem for mean-field BDSDEs. We only consider one-dimensional mean-field BDSDEs, i.e., n=1.

We consider the following mean-field BDSDEs: (0tT)

(3.9)
(3.10)

Theorem 3.2 (Comparison theorem)

Assume mean-field BDSDEs (3.9) and (3.10) satisfy the conditions of Theorem  3.1. Let ( Y 1 , Z 1 ) and ( Y 2 , Z 2 ) be the solutions of mean-field BDSDEs (3.9) and (3.10), respectively. Moreover, for the two generators of f 1 and f 2 , we suppose:

  1. (i)

    One of the two generators is independent of z .

  2. (ii)

    One of the two generators is nondecreasing in y .

Then if ξ 1 ξ 2 , a.s., f 1 (t, y , z ,y,z) f 2 (t, y , z ,y,z), a.s., there also holds that Y t 1 Y t 2 , a.s. t[0,T].

Remark 3.2 The conditions (i) and (ii) of Theorem 3.2 are, in particular, satisfied if they hold for the same generator f j (j=1,2), but also if (i) is satisfied by one generator and (ii) by the other one.

Proof Without loss of generality, we suppose that (i) is satisfied by f 1 and (ii) by f 2 . For notational simplicity, we set ξ ¯ := ξ 1 ξ 2 , ( Y ¯ , Z ¯ ):=( Y 1 Y 2 , Z 1 Z 2 ), then

Y ¯ t = ξ ¯ + t T E [ f 1 ( s , ( Y s 1 ) , Y s 1 , Z s 1 ) f 2 ( s , ( Y s 2 ) , ( Z s 2 ) , Y s 2 , Z s 2 ) ] d s + t T E [ g ( s , ( Y s 1 ) , ( Z s 1 ) , Y s 1 , Z s 1 ) g ( s , ( Y s 2 ) , ( Z s 2 ) , Y s 2 , Z s 2 ) ] d B s t T Z ¯ s d W s , 0 t T .

By Itô’s formula applied to | Y ¯ t + | 2 and noting that ξ 1 ξ 2 , it easily follows that

(3.11)

Since f 1 (t, y , z ,y,z) f 2 (t, y , z ,y,z) a.s. and f 2 is nondecreasing in y , we have

Then we have

Δ : = 2 E [ Y ¯ s + ( E [ f 1 ( s , ( Y s 1 ) , Y s 1 , Z s 1 ) f 2 ( s , ( Y s 2 ) , ( Z s 2 ) , Y s 2 , Z s 2 ) ] ) ] ( 2 λ ¯ N + + L N 1 2 α + 2 | λ N | ) E [ | Y ¯ s + | 2 ] + ( 1 2 α ) E [ I { Y ¯ s > 0 } | Z ¯ s | 2 ] .
(3.12)

With the assumption (A1), we obtain

(3.13)

Combining (3.12), (3.13) with (3.11) yields

E [ | Y ¯ t + | 2 ] ( 2 λ ¯ N + + L N 1 2 α + 2 | λ N | + 2 L ) t T E [ ( Y ¯ s + ) 2 ] ds.

By Gronwall’s inequality, it follows that

E [ | Y ¯ t + | 2 ] =0,t[0,T],

that is, Y t 1 Y t 2 , P-a.s., t[0,T]. □

4 Decoupled mean-field forward-backward doubly SDEs

In this section, we study the decoupled mean-field forward-backward doubly stochastic differential equations. First, we recall some results of Buckdahn, Li and Peng [15] on McKean-Vlasov SDEs. Given continuous functions b: Ω ¯ ×[0,T]× R d × R d R d and σ: Ω ¯ ×[0,T]× R d × R d R d × d which are supposed to satisfy the following conditions:

Assumption 4.1 (i) b(t,0,0) and σ(t,0,0) are F ¯ t -measurable continuous processes and there exists some constant L>0 such that

| b ( t , x , x ) | + | σ ( t , x , x ) | L ( 1 + | x | ) ,a.s., for all 0tT,x, x R d ;
  1. (ii)

    b and σ are Lipschitz in x, x , i.e., there is some constant L>0 such that

For any x 0 R d , we consider the following SDE parameterized by the initial condition (t,ζ)[0,T]× L 2 (Ω, F t ,P; R d ):

{ d X s t , ζ = E [ b ( s , ( X s 0 , x 0 ) , X s t , ζ ) ] d s + E [ σ ( s , ( X s 0 , x 0 ) , X s t , ζ ) ] d W s , s [ t , T ] , X t t , ζ = ζ .
(4.1)

From the result about Eq. (5.1) in [15], we know that under Assumption 4.1, SDE (4.1) has a unique strong solution, and we can obtain that X T t , ζ has a continuous version with the following well-known standard estimates.

Proposition 4.1 p2, there exists C p R + such that, for all t[0,T] and ζ, ζ L p (Ω, F t ,P; R d ),

(4.2)

for all δ>0 with t+δT.

Now, let f(t, x ,x, y ,y, z ,z), g(t, x ,x, y ,y, z ,z) and Φ( x ,x) be real-valued functions and satisfy the following conditions.

Assumption 4.2 (i) Φ: Ω ¯ × R d × R d R n is an F ¯ T B( R d )-measurable random variable, f: Ω ¯ ×[0,T]× R d × R d × R n × R n × R n × d × R n × d R n and g: Ω ¯ ×[0,T]× R d × R d × R n × R n × R n × d × R n × d R n × l are two measurable processes such that f(t, x ,x, y ,y, z ,z), g(t, x ,x, y ,y, z ,z) are F ¯ t -measurable, for all ( x ,x, y ,y, z ,z) R d × R d × R n × R n × R n × d × R n × d .

  1. (ii)

    For all 0tT, x 1 , x 1 , x 2 , x 2 R d , y 1 , y 1 , y 2 , y 2 R n , z 1 , z 1 , z 2 , z 2 R n × d , there exist constants L>0, λ 1 , λ 2 R and 0<α< 1 2 such that

  2. (iii)

    f, g and Φ satisfy a linear growth condition, i.e., there exists some L>0 such that, a.s., for all x ,x R d

    | Φ ( x , x ) | + | f ( t , x , x , 0 , 0 , 0 , 0 ) | + | g ( t , x , x , 0 , 0 , 0 , 0 ) | L ( 1 + | x | + | x | ) .

Next, we investigate the solution of the following BDSDE:

{ d Y s t , ζ = E [ f ( s , ( X s 0 , x 0 ) , X s t , ζ , ( Y s 0 , x 0 ) , Y s t , ζ , ( Z s 0 , x 0 ) , Z s t , ζ ) ] d s d Y s t , ζ = + E [ g ( s , ( X s 0 , x 0 ) , X s t , ζ , ( Y s 0 , x 0 ) , Y s t , ζ , ( Z s 0 , x 0 ) , Z s t , ζ ) ] d B s Z s t , ζ d W s , Y T t , ζ = E [ Φ ( ( X T 0 , x 0 ) , X T t , ζ ) ] , s [ t , T ] .
(4.3)

Firstly, we study the case (t,ζ)=(0, x 0 ). From Theorem 2.1, we know that there exists a unique solution ( Y 0 , x 0 , Z 0 , x 0 ) S F 2 (0,T; R n )× H F 2 (0,T; R n × d ) to the mean-field BDSDE (4.3). Once we have ( Y 0 , x 0 , Z 0 , x 0 ), Eq. (4.3) becomes a classical BDSDE with coefficients

and Φ ˜ (ω, X T t , ζ (ω))= E [Φ(,ω, ( X T 0 , x 0 ) , X T t , ζ )] L 2 (Ω, F T ,P; R n ). Then due to Theorem 2.2 in [12], we obtain that there exists a unique solution ( Y t , ζ , Z t , ζ ) S F 2 ([0,T]; R n )× H F 2 (0,T; R n × d ) to Eq. (4.3).

For BDSDE (4.3), we give the following proposition.

Proposition 4.2 For any t[0,T] and ζ, ζ L 2 (Ω, F t ,P; R d ), there exists a constant C>0 such that

(4.4)
(4.5)

Proof Combining classical BDSDE estimates (see the proof of Theorem 2.1 in Pardoux and Peng [6]) with the techniques presented in Theorem 3.1, we can get the proof easily. □

5 Mean-field BDSDEs and McKean-Vlasov SPDEs

We now pay attention to investigation of the following system of quasilinear backward stochastic partial differential equations which are called McKean-Vlasov SPDEs: for any (t,x)[0,T]× R d ,

u ( t , x ) = E [ Φ ( X T 0 , x 0 , x ) ] + t T L u ( s , x ) d s + t T E [ f ( s , X s 0 , x 0 , x , u ( s , X s 0 , x 0 ) , u ( s , x ) , σ ˆ u ( s , X s 0 , x 0 ) , σ ˆ u ( s , x ) ) ] d s + t T E [ g ( s , X s 0 , x 0 , x , u ( s , X s 0 , x 0 ) , u ( s , x ) , σ ˆ u ( s , X s 0 , x 0 ) , σ ˆ u ( s , x ) ) ] d B s ,
(5.1)

with σ ˆ :=E[σ(s, X s 0 , x 0 ,x)]= E [σ(s, ( X s 0 , x 0 ) ,x)], and u: R + × R d R n ,

Lu=( L u 1 L u n )

with

L:= 1 2 i , j = 1 d a i j 2 x i x j + i = 1 d E [ b i ( t , X t 0 , x 0 , x ) ] x i ,

where

a:=( a i , j )= ( E [ σ ( t , X t 0 , x 0 , x ) ] E [ σ ( t , X t 0 , x 0 , x ) ] ) .

Note that

In fact, Eq. (5.1) is a new kind of nonlocal SPDE because of the mean-field term. Here, the functions b, σ, f, g and Φ are supposed to satisfy Assumption 4.1 and Assumption 4.2 respectively, and X 0 , x 0 is the solution of the mean-field SDE (4.1) with (t,ξ)=(0, x 0 ).

Now, we give the main theorem of this section.

Theorem 5.1 Suppose that Assumption  4.1 and Assumption  4.2 hold. Let {u(t,x);0tT,x R d } be a F t , T B -measurable random field such that u(t,x) satisfies Eq. (5.1) and for each (t,x), u C 0 , 2 ([0,T]× R d ; R n ) a.s. Moreover, we assume that f,gC([0,T]× R d × R d × R n × R n × R n × d × R n × d ) for a.s. ω ¯ Ω ¯ .

Then we have u(t,x)= Y t t , x , where { ( Y s t , x , Z s t , x ) ; t s T } t 0 , x R d is the unique solution of the mean-field BDSDEs (4.3) and

Y s t , x =u ( s , X s t , x ) , Z s t , x = E [ σ ( s , ( X s 0 , x 0 ) , X s t , x ) ] u ( s , X s t , x ) .
(5.2)

Proof It suffices to show that {u(s, X s t , x ), E [ σ ( s , ( X s 0 , x 0 ) , X s t , x ) ] u(s, X s t , x );0ts} solves the mean-field BDSDE (4.3). To simplify the notation, we define

φ ˆ ( s , X s 0 , x 0 , x ) E [ φ ( s , X s 0 , x 0 , x , u ( s , X s 0 , x 0 ) , u ( s , x ) , σ ˆ u ( s , X s 0 , x 0 ) , σ ˆ u ( s , x ) ) ] , for  φ = f , g .

According to our notations introduced in Section 2, we know that

φ ˆ ( s , X s 0 , x 0 , x ) = E [ φ ( ω , s , X s 0 , x 0 ( ω ) , x , u ( s , X s 0 , x 0 ( ω ) ) , u ( s , x ) , σ ˆ u ( s , X s 0 , x 0 ( ω ) ) , σ ˆ u ( s , x ) ) ] .

Let t= t 0 < t 1 < t 2 << t n =T and λ:=max{ t i + 1 t i }0. For each t i [t,T), applying Itô’s formula to u( t i , X t i t , x ) and noticing that u satisfies Eq. (5.1), we get

The condition that u C 0 , 2 ([0,T]× R d ; R n ) and the continuity of f and g are adopted in the last equation.

Then we have

u ( t , x ) u ( T , X T t , x ) = i = 0 n 1 [ u ( t i , X t i t , x ) u ( t i + 1 , X t i + 1 t , x ) ] = i = 0 n 1 t i t i + 1 f ˆ ( s , X s 0 , x 0 , X s t , x ) d s + i = 0 n 1 t i t i + 1 g ˆ ( s , X s 0 , x 0 , X s t , x ) d B s i = 0 n 1 t i t i + 1 E [ σ ( s , ( X s 0 , x 0 ) , X s t , x ) ] u ( s , X s t , x ) d W s = t T f ˆ ( s , X s 0 , x 0 , X s t , x ) d s + t T g ˆ ( s , X s 0 , x 0 , X s t , x ) d B s t T E [ σ ( s , ( X s 0 , x 0 ) , X s t , x ) ] u ( s , X s t , x ) d W s .

So, Y s t , x :=u(s, X s t , x ), Z s t , x := E [ σ ( s , ( X s 0 , x 0 ) , X s t , x ) ] u(s, X s t , x ) solves the mean-field BDSDE (4.3). The proof is now complete. □

Remark 5.1 Formula (5.2) generalizes the stochastic Feynman-Kac formula for SPDEs of the mean-field type.

References

  1. Bossy M: Some stochastic particle methods for nonlinear parabolic PDEs. ESAIM Proc. 2005, 15: 18-57.

    MathSciNet  Google Scholar 

  2. Méléard S: Asymptotic behaviour of some interacting particle systems; McKean-Vlasov and Boltzmann models. Lectures Notes in Math. 1627. In Probabilistic Models for Nonlinear PDE’s. Edited by: Talay D, Tubaro L. Springer, Berlin; 1996:42-95.

    Google Scholar 

  3. Bossy M, Talay D: A stochastic particle method for the McKean-Vlasov and the Burgers equation. Math. Comput. 1997, 217(66):157-192.

    Article  MathSciNet  Google Scholar 

  4. Talay D, Vaillant O: A stochastic particle method with random weights for the computation of statistical solutions of McKean-Vlasov equations. Ann. Appl. Probab. 2003, 13(1):140-180.

    Article  MathSciNet  Google Scholar 

  5. Lasry JM, Lions PL: Mean field games. Jpn. J. Math. 2007, 2: 229-260. 10.1007/s11537-007-0657-8

    Article  MathSciNet  Google Scholar 

  6. Pardoux E, Peng SG: Backward doubly stochastic differential equations and systems of quasilinear SPDEs. Probab. Theory Relat. Fields 1994, 98: 209-227. 10.1007/BF01192514

    Article  MathSciNet  Google Scholar 

  7. Bally V, Matoussi A: Weak solutions for SPDEs and backward doubly stochastic differential equations. J. Theor. Probab. 2001, 14(1):125-164. 10.1023/A:1007825232513

    Article  MathSciNet  Google Scholar 

  8. Boufoussi B, Casteren J, Mrhardy N: Generalized backward doubly stochastic differential equations and SPDEs with nonlinear Neumann boundary conditions. Bernoulli 2007, 13(2):423-446. 10.3150/07-BEJ5092

    Article  MathSciNet  Google Scholar 

  9. Han YC, Peng SG, Wu Z: Maximum principle for backward doubly stochastic control systems with applications. SIAM J. Control Optim. 2010, 48(7):4224-4241. 10.1137/080743561

    Article  MathSciNet  Google Scholar 

  10. N’zi M, Owo J: Backward doubly stochastic differential equations with discontinuous coefficients. Stat. Probab. Lett. 2009, 79: 920-926. 10.1016/j.spl.2008.11.011

    Article  MathSciNet  Google Scholar 

  11. Shi YF, Gu YL, Liu K: Comparison theorems of backward doubly stochastic differential equations and applications. Stoch. Anal. Appl. 2005, 23: 97-110. 10.1081/SAP-200044444

    Article  MathSciNet  Google Scholar 

  12. Wu Z, Zhang F: BDSDEs with locally monotone coefficients and Sobolev solutions for SPDEs. J. Differ. Equ. 2011, 251: 759-784. 10.1016/j.jde.2011.05.017

    Article  Google Scholar 

  13. Zhang Q, Zhao HZ: Stationary solutions of SPDEs and infinite horizon BDSDEs. J. Funct. Anal. 2007, 252: 171-219. 10.1016/j.jfa.2007.06.019

    Article  MathSciNet  Google Scholar 

  14. Buckdahn R, Djehiche B, Li J, Peng SG: Mean-field backward stochastic differential equations: a limit approach. Ann. Probab. 2009, 37(4):1524-1565. 10.1214/08-AOP442

    Article  MathSciNet  Google Scholar 

  15. Buckdahn R, Li J, Peng SG: Mean-field backward stochastic differential equations and related partial differential equations. Stoch. Process. Appl. 2009, 119(10):3133-3154. 10.1016/j.spa.2009.05.002

    Article  MathSciNet  Google Scholar 

Download references

Acknowledgements

The author would like to thank the editor and anonymous referees for their constructive and insightful comments on improving the quality of this revision, and to thank professor Zhen Wu for many helpful suggestions.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Ruimin Xu.

Additional information

Competing interests

The author declares that he has no competing interests.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 2.0 International License (https://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Xu, R. Mean-field backward doubly stochastic differential equations and related SPDEs. Bound Value Probl 2012, 114 (2012). https://doi.org/10.1186/1687-2770-2012-114

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1687-2770-2012-114

Keywords