arxiv: v2 [math.na] 1 Aug 2016

Transcription

1 STOCHASTIC C-STABILITY AND B-CONSISTENCY OF EXPLICIT AND IMPLICIT MILSTEIN-TYPE SCHEMES WOLF-JÜRGEN BEYN, ELENA ISAAK, AND RAPHAEL KRUSE arxiv: v [math.na] 1 Aug 016 Abstract. This paper focuses on two variants of the Milstein scheme, namely the split-step backward Milstein method and a newly proposed projected Milstein scheme, applied to stochastic differential equations which satisfy a global monotonicity condition. In particular, our assumptions include equations with super-linearly growing drift and diffusion coefficient functions and we show that both schemes are mean-square convergent of order 1. Our analysis of the error of convergence with respect to the mean-square norm relies on the notion of stochastic C-stability and B-consistency, which was set up and applied to Euler-type schemes in [Beyn, Isaak, Kruse, J. Sci. Comp., 015]. As a direct consequence we also obtain strong order 1 convergence results for the split-step backward Euler method and the projected Euler-Maruyama scheme in the case of stochastic differential equations with additive noise. Our theoretical results are illustrated in a series of numerical experiments. 1. Introduction More than four decades ago Grigori N. Milstein proposed a new numerical method for the approximate integration of stochastic ordinary differential equations (SODEs) in [15] (see [16] for an English translation). This scheme is nowadays called the Milstein method and offers a higher order of accuracy than the classical Euler-Maruyama scheme. In fact, G. N. Milstein showed that his method converges with order 1 to the exact solution with respect to the root mean square norm under suitable conditions on the coefficient functions of the SODE while the Euler-Maruyama scheme is only convergent of order 1, in general. In its simplest form, that is for scalar stochastic differential equations driven by a scalar Wiener process W, the Milstein method is given by the recursion X h (t+h) = X h (t)+hf(t,x h (t))+g(t,x h (t)) h W(t) (1) + 1 ( g x g) (t,x h (t)) ( h W(t) h ), where h denotes the step size, h W(t) = W(t+h) W(t) is the stochastic increment, and f and g are the drift and diffusion coefficient functions of the underlying SODE (Equation (3) below shows the SODE in the full generality considered in this paper). Since the derivation of the Milstein method in [15] relies on an iterated application of the Itō formula, the error analysis requires the boundedness and continuity of the coefficient functions f and g and their partial derivatives up to the fourth order. Similar conditions also appear in the standard literature on this topic [9, 17, 18]. 010 Mathematics Subject Classification. 65C30, 65L0. Key words and phrases. stochastic differential equations, global monotonicity condition, splitstep backward Milstein method, projected Milstein method, mean-square convergence, strong convergence, C-stability, B-consistency. 1

2 W.-J. BEYN, E. ISAAK, AND R. KRUSE In more recent publications these conditions have been relaxed: For instance in [10] it is proved that the strong order 1 result for the scheme (1) stays true if the coefficient functions are only two times continuously differentiable with bounded partial derivatives, provided the exact solution has sufficiently high moments and the mapping x ( g x g)(t,x) is globally Lipschitz continuous for every t [0,T]. On the other hand, from the results in [7] it follows that the explicit Euler-Maruyama method is divergent in the strong and weak sense if the coefficient functions are super-linearly growing. Since the same reasoning also applies to the Milstein scheme (1) it is necessary to consider suitable variants in this situation. One possibility to treat super-linearly growing coefficient functions is proposed in [3]. Here the authors combine the Milstein scheme with the taming strategy from [8]. This allows to prove the strong convergence rate 1 in the case of SODEs whose drift coefficient functions satisfy a one-sided Lipschitz condition. The same approach is used in [1], where the authors consider SODEs driven by Lèvy noise. However, both papers still require that the diffusion coefficient functions are globally Lipschitz continuous. This is not needed for the implicit variant of the Milstein scheme considered in [6], where the strong convergence result also applies to certain SODEs with superlinearly growing diffusion coefficient functions. However, the authors only consider scalar SODEs and did not determine the order of convergence. The first result bypassing all these restrictions is found in [5], which deals with an explicit first order method based on a variant of the taming idea. A more recent result based on the taming strategy is also given in [13]. In this paper we propose two further variants of the Milstein scheme which apply to multi-dimensional SODEs of the form (3). First, we follow an idea from [] and study the projected Milstein method which consists of the standard explicit Milstein scheme together with a nonlinear projection onto a sphere whose radius is expanding with a negative power of the step size. The second scheme is a Milstein-type variant of the split-step backward Euler scheme (see [5]) termed split-step backward Milstein method. For both schemes we prove the optimal strong convergence rate 1 in the following sense: Let X: [0,T] Ω R d and X h : {t 0,t 1,...,t N } Ω R d denote the exact solution and its numerical approximation with corresponding step size h. Then, there exists a constant C independent of h such that () max X(t n) X h (t n ) L n {1,...,N} (Ω;R d ) Ch, where t n = nh. For the proof we essentially impose the global monotonicity condition (4) and certain local Lipschitz assumptions on the first order derivatives of the drift and diffusion coefficient functions. For a precise statement of all our assumptions and the two convergence results we refer to Assumption.1 and Theorems. and.3 below. Together with the result on the balanced scheme found in [5], these theorems are the first results which determine the optimal strong convergence rate for some Milstein-type schemes without any linear growth or global Lipschitz assumption on the diffusion coefficient functions and for multi-dimensional SODEs. The remainder of this paper is organized as follows: In Section we introduce the projected Milstein method and the split-step backward Milstein scheme in full detail. We state all assumptions and the convergence results, which are then proved in later sections. Further, we apply the convergence results to SODEs with additive

3 STOCHASTIC C-STAB. AND B-CONS. OF MILSTEIN-TYPE SCHEMES 3 noise for which the Milstein-type schemes coincide with the corresponding Eulertype scheme. The proofs follow the same steps as the error analysis in []. In order to keep this paper as self-contained as possible we briefly recall the notions of C-stability and B- consistency and the abstract convergence theorem from [] in Section 3. Then, in the following four sections we verify that the two considered Milstein-type schemes are indeed stable and consistent in the sense of Section 3. Finally, in Section 8 we report on a couple of numerical experiments which illustrate our theoretical findings. Note that both examples include non-globally Lipschitz continuous coefficient functions, which are not covered by the standard results found in [9, 17].. Assumptions and main results This section contains a detailed description of our assumptions on the stochastic differential equation, under which our strong convergence results hold. Further, we introduce the projected Milstein method and the split-step backward Milstein scheme and we state our main results. Our starting point is the stochastic ordinary differential equation (3) below. We apply the same notation as in [] and we fix d,m N, T (0, ), and a filtered probability space (Ω,F,(F t ),P) satisfying the usual conditions. By X: [0,T] Ω R d we denote a solution to the SODE (3) dx(t) = f(t,x(t))dt+ X(0) = X 0. g r (t,x(t))dw r (t), t [0,T], Here f: [0,T] R d R d stands for the drift coefficient function, while g r : [0,T] R d R d, r = 1,...,m, are the diffusion coefficient functions. By W r : [0,T] Ω R, r = 1,...,m, we denote an independent family of real-valued standard (F t ) -Brownian motions on (Ω,F,P). For a sufficiently large p [, ) the initial condition X 0 is assumed to be an element of the space L p (Ω,F 0,P;R d ). Let us fix some further notation: We write, and for the Euclidean inner product and the Euclidean norm on R d, respectively. Further, we denote by L(R d ) = L(R d,r d ) the set of all bounded linear operators on R d endowed with the matrix norm L(Rd ) induced by the Euclidean norm. For a sufficiently smooth mapping f: [0,T] R d R d and a given t [0,T] we denote by f x (t,x) L(Rd ) the Jacobian matrix of the mapping R d x f(t,x) R d. Having established this we formulate the conditions on the drift and the diffusion coefficient functions: Assumption.1. The mappings f: [0,T] R d R d and g r : [0,T] R d R d, r = 1,...,m, are continuously differentiable. Further, there exist L (0, ) and η ( 1, ) such that for all t [0,T] and x 1,x R d it holds (4) m f(t,x1 ) f(t,x ),x 1 x +η g r (t,x 1 ) g r (t,x ) L x 1 x. In addition, there exists q [, ) such that f x (t,x 1) f x (t,x ) L(Rd ) L( 1+ x 1 + x ) q (5) x1 x

4 4 W.-J. BEYN, E. ISAAK, AND R. KRUSE and, for every r = 1,...,m, gr x (t,x 1) gr x (t,x ) L(Rd ) L( 1+ x 1 + x ) q 3 (6) x 1 x for all t [0,T] and x 1,x R d. Moreover, it holds f t (t,x) ( ) q, L 1+ x g r t (t,x) ( ) q+1 (7) L 1+ x, for all t [0,T], x R d, and all r = 1,...,m. First we note that Assumption.1 is slightly weaker than the conditions imposed in [5, Lemma 4.] in terms of smoothness requirements on the coefficient functions. Further, we recall that Equation (4) is often termed global monotonicity condition in the literature. It is easy to check that Assumption.1 is satisfied (with q = 3) iff and g r and all their first order partial derivatives are globally Lipschitz continuous. However, Assumption.1 includes several SODEs which cannot be treated by the standard results found in [9, 17]. We refer to Section 8 for two more concrete examples. For a possibly enlarged L the following estimates are an immediate consequence of Assumption.1 and the mean value theorem: For all t,t 1,t [0,T] and x,x 1,x R d it holds (8) f(t,x) L ( 1+ x ) q, (9) f x (t,x) L(R d ) L( 1+ x ) q 1, (10) f(t1,x) f(t,x) L ( 1+ x ) q t1 t, (11) and, for all r = 1,...,m, (1) (13) (14) (15) f(t,x 1 ) f(t,x ) L ( 1+ x 1 + x ) q 1 x1 x, g r (t,x) L ( 1+ x ) q+1, gr x (t,x) L(R d ) L( 1+ x ) q 1, g r (t 1,x) g r (t,x) L ( 1+ x ) q+1 t 1 t, g r (t,x 1 ) g r (t,x ) L ( 1+ x 1 + x ) q 1 x 1 x. Thus, Assumption.1 implies [, Assumption.1] and all results of that paper also hold true in the situation considered here. Note that in this paper we use the weights (1+ x ) p instead of 1+ x p as in []. For p 0 this makes no difference, however in condition (6) we may have p = q 3 < 0 if q < 3, so that Lipschitz constants actually decrease at infinity. In the following it will be convenient to introduce the abbreviation (16) g r1,r (t,x) := gr1 (t,x), t [0,T], x R d, x (t,x)gr for r 1,r = 1,...,m. As above, one easily verifies under Assumption.1 that the mappings g r1,r satisfy (for a possibly larger L) the polynomial growth bound g r 1,r (t,x) ( ) q (17) L 1+ x as well as the local Lipschitz bound g r 1,r (t,x 1 ) g r1,r (t,x ) ( L 1+ x1 + x ) q 1 (18) x1 x

5 STOCHASTIC C-STAB. AND B-CONS. OF MILSTEIN-TYPE SCHEMES 5 for all x,x 1,x R d, t [0,T], and r 1,r = 1,...,m. For this conclusion to hold in case q < 3, it is essential to use the modified weight function in (6). We say that an almost surely continuous and (F t ) -adapted stochastic process X: [0,T] Ω R d is a solution to (3) if it satisfies P-almost surely the integral equation t t (19) X(t) = X 0 + f(s,x(s))ds+ g r (s,x(s))dw r (s) 0 for all t [0,T]. It is well-known that Assumption.1 is sufficient to ensure the existence of a unique solution to (3), see for instance [11], [14, Chap..3] or the SODE chapter in [0, Chap. 3]. In addition, the exact solution has finite p-th moments, that is (0) sup X(t) Lp (Ω;R d ) <, if the following global coercivity condition is satisfied: There exist C (0, ) and p [, ) such that p 1 (1) f(t,x),x + g r (t,x) C ( 1+ x ) for all x R d, t [0,T]. A proof is found, for example, in [14, Chap..4]. For the formulation of the numerical methods we recall the following terminology from []: By h (0,T] we denote an upper step size bound. Then, for every N N we say that h = (h 1,...,h N ) (0,h] N is a vector of (deterministic) step sizes if N i=1 h i = T. Every vector of step sizes h induces a set of temporal grid points T h given by T h := { t n := where = 0. For short we write n i=1 0 } h i : n = 0,...,N, h := max i {1,...,N} h i for the maximal step size in h. Moreover, we recall from [9, 17] the following notation for the stochastic increments: Let t,s [0,T] with s < t. Then we define () for r {1,...,m} and, similarly, (3) I s,t (r 1,r ) := t I s,t (r) := t s τ1 s s dw r (τ), dw r1 (τ )dw r (τ 1 ), where r 1,r {1,...,m}. The joint family of the iterated stochastic integrals (I s,t (r 1,r ) )m r 1,r =1 is not easily generated on a computer. Besides special cases such as commutative noise one relies on an additional approximation method from e.g. [3, 1, 4]. We also refer to the corresponding discussion in [9, Chap. 10.3]. The first numerical scheme, which we study in this paper, is an explicit onestep scheme and termed projected Milstein method (PMil). It is the Milstein-type counterpart of the projected Euler-Maruyama method form [] and consists of the

6 6 W.-J. BEYN, E. ISAAK, AND R. KRUSE standard Milstein scheme and a projection onto a ball in R d whose radius is expanding with a negative power of the step size. To be more precise, let h (0,h] N, N N, be an arbitrary vector of step sizes with upper step size bound h = 1. For a given parameter α (0, ) the PMil method is determined by the recursion (4) X PMil h X PMil h (t i ) = min ( 1,h α i X PMil h (t i 1 ) 1 ) X PMil h (t i 1 ), (t i ) = X PMil h (t i )+h i f(t i 1,X PMil h (t i ))+ + r 1,r =1 g r1,r (t i 1,X PMil h g r (t i 1,X PMil h (t i ))I ti 1,ti (r) (t i ))I ti 1,ti, for 1 i N, (r,r 1) where Xh PMil (0) := X 0. The results of Section 4 indicate that the parameter value for α is optimally chosen by setting α = 1 (q 1) in dependence of the growth rate q appearing in Assumption.1. One aim of this paper is the proof of the following strong convergence result for the PMil method. It follows directly from Theorems 4.4 and 5.1 as well as Theorem 3.5. Theorem.. Let Assumption.1 be satisfied with polynomial growth rate q [, ). If the exact solution X to (3) satisfies sup τ [0,T] X(τ) L 8q 6 (Ω;R d ) <, then the projected Milstein method (4) with parameter value α = 1 (q 1) and with arbitrary upper step size bound h (0,1] is strongly convergent of order γ = 1. Next, we come to the second numerical scheme, which is called split-step backward Milstein method (SSBM). For a suitable upper step size bound h (0,T] and a given vector of step sizes h = (h 1,...,h N ) (0,h] N, N N, this method is defined by setting Xh SSBM (0) = X 0 and by the recursion (5) X SSBM h (t i ) = X SSBM h Xh SSBM (t i ) = X SSBM + (t i 1 )+h i f(t i,x SSBM h (t i )), h (t i )+ r 1,r =1 g r (t i,x SSBM h (t i ))I ti 1,ti (r) g r1,r (t i,x SSBM h (t i ))I ti 1,ti (r,,r 1) for every i = 1,...,N. Let us note that the recursion defining the SSBM method evaluates the diffusion coefficient functions g r at time t i in the i-th step. This phenomenon was already apparent in the definition of the split-step backward Euler method in []. It turns out that by this modification we avoid some technical issues in the proofs as condition (6) is applied to f and g r, r = 1,...,m, simultaneously at the same point t [0,T] in time. Compare also with the inequality (50) further below. It is shown in Section 6 that the SSBM scheme is a well-defined stochastic onestep method under Assumption.1. The second main result of this paper is the proof of the following strong convergence result:

7 STOCHASTIC C-STAB. AND B-CONS. OF MILSTEIN-TYPE SCHEMES 7 Theorem.3. Let Assumption.1 be satisfied with L (0, ) and q [, ). In addition, we assume that there exist η 1 (1, ) and η (0, ) such that it holds (6) m f(t,x1 ) f(t,x ),x 1 x +η1 g r (t,x 1 ) g r (t,x ) +η m r 1,r =1 g r 1,r (t,x 1 ) g r1,r (t,x ) L x 1 x for all x 1,x R d. If the solution X to (3) satisfies sup τ [0,T] X(τ) L 6q 4 (Ω;R d ) <, then the split-step backward Milstein method (5) with arbitrary upper step size bound h (0,max(L 1, η η 1 )) is strongly convergent of order γ = 1. As we show below this theorem follows directly from Theorem 3.5 together with Theorems 6.3 and 7.1. Note that (6) is more restrictive than the global monotonicity condition (4) if the mappings g r1,r are not globally Lipschitz continuous for all r 1,r = 1,...,m. In the remainder of this section we briefly summarize the corresponding convergence results in the case of stochastic differential equations with additive noise, that is if the mappings g r, r = 1,...,m, do not depend explicitly on the state of X. In this case it is well-known that Milstein-type schemes coincide with their Euler-type counterparts. To be more precise, we consider the solution X: [0,T] Ω R d to an SODE of the form dx(t) = f(t,x(t))dt+ g r (t)dw r (t), t [0,T], (7) X(0) = X 0. In this case, the conditions on the drift coefficient function f: [0,T] R d R d and the diffusion coefficient functions g r : [0,T] R d, r = 1,...,m, in Assumption.1 simplify to Assumption.4 (Additive noise). The coefficient functions f: [0,T] R d R d and g r : [0,T] R d, r = 1,...,m, are continuously differentiable, and there exist constants L (0, ), q [, ) such that for all t [0,T] and x,x 1,x R d the following properties hold f(t,x1 ) f(t,x ),x 1 x L x1 x, f t (t,x) L ( 1+ x ) q, f x (t,x 1) f x (t,x ) L(R d ) L( 1+ x 1 + x ) q x1 x. Under this assumption it directly follows that g r1,r 0 for all r 1,r = 1,...,m for the coefficient functions defined in (16). Consequently, the PMil method and the SSBM scheme coincide with the PEM method and the SSBE scheme from [], respectively. Let us note that Assumption.4 implies the global coercivity condition (1) for every p [, ). Consequently, under Assumption.4 the exact solution to (7) has finite p-th moments for every p [, ). From this and Theorems. and.3 we directly obtain the following convergence result: Corollary.5. Let Assumption.4 be satisfied with L (0, ) and q [, ). Then it holds that

8 8 W.-J. BEYN, E. ISAAK, AND R. KRUSE (i) the projected Euler-Maruyama method with α = 1 (q 1) and arbitrary upper step size bound h (0,1] is strongly convergent of order γ = 1. (ii) the split-step backward Euler method with arbitrary upper step size bound h (0,L 1 ) is strongly convergent of order γ = A reminder on stochastic C-stability and B-consistency In this section we give a brief overview of the notions of stochastic C-stability and B-consistency introduced in []. We also state the abstract convergence theorem, which, roughly speaking, can be summarized by stoch. C-stability + stoch. B-consistency Strong convergence. We first recall some additional notation from []: For an arbitrary upper step size bound h (0,T] we define the set T := T(h) [0,T) (0,h] to be T(h) := { (t,) [0,T) (0,h] : t+ T }. Further, for a given vector of step sizes h (0,h] N, N N, we denote by G (T h ) the space of all adapted and square integrable grid functions, that is G (T h ) := { Z: T h Ω R d : Z(t n ) L (Ω,F tn,p;r d ) for all n = 0,1,...,N }. The next definition describes the abstract class of stochastic one-step methods which we consider in this section. Definition 3.1. Let h (0,T] be an upper step size bound and Ψ: R d T Ω R d be a mapping satisfying the following measurability and integrability condition: For every (t,) T and Z L (Ω,F t,p;r d ) it holds (8) Ψ(Z,t,) L (Ω,F t+,p;r d ). Then, for every vector of step sizes h (0,h] N, N N, we say that a grid function X h G (T h ) is generated by the stochastic one-step method (Ψ,h,ξ) with initial condition ξ L (Ω,F 0,P;R d ) if (9) X h (t i ) = Ψ(X h (t i 1 ),t i 1,h i ), 1 i N, X h (t 0 ) = ξ. We call Ψ the one-step map of the method. For the formulation of the next definition we denote by E[Y F t ] the conditional expectation of a random variable Y L 1 (Ω;R d ) with respect to the sigma-field F t. Note that if Y is square integrable, then E[Y F t ] coincides with the orthogonal projection onto the closed subspace L (Ω,F t,p;r d ). By (id E[ F t ]) we denote the associated projector onto the orthogonal complement. Definition 3.. A stochastic one-step method (Ψ, h, ξ) is called stochastically C- stable (with respect to the norm in L (Ω;R d )) if there exist a constant C stab and a parameter value ν (1, ) such that for all (t,) T and all random variables Y,Z L (Ω,F t,p;r d ) it holds E [ Ψ(Y,t,) Ψ(Z,t,) Ft ] L (Ω;R d ) (30) +ν ( id E[ F t ] )( Ψ(Y,t,) Ψ(Z,t,) ) L (Ω;R d ) ( 1+C stab ) Y Z L (Ω;R d ).

9 STOCHASTIC C-STAB. AND B-CONS. OF MILSTEIN-TYPE SCHEMES 9 A first consequence of the notion of stochastic C-stability is the following a priori estimate: Let (Ψ, h, ξ) be a stochastically C-stable one-step method. If there exists a constant C 0 such that for all (t,) T it holds (31) (3) E [ ] Ψ(0,t,) F L t C (Ω;R d ) 0, ( id E [ ]) F t Ψ(0,t,) L (Ω;R d ) C 0 1, then there exists a positive constant C with max X h(t n ) L n {0,...,N} (Ω;R d ) e CT( )1 ξ L (Ω;R d ) +CC 0T, for every vector of step sizes h (0,h] N, N N, where X h denotes the grid function generated by (Ψ,h,ξ) with step sizes h. A proof for this result is found in [, Cor. 3.6]. Definition 3.3. A stochastic one-step method (Ψ, h, ξ) is called stochastically B- consistent of order γ > 0 to (3) if there exists a constant C cons such that for every (t,) T it holds [ ] (33) E X(t+) Ψ(X(t),t,) Ft L (Ω;R d ) C cons γ+1 and (34) ( id E[ F t ] )( X(t+) Ψ(X(t),t,) ) L (Ω;R d ) C cons γ+1, where X: [0,T] Ω R d denotes the exact solution to (3). Finally, it remains to give our definition of strong convergence. Definition 3.4. A stochastic one-step method (Ψ, h, ξ) converges strongly with order γ > 0 to the exact solution of (3) if there exists a constant C such that for every vector of step sizes h (0,h] N, N N, it holds max n {0,...,N} X h (t n ) X(t n ) L (Ω;R d ) C h γ. Here X denotes the exact solution to (3) and X h G (T h ) is the grid function generated by (Ψ,h,ξ) with step sizes h (0,h] N. We close this section with the following abstract convergence theorem, which is proved in [, Theorem 3.7]. Theorem 3.5. Let the stochastic one-step method (Ψ, h, ξ) be stochastically C- stable and stochastically B-consistent of order γ > 0. If ξ = X 0, then there exists a constant C depending on C stab, C cons, T, h, and ν such that for every vector of step sizes h (0,h] N, N N, it holds max n {0,...,N} X(tn ) X h (t n ) L (Ω;R d ) C h γ, where X denotes the exact solution to (3) and X h the grid function generated by (Ψ, h, ξ) with step sizes h. In particular, (Ψ, h, ξ) is strongly convergent of order γ.

10 10 W.-J. BEYN, E. ISAAK, AND R. KRUSE 4. C-stability of the projected Milstein method In this section we prove that the projected Milstein (PMil) method defined in (4) is stochastically C-stable. Throughout this section we assume that Assumption.1 is satisfied with growth rate q [, ). First, we choose an arbitrary upper step size bound h (0,1] and a parameter value α (0, ). Later it will turn out to be optimal to set α = 1 (q 1) in dependence of the growth q in Assumption.1. For the definition of the one-step map of the PMil method it is convenient to introduce the following short hand notation: For every (0,h], we denote the projection of x R d onto the ball of radius α by (35) x := min(1, α x 1 )x. Then, the one-step map Ψ PMil : R d T Ω R d is given by (36) Ψ PMil (x,t,) := x +f(t,x )+ g r (t,x )I t,t+ (r) + r 1,r =1 g r1,r (t,x )I t,t+ (r,r 1) for every x R d and (t,) T. Recall () and (3) for the definition of the stochastic increments. First, we check that the PMil method is a stochastic one-step method in the sense of Definition 3.1. At the same time we verify that the one step map satisfies conditions (31) and (3). Proposition 4.1. Let the functions f and g r, r = 1,...,m, satisfy Assumption.1 with L (0, ) and let h (0,1]. For every initial value ξ L (Ω;F 0,P;R d ) and for every α (0, ) it holds that (Ψ PMil,h,ξ) is a stochastic one-step method. In addition, there exists a constant C 0 only depending on L and m such that [ ] (37) E Ψ PMil (0,t,) F L t (Ω;R d ) C 0, (38) ( id E[ F t ] ) Ψ PMil (0,t,) L (Ω;R d ) C 0 1 for all (t,) T. Proof. We first verify thatψ PMil satisfies (8). For this let us fix arbitrary(t,) T and Z L (Ω,F t,p;r d ). Then, the continuity and boundedness of the mapping R d x x = min(1, α x 1 )x R d yields Z L (Ω,F t,p;r d ). Consequently, by the smoothness of the coefficient functions and by (8), (1), and (17) it follows that f(t,z ), g r1 (t,z ), g r1,r (t,z ) L (Ω,F t,p;r d ) for every r 1,r = 1,...,m. Therefore, Ψ PMil (Z,t,): Ω R d is an F t+ /B(R d )- measurable random variable satisfying condition (8). It remains to show (37) and (38). From (8) we get immediately that E [ Ψ PMil (0,t,) F t ] L (Ω;R d ) = f(t,0) L.

11 STOCHASTIC C-STAB. AND B-CONS. OF MILSTEIN-TYPE SCHEMES 11 Next, recall that the stochastic increments (I t,t+ (r) ) m and (I t,t+ (r 1,r ) )m r 1,r =1 are pairwise uncorrelated. Therefore, we obtain that ( id E[ F t ] ) Ψ PMil (0,t,) L (Ω;R d ) = = g r (t,0)i t,t+ (r) + g r (t,0) + r 1,r =1 r 1,r =1 g r1,r (t,0)i t,t+ (r 1,r ) L (Ω;R d ) g r 1,r (t,0) L m + 1 L m, where the last step follows from (1) and (17). Since h 1 this verifies (38). The next result is concerned with the projection onto the ball of radius α. The proof is found in [, Lem. 6.]. Lemma 4.. For every α (0, ) and (0,1] the mapping R d x x R d defined in (35) is globally Lipschitz continuous with Lipschitz constant 1. In particular, it holds x (39) 1 x x1 x for all x 1,x R d. The following inequality (40) follows from the global monotonicity condition (4) and plays an import role in the stability analysis of the PMil method. The proof is given in [, Lem. 6.3]. Lemma 4.3. Let the functions f and g r, r = 1,...,m, satisfy Assumption.1 with L (0, ), q [, ), and η ( 1, ). Consider the mapping Rd x x R d 1 defined in (35) with α (0, (q 1) ] and (0,1]. Then there exists a constant C only depending on L such that (40) for all x 1,x R d. x 1 x +(f(t,x 1) f(t,x )) +η (1+C) x 1 x g r (t,x 1) g r (t,x )) The next theorem verifies that the PMil method is stochastically C-stable. Theorem 4.4. Let the functions f and g r, r = 1,...,m, satisfy Assumption.1 with L (0, ), q [, ), and η ( 1, ). Further, let h (0,1]. Then, for every ξ L (Ω,F 0,P;R d ) the projected Milstein method (Ψ PMil,h,ξ) with α = 1 (q 1) is stochastically C-stable. Proof. Let (t,) T and Y,Z L (Ω,F t,p;r d ) be arbitrary. By recalling (36) we obtain E [ Ψ PMil (Y,t,) Ψ PMil (Z,t,) F t ] = Y +f(t,y ) (Z +f(t,z )) and ( id E[ Ft ] )( Ψ PMil (Y,t,) Ψ PMil (Z,t,) ) = ( g r (t,y ) g r (t,z ) ) I t,t+ (r) + r 1,r =1 ( g r 1,r (t,y ) g r1,r (t,z ) ) I t,t+ (r,r 1).

12 1 W.-J. BEYN, E. ISAAK, AND R. KRUSE In order to verify (30) with ν = η (1, ) let us note that the stochastic increments are pairwise uncorrelated and independent of Y and Z. Hence it follows E [ Ψ PMil (Y,t,) Ψ PMil (Z,t,) F t ] L (Ω;R d ) +ν ( id E[ F t ] )( Ψ PMil (Y,t,) Ψ PMil (Z,t,) ) L (Ω;R d ) = Y +f(t,y ) (Z +f(t,z )) L (Ω;R d ) +ν g r (t,y ) g r (t,z ) L (Ω;R d ) + 1 ν m r 1,r =1 g r 1,r (t,y ) g r1,r (t,z ) L (Ω;R d ). An application of Lemma 4.3 with ν = η shows that the first two terms are dominated by E[ Y +f(t,y ) (Z +f(t,z )) +ν g r (t,y ) g r (t,z ) ] (1+C) Y Z L (Ω;R d ). In addition, applications of (18) and (39) yield g r 1,r (t,x 1 (t,x ) gr1,r ) ( L 1+ x 1 + x ) q 1 x 1 x L ( 1+ α) q 1 x 1 x, where we made use of the fact that x 1, x α. Since α(q 1) = 1 and (0,1] it follows 1 (1+ α ) q 1 3 q 1 and, therefore, ν m r 1,r =1 g r1,r (t,y ) g r1,r (t,z ) L (Ω;R d ) 3(q 1) νm L Y Z L (Ω;R d ). This completes the proof. 5. B-consistency of the projected Milstein method In this section we show that the PMil method is stochastically B-consistent of order γ = 1. To be more precise, we prove the following result: Theorem 5.1. Let f and g r, r = 1,...,m, satisfy Assumption.1 with L (0, ) and q [, ). Let h (0,1] be arbitrary. If the exact solution X to (3) satisfies sup τ [0,T] X(τ) L 8q 6 (Ω;R d ) <, then the projected Milstein method (Ψ PMil,h,X 0 ) with α = 1 (q 1) is stochastically B-consistent of order γ = 1. In preparation for the proof of Theorem 5.1 we introduce several more technical lemmas. The first is cited from [, Lemma 6.5]. It formalizes a method of proof already found in [5, Theorem.]. Lemma 5.. For arbitrary α (0, ) and (0,1] consider the mapping R d x x R d defined in (35). Let L (0, ), κ [1, ) and let ϕ: R d R d be a measurable mapping which satisfies ϕ(x) L ( 1+ x κ)

13 STOCHASTIC C-STAB. AND B-CONS. OF MILSTEIN-TYPE SCHEMES 13 for all x R d. For some p (, ) let Y L pκ (Ω;R d ). Then there exists a constant C only depending on L and p with ϕ(y) ϕ(y ) L (Ω;R d ) C( 1+ Y 1 pκ ) 1 L pκ (Ω;R d ) α(p )κ. The proof of consistency also depends on the Hölder continuity of the exact solution to (3) with respect to the norm in L p (Ω;R d ) for some p [, ). A proof is given in [, Proposition 5.4]. Proposition 5.3. Let f and g r, r = 1,...,m, satisfy Assumption.1 with L (0, ) and q [, ). For every p [, ) there exists a constant C = C(L,q,p) such that (41) X(t 1 ) X(t ) Lp (Ω;R d ) C( 1+ sup X(t) q L pq (Ω;R )) t1 t d 1 holds for all t 1,t [0,T] and for every solution X of (3) satisfying sup X(t) L pq (Ω;R d ) <. The next auxiliary result combines the Hölder estimates with growth functions. Proposition 5.4. Let q 1 0,q > 0 and consider an R d -valued process X(t),t [0,T] satisfying (41) for pq = (qq +q 1 ) and sup X(t) L pq (Ω;R d ) <. Then there exists a constant C such that for all 0 t 1 t T (4) ( 1+ X(t 1 ) + X(t ) ) q 1 X(t1 ) X(t ) q L (Ω;R) C ( 1+ sup X(t) pq L pq (Ω;R )) t1 t d q. Proof. We apply a Hölder estimate with arbitrary ν > 1, ν = ν ν 1 ( 1+ X(t 1 ) + X(t ) ) q 1 X(t1 ) X(t ) q L (Ω;R) and use (41), C 1+ X(t 1 ) + X(t ) q1 L ν q 1(Ω;R) X(t ) X(t 1 ) q L νq (Ω;R d ) C ( )( 1+ sup X(t) q1 L ν q 1(Ω;Rd ) 1+ sup X(t) qq L νqq (Ω;R )) t1 t d q. The norms are balanced if we choose ν q 1 = qνq which leads to ν = 1 + q1 qq and νqq = (qq +q 1 ). This shows our assertion for q 1 > 0. In case q 1 = 0 it is enough to just apply (41) with p = q. The following lemma is quoted from [, Lemma 5.5]. Lemma 5.5. Let Assumption.1 be satisfied by f and g r, r = 1,...,m, with L (0, ) and q [, ). Further, let the exact solution X to the SODE (3) satisfy sup X(t) L 4q (Ω;R d ) <. Then, there exists a constant C such that for all t 1,t,s [0,T] with 0 t 1 s t T it holds t f(τ,x(τ)) f(s,x(t1 )) L (Ω;R d ) dτ t 1 C ( 1+ sup X(t) q 1 L 4q (Ω;R )) t1 t d 3. The order of convergence indicated by Lemma 5.5 can be increased if we insert the conditional expectation with respect to the σ-field F t1 :

14 14 W.-J. BEYN, E. ISAAK, AND R. KRUSE Lemma 5.6. Let Assumption.1 be satisfied by f and g r, r = 1,...,m with L (0, ) and q [, ). Further, let the exact solution X to the SODE (3) satisfy sup X(t) L 6q 4 (Ω;R d ) <. Then, there exists a constant C such that for all t 1,t,s [0,T] with 0 t 1 s t T it holds t [ ] E f(τ,x(τ)) f(s,x(t1 )) F L t1 (Ω;R d ) dτ t 1 C ( 1+ sup X(t) 3q L 6q 4 (Ω;R d ) ) t1 t. Proof. Since E[Y F t1 ] L (Ω;R d ) Y L (Ω;R d ) for all Y L (Ω;R d ) the integrand is estimated by E [ f(τ,x(τ)) f(s,x(t1 )) F t1 ] L (Ω;R d ) f(τ,x(τ)) f(s,x(τ)) L (Ω;R d ) + E [ f(s,x(τ)) f(s,x(t1 )) F t1 ] L (Ω;R d ) for every τ [t 1,t ]. From (10) it follows that f(τ,x(τ)) f(s,x(τ)) L (Ω;R d ) L 1+ X(τ) q L q (Ω;R) τ s (43) C ( 1+ sup X(t) q L q (Ω;R )) t t d 1, which after integrating over τ, yields the desired estimate since q 6q 4 for q 1. Next, from the mean value theorem we obtain f(s,x(τ)) f(s,x(t 1 )) = f x (s,x(t 1)) ( X(τ) X(t 1 ) ) +R f, where the remainder term R f is given by R f = 1 0 ( f ( s,x(t1 )+ρ(x(τ) X(t 1 )) ) f ) x x (s,x(t 1)) dρ ( X(τ) X(t 1 ) ). Using the SODE (3) we obtain [ f E x (s,x(t 1)) ( X(τ) X(t 1 ) ) ] Ft1 [ f τ = E x (s,x(t 1)) After taking the L -norm and inserting (8) and (9) we arrive at [ f E x (s,x(t 1)) ( X(τ) X(t 1 ) ) ] L Ft1 (Ω;R d ) (44) τ L(1+ X(t1 ) ) q 1 ( ) q L 1+ X(σ) L (Ω;R) dσ t 1 C ( 1+ sup X(t) q 1 L 4q (Ω;R )) τ t1. d f(σ,x(σ))dσ F t1 ]. t 1 Hence, we also obtain the desired estimate for this term after integrating over τ. Finally, we have to estimate the L -norm of the remainder term R f. For this we make use of (5) and get (45) R f 1 0 L ( 1+ X(t 1 )+ρ(x(τ) X(t 1 )) + X(t 1 ) ) q dρ X(τ) X(t1 ) C ( 1+ X(t1 ) + X(τ) ) q X(τ) X(t1 )

15 STOCHASTIC C-STAB. AND B-CONS. OF MILSTEIN-TYPE SCHEMES 15 for a constant C only depending on L and q. Applying Proposition 5.4 with q =,q 1 = q yields (qq +q 1 ) = 6q 4 and therefore, (46) R f L (Ω;R d ) C ( 1+ sup X(t) 3q L 6q 4 (Ω;R )) τ t1. d The next lemma contains the corresponding estimate for the stochastic integral. Lemma 5.7. Let Assumption.1 be satisfied by f and g r, r = 1,...,m with L (0, ) and q [, ). Further, let the exact solution X to (3) satisfy sup X(t) L 6q 4 (Ω;R d ) <. Then, there exists a constant C such that for all r = 1,...,m and t 1,t,s [0,T] with 0 t 1 s t T it holds t g r (τ,x(τ)) g r (s,x(t 1 ))dw r (τ) g r,r (s,x(t 1 ))I t1,t (r,r) t L 1 (Ω;R d ) C ( 1+ sup X(t) 3q r =1 L 6q 4 (Ω;R d )) t1 t 3. Proof. Let us fixr = 1,...,m arbitrary. We first consider the square of thel -norm and by recalling () and (3) we get [ t ( ) E g r (τ,x(τ)) g r (s,x(t 1 )) g r,r (s,x(t 1 ))I t1,τ (r ) dw r (τ) ] t 1 r =1 t [ g = r (τ,x(τ)) g r (s,x(t 1 )) g r,r (s,x(t 1 ))I t1,τ ]dτ t 1 E by an application of the Itō isometry. Thus, the assertion is proved if there exists a constant C independent of τ, t 1, t, and s such that Γ(τ) := g r (τ,x(τ)) g r (s,x(t 1 )) g r,r (s,x(t 1 ))I t1,τ r =1 r =1 C ( 1+ sup X(t) 3q L 6q 4 (Ω;R )) t1 t d for every τ [t 1,t ]. For this we first estimate Γ(τ) by (r ) (r ) L (Ω;R d ) Γ(τ) g r (τ,x(τ)) g r (s,x(τ)) L (Ω;R d ) + g r (s,x(τ)) g r (s,x(t 1 )) g r,r (s,x(t 1 ))I t1,τ. L (Ω;R d ) r =1 In the same way as in (43) one shows for the first term g r (τ,x(τ)) g r (s,x(τ)) L (Ω;R d ) C( 1+ sup X(t) q+1 L q+1 (Ω;R )) t t d 1 and notes q +1 6q 4. Next, we again apply the mean value theorem g r (s,x(τ)) g r (s,x(t 1 )) = gr x (s,x(t 1)) ( X(τ) X(t 1 ) ) +R g, where this time the remainder term R g is given by R g := 1 0 ( g r x ( s,x(t1 )+ρ(x(τ) X(t 1 )) ) ) gr x (s,x(t 1)) dρ ( X(τ) X(t 1 ) ). (r )

16 16 W.-J. BEYN, E. ISAAK, AND R. KRUSE Using the condition (6) we get R g C ( 1+ X(t 1 ) + X(τ) ) q 1 X(τ) X(t1 ), where q 1 = (q 3)+. Therefore, Proposition 5.4 applies with q = and leads to R g L (Ω;R d ) C ( 1+ sup X(t) 3q L 6q 4 (Ω;R )) τ t1, d since (qq + q 1 ) = max(5q 3,4q) 6q 4 for q. It remains to give a corresponding estimate for Γ (τ) := gr x (s,x(t 1)) ( X(τ) X(t 1 ) ) g r,r (s,x(t 1 ))I t1,τ (r ). L (Ω;R d ) After inserting (19) we finally arrive at the two terms τ Γ (τ) gr x (s,x(t 1)) f(σ,x(σ))dσ t L 1 (Ω;R d ) τ + gr x (s,x(t 1)) g r (σ,x(σ))dw r (σ) g r,r (s,x(t 1 ))I t1,τ (r ). t L 1 (Ω;R d ) r =1 Using (13), the first term is estimated analogously to (44), τ gr x (s,x(t 1)) f(σ,x(σ))dσ C( 1+ sup X(t) 3q 1 t 1 L (Ω;R d ) L 3q 1 (Ω,R )) τ t1. d r =1 For the second term we insert (16) and () and obtain from Itō s isometry τ gr x (s,x(t 1)) g r (σ,x(σ))dw r (σ) g r,r (s,x(t 1 ))I t1,τ (r ) r =1 t L 1 (Ω;R d ) τ g r = r =1 t 1 x (s,x(t 1)) ( g r (σ,x(σ)) g r (s,x(t 1 )) ) dw r (σ) ( τ = gr x (s,x(t 1)) ( g r (σ,x(σ)) g r (s,x(t 1 )) ) ) 1 dσ. L (Ω;R d ) r =1 t 1 Now, it follows from (14) and (15) that g r (σ,x(σ)) g r (s,x(t 1 )) g r (σ,x(σ)) g r (σ,x(t 1 )) + g r (σ,x(t 1 )) g r (s,x(t 1 )) L (Ω;R d ) L ( 1+ X(t 1 ) + X(σ) ) q 1 X(σ) X(t 1 ) +L ( 1+ X(t 1 ) ) q+1 σ s. Hence, the growth estimate (13) and Proposition 5.4 with q 1 = q 1,q = 1 yield gr x (s,x(t 1)) ( g r (σ,x(σ)) g r (s,x(t 1 )) ) L (Ω;R d ) L ( 1+ X(t 1 ) ) q 1 g r (σ,x(σ)) g r (s,x(t 1 )) L (Ω,R) L ( 1+ X(t 1 ) + X(σ) ) q 1 X(σ) X(t1 ) L (Ω,R) +L ( 1+ X(t 1 ) ) q L (Ω,R) σ s (1+ C( sup X(t) q 1 ) L 4q (Ω,R d ) σ t1 1 + ( 1+ sup X(t) q ) ) L q (Ω,R d ) σ s.

17 STOCHASTIC C-STAB. AND B-CONS. OF MILSTEIN-TYPE SCHEMES 17 To sum up, we have shown Γ (τ) C ( 1+ sup X(t) q 1 L 4q (Ω;R )) t1 t d. Since 4q 6q 4, this completes the proof. The proof shows that it is sufficient to have bounds for moments of order max(5q 3,4q) instead of 6q 4. However, in view of the weaker estimate in Lemma 5.6, this does not improve the result of Theorem 5.1. Now we are well-prepared for the proof of Theorem 5.1. Proof of Theorem 5.1. We first verify (33) with γ = 1 for the PMil method. For this, let (t,) T be arbitrary. After inserting (19) and (36) we obtain [ E X(t+) Ψ PMil (X(t),t,) ] Ft L (Ω;R d ) [ t+ ] = E X(t)+ f(τ,x(τ))dτ X (t) f(t,x L (t)) F t t (Ω;R d ) X(t) X (t) L (Ω;R d ) + f(t,x(t)) f(t,x (t)) L (Ω;R d ) + t+ t [ ] E f(τ,x(τ)) f(t,x(t)) Ft L (Ω;R d ) dτ. By applying Lemma 5. with ϕ = id, κ = 1, p = 8q 6, and α = 1 (q 1) X(t) X (t) L (Ω;R d ) C( 1+ X(t) 4q 3 L 8q 6 (Ω;R )), d we obtain since 1 α(p )κ =. Similarly, we estimate the second term by Lemma 5. with ϕ = f(t, ), κ = q, and p = 6 4 q. Since in this case 1 α(p )κ = 1 we get (47) f(t,x(t)) f(t,x (t)) L (Ω;R d ) C( 1+ X(t) 3q L 6q 4 (Ω;R d )). The last term is estimated by Lemma 5.6 with t 1 = s = t and t = t+, t+ t [ ] E f(τ,x(τ)) f(t,x(t)) Ft L (Ω;R d ) dτ C ( 1+ sup X(t) 3q L 6q 4 (Ω;R )). d This completes the proof of (33). For the proof of (34) we first insert (19) and (36). Then, in the same way as above we obtain the following four terms (id E[ Ft ]) ( X(t+) Ψ PMil (X(t),t,) ) L (Ω;R d ) + ( g r (t,x(t)) g r (t,x (t)) ) I t,t+ (r) r,r =1 L (Ω;R d ) ( g r,r (t,x(t)) g r,r (t,x (t)) ) I t,t+ + (id E[ F t ]) + t+ t t+ t f(τ,x(τ))dτ L (Ω;R d ) (r,r ) L (Ω;R d ) g r (τ,x(τ)) g r (t,x(t))dw r (τ) g r,r (t,x(t))i t,t+. L (Ω;R d ) r =1 (r,r)

18 18 W.-J. BEYN, E. ISAAK, AND R. KRUSE Since the stochastic increment I t,t+ (r) is independent of F t it directly follows that ( g r (t,x(t)) g r (t,x (t)) ) I t,t+ (r) L (Ω;R d ) = 1 g r (t,x(t)) g r (t,x (t)) L (Ω;R d ). Then, we apply Lemma 5. with ϕ = g r (t, ), κ = q+1 16, and p = 10 q+1. As above, this yields 1 α(p )κ = 1 and we get 1 g r (t,x(t)) g r (t,x (t)) L (Ω;R d ) C( 1+ sup X(t) 5 q 3 L 5q 3 (Ω;R d ) for every r = 1,..., m. In the same way we obtain for the second term ( g r,r (t,x(t)) g r,r (t,x (t)) ) I t,t+ (r,r ) L (Ω;R d ) = 1 g r,r (t,x(t)) g r,r (t,x (t)) L (Ω;R d ). ) 3 Then, a further application of Lemma 5. with ϕ = g r,r (t, ), κ = q, and p = 4 q gives g r,r (t,x(t)) g r,r (t,x (t)) C( L 1+ sup X(t) q 1 ) 3 (Ω;R d ) L 4q (Ω;R d ) for every r,r = 1,...,m, since in this case 1 α(p )κ = 1. Next, since f(t,x(t)) is F t -measurable it follows for the third term that t+ (id E[ F t ]) f(τ,x(τ))dτ = (id E[ F t ]) t t+ t L (Ω;R d ) f(τ,x(τ)) f(t,x(t))dτ. L (Ω;R d ) By making use of (id E[ F t ])Y L (Ω;R d ) Y L (Ω;R d ) one directly deduces the desired estimate from Lemma 5.5. Finally, the last term is estimated by Lemma C-stability of the split-step backward Milstein method In this section we verify that Assumption.1 and condition (6) are sufficient for the C-stability of the split-step backward Milstein method. The results of Proposition 6.1 below are needed in order to show that the SSBM method is a well-defined one-step method in the sense of Definition 3.1. Further, the inequality (50) plays a key role in the proof of the C-stability of the SSBM method and generalizes a similar estimate for the split-step backward Euler method from [, Corollary 4.]. Proposition 6.1. Let the functions f: [0,T] R d R d and g r : [0,T] R d R d, r = 1,...,m, satisfy Assumption.1 and condition (6) with L (0, ), η 1 (1, ), and η (0, ). Let h (0,L 1 ) be given and define for every (0,h] the mapping F : [0,T] R d R d by F (t,x) = x f(t,x). Then, the mapping R d x F (t,x) R d is a homeomorphism for every t [0,T]. In addition, the inverse F 1 (t, ): R d R d satisfies F 1 (48) (t,x 1 ) F 1 (t,x ) (1 L) 1 x 1 x, F 1 (t,x) (1 L) 1 (49) ( L + x ),

19 STOCHASTIC C-STAB. AND B-CONS. OF MILSTEIN-TYPE SCHEMES 19 for every x,x 1,x R d and t [0,T]. Moreover, there exists a constant C 1 only depending on L and h such that (50) F 1 (t,x 1 ) F 1 (t,x ) +η 1 +η r 1,r =1 g r (t,f 1 (t,x 1 )) g r (t,f 1 (t,x )) g r 1,r (t,f 1 (t,x 1 )) g r1,r (t,f 1 (t,x )) (1+C 1 ) x1 x for every x 1,x R d and t [0,T]. Proof. The first part is a direct consequence of the Uniform Monotonicity Theorem (see for instance, [19, Chap.6.4], [, Theorem C.]). The estimates (48) and (49) are standard and a proof is found, for example, in [, Sec. 4]. Regarding (50) it first follows from (6) that F (t,x 1 ) F (t,x ),x 1 x = x 1 x f(t,x 1 ) f(t,x ),x 1 x (1 L) x 1 x +η 1 g r (t,x 1 ) g r (t,x ) +η r 1,r =1 g r 1,r (t,x 1 ) g r1,r (t,x ) for all x 1,x R d. For some y 1,y R d we substitute x 1 = F 1 (t,y 1 ) and x = F 1 (t,y ) into this inequality. Then, after some rearranging we obtain F 1 (t,y 1 ) F 1 (t,y ) +η 1 g r (t,f 1 (t,y 1 )) g r (t,f 1 (t,y )) +η r 1,r =1 g r 1,r (t,f 1 (t,y 1 )) g r1,r (t,f 1 (t,y )) y 1 y,f 1 (t,y 1 ) F 1 (t,y ) +L F 1 (t,y 1 ) F 1 (t,y ). Next, as in the proof of [, Corollary 4.] we apply the Cauchy-Schwarz inequality and (48). This yields y1 y,f 1 (t,y 1 ) F 1 (t,y ) +L F 1 (t,y 1 ) F 1 (t,y ) y 1 y F 1 (t,y 1 ) F 1 (t,y ) +L F 1 (t,y 1 ) F 1 (t,y ) (1 L) 1( 1+(1 L) 1 L ) y 1 y = (1 L) y 1 y for all y 1,y R d. Finally, note that b() = (1 L) is a convex function, hence for all [0,h], (1 L) 1+C 1, with C 1 = b(h) b(0) h and inequality (50) is verified. = L( Lh)(1 Lh), Proposition 6.1 ensures that the implicit step of the SSBM method admits a unique solution if f satisfies Assumption.1 with one-sided Lipschitz constant L. To be more precise, for a given h (0,L 1 ) let us consider an arbitrary vector

20 0 W.-J. BEYN, E. ISAAK, AND R. KRUSE of step sizes h (0,h] N, N N. Then, it follows from Proposition 6.1 that the nonlinear equations X SSBM h (t i ) = Xh SSBM (t i 1 )+h i f(t i,x SSBM h (t i )), 1 i N, are uniquely solvable. Further, there exists a homeomorphism F hi (t i, ): R d R d such that X SSBM h (t i ) = F 1 h i (t i,xh SSBM (t i 1 )). Therefore, the one-step map Ψ SSBM : R d T Ω R d of the split-step backward Milstein method is given by Ψ SSBM (x,t,) = F 1 (t+,x)+ g r (t+,f 1 (t+,x))i t,t+ (r) (51) + r 1,r =1 g r1,r (t+,f 1 (t+,x))i t,t+ (r,r 1) for every x R d and (t,) T, where the stochastic increments are defined in () and (3). Next, we verify that Ψ SSBM satisfies condition (8) as well as (31) and (3). Proposition 6.. Let the functions f and g r, r = 1,...,m, satisfy Assumption.1 and condition (6) with L (0, ), q [, ), η 1 (1, ), and η (0, ). For every h (0,max(L 1, η η 1 )) and initial value ξ L (Ω;F 0,P;R d ) it holds that (Ψ SSBM,h,ξ) is a stochastic one-step method. In addition, there exists a constant C 0 depending on L, q, m, and h, such that (5) (53) for all (t,) T. E [ Ψ SSBM (0,t,) F t ] L (Ω;R d ) C 0, ( id E[ F t ] ) Ψ SSBM (0,t,) L (Ω;R d ) C 0 1 Proof. Regarding the first assertion we show that Ψ SSBM satisfies (8). For this we fix arbitrary (t,) T and Z L (Ω,F t,p;r d ). Then, we obtain from Proposition 6.1 that the mapping F 1 (t+, ): R d R d is a homeomorphism satisfying the linear growth bound (49). Hence, we have F 1 (t+,z) L (Ω,F t,p;r d ). Consequently, by the continuity of g r and g r1,r the mappings and Ω ω g r (t+,f 1 (t+,z(ω))) R d Ω ω g r1,r (t+,f 1 (t+,z(ω))) R d aref t /B(R d )-measurable for everyr,r 1,r = 1,...,m. Hence,Ψ SSBM (Z,t,): Ω R d is an F t+ /B(R d )-measurable random variable. Next, we show that Ψ SSBM (Z,t,) is square integrable. First, it follows from (49) that E [ Ψ SSBM (Z,t,) F t ] L (Ω;R d ) = F 1 (t+,z) L (Ω;R d ) (1 L) 1( L + Z L (Ω;R d )). In particular, if Z = 0 L (Ω;R d ) we get E [ Ψ SSBM (0,t,) F t ] L (Ω;R d ) (1 Lh) 1 L,

21 STOCHASTIC C-STAB. AND B-CONS. OF MILSTEIN-TYPE SCHEMES 1 which is (5). Further, since the stochastic incrementsi t,t+ (r) andi t,t+ (r 1,r ) are pairwise uncorrelated and satisfy E[ I t,t+ (r) ] = and E[ I t,t+ (r 1,r ) ] = 1 we obtain ( id E[ F t ] ) Ψ SSBM (0,t,) L (Ω;R d ) = g r (t+,f 1 (t+,0))i t,t+ + = r 1,r =1 (r) L (Ω;R d ) g r1,r (t+,f 1 (t+,0))i t,t+ g r (t+,f 1 (t+,0)) + 1 Then, applications of (1) and (49) yield (r,r 1) L (Ω;R d ) r 1,r =1 g r1,r (t+,f 1 (t+,0)). g r (t+,f 1 (t+,0)) L ( 1+ F 1 (t+,0) ) q+1 L ( 1+(1 Lh) 1 Lh ) q+1 and, similarly, by (17) g r 1,r (t+,f 1 (t+,0)) L ( 1+ F 1 (t+,0) ) q L ( 1+(1 Lh) 1 Lh ) q. Therefore, there exists a constant C 0 depending on L, q, m, and h, such that (53) is satisfied. In particular, this proves that Ψ SSBM (0,t,) L (Ω,F t+,p;r d ). Next, for arbitrary Z L (Ω;F t,p;r d ) the same arguments as above yield Ψ SSBM (Z,t,) Ψ SSBM (0,t,) L (Ω;R d ) = F 1 (t+,z) F 1 (t+,0) L (Ω;R d ) + g r (t+,f 1 (t+,z)) g r (t+,f 1 (t+,0)) L (Ω;R d ) + r 1,r =1 g r 1,r (t+,f 1 (t+,z)) g r1,r (t+,f 1 (t+,0)) L (Ω;R d ). Note that η 1 > 1 and h η. Thus, the inequality (50) is applicable and we obtain Ψ SSBM (Z,t,) Ψ SSBM (0,t,) L (Ω;R d ) (1+C 1) Z L (Ω;R d ). Hence Ψ SSBM (Z,t,) L (Ω,F t+,p;r d ). Theorem 6.3. Let the functions f and g r, r = 1,...,m, satisfy Assumption.1 and condition (6) with L (0, ), η 1 (1, ), and η (0, ). Further, let h (0,max(L 1, η η 1 )). Then, for every ξ L (Ω,F 0,P;R d ) the SSBM scheme (Ψ SSBM,h,ξ) is stochastically C-stable. Proof. Let (t,) T be arbitrary. For every Y L (Ω,F t,p;r d ) we have E [ Ψ SSBM (Y,t,) F t ] = F 1 (t+,y)

22 W.-J. BEYN, E. ISAAK, AND R. KRUSE and ( id E[ Ft ] ) Ψ SSBM (Y,t,) = + g r (t+,f 1 (t+,y))i t,t+ (r) r 1,r =1 g r1,r (t+,f 1 (t+,y))i t,t+ (r,r 1). For the computation of the L -norm, we make use of the facts that the stochastic increments are independent of F t and pairwise uncorrelated. Further, since E[ I t,t+ (r) ] = and E[ I t,t+ (r 1,r ) ] = 1 it follows for (30) with ν = η 1 [ ] E Ψ SSBM (Y,t,) Ψ SSBM (Z,t,) F t L (Ω;R d ) +ν ( id E[ F t ] )( Ψ SSBM (Y,t,) Ψ SSBM (Z,t,) ) L (Ω;R d ) = E [ F 1 (t+,y) F 1 (t+,z) ] +η 1 E [ g r (t+,f 1 (t+,y)) g r (t+,f 1 (t+,z)) ] + 1 η 1 m r 1,r =1 E [ g r1,r (t+,f 1 (t+,y)) g r1,r (t+,f 1 (t+,z)) ]. Due to 1 η 1 1 η 1h < η an application of inequality (50) yields E [ Ψ SSBM (Y,t,) Ψ SSBM (Z,t,) F t ] L (Ω;R d ) +ν ( id E[ F t ] )( Ψ SSBM (Y,t,) Ψ SSBM (Z,t,) ) L (Ω;R d ) (1+C 1 ) Y Z L (Ω;R d ), which is the C-stability condition (30) with ν = η 1 (1, ). 7. B-consistency of the split-step backward Milstein method This section is devoted to the proof of the following result, which is concerned with the B-consistency of the SSBM method. Theorem 7.1. Let the functions f and g r, r = 1,...,m, satisfy Assumption.1 with L (0, ) and q [, ). Let h (0,L 1 ). If the exact solution X to (3) satisfies sup τ [0,T] X(τ) L 6q 4 (Ω;R d ) <, then the split-step backward Milstein method (Ψ SSBM,h,X 0 ) is stochastically B-consistent of order γ = 1. For the proof we recall some estimates of the homeomorphism F 1 from [, Lemma 4.3], which will be useful for the estimate of the local truncation error. Lemma 7.. Consider the same situation as in Proposition 6.1. Then there exist constants C, C 3 only depending on L, h and q such that for every (0,h] the inverse F 1 (t, ): R d R d satisfies the estimates (54) (55) for every x R d and t [0,T]. F 1 (t,x) x C ( 1+ x q ), F 1 (t,x) x f(t,x) C 3 ( 1+ x q 1 )

23 STOCHASTIC C-STAB. AND B-CONS. OF MILSTEIN-TYPE SCHEMES 3 Proof of Theorem 7.1. The proof follows the same steps as the proof of [, Theorem 5.7]. Let us fix arbitrary (t,) T. Then, by inserting (19) and (51) we obtain the following representation of the local truncation error X(t+) Ψ SSBM (X(t),t,) = X(t)+f(t+,X(t)) F 1 (t+,x(t)) + t+ t ( m r 1,r =1 ( f(τ,x(τ)) f(t+,x(t)) ) dτ t+ t r 1,r =1 ( g r (τ,x(τ)) g r (t+,x(t)) ) dw r (τ) ) g r1,r (t+,x(t))i t,t+ (r,r 1) ( g r (t+,x(t)) g r (t+,f 1 (t+,x(t))) ) I t,t+ (r) ( g r 1,r (t+,x(t)) g r1,r (t+,f 1 (t+,x(t))) ) I t,t+ (r,r 1) =: T 1 +T +T 3 +T 4 +T 5. We discuss the five terms separately. It is already shown in the proof of [, Theorem 5.7] that by applying (55) the L -norm of the term T 1 is dominated by (56) T 1 L (Ω;R d ) = X(t)+f(t+,X(t)) F 1 (t+,x(t)) L (Ω;R d ) C 31+ X(t) q 1 L (Ω;R) ) C 3 (1+ sup X(τ) q 1. L 4q (Ω;R d ) τ [0,T] Moreover, if we consider the conditional expectation of the term T with respect to F t, then after taking the L -norm we arrive at [ ] [ t+ E T F L E t (Ω;R d ) = t ( f(τ,x(τ)) f(t+,x(t)) ) dτ Ft ] L (Ω;R d). Hence, an application of Lemma 5.6 with t 1 = t and s = t+ yields [ ] E T F L t (Ω;R d ) C( 1+ sup X(t) 3q (57) L 6q 4 (Ω;R )). d Since E[T 1 F t ] = T 1 and E[T i F t ] = 0 for i {3,4,5} we get from (56) and (57) E [ X(t+) Ψ SSBM (X(t),t,) Ft ] L (Ω;R d ) T 1 L (Ω;R d ) + [ ] E T F L t (Ω;R d ) C ( 1+ sup X(t) 3q L 6q 4 (Ω;R )). d This proves (33) with γ = 1 and it remains to show (34). For this we estimate ( id E[ F t ] )( X(t+) Ψ SSBM (X(t),t,) ) L (Ω;R d ) 5 ( id E[ F t ] ) T L i (Ω;R d ). i=1 Then, note that (id E[ F t ])T 1 L (Ω;R d ) = 0 since T 1 is F t -measurable. Further, by making use of the fact that (id E[ F t ])Y L (Ω;R d ) Y L (Ω;R d ) for all

24 4 W.-J. BEYN, E. ISAAK, AND R. KRUSE Y L (Ω;R d ) we get ( id E[ F t ] ) T L (Ω;R d ) T L (Ω;R d ). After inserting T it follows from Lemma 5.5 that T L (Ω;R d ) t+ t f(τ,x(τ)) f(t+,x(t)) L (Ω;R d ) dτ C ( 1+ sup X(t) q 1 L 4q (Ω;R d ) Regarding the term T 3 we first couple the summation indices r = r 1. Then, the triangle inequality yields t+ ( T 3 L (Ω;R d ) g r (τ,x(τ)) g r (t+,x(t)) ) dw r (τ) t ) 3 g r,r (t+,x(t))i t,t+. L (Ω;R d ) r =1 (r,r) Hence, we are in the situation of Lemma 5.7 with t 1 = t and t = s = t+ and we obtain T 3 L (Ω;R d ) C ( 1+ sup X(t) 3q ) 3 L 6q 4 (Ω;R d ). The L -norm estimates of the remaining terms T 4 and T 5 follow the same line of arguments as the last part of the proof of [, Theorem 5.7]. For instance, the term T 5 is estimated as follows: From (18), (49), and (54) we obtain g r1,r (t+,x(t)) g r1,r (t+,f 1 (t+,x(t))) L ( 1+ X(t) + F 1 (t+,x(t)) ) q 1 X(t) F 1 (t+,x(t)) C L ( 1+ X(t) +(1 L) 1 (L + X(t) ) ) q 1( 1+ X(t) q ) C ( 1+ X(t) q 1), for a constant C only depending on C, L, q, and h. Therefore, ( g r 1,r (t+,x(t)) g r1,r (t+,f 1 (t+,x(t))) ) I t,t+ r 1,r =1 = 1 m r 1,r =1. (r,r 1) L (Ω;R d ) g r 1,r (t+,x(t)) g r1,r (t+,f 1 (t+,x(t))) L (Ω;R d ) Cm ( ) 1+ sup X(τ) q 1 4, L 4q (Ω;R d ) τ [0,T] which is the desired estimate of T 5. The corresponding estimate of T 4 reads ( ) T 4 L (Ω;R d ) Cm 1+ sup X(τ) 3 q 1 3 L 3q 1 (Ω;R d ) τ [0,T] and is obtained in the same way as for T 5 but with (15) in place of (18). Altogether, this completes the proof of (34) with γ = Numerical experiments In this section we perform several numerical experiments which illustrate the preceding theory for two characteristic examples.