BSDE with quadratic growth and unbounded terminal value

Philippe Briand

Ying Hu

IRMAR, Université Rennes 1, 35 042 RENNES Cedex, FRANCE philippe.briand@univ-rennes1.fr ying.hu@univ-rennes1.fr

1. Introduction. In this paper we are concerned with real valued backward stochastic differential equations – BSDEs for short in the remaining – Y t = ξ + ∫ t T f ( s , Y s , Z s ) d s − ∫ t T Z s ⋅ d B s , 0 ≤ t ≤ T   where ( B t ) t ≥ 0   is a standard brownian motion. Such equations have been extensively studied since the first paper of E. Pardoux and S. Peng [PP90] . The full list of contributions is too long to give and we will only quote results in our framework.
Our setting is mainly the following : the generator, namely the function f   , is of quadratic growth in the variable z   and the terminal condition, the random variable ξ   , will not be bounded. BSDEs with quadratic growth have been first studied by Magdalena Kobylanski in her PhD (see [Kob97, Kob00] ) and then by Jean-Pierre Lepeltier and Jaime San Martin in [LSM98] . We should point out that BSDEs with quadratic growth in the variable z   have found applications in control and finance, see, e.g., Bismut [Bis78] , El Karoui, Rouge [EKR00] , Hu, Imkeller, Muller [HIM05] , . . .
All the results on BSDEs with quadratic growth require that the terminal condition ξ   is a bounded random variable. The boundedness of the terminal condition appears, from the point of view of the applications, to be restrictive and, moreover, from a theoretical point of view, is not necessary to obtain a solution. Indeed, let us consider the following well known equation :
Y t = ξ + 1 2 ∫ t T | Z s | 2 d s − ∫ t T Z s d B s , 0 ≤ t ≤ T ;   the change of variables P t = e Y t   , Q t = e Y t Z t   , leads to the equation Ph. Briand, Draft 25 Mars 2005 P t = e ξ − ∫ t T Q s d B s   which has a solution as soon as e ξ   is integrable.
On this simple example we see that the existence of exponential moments of the terminal condition is sufficient to construct a solution to our BSDE. Our paper will be focused on the theoretical study of these BSDEs but with unbounded terminal value with only exponential moments.
To fill the gap between boundedness and existence of exponential moments, we will use an approach based upon a localization procedure together with a priori bounds. Let us quickly explain how it works on a simple example. Let f : R × R d ⟶ R   be a continuous function and ξ   be a nonnegative terminal condition such that | f ( y , z ) | ≤ 1 2 | z | 2 , E [ e ξ ] < ∞ ,   and let us try to construct a solution to the BSDE Y t = ξ + ∫ t T f ( Y s , Z s ) d s − ∫ t T Z s ⋅ d B s , 0 ≤ t ≤ T .   As mentioned before, BSDEs with quadratic growth in the variable z   can be solved when the terminal solution is bounded. That is why we introduce ( Y n , Z n )   as the minimal solution to the BSDE Y t n = ξ ∧ n + ∫ t T f ( Y s n , Z s n ) d s − ∫ t T Z s n ⋅ d B s ,   and of course we want to pass to the limit when n → ∞   in this equation.
The process Y n   is known to be bounded but the estimate depends on ∥ ξ ∧ n ∥ ∞   and thus is far from being useful when ξ   is not bounded. The first step of our approach consists in finding an estimation for Y n   independent of n   . In this example, we can use the explicit formula mentioned before to show that 0 ≤ − ln E ( e − ( ξ ∧ n ) | ℱ t ) ≤ Y t n ≤ ln E ( e ξ ∧ n | ℱ t ) ≤ ln E ( e ξ | ℱ t ) .   With these inequalities in hands, we introduce the stopping time τ k = inf { t ∈ [ 0 , T ] : ln E ( e ξ | ℱ t ) ≥ k } ∧ T   and instead of working on the time interval [ 0 , T ]   we will restrict ourselves to [ 0 , τ k ]   by considering the BSDE Y t ∧ τ k n = Y τ k n + ∫ t ∧ τ k T ∧ τ k f ( Y s n , Z s n ) d s − ∫ t ∧ τ k T ∧ τ k Z s n ⋅ d B s , 0 ≤ t ≤ T .   By construction, we have sup n sup t ‖ Y t ∧ τ k n ‖ ∞ ≤ k   . This last property together with the fact that the sequence ( Y n ) n ≥ 1   is nondecreasing allows us, with the help of a result of Kobylanski, to pass to the limit when n → ∞   , k   being fixed and then to send k   to infinity to get a solution.
The rest of the paper is organized as follows. Next section is devoted to the notations we use during this text. In Section  3 , we claim our main result that we prove in Section  4 .
Section  5 is devoted to some additional results on BSDEs with quadratic growth in z   .
Finally, in the last section, we apply the same approach to study BSDEs with terminal value in L 1   .
2. Notations. For the remaining of the paper, let us fix a nonnegative real number T > 0   .
First of all, B = { B t } t ≥ 0   is a standard brownian motion with values in R d   defined on some complete probability space ( Ω , ℱ , P )   . { ℱ t } t ≥ 0   is the augmented natural filtration of B   which satisfies the usual conditions. In this paper, we will always use this filtration. P   denotes the sigma-field of predictable subsets of [ 0 , T ] × Ω   .
As mentioned in the introduction, a BSDE is an equation of the following type
Y t = ξ + ∫ t T f ( s , Y s , Z s ) d s − ∫ t T Z s ⋅ d B s , 0 ≤ t ≤ T . (1)
f   is called the generator and ξ   the terminal condition.
Let us recall that a generator is a random function f : [ 0 , T ] × Ω × R × R d ⟶ R   which is measurable with respect to P × ℬ ( R ) × ℬ ( R d )   and a terminal condition is simply a real ℱ T   –measurable random variable.
By a solution to the BSDE  1 we mean a pair ( Y , Z ) = { ( Y t , Z t ) } t ∈ [ 0 , T ]   of predictable processes with values in R × R d   such that P   –a.s., t ↦ Y t   is continuous, t ↦ Z t   belongs to L 2 ( 0 , T )   , t ↦ f ( t , Y t , Z t )   belongs to L 1 ( 0 , T )   and P   –a.s.
Y t = ξ + ∫ t T f ( s , Y s , Z s ) d s − ∫ t T Z s ⋅ d B s , 0 ≤ t ≤ T .   We will use the notation BSDE( ξ , f   ) to say that we consider the BSDE whose generator is f   and whose terminal condition is ξ   ; ( Y f ( ξ ) , Z f ( ξ ) )   means a solution to the BSDE( ξ , f   ).
( Y f ( ξ ) , Z f ( ξ ) )   is said to be minimal if P   -a.s., for each t ∈ [ 0 , T ]   , Y t f ( ξ ) ≤ Y t g ( ζ )   whenever P   –a.s. ξ ≤ ζ   and f ( t , y , z ) ≤ g ( t , y , z )   for all ( t , y , z )   . ( Y f ( ξ ) , Z f ( ξ ) )   is said to be minimal in some space ℬ   if it belongs to this space and the previous property holds true as soon as ( Y g ( ζ ) , Z g ( ζ ) ) ∈ ℬ   .
For any real p > 0   , S p   denotes the set of real-valued, adapted and càdlàg processes { Y t } t ∈ [ 0 , T ]   such that ∥ Y ∥ S p : = E [ sup t ∈ [ 0 , T ] | Y t | p ] 1 ∧ 1 / p < + ∞ .   If p ≥ 1   , ∥ ⋅ ∥ S p   is a norm on S p   and if p ∈ ( 0 , 1 )   , ( X , X ′ ) ↦ ∥ X − X ′ ∥ S p   defines a distance on S p   . Under this metric, S p   is complete.
M p   denotes the set of (equivalent classes of ) predictable processes { Z t } t ∈ [ 0 , T ]   with values in R d   such that ∥ Z ∥ M p : = E [ ( ∫ 0 T | Z s | 2 d s ) p / 2 ] 1 ∧ 1 / p < + ∞ .   For p ≥ 1   , M p ( R n )   is a Banach space endowed with this norm and for p ∈ ( 0 , 1 )   , M p   is a complete metric space with the resulting distance.
We set S = ∪ p > 1 S p   , M = ∪ p > 1 M p   and denote by S ∞   the set of predictable bounded processes. Finally, let us recall that a continuous process { Y t } t ∈ [ 0 , T ]   belongs to the class (D) if the family { Y τ : τ stopping time bounded by T }   is uniformly integrable.
3. Quadratic BSDEs. In this section, we consider BSDE( ξ , f   ) when the generator f   has a linear growth in y   and a quadratic growth in z   . We denote  H1 the assumption: there exist α ≥ 0   , β ≥ 0   and γ > 0   such that P   –a.s.
∀ t ∈ [ 0 , T ] , ( y , z ) ↦ f ( t , y , z ) is continuous , ∀ ( t , y , z ) ∈ [ 0 , T ] × R × R d , | f ( t , y , z ) | ≤ α + β | y | + γ 2 | z | 2 . (H1)
Concerning the terminal condition ξ   , we will assume that
E [ e γ e β T | ξ | ] < + ∞ . (H2)
We will use also a stronger assumption on the integrability of ξ   namely
∃ λ > γ e β T , E [ e λ | ξ | ] < + ∞ . (H3)
It is clear that we can assume without loss of generality that α ≥ β / γ   .
As we explained in the introduction, our method relies heavily on a priori estimate. To obtain such estimations, we will use the change of variable P t = e γ Y t   , Q t = γ e γ Y t Z t   ; if ( Y , Z )   is a solution to the BSDE( ξ , f   ), ( P , Q )   solves the BSDE P t = e γ ξ + ∫ t T F ( s , P s , Q s ) d s − ∫ t T Q s ⋅ d B s , 0 ≤ t ≤ T ,   with the function F   defined by
F ( s , p , q ) = 1 p > 0 ( γ p f ( s , ln p γ , q γ p ) − 1 2 | q | 2 p ) . (2)
In view of the growth of the generator f   , we have F ( s , p , q ) ≤ 1 p > 0 p ( α γ + β | ln p | )   . For notational convenience, we denote by H   the function ∀ p ∈ R , H ( p ) = p ( α γ + β ln p ) 1 [ 1 , + ∞ ) ( p ) + γ α 1 ( − ∞ , 1 ) ( p ) .   It is straightforward to check that, since α ≥ β / γ   , H   is convex and locally Lipschitz continuous and that, for any real p > 0   , p ( α γ + β | ln p | ) ≤ H ( p )   . Thus we deduce the inequality
∀ s ∈ [ 0 , T ] , ∀ p ∈ R , ∀ q ∈ R d , F ( s , p , q ) ≤ H ( p ) . (3)
To get an upper bound for Y t   , the idea is to compare more or less P t   with φ t ( ξ )   where, for any real z   , { φ t ( z ) } 0 ≤ t ≤ T   stands for the solution to the differential equation
φ t = e γ z + ∫ t T H ( φ s ) d s , 0 ≤ t ≤ T . (4)
Using the convexity of H   , we will able to prove that P t ≤ E ( φ t ( ξ ) | ℱ t ) , Y t ≤ 1 γ ln E ( φ t ( ξ ) | ℱ t ) .   Before proving this result rigorously, let us recall that the differential equation  4 can be solved easily. Indeed, we have, for any z ≥ 0   , φ t ( z ) = exp ( γ α e β ( T − t ) − 1 β ) exp ( z γ e β ( T − t ) ) , if β > 0 ,   and φ t ( z ) = e γ α ( T − t ) e γ z   if β = 0   .
Let us consider the case where z < 0   . If e γ z + T γ α ≤ 1   then the solution is φ t = e γ z + γ α ( T − t )   and otherwise there exists 0 < S < T   such that e γ z + γ α ( T − S ) = 1   and φ t = [ e γ z + γ α ( T − t ) ] 1 ( S , T ] ( t ) + exp ( γ α e β ( S − t ) − 1 β ) 1 [ 0 , S ] ( t ) .   It is plain to check that t ↦ φ t ( z )   is decreasing and that z ↦ φ t ( z )   is increasing and continuous.
Lemma 1. Let the assumption  H1 hold and let ξ   be a bounded ℱ T   –measurable random variable.
If ( Y , Z )   is a solution to the BSDE( ξ , f   ) in S ∞ × M 2   then − 1 γ ln E ( φ t ( − ξ ) | ℱ t ) ≤ Y t ≤ 1 γ ln E ( φ t ( ξ ) | ℱ t ) .  
  • Proof. Let us set Φ t = E ( φ t ( ξ ) | ℱ t )   . We have Φ t = E ( e γ ξ + ∫ t T H ( φ s ( ξ ) ) d s | ℱ t ) = E ( e γ ξ + ∫ t T E ( H ( φ s ( ξ ) ) | ℱ s ) d s | ℱ t ) .   Thus writing the bounded brownian martingale E ( e γ ξ + ∫ 0 T E ( H ( φ s ( ξ ) ) | ℱ s ) d s | ℱ t ) = E [ e γ ξ + ∫ 0 T E ( H ( φ s ( ξ ) ) | ℱ s ) d s ] + ∫ 0 t Ψ s ⋅ d B s   ( Φ , Ψ )   solves the BSDE Φ t = e γ ξ + ∫ t T E ( H ( φ s ( ξ ) ) | ℱ s ) d s − ∫ t T Ψ s ⋅ d B s .   On the other hand, if ( Y , Z ) ∈ S ∞ × M 2   is a solution of  1 , setting as before P t = e γ Y t   , Q t = γ e γ Y t Z t   , we have P t = e γ ξ + ∫ t T F ( s , P s , Q s ) d s − ∫ t T Q s ⋅ d B s ,   with F   defined by  2 .
    It follows that Φ t − P t = ∫ t T ( H ( Φ s ) − H ( P s ) ) d s + ∫ t T R s d s − ∫ t T ( Ψ s − Q s ) ⋅ d B s   where, in view of the inequality  3 and since H   is convex, R s = E ( H ( φ s ( ξ ) ) | ℱ s ) − H ( E ( φ s ( ξ ) | ℱ s ) ) + H ( P s ) − F ( s , P s , Q s )   is a nonnegative process.
    H   is only locally Lipschitz but since Φ   and P   are bounded we can apply the comparison theorem to get P t ≤ Φ t   and Y t ≤ 1 γ ln Φ t   .
    Finally, since the function − f ( t , − y , − z )   still satisfies the assumption  H1 , we get also the inequality − Y t ≤ 1 γ ln E ( φ t ( − ξ ) | ℱ t )   .
We are now in position to prove that under the assumptions described before the BSDE  1 has at least a solution.
Theorem 2. Let the assumptions  H1 and  H2 hold. Then the BSDE  1 has at least a solution ( Y , Z )   such that :
− 1 γ ln E ( φ t ( − ξ ) | ℱ t ) ≤ Y t ≤ 1 γ ln E ( φ t ( ξ ) | ℱ t ) . (5)
If moreover,  H3 holds, then Z   belongs to M 2   .
  • Proof of the last part of Theorem  2 . If ( Y , Z )   is a solution to the BSDE  1 such that the inequalities  5 hold, then | Y t | ≤ 1 γ ln E ( φ 0 ( | ξ | ) | ℱ t )   and, under the assumption  H3 , we deduce that, for some p > 1   , E [ sup t ∈ [ 0 , T ] e p γ | Y t | ] < + ∞ .   For n ≥ 1   , let τ n   be the following stopping time τ n = inf { t ≥ 0 : ∫ 0 t e 2 γ | Y s | | Z s | 2 d s ≥ n } ∧ T ,   and let us consider the function from R +   into itself defined by u ( x ) = 1 γ 2 ( e γ x − 1 − γ x ) .   x ↦ u ( | x | )   is C 2   and we have from Itô's formula, with the notation sgn ( x ) = − 1 x ≤ 0 + 1 x > 0   ,
    u ( | Y 0 | ) = u ( | Y t ∧ τ n | ) + ∫ 0 t ∧ τ n ( u ′ ( | Y s | ) sgn ( Y s ) f ( s , Y s , Z s ) − 1 2 u ′ ′ ( | Y s | ) | Z s | 2 ) d s
    − ∫ 0 t ∧ τ n u ′ ( | Y s | ) sgn ( Y s ) Z s ⋅ d B s .
    It follows from  H1 since u ′ ( x ) ≥ 0   for x ≥ 0   that
    u ( | Y 0 | ) ≤ u ( | Y t ∧ τ n | ) + ∫ 0 t ∧ τ n u ′ ( | Y s | ) ( α + β | Y s | ) d s − ∫ 0 t ∧ τ n u ′ ( | Y s | ) sgn ( Y s ) Z s ⋅ d B s
    − 1 2 ∫ 0 t ∧ τ n ( u ′ ′ ( | Y s | ) − γ u ′ ( | Y s | ) ) | Z s | 2 d s .
    Moreover, we have ( u ′ ′ − γ u ′ ) ( x ) = 1   for x ≥ 0   and, taking expectation of the previous inequality, we get 1 2 E [ ∫ 0 T ∧ τ n | Z s | 2 d s ] ≤ E [ 1 γ 2 sup t ∈ [ 0 , T ] e γ | Y t | + 1 γ ∫ 0 T e γ | Y s | ( α + β | Y s | ) d s ]   Fatou's lemma together with the fact that e γ | Y t | ∈ S p   gives the result.
4. Proof of Theorem  2 .
Let us first construct a solution to the BSDE  1 in the case where ξ   is nonnegative.
For each n ∈ N *   , we set ξ n = ξ ∧ n   . Then it is known from [Kob00,Theorem2.3] that the BSDE Y t n = ξ n + ∫ t T f ( s , Y s n , Z s n ) d s − ∫ t T Z s n ⋅ d B s , 0 ≤ t ≤ T   has a minimal solution ( Y n , Z n )   in S ∞ × M 2   . Lemma  1 implies the inequalities − 1 γ ln E ( φ t ( − ξ n ) | ℱ t ) ≤ Y t n ≤ 1 γ ln E ( φ t ( ξ n ) | ℱ t ) .   Since we consider only minimal solutions, we have, ∀ t ∈ [ 0 , T ] , Y t n ≤ Y t n + 1 .   We define Y = sup n ≥ 1 Y n   . Since 0 ≤ φ t ( ξ n ) ≤ φ 0 ( | ξ | )   and 0 ≤ φ t ( − ξ n ) ≤ φ 0 ( | ξ | )   , we deduce from the dominated convergence theorem, noting that the random variable φ 0 ( | ξ | )   is integrable by  H 2     , that − 1 γ ln E ( φ t ( − ξ ) | ℱ t ) ≤ Y t ≤ 1 γ ln E ( φ t ( ξ ) | ℱ t ) .   In particular, we have lim t → T Y t = ξ = Y T   . Indeed, for each S < T   , limsup t → T Y t ≤ limsup t → T 1 γ ln E ( φ t ( ξ ) | ℱ t ) ≤ lim t → T 1 γ ln E ( φ S ( ξ ) | ℱ t ) = 1 γ ln φ S ( ξ ) ,   and lim S → T 1 γ ln φ S ( ξ ) = ξ   . We can do the same for liminf   .
Let us introduce the following stopping time :
τ k = inf { t ∈ [ 0 , T ] : 1 γ ln E ( φ 0 ( | ξ | ) | ℱ t ) ≥ k } ∧ T .   Then ( Y k n , Z k n ) : = ( Y t ∧ τ k n , Z t n 1 t ≤ τ k )   satisfies the following BSDE Y k n = ξ k n + ∫ t T 1 s ≤ τ k f ( s , Y k n ( s ) , Z k n ( s ) ) d s − ∫ t T Z k n ( s ) ⋅ d B s ,   where of course ξ k n = Y k n ( T ) = Y τ k n   .
We are going to pass to the limit when n   tends to + ∞   for k   fixed in this last equation.
The key point is that Y k n   is increasing in n   and remains bounded by k   . At this stage, let us mention a mere generalization of Proposition 2.4 in [Kob00] .
Lemma 3 ([Kob00] ). Let ( ξ n ) n ≥ 1   be a sequence of ℱ T   –measurable bounded random variables and ( f n ) n ≥ 1   be a sequence of generators which are continuous with respect to ( y , z )   .
We assume that ( ξ n ) n ≥ 1   converges P   –a.s. to ξ   , that ( f n ) n ≥ 1   converges locally uniformly in ( y , z )   to the generator f   , and also that
  • 1. sup n ≥ 1 ∥ ξ n ∥ ∞ < + ∞   ;
  • 2. sup n ≥ 1 | f n ( t , y , z ) |   satisfies the inequality in  H1 .
If for each n ≥ 1   , the BSDE( ξ n , f n   ) has a solution in S ∞ × M 2   , such that ( Y f n ( ξ n ) ) n ≥ 1   is nondecreasing (respectively nonincreasing), then P   –a.s. ( Y t f n ( ξ n ) ) n ≥ 1   converges uniformly on [ 0 , T ]   to Y t = sup n ≥ 1 Y t f n ( ξ n )   (respectively Y t = inf n ≥ 1 Y t f n ( ξ n )   ), ( Z f n ( ξ n ) ) n ≥ 1   converges to some Z   in M 2   and ( Y , Z )   is a solution to BSDE( ξ , f   ) in S ∞ × M 2   .
  • Proof. It follows from Lemma  1 that there exists r > 0   such that, P   –a.s.
    ∀ n ≥ 1 , ∀ t ∈ [ 0 , T ] , | Y t f n ( ξ n ) | ≤ r .   Let us consider the continuous function ρ ( x ) = x r / max ( r , | x | )   . Since ρ ( x ) = x   for | x | ≤ r   , ( Y f n ( ξ n ) , Z f n ( ξ n ) )   solves the BSDE( ξ n , g n   ) where g n ( t , y , z ) = f n ( t , ρ ( y ) , z )   . Obviously, we have, for each n ≥ 1   , | g n ( t , y , z ) | ≤ α + β r + γ 2 | z | 2 ,   and thus we can apply the result of Kobylanski.
Setting Y k ( t ) = sup n Y k n ( t )   , it follows from the previous lemma that there exists a process Z k ∈ M 2   such that lim n Z k n = Z k   in M 2   and ( Y k , Z k )   solves the BSDE
Y k ( t ) = ξ k + ∫ t T 1 s ≤ τ k f ( s , Y k ( s ) , Z k ( s ) ) d s − ∫ t T Z k ( s ) ⋅ d B s , (6)
where ξ k = sup n Y τ k n   .
But τ k ≤ τ k + 1   , and thus we get, coming back to the definition of Y k   , Z k   and Y   , Y t ∧ τ k = Y k + 1 ( t ∧ τ k ) = Y k ( t ) , Z k + 1 ( t ) 1 t ≤ τ k = Z k ( t ) .   As τ k → T   and the Y k   's are continuous processes we deduce in particular that Y   is continuous on [ 0 , T )   . On the other hand, as mentioned before lim t → T Y t = ξ   and Y T   is equal to ξ   by construction. Thus Y   is a continuous process on the closed interval [ 0 , T ]   .
Then we define Z   on ( 0 , T )   by setting :
Z t = Z k ( t ) , if t ∈ ( 0 , τ k ) .   From  6 , ( Y , Z )   satisfies:
Y t ∧ τ k = Y τ k + ∫ t ∧ τ k τ k f ( s , Y s , Z s ) d s − ∫ t ∧ τ k τ k Z s ⋅ d B s . (7)
Finally, we have
P ( ∫ 0 T | Z s | 2 d s = ∞ ) = P ( ∫ 0 T | Z s | 2 d s = ∞ , τ k = T ) + P ( ∫ 0 T | Z s | 2 d s = ∞ , τ k < T )
≤ P ( ∫ 0 τ k | Z k ( s ) | 2 d s = ∞ ) + P ( τ k < T ) ,
and we deduce that, P   –a.s.
∫ 0 T | Z s | 2 d s < ∞ .   By sending k   to infinity in  7 , we deduce that ( Y , Z )   is a solution of  1 .
Let us explain quickly how to extend this construction to the general case. Let us fix n ∈ N *   and p ∈ N *   and set ξ n , p = ξ + ∧ n − ξ − ∧ p   . Let us consider, ( Y n , p , Z n , p )   the minimal bounded solution to the BSDE Y t n , p = ξ n , p + ∫ t T f ( s , Y s n , p , Z s n , p ) d s − ∫ t T Z s n , p ⋅ d B s , 0 ≤ t ≤ T   which satisfies − 1 γ ln E ( φ t ( − ξ n , p ) | ℱ t ) ≤ Y t n , p ≤ 1 γ ln E ( φ t ( ξ n , p ) | ℱ t ) .   We have, ∀ t ∈ [ 0 , T ] , Y t n , p + 1 ≤ Y t n , p ≤ Y t n + 1 , p ,   and we define Y p = sup n ≥ 1 Y n , p   so that Y t p + 1 ≤ Y t p   and Y t = inf p ≥ 1 Y t p   .
By the dominated convergence theorem, we have − 1 γ ln E ( φ t ( − ξ ) | ℱ t ) ≤ Y t ≤ 1 γ ln E ( φ t ( ξ ) | ℱ t ) ,   and in particular, we have lim t → T Y t = ξ = Y T   .
( Y t ∧ τ k n , p , Z t n , p 1 t ≤ τ k )   solves the BSDE Y t ∧ τ k n , p = Y τ k n , p + ∫ t T 1 s ≤ τ k f ( s , Y s n , p , Z s n , p ) d s − ∫ t T Z s n , p 1 s ≤ τ k ⋅ d B s .   But, once again Y t ∧ τ k n , p   is increasing in n   and decreasing in p   and remains bounded by k   . Arguing as before, setting Y k ( t ) = inf p sup n Y t ∧ τ k n , p   , there exists a process Z k   such that lim p lim n Z n , p ( s ) 1 s ≤ τ k = Z k ( s )   and ( Y k , Z k )   still solves the BSDE  6 . The rest of the proof is unchanged.
5. Additional results on quadratic BSDEs. 5.1. Minimal solution. In this section, we give some complements on BSDEs with quadratic growth in z   .
Proposition 4. Let  H1 holds and assume moreover that there exists an integer r ≥ 0   such that P   –a.s.
f ( t , y , z ) ≥ − r ( 1 + | y | + | z | ) .   Let us assume also that  H3 holds for ξ +   and that, for some p > 1   , ξ − ∈ L p   .
Then BSDE( ξ , f   ) has a minimal solution in S   .
  • Proof. For each n ≥ r   , let us consider the function f n ( t , y , z ) = inf { f ( t , p , q ) + n | p − y | + n | q − z | : ( p , q ) ∈ Q 1 + d } .   Then f n   is well defined and it is globally Lipschitz continuous with constant n   . Moreover ( f n ) n ≥ r   is increasing and converges pointwise to f   . Dini's theorem implies that the convergence is also uniform on compact sets. We have also, for all n ≥ r   , − r ( 1 + | y | + | z | ) ≤ f n ( t , y , z ) ≤ f ( t , y , z )   Let ( Y n , Z n )   be the unique solution in S p × M p   to BSDE( ξ , f n   ). It follows from the classical comparison theorem that Y t r ≤ Y t n ≤ Y t n + 1 .   Let us prove that Y t n ≤ 1 γ ln E ( φ t ( ξ ) | ℱ t )   . To do this let us recall that, since f n   is Lipschitz, Y t n = lim m → + ∞ Y t f n ( ξ m )   where ξ m = ξ 1 | ξ | ≤ m   . Moreover E ( φ t ( ξ m ) | ℱ t ) ⟶ E ( φ t ( ξ ) | ℱ t )   a.s.
    since sup m ≥ 1 | φ t ( ξ m ) | ≤ φ 0 ( ξ + )   which is integrable. Thus we have only to prove that Y t f n ( ξ m ) ≤ 1 γ ln E ( φ t ( ξ m ) | ℱ t )   . We keep the notations of the beginning of Section  4 .
    ( Φ , Ψ )   is solution to the BSDE Φ t = e γ ξ m + ∫ t T H ( Φ s ) d s + ∫ t T Γ s d s − ∫ t T Ψ s ⋅ d B s ,   where Γ s = E ( H ( φ s ( ξ ) ) | ℱ s ) − H ( Φ s )   is a nonnegative process since H   is convex.
    It follows by setting U t = 1 γ ln Φ t   , V t = Ψ t γ Φ t   that ( U , V )   solves the BSDE U t = ξ m + ∫ t T g ( s , U s , V s ) d s − ∫ t T V s ⋅ d B s   where we have set g ( s , u , v ) = ( α + β u ) 1 u ≥ 0 + α e γ | u | 1 u < 0 + γ 2 | v | 2 + C s   with C s = 1 γ e − γ U s Γ s   .
    Since the process C   is still nonnegative, we have the inequalities f n ( t , u , v ) ≤ f ( t , u , v ) ≤ g ( t , u , v )   taking into account the fact that α γ ≥ β   .
    Since f n   is Lipschitz continuous and ( Y f n ( ξ m ) − U ) +   belongs to S   , we can apply the extended comparison theorem (see Proposition  5 ) to get, for each m ≥ 1   , Y t f n ( ξ m ) ≤ U t   and thus the inequality we want to obtain.
    We set Y = sup n ≥ r Y n   and, for k ≥ 1   , τ k = inf { t ∈ [ 0 , T ] : max ( 1 γ ln E ( φ t ( ξ ) | ℱ t ) , − Y t r ) ≥ k } ∧ T .   Arguing as in the proof of Theorem  2 , we construct a process Z   such that ( Y , Z )   solves BSDE( ξ , f   ).
    Let us show that this solution is minimal in S   . Let ( Y ′ , Z ′ )   be a solution to the BSDE( ξ ′ , f ′   ) where ξ ≤ ξ ′   and f ≤ f ′   . It is enough to check that Y n ≤ Y ′   to prove that Y ≤ Y ′   . But this is a direct consequence of Proposition  5 .
To be complete, let us claim and prove the extended comparison theorem that we used in the proof of the previous result.
Proposition 5. Let ( Y , Z )   be a solution to BSDE( ξ , f   ) and ( Y ′ , Z ′ )   be a solution to BSDE( ξ ′ , f ′   ). We assume that ξ ≤ ξ ′   and that f   satisfies, for some constants μ   and λ   , P   –a.s.
( y − y ′ ) ⋅ ( f ( t , y , z ) − f ( t , y ′ , z ) ) ≤ μ | y − y ′ | 2 ;
| f ( t , y , z ) − f ( t , y , z ′ ) | ≤ λ | z − z ′ | ;
If ( Y − Y ′ ) +   belongs to S   , then P   –a.s. Y t ≤ Y t ′   .
  • Proof. Let us fix n ∈ N *   and denote τ n   the stopping time τ n = inf { t ∈ [ 0 , T ] : ∫ 0 t ( | Z s | 2 + | Z s ′ | 2 ) d s ≥ n } ∧ T .   Tanaka's formula leads to the equation, setting U t = Y t − Y t ′   , V t = Z t − Z t ′   ,
    e μ ( t ∧ τ n ) U t ∧ τ n + ≤ e μ τ n U τ n + − ∫ t ∧ τ n τ n e μ s 1 U s > 0 V s ⋅ d B s + ∫ t ∧ τ n τ n e μ s { 1 U s > 0 ( f ( s , Y s , Z s ) − f ′ ( s , Y s ′ , Z s ′ ) ) − μ U s + } d s . (8)
    First of all, we write f ( s , Y s , Z s ) − f ′ ( s , Y s ′ , Z s ′ ) = f ( s , Y s , Z s ) − f ( s , Y s ′ , Z s ) + f ( s , Y s ′ , Z s ) − f ′ ( s , Y s ′ , Z s ′ )   and we deduce, using the monotonicity of f   in y   that 1 U s > 0 ( f ( s , Y s , Z s ) − f ′ ( s , Y s ′ , Z s ′ ) ) − μ U s + ≤ 1 U s > 0 ( f ( s , Y s ′ , Z s ) − f ′ ( s , Y s ′ , Z s ′ ) ) .   But f ( s , Y s ′ , Z s ′ ) − f ′ ( s , Y s ′ , Z s ′ )   is nonpositive so that 1 U s > 0 ( f ( s , Y s , Z s ) − f ′ ( s , Y s ′ , Z s ′ ) ) − μ U s + ≤ 1 U s > 0 ( f ( s , Y s ′ , Z s ) − f ( s , Y s ′ , Z s ′ ) ) .   Finally we set β s = ( f ( s , Y s ′ , Z s ) − f ( s , Y s ′ , Z s ′ ) ) V s | V s | 2   which is a process bounded by λ   .
    Coming back to  8 , we obtain the following inequality e μ ( t ∧ τ n ) U t ∧ τ n + ≤ e μ τ n U τ n + + ∫ t ∧ τ n τ n e μ s 1 U s > 0 β s ⋅ V s d s − ∫ t ∧ τ n τ n e μ s 1 U s > 0 V s ⋅ d B s   By Girsanov's theorem, we deduce that E * [ e μ ( t ∧ τ n ) U t ∧ τ n + ] ≤ E * [ e μ τ n U τ n + ] ,   where P *   is the probability measure on ( Ω , ℱ T )   whose density with respect to P   is D T = exp { ∫ 0 T β s ⋅ d B s − 1 2 ∫ 0 T | β s | 2 d s } ;   it is worth noting that, since β   is a bounded process, D T   has moments of all order.
    Since we know that U +   belongs to S   , we can easily send n   to infinity to get E * [ e μ t U t + ] ≤ 0 .   Thus U t ≤ 0   P *   –a.s. and since P *   is equivalent to P   on ( Ω , ℱ T )   , Y t ≤ Y t ′   P   –a.s..
5.2. One extension. In this paragraph, we explain how we can extend our results to a more general setting allowing a superlinear growth of the generator in the variable y   as in the work [LSM98] .
Let h : R + ⟶ R +   be a nondecreasing convex C 1   function with h ( 0 ) > 0   such that ∫ 0 + ∞ d u h ( u ) = + ∞ .   We denote  H1' the assumption: there exists γ > 0   such that P   –a.s.
∀ t ∈ [ 0 , T ] , ( y , z ) ↦ f ( t , y , z ) is continuous , ∀ ( t , y , z ) ∈ [ 0 , T ] × R × R d , | f ( t , y , z ) | ≤ h ( | y | ) + γ 2 | z | 2 , sup y > 0 e − γ y h ( y ) < + ∞ . (H1')
Let us point out that the previous setting, namely the linear growth condition, corresponds to h ( y ) = α + β y   but we can also have a superlinear growth in y   ; for instance, we can take h ( y ) = α ( y + e ) ln ( y + e )   .
Before giving our integrability condition for the terminal value ξ   , let us explain what is the first modification we have to do. We consider only the case where h   is not constant.
According to the third point of  H1' , let us denote by c = sup p ∈ ( 0 , 1 ) γ p h ( − ln p γ )   and let us define p 0 = inf { p ≥ 1 : γ p h ( ln p γ ) ≥ c } .   We define finally H ( p ) = γ p h ( ln p γ ) 1 p ≥ p 0 + c 1 p < p 0 .   Then H   is convex and we have the following result.
Lemma 6. Let z ∈ R   . The differential equation φ t = e γ z + ∫ t T H ( φ s ) d s , 0 ≤ t ≤ T ,   has a unique continuous solution { φ t ( z ) } 0 ≤ t ≤ T   which is decreasing. Moreover, for each t ∈ [ 0 , T ]   , the map z ↦ φ t ( z )   is increasing and continuous.
  • Proof. φ t   is solution if and only if u t = ln φ t / γ   is a solution of the differential equation u t ′ = − θ ( u t ) , 0 ≤ t ≤ T , u T = z ≥ 0 ,   where θ ( x ) = h ( x ) 1 x ≥ ln p 0 γ + c γ e − γ x 1 x < ln p 0 γ   . Let us consider the function Θ   defined by Θ ( x ) = ∫ − ∞ x 1 θ ( u ) d u , x ∈ R .   Since θ   is positive, Θ   is an increasing bijection from R   onto ( 0 , ∞ )   of class C 1   . It's plain to check that the unique solution to the previous differential equation is Θ − 1 ( T − t + Θ ( z ) )   since for any solution we have Θ ( u t ) ′ = − 1   . Thus φ t = e γ Θ − 1 ( T − t + Θ ( z ) )   and the proof of the lemma is complete.
We are now in position to give our second assumption.
φ 0 ( | ξ | ) is integrable. (H2')
Exactly as in the linear case, we can prove the following existence result.
Theorem 7. Let assumptions  H1' and  H2' hold. Then the BSDE  1 has at least a solution ( Y , Z )   such that :
− 1 γ ln E ( φ t ( ξ ) | ℱ t ) ≤ Y t ≤ 1 γ ln E ( φ t ( ξ ) | ℱ t ) .  
6. BSDEs in L 1   .
In this section, we use the method developed before to construct solutions to BSDEs when the data are only integrable. BSDEs with integrable data have been studied in [BDH + 03] and we show that, in the one dimensional case, we can extend the result quoted before.
Let us recall the framework of [BDH + 03] : assumption (A) holds true for the random function f   if there exist constants μ ∈ R   , λ ≥ 0   , δ ≥ 0   and α ∈ ( 0 , 1 )   such that
( y − y ′ ) ⋅ ( f ( t , y , z ) − f ( t , y ′ , z ) ) ≤ μ | y − y ′ | 2 ;
| f ( t , y , z ) − f ( t , y , z ′ ) | ≤ λ | z − z ′ | ;
y ⟶ f ( t , y , z ) is continuous ;
for each r > 0 , ψ r ( t ) : = sup | y | ≤ r | f ( t , y , 0 ) − f ( t , 0 , 0 ) | ∈ L 1 ( ( 0 , T ) × Ω ) ;
| f ( t , y , z ) − f ( t , y , 0 ) | ≤ δ ( g t + | y | + | z | ) α ,
where moreover the progressively measurable processes ( | f ( t , 0 , 0 ) | )   and ( g t )   and the terminal condition ξ   satisfy E [ | ξ | + ∫ 0 T ( | f ( s , 0 , 0 ) | + g s ) d s ] < + ∞ .   Let us recall the following result.
Lemma 8 ([BDH + 03] ). Under the assumption (A), BSDE  1 has a unique solution ( Y , Z )   such that Y   is of class (D) and Z ∈ M β   for some β > α   . Moreover ( Y , Z ) ∈ S β × M β   for each β ∈ ( 0 , 1 )   .
We should point out that the result holds true in all dimension not only in the real case.
The last assumption on the generator f   does not seem to be very natural; it would be better to have a condition of the type | f ( t , y , z ) | ≤ c ( 1 + | y | + | z | α )   and the remaining of this section is devoted to the construction of a solution under this assumption. But before, we state a comparison result for BSDEs under assumption (A).
Lemma 9. Let (A) holds. Then, if ξ ≤ ξ ′   belong to L 1   , Y t f ( ξ ) ≤ Y t f ( ξ ′ )   .
  • Proof. Let us fix n ∈ N *   and denote τ n   the stopping time τ n = inf { t ∈ [ 0 , T ] : ∫ 0 t ( | Z s | 2 + | Z s ′ | 2 ) d s ≥ n } ∧ T .   Tanaka's formula leads to the equation, setting U t = Y t − Y t ′   , V t = Z t − Z t ′   ,
    e μ ( t ∧ τ n ) U t ∧ τ n + ≤ e μ τ n U τ n + − ∫ t ∧ τ n τ n e μ s 1 U s > 0 V s ⋅ d B s + ∫ t ∧ τ n τ n e μ s { 1 U s > 0 ( f ( s , Y s , Z s ) − f ( s , Y s ′ , Z s ′ ) ) − μ U s + } d s .
    We deduce from the previous inequality, using the monotonicity of f   in y   that 1 U s > 0 ( f ( s , Y s , Z s ) − f ( s , Y s ′ , Z s ′ ) ) − μ U s + ≤ 1 U s > 0 ( f ( s , Y s ′ , Z s ) − f ( s , Y s ′ , Z s ′ ) )   and taking into account the last condition on f   , the right hand side of the previous inequality is bounded from above by X s : = 2 δ ( g s + | Y s ′ | + | Z s | + | Z s ′ | ) α .   It follows that e μ ( t ∧ τ n ) U t ∧ τ n + ≤ e μ τ n U τ n + + ∫ 0 T e μ s X s d s − ∫ t ∧ τ n τ n e μ s 1 U s > 0 V s ⋅ d B s   and thus that e μ ( t ∧ τ n ) U t ∧ τ n + ≤ E ( e μ τ n U τ n + + ∫ 0 T e μ s X s d s | ℱ t ) .   Since Y   and Y ′   belongs to the class (D), we can send n   to ∞   in the previous inequality (see [BDH + 03] for details) to get e μ t U t + ≤ E ( ∫ 0 T e μ s X s d s | ℱ t ) .   As a byproduct, we deduce that U +   belongs to S p   as soon as α p < 1   . Thus we can choose p > 1   such that α p < 1   .
    Since U +   belongs to S   , we can apply Proposition  5 to conclude the proof.
From now on we assume that f   is continuous and satisfies, for some constants c ≥ 0   and α ∈ ( 0 , 1 )   ,
| f ( t , y , z ) | ≤ c ( 1 + | y | + | z | α ) . (H4)
Moreover we will suppose that ξ ∈ L 1   . This last assumption is denoted by (H5).
Theorem 10. Let  H4 and (H5) hold. Then the BSDE  1 has a solution ( Y , Z )   such that Y   belongs to the class (D). Moreover ( Y , Z )   belongs to S β × M β   for all 0 < β < 1   .
  • Proof. To prove this result, we use the same approach as in the case of quadratic generators.
    So let us fix n ∈ N *   and p ∈ N *   and set ξ n , p = ξ + ∧ n − ξ − ∧ p   . Since f   is assumed to be a continuous map, we can consider, according to [LSM97] , ( Y f ( ξ n , p ) , Z f ( ξ n , p ) )   as the minimal solution to the BSDE Y t = ξ n , p + ∫ t T f ( s , Y s , Z s ) d s − ∫ t T Z s ⋅ d B s , 0 ≤ t ≤ T .   Since we are dealing with minimal solutions, we have Y t f ( ξ n , p + 1 ) ≤ Y t f ( ξ n , p ) ≤ Y t f ( ξ n + 1 , p ) ,   and we set Y t = inf p ≥ 1 sup n ≥ 1 Y t f ( ξ n , p )   .
    In order to apply the method described before, we have to find an upper bound independent of ( n , p )   for | Y f ( ξ n , p ) |   . For this let us observe that | f ( t , y , z ) | ≤ g ( y , z ) : = 2 c ( 1 + | y | + | z | α ∧ | z | ) .   This function g   is globally Lipschitz continuous so that we have from the classical comparison theorem Y f ( ξ n , p ) ≤ Y g ( ξ n , p )   and since we have the same inequality for − f ( t , − y , − z )   , − Y f ( ξ n , p ) ≤ Y g ( − ξ n , p ) .   But the function g   satisfies also the assumption (A) and thus form the comparison theorem in the integrable framework – Lemma  9 –, we deduce that | Y t f ( ξ n , p ) | ≤ Y t g ( | ξ | ) .   Let us define Y t = inf p ≥ 1 sup n ≥ 1 Y t f ( ξ n , p )   and for each k ≥ 1   , τ k = inf { t ∈ [ 0 , T ] : Y t g ( | ξ | ) ≥ k } ∧ T .   Exactly as in the proof of Theorem  2 , we construct a process Z   such that ( Y , Z )   solves the BSDE  1 .
    To conclude the proof, let us observe that since | Y t | ≤ Y t g ( | ξ | )   , Y   belongs to the class (D) and to S β   for each β ∈ ( 0 , 1 )   . It follows from [BDH + 03,Lemma3.1] that Z   belongs to M β   for β ∈ ( 0 , 1 )   .
References.

  1. Ph. Briand, B. Delyon, Y. Hu, É. Pardoux, and L. Stoica, L p   solutions of backward stochastic differential equations, Stochastic Process. Appl. 108 (2003), 109–129.
  2. J.-M. Bismut, Contrôle des systèmes linéaires quadratiques: applications de l'intégrale stochastique, Séminaire de Probabilités, XII (Univ. Strasbourg, Strasbourg, 1976/1977), Lecture Notes in Math., vol. 649, Springer, Berlin, 1978, pp. 180–264.
  3. N. El Karoui and R. Rouge, Pricing via utility maximization and entropy, Math. Finance 10 (2000), no. 2, 259–276, INFORMS Applied Probability Conference (Ulm, 1999).
  4. Y. Hu, P. Imkeller, and M. Müller, Utility maximization in incomplete markets, Ann. Appl. Probab. (2005), To appear.
  5. M. Kobylanski, Résultats d'existence et d'unicité pour des équations différentielles stochastiques rétrogrades avec des générateurs à croissance quadratique, C. R. Acad. Sci. Paris Sér. I Math. 324 (1997), no. 1, 81–86.
  6. , Backward stochastic differential equations and partial differential equations with quadratic growth, Ann. Probab. 28 (2000), no. 2, 558–602.
  7. J.-P. Lepeltier and J. San Martin, Backward stochastic differential equations with continuous coefficients, Statist. Probab. Lett. 32 (1997), no. 4, 425–430.
  8. , Existence for BSDE with superlinear-quadratic coefficient, Stochastics Stochastics Rep. 63 (1998), no. 3-4, 227–240.
  9. E. Pardoux and S. Peng, Adapted solution of a backward stochastic differential equation, Systems Control Lett. 14 (1990), no. 1, 55–61.