Markovian Approximation and Linear Response Theory for Classical Open Systems
|
|
- Sheila Garrison
- 5 years ago
- Views:
Transcription
1 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 2 28/11/2011 G.A. Pavliotis (MICMAC-IC) Open Classical Systems 1 / 67
2 Contents The Langevin equation. Classical open systems and the Kac-Zwanzig model. Markovian approximation. Derivation of the Lagevin equation. Linear response theory for the Langevin equation. The Green-Kubo formalism. G.A. Pavliotis (MICMAC-IC) Open Classical Systems 2 / 67
3 Consider the Langevin equation q = V (q) γ q + 2γβ 1Ẇ, (1) with V (q) being a confining potential (later on we will also consider periodic potentials), γ is the friction coefficient and β the inverse temperature. Write it as the first order system dq t = p t dt, (2a) dp t = V (q t ) dt γp t dt + 2γβ 1 dw t. (2b) G.A. Pavliotis (MICMAC-IC) Open Classical Systems 3 / 67
4 The process (q t, p t ) is a Markov process on R d R d (or T d R d for periodic potentials) with generator L = p q q V (q) p + γ( p p + β 1 p ). (3) The Fokker-Planck operator (L 2 (R 2d ) adjoint of the generator) is L = p q + q V p +γ ( p (p ) + β 1 p ) (4) The law of the process at time t (distribution function) ρ(q, p, t) is the solution of the Fokker-Planck equation ρ t = L ρ, ρ t=0 = ρ 0. (5) G.A. Pavliotis (MICMAC-IC) Open Classical Systems 4 / 67
5 Introduce the notation X q,p t := ( q t, p t ; q 0 = q, p 0 = p ). The expectation of an observable u(q, p, t) = Ef (X q,p t ), can be computed by solving the backward Kolmogorov equation u t = Lu, u t=0 = f (q, p). (6) G.A. Pavliotis (MICMAC-IC) Open Classical Systems 5 / 67
6 For sufficiently nice potentials the Langevin dynamics (2) is ergodic with respect to the Gibbs measure µ(dqdp) = 1 Z e βh(q,p) dqdp =: ρ (q, p) dpdq, Z = e βh(q,p) dqdp. R 2d The density ρ can be obtained as the solution of the stationary Fokker-Planck equation L ρ = 0. (7a) (7b) Furthermore, we have exponentially fast convergence to equilibrium ρ(t, ) ρ Ce λt, (8) for positive constants C, λ. See the lecture notes by F. Nier. G.A. Pavliotis (MICMAC-IC) Open Classical Systems 6 / 67
7 Why use the Langevin equation? As a thermostat, i.e. to sample from the distribution (7) and to also compute time dependent quantities such as correlation functions. As a reduced description of a more complicated system. Think of the physical model of Brownian motion: colloid particle interacting with the molecules of the surrounding fluid. We obtain a closed (stochastic) equation for the dynamics of the Brownian particle by integrating out the fluid molecules. One of the goals of these lectures is to derive the Langevin equation from a more basic model ( first principles ). We also want to give a miscroscopic definition of the friction coefficient γ. G.A. Pavliotis (MICMAC-IC) Open Classical Systems 7 / 67
8 It is often useful (also in experimental set-ups) to probe a system which is at equilibrium by adding a weak external forcing. The perturbed Langevin dynamics (assumed to be at equilibrium at t = 0) is q ɛ = V (q ɛ ) + ɛf(t) γ q ɛ + 2γβ 1Ẇ. (9) We are interested in understanding the dynamics in the weak forcing regime ɛ 1. We will do this by developing linear response theory (first order perturbation theory). We will see that the response of the system to the external forcing is related to the fluctuations at equilibrium. A result of this analysis is the derivation of Green-Kubo formulas for transport coefficients. For the diffusion coefficient we have (Einstein relation/green-kubo) D = β 1 µ = + V where µ = lim F 0 F, where V is the effective drift. 0 p t p 0 eq dt, (10) G.A. Pavliotis (MICMAC-IC) Open Classical Systems 8 / 67
9 REFERENCES K. Lindenberg and B.J. West The nonequilibrium statistical mechanics of open and closed systems, R.M. Mazo Brownian motion, R. Zwanzig nonequilibrium statistical mechanics, R. Kubo, M. Toda, N. Hashitsume Statistical Physics II. Nonequilibrium statistical mechanics, L. Rey-Bellet Open Classical Systems. In Open Quantum Systems II, D. Givon, R. Kupferman, A.M. Stuart, Extracting macroscopic dynamics: model problems and algorithms, Nonlinearity, 17(6), R55 R127, (2004). Lecture notes, G.P. G.A. Pavliotis (MICMAC-IC) Open Classical Systems 9 / 67
10 DERIVATION OF THE LANGEVIN EQUATION G.A. Pavliotis (MICMAC-IC) Open Classical Systems 10 / 67
11 We will consider the dynamics of a Brownian particle ( small system ) in contact with its environment ( heat bath ). We assume that the environment is at equilibrium at time t = 0. We assume that the dynamics of the full system is Hamiltonian: H(Q, P; q, p) = H BP (Q, P) + H(q, p) + λh I (Q, q), (11) where Q, P are the position and momentum of the Brownian particle, q, p of the environment (they are actually fields) and λ controls the strength of the coupling between the Brownian particle and the environment. G.A. Pavliotis (MICMAC-IC) Open Classical Systems 11 / 67
12 Let H denote the phase space of the full dynamics and X := (P, Q, p, q). Our goal is to show that in some appropriate asymptotic limit we can obtain a closed equation for the dynamics of the Brownian particle. In other words, we need to find a projection P : H R 2n so that PX satisfies an equation that is of the Langevin type (2). We can use projection operator techniques at the level of the Liouville equation ρ = {H, ρ} (12) t associated to the Hamiltonian dynamics (11) Mori-Zwanzig formalism. Projection operator techniques tend to give formal results. Additional approximation/asymptotic calculations are needed. G.A. Pavliotis (MICMAC-IC) Open Classical Systems 12 / 67
13 We will consider a model for which we can obtain the Langevin equation (2) from the Hamiltonian dynamics (11) in a quite explicit way. We will make two basic assumptions: 1 The coupling between the Brownian particle and the environment is linear. 2 The Hamiltonian describing the environment is quadratic in positions (and momenta). In addition, the environment is much larger than the Brownian particle (i.e. infinite dimensional) and at equilibrium at time t = 0. G.A. Pavliotis (MICMAC-IC) Open Classical Systems 13 / 67
14 We have 4 different levels of description: 1 The full Hamiltonian dynamics (11). 2 The generalized Langevin equation (GLE) that we will obtain after eliminating the heat bath variables: t Q = V (Q) γ(t s) Q(s) ds + F(t). (13) 0 3 The Langevin equation (1) that we will obtain in the rapid decorrelation limit 4 The overdamped Langevin (Smoluchowski) dynamics that we obtain in the high friction limit q = V (q) + 2β 1Ẇ. (14) Each reduced level of description includes less information but is easier to analyze. For example, we can prove exponentially fast convergence to equilibrium for (1) and (14), something that is not known, without additional assumptions, for (11) and (13). G.A. Pavliotis (MICMAC-IC) Open Classical Systems 14 / 67
15 The Hamiltonian of the Brownian particle (in one dimension, for simplicity) is described by the Hamiltonian (notice change of notation) H BP = 1 2 p2 + V (q), (15) where V is a confining potential. The environment is modeled as a linear the wave equation (with infinite energy): 2 t φ(t, x) = 2 x φ(t, x). (16) The Hamiltonian of this system (which is quadratic) is H HB (φ, π) = ( x φ 2 + π(x) 2) dx. (17) where π(x) denotes the conjugate momentum field. The initial conditions are distributed according to the Gibbs measure (which in this case is a Gaussian measure) at inverse temperature β, which we formally write as µ β = Z 1 e βh HB(φ,π) dφdπ. (18) G.A. Pavliotis (MICMAC-IC) Open Classical Systems 15 / 67
16 Under this assumption on the initial conditions, typical configurations of the heat bath have infinite energy. In this way, the environment can pump enough energy into the system so that non-trivial fluctuations emerge. The coupling between the particle and the field is linear: H I (Q, φ) = q x φ(x)ρ(x) dx, (19) where The function ρ(x) models the coupling between the particle and the field. The coupling (19) can be thought of as the first term in a Taylor series expansion of a nonlinear, nonlocal coupling (dipole coupling approximation). The Hamiltonian of the particle-field model is H(Q, P, φ, π) = H BP (P, Q) + H(φ, π) + λh I (Q, φ). (20) G.A. Pavliotis (MICMAC-IC) Open Classical Systems 16 / 67
17 Let us first consider a caricature of (20) where there is only one oscillator in the environment : [( H(Q, P, q, p) = P2 p V (Q) ) ] 2 ω2 q 2 λqq, (21) Hamilton s equations of motion are: Q + V (Q) = λq, (22a) q + ω 2 (q λq) = 0. (22b) We can solve (22b) using the variation of constants formula. Set z = (q p) T, p = q. Eqn (22b) can be written as where A = dz dt ( 0 1 ω 2 0 = Az + λh(t), (23) ) and h(t) = ( 0 Q(t) ) G.A. Pavliotis (MICMAC-IC) Open Classical Systems 17 / 67
18 The solution of (23) is t z(t) = e At z(0) + λ e A(t s) h(s) ds. 0 We calculate e At = cos(ωt)i + 1 sin(ωt)a, (24) ω Consequently: where I stands for the 2 2 identity matrix. From this we obtain q(t) = q(0) cos(ωt) + p(0) ω sin(ωt) +λ 1 ω t 0 sin(ω(t s))q(s) ds. (25) G.A. Pavliotis (MICMAC-IC) Open Classical Systems 18 / 67
19 Now we substitute (25) into (22a) to obtain a closed equation that describes the dynamics of the Brownian particle: where t Q = V (Q) λ 2 D(t s)q(s) ds + λf(t), (26) 0 D(t) = 1 sin(ωt), (27) ω and F(t) = q(0) cos(ωt) + p(0) sin(ωt). (28) ω Since q(0) and p(0) are Gaussian random variables with mean 0 and q(0) 2 = β 1 ω 2, p(0) 2 = β 1, q(0)p(0) = 0 we have Consequently, F(t)F(s) = β 1 ω 2 cos(ω(t s)) =: β 1 C(t s). D(t) = Ċ(t). (29) This is a form of the fluctuation-dissipation theorem. G.A. Pavliotis (MICMAC-IC) Open Classical Systems 19 / 67
20 We perform an integration by parts in (25): q(t) = (q(0) λω ) 2 Q(0) cos(ωt) + p(0) ω sin(ωt) +λ 1 ω 2 Q(t) λ 1 t ω 2 cos(ω(t s)) Q(s) ds. We substitute this in equation (22a) to obtain the Generalized Langevin Equation (GLE) 0 t Q = V eff (Q) λ 2 γ(t s) Q(s) ds + λf(t), (30) where V eff (Q) = V (Q) λ2 2ω 2 Q 2, 0 γ(t) = 1 cos(ωt), (31a) ω2 [ ( F(t) = q(0) λ ) ω 2 Q(0) cos(ωt) + p(0) ] ω sin(ωt). (31b) G.A. Pavliotis (MICMAC-IC) Open Classical Systems 20 / 67
21 The same calculation can be done for an arbitrary number of harmonic oscillators in the environment. When writing the dynamics of the Brownian particle in the form (30) we need to introduce an effective potential. We have assumed that the initial conditions of the Brownian particle are deterministic, independent of the initial distribution of the environment. i.e. the environment is initially at equilibrium in the absence of the Brownian particle. We can also assume that the environment is initially in equilibrium in the presence of the distinguished particle, i.e. that the initial positions and momenta of the heat bath particles are distributed according to a Gibbs distribution, conditional on the knowledge of {Q 0, P 0 }: µ β (dpdq) = Z 1 e βh eff(q,p,q) dqdp, (32) where H eff (q, p, Q) = [ p ω2 (q λω 2 Q ) 2 ]. (33) G.A. Pavliotis (MICMAC-IC) Open Classical Systems 21 / 67
22 This assumption implies that q(0) = λ ω 2 Q 0 + β 1 ω 2 ξ, p(0) = β 1 η, (34) where the ξ, η are independent N (0, 1) random variables. This assumption ensures that the forcing term in (30) is mean zero. The fluctuation-dissipation theorem takes the form F(t)F(s) = β 1 γ(t s). (35) G.A. Pavliotis (MICMAC-IC) Open Classical Systems 22 / 67
23 We can perform the same calculation for the coupled particle-field model (20). We obtain where Q = V (Q) t 0 A = D(t s)q(s) ds + φ 0, e Lt α, (36) ( 0 1 x 2 0 ), f, h = ( x f 1 x h 1 + f 2 h 2 ) dx, f = (f 1, f 2 ) and ˆα(k) = ( ik ˆρ(k)/k 2, 0 ) in Fourier space. Furthermore D(t) = e Lt Lα, α = Ċ(t). C(t) is the covariance of the Gaussian noise process In particular: F(t) = φ 0, e Lt α. E(F(t)F(s)) = β 1 C(t s) = β 1 ˆρ(k) 2 e ikt dk. G.A. Pavliotis (MICMAC-IC) Open Classical Systems 23 / 67
24 The spectral density of the autocorrelation function of the noise process is the square of the Fourier transform of the density ρ(x) which controls the coupling between the particle and the environment. The GLE (36) is equivalent to the original infinite dimensional Hamiltonian system with random initial conditions. Proving ergodicity, convergence to equilibrium etc for (36) is equivalent to proving ergodicity and convergence to equilibrium for the infinite dimensional Hamiltonian system. G.A. Pavliotis (MICMAC-IC) Open Classical Systems 24 / 67
25 For general coupling functions ρ(x) the GLE (36) describes a non-markovian system. It can be represented as a Markovian system only if we add an infinite number of additional variables. However, for appropriate choices of the coupling function ρ(x) the GLE (36) is equivalent to a Markovian process in a finite dimensional extended phase space. Definition We will say that a stochastic process X t is quasi-markovian if it can be represented as a Markovian stochastic process by adding a finite number of additional variables: There exists a stochastic process Y t so that {X t, Y t } is a Markov process. G.A. Pavliotis (MICMAC-IC) Open Classical Systems 25 / 67
26 Proposition Assume that p(k) = M m=1 c m( ik) m is a polynomial with real coefficients and roots in the upper half plane then the Gaussian process with spectral density p(k) 2 is the solution of the linear SDE ( p ( i d dt ) x t ) = dw t. (37) dt Proof. The solution of (37) is x t = t k(t s) dw (s), k(t) = 1 2π e ikt 1 p(k) dk. We compute E(x t x s ) = e ik(t s) 1 p(k) 2 dk. G.A. Pavliotis (MICMAC-IC) Open Classical Systems 26 / 67
27 Example Take p(k) (ik + α). The spectral density is ˆρ(k) 2 = α π 1 π 2 + α 2. The autocorrelation function is C(t) = α e ikt 1 π k 2 + α 2 dk = e α t. The linear SDE is dx t = αx t dt + 2αdW t. G.A. Pavliotis (MICMAC-IC) Open Classical Systems 27 / 67
28 Consider the GLE t Q = V (Q) λ 2 γ(t s) Q(s) ds + λf(t), (38) with F (t)f(s) = β 1 γ(t s) = β 1 e α(t s). F(t) is the stationary Ornstein-Uhlenbeck process: df dt = αf + 2β 1 α dw dt We can rewrite (38) as a system of SDEs: dq dt dp dt dz dt = P, where Z (0) N (0, β 1 ). 0 = V (Q) + λz,, with F(0) N (0, β 1 ). (39) = αz λp + 2αβ 1 dw dt G.A. Pavliotis (MICMAC-IC) Open Classical Systems 28 / 67,
29 The process {Q(t), P(t), Z (t)} R 3 is Markovian. It is a degenerate Markov process: noise acts directly only on one of the 3 degrees of freedom. The generator of this process is L = p q + ( V (q) + λz ) p ( αz + λp ) z + αβ 1 2 z. It is a hypoelliptic and hypocoercive operator. G.A. Pavliotis (MICMAC-IC) Open Classical Systems 29 / 67
30 Rescale λ λ ɛ, α α ɛ 2 Eqn. (40) becomes dq ɛ = p ɛ dt, (41a) dp ɛ = V (q ɛ ) dt + λ ɛ zɛ dt, (41b) dz ɛ = α ɛ 2 zɛ dt λ 2αβ 1 ɛ pɛ dt + ɛ 2 dw, (41c) In the limit as ɛ 0 we obtain the Langevin equation for q ɛ (t), p ɛ (t). G.A. Pavliotis (MICMAC-IC) Open Classical Systems 30 / 67
31 Proposition Let {q ɛ (t), p ɛ (t), z ɛ (t)} on R 3 be the solution of (41) with V (q) C (T) with stationary initial conditions. Then {q ɛ (t), p ɛ (t)} converge weakly to the solution of the Langevin equation dq = p dt, dp = V (q) dt γp dt + 2γβ 1 dw, (42) where the friction coefficient is given by the formula γ = λ2 α. (43) G.A. Pavliotis (MICMAC-IC) Open Classical Systems 31 / 67
32 We have derived the (generalized) Langevin equation from a particle + field model where the field is at equilibrium at t = 0. The model that we can considered is very specific since the field and the coupling are linear. The GLE can be derived (at least formally) for more general Hamiltonian systems using the Mori-Zwanzig formalism (projection operator techniques for the Liouville equation). The Markovian approximation of the GLE leads to a finite dimensional hypoelliptic diffusion. Conjecture (F. Nier): Assume that L = V + β 1 has compact resolvent. Then so does L == p q q V p + γ( p p + β 1 p ). Formulate a similar conjecture for all Markovian approximations of the GLE. G.A. Pavliotis (MICMAC-IC) Open Classical Systems 32 / 67
33 LINEAR RESPONSE THEORY G.A. Pavliotis (MICMAC-IC) Open Classical Systems 33 / 67
34 Let X t denote a stationary dynamical system with state space X and invariant measure µ(dx) = f (x) dx. We probe the system by adding a time dependent forcing ɛf(t) with ɛ 1 at time t 0. Goal: calculate the distribution function f ɛ (x, t) of the perturbed systems Xt ɛ, ɛ 1, in particular in the long time limit t +. We can then calculate expectation value of observables as well as correlation functions. G.A. Pavliotis (MICMAC-IC) Open Classical Systems 34 / 67
35 We assume that the distribution function f ɛ (x, t) satisfies a linear kinetic equation e.g. the Liouville or the Fokker-Planck equation: f ɛ t = L ɛ f ɛ, (44a) f ɛ t=t0 = f. (44b) The choice of the initial conditions reflects the fact that at t = t 0 the system is at equilibrium. Since f is the unique equilibrium distribution, we have that L 0 f = 0. (45) G.A. Pavliotis (MICMAC-IC) Open Classical Systems 35 / 67
36 The operator L ɛ can be written in the form L ɛ = L 0 + ɛl 1, (46) where cl 0 denotes the Liouville or Fokker-Planck operator of the unperturbed system and L 1 is related to the external forcing. We will assume that L 1 is of the form L 1 where D is some linear (differential) operator. = F(t) D, (47) G.A. Pavliotis (MICMAC-IC) Open Classical Systems 36 / 67
37 Example (A deterministic dynamical system). Let X t be the solution of the differential equation dx t = h(x t ), (48) dt on a (possibly compact) state space X. We add a weak time dependent forcing to obtain the dynamics dx t dt = h(x t ) + ɛf(t). (49) We assume that the unperturbed dynamics has a unique invariant distribution f which is the solution of the stationary Liouville equation ) (h(x)f = 0, (50) equipped with appropriate boundary conditions. G.A. Pavliotis (MICMAC-IC) Open Classical Systems 37 / 67
38 Example The operator L ɛ in (46) has the form L ɛ ( ) = h(x) ɛf (t). In this example, the operator D in (47) is D =. A particular case of a deterministic dynamical system of the form (48), and the most important in statistical mechanics, is that of an N-body Hamiltonian system. G.A. Pavliotis (MICMAC-IC) Open Classical Systems 38 / 67
39 Example (A stochastic dynamical system). Let X t be the solution of the stochastic differential equation dx t = h(x t ) dt + σ(x t ) dw t, (51) on R d, where σ(x) is a positive semidefinite matrix and where the Itô interpretation is used. We add a weak time dependent forcing to obtain the dynamics dx t = h(x t ) dt + ɛf(t) dt + σ(x t ) dw t. (52) We assume that the unperturbed dynamics has a unique invariant distribution f which is the solution of the stationary Fokker-Planck equation ) (h(x)f + 1 ) 2 D2 : (Σ(x)f = 0, (53) where Σ(x) = σ(x)σ T (x). G.A. Pavliotis (MICMAC-IC) Open Classical Systems 39 / 67
40 Example The operator L ɛ in (46) has the form L ɛ ( ) = h(x) + 1 ( ) 2 D2 : Σ(x) ɛf(t). As in the previous example, the operator D in (47) is D =. A particular case of Example 5 is the Langevin equation: q = V (q) + ɛf(t) γ q + 2γβẆ. (54) We can also add a forcing term to the noise in (51), i.e. we can consider a time-dependent temperature. For the Langevin equation the perturbed dynamics is q = V (q) + ɛf(t) γ q + 2γβ 1 (1 + ɛt (t))ẇ, (55) with 1 + ɛt (t) > 0. in this case the operator L 1 is L 1 = F(t) p + γβ 1 T (t) p, where p = q. G.A. Pavliotis (MICMAC-IC) Open Classical Systems 40 / 67
41 We proceed with the analysis of (44). We look for a solution in the form of a power series expansion in ɛ: f ɛ = f 0 + ɛf We substitute this into (44a) and use the initial condition (44b) to obtain the equations f 0 t = L 0 f 0, f 0 t=0 = f, (56a) f 1 = L 0 t 1 + L 1 f 0, f t=0 1 = 0. (56b) The only solution to (56a) is f 0 = f. We use this into (56b) and use (47) to obtain f 1 t = L 0 f 1 + F(t) Df, f 1 t=0 = 0. G.A. Pavliotis (MICMAC-IC) Open Classical Systems 41 / 67
42 We use the variation of constants formula to solve this equation: f 1 (t) = t t 0 e L 0 (t s) F(s) Df ds. (57) Now we can calculate the deviation in the expectation value of an observable due to the external forcing: Let eq and denote the expectation value with respect to f and f ɛ, respectively. Let A( ) be an observable and denote by A(t) the deviation of its expectation value from equilibrium, to leading order: A(t) := A(X t ) A(X t ) eq = A(x) ( f ɛ (x, t) f eq (x) ) dx ( t ) = ɛ A(x) e L (t s) 0 F(s) Df ds dx. t 0 G.A. Pavliotis (MICMAC-IC) Open Classical Systems 42 / 67
43 Assuming now that we can interchange the order of integration we can rewrite the above formula as ( t ) A(t) = ɛ A(x) e L (t s) 0 F(s) Df ds dx t = ɛ =: ɛ t 0 t ( t 0 ) A(x)e L 0 (t s) Df dx ds t 0 R L0,A(t s)f(s) ds, (58) where we have defined the response function R L0,A(t) = A(x)e L 0 t Df dx (59) We set now the lower limit of integration in (58) to be t 0 = (define R L0,A(t) in (59) to be 0 for t < 0) and assume that we can extend the upper limit of integration to + to write A(t) = ɛ + R L0,A(t s)f (s) ds. (60) G.A. Pavliotis (MICMAC-IC) Open Classical Systems 43 / 67
44 As expected (since we have used linear perturbation theory), the deviation of the expectation value of an observable from its equilibrium value is a linear function of the forcing term. (60) has the form of the solution of a linear differential equation with R L0,A(t) playing the role of the Green s function. If we consider a delta-like forcing at t = 0, F(t) = δ(t), then the above formula gives A(t) = ɛr L0,A(t). Thus, the response function gives the deviation of the expectation value of an observable from equilibrium for a delta-like force. Consider a constant force that is exerted to the system at time t = 0, F(t) = FΘ(t) where Θ(t) denotes the Heaviside step function. For this forcing (58) becomes A(t) = ɛf t 0 R L0,A(t s) ds. (61) G.A. Pavliotis (MICMAC-IC) Open Classical Systems 44 / 67
45 There is a close relation between the response function (59) and stationary autocorrelation functions. Let X t be a stationary Markov process in R d with generator L and invariant distribution f. Let A( ) and B( ) be two observables. The stationary autocorrelation function A(X t )B(X 0 ) eq can be calculated as follows κ A,B (t) := A(X t )B(X 0 ) eq = A(x)B(x 0 )p(x, t x 0, 0)f (x 0 ) dxdx 0 = A(x)B(x 0 )e L t δ(x x 0 )f (x 0 ) dxdx 0 = e Lt A(x)B(x 0 )δ(x x 0 )f (x 0 ) dxdx 0 = e Lt A(x)B(x)f (x) dx, where L acts on functions of x. G.A. Pavliotis (MICMAC-IC) Open Classical Systems 45 / 67
46 Thus we have established the formula κ A,B (t) = S t A(x), B(x) f, (62) where S t denotes the semigroup generated by L and, f denotes the L 2 inner product weighted by the invariant distribution of the diffusion process. Consider now the particular choice B(x) = f 1 Df. We combine (59) and (62) to deduce κ A,f 1 Df (t) = R L0,A(t). (63) This is a version of the fluctuation-dissipation theorem. G.A. Pavliotis (MICMAC-IC) Open Classical Systems 46 / 67
47 Example Consider the Langevin dynamics with an external forcing F. dq = p dt, dp = V (q) dt + F dt γp dt + 2γβ 1 dw t. We have D = p and We use (63) with A = p: B = f 1 Df = βp. β p(t)p(0) eq = R L0,p(t). G.A. Pavliotis (MICMAC-IC) Open Classical Systems 47 / 67
48 Example (continued) When the potential is harmonic, V (q) = 1 2 ω2 0 q2, we can compute explicitly the response function and, consequently, the velocity autocorrelation function at equilibrium: and Consequently: R q,l (t) = 1 ω 1 e γt 2 sin(ω1 t), ω 1 = R p,l (t) = e γt 2 p(t)p(0) eq = β 1 e γt 2 ω 2 0 γ2 4 ( cos(ω 1 t) γ ) sin(ω 1 t). 2ω 1 ( cos(ω 1 t) γ ) sin(ω 1 t). 2ω 1 Similar calculations can be done for all linear SDEs. G.A. Pavliotis (MICMAC-IC) Open Classical Systems 48 / 67
49 Example Consider the Langevin dynamics with a perturbation in the temperature dq = p dt, dp = V (q) dt + F dt γp dt + 2γβ 1 (1 + F) dw t. We have D = γβ 1 2 p and B = f 1 Df = γβ(p 2 β 1 ). Let H(p, q) = p 2 /2 + V (q) denote the total energy. We have f 1 L 0 H(p, q)f = L 0 H(p, q) = γ( p 2 + β 1 ). Consequently (using (62)): κ A,f 1 Df (t) = β d dt κ A,H(t). For A = H we obtain R H,L (t) = β d dt H(t)H(0) eq. G.A. Pavliotis (MICMAC-IC) Open Classical Systems 49 / 67
50 We calculate the long time limit of A(t) when the external forcing is a step function. The following formal calculations can be justified e.g. for reversible diffusions using functional calculus. t 0 R L,A (t s) ds = = = = = t 0 A(x)e L 0(t s) Df dx ds t ( ) e L(t s) A(x) Df ds dx 0 ( t ) e L 0t e L( s) dsa(x) Df dx 0 ( ( e L0t ( L 0 ) 1 e L0( t) I) ((I e L 0 t ) ( L) 1 A(x)) Df dx ) A(x) Df dx G.A. Pavliotis (MICMAC-IC) Open Classical Systems 50 / 67
51 Assuming now that lim t + e Lt = 0 (again, think of reversible diffusions) we have that D := t lim R L,A (t s) ds = t + 0 ( L) 1 A(x)Df dx. Using this in (61) and relabeling ɛf F we obtain lim lim A(t) F 0 t + F = ( L) 1 A(x)Df dx. (64) G.A. Pavliotis (MICMAC-IC) Open Classical Systems 51 / 67
52 Remark At least formally, we can interchange the order with which we take the limits in (64). Formulas of the form (64) enable us to calculate transport coefficients, such as the diffusion coefficient. We can rewrite the above formula in the form lim lim A(t) F 0 t + F = φdf dx, (65) where φ is the solution of the Poisson equation (when A(x) is mean zero) Lφ = A(x), (66) equipped with appropriate boundary conditions. This is precisely the formalism that we obtain using homogenization theory. G.A. Pavliotis (MICMAC-IC) Open Classical Systems 52 / 67
53 Consider the Langevin dynamics in a periodic or random potential. dq t = p t dt, dp t = V (q t ) dt γp t dt + 2γβ 1 dw. From Einstein s formula (10) we have that the diffusion coefficient is related to the mobility according to D = β 1 lim where we have used p t eq = 0. lim F 0 t + p t F We use now (65) with A(t) = p t, D = p, f = 1 Z e βh(q,p) to obtain D = φpf dpdq = Lφ, φ f, (67) which is the formula obtained from homogenization theory (e.g. Kipnis and Varadhan 1985, Rodenhausen 1989). G.A. Pavliotis (MICMAC-IC) Open Classical Systems 53 / 67
54 Consider the unperturbed dynamics dq = p dt, dp = V (q) dt γp dt + 2γβ 1 dw, (68) where V is a smooth periodic potential. Then the rescaled process q ɛ (t) := ɛq(t/ɛ 2 ), Converges to a Brownian motion with diffusion coefficient D = (( L) 1 p)pf dpdq. (69) A similar result can be proved when V is a random potential. G.A. Pavliotis (MICMAC-IC) Open Classical Systems 54 / 67
55 Consider now the perturbed dynamics dq F = p F dt, dp F = V (q F ) dt + F dt γp F dt + 2γβ 1 dw. (70) At least for F sufficiently small, the dynamics is ergodic on T R with a smooth invariant density f F (p, q) which is a differentiable function of F. The external forcing induces an effective drift V F = pf F dpdq. (71) The mobility is defined as µ := d df V F. (72) F =0 G.A. Pavliotis (MICMAC-IC) Open Classical Systems 55 / 67
56 Theorem (Rodenhausen J. Stat. Phys. 1989) The mobility is well defined by (72) and D = β 1 µ. (73) Proof. Ergodic theory for hypoelliptic diffusions. Study of the Poisson and stationary Fokker-Planck equations. Girsanov s formula. G.A. Pavliotis (MICMAC-IC) Open Classical Systems 56 / 67
57 Upon combining (63) with (64) we obtain D = t lim t + 0 κ A,f 1 Df (t s) ds. (74) Thus, a transport coefficient can be computed in terms of the time integral of an appropriate autocorrelation function. This is an example of the Green-Kubo formula G.A. Pavliotis (MICMAC-IC) Open Classical Systems 57 / 67
58 We can obtain a more general form of the Green-Kubo formalism. We define the generalized drift and diffusion coefficients as follows : V f 1 ( ) (x) = lim h 0 h E f (X h ) f (X 0 ) X 0 = x = Lf (75) and D f,g (x) := 1 ( ) lim h 0 h E (f (X t+h ) f (X t ))((g(x t+h ) g(x t )) X t = x = L(fg)(x) (glf )(x) (f Lg)(x), (76) where f, g are smooth functions (in fact, all we need is f, g D(L) and fg D(L 0 )). The equality in (75) follows from the definition of the generator of a diffusion process. We will prove the equality in (75). Sometimes D f,g (x) is called the opérateur carré du champ. G.A. Pavliotis (MICMAC-IC) Open Classical Systems 58 / 67
59 To prove (76), notice first that, by stationarity, it is sufficient to prove it at t = 0. Furthermore (f (X h ) f (X 0 )) (g(x h ) g(x 0 )) = (fg)(x h ) (fg)(x 0 ) Consequently: f (X 0 )(g(x h ) g(x 0 )) g(x 0 )(f (X h ) f (X 0 )). D f,g 1 ( ) (x) := lim h 0 h E (f (X h ) f (X 0 ))((g(x h ) g(x 0 )) X 0 = x 1 ( ) = lim h 0 h E (fg)(x h ) (fg)(x 0 ) X 0 = x 1 ( ) lim h 0 h E g(x h ) g(x 0 ) X 0 = x f (x) 1 ( ) lim h 0 h E f (X h ) f (X 0 ) X 0 = x g(x) = L(fg)(x) (glf )(x) (f Lg)(x). G.A. Pavliotis (MICMAC-IC) Open Classical Systems 59 / 67
60 We have the following result. Theorem (The Green-Kubo formula) Let X t be a stationary reversible diffusion process with state space X, generator L, invariant measure µ(dx) and let V f (x), D f,g (x). Then 1 2 D f,g µ(dx) = 0 ( ) E V f (X t )V g (X 0 ) dt. (77) G.A. Pavliotis (MICMAC-IC) Open Classical Systems 60 / 67
61 Proof. Let (, ) µ denote the inner product in L 2 (X, µ). We note that 1 D f,g µ(dx) = ( Lf, g) µ, (78) 2 In view of (78), formula (77) becomes ( Lf, g) µ = 0 ( ) E V f (X t )V g (X 0 ) dt. (79) Now we use (62), together with the identity 0 elt dt = ( L) 1 to obtain 0 κ A,B (t) dt = (( L) 1 A, B) µ. We set now A = V f = Lf, B = V g = Lg in the above formula to deduce (79) from which (77) follows. G.A. Pavliotis (MICMAC-IC) Open Classical Systems 61 / 67
62 Remark The above formal calculation can also be performed in the nonreversible case to obtain 1 ( ) D f,f µ(dx) = E V f (X t )V f (X 0 ) dt. 2 In the reversible case (i.e. L being a selfadjoint operator in H := L 2 (X, µ)) the formal calculations can be justified using functional calculus. For the reversible diffusion process 0 dx t = V (X t ) dt + 2β 1 dw t and under appropriate assumptions on the potential, the generator L has compact resolvent and its eigenfunctions form an orthonormal basis in H. In this case the proof of (11) becomes quite simple. G.A. Pavliotis (MICMAC-IC) Open Classical Systems 62 / 67
63 The spectral representation of L is from which it follows that e Lt = 0 L = 0 λ de λ, e λt de λ and e Lt L = 0 e λt λ. Remember that V f = Lf We calculate (with S t = e Lt ) E[V f (X t )V g (X 0 )] = (S t Lf, Lg) µ = 0 λ 2 e λt d(e λ f, g) µ. (80) From Fubini s theorem and Cauchy-Schwarz it follows that 0 0 with D L (f ) := ( Lf, f ) µ. λ 2 e λt d(e λ f, g) µ dt D L (f ) 1/2 D L (g) 1/2 < + G.A. Pavliotis (MICMAC-IC) Open Classical Systems 63 / 67
64 We use now (80) and Fubini s theorem to calculate + 0 E[V f (X t )V g (X 0 )] dt = = t (S t Lf, Lg) µ dt λ 2 e λt d(e λ f, g) µ dt = ( λ)d(e λ f, g) µ ( 0 ) 0 = ( λ)de λ f, de λ g µ = ( Lf, g) µ = 1 D f,g (x)µ(dx). 2 G.A. Pavliotis (MICMAC-IC) Open Classical Systems 64 / 67
65 Example Consider the diffusion process dx t = b(x t ) dt + σ(x t ) dw t, (81) where the noise is interpreted in the Stratonovich sense. The generator is L = b(x) + 1 (A(x) ), (82) 2 where A(x) = (σσ T )(x). We assume that the diffusion process has a unique invariant distribution which is the solution of the stationary Fokker-Planck equation L ρ = 0. (83) When A is strictly positive definite ergodicity follows from PDE arguments. G.A. Pavliotis (MICMAC-IC) Open Classical Systems 65 / 67
66 Example (Continued) The stationary process X t (i.e. X 0 ρ(x)dx) is reversible provided that Let f = x i, g = x j. We calculate b(x) = 1 A(x) log ρ(x). (84) 2 V x i (x) = Lx i = b i ka ik, i = 1,... d. (85) We use the detailed balance condition (84) and (78) to calculate ( 1 D x i,x j µ(dx) = = b i (x) + 1 ) 2 2 ka ik (x) x j ρ(x) dx = 1 (A ik k ρ(x) + k A ik (x)ρ(x)) x j dx 2 = 1 A ij (x)ρ(x) dx. 2 G.A. Pavliotis (MICMAC-IC) Open Classical Systems 66 / 67
67 Example (Continued) The Green-Kubo formula (77) gives: 1 2 A ij (x)ρ(x) dx = + where the drift V x i (x) is given by (84). 0 ( ) E V x i (X t )V x j (X 0 ) dt, (86) G.A. Pavliotis (MICMAC-IC) Open Classical Systems 67 / 67
Linear 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 informationNon-Equilibrium Statistical Mechanics. G.A. Pavliotis Department of Mathematics Imperial College London London SW7 2AZ, UK
Non-Equilibrium Statistical Mechanics G.A. Pavliotis Department of Mathematics Imperial College London London SW7 2AZ, UK June 2, 212 2 Contents 1 Derivation of the Langevin Equation 1 1.1 Open Classical
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 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 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 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 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 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 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 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 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 Langevin Equations
1 Brownian Motion and Langevin Equations 1.1 Langevin Equation and the Fluctuation- Dissipation Theorem The theory of Brownian motion is perhaps the simplest approximate way to treat the dynamics of nonequilibrium
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 information10. Zwanzig-Mori Formalism
University of Rhode Island DigitalCommons@URI Nonequilibrium Statistical Physics Physics Course Materials 0-9-205 0. Zwanzig-Mori Formalism Gerhard Müller University of Rhode Island, gmuller@uri.edu Follow
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 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 informationComputing ergodic limits for SDEs
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
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 informationM4A42 APPLIED STOCHASTIC PROCESSES
M4A42 APPLIED STOCHASTIC PROCESSES G.A. Pavliotis Department of Mathematics Imperial College London, UK LECTURE 1 12/10/2009 Lectures: Mondays 09:00-11:00, Huxley 139, Tuesdays 09:00-10:00, Huxley 144.
More informationA path integral approach to the Langevin equation
A path integral approach to the Langevin equation - Ashok Das Reference: A path integral approach to the Langevin equation, A. Das, S. Panda and J. R. L. Santos, arxiv:1411.0256 (to be published in Int.
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 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 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 informationCorrections to Einstein s relation for Brownian motion in a tilted periodic potential
Corrections to Einstein s relation for Brownian motion in a tilted periodic potential J.C. Latorre Department of Mathematics and Computer Science, Freie Universität Berlin, Germany G.A. Pavliotis Department
More informationG : Statistical Mechanics
G25.2651: Statistical Mechanics Notes for Lecture 15 Consider Hamilton s equations in the form I. CLASSICAL LINEAR RESPONSE THEORY q i = H p i ṗ i = H q i We noted early in the course that an ensemble
More informationNPTEL
NPTEL Syllabus Nonequilibrium Statistical Mechanics - Video course COURSE OUTLINE Thermal fluctuations, Langevin dynamics, Brownian motion and diffusion, Fokker-Planck equations, linear response theory,
More information10. Zwanzig-Mori Formalism
University of Rhode Island DigitalCommons@URI Nonequilibrium Statistical Physics Physics Course Materials 205 0. Zwanzig-Mori Formalism Gerhard Müller University of Rhode Island, gmuller@uri.edu Creative
More informationTable of Contents [ntc]
Table of Contents [ntc] 1. Introduction: Contents and Maps Table of contents [ntc] Equilibrium thermodynamics overview [nln6] Thermal equilibrium and nonequilibrium [nln1] Levels of description in statistical
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 informationF r (t) = 0, (4.2) F (t) = r= r a (t)
Chapter 4 Stochastic Equations 4.1 Langevin equation To explain the idea of introduction of the Langevin equation, let us consider the concrete example taken from surface physics: motion of an atom (adatom)
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 informationCHAPTER V. Brownian motion. V.1 Langevin dynamics
CHAPTER V Brownian motion In this chapter, we study the very general paradigm provided by Brownian motion. Originally, this motion is that a heavy particle, called Brownian particle, immersed in a fluid
More informationDerivation of the GENERIC form of nonequilibrium thermodynamics from a statistical optimization principle
Derivation of the GENERIC form of nonequilibrium thermodynamics from a statistical optimization principle Bruce Turkington Univ. of Massachusetts Amherst An optimization principle for deriving nonequilibrium
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 informationSeparation of Variables in Linear PDE: One-Dimensional Problems
Separation of Variables in Linear PDE: One-Dimensional Problems Now we apply the theory of Hilbert spaces to linear differential equations with partial derivatives (PDE). We start with a particular example,
More informationPath integrals for classical Markov processes
Path integrals for classical Markov processes Hugo Touchette National Institute for Theoretical Physics (NITheP) Stellenbosch, South Africa Chris Engelbrecht Summer School on Non-Linear Phenomena in Field
More informationFluctuation theorem in systems in contact with different heath baths: theory and experiments.
Fluctuation theorem in systems in contact with different heath baths: theory and experiments. Alberto Imparato Institut for Fysik og Astronomi Aarhus Universitet Denmark Workshop Advances in Nonequilibrium
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 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 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 information16. Working with the Langevin and Fokker-Planck equations
16. Working with the Langevin and Fokker-Planck equations In the preceding Lecture, we have shown that given a Langevin equation (LE), it is possible to write down an equivalent Fokker-Planck equation
More informationBrownian Motion: Fokker-Planck Equation
Chapter 7 Brownian Motion: Fokker-Planck Equation The Fokker-Planck equation is the equation governing the time evolution of the probability density of the Brownian particla. It is a second order differential
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 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 informationNon-equilibrium phenomena and fluctuation relations
Non-equilibrium phenomena and fluctuation relations Lamberto Rondoni Politecnico di Torino Beijing 16 March 2012 http://www.rarenoise.lnl.infn.it/ Outline 1 Background: Local Thermodyamic Equilibrium 2
More informationMicroscopic Hamiltonian dynamics perturbed by a conservative noise
Microscopic Hamiltonian dynamics perturbed by a conservative noise CNRS, Ens Lyon April 2008 Introduction Introduction Fourier s law : Consider a macroscopic system in contact with two heat baths with
More information08. Brownian Motion. University of Rhode Island. Gerhard Müller University of Rhode Island,
University of Rhode Island DigitalCommons@URI Nonequilibrium Statistical Physics Physics Course Materials 1-19-215 8. Brownian Motion Gerhard Müller University of Rhode Island, gmuller@uri.edu Follow this
More informationConvergence to equilibrium of Markov processes (eventually piecewise deterministic)
Convergence to equilibrium of Markov processes (eventually piecewise deterministic) A. Guillin Université Blaise Pascal and IUF Rennes joint works with D. Bakry, F. Barthe, F. Bolley, P. Cattiaux, R. Douc,
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 informationStatistical mechanics of random billiard systems
Statistical mechanics of random billiard systems Renato Feres Washington University, St. Louis Banff, August 2014 1 / 39 Acknowledgements Collaborators: Timothy Chumley, U. of Iowa Scott Cook, Swarthmore
More information5.74 Introductory Quantum Mechanics II
MIT OpenCourseWare http://ocw.mit.edu 5.74 Introductory Quantum Mechanics II Spring 9 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms. Andrei Tokmakoff,
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 informationStatistical Mechanics of Active Matter
Statistical Mechanics of Active Matter Umberto Marini Bettolo Marconi University of Camerino, Italy Naples, 24 May,2017 Umberto Marini Bettolo Marconi (2017) Statistical Mechanics of Active Matter 2017
More informationStochastic equations for thermodynamics
J. Chem. Soc., Faraday Trans. 93 (1997) 1751-1753 [arxiv 1503.09171] Stochastic equations for thermodynamics Roumen Tsekov Department of Physical Chemistry, University of Sofia, 1164 Sofia, ulgaria The
More informationHierarchical Modeling of Complicated Systems
Hierarchical Modeling of Complicated Systems C. David Levermore Department of Mathematics and Institute for Physical Science and Technology University of Maryland, College Park, MD lvrmr@math.umd.edu presented
More informationHAMILTON S PRINCIPLE
HAMILTON S PRINCIPLE In our previous derivation of Lagrange s equations we started from the Newtonian vector equations of motion and via D Alembert s Principle changed coordinates to generalised coordinates
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 informationIntertwinings for Markov processes
Intertwinings for Markov processes Aldéric Joulin - University of Toulouse Joint work with : Michel Bonnefont - Univ. Bordeaux Workshop 2 Piecewise Deterministic Markov Processes ennes - May 15-17, 2013
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 informationSTOCHASTIC PROCESSES IN PHYSICS AND CHEMISTRY
STOCHASTIC PROCESSES IN PHYSICS AND CHEMISTRY Third edition N.G. VAN KAMPEN Institute for Theoretical Physics of the University at Utrecht ELSEVIER Amsterdam Boston Heidelberg London New York Oxford Paris
More informationDecay rates for partially dissipative hyperbolic systems
Outline Decay rates for partially dissipative hyperbolic systems Basque Center for Applied Mathematics Bilbao, Basque Country, Spain zuazua@bcamath.org http://www.bcamath.org/zuazua/ Numerical Methods
More informationAnomalous Transport and Fluctuation Relations: From Theory to Biology
Anomalous Transport and Fluctuation Relations: From Theory to Biology Aleksei V. Chechkin 1, Peter Dieterich 2, Rainer Klages 3 1 Institute for Theoretical Physics, Kharkov, Ukraine 2 Institute for Physiology,
More informationDensities for the Navier Stokes equations with noise
Densities for the Navier Stokes equations with noise Marco Romito Università di Pisa Universitat de Barcelona March 25, 2015 Summary 1 Introduction & motivations 2 Malliavin calculus 3 Besov bounds 4 Other
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 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 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 informationExample 4.1 Let X be a random variable and f(t) a given function of time. Then. Y (t) = f(t)x. Y (t) = X sin(ωt + δ)
Chapter 4 Stochastic Processes 4. Definition In the previous chapter we studied random variables as functions on a sample space X(ω), ω Ω, without regard to how these might depend on parameters. We now
More informationProblem reduction, memory and renormalization
Problem reduction, memory and renormalization Panos Stinis Northwest Institute for Advanced Computing Pacific Northwest National Laboratory A simple example The linear differential system for x(t) and
More informationFast-slow systems with chaotic noise
Fast-slow systems with chaotic noise David Kelly Ian Melbourne Courant Institute New York University New York NY www.dtbkelly.com May 1, 216 Statistical properties of dynamical systems, ESI Vienna. David
More informationReflected Brownian Motion
Chapter 6 Reflected Brownian Motion Often we encounter Diffusions in regions with boundary. If the process can reach the boundary from the interior in finite time with positive probability we need to decide
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 informationIntroduction to asymptotic techniques for stochastic systems with multiple time-scales
Introduction to asymptotic techniques for stochastic systems with multiple time-scales Eric Vanden-Eijnden Courant Institute Motivating examples Consider the ODE {Ẋ = Y 3 + sin(πt) + cos( 2πt) X() = x
More informationDamped Harmonic Oscillator
Damped Harmonic Oscillator Wednesday, 23 October 213 A simple harmonic oscillator subject to linear damping may oscillate with exponential decay, or it may decay biexponentially without oscillating, or
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 informationThis is a Gaussian probability centered around m = 0 (the most probable and mean position is the origin) and the mean square displacement m 2 = n,or
Physics 7b: Statistical Mechanics Brownian Motion Brownian motion is the motion of a particle due to the buffeting by the molecules in a gas or liquid. The particle must be small enough that the effects
More informationWeak convergence and large deviation theory
First Prev Next Go To Go Back Full Screen Close Quit 1 Weak convergence and large deviation theory Large deviation principle Convergence in distribution The Bryc-Varadhan theorem Tightness and Prohorov
More informationCanonical transformations (Lecture 4)
Canonical transformations (Lecture 4) January 26, 2016 61/441 Lecture outline We will introduce and discuss canonical transformations that conserve the Hamiltonian structure of equations of motion. Poisson
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 informationSimulation methods for stochastic models in chemistry
Simulation methods for stochastic models in chemistry David F. Anderson anderson@math.wisc.edu Department of Mathematics University of Wisconsin - Madison SIAM: Barcelona June 4th, 21 Overview 1. Notation
More informationQUANTUM MARKOVIAN KINETIC EQUATION FOR HARMONIC OSCILLATOR. Boris V. Bondarev
QUANTUM MARKOVIAN KINETIC EQUATION FOR HARMONIC OSCILLATOR Boris V. Bondarev Moscow Aviation Institute, Volokolamsk road, 4, 15871, Moscow, Russia E-mail: bondarev.b@mail.ru arxiv:130.0303v1 [physics.gen-ph]
More informationPhysics 5153 Classical Mechanics. Canonical Transformations-1
1 Introduction Physics 5153 Classical Mechanics Canonical Transformations The choice of generalized coordinates used to describe a physical system is completely arbitrary, but the Lagrangian is invariant
More informationfor changing independent variables. Most simply for a function f(x) the Legendre transformation f(x) B(s) takes the form B(s) = xs f(x) with s = df
Physics 106a, Caltech 1 November, 2018 Lecture 10: Hamiltonian Mechanics I The Hamiltonian In the Hamiltonian formulation of dynamics each second order ODE given by the Euler- Lagrange equation in terms
More informationLecture 6: Multiple Model Filtering, Particle Filtering and Other Approximations
Lecture 6: Multiple Model Filtering, Particle Filtering and Other Approximations Department of Biomedical Engineering and Computational Science Aalto University April 28, 2010 Contents 1 Multiple Model
More informationStatistical methods in atomistic computer simulations
Statistical methods in atomistic computer simulations Prof. Michele Ceriotti, michele.ceriotti@epfl.ch This course gives an overview of simulation techniques that are useful for the computational modeling
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 informationActive Matter Lectures for the 2011 ICTP School on Mathematics and Physics of Soft and Biological Matter Lecture 3: Hydrodynamics of SP Hard Rods
Active Matter Lectures for the 2011 ICTP School on Mathematics and Physics of Soft and Biological Matter Lecture 3: of SP Hard Rods M. Cristina Marchetti Syracuse University Baskaran & MCM, PRE 77 (2008);
More informationIntroduction to Random Diffusions
Introduction to Random Diffusions The main reason to study random diffusions is that this class of processes combines two key features of modern probability theory. On the one hand they are semi-martingales
More informationStrong uniqueness for stochastic evolution equations with possibly unbounded measurable drift term
1 Strong uniqueness for stochastic evolution equations with possibly unbounded measurable drift term Enrico Priola Torino (Italy) Joint work with G. Da Prato, F. Flandoli and M. Röckner Stochastic Processes
More informationHypoelliptic multiscale Langevin diffusions and Slow fast stochastic reaction diffusion equations.
Hypoelliptic multiscale Langevin diffusions and Slow fast stochastic reaction diffusion equations. Wenqing Hu. 1 (Joint works with Michael Salins 2 and Konstantinos Spiliopoulos 3.) 1. Department of Mathematics
More informationClassical Statistical Mechanics: Part 1
Classical Statistical Mechanics: Part 1 January 16, 2013 Classical Mechanics 1-Dimensional system with 1 particle of mass m Newton s equations of motion for position x(t) and momentum p(t): ẋ(t) dx p =
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 information4. The Green Kubo Relations
4. The Green Kubo Relations 4.1 The Langevin Equation In 1828 the botanist Robert Brown observed the motion of pollen grains suspended in a fluid. Although the system was allowed to come to equilibrium,
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 informationExtracting Macroscopic Stochastic Dynamics: Model Problems
Extracting Macroscopic Stochastic Dynamics: Model Problems WILHELM HUISINGA Free University CHRISTOF SCHÜTTE Free University AND ANDREW M. STUART University of Warwick Abstract The purpose of this work
More informationLecture 4: Stochastic Processes (I)
Miranda Holmes-Cerfon Applied Stochastic Analysis, Spring 215 Lecture 4: Stochastic Processes (I) Readings Recommended: Pavliotis [214], sections 1.1, 1.2 Grimmett and Stirzaker [21], section 9.4 (spectral
More informationFrom a Mesoscopic to a Macroscopic Description of Fluid-Particle Interaction
From a Mesoscopic to a Macroscopic Description of Fluid-Particle Interaction Carnegie Mellon University Center for Nonlinear Analysis Working Group, October 2016 Outline 1 Physical Framework 2 3 Free Energy
More informationLecture 6: Bayesian Inference in SDE Models
Lecture 6: Bayesian Inference in SDE Models Bayesian Filtering and Smoothing Point of View Simo Särkkä Aalto University Simo Särkkä (Aalto) Lecture 6: Bayesian Inference in SDEs 1 / 45 Contents 1 SDEs
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 informationSDE Coefficients. March 4, 2008
SDE Coefficients March 4, 2008 The following is a summary of GARD sections 3.3 and 6., mainly as an overview of the two main approaches to creating a SDE model. Stochastic Differential Equations (SDE)
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 information