Computing ergodic limits for SDEs
|
|
- Peter Shields
- 6 years ago
- Views:
Transcription
1 Computing ergodic limits for SDEs M.V. Tretyakov School of Mathematical Sciences, University of Nottingham, UK Talk at the workshop Stochastic numerical algorithms, multiscale modelling and high-dimensional data analytics, ICERM, 21st July 2016
2 Plan of the talk ϕ erg = 1 t ϕ(x)ρ(x) dx = lim Eϕ(X (t)) = lim ϕ(x (s))ds a.s. t t t 0
3 Plan of the talk ϕ erg = 1 t ϕ(x)ρ(x) dx = lim Eϕ(X (t)) = lim ϕ(x (s))ds a.s. t t t 0 Introduction Examples of Langevin-type equations and stochastic gradient systems Geometric integrators for Langevin equations and the gradient system Non-Markovian scheme for stochastic gradient systems [Davidchack, Ouldridge&T. J Chem Phys 2015] [Leimkuhler, Matthews,T 2014]
4 Introduction Hamiltonian H(r, p) canonical ensemble (NVT ) ρ(r, p) exp( βh(r, p)), where β = 1/(k B T ) > 0 is an inverse temperature. ϕ erg = 1 t ϕ(r)ρ(r, p) drdp = lim Eϕ(R(t)) = lim ϕ(r(s))ds a.s. t t t 0
5 Rigid Body Dynamics Following [Miller III et al J. Chem. Phys., 2002] H(r, p, q, π) = pt p 2m + n j=1 l=1 3 V l (q j, π j ) + U(r, q), (1) where r = (r 1T,..., r nt ) T R 3n, and p = (p 1T,..., p nt ) T R 3n are the center-of-mass coordinates and momenta; q = (q 1T,..., q nt ) T R 4n, q j = (q j 0, qj 1, qj 2, qj 3 )T R 4, q j = 1, are the rotational coordinates in the quaternion representation, π = (π 1T,..., π nt ) T R 4n, are the angular momenta; 3 V l (q, π) = 1 8 l=1 3 l=1 1 I l [ π T S l q ] 2 = 1 8 πt S(q)DS T (q)π, q, π R 4, l = 1, 2, 3, I l the principal moments of inertia and the constant 4-by-4 matrices S l : S 1 q = ( q 1, q 0, q 3, q 2 ) T, S 2 q = ( q 2, q 3, q 0, q 1 ) T, S 3 q = ( q 3, q 2, q 1, q 0 ) T, S 0 = diag(1, 1, 1, 1) T. S(q) = [S 0 q, S 1 q, S 2 q, S 3 q], D = diag(0, 1/I 1, 1/I 2, 1/I 3 ).
6 Langevin thermostat for Rigid Body Dynamics dr j = Pj m dt, Rj (0) = r j, dp j = f j (R, Q)dt γp j dt + 2mγ β dw j (t), P j (0) = p j, (2) dq j = 1 4 S(Qj )DS T (Q j )Π j dt, Q j (0) = q j, q j = 1, (3) 3 dπ j = l=1 1 I l ( Π j T S l Q j) S l Π j dt + F j (R, Q)dt ΓJ(Q j )Π j dt 2MΓ β 3 l=1 S l Q j dw j l (t), Πj (0) = π j, q j T π j = 0, j = 1,..., n, where f j (r, q) = r j U(r, q) R 3, F j (r, q) = q j U(r, q) T q j S 3, (w T, W T ) T = (w 1 T,..., w n T, W 1 T,..., W n T ) T is a (3n + 3n)-dimensional standard Wiener process; γ 0 and Γ 0 are the friction coefficients for the translational and rotational motions, and 3 J(q) = MS(q)DS T (q)/4, M = 4/ 1/I l. Davidchack, Ouldridge&T. J Chem Phys 2015 l=1
7 Langevin-type equations The Ito interpretation of the SDEs (2) (3) coincides with its Stratonovich interpretation. The solution of (2) (3) preserves the quaternion length Q j (t) = 1, j = 1,..., n, for all t 0. (4) The solution of (2) (3) automatically preserves the constraint: Q j T (t)π j (t) = 0, j = 1,..., n, for t 0 (5) Assume that the solution X (t) = (R T (t), P T (t), Q T (t), Π T (t)) T of (2) (3) is an ergodic process on D = {x = (r T, p T, q T, π T ) T R 14n : q j = 1, q j T π j = 0, j = 1,..., n}. Then it can be shown that the invariant measure of X (t) is Gibbsian with the density ρ(r, p, q, π) on D: ρ(r, p, q, π) exp( βh(r, p, q, π)) (6) Davidchack, Ouldridge&T. J Chem Phys 2015
8 Langevin equations and quasi-symplectic integrators dr j = Pj m dt, Rj (0) = r j, dp j = f j (R, Q)dt γp j dt + 2mγ β dw j (t), P j (0) = p j, (9) dq j = 1 4 S(Qj )DS T (Q j )Π j dt, Q j (0) = q j, q j = 1, (10) 3 dπ j = l=1 1 I l ( Π j T S l Q j) S l Π j dt + F j (R, Q)dt ΓJ(Q j )Π j dt 2MΓ β 3 l=1 S l Q j dw j l (t), Πj (0) = π j, q j T π j = 0, j = 1,..., n,
9 Langevin equations and quasi-symplectic integrators dr j = Pj m dt, Rj (0) = r j, dp j = f j (R, Q)dt γp j dt + 2mγ β dw j (t), P j (0) = p j, (9) dq j = 1 4 S(Qj )DS T (Q j )Π j dt, Q j (0) = q j, q j = 1, (10) 3 dπ j = l=1 1 I l ( Π j T S l Q j) S l Π j dt + F j (R, Q)dt ΓJ(Q j )Π j dt 2MΓ β 3 l=1 S l Q j dw j l (t), Πj (0) = π j, q j T π j = 0, j = 1,..., n, Let D 0 R d, d = 14n, be a domain with finite volume. The transformation x = (r, p, q, π) X (t) = X (t; x) = (R(t; x), P(t; x), Q(t; x), Π(t; x)) maps D 0 into the domain D t.
10 Langevin equations and quasi-symplectic integrators V t = = dx 1... dx d (7) D t D(X 1,..., X d ) D(x 1,..., x d ) dx 1... dx d. D 0 The Jacobian J is equal to J = D(X 1,..., X d ) D(x 1,..., x d ) = exp ( n(3γ + Γ) t). (8) [Milstein&T. IMA J. Numer. Anal (also Springer 2004)], [Davidchack, Ouldridge&T. J Chem Phys 2015]
11 The stochastic gradient system It is easy to verify that exp( βh(r, p, q, π))dpdπ (9) D mom exp( βu(r, q)) =: ρ(r, q), where (r T, q T ) T D = {(r T, q T ) T R 7n : q j = 1} and the domain of conjugate momenta D mom = {(p T, π T ) T R 7n : q T π = 0}.
12 The stochastic gradient system It is easy to verify that exp( βh(r, p, q, π))dpdπ (9) D mom exp( βu(r, q)) =: ρ(r, q), where (r T, q T ) T D = {(r T, q T ) T R 7n : q j = 1} and the domain of conjugate momenta D mom = {(p T, π T ) T R 7n : q T π = 0}. We introduce the gradient system in the form of Stratonovich SDEs: dr = υ m f(r, Q)dt + 2υ dw(t), R(0) = r, (10) mβ dq j = Υ M F j 2Υ 3 (R, Q)dt + S l Q j dw j l (t), (11) Mβ l=1 Q j (0) = q j, q j = 1, j = 1,..., n, where indicates the Stratonovich form of the SDEs, parameters υ > 0 and Υ > 0 control the speed of evolution of the gradient system (10) (11), f = (f 1 T,..., f n T ) T and the rest of the notation is as in (2) (3). [Davidchack, Ouldridge&T. J Chem Phys 2015]
13 The gradient thermostat for rigid body dynamics This new gradient thermostat possesses the following properties. As in the case of (2) (3), the solution of (10) (11) preserves the quaternion length (4). Assume that the solution X (t) = (R T (t), Q T (t)) T D of (10) (11) is an ergodic process. Then, by the usual means of the stationary Fokker-Planck equation, one can show that its invariant measure is Gibbsian with the density ρ(r, q) from (9).
14 Langevin integrators Davidchack, Ouldridge&T. J Chem Phys 2015 For simplicity we use a uniform time discretization of a time interval [0, T ] with the step h = T /N. Goal: to construct integrators quasi-symplectic preserve Q j (t k ) = 1, j = 1,..., n, for all t 0 automatically preserve Q j T (t k ) Π j (t k ) = 0, of weak order 2 To this end: j = 1,..., n, for t 0 automatically stochastic numerics+splitting techniques [see e.g. Milstein&T, Springer 2004] the deterministic symplectic integrator from [Miller III et al J. Chem. Phys., 2002]
15 Langevin A integrator The first integrator is based on splitting the Langevin system in dr j = Pj m dt, Rj (0) = r j, (12) dp j = f j 2mγ (R, Q)dt + β dw j (t), dq j = 1 4 S(Qj )DS T (Q j )Π j dt, (13) 3 dπ j = l=1 1 I l ( Π j T S l Q j) S l Π j dt + F j (R, Q)dt 2MΓ β 3 l=1 S l Q j dw j l (t), j = 1,..., n,
16 Langevin A integrator The first integrator is based on splitting the Langevin system in dr j = Pj m dt, Rj (0) = r j, (12) dp j = f j 2mγ (R, Q)dt + β dw j (t), dq j = 1 4 S(Qj )DS T (Q j )Π j dt, (13) 3 dπ j = l=1 1 I l ( Π j T S l Q j) S l Π j dt + F j (R, Q)dt 2MΓ β 3 l=1 S l Q j dw j l (t), j = 1,..., n, and the deterministic system of linear differential equations ṗ = γp, π j = ΓJ(q j )π j, j = 1,..., n. (14)
17 Langevin B integrator 2mγ dp I = γp I dt + β dw(t), dπ j I = ΓJ(q)Π j I dt + 2MΓ 3 S l qdw j l β (t); l=1 (15) dr II = P II m dt, dp II = f(r II, Q II )dt, dq j II = 1 4 S(Qj II )DS T (Q j II )Πj II dt, (16) dπ j II = F j (R II, Q II )dt ] [(Π j II 4 I )T S l Q j II S l Π j II dt, j = 1,..., n. l l=1
18 Langevin B integrator 2mγ dp I = γp I dt + β dw(t), dπ j I = ΓJ(q)Π j I dt + 2MΓ 3 S l qdw j l β (t); l=1 (15) dr II = P II m dt, dp II = f(r II, Q II )dt, dq j II = 1 4 S(Qj II )DS T (Q j II )Πj II dt, (16) dπ j II = F j (R II, Q II )dt ] [(Π j II 4 I )T S l Q j II S l Π j II dt, j = 1,..., n. l l=1 The SDEs (15) have the exact solution: 2mγ t P I (t) = P I (0) exp( γt) + β 0 Π j I (t) = exp( ΓJ(q)t)Πj I (0) + 2MΓ β exp( γ(t s))dw(s), (17) 3 l=1 t 0 exp( ΓJ(q)(t s))dw j l (s). To construct the method: half a step of (17) + one step of a symplectic method for (16) + half a step of (17).
19 Langevin C integrator Based on the same spliting (15) and (16) as Langevin B, i.e., the determinisitic Hamiltonian system + OU. To construct the method: half a step of a symplectic method for (16) + step of (17) + half a step of a symplectic method for (16).
20 Langevin C integrator Based on the same spliting (15) and (16) as Langevin B, i.e., the determinisitic Hamiltonian system + OU. To construct the method: half a step of a symplectic method for (16) + step of (17) + half a step of a symplectic method for (16). Various splittings are compared for a translational Langevin thermostat in [Leimkuhler&Matthews 2013]
21 Numerical experiments Davidchack, Handel&T. J Chem Phys 2009 and Davidchack, Ouldridge&T. J Chem Phys 2015 Two objectives for the experiments: the dependence of Langevin dynamics properties on the choice of parameters γ and Γ errors of the numerical schemes. TIP4P rigid model of water (Jorgensen et. al J. Chem. Phys. 1983) A h : = EA( X ) = EA(X ) + E A h p + O(h p+1 ) = A 0 + E A h p + O(h p+1 ) p = 2 for Langevin integrators and p = 1 for the gradient thermostat integrator Talay&Tubaro Stoch.Anal.Appl. 1990
22 Accuracy of integrators Langevin A (left), Langevin B (centre), and Langevin C (right) with γ = 5 ps 1 and Γ = 10 ps 1. Error bars denote estimated 95% confidence intervals.
23 Geometric integrator for the gradient thermostat The main idea is to rewrite the components Q j of the solution to (10) (11) in the form Q j (t) = exp(y j (t))q j (0) and then solve numerically the SDEs for the 4 4-matrices Y j (t)
24 Geometric integrator for the gradient thermostat The main idea is to rewrite the components Q j of the solution to (10) (11) in the form Q j (t) = exp(y j (t))q j (0) and then solve numerically the SDEs for the 4 4-matrices Y j (t) Lie group methods [Hairer, Lubich, Wanner; Springer, 2002] To this end, we introduce the 4 4 skew-symmetric matrices: F j (r, q) = F j (r, q)q j T q j (F j (r, q)) T, j = 1,..., n. Note that F j (r, q)q j = F j (r, q) under q j = 1 and the equations (11) can be written as dq j = Υ M F j(r, Q)Q j 2Υ 3 dt+ S l Q j dw j l Mβ (t), Qj (0) = q j, q j = 1. One can show that l=1 Y j (t + h) = h Υ M F j(r(t), Q(t)) + 2Υ Mβ 3 ( l=1 W j l + terms of higher order. (18) (t + h) W j l (t) ) S l
25 Geometric integrator for the gradient thermostat R 0 = r, Q 0 = q, q j = 1, j = 1,..., n, (19) R k+1 = R k + h υ m f(r k, Q k ) + 2υ h mβ ξ k, Y j k = h Υ M F j(r k, Q k ) + 2Υ 3 h η j,l k Mβ S l, Q j k+1 = exp(y j k )Qj k, j = 1,..., n, where ξ k = (ξ 1,k,..., ξ 3n,k ) T and ξ i,k, i = 1,..., 3n, η j,l k, l = 1, 2, 3, j = 1,..., n, are i.i.d. random variables with the same law l=1 P(θ = ±1) = 1/2. (20) Proposition. The numerical scheme (19) (20) for (10) (11) preserves the length of quaternions, i.e., Q j k = 1, j = 1,..., n, for all k, and it is of weak order one. Davidchack, Ouldridge&T. J Chem Phys 2015
26 Stochastic gradient system dx = a(x)dt + σdw, X(0) = X 0, (21) a(x) := U(x), x R d, (22) σ = 2/β, d = 3n, and w(t) is a standard d-dimensional Wiener process.
27 Stochastic gradient system dx = a(x)dt + σdw, X(0) = X 0, (21) a(x) := U(x), x R d, (22) σ = 2/β, d = 3n, and w(t) is a standard d-dimensional Wiener process. We use the following notation for the solution of (21): X(t) = X t0,x(t) when X(t 0 ) = x, t t 0, and also we will write X x (t) when t 0 = 0.
28 Stochastic gradient system dx = a(x)dt + σdw, X(0) = X 0, (21) a(x) := U(x), x R d, (22) σ = 2/β, d = 3n, and w(t) is a standard d-dimensional Wiener process. We use the following notation for the solution of (21): X(t) = X t0,x(t) when X(t 0 ) = x, t t 0, and also we will write X x (t) when t 0 = 0. The process X(t) is exponentially ergodic [Hasminskii 1980] if for any x R d and any function ϕ with a polynomial growth there are C(x) > 0 and λ > 0 such that Eϕ(X x (t)) ϕ erg C(x)e λt, t 0,
29 Stochastic gradient system dx = a(x)dt + σdw, X(0) = X 0, (21) a(x) := U(x), x R d, (22) σ = 2/β, d = 3n, and w(t) is a standard d-dimensional Wiener process. We use the following notation for the solution of (21): X(t) = X t0,x(t) when X(t 0 ) = x, t t 0, and also we will write X x (t) when t 0 = 0. The process X(t) is exponentially ergodic [Hasminskii 1980] if for any x R d and any function ϕ with a polynomial growth there are C(x) > 0 and λ > 0 such that Eϕ(X x (t)) ϕ erg C(x)e λt, t 0, The solution X(t) of (21) is exponentially ergodic and ρ(x) exp ( 2σ ) 2 U(x) if there exist c 0 R and c 1 > 0 such that (x, a(x)) c 0 c 1 x 2. (23) [Hasminskii (1980), Mattingly, Stuart, Higham (2002)]
30 Stochastic gradient system The Euler scheme: X k+1 = X k + ha(x k ) + σ hξ k+1, (24) where ξ k = (ξ 1 k,..., ξ d k ) and ξ i k, i = 1,..., d, k = 1,..., are i.i.d. random variables with the law N (0, 1).
31 Stochastic gradient system The Euler scheme: X k+1 = X k + ha(x k ) + σ hξ k+1, (24) where ξ k = (ξ 1 k,..., ξ d k ) and ξ i k, i = 1,..., d, k = 1,..., are i.i.d. random variables with the law N (0, 1). Heun s scheme: ˆX k+1 = X k + ha(x k ) + σ hξ k+1, X k+1 = X k + h [ a( 2 ˆX ] k+1 ) + a(x k ) + σ hξ k+1. (25)
32 Stochastic gradient system The Euler scheme: X k+1 = X k + ha(x k ) + σ hξ k+1, (24) where ξ k = (ξ 1 k,..., ξ d k ) and ξ i k, i = 1,..., d, k = 1,..., are i.i.d. random variables with the law N (0, 1). Heun s scheme: ˆX k+1 = X k + ha(x k ) + σ hξ k+1, X k+1 = X k + h [ a( 2 ˆX ] k+1 ) + a(x k ) + σ hξ k+1. (25) Weak convergence: Eϕ(X (T )) Eϕ(X N ) Kh p (26) N = T /h; Euler p = 1, Heun p = 2 [e.g. Milstein, T., Springer 2004]
33 Stochastic gradient system The Euler scheme: X k+1 = X k + ha(x k ) + σ hξ k+1, (24) where ξ k = (ξ 1 k,..., ξ d k ) and ξ i k, i = 1,..., d, k = 1,..., are i.i.d. random variables with the law N (0, 1). Heun s scheme: ˆX k+1 = X k + ha(x k ) + σ hξ k+1, X k+1 = X k + h [ a( 2 ˆX ] k+1 ) + a(x k ) + σ hξ k+1. (25) Weak convergence: Eϕ(X (T )) Eϕ(X N ) Kh p (26) N = T /h; Euler p = 1, Heun p = 2 [e.g. Milstein, T., Springer 2004] Ergodic limits [Talay (1990); Talay, Tubaro (1990); Mattingly, Stuart, Higham (2002); Milstein, T (2007); Mattingly, Stuart, T. (2010)]: ϕ erg Eϕ(X N ) Kh p + Ce λt (27)
34 Non-Markovian scheme X k+1 = X k + ha(x k ) + σ h 2 (ξ k + ξ k+1 ), (28) where ξ k = (ξ 1 k,..., ξ i k) defined on (Ω, P, F) and ξ i k, i = 1,..., d, k = 1,..., are i.i.d. random variables with the law N (0, 1) [Leimkuhler, Matthews (2013)]
35 Example Let a(x) = αx with α > 0, then X(t) from (21) is the Ornstein-Uhlenbeck process, which is Gaussian with EX x (t) = xe αt, Cov(X x (s), X x (t)) = σ2 2α (e α(t s) e α(t+s) ) for s t and Var(X x (T )) = σ2 2α (1 e 2αT ).
36 Example Let a(x) = αx with α > 0, then X(t) from (21) is the Ornstein-Uhlenbeck process, which is Gaussian with EX x (t) = xe αt, Cov(X x (s), X x (t)) = σ2 2α (e α(t s) e α(t+s) ) for s t and Var(X x (T )) = σ2 2α (1 e 2αT ). For the Euler scheme (24): EX N = x 0 (1 αh) N = x 0 e αt (1 + O(h)), Var(X N ) = σ2 1 (1 αh) 2N 2α 1 + αh = σ2 2α (1 e 2αT ) σ2 2 h + e 2αT O(h) + O(h 2 ), αh < 1, where O(h p ) Kh with K > 0 independent of T.
37 Example Let a(x) = αx with α > 0, then X(t) from (21) is the Ornstein-Uhlenbeck process, which is Gaussian with EX x (t) = xe αt, Cov(X x (s), X x (t)) = σ2 2α (e α(t s) e α(t+s) ) for s t and Var(X x (T )) = σ2 2α (1 e 2αT ). For the Euler scheme (24): EX N = x 0 (1 αh) N = x 0 e αt (1 + O(h)), Var(X N ) = σ2 1 (1 αh) 2N 2α 1 + αh = σ2 2α (1 e 2αT ) σ2 2 h + e 2αT O(h) + O(h 2 ), αh < 1, where O(h p ) Kh with K > 0 independent of T. For the scheme (28): EX N = x 0 (1 αh) N = x 0 e αt (1 + O(h)), ] Var(X N ) = [1 σ2 (1 αh)2n = σ2 2α 1 αh 2α (1 e 2αT ) + e 2αT O(h).
38 Assumptions Assumption 1 There exist c 0 R and c 1 > 0 such that (x, a(x)) c 0 c 1 x 2. Assumption 2 The potential U(x) C 7 (R d ), its first-order derivatives grow not faster than a linear function at infinity and higher derivatives are bounded. The function ϕ(x) C 6 (R d ) and it and its derivatives grow not faster than a polynomial function at infinity. The most restrictive condition in Assumption 2 is the requirement for a(x) = U to be globally Lipschitz: a(x) 2 K(1 + x 2 ), (29) where K > 0 is independent of x R d, which can be relaxed via [Milstein, T. (2005) and (2007): the concept of rejecting exploding trajectories]
39 Introduce the operator L L := t + L, where L The function L := d i=1 a i (x) x i + σ2 2 d i=1 2 ( x i ) 2. (30) u(t, x) = Eϕ(X t,x (T )) (31) satisfies the Cauchy problem for the backward Kolmogorov equation Lu = 0, (32) u(t, x) = ϕ(x).
40 Main theorem Theorem (1) Let Assumptions 1-2 hold. Then the scheme (28) is first order weakly convergent and for all sufficiently small h > 0 its error has the form + σ2 2 Eϕ(X x (T )) Eϕ(X N ) = C 0 (T, x)h + C(T, x)h 2, (33) T C 0 (T, x) = E B 0 (t, X x (t))dt, (34) 0 B 0 (t, x) = 1 d 2 d i,j=1 i,j=1 a i (x) 2 u(t, x) x j x i x j + σ2 2 a j (x) ai (x) u(t, x) x j x i d 2 a i (x) u(t, x) ( x j ) 2 x i, i,j=1 C(T, x) K(1 + x κ e λt ), for some K > 0, κ N and λ > 0 independent of h and T.
41 Corollary Theorem 1: Eϕ(X x (T )) Eϕ(X N ) = C 0 (T, x)h + C(T, x)h 2, T C 0 (T, x) = E B 0 (t, X x (t))dt. 0
42 Corollary Theorem 1: Eϕ(X x (T )) Eϕ(X N ) = C 0 (T, x)h + C(T, x)h 2, T C 0 (T, x) = E B 0 (t, X x (t))dt. 0 Theorem (2) Let Assumptions 1-2 hold. Then the coefficient C 0 (T, x) from (34) goes to zero as T : C 0 (T, x) K(1 + x κ )e λt (35) for some constants K > 0, κ N and λ > 0, i.e., over a long integration time the scheme (28) is of order two up to exponentially small correction.
43 Sketch of the proof C 0 (T, x) = = T 0 T 0 T + 0 EB 0 (t, X x (t))dt = T B 0 (t, y)ρ(y)dydt R d R d B 0 (t, y)[p(t, x, y) ρ(y)]dydt, 0 R d B 0 (t, y)p(t, x, y)dydt (36) where p(t, x, y) is the transition density for (21) and ρ(y) is the invariant density.
44 Sketch of the proof C 0 (T, x) = = T 0 T 0 T + 0 EB 0 (t, X x (t))dt = T B 0 (t, y)ρ(y)dydt R d R d B 0 (t, y)[p(t, x, y) ρ(y)]dydt, 0 R d B 0 (t, y)p(t, x, y)dydt (36) where p(t, x, y) is the transition density for (21) and ρ(y) is the invariant density. B 0 (t, y) exp ( 2σ ) R 2 U(y) dy = 0. (37) d
45 Discussion 1 We emphasize that the fact that the average of B 0 (t, x) with respect to the invariant measure is equal to zero is the reason why the scheme (28) is second order accurate in approximating ergodic limits.
46 Discussion 1 We emphasize that the fact that the average of B 0 (t, x) with respect to the invariant measure is equal to zero is the reason why the scheme (28) is second order accurate in approximating ergodic limits. 2 In the case of the Euler scheme (24) we get the same error expansion as (33) for the scheme (28) but with a different B 0 (t, x) = B0 E (t, x) (see [Milstein, T. (2004)]): B0 E (t, x) = 1 d a j u u 2 x j ai x i + σ2 2 +σ 2 d i,j=1 d i,j i,j=1 2 a j u ( x i ) 2 x j + σ2 2 a j 2 u x i x j x i + σ4 6 d a i 3 u i,j=1 x i ( x j ) 2 d 4 u ( x i ) 2 ( x j ) 2. i,j=1 The average of B0 E (t, x) with respect to the invariant measure is not equal to zero and, consequently, the Euler scheme (24) approximates ergodic limits with order one the same order as its weak convergence over a finite time interval.
47 Discussion 3 Let a one-step weak approximation X t,x (t + h) of the solution X t,x (t + h) of (21) generate a method of order p. The global error of the method: R : = Eϕ(X x (T )) Eϕ( X x (T )) (38) = C 0 (T, x)h p + + C n (T, x)h p+n + O(h p+n+1 ), where n N and the functions C 0 (T, x),..., C n (T, x) are independent of h which can be presented in the form C i (T, x) = T 0 EB i (s, X x (s))ds. One can deduce from the proof of Theorem 2 that if the averages of B i (s, x) 0 i n, with respect to the invariant measure are equal to zero then in the limit of T the scheme has p + n order of accuracy in h.
48 Discussion 3 Let a one-step weak approximation X t,x (t + h) of the solution X t,x (t + h) of (21) generate a method of order p. The global error of the method: R : = Eϕ(X x (T )) Eϕ( X x (T )) (38) = C 0 (T, x)h p + + C n (T, x)h p+n + O(h p+n+1 ), where n N and the functions C 0 (T, x),..., C n (T, x) are independent of h which can be presented in the form C i (T, x) = T 0 EB i (s, X x (s))ds. One can deduce from the proof of Theorem 2 that if the averages of B i (s, x) 0 i n, with respect to the invariant measure are equal to zero then in the limit of T the scheme has p + n order of accuracy in h. [Abdulle, Vilmart, Zygalakis ]
49 Numerical experiments Anharmonic scalar model: the one-dimensional potential energy U(x) = cos(x) L 2 error: (ˆρ i ρ i ) 2 i tropy error 10-1 L 2 error uyama 10-2 Euler-Maruyama ethod didate Scheme (1.7) Error Heun s Method Candidate Scheme (1.7) Timestep Timestep Figure: The error in computed distributions is plotted for each scheme.
50 Error in finite time Error at time t t=1.92 t=3.84 t=5.76 t= Timestep Timestep Timestep Timestep 10-1 Error at time t t=0.96 t=2.88 t=4.80 t=6.72 t= Timestep Timestep Timestep Timestep Timestep 10-1 Error using h = 0.16 L error Euler-Maruyama Heun s Method Candidate Scheme (1.7) Time Figure: The lower plot shows the error in the distribution after time t, as computed using each scheme at h = In the plots at the top, we compare the error growth with respect to stepsize h at multiples of t = 0.96.
51 Conclusions Two new thermostats, one of Langevin type and one of gradient (Brownian) type, for rigid body dynamics are introduced. As in the deterministic case, it is important to preserve structural properties of stochastic systems for accurate long term simulations. Current work includes stochastic rigid body dynamics with hydrodynamic interactions. Development of more efficient methods for stochastic gradient systems.
52 Conclusions Two new thermostats, one of Langevin type and one of gradient (Brownian) type, for rigid body dynamics are introduced. As in the deterministic case, it is important to preserve structural properties of stochastic systems for accurate long term simulations. Current work includes stochastic rigid body dynamics with hydrodynamic interactions. Development of more efficient methods for stochastic gradient systems. THANK YOU!
Qualitative behaviour of numerical methods for SDEs and application to homogenization
Qualitative behaviour of numerical methods for SDEs and application to homogenization K. C. Zygalakis Oxford Centre For Collaborative Applied Mathematics, University of Oxford. Center for Nonlinear Analysis,
More informationPROBABILITY: LIMIT THEOREMS II, SPRING HOMEWORK PROBLEMS
PROBABILITY: LIMIT THEOREMS II, SPRING 218. HOMEWORK PROBLEMS PROF. YURI BAKHTIN Instructions. You are allowed to work on solutions in groups, but you are required to write up solutions on your own. Please
More informationPROBABILITY: LIMIT THEOREMS II, SPRING HOMEWORK PROBLEMS
PROBABILITY: LIMIT THEOREMS II, SPRING 15. HOMEWORK PROBLEMS PROF. YURI BAKHTIN Instructions. You are allowed to work on solutions in groups, but you are required to write up solutions on your own. Please
More informationDiscretization of SDEs: Euler Methods and Beyond
Discretization of SDEs: Euler Methods and Beyond 09-26-2006 / PRisMa 2006 Workshop Outline Introduction 1 Introduction Motivation Stochastic Differential Equations 2 The Time Discretization of SDEs Monte-Carlo
More informationWORD SERIES FOR THE ANALYSIS OF SPLITTING SDE INTEGRATORS. Alfonso Álamo/J. M. Sanz-Serna Universidad de Valladolid/Universidad Carlos III de Madrid
WORD SERIES FOR THE ANALYSIS OF SPLITTING SDE INTEGRATORS Alfonso Álamo/J. M. Sanz-Serna Universidad de Valladolid/Universidad Carlos III de Madrid 1 I. OVERVIEW 2 The importance of splitting integrators
More informationThe concentration of a drug in blood. Exponential decay. Different realizations. Exponential decay with noise. dc(t) dt.
The concentration of a drug in blood Exponential decay C12 concentration 2 4 6 8 1 C12 concentration 2 4 6 8 1 dc(t) dt = µc(t) C(t) = C()e µt 2 4 6 8 1 12 time in minutes 2 4 6 8 1 12 time in minutes
More informationExercises. T 2T. e ita φ(t)dt.
Exercises. Set #. Construct an example of a sequence of probability measures P n on R which converge weakly to a probability measure P but so that the first moments m,n = xdp n do not converge to m = xdp.
More informationA Short Introduction to Diffusion Processes and Ito Calculus
A Short Introduction to Diffusion Processes and Ito Calculus Cédric Archambeau University College, London Center for Computational Statistics and Machine Learning c.archambeau@cs.ucl.ac.uk January 24,
More informationThermostatic Controls for Noisy Gradient Systems and Applications to Machine Learning
Thermostatic Controls for Noisy Gradient Systems and Applications to Machine Learning Ben Leimkuhler University of Edinburgh Joint work with C. Matthews (Chicago), G. Stoltz (ENPC-Paris), M. Tretyakov
More informationNumerical Integration of SDEs: A Short Tutorial
Numerical Integration of SDEs: A Short Tutorial Thomas Schaffter January 19, 010 1 Introduction 1.1 Itô and Stratonovich SDEs 1-dimensional stochastic differentiable equation (SDE) is given by [6, 7] dx
More informationHomogenization with stochastic differential equations
Homogenization with stochastic differential equations Scott Hottovy shottovy@math.arizona.edu University of Arizona Program in Applied Mathematics October 12, 2011 Modeling with SDE Use SDE to model system
More informationLangevin Methods. Burkhard Dünweg Max Planck Institute for Polymer Research Ackermannweg 10 D Mainz Germany
Langevin Methods Burkhard Dünweg Max Planck Institute for Polymer Research Ackermannweg 1 D 55128 Mainz Germany Motivation Original idea: Fast and slow degrees of freedom Example: Brownian motion Replace
More informationLecture 4: Introduction to stochastic processes and stochastic calculus
Lecture 4: Introduction to stochastic processes and stochastic calculus Cédric Archambeau Centre for Computational Statistics and Machine Learning Department of Computer Science University College London
More informationAlgorithms for Ensemble Control. B Leimkuhler University of Edinburgh
Algorithms for Ensemble Control B Leimkuhler University of Edinburgh Model Reduction Across Disciplines Leicester 2014 (Stochastic) Ensemble Control extended ODE/SDE system design target ensemble SDE invariant
More informationI forgot to mention last time: in the Ito formula for two standard processes, putting
I forgot to mention last time: in the Ito formula for two standard processes, putting dx t = a t dt + b t db t dy t = α t dt + β t db t, and taking f(x, y = xy, one has f x = y, f y = x, and f xx = f yy
More informationTwo recent works on molecular systems out of equilibrium
Two recent works on molecular systems out of equilibrium Frédéric Legoll ENPC and INRIA joint work with M. Dobson, T. Lelièvre, G. Stoltz (ENPC and INRIA), A. Iacobucci and S. Olla (Dauphine). CECAM workshop:
More informationBrownian Motion. 1 Definition Brownian Motion Wiener measure... 3
Brownian Motion Contents 1 Definition 2 1.1 Brownian Motion................................. 2 1.2 Wiener measure.................................. 3 2 Construction 4 2.1 Gaussian process.................................
More informationELEMENTS OF PROBABILITY THEORY
ELEMENTS OF PROBABILITY THEORY Elements of Probability Theory A collection of subsets of a set Ω is called a σ algebra if it contains Ω and is closed under the operations of taking complements and countable
More informationStochastic Gradient Descent in Continuous Time
Stochastic Gradient Descent in Continuous Time Justin Sirignano University of Illinois at Urbana Champaign with Konstantinos Spiliopoulos (Boston University) 1 / 27 We consider a diffusion X t X = R m
More informationGSHMC: An efficient Markov chain Monte Carlo sampling method. Sebastian Reich in collaboration with Elena Akhmatskaya (Fujitsu Laboratories Europe)
GSHMC: An efficient Markov chain Monte Carlo sampling method Sebastian Reich in collaboration with Elena Akhmatskaya (Fujitsu Laboratories Europe) 1. Motivation In the first lecture, we started from a
More informationKolmogorov equations in Hilbert spaces IV
March 26, 2010 Other types of equations Let us consider the Burgers equation in = L 2 (0, 1) dx(t) = (AX(t) + b(x(t))dt + dw (t) X(0) = x, (19) where A = ξ 2, D(A) = 2 (0, 1) 0 1 (0, 1), b(x) = ξ 2 (x
More informationMathematical Methods for Neurosciences. ENS - Master MVA Paris 6 - Master Maths-Bio ( )
Mathematical Methods for Neurosciences. ENS - Master MVA Paris 6 - Master Maths-Bio (2014-2015) Etienne Tanré - Olivier Faugeras INRIA - Team Tosca November 26th, 2014 E. Tanré (INRIA - Team Tosca) Mathematical
More informationp 1 ( Y p dp) 1/p ( X p dp) 1 1 p
Doob s inequality Let X(t) be a right continuous submartingale with respect to F(t), t 1 P(sup s t X(s) λ) 1 λ {sup s t X(s) λ} X + (t)dp 2 For 1 < p
More informationTheoretical Tutorial Session 2
1 / 36 Theoretical Tutorial Session 2 Xiaoming Song Department of Mathematics Drexel University July 27, 216 Outline 2 / 36 Itô s formula Martingale representation theorem Stochastic differential equations
More informationNumerical solution of second-order stochastic differential equations with Gaussian random parameters
Journal of Linear and Topological Algebra Vol. 2, No. 4, 23, 223-235 Numerical solution of second-order stochastic differential equations with Gaussian random parameters R. Farnoosh a, H. Rezazadeh a,,
More informationHJB equations. Seminar in Stochastic Modelling in Economics and Finance January 10, 2011
Department of Probability and Mathematical Statistics Faculty of Mathematics and Physics, Charles University in Prague petrasek@karlin.mff.cuni.cz Seminar in Stochastic Modelling in Economics and Finance
More informationStochastic differential equation models in biology Susanne Ditlevsen
Stochastic differential equation models in biology Susanne Ditlevsen Introduction This chapter is concerned with continuous time processes, which are often modeled as a system of ordinary differential
More informationStochastic Differential Equations
Chapter 5 Stochastic Differential Equations We would like to introduce stochastic ODE s without going first through the machinery of stochastic integrals. 5.1 Itô Integrals and Itô Differential Equations
More informationEstimating Accuracy in Classical Molecular Simulation
Estimating Accuracy in Classical Molecular Simulation University of Illinois Urbana-Champaign Department of Computer Science Institute for Mathematics and its Applications July 2007 Acknowledgements Ben
More informationNonlinear representation, backward SDEs, and application to the Principal-Agent problem
Nonlinear representation, backward SDEs, and application to the Principal-Agent problem Ecole Polytechnique, France April 4, 218 Outline The Principal-Agent problem Formulation 1 The Principal-Agent problem
More informationHeat bath models for nonlinear partial differential equations Controlling statistical properties in extended dynamical systems
Heat bath models for nonlinear partial differential equations Controlling statistical properties in extended dynamical systems Ben Leimkuhler University of Edinburgh Funding: NAIS (EPSRC) & NWO MIGSEE
More informationStochastic Variational Integrators
Stochastic Variational Integrators Jeremy Schmitt Variational integrators are numerical geometric integrator s derived from discretizing Hamilton s principle. They are symplectic integrators that exhibit
More informationSpeeding up Convergence to Equilibrium for Diffusion Processes
Speeding up Convergence to Equilibrium for Diffusion Processes G.A. Pavliotis Department of Mathematics Imperial College London Joint Work with T. Lelievre(CERMICS), F. Nier (CERMICS), M. Ottobre (Warwick),
More informationHandbook of Stochastic Methods
C. W. Gardiner Handbook of Stochastic Methods for Physics, Chemistry and the Natural Sciences Third Edition With 30 Figures Springer Contents 1. A Historical Introduction 1 1.1 Motivation I 1.2 Some Historical
More informationWiener Measure and Brownian Motion
Chapter 16 Wiener Measure and Brownian Motion Diffusion of particles is a product of their apparently random motion. The density u(t, x) of diffusing particles satisfies the diffusion equation (16.1) u
More informationNUMERICAL METHODS FOR SECOND-ORDER STOCHASTIC DIFFERENTIAL EQUATIONS
SIAM J SCI COMPUT Vol 9, No, pp 45 64 c 7 Society for Industrial and Applied Mathematics NUMERICAL METHODS FOR SECOND-ORDER STOCHASTIC DIFFERENTIAL EQUATIONS KEVIN BURRAGE, IAN LENANE, AND GRANT LYTHE
More informationDISCRETE-TIME STOCHASTIC MODELS, SDEs, AND NUMERICAL METHODS. Ed Allen. NIMBioS Tutorial: Stochastic Models With Biological Applications
DISCRETE-TIME STOCHASTIC MODELS, SDEs, AND NUMERICAL METHODS Ed Allen NIMBioS Tutorial: Stochastic Models With Biological Applications University of Tennessee, Knoxville March, 2011 ACKNOWLEDGEMENT I thank
More informationLecture 12: Detailed balance and Eigenfunction methods
Miranda Holmes-Cerfon Applied Stochastic Analysis, Spring 2015 Lecture 12: Detailed balance and Eigenfunction methods Readings Recommended: Pavliotis [2014] 4.5-4.7 (eigenfunction methods and reversibility),
More information1. Stochastic Processes and filtrations
1. Stochastic Processes and 1. Stoch. pr., A stochastic process (X t ) t T is a collection of random variables on (Ω, F) with values in a measurable space (S, S), i.e., for all t, In our case X t : Ω S
More informationIntroduction. Stochastic Processes. Will Penny. Stochastic Differential Equations. Stochastic Chain Rule. Expectations.
19th May 2011 Chain Introduction We will Show the relation between stochastic differential equations, Gaussian processes and methods This gives us a formal way of deriving equations for the activity of
More informationMean-square Stability Analysis of an Extended Euler-Maruyama Method for a System of Stochastic Differential Equations
Mean-square Stability Analysis of an Extended Euler-Maruyama Method for a System of Stochastic Differential Equations Ram Sharan Adhikari Assistant Professor Of Mathematics Rogers State University Mathematical
More informationState Space Representation of Gaussian Processes
State Space Representation of Gaussian Processes Simo Särkkä Department of Biomedical Engineering and Computational Science (BECS) Aalto University, Espoo, Finland June 12th, 2013 Simo Särkkä (Aalto University)
More informationIntroduction to nonequilibrium physics
Introduction to nonequilibrium physics Jae Dong Noh December 18, 2016 Preface This is a note for the lecture given in the 2016 KIAS-SNU Physics Winter Camp which is held at KIAS in December 17 23, 2016.
More informationLecture 12: Detailed balance and Eigenfunction methods
Lecture 12: Detailed balance and Eigenfunction methods Readings Recommended: Pavliotis [2014] 4.5-4.7 (eigenfunction methods and reversibility), 4.2-4.4 (explicit examples of eigenfunction methods) Gardiner
More informationLecture 4: Numerical Solution of SDEs, Itô Taylor Series, Gaussian Process Approximations
Lecture 4: Numerical Solution of SDEs, Itô Taylor Series, Gaussian Process Approximations Simo Särkkä Aalto University Tampere University of Technology Lappeenranta University of Technology Finland November
More informationNumerical methods for solving stochastic differential equations
Mathematical Communications 4(1999), 251-256 251 Numerical methods for solving stochastic differential equations Rózsa Horváth Bokor Abstract. This paper provides an introduction to stochastic calculus
More informationStochastic Volatility and Correction to the Heat Equation
Stochastic Volatility and Correction to the Heat Equation Jean-Pierre Fouque, George Papanicolaou and Ronnie Sircar Abstract. From a probabilist s point of view the Twentieth Century has been a century
More informationOn the controllability of conservative systems
On the controllability of conservative systems Martin Hairer Department of Mathematics, The University of Warwick, Coventry CV4 7AL, United Kingdom. Email: hairer@maths.warwick.co.uk Abstract We consider
More informationStochastic Calculus February 11, / 33
Martingale Transform M n martingale with respect to F n, n =, 1, 2,... σ n F n (σ M) n = n 1 i= σ i(m i+1 M i ) is a Martingale E[(σ M) n F n 1 ] n 1 = E[ σ i (M i+1 M i ) F n 1 ] i= n 2 = σ i (M i+1 M
More informationAn Introduction to Malliavin calculus and its applications
An Introduction to Malliavin calculus and its applications Lecture 3: Clark-Ocone formula David Nualart Department of Mathematics Kansas University University of Wyoming Summer School 214 David Nualart
More informationLecture 4: Numerical Solution of SDEs, Itô Taylor Series, Gaussian Approximations
Lecture 4: Numerical Solution of SDEs, Itô Taylor Series, Gaussian Approximations Simo Särkkä Aalto University, Finland November 18, 2014 Simo Särkkä (Aalto) Lecture 4: Numerical Solution of SDEs November
More informationThe Wiener Itô Chaos Expansion
1 The Wiener Itô Chaos Expansion The celebrated Wiener Itô chaos expansion is fundamental in stochastic analysis. In particular, it plays a crucial role in the Malliavin calculus as it is presented in
More informationDynamical systems with Gaussian and Levy noise: analytical and stochastic approaches
Dynamical systems with Gaussian and Levy noise: analytical and stochastic approaches Noise is often considered as some disturbing component of the system. In particular physical situations, noise becomes
More informationFinite element approximation of the stochastic heat equation with additive noise
p. 1/32 Finite element approximation of the stochastic heat equation with additive noise Stig Larsson p. 2/32 Outline Stochastic heat equation with additive noise du u dt = dw, x D, t > u =, x D, t > u()
More informationLinear Response Theory, The Green-Kubo formula and Langevin Dynamics
Linear Response Theory, The Green-Kubo formula and Langevin Dynamics G.A. Pavliotis Department of Mathematics Imperial College London 07/07/2015 G.A. Pavliotis (IC London) Linear Response Theory 1 / 64
More informationGaussian Processes. 1. Basic Notions
Gaussian Processes 1. Basic Notions Let T be a set, and X : {X } T a stochastic process, defined on a suitable probability space (Ω P), that is indexed by T. Definition 1.1. We say that X is a Gaussian
More informationLecture 1: Pragmatic Introduction to Stochastic Differential Equations
Lecture 1: Pragmatic Introduction to Stochastic Differential Equations Simo Särkkä Aalto University, Finland (visiting at Oxford University, UK) November 13, 2013 Simo Särkkä (Aalto) Lecture 1: Pragmatic
More informationInterest Rate Models:
1/17 Interest Rate Models: from Parametric Statistics to Infinite Dimensional Stochastic Analysis René Carmona Bendheim Center for Finance ORFE & PACM, Princeton University email: rcarmna@princeton.edu
More informationHandbook of Stochastic Methods
Springer Series in Synergetics 13 Handbook of Stochastic Methods for Physics, Chemistry and the Natural Sciences von Crispin W Gardiner Neuausgabe Handbook of Stochastic Methods Gardiner schnell und portofrei
More informationPhysics 106a, Caltech 13 November, Lecture 13: Action, Hamilton-Jacobi Theory. Action-Angle Variables
Physics 06a, Caltech 3 November, 08 Lecture 3: Action, Hamilton-Jacobi Theory Starred sections are advanced topics for interest and future reference. The unstarred material will not be tested on the final
More informationIntroduction to the Numerical Simulation of Stochastic Differential Equations with Examples
Introduction to the Numerical Simulation of with Examples Prof. Michael Mascagni Department of Computer Science Department of Mathematics Department of Scientific Computing Florida State University, Tallahassee,
More informationStochastic Numerical Analysis
Stochastic Numerical Analysis Prof. Mike Giles mike.giles@maths.ox.ac.uk Oxford University Mathematical Institute Stoch. NA, Lecture 3 p. 1 Multi-dimensional SDEs So far we have considered scalar SDEs
More informationPoisson Jumps in Credit Risk Modeling: a Partial Integro-differential Equation Formulation
Poisson Jumps in Credit Risk Modeling: a Partial Integro-differential Equation Formulation Jingyi Zhu Department of Mathematics University of Utah zhu@math.utah.edu Collaborator: Marco Avellaneda (Courant
More informationStochastic contraction BACS Workshop Chamonix, January 14, 2008
Stochastic contraction BACS Workshop Chamonix, January 14, 2008 Q.-C. Pham N. Tabareau J.-J. Slotine Q.-C. Pham, N. Tabareau, J.-J. Slotine () Stochastic contraction 1 / 19 Why stochastic contraction?
More information5 Applying the Fokker-Planck equation
5 Applying the Fokker-Planck equation We begin with one-dimensional examples, keeping g = constant. Recall: the FPE for the Langevin equation with η(t 1 )η(t ) = κδ(t 1 t ) is = f(x) + g(x)η(t) t = x [f(x)p
More informationLangevin Dynamics of a Single Particle
Overview Langevin Dynamics of a Single Particle We consider a spherical particle of radius r immersed in a viscous fluid and suppose that the dynamics of the particle depend on two forces: A drag force
More informationWeek 9 Generators, duality, change of measure
Week 9 Generators, duality, change of measure Jonathan Goodman November 18, 013 1 Generators This section describes a common abstract way to describe many of the differential equations related to Markov
More informationOn the construction of discrete gradients
On the construction of discrete gradients Elizabeth L. Mansfield School of Mathematics, Statistics and Actuarial Science, University of Kent, Canterbury, CT2 7NF, U.K. G. Reinout W. Quispel Department
More informationMarkovian Approximation and Linear Response Theory for Classical Open Systems
Markovian Approximation and Linear Response Theory for Classical Open Systems G.A. Pavliotis INRIA-MICMAC Project Team and Department of Mathematics Imperial College London Lecture 1 09/11/2011 Lecture
More informationHigher order weak approximations of stochastic differential equations with and without jumps
Higher order weak approximations of stochastic differential equations with and without jumps Hideyuki TANAKA Graduate School of Science and Engineering, Ritsumeikan University Rough Path Analysis and Related
More informationEdinburgh Research Explorer
Edinburgh Research Explorer Adaptive Thermostats for Noisy Gradient Systems Citation for published version: Leimkuhler, B & Shang, X 2016, 'Adaptive Thermostats for Noisy Gradient Systems' SIAM Journal
More informationAdaptive timestepping for SDEs with non-globally Lipschitz drift
Adaptive timestepping for SDEs with non-globally Lipschitz drift Mike Giles Wei Fang Mathematical Institute, University of Oxford Workshop on Numerical Schemes for SDEs and SPDEs Université Lille June
More informationModified Equations for Variational Integrators
Modified Equations for Variational Integrators Mats Vermeeren Technische Universität Berlin Groningen December 18, 2018 Mats Vermeeren (TU Berlin) Modified equations for variational integrators December
More informationThe Smoluchowski-Kramers Approximation: What model describes a Brownian particle?
The Smoluchowski-Kramers Approximation: What model describes a Brownian particle? Scott Hottovy shottovy@math.arizona.edu University of Arizona Applied Mathematics October 7, 2011 Brown observes a particle
More informationUniformly Uniformly-ergodic Markov chains and BSDEs
Uniformly Uniformly-ergodic Markov chains and BSDEs Samuel N. Cohen Mathematical Institute, University of Oxford (Based on joint work with Ying Hu, Robert Elliott, Lukas Szpruch) Centre Henri Lebesgue,
More informationEffective dynamics for the (overdamped) Langevin equation
Effective dynamics for the (overdamped) Langevin equation Frédéric Legoll ENPC and INRIA joint work with T. Lelièvre (ENPC and INRIA) Enumath conference, MS Numerical methods for molecular dynamics EnuMath
More informationBrownian motion and the Central Limit Theorem
Brownian motion and the Central Limit Theorem Amir Bar January 4, 3 Based on Shang-Keng Ma, Statistical Mechanics, sections.,.7 and the course s notes section 6. Introduction In this tutorial we shall
More informationA Concise Course on Stochastic Partial Differential Equations
A Concise Course on Stochastic Partial Differential Equations Michael Röckner Reference: C. Prevot, M. Röckner: Springer LN in Math. 1905, Berlin (2007) And see the references therein for the original
More informationPart II. Classical Dynamics. Year
Part II Year 28 27 26 25 24 23 22 21 20 2009 2008 2007 2006 2005 28 Paper 1, Section I 8B Derive Hamilton s equations from an action principle. 22 Consider a two-dimensional phase space with the Hamiltonian
More informationBoltzmann-Gibbs Preserving Langevin Integrators
Boltzmann-Gibbs Preserving Langevin Integrators Nawaf Bou-Rabee Applied and Comp. Math., Caltech INI Workshop on Markov-Chain Monte-Carlo Methods March 28, 2008 4000 atom cluster simulation Governing Equations
More informationEXPONENTIAL MEAN-SQUARE STABILITY OF NUMERICAL SOLUTIONS TO STOCHASTIC DIFFERENTIAL EQUATIONS
London Mathematical Society ISSN 1461 157 EXPONENTIAL MEAN-SQUARE STABILITY OF NUMERICAL SOLUTIONS TO STOCHASTIC DIFFERENTIAL EQUATIONS DESMOND J. HIGHAM, XUERONG MAO and ANDREW M. STUART Abstract Positive
More informationContinuous dependence estimates for the ergodic problem with an application to homogenization
Continuous dependence estimates for the ergodic problem with an application to homogenization Claudio Marchi Bayreuth, September 12 th, 2013 C. Marchi (Università di Padova) Continuous dependence Bayreuth,
More informationStochastic Differential Equations Related to Soft-Edge Scaling Limit
Stochastic Differential Equations Related to Soft-Edge Scaling Limit Hideki Tanemura Chiba univ. (Japan) joint work with Hirofumi Osada (Kyushu Unv.) 2012 March 29 Hideki Tanemura (Chiba univ.) () SDEs
More informationANOMALOUS TRANSPORT IN RANDOM MEDIA: A ONE-DIMENSIONAL GAUSSIAN MODEL FOR ANOMALOUS DIFFUSION
THESIS FOR THE DEGREE OF MASTER OF SCIENCE ANOMALOUS TRANSPORT IN RANDOM MEDIA: A ONE-DIMENSIONAL GAUSSIAN MODEL FOR ANOMALOUS DIFFUSION ERIK ARVEDSON Department of Theoretical Physics CHALMERS UNIVERSITY
More informationGaussian processes for inference in stochastic differential equations
Gaussian processes for inference in stochastic differential equations Manfred Opper, AI group, TU Berlin November 6, 2017 Manfred Opper, AI group, TU Berlin (TU Berlin) inference in SDE November 6, 2017
More informationStochastic Modelling in Climate Science
Stochastic Modelling in Climate Science David Kelly Mathematics Department UNC Chapel Hill dtbkelly@gmail.com November 16, 2013 David Kelly (UNC) Stochastic Climate November 16, 2013 1 / 36 Why use stochastic
More informationSmoluchowski Diffusion Equation
Chapter 4 Smoluchowski Diffusion Equation Contents 4. Derivation of the Smoluchoswki Diffusion Equation for Potential Fields 64 4.2 One-DimensionalDiffusoninaLinearPotential... 67 4.2. Diffusion in an
More informationLecture 21: Physical Brownian Motion II
Lecture 21: Physical Brownian Motion II Scribe: Ken Kamrin Department of Mathematics, MIT May 3, 25 Resources An instructie applet illustrating physical Brownian motion can be found at: http://www.phy.ntnu.edu.tw/jaa/gas2d/gas2d.html
More informationUnbiased simulation of stochastic differential equations using parametrix expansions
Submitted to Bernoulli Unbiased simulation of stochastic differential equations using parametrix expansions PATRIK ANDERSSON and ARTURO KOHATSU-HIGA Department of Statistics, Uppsala University, Box 53,
More information1 Introduction. 2 Diffusion equation and central limit theorem. The content of these notes is also covered by chapter 3 section B of [1].
1 Introduction The content of these notes is also covered by chapter 3 section B of [1]. Diffusion equation and central limit theorem Consider a sequence {ξ i } i=1 i.i.d. ξ i = d ξ with ξ : Ω { Dx, 0,
More informationanalysis for transport equations and applications
Multi-scale analysis for transport equations and applications Mihaï BOSTAN, Aurélie FINOT University of Aix-Marseille, FRANCE mihai.bostan@univ-amu.fr Numerical methods for kinetic equations Strasbourg
More informationLocal vs. Nonlocal Diffusions A Tale of Two Laplacians
Local vs. Nonlocal Diffusions A Tale of Two Laplacians Jinqiao Duan Dept of Applied Mathematics Illinois Institute of Technology Chicago duan@iit.edu Outline 1 Einstein & Wiener: The Local diffusion 2
More informationStochastic continuity equation and related processes
Stochastic continuity equation and related processes Gabriele Bassi c Armando Bazzani a Helmut Mais b Giorgio Turchetti a a Dept. of Physics Univ. of Bologna, INFN Sezione di Bologna, ITALY b DESY, Hamburg,
More informationLECTURES ON MEAN FIELD GAMES: II. CALCULUS OVER WASSERSTEIN SPACE, CONTROL OF MCKEAN-VLASOV DYNAMICS, AND THE MASTER EQUATION
LECTURES ON MEAN FIELD GAMES: II. CALCULUS OVER WASSERSTEIN SPACE, CONTROL OF MCKEAN-VLASOV DYNAMICS, AND THE MASTER EQUATION René Carmona Department of Operations Research & Financial Engineering PACM
More informationLONG TERM BEHAVIOUR OF LARGE MECHANICAL SYSTEMS WITH RANDOM INITIAL DATA
Stochastics and Dynamics c World Scientific Publishing Company LONG TERM BEHAVIOUR OF LARGE MECHANICAL SYSTEMS WITH RANDOM INITIAL DATA R. KUPFERMAN Institute of Mathematics, The Hebrew University Jerusalem
More informationOn continuous time contract theory
Ecole Polytechnique, France Journée de rentrée du CMAP, 3 octobre, 218 Outline 1 2 Semimartingale measures on the canonical space Random horizon 2nd order backward SDEs (Static) Principal-Agent Problem
More informationGillespie s Algorithm and its Approximations. Des Higham Department of Mathematics and Statistics University of Strathclyde
Gillespie s Algorithm and its Approximations Des Higham Department of Mathematics and Statistics University of Strathclyde djh@maths.strath.ac.uk The Three Lectures 1 Gillespie s algorithm and its relation
More informationLAN property for sde s with additive fractional noise and continuous time observation
LAN property for sde s with additive fractional noise and continuous time observation Eulalia Nualart (Universitat Pompeu Fabra, Barcelona) joint work with Samy Tindel (Purdue University) Vlad s 6th birthday,
More informationConvoluted Brownian motions: a class of remarkable Gaussian processes
Convoluted Brownian motions: a class of remarkable Gaussian processes Sylvie Roelly Random models with applications in the natural sciences Bogotá, December 11-15, 217 S. Roelly (Universität Potsdam) 1
More informationNumerical treatment of stochastic delay differential equations
Numerical treatment of stochastic delay differential equations Evelyn Buckwar Department of Mathematics Heriot-Watt University, Edinburgh Evelyn Buckwar Roma, 16.10.2007 Overview An example: machine chatter
More information