arxiv: v2 [math.pr] 19 Aug 2016

Transcription

1 Strong convergence of the ymmetrized Miltein cheme for ome CV-like SD Mireille Boy and Héctor Olivero arxiv:58458v mathpr 9 Aug 6 TOSCA Laboratory, INRIA Sophia Antipoli Méditerranée, France Departamento de Ingeniería Matemática, Univeridad de Chile, Chile March, 8 Abtract In thi paper we tudy the rate of convergence of a ymmetrized verion of the Miltein cheme applied to the olution of the one dimenional SD X t = x + bx d + σ X α dw, x >, σ >, α, Auming b/σ big enough, and b mooth, we prove a trong rate of convergence of order one, recovering the claical reult of Miltein for SD with mooth diffuion coefficient In contrat with other recent reult, our proof doe not relie on Lamperti tranformation, and it can be applied to a wide cla of drift function On the downide, our hypothei on the critical parameter value b/σ i more retrictive than other available in the literature Some numerical experiment and comparion with variou other cheme complement our theoretical analyi that alo applie for the imple projected Miltein cheme with ame convergence rate Introduction and main reult The Miltein cheme wa introduced by Miltein in 6 for one dimenional Stochatic Differential quation SD having mooth diffuion coefficient Introducing an appropriated correction term, thi cheme ha better convergence rate for the trong error than the claical uler-maruyama cheme Typically, when the drift and diffuion coefficient of one dimenional SD are twice continuouly differentiable with bounded derivative, the Miltein cheme i of order one for trong error ee eg Talay intead of one-half for the uler-maruyama cheme Thi well-know fact produce remark on blog, internet forum, and oftware package that ometime recommend to ue the Miltein cheme for contant elaticity of variance CV model in finance, or it extenion with tochatic volatility a SABR model, ee eg Delbaen and Shirakawa 8 and Lion and Muiela 5 for a dicuion on the weak exitence of uch model; CV are popular tochatic volatility model of the form dx t = µx t dt + σx γ t dw t with < γ < But the intereting fact in thi tory i that the rate of convergence of the Miltein cheme, for uch family of procee with < γ < i not yet well tudied, to the bet of our knowledge mireilleboy@inriafr holivero@dimuchilecl Financially upported by: Proyecto Meceup UCH67, the Dirección de Potgrado y Potítulo de la Vicerrectoría de Aunto Académico de la Univeridad de Chile, the Intituto Francé de Chile - mbajada de Francia en Chile, and the Center for Mathematical Modeling CMM

2 In thi paper we etablih a rate of convergence reult for a ymmetrized verion of the Miltein cheme applied to the olution of the one dimenional SD X t = x + bx d + σ X α dw, where x >, σ > and α < Of coure quation doe not atifie the hypothei to apply the claical reult of Miltein 6 In particular, the diffuion coefficient i only Hölder continuou wherea the claical hypothei i to have a C diffuion coefficient The main picture of our convergence rate reult i that Miltein cheme tay of order one in the cae of quation, but ome attention mut be paid to the value of b, α and σ There exit in the literature other trategie for the dicretization of the olution to There are ome reult baed on the Lamperti tranformation of the equation, for example, by Alfoni,, and by Chaagneux, Jacquier and Mihaylov 6 And alo, there ome are reult where the equation i dicretized directly, a in Berkaoui, Boy and Diop 3 or in Kahl and Jackël In the numerical experiment ection, we compare the ymmetrized Miltein cheme with a election of cheme propoed in the aforementioned reference We alo experiment the ymmetrized Miltein cheme in a multilevel Monte Carlo application and we compare with other cheme In the whole paper, we work under the following bai-hypothei: Hypothei The power parameter α in the diffuion coefficient of quation belong to, The drift coefficient b i Lipchitz with contant K >, and i uch that b > Hypothei i a claical aumption to enure a unique trong olution valued in R + We aume it in all the forthcoming reult of the paper, without recall it explicitly To tate the convergence reult ee Theorem 6, another Hypothei 5 will be added and dicued, that in particular contrain the value α, b and σ The ymmetrized Miltein cheme To complete our tak we follow the idea of Berkaoui, Boy and Diop in 3 who analyze the rate of convergence of the trong error for the ymmetrized uler cheme applied to quation Although, wherea they utilize an argument of change of time, we conider firt a weighted L p Ω-error for which we prove a convergence reult, and then we utilize thi reult to prove the convergence of the actual L p Ω-error We conider x >, T >, and N N We define the contant tep ize = T/N and t k = k Over thi dicretization of the interval, T we define the Symmetrized Miltein Scheme SMS X tk, k =,, N by x, for k =, X tk = X t k + bx tk + σx α t k W tk W tk + ασ for k =,, N Xα t k Wtk W tk, In the following, we ue the time continuou verion of the SMS, X t, t T atifying X t = X ηt + bx ηt t ηt + σx α ηtw t W ηt + ασ Xα ηt Wt W ηt t ηt, where ηt = up k {,,N} {t k : t k t} We alo introduce the increment proce Z t, t T defined by Z t = X ηt + bx ηt t ηt + σx α ηtw t W ηt + ασ Xα ηt Wt W ηt t ηt, 3 o that X t = Z t Thank to Tanaka Formula, the emi-martingale decompoition of X t i given by X t = x + where gnx = x gnz bx η d + L t X + gnz σx α η + ασ X α η W W η dw 4

3 Moment upper bound etimation for X and X We ummarize ome fact about the proce X t, t T, the proof of which can be found in Boy and Diop 5 Lemma For any q, there exit a poitive contant C depending on q, but alo on the parameter b, K, σ, α and T uch that, for any x >, X q t up t T C + x q 5 When < α <, for any q >, up t T X q t C + x q 6 When α = b, for any q uch that < q < σ, up t T X q t Cx q 7 Lemma 3 Let X t, t T be the olution of with < α < For all µ, there exit a poitive contant CT, µ, increaing in µ and T, depending alo on b, σ, α and x uch that exp µ T d X α When α = σ σ, the inequality 8 hold if b > and µ 8 b σ CT, µ 8 Notice that the condition b > σ / i alo impoed by the Feller tet in the cae α = for the trict poitivity of X, that allow to rewrite quation a X t = x + bx d + σ X dw Uing the emimartingale repreentation 4, we prove the following Lemma regarding the exitence of moment of any order for X t Lemma 4 For any q, there exit a poitive contant C depending on q, but alo on the parameter b, K, σ, α and T uch that for any x >, X q t C + x q up t,t The proof of thi lemma i baed on the Lipchitz property of b and claical combination of Itô formula and Young Inequality For the ake of completene, we give a hort proof in the Appendix Strong rate of convergence The main reult of thi work i the trong convergence at rate one of the SMS X to the exact proce X The convergence hold in L p for p To tate it, we add to Hypothei the following For any x in R +, we denote x the rounded up integer Hypothei 5 i Let p To control the L p Ω-norm of the error, if α > we aume b > α α σ Wherea for α = we aume b > 3p + σ / ii The drift coefficient b i of cla C R, and b ha polynomial growth 3

4 We now tate our main theorem To lighten the notation, we conider for α, b σ α := b α ασ, Kα := K + ασ α α α α and we extend thi definition to α = taking limit So, b σ = lim α b σ α = b σ /4, and K = lim Kα = K Notice that lim α Kα = K + σ /, and ince Kα i continuou on α,, we have that Kα i bounded Thi i epecially important in the definition of max α bellow, becaue tell u that α max α i trictly poitive and bounded on, Theorem 6 Aume Hypothee and 5 Define a maximum tep ize max α a x max α = αb σ α 4αKα, for α, 4K x, for α = Let X t, t T be the proce defined on and X t, t T the ymmetrized Miltein cheme given in Then for any p that allow Hypothee 5, there exit a contant C depending on p, T, b, α, σ, K, and x, but not on, uch that for all max α, 9 up X t X t p p C t T About Hypothei 5 Notice that for α >, Aumption i doe not depend on p and become eaier to fulfill a α increae On the other hand, for α =, Aumption i depend on p in a unpleaant manner However, a we will ee later in Section 4 ee Table, thi kind of dependence in p i expected, and imilar condition are aked in the literature for other approximation cheme in order to obtain imilar rate of convergence reult Alo, notice that i i a ufficient condition: in the numerical experiment we till oberve a rate of convergence of order one for parameter that do not atify it, but we alo oberve that for parameter uch that b σ, although the convergence occur, it doe in a ublinear fahion On the other hand, Aumption ii i the claical requirement for the trong convergence of the Miltein cheme A we will ee later in the proof of the main theorem, with the help of the Itô formula, thi hypothei let u conclude that X X p bx η bx C up X u X u p + p u intead of X X p bx η bx C up X u X u p + p, u which i the claical bound obtained for the uler-maruyama cheme under a Lipchitz condition for a drift b The ret of the paper i organized a follow In Section we tate ome preliminary reult on the cheme which will be building block in the proof of Theorem 6 Section 3 i devoted to the proof of the convergence rate The main idea i firt to introduce a weight proce in the L 4p Ω-error Wet get the rate of convergence for thi weighted error proce, and we ue thi intermediate bound to control the L p Ω-error, from where we finally control the L p Ω-error Alo, a a byproduct we obtain the order one convergence of the Projected Miltein Scheme ee 33 In ection 4, we diplay ome numerical experiment to how the effectivene of the theoretical rate of convergence of the cheme, but alo to tet Hypothee 5-i on a et of parameter In thi ection we alo how how the incluion of the SMS in a Multilevel Monte Carlo framework could help to optimize the computational time of weak approximation of aet valuation In Section 5 we preent the proof of the preliminary reult on the cheme Finally we have included in a mall appendix a couple of proof to make thi paper elf contained 4

5 Some preliminary reult for X Thi hort ection i devoted to tate ome reult about the behavior of X, their proof are potponed to Section 5 All thee reult hold under Hypothei 5-i which i in fact tronger than what we need here So, we preent the next lemma with their minimal hypothee till auming Hypothei Lemma Local error For any x >, for any p, there exit a poitive contant C, depending on p, T, the parameter of the model b, K, σ, α, but not on uch that up X t X ηt p C p t T By contruction the cheme X i nonnegative, but a key point of the convergence proof reide in the analyi of the behavior of X or Z viiting the point The next Lemma how that although Z t i not alway poitive, the probability of Z t being negative i actually very mall under uitable hypothee Lemma For α,, if b > α α σ, and Kα, then there exit a poitive contant γ, depending on the parameter of the model, but not on, uch that In particular, up P k=,,n inf Z t k < t k+ exp γ up P Z t C exp γ t T To prove Lemma, it i neceary to etablih before the following one, which although technical, give ome intuition about the difference between the SMS and the Symmetrized uler cheme preented in 3 Lemma 3 For α,, if b > α α σ, we et xα = bσα Kα > Then for all t, T, and for all ρ,, P Z t ρb σ α, X ηt < ρ xα = Roughly peaking, from thi lemma we ee that when Z ηt >, Z t become negative only when Z t Z ηt > ρ xα But oberve that only the left-hand ide of thi inequality depend on, and it expectation decreae to zero proportionally to, according to Lemma Now impoing mall enough, we prove an explicit bound for the local time moment of X Lemma 4 For α,, if b > α α σ and Kα x αb σα, then there exit poitive contant C and γ > depending on α, b, K, and σ but not in uch that L T X C exp γ We end thi ection with another key preliminary reult, which i the convergence rate of order for the corrected local error Although the claical local error i of order /, a tated in Lemma, the local error een by the diffuion coefficient function, corrected with the Miltein term tay of order Lemma 5 Corrected local error proce Let u fix p, and α, For α >, aume b > α α σ, wherea for α =, aume b > 3p + σ / Then, there exit C >, depending on the parameter of the model but not in, uch that for all max α, the Corrected Local rror atifie up t T σx α t σx α ηt ασ X α ηt W t W ηt p C p 5

6 3 Proof of the main Theorem 6 The proof of Theorem 6 i built in everal tep Firt, we work with the L p Ω-norm of the error, for p, then at the lat tep of the proof we go back to the L p Ω-norm for p In what follow we denote t := X t X t and o that Σ t := gnz t σx α ηt + ασ X α ηt W t W ηt σxt α d t = gnz t bx ηt bx t dt + dl t X + Σ t dw t Alo, to make the notation lighter, we will denote the Corrected Local rror by 3 The Weighted rror D t X := σx α t σx α ηt ασ X α ηt W t W ηt Before to prove the main theorem, we etablih in the Propoition 3 the convergence of a weighted error For p, let u conider β t, t T, defined by and the Weight Proce Γ t, t T defined by β t = p b + p4p + 8α p4p σ, 3 Γ t = exp X α t β d 3 The Weight Proce i adapted, almot urely poitive, and bounded by It path are non increaing and hence ha bounded variation, and alo atifie dγ t = β t Γ t dt The proce Γ t t can be een a an integrating factor in the ene of linear firt order OD ee for example 9 When we apply the Itô Lemma to Γ t 4p t, intead of 4p t alone, we can remove a very annoying term that appear in the righthand ide bound ee the proof of Lemma 33 The exponential weight in a L p Ω-norm i a ueful tool to obtain a priori bound a an example, for the exitence and uniquene of the olution of a Backward SD, it i introduced a norm with exponential weight, uch that the operator aociated to the BSD i contractive under thi new norm ee proof of Theorem in 7 In the ame way, here we introduce thi exponential weight to get the following a priori error-bound, which will allow u to prove Theorem 6 Propoition 3 Weighted rror Under the hypothee of Theorem 6, for p and α,, there exit a contant C not depending on uch that for all max α Remark 3 up t T Γ t 4p t C 4p 33 i Since we are going to work firt with the L p Ω-norm of the error, when α =, Hypothei 5-i become in particular b > 3p + σ = 34p + σ, b σ > p

7 ii From Lemma, the proce β ha polynomial moment of any order for α >, and when α =, there exit C uch that β q t < C, for all < q < b/σ From the previou point in thi Remark, it follow that the proce β ha moment at leat up to order p + iii From the definition β in 3, and due to Lemma 3, there exit a contant C uch that Γ q T < C for all q > when α >, wherea for α =, the q-th negative moment of the weight proce i finite, a oon a p4p σ q σ 8 b σ Notice that, thank to point i in thi Remark, a ufficient condition uch that thi lat inequality hold i 6p4p q p + We cut the proof of Propoition 3 in two technical lemma Lemma 33 Under the hypothee of Theorem 6, for p and α,, there exit a contant C not depending on uch that for all max α Γ t 4p t 4p Γ 4p bxη bx d + 4p b up u Γ u 4p u d + C 4p 35 Lemma 34 Under the hypothee of Theorem 6, for p and α,, there exit a contant C not depending on uch that for all max α, and for any, T Γ 4p bx η bx C up u Γ u 4p u + C 4p 36 A we will ee oon, we prove 36 with the help of the Itô formula applied to b, and here i where we need b of cla C required in Hypothee 5-ii Proof of Propoition 3 Thank to Lemma 33 and 34, we have Γ t 4p t C and ince the right-hand ide i increaing, it follow up t Γ 4p C up u up u from where we conclude the reult thank to Gronwall Inequality Now we preent the proof of the technical lemma Proof of Lemma 33 By the integration by part formula, Γ t 4p t = 4p + p4p + 4p Γ u 4p u d + C 4p Γ u 4p u d + C 4p, Γ 4p { gnz bx η bx } d Γ 4p Σ d Γ 4p dl X β Γ 4p d 7

8 Thank to Lemma 4 and the control in the moment of the exact proce in Lemma we have Γ 4p dl X X 4p dl X up X 8p T L T X C 4p On the other hand, with gnx = {x<}, calling W = W W η, we get for all t Σ σx α σx α + σx α σx α η ασ X α η W + R Σ, {Z<} where we put aide all the term multiplied by {Z<} in R Σ :=4 σx α η + ασ X α η W { σx α η + ασ X α η W + } σx α σx α + σx α σx α η ασ X α η W So, from the previou computation, the Lipchitz property of b, and Young Inequality, we conclude Γ t 4p t 4p where R = 4p4p 4p +4p b Γ 4p bxη bx d +4p4p +4p 4p +4p Γ η 4p η d + 4p b Γ 4p σx α σx α D X 4p d + β Γ 4p d Γ 4p d d Γ 4p d R {Z<} d + C 4p R Σ + 8p 4p bx η, and from Lemma we have R {Z<} C 4p Since max α and Remark 3-i, we can apply Lemma 5 o D X 4p C 4p Introducing thee etimation in the previou computation, we have Γ t 4p t 4p Γ 4p bxη bx d + 4p b Γ η 4p η d + 4p b + 4p 4p + 4p4p β Γ 4p d Γ 4p σx α σx α + C 4p Γ 4p d d 8

9 Now we ue the particular form of the weight proce Since for all α, for all x, y, we have x α y α x α + y α α x y, 37 Γ 4p σx α σx α and then, from the definition of β in 3, we conclude from where Γ t 4p t 4p Γ t 4p t 4p Γ 4p d Γ 4p bxη bx d + 4p b Γ 4p bxη bx d + 4p b 4α σ X α d, Γ η 4p η d + C 4p up u Γ u 4p u d + C 4p Proof of Lemma 34 By integration by part Γ 4p bx η bx = + η η Γ u u 4p bx η bx u β u Γ u du Γ ud 4p u bx η bx u Applying Hölder Inequality to the firt term in the right-hand ide we have Γ uu 4p bx η bx u β u du Γ u u 4p 4p bx η bx u 4p βu 4p 4p du η Recalling the Remark 3-ii, we have that β 8p u i finite, o applying Lemma, η 38 Then, bx η bx u 4p βu 4p bx η bx u 8p βu 8p C p η Γ u 4p u bx η bx u β u du C up Γ uu 4p u 4p 3/ Applying the Itô Formula to the econd term in the right-hand ide of 38, and taking expectation we get Γ ud u 4p bx η bx u η = σ Γ uu 4p η +4p η + σ 4p 4p 4p σ + 4p η η =: I + I + I 3 + I 4 + I 5 b X u bx u + σ b X u Xu α du Γ uu 4p bx η bx u { gnz u bx η bx u } du η Γ uu 4p 3 bx η bx u Σ udu Γ uu 4p b X u Xu α Σ u du Γ uu 4p bx η bx u dl ux 39 9

10 By the finitene of the moment of X, the linear growth of b, and the polynomial growth of b, applying Holder inequality, we have I C Γ u u 4p 4p b X u bx u + σ η b X u Xu α 4p 4p C up Γ uu 4p 4p u For the bound of I, ince b i Lipchitz, and gnx = {x<}, we have I C + C η Γ uu 4p X η X u X η X u du η C up Γ uu 4p u Γ u 4p u X η X u du + p + C η up Γ uu 4p u R u du {Zu<} du 4p 3/ + C 4p Where again, all the term multiplied by {Zu<} are putted in the ret R u, and the expectation of the product i bounded with Lemma In a imilar way for the bound of I 3, decompoing Σ u with gnx = {x<}, I 3 C + C η Γ uu 4p 3 X η X u D X du η Γ uu 4p 3 For the firt term in the right-hand ide we have Γ u 4p 3 u X η X u σx u σx α u du + η R u 3 du {Zu<} X η X u D X Γ uu 4p 3 4p X η X u 4p 4p D X 4p p 3 4p 5/, up Γ uu 4p u due to the bound for the increment of the exact proce and Lemma 5 For the econd term, applying 37, and noting that from Remark 3-ii, the exact proce ha negative moment up to of order p + Γ uu 4p 3 X η X u σx u σxu α CΓ uu 4p X η X u X α u C Γ uu 4p 4p X η X u 8p /8p up Γ uu 4p u 4p / X 6 αp u /8p We control the third term in the right-hand ide in the bound for I 3 uing again Lemma, o I 3 C up Γ uu 4p u 3 4p 7/ + up Γ uu 4p u 4p 3/ + C 4p

11 Now we bound I 4 I 4 C + C η η Γ uu 4p D u Xb X u Xu α du Γ u 4p u σx u σx α u b X u X α u du + η R u 4 du {Zu<} We control the firt term in the right-hand ide uing Hölder inequality, Lemma 5 and the control in the moment of the exact proce for all u Γ u 4p u D u Xb X u Xu α Γ u 4p C u upγ uu 4p u p D X 4p 4p b X u 4p Xu 4pα 4p p For the econd term in the right-hand ide of the bound for I 4, we ue one more time 37, and the exitence of negative moment of the exact proce X, and then Γ uu 4p σxu σxu α b X u Xu α CΓ uu 4p X α u b X u Xu α C up Γ uu 4p u To control the third term in the right-hand ide of the bound for I 4 we ue Lemma jut a before So Finally, I 4 C up Γ uu 4p u 4p I 5 = C up u C η p + C up Γ uu 4p u + X 4p u up u the lat inequality come from Lemma and 4 Putting all the lat calculation in 38 we obtain Γ 4p bx η bx C up u + C 4p + C 4p Γ uxu 4p bx η bx u dl ux + X 4p u + C 4p L T X L T X C 4p, Γ uu 4p 4p 3/ + C up Γ uu 4p u p + C Applying Young Inequality in all term in the right, we get the deired inequality 36 up Γ uu 4p u up Γ uu 4p u Remark 35 Proceeding a in the Proof of Lemma 34, if p = or p 3/, i not difficult to prove p bx η bx C up u 4p 4p 3 4p 7/ p u + C p 3

12 Notice that 3 i imilar to 36 with the proce Γ = The retriction on p come from the following obervation Applying the Itô Lemma to the function x, y x p y for p, with the couple of procee, bx η bx, in the analogou of the identity 39 with Γ =, it will appear a term of the form η u p 3 bx η bx u Σ udu The retriction p 3/ avoid the ituation where p 3 i negative 3 Proof of Theorem 6 We tart by controlling the L p -error Firt, aume p = or p 3/ By the Itô formula we have p t = p p A we have een before, { gnz bx η bx } d + p p p t p + p where R = 4pp p p dl X + pp dl X C p, and gnx = {x<}, o p bxη bx η d + 8pp + 8pp p bxη bx d σx α η + ασ X α η p σx α σx α d p D X d + R {Z<} d + Cp, p Σ d 3 W W η + 8p p bx η If we ue the Lipchitz property of b, and Young inequality in the firt term in the right of 3, Lemma 5 in the fourth one, and Lemma in the fifth one, we have p t C up u + 8pp And, according to Remark 35, we have p p u d + p p bx η bx d p σx α σx α d + C p bx η bx C up u p u + C p 3 On the other hand, uing again 37, we have p and applying Cauchy-Schwartz inequality, σx α σx α C p Γ p X α X α Γ Γ 4p = CΓ p X 4 α X α Γ Γ, The firt term in the right-hand ide i the weight error controlled by Propoition 3 To control the econd one, let u recall Remark 3 For α >, the exact proce and the weight proce Γ have negative moment of any order,

13 therefore the econd term in the lat inequality i bounded by a contant On the other hand, for α = we need a finer analyi From the econd point in Remark 3, the p + -th negative moment of the exact proce i finite, and ince 6p4p 6p + < p +, 6p according with the third point of Remark 3, the 6p + /6p-th negative moment of the weight proce Γ i alo finite Therefore, when α =, and then in any cae X 4 α Γ = X Γ Γ p Introducing all the lat computation in 3 we get p t C X α up u X p+ 6p+ Γ C p Γ u p d + C p 6p+ 6p 6p 6p+ C, Since the right-hand ide i increaing, thank to Gronwall Inequality we have, for p = or p 3/, We extend to p, 3/, thank to Jenen inequality up p p t C t T p t 6 p t 6 C, and we conclude that p t p C for all p atifying Remark 3-i Now we control the L p -error For α = and p, denoting p = p, Hypothei 5-i give up p t p = up p p t t T t T b > 3p + σ = 34p + σ C Since 34 in Remark 3 i atified, we can control the L p Ω-norm of the error and then up p t p = up p p t t T t T C If p,, Hypothei 5-i i b > 5 σ, which i enough to bound the L Ω-norm of the error, and then from Jenen inequality up t T p p up t T t C The cae α > / i eaier Since the Hypothei in the parameter for thi cae doe not depend on p, we can conclude for any p from Jenen inequality and the control for the L p Ω-norm of the error 3

14 Remark 36 Let u mention an example of extenion of our convergence reult, baed on imple tranformation method: conider the 3/-model, namely the olution of r t = r + c r c r d + Applying the Itô Formula to v t = fr t, with fx = x, we have v t = v + c + c 3 c c v d + c 3 r 3/ dw c 3 v / db, where B = W i a Brownian motion We can approximate v with the SMS v, and then define r t := / v t Then we can deduce the trong convergence with rate one of r t to r t from our previou reult Tranformation method can be ued in a more exhautive manner, in the context of CV-like SD and we refer to 6 for approximation reult and example, uing thi approach 33 Strong Convergence of the Projected Miltein Scheme The Projected Miltein Scheme PMS i defined by X tk = Xtk + b X + σ X α W ασ Xα t tk tk t k W tk + k where for all x R, x + = max, x The continuou time verion of the PMS i given by X t = Xηt + b X ηt + σ X ηt α W t W ηt + ασ Xα ηt Wtk W tk +, X = x Wt W ηt + 33 Notice that for all t, T, X t X t, then the poitive moment of the PMS are bounded ee Lemma 4 To obtain a trong convergence rate for the PMS, we firt how that the PMS and the SMS coincide with a large probability Lemma 37 Let u conider the topping time τ = inf{ : X X } Aume that b > α α σ and /Kα Then for any p, P τ T C p Proof Notice that τ i almot urely trictly poitive becaue both cheme tart from the ame determinitic initial condition x On the other hand P τ T = N and according to Lemma P τ t i, t i+ = P inf Z, X tk = X tk t k < t k+ i= P τ t k, t k+, P inf Z t k < t k+ exp γ So, P τ T T exp γ Since for any p, there exit a contant C p uch that exp γ// C p p, we have P τ T C p 4

15 Corollary 38 Aume Hypothee and 5 Conider a maximum tep ize max α defined in Let X t, t T be the proce defined on and X t, t T the Projecter Miltein cheme given in 33 Then for any p that allow Hypothee 5, there exit a contant C depending on p, T, b, α, σ, K, and x, but not on, uch that for all max α, Proof Notice that for all t, T, with τ = inf{ : X X }, up X t X t p p C 34 t T X t X t p = X t X t p {τ T } + X t X t p {τ>t } X t X t p P τ T + X t X t p C p where the lat inequality come from Lemma 37 and Theorem 6 4 Numerical xperiment and Concluion We tart thi ection with the analyi of two numerical experiment The firt one aim to tudy empirically the trong rate of convergence of the SMS in comparion with other cheme propoed the literature The econd one aim to tudy the impact of including the SMS in a Multilevel Monte Carlo application 4 mpirical tudy of the trong rate of converge In thi experiment we compute the error of the cheme a a function of the tep ize for different value of the parameter α and σ For α > we compare the SMS with the Symmetrized uler Scheme SS introduced in 3, and with the Balanced Miltein Scheme BMS preented in 3 Wherea for α =, in addition to the aforementioned cheme, we will alo compare SMS with the Modified uler Scheme MS propoed in 6, and with the Alfoni Implicit Scheme AIS propoed in Let u firt, hortly review thoe different cheme Alfoni Implicit Scheme AIS Propoed in, the AIS can be applied to equation when the drift i a linear function A priori, the AIS can be applied for α,, but i relevant to oberve that only when α =, the AIS i in fact an explicit cheme alo know a drift-implicit quare-root uler approximation wherea in any other cae i not Thi implie that in order to compute the AIS for α >, at each time tep it i neceary to olve numerically a non linear equation Thi extra tep in the implementation of the cheme bring quetion about the impact of the error of thi ubroutine in the error of the cheme, and about the computing performance of the cheme Since thi quetion are beyond the cope of the preent work, we include the AIS in the comparion only in the Cox Ingerol Ro CIR cae linear drift and α = In thi context, the AIS can be ue only if σ > 4b, for other value of the parameter the AIS i not defined In term of convergence, when α =, according to Theorem in Alfoni, the AIS converge in the L p Ω-norm, for p, at rate when 3p/4σ < b When α >, the AIS ee Section 3 of converge a oon a b >, at rate to the exact olution Balanced Miltein Scheme BMS The BMS wa introduced by Kahl and Schurz in 3, and although it convergence it i not proven for quation ee Remark 5 in 3, numerical experiment how a competitive behavior ee Alo, the BMS can be eaily implemented for α,, o we decide to include it in our numerical comparion 5

16 Modified uler Scheme MS Introduced in 6, the MS can be applied to the quation for α, when the drift ha the form bx = µ x µ xx for µ and µ uitable function The rate of convergence in the L Ω-norm of the MS depend on the parameter For α = the rate i if σ i big enough compared with b, and it i ρ < in other cae When α >, the MS converge at rate a oon a b > ee Propoition 4 in 6 When α >, the implementation of the MS require ome extra tuning which i not explicitly given in 6 ee Remark 5, o we implement the MS only for α = Symmetryzed uler Scheme SS The SS, introduced in 3, i an explicit cheme which can be apply to the equation for α, and any Lipchitz drift function b It ha the weaket hypothei over b of all the cheme dicued in thi paper If α =, according to Theorem in 3, the rate of convergence of the SS i under uitable condition for b, σ and K When α > the SS converge at rate a oon a b > ee Theorem in 3 We ummarize the theoretical analyi of the cheme above in Table for α =, and in Table for α > Scheme Norm Drift Convergence Condition SMS L p, p AIS BMS b Lipchitz, b C b with polynomial growth b > 3 p + σ Theoretical rate L p, p, 4b 3σ bx = a bx b > 3 4 pσ bx = µ x µ xx b > 5σ MS 6 L µ i Cb C b, µ b > 3σ undetermined µ, µ b > σ 6, σ 8 b > σ K p + σ, SS 3 L p, p b Lipchitz Kq = K6q 4σ 8p b+σ Table : Summary of the condition over the parameter for the convergence of the different cheme for α = Simulation etup In our imulation we conider a time horizon T =, and x = In order to include a many cheme a poible we conider for all imulation a linear drift bx = x To meaure the error of each cheme, we etimate it L Ω-norm for which a theoretical rate i propoed for all the elected cheme Let T SMS BMS SS MS AIS, T, T, T, and T be the L Ω-norm of the error for the SMS, BMS, SS, MS and AIS repectively To etimate thee quantitie, we conider a a reference olution the AIS approximation for = max α/ when α =, and the SMS for = maxα/ when α > Then for each { } max α n, n =, 9, we etimate T by computing 5 4 trajectorie of the correponding cheme, and comparing them with the reference olution The reult of thee imulation are reported in Figure α = and α > The graph 6

17 Scheme Norm Drift Convergence Condition SMS L p, p AIS BMS b Lipchitz, b C b with polynomial growth Theoretical rate b > α α σ L p, p, 4b 3σ bx = a bx b > MS L µ i Cb C b, µ bx = µ x µ xx µ, µ undetermined b > SS L p, p b Lipchitz b > Table : Summary of the condition over the parameter for the convergence of the different cheme when α > plot the Log T in term of the Log, and we have added the plot of the identity map to erve a reference for rate of order The cheme with a lope maller than the lope of the reference line have an order of convergence maller than one To obtain a more quantitative comparion of the cheme, we alo perform a regreion analyi on the model log T = ρ log + C Notice that ˆρ, the etimated value for ρ, correpond to the empirical rate of convergence of the different cheme We preent the reult of thi regreion analyi in Table 3 and 4 mpirical reult for α = Figure and Table 3 preent the reult for the CIR cae From Table, we oberve that we can ditinguih five cae for the parameter The firt cae σ = i uch that b > 6σ : the SMS, the AIS, and the MS have a theoretical rate of convergence equal to, wherea the SS ha a theoretical rate of convergence equal to In Figure a, we oberve that the graph of the SMS, the AIS, the BMS, and the MS eem parallel to the reference line, which i expected, while the SS ha a maller lope Thi i alo confirm in the firt line of the Table 3, where we oberve that the empirical rate of convergence are cloe to the theoretical one Notice that the BMS ha a competitive empirical rate of convergence, although the theoretical one i not known The econd cae σ = 4 i uch that b 5 σ, 6σ, now only the AIS and the MS have a theoretical rate of convergence equal to However, how we can ee in Figure b the SMS till how a linear behavior in thi cae Recall that the condition over the parameter i a ufficient condition and we believe that could be improved Notice that alo the BMS how a linear behavior In the econd line of Table 3, we can oberve that empirical rate of convergence are cloe to one for all the cheme but the SS In Figure c, we illutrate the third cae σ = 65 and then b 3σ /, 5σ / In thi cae, only the AIS ha a theoretical rate of convergence equal to one For the MS i, but in the graphic we till oberve a linear behavior for the MS, and alo for the SMS and the BMS Thi i confirm in the third line of Table 3 The fourth cae i σ = 9 and b σ, 3σ /, which we diplay in Figure d For thi value of the parameter the theoretical rate of convergence i known only for the AIS and the MS Neverthele, we oberve in the graph and in the fourth line of Table 3 that all the cheme eem to reach their optimal convergence rate Finally, the fifth cae i σ = 36 and then b < σ In thi cae all the cheme have a ublinear behavior a we can ee in Figure d and the fifth line of Table 3 Thi cae illutrate the neceity of ome condition over the parameter of the model to obtain the optimal rate of convergence for the SMS mpirical reult for α > α = 7 In Figure and Table 4 we preent the reult of the imulation for α = 6, and 7

18 σ Oberved L Ω convergence rate ˆρ and it R value SMS AIS BMS MS SS ˆρ R ˆρ R ˆρ R ˆρ R ˆρ R % 6 999% % % % % % % % % % % 999% % % % 5 999% % % % % % % % % Table 3: mpirical rate of convergence ˆρ for the L Ω-error of the cheme when α = for different value of σ In thee cae, it can be oberved in numerical experiment that the MS need maller to achieve it theoretical order one convergence rate, unle one tune the projection operator in the manner of Remark 5 in 6 Since thi tuning i not explicitly given we do not include the MS in thee imulation Parameter Oberved L Ω convergence rate ˆρ and it R value α σ SMS BMS SS ˆρ R ˆρ R ˆρ R % % % % % % % % % % 9 997% % % % % % % 64 99% Table 4: mpirical rate of convergence ˆρ for the L Ω-error, when α > for different value of α and σ We have oberved in the numerical experiment three cae for the parameter The firt one i when b > α α σ σ = 49 and σ = 64 In thi cae, Theorem 6 hold and we oberve the order one convergence ee Figure a, b, and firt and fourth row in Table 4 The econd cae i when the parameter do not atify b > α α σ σ = 539 and σ = 8, and then we can not apply Theorem 6, but in the numerical imulation we till oberve the order one convergence ee Figure c, d, and econd and fifth row in Table 4 Finally the third cae, i when σ b, and then we do not oberve a linear convergence anymore ee Figure e, f, and third and ix row in Table 4 Notice that in the three cae the SMS perform better than the BMS, pecially when σ grow ee Table 4 8

19 The econd and third cae how u that ome retriction ha to be impoe on the parameter to oberve the convergence of order one But our retriction, although ufficient, it eem to be too trong, pecially for α cloe to one 4 Application of the SMS in Multilevel Monte Carlo We continue thi ection by teting the SMS in the context of a multilevel Monte Carlo application widely ued nowaday in computational finance ee eg and reference therein Multilevel Monte Carlo i an efficient technique introduced by Gile 9 to decreae the computational complexity of an etimator combining Monte Carlo imulation and time dicretiation cheme for a given threhold in the accuracy For detail we refer to 9,, For thi experiment, we conider the claical but non trivial tet-cae of the Zero Coupon Bound ZCB pricing of maturity T, T B, T = exp r d, under the hypothei that the hort interet rate dynamic r t, t i modeled with a CIR proce α = and bx = a bx : dr t = a br t dt + σ r t dw t In thi context, the price of the ZCB admit a wellknown cloed-form olution given by ee eg 7, 4 B, T = AT e BT r, AT = λe b+λt/ λ + be λt + λ a σ, BT = e λt λ + be λt + λ where r i the initial value of the interet rate, and λ = b + σ Let B l a dicrete-time weak approximation of B, T with the time tep l We conider the L-level Monte Carlo etimator : Ŷ T = N N i= B i + L N l N l l= i= For a targeted mean-quare error ɛ on the computation of the quantity B, T ŶT B, T = Oɛ, Bi l B i l one can chooe the following a priori parametrization of the MLMC method in order to minimize the computational log ɛ time complexity ee 9,, : we ue the etimation L = ; from one level to the next one the time tep i log divided by, l = l+ ; the number of trajectorie to imulate i etimated with Gile formula 9 N l = ɛ Vl l L l= V l / l, with V l = Var B l B l A an etimator for the bia variance, the trong rate of convergence of the dicretiation cheme enter a a key ingredient in the N l a priori etimation A cheme with a reduced trong bia will then allow a maller N l We apply the MLMC computation for the SMS, the PMS, the AIS and the BMS We ummarize the reult of the performance comparion between the four cheme in Table 5 The computation have been run uing a initial interet rate r =, the maturity of the bond T =, the drift parameter a = b =, the volatility σ = For the MLMC imulation, we fix a minimum number of trajectorie equal to 5, and a minimal number of level equal to 6 In Table 5, we give the meaure of the CPU time for a et of three decreaing targeted error, a long a the effective meaured error and the total number of imulated trajectorie A expected, the required threhold error ha been reached by the MLMC trategy A alo expected ee Gile, Miltein cheme perform better than their uler verion Finally, a ɛ decreae, the SMS clearly perform better than hi PMS verion 9

20 - Parameter: α = 5, σ =, max = 5 - Parameter: α = 5, σ = 4, max = logrror log Ref SMS SS BMS MS AIS logrror log Ref SMS SS BMS MS AIS - a Parameter in cae : b > 6σ Parameter: α = 5, σ = 65, max = 5 - b Parameter in cae : b 5σ /, 6σ Parameter: α = 5, σ = 9, max = logrror log Ref SMS SS BMS MS AIS logrror log Ref SMS SS BMS MS AIS c Parameter in cae 3: b 3σ /, 5σ / d Parameter in cae 4: b σ, 3σ / Parameter: α = 5, σ = 36, max = 5 - logrror log Ref SMS SS BMS MS AIS e Parameter in cae 5: b < σ Figure : Step ize veru the etimated L Ω-trong error for the CIR Proce Log-Log cale The identity map erve a a reference line of rate one

21 Parameter: α = 6, σ = 49, max = 43 Parameter: α = 7, σ = 64, max = - - logrror log Ref SMS SS BMS logrror log Ref SMS SS BMS a Parameter in cae : α = 6 Parameter: α = 6, σ = 539, max = b Parameter in cae : α = 7 Parameter: α = 7, σ = 8, max = - - logrror log Ref SMS SS BMS logrror log Ref SMS SS BMS c Parameter in cae : α = 6 Parameter: α = 6, σ = 44, max = d Parameter in cae : α = 7 Parameter: α = 7, σ = 5, max = - - logrror - logrror log SMS Ref SS BMS log Ref SMS SS BMS e Parameter in cae 3: α = 6 f Parameter in cae 3: α = 7 Figure : Step ize veru the etimated L Ω-error for α > and different value for σ Log-Log cale The identity map erve a a reference line of rate one

22 ɛ = e-3 L = 9, L = / SMS PMS AIS BMS CPU time N + + N L oberved error e e e e-4 ɛ = e-4 L = 3, L = / 4 SMS PMS AIS BMS CPU time N + + N L oberved error e e e e-6 ɛ = e-5 L = 6, L = / 7 SMS PMS AIS BMS CPU time N + + N L oberved error e e e e-7 Table 5: CPU time to achieve the target error for the different cheme The oberved error i ŶT B, T 43 Concluion In thi paper we have recovered the claical rate of convergence of the Miltein cheme in a context of non mooth diffuion coefficient, although we have to impoe ome retriction over the parameter of the SD to enure the theoretical order one of convergence Typically, if the quotient b/σ i big enough we will oberve the optimal convergence rate In the numerical imulation we have oberved that, depite the fact it i neceary to impoe ome retriction over the parameter of SD to obtain the order one convergence, Hypothei 5 eem to be not optimal, pecially for α = Alo, through numerical imulation, we have oberved that the ue of SMS could improve the computation time in a Multilevel Monte Carlo framework, at leat a well a the CIR pecialized one-order cheme Although our reult eem more retrictive in term of hypothee on the et of parameter, in particular if we compare SMS with Lamperti tranformation-baed cheme ee the recent work in and 6, the SMS can be applied to a more general cla of drift function and in variou context It i thu a ueful complement of the exiting literature

23 5 Proof for preliminary lemma 5 On the Local rror of the SMS Proof of Lemma From the definition of X, and the algebraic inequality for poitive real number a + + a n p n p a p + + ap n we have Xt X ηt p 3 p bx ηt p t ηt p + σ p X αp ηtw t W ηt p + αp σ 4p p X α p ηt Wt W ηt t ηt p Thank to the linear growth of b, Lemma 4 and the propertie of the Brownian Motion it i quite eay to conclude the exitence of a contant C uch that X t X ηt p C p, from where the reult follow 5 On the Probability of SMS being cloe to zero From b σ α and Kα defined in 9, let u recall the notation xα := b σα Kα introduced in Lemma 3 A b σ α > under Hypothei 5-i, xα i bounded away from In particular, lim α xα = b σ /4, wherea lim xα = K α b K + σ / Proof of Lemma 3 Denoting W = W W η, and = η, we have for all, T, Z = ασ Xα η W + σx α η W + X η + bx η ασ Xα η From the Lipchitz property of b and the following bound for any x > we have So, x α 4 α + α α α α x, 5 Z ασ Xα η W + σx α η W + X η + b σ α KαX η 5 P Z ρb σ α, X η < ρ xα ασ P Xα η W + σx α η W + X η + ρb σ α KαX η, Xηt < ρ xα From the independence of W with repect to F η, if we denote by N a tandard Gauian variable, we have P Z ρb σ α, X η < ρ xα ασ P xα N + σ x α N + x + ρb σ α Kαx {Xη <ρ xα} x=xη 3

24 Notice that in the right-hand ide we have a quadratic polynomial of a tandard Gauian random variable Let u compute it dicriminant: x, α = σ x α ασ x α x + ρb σ α Kαx = α σ x α ασ x α ρb σ α Kαx Since b σ α >, we have x, α < for all α,, and x ρ xα So, for all we have taking = we conclude on the Lemma P Z ρb σ α, X η < ρ xα =, Proof of Lemma We have P inf Z = P inf Z, X tk xα t k <t k+ t k <t k+ + P inf Z, X tk < xα t k <t k+ 53 We tart with the econd term in the right hand of the lat inequality By continuity of the path of Z and Lemma 3, we have P inf Z, X tk < xα = P Z, X tk < xα = t k <t k+ On the other hand, from 5 we have Q t k,t k+ Z σx α η W + Kα X η + b σ α 54 Then P inf t k <t k+ Z, X tk xα α Xη bσ α KαX η P inf + t k < t k+ σ σx α η = ψx tk {Xtk xα}, + W, X tk xα where the lat equality hold thank to the Markov Property of the Brownian motion, for ψx = P inf <u x α σ + b σα Kαx σx α u + B u, where B t denote a Brownian Motion independent of W t If B µ t, t T i a Brownian motion with drift µ, tarting at y, then for all y y, we have ee 4: P inf < t Bµ y = erfc y y + µ t t + exp µy y y y erfc µ t, t 55 where for z R, erfc z = /π z exp u / du In our cae y = x α /σ, µ = b σ α Kαx/σx α, and y = Then ψx = Kαx + erfc bσ α σx α + exp b σα Kαxx σ x α erfc x bσ α Kαx σx α 4

25 Since /Kα, for any x xα, the argument in the erfc function in the lat equality are both poitive, and then recalling that for all z > erfcz exp z, we obtain So, for all x xα Then P ψx exp Kαx + b σα σ x α + exp b σα Kαxx σ x α exp x b σα Kαx σ x α exp Kαx + b σα σ x α inf t k <t k+ Z, X tk ψx exp Kα x α σ xα exp and finally, chooing γ = xα α /8σ, we get exp P inf Z t k <t k+ Kα X α t k σ Kα xα α σ exp γ, {Xtk xα} 53 On the Local Time of the SMS at Zero The Stopping Time Θ α, α < In what follow, we conider Θ α = inf { > : X < αb σ α } 56 Lemma 5 Aume b > α α σ, and /Kα x / αb σ α Then there exit a poitive contant γ depending on α, b, K and σ but not on uch that PΘ α T T γ exp 57 Proof Firt, notice that the condition < x / αb σ α enure that the topping time Θ α i almot urely trictly poitive To enlighten the notation along thi proof, let u call l σ α := αb σ α, and ζ k = inf tk < t k+ Z We plit the proof in three tep Step Let u prove that for a uitable function ψ : R, and the et A k = {X tk > xα α} F tk : P Θ α T N k= ψx tk Ak 58 5

26 Indeed, But, for each k =,, N P Θ α T N k= P ζ k l σ α, X tk > l σ α P ζ k l σ α, X tk > l σ α = P ζ k l σ α, X tk > l σ α, X tk < xα α Since l σ α = αb σ α xα α, we have P ζ k l σ α, X tk > l σ α, X tk < xα α + P ζ k l σ α, X tk > l σ α, X tk xα α P ζ k l σ α, X tk < xα α P Z l σ α, X tk < xα α =, Q t k,t k+ thank to Lemma 3 On the other hand, we have P ζ k l σ α, X tk > l σ α, X tk xα α = P ζ k l σ α, X tk xα α α Xη bσ α KαX η P inf + t k < t k+ σ σx α η = ψx tk {Xtk > xα α}, + W l σα σx α, X tk xα α η where the inequality come from 5, and the lat equality hold thank to the Markov Property of the Brownian motion, for x α ψx = P inf + b σα Kαx <u σ σx α u + B u l σα σx α, where B t denote a Brownian Motion independent of W t Summarizing P ζ k l σ α, X tk > l σ α ψx tk {Xtk > xα α}, and we have 58 for A k = { X tk > xα α } Step Let u prove that for all x xα α: ψx exp Kα xα α α σ 59 Applying again 55, we have ψx = Kαx + erfc αbσ α σx α + exp b σα Kαxx αb σ α σ x α erfc + Kαx αbσ α σx α =: Ax + Bx 6

27 Since /Kα, and erfcz exp z for all z > we have Ax exp Kαx + αb σ α σ x α exp Kα x α σ On the other hand, for x xα α, and /Kα, it follow x > αb σ α/ + K, o the argument of the function erfc in B i poitive, and then Bx exp b σα Kαxx αb σ α σ x α exp + Kαx αb σ α σ x α So = exp Kαx + αb σ α σ x α exp Kα x α σ ψx = Ax + Bx exp Kα x α σ, and ince the right-hand ide i decreaing on x, we have 59 Step 3 Let u conclude Putting together 58 and 59 we have P Θ α T N k= N k= with γ = xα α α /8σ ψx tk {Xtk > xα α} exp Kα xα α α σ C exp γ P X tk > xα α Proof of Lemma 4 From 4, tandard argument how that L T X4 CT On the other hand, thank to Corollary VI9 on Revuz and Yor 8, p, we have almot urely L T Θα T Θ α X = lim,ε X d X =, ε ε becaue for ε < αb σ α, and T Θ α,,ε X = a Now, ince we can conclude that L T X = L T X {Θα<T } + L T Θ α X {T Θα} = L T X {Θα<T }, L T X = L T X {Θα<T } L T X4 P Θ α < T C exp γ 7

28 54 On the negative moment of the topped increment proce Z t Θα To prove Lemma 5 ee ection 55 below, we need to control the negative moment of the topped increment proce {Z t Θα } t T Thi i the object of the following lemma, that can be ummarize in the following Lemma 5 Let q Let Θ α be the topping time defined in 56 Let u aume max α Moreover, let u aume b > α α when α,, and b > 3 σ q + when α = Then there exit a contant C depending on b, σ, α, T and q but not on, uch that t, T, Z q t Θ α C + x q xitence of Negative Moment Cae α = The proof of the exitence of Negative Moment of Z t Θα ha two part Firt we tudy the quotient X η /Z, and then we proof the main reult of the ection Lemma 53 For α =, and /4K x we have up P T Z X η C 5 8σ bσ/ 5 To prove thi lemma, we need the following auxiliary reult, the proof of which i potponed in Appendix A a a traightforward adaptation of the Lemma 36 in 5 Lemma 54 Aume Hypothei hold, and b > σ /4 Aume alo that /4K x Then, for any γ there exit a contant C depending on the parameter b, K, σ, x, T, and alo on γ, uch that Proof of Lemma 53 We tart by proving up exp X t k k=,,n γσ C up P T Z X η Indeed, if we call = η, and W = W W η, then P Z X η x σ bσ/ γ up exp X t k k=,,n γσ 5 P σ X η W + b σ / + K X η X η P W b σ/ + K X η σ X η F η exp b σ/ + K X η σ X η exp K X η 8σ From here, the bound 5 follow eaily, and then we conclude uing Lemma 54 8

29 Lemma 55 Let Θ be the topping time defined in 56, and q If max /, and b > 3 σ q + 5 Then there exit a contant C depending on b, σ, α, T and q but not on, uch that t, T, Z q t Θ C + x q Proof Let u call W := W W η, and := η By Ito formula Z q t Θ = t Θ bx η x q q Z q+ d qq + t Θ + But, Indeed, Z q+ σ X η + σ W d 53 σ X η + σ W σ X η + σ Z, P a 54 σ X η + σ W =σ X η + σ X η W + σ 4 W =σ X η + Z X η + bx η σ 4 But, thank to the Lipchitz property of b, X η + bx η σ 4 X η + b KX η σ 4 = b σ / + K X η, ince /K, and b σ / > So we have 54 Introducing 54 in 53, and uing bx b Kx, we have Z q t Θ t Θ b t Θ X η x q q Z q+ d + qk Z q+ d 55 qq + t Θ + σ { } Xη + Z d Since X η Z Z q+ X η Z {Z X η /} +, and applying Hölder Inequality for ome ε >, we have Z q t Θ t Θ b x q q Z q+ d 3qq + t Θ + σ + qk Θ Z q Z q+ + C T q+ X /ε ε η P Z X η / ε d d d 9

30 Since b > 3σ q + /, we have 5b σ //8σ > q +, o chooing ε = q/q +, and applying Lemma 53 we have P Z X η / ε C q+, and then Z q t Θ t Θ x q + qk Z q d 3q + t Θ + q σ b Since from the Hypothee, the third term in the right-hand ide i negative, we can conclude thank to Gronwall Lemma Z q+ d + C xitence of Negative Moment Cae α > Lemma 56 For α,, if b > α α σ and /4αKα, there exit γ > uch that up P Z X η exp γ 56 T α Proof Let u call W := W W η, := η, and qx η, W = ασ Xα η W + σx α η W + X η α + b σ α KαX η Notice that for fix x R, qx, i a quadratic polynomial Uing 5, we have P Z X η P qx η, W, X η xα α + P qx η, W, X η xα, where recall, xα = b σ α/kα But P qx η, W, X η xα = P qx, N {Xη x=xη xα}, where N tand for a tandard Gauian random variable A in the Lemma 3, we have a quadratic polynomial in N, it dicriminant i x = σ x α ασ x α α + b σα Kαx = ασ x α b σ α Kαx, o if x xα, < and the quadratic form in N ha not real root, and in particular i non negative almot urely Then P qx η, W, X η xα = On the other hand, P qx η, W, X η xα P N b σα + α Kα x σx α {Xη xα}, x=xη and ince /4αKα we can apply the exponential bound for Gauian tail and get α P qx η, W, X η xα exp Kα x α σ {Xη xα} x=xη We conclude by taking γ = xα α /6σ 3

31 Lemma 57 Let Θ α be the topping time defined in 56 Let u aume for α,, b > α α, and max α, then for all q, there exit a contant C depending on b, σ, α, T and p but not on, uch that t, T, Z q t Θ α C + x q Proof Let u call W := W W η By Ito formula and the Lipchitz property of b, Z q t Θ α t Θα x q q b t Θα d + qk Z q+ + qq + t Θα X η Z q+ d Z q+ σx α η + ασ X α η W d 57 Following the ame idea to prove 54, for all, t we can eaily prove that almot urely σx α η + ασ X α η W σ X α η + ασ X α η Z Introducing thi bound in the previou inequality, we have Z q t Θ α t Θα b t Θα x q q Z q+ d + qk qq + t Θα + σ Z q+ X η Z q+ d { X α η + αx α η Z } d 58 ince for r {, α, α}, we get Xη Z r Xη Z r r α {Z X η + α} α Z q t Θ α t Θα b x q q Z q+ d + α t Θα α qk Z q d qq + + σ αα+ t Θα α α d Z q+ α { } t Θα Xη + C Z q+ + Xα η Z q+ + Xα η Z q+ {Z Xη α}d The lat term in the previou inequality i bounded becaue of the definition of Θ α and the Lemma 56 Indeed, { t Θ α Xη Z q+ C q+ + Xα η Z q+ T + Xα η Z q+ C q+ exp γ } {Z X η α}d X η + X α η + X α η C P Z X η d α 3