Effective dynamics for the (overdamped) Langevin equation
|
|
- Julian Mosley
- 5 years ago
- Views:
Transcription
1 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 conference, Leicester, Sept 5-9, 2011 p. 1
2 Molecular simulation Some quantities of interest in molecular dynamics: thermodynamical averages wrt Gibbs measure: Φ = Φ(X) dµ, dµ = Z 1 exp( βv (X)) dx, X R n R n or dynamical quantities: diffusion coefficients rate constants residence times in metastable basins EnuMath conference, Leicester, Sept 5-9, 2011 p. 2
3 Molecular simulation Some quantities of interest in molecular dynamics: thermodynamical averages wrt Gibbs measure: Φ = Φ(X) dµ, dµ = Z 1 exp( βv (X)) dx, X R n R n or dynamical quantities: diffusion coefficients rate constants residence times in metastable basins In practice, quantities of interest often depend on a few variables. Reduced description of the system, that still includes some dynamical information? EnuMath conference, Leicester, Sept 5-9, 2011 p. 2
4 Example of a biological system Courtesy of Chris Chipot, CNRS Nancy EnuMath conference, Leicester, Sept 5-9, 2011 p. 3
5 Reference dynamics We are interested in dynamical properties. Two possible choices for the reference dynamics of the system: overdamped Langevin equation: dx t = V (X t ) dt + 2β 1 dw t, X t R n Langevin equation (with masses set to 1): dx t = P t dt, X t R n, dp t = V (X t ) dt γp t dt + 2γβ 1 dw t, P t R n. For both dynamics, 1 T T 0 Φ(X t ) dt R n Φ(X) dµ, dµ = Z 1 exp( βv (X)) dx. We will mostly argue with overdamped Langevin, and next turn to Langevin. EnuMath conference, Leicester, Sept 5-9, 2011 p. 4
6 Metastability and reaction coordinate dx t = V (X t ) dt + 2β 1 dw t, X t position of all atoms in practice, the dynamics is metastable: the system stays a long time in a well of V before jumping to another well: we assume that wells are fully described through a well-chosen reaction coordinate ξ : R n R ξ(x) may e.g. be a particular angle in the molecule. Quantity of interest: path t ξ(x t ). EnuMath conference, Leicester, Sept 5-9, 2011 p. 5
7 Our aim dx t = V (X t ) dt + 2β 1 dw t in R n Given a reaction coordinate ξ : R n R, propose a dynamics z t that approximates ξ(x t ). preservation of equilibrium properties: when X dµ, then ξ(x) is distributed according to exp( βa(z)) dz, where A is the free energy. The dynamics z t should be ergodic wrt exp( βa(z)) dz. recover in z t some dynamical information included in ξ(x t ). Related approaches: Mori-Zwanzig formalism, asymptotic expansion of the generator (Papanicolaou,... ), averaging principle for SDE (Pavliotis and Stuart, Hartmann,... ), effective dynamics using Markov state models (Schuette and Sarich, Lu),... EnuMath conference, Leicester, Sept 5-9, 2011 p. 6
8 A super-simple case: ξ(x,y) = x Consider the dynamics in two dimensions: X = (x, y) R 2, dx t = x V (x t, y t ) dt + 2β 1 dwt x, dy t = y V (x t, y t ) dt + 2β 1 dw y t, and assume that ξ(x, y) = x. EnuMath conference, Leicester, Sept 5-9, 2011 p. 7
9 A super-simple case: ξ(x,y) = x Consider the dynamics in two dimensions: X = (x, y) R 2, dx t = x V (x t, y t ) dt + 2β 1 dwt x, dy t = y V (x t, y t ) dt + 2β 1 dw y t, and assume that ξ(x, y) = x. Let ψ(t, X) be the density of X at time t: for any B R 2, P(X t B) = ψ(t, X)dX Introduce the mean of the drift over all configurations satisfying ξ(x) = z: x V (z, y) ψ(t, z, y) dy b(t, z) := R = E [ x V (X) ξ(x t ) = z] ψ(t, z, y) dy R and consider dz t = b(t, z t ) dt + 2β 1 db t Then, for any t, the law of z t is equal to the law of x t (Gyongy 1986) B EnuMath conference, Leicester, Sept 5-9, 2011 p. 7
10 Making the approach practical b(t, z) = R x V (z, y) ψ(t, z, y) dy = E [ x V (X) ξ(x t ) = z] b(t, z) is extremely difficult to compute... Need for approximation: EnuMath conference, Leicester, Sept 5-9, 2011 p. 8
11 Making the approach practical b(t, z) = R x V (z, y) ψ(t, z, y) dy = E [ x V (X) ξ(x t ) = z] b(t, z) is extremely difficult to compute... Need for approximation: b(z) := x V (z, y) ψ (z, y) dy = E µ [ x V (X) ξ(x) = z] R with ψ (x, y) = Z 1 exp( βv (x, y)). Effective dynamics: dz t = b(z t ) dt + 2β 1 db t Idea: b(t, x) b(x) if the equilibrium in each manifold Σ x = {(x, y), y R} is quickly reached: x t is much slower than y t. EnuMath conference, Leicester, Sept 5-9, 2011 p. 8
12 The general case: X R n and arbitrary ξ dx t = V (X t ) dt + 2β 1 dw t, ξ : R n R From the dynamics on X t, we obtain (chain rule) d [ξ(x t )] = ( V ξ + β 1 ξ ) (X t ) dt + 2β 1 ξ (X t ) db t where B t is a 1D brownian motion. EnuMath conference, Leicester, Sept 5-9, 2011 p. 9
13 The general case: X R n and arbitrary ξ dx t = V (X t ) dt + 2β 1 dw t, ξ : R n R From the dynamics on X t, we obtain (chain rule) d [ξ(x t )] = ( V ξ + β 1 ξ ) (X t ) dt + 2β 1 ξ (X t ) db t where B t is a 1D brownian motion. Introduce the average of the drift and diffusion terms: ( V b(z) := ξ + β 1 ξ ) (X) ψ (X) δ ξ(x) z dx σ 2 (z) := ξ(x) 2 ψ (X) δ ξ(x) z dx Eff. dyn.: dz t = b(z t ) dt + 2β 1 σ(z t ) db t The approximation makes sense if, in the manifold Σ z = {X R n, ξ(x) = z}, X t quickly reaches equilibrium. ξ(x t ) much slower than evolution of X t in Σ z. EnuMath conference, Leicester, Sept 5-9, 2011 p. 9
14 Some remarks Effective dynamics: dz t = b(z t ) dt + 2β 1 σ(z t ) db t OK from the statistical viewpoint: the dynamics is ergodic wrt exp( βa(z))dz. Using different arguments, this dynamics has been obtained by [E and Vanden-Eijnden 2004], and [Maragliano, Fischer, Vanden-Eijnden and Ciccotti, 2006]. In the following, we will numerically assess its accuracy derive error bounds In practice, we pre-compute b(z) and σ(z) for values on a grid (remember z is scalar), and next linearly interpolate between these values. EnuMath conference, Leicester, Sept 5-9, 2011 p. 10
15 Dimer in solution: comparison of residence times solvent-solvent, solvent-monomer: truncated LJ on r = x i x j : ( ) σ 12 σ6 V WCA (r) = 4ε 2 if r σ, 0 otherwise (repulsive potential) r12 r 6 monomer-monomer: double well on r = x 1 x 2 Reaction coordinate: the distance between the two monomers EnuMath conference, Leicester, Sept 5-9, 2011 p. 11
16 Dimer in solution: comparison of residence times solvent-solvent, solvent-monomer: truncated LJ on r = x i x j : ( ) σ 12 σ6 V WCA (r) = 4ε 2 if r σ, 0 otherwise (repulsive potential) r12 r 6 monomer-monomer: double well on r = x 1 x 2 Reaction coordinate: the distance between the two monomers β Reference dyn. Effective dyn ± ± ± ± 0.04 EnuMath conference, Leicester, Sept 5-9, 2011 p. 11
17 Accuracy assessment: some background materials EnuMath conference, Leicester, Sept 5-9, 2011 p. 12
18 Accuracy assessment: some background materials dx t = V (X t ) dt + 2β 1 dw t Let ψ(t, X) be the probability distribution function of X t : P (X t B) = ψ(t, X) dx B Under mild assumptions, ψ(t, X) converges to ψ (X) = Z 1 exp( βv (X)) exponentially fast: ψ(t, ) H (ψ(t, ) ψ ) := ψ(t, ) ln C exp( 2ρt) ψ Relative entropy is interesting because ψ(t, ) ψ 2 L 1 2H (ψ(t, ) ψ ). The larger ρ is, the faster the convergence to equilibrium. Remark: ρ is the Logarithmic Sobolev inequality constant of ψ. EnuMath conference, Leicester, Sept 5-9, 2011 p. 12
19 A convergence result dx t = V (X t ) dt + 2β 1 dw t, consider ξ(x t ) Let ψ exact (t, z) be the probability distribution function of ξ(x t ): P (ξ(x t ) I) = ψ exact (t, z) dz I On the other hand, we have introduced the effective dynamics dz t = b(z t ) dt + 2β 1 σ(z t ) db t Let φ eff (t, z) be the probability distribution function of z t. Introduce the error E(t) := R ψ exact (t, ) ln ψ exact(t, ) φ eff (t, ) We would like ψ exact φ eff, e.g. E small... EnuMath conference, Leicester, Sept 5-9, 2011 p. 13
20 Decoupling assumptions Σ z = {X R n, ξ(x) = z}, dµ z exp( βv (X)) δ ξ(x) z assume that the Gibbs measure restricted to Σ z satisfy a Logarithmic Sobolev inequality with a large constant ρ, uniform in z (no metastability in Σ z ). EnuMath conference, Leicester, Sept 5-9, 2011 p. 14
21 Decoupling assumptions Σ z = {X R n, ξ(x) = z}, dµ z exp( βv (X)) δ ξ(x) z assume that the Gibbs measure restricted to Σ z satisfy a Logarithmic Sobolev inequality with a large constant ρ, uniform in z (no metastability in Σ z ). assume that the coupling between the dynamics of ξ(x t ) and the dynamics in Σ z is weak: If ξ(x, y) = x, we request xy V to be small. In the general case, recall that the free energy derivative reads A (z) = We assume that max Σz F κ. Σ z F(X)dµ z EnuMath conference, Leicester, Sept 5-9, 2011 p. 14
22 Decoupling assumptions Σ z = {X R n, ξ(x) = z}, dµ z exp( βv (X)) δ ξ(x) z assume that the Gibbs measure restricted to Σ z satisfy a Logarithmic Sobolev inequality with a large constant ρ, uniform in z (no metastability in Σ z ). assume that the coupling between the dynamics of ξ(x t ) and the dynamics in Σ z is weak: If ξ(x, y) = x, we request xy V to be small. In the general case, recall that the free energy derivative reads A (z) = We assume that max Σz F κ. Σ z F(X)dµ z assume that ξ is close to a constant on each Σ z, e.g. λ = max ξ 2 (X) σ 2 (ξ(x)) X σ 2 (ξ(x)) is small EnuMath conference, Leicester, Sept 5-9, 2011 p. 14
23 Error estimate E(t) = error = R ψ exact (t, ) ln ψ exact(t, ) φ eff (t, ) Under the above assumptions, for all t 0, E(t) C(ξ, Initial Cond.) (λ + β2 κ 2 ) Hence, if the coarse variable ξ is such that ρ is large (fast ergodicity in Σ z ), κ is small (small coupling between dynamics in Σ z and on z t ), λ is small ( ξ is close to a constant on each Σ z ), then the effective dynamics is accurate: ρ 2 at any time, law of ξ(x t ) law of z t. Remark: this is not an asymptotic result, and this holds for any ξ. This bound may be helpful, in some cases, to discriminate between various reaction coordinates. EnuMath conference, Leicester, Sept 5-9, 2011 p. 15
24 Rough estimation in a particular case Standard expression in MD: V ε (X) = V 0 (X) + 1 ε q2 (X) : q fast direction E(t) = error = ψ exact (t, ) ln ψ exact(t, ) φ eff (t, ) R If ξ q = 0, then the direction ξ is decoupled from the fast direction q, hence ξ is indeed a slow variable, and it turns out that E(t) = O(ε). If ξ q 0, then the variable ξ does not contain all the slow motion, and bad scale separation: E(t) = O(1), hence the laws of ξ(x t ) and of z t are not close one to each other. The condition ξ q = 0 seems important to obtain good accuracy. EnuMath conference, Leicester, Sept 5-9, 2011 p. 16
25 Tri-atomic molecule B B A C A C V (X) = 1 2 k 2 (r AB l eq ) k 2 (r BC l eq ) 2 + k 3 W DW (θ ABC ), k 2 k 3, where W DW is a double well potential. Reaction coordinate Orthogonality Reference Residence time condition residence time using Eff. Dyn. ξ 1 (X) = θ ABC yes ± ± ξ 2 (X) = rac 2 no ± ± EnuMath conference, Leicester, Sept 5-9, 2011 p. 17
26 Extension to the Langevin equation We have considered until now the overdamped Langevin equation: dx t = V (X t ) dt + 2β 1 dw t Consider now the Langevin equation: dx t = P t dt dp t = V (X t ) dt γp t dt + 2γβ 1 dw t Can we extend the numerical strategy? Joint work with F. Galante and T. Lelièvre. EnuMath conference, Leicester, Sept 5-9, 2011 p. 18
27 Building the effective dynamics - 1 We compute dx t = P t dt, dp t = V (X t ) dt γp t dt + 2γβ 1 dw t d [ξ(x t )] = ξ(x t ) P t dt We introduce the coarse-grained velocity and have (chain rule) v(x, P) = ξ(x) P R d [v(x t, P t )] = [ Pt T 2 ξ(x t )P t ξ(x t ) T V (X t ) ] dt γv(x t, P t ) dt + 2γβ 1 ξ(x t ) db t where B t is a 1D Brownian motion. We wish to write a closed equation on ξ t = ξ(x t ) and v t = v(x t, P t ). Introduce D(X, P) = P T 2 ξ(x)p ξ(x) T V (X). EnuMath conference, Leicester, Sept 5-9, 2011 p. 19
28 Building the effective dynamics - 2 Without any approximation, we have obtained dξ t = v t dt, dv t = D(X t, P t ) dt γv t dt + 2γβ 1 ξ(x t ) db t To close the system, we introduce the conditional expectations with respect to the equilibrium measure µ(x, P) = Z 1 exp[ β ( V (X) + P T P/2 ) ]: Effective dynamics: D eff (ξ 0, v 0 ) = E µ ( D(X, P) ξ(x) = ξ0, v(x, P) = v 0 ) σ 2 (ξ 0, v 0 ) = E µ ( ξ 2 (X) ξ(x) = ξ 0, v(x, P) = v 0 ) dξ t = v t dt, dv t = D eff (ξ t, v t ) dt γv t dt + 2γβ 1 σ(ξ t, v t ) db t Again, this dynamics is consistent with the equilibrium properties: it preserves the equilibrium measure exp[ βa(ξ 0, v 0 )] dξ 0 dv 0. EnuMath conference, Leicester, Sept 5-9, 2011 p. 20
29 Numerical results: the tri-atomic molecule B B A C A C V (X) = 1 2 k 2 (r AB l eq ) k 2 (r BC l eq ) 2 + k 3 W DW (θ ABC ), k 2 k 3, where W DW is a double well potential. Reaction coordinate: ξ(x) = θ ABC. Inverse temp. Reference Eff. dyn. β = ± ± β = ± ± 1.25 Excellent agreement on the residence times in the well. EnuMath conference, Leicester, Sept 5-9, 2011 p. 21
30 Pathwise accuracy (ongoing work with T. Lelièvre and S. Olla) Going back to the overdamped case, with ξ(x, y) = x, can we get pathwise accuracy, e.g. [ ] E x t z t 2 C(T)? ρ sup 0 t T EnuMath conference, Leicester, Sept 5-9, 2011 p. 22
31 Pathwise accuracy (ongoing work with T. Lelièvre and S. Olla) Going back to the overdamped case, with ξ(x, y) = x, can we get pathwise accuracy, e.g. [ ] E x t z t 2 C(T)? ρ sup 0 t T z t is the effective dynamics trajectory: dz t = b(z t ) dt + 2β 1 db t (x t, y t ) is the exact trajectory: dx t = x V (x t, y t ) dt + 2β 1 db t = b(x t ) dt + e(x t, y t ) dt + 2β 1 db t By construction, b(x) = E µ [ x V (X) ξ(x) = x], hence x, E µ [e(x) ξ(x) = x] = 0. EnuMath conference, Leicester, Sept 5-9, 2011 p. 22
32 Poisson equation Let L x be the Fokker-Planck operator corresponding to a simple diffusion in y, at fixed x: For any x, L x φ = div y (φ y V ) + β 1 y φ E µ [e(x) ξ(x) = x] = 0 = u(x, ) s.t. (L x ) u(x, ) = e(x, ) Assume that L x satisfies a (large) spectral gap ρ 1. Then, u C/ρ 1. We have not been able to show a bound on E [ sup 0 t T x t z t 2], but we have shown that ( P sup 0 t T x t z t cρ α ) C ln(ρ) for any 0 α < 1/2. This somewhat explains the good results we observe on the residence times. EnuMath conference, Leicester, Sept 5-9, 2011 p. 23
33 Numerical illustration: the tri-atomic molecule B A C V (X) = 1 2 k 2 (r AB l eq ) k 2 (r BC l eq ) 2 + k 3 W(θ ABC ), k 2 k θ ABC (X t ) z t t EnuMath conference, Leicester, Sept 5-9, 2011 p. 24
34 Conclusions We have proposed a natural way to obtain a closed equation on ξ(x t ). Encouraging numerical results and rigorous error bounds (marginals at time t and in probability on the paths). Once a reaction coordinate ξ has been chosen, computing the drift and diffusion functions b(z) and σ(z) is as easy/difficult as computing the free energy derivative A (z). The approach can be extended to the Langevin equation. FL, T. Lelièvre, Nonlinearity 23, FL, T. Lelièvre, Springer LN Comput. Sci. Eng., vol. 82, 2011, arxiv FL, T. Lelièvre, S. Olla, in preparation F. Galante, FL, T. Lelièvre, in preparation EnuMath conference, Leicester, Sept 5-9, 2011 p. 25
Energie libre et dynamique réduite en dynamique moléculaire
Energie libre et dynamique réduite en dynamique moléculaire Frédéric Legoll Ecole des Ponts ParisTech et INRIA http://cermics.enpc.fr/ legoll Travaux en collaboration avec T. Lelièvre (Ecole des Ponts
More informationNumerical methods in molecular dynamics and multiscale problems
Numerical methods in molecular dynamics and multiscale problems Two examples T. Lelièvre CERMICS - Ecole des Ponts ParisTech & MicMac project-team - INRIA Horizon Maths December 2012 Introduction The aim
More informationTensorized Adaptive Biasing Force method for molecular dynamics
Tensorized Adaptive Biasing Force method for molecular dynamics Virginie Ehrlacher 1, Tony Lelièvre 1 & Pierre Monmarché 2 1 CERMICS, Ecole des Ponts ParisTech & MATHERIALS project, INRIA 2 Laboratoire
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 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 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 informationSDEs in large dimension and numerical methods Part 1: Sampling the canonical distribution
SDEs in large dimension and numerical methods Part 1: Sampling the canonical distribution T. Lelièvre CERMICS - Ecole des Ponts ParisTech & Matherials project-team - INRIA RICAM Winterschool, December
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 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 informationStochastic differential equations in neuroscience
Stochastic differential equations in neuroscience Nils Berglund MAPMO, Orléans (CNRS, UMR 6628) http://www.univ-orleans.fr/mapmo/membres/berglund/ Barbara Gentz, Universität Bielefeld Damien Landon, MAPMO-Orléans
More informationEntropic structure of the Landau equation. Coulomb interaction
with Coulomb interaction Laurent Desvillettes IMJ-PRG, Université Paris Diderot May 15, 2017 Use of the entropy principle for specific equations Spatially Homogeneous Kinetic equations: 1 Fokker-Planck:
More informationA mathematical framework for Exact Milestoning
A mathematical framework for Exact Milestoning David Aristoff (joint work with Juan M. Bello-Rivas and Ron Elber) Colorado State University July 2015 D. Aristoff (Colorado State University) July 2015 1
More informationA Lévy-Fokker-Planck equation: entropies and convergence to equilibrium
1/ 22 A Lévy-Fokker-Planck equation: entropies and convergence to equilibrium I. Gentil CEREMADE, Université Paris-Dauphine International Conference on stochastic Analysis and Applications Hammamet, Tunisia,
More informationMathematical analysis of an Adaptive Biasing Potential method for diffusions
Mathematical analysis of an Adaptive Biasing Potential method for diffusions Charles-Edouard Bréhier Joint work with Michel Benaïm (Neuchâtel, Switzerland) CNRS & Université Lyon 1, Institut Camille Jordan
More informationOn some relations between Optimal Transport and Stochastic Geometric Mechanics
Title On some relations between Optimal Transport and Stochastic Geometric Mechanics Banff, December 218 Ana Bela Cruzeiro Dep. Mathematics IST and Grupo de Física-Matemática Univ. Lisboa 1 / 2 Title Based
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 informationMonte Carlo methods in molecular dynamics.
T. Lelièvre, Rencontres EDP/Proba, Octobre 2009 p. 1 Monte Carlo methods in molecular dynamics. Tony Lelièvre CERMICS, Ecole Nationale des Ponts et Chaussées, France, MicMac team, INRIA, France. http://cermics.enpc.fr/
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 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 informationDIFFUSION MAPS, REDUCTION COORDINATES AND LOW DIMENSIONAL REPRESENTATION OF STOCHASTIC SYSTEMS
DIFFUSION MAPS, REDUCTION COORDINATES AND LOW DIMENSIONAL REPRESENTATION OF STOCHASTIC SYSTEMS R.R. COIFMAN, I.G. KEVREKIDIS, S. LAFON, M. MAGGIONI, AND B. NADLER Abstract. The concise representation of
More informationLecture No 1 Introduction to Diffusion equations The heat equat
Lecture No 1 Introduction to Diffusion equations The heat equation Columbia University IAS summer program June, 2009 Outline of the lectures We will discuss some basic models of diffusion equations and
More informationAsymptotic properties of maximum likelihood estimator for the growth rate for a jump-type CIR process
Asymptotic properties of maximum likelihood estimator for the growth rate for a jump-type CIR process Mátyás Barczy University of Debrecen, Hungary The 3 rd Workshop on Branching Processes and Related
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 informationLarge deviations and averaging for systems of slow fast stochastic reaction diffusion equations.
Large deviations and averaging for systems of slow fast stochastic reaction diffusion equations. Wenqing Hu. 1 (Joint work with Michael Salins 2, Konstantinos Spiliopoulos 3.) 1. Department of Mathematics
More informationA mathematical introduction to coarse-grained stochastic dynamics
A mathematical introduction to coarse-grained stochastic dynamics Gabriel STOLTZ gabriel.stoltz@enpc.fr (CERMICS, Ecole des Ponts & MATHERIALS team, INRIA Paris) Work also supported by ANR Funding ANR-14-CE23-0012
More informationOn rough PDEs. Samy Tindel. University of Nancy. Rough Paths Analysis and Related Topics - Nagoya 2012
On rough PDEs Samy Tindel University of Nancy Rough Paths Analysis and Related Topics - Nagoya 2012 Joint work with: Aurélien Deya and Massimiliano Gubinelli Samy T. (Nancy) On rough PDEs Nagoya 2012 1
More informationMalliavin Calculus in Finance
Malliavin Calculus in Finance Peter K. Friz 1 Greeks and the logarithmic derivative trick Model an underlying assent by a Markov process with values in R m with dynamics described by the SDE dx t = b(x
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 informationShort-time expansions for close-to-the-money options under a Lévy jump model with stochastic volatility
Short-time expansions for close-to-the-money options under a Lévy jump model with stochastic volatility José Enrique Figueroa-López 1 1 Department of Statistics Purdue University Statistics, Jump Processes,
More informationNonparametric Drift Estimation for Stochastic Differential Equations
Nonparametric Drift Estimation for Stochastic Differential Equations Gareth Roberts 1 Department of Statistics University of Warwick Brazilian Bayesian meeting, March 2010 Joint work with O. Papaspiliopoulos,
More informationOn different notions of timescales in molecular dynamics
On different notions of timescales in molecular dynamics Origin of Scaling Cascades in Protein Dynamics June 8, 217 IHP Trimester Stochastic Dynamics out of Equilibrium Overview 1. Motivation 2. Definition
More informationImportance splitting for rare event simulation
Importance splitting for rare event simulation F. Cerou Frederic.Cerou@inria.fr Inria Rennes Bretagne Atlantique Simulation of hybrid dynamical systems and applications to molecular dynamics September
More informationin Bounded Domains Ariane Trescases CMLA, ENS Cachan
CMLA, ENS Cachan Joint work with Yan GUO, Chanwoo KIM and Daniela TONON International Conference on Nonlinear Analysis: Boundary Phenomena for Evolutionnary PDE Academia Sinica December 21, 214 Outline
More information2012 NCTS Workshop on Dynamical Systems
Barbara Gentz gentz@math.uni-bielefeld.de http://www.math.uni-bielefeld.de/ gentz 2012 NCTS Workshop on Dynamical Systems National Center for Theoretical Sciences, National Tsing-Hua University Hsinchu,
More informationConvergence to equilibrium for rough differential equations
Convergence to equilibrium for rough differential equations Samy Tindel Purdue University Barcelona GSE Summer Forum 2017 Joint work with Aurélien Deya (Nancy) and Fabien Panloup (Angers) Samy T. (Purdue)
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 informationSMAI-JCM SMAI Journal of Computational Mathematics
SMAI-JCM SMAI Journal of Computational Mathematics Long-time convergence of an adaptive biasing force method: Variance reduction by Helmholtz projection Houssam Alrachid & Tony Lelièvre Volume (05, p.
More informationNONLOCAL DIFFUSION EQUATIONS
NONLOCAL DIFFUSION EQUATIONS JULIO D. ROSSI (ALICANTE, SPAIN AND BUENOS AIRES, ARGENTINA) jrossi@dm.uba.ar http://mate.dm.uba.ar/ jrossi 2011 Non-local diffusion. The function J. Let J : R N R, nonnegative,
More informationEntropy-dissipation methods I: Fokker-Planck equations
1 Entropy-dissipation methods I: Fokker-Planck equations Ansgar Jüngel Vienna University of Technology, Austria www.jungel.at.vu Introduction Boltzmann equation Fokker-Planck equations Degenerate parabolic
More informationFirst-Order ODE: Separable Equations, Exact Equations and Integrating Factor
First-Order ODE: Separable Equations, Exact Equations and Integrating Factor Department of Mathematics IIT Guwahati REMARK: In the last theorem of the previous lecture, you can change the open interval
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 informationSPDEs, criticality, and renormalisation
SPDEs, criticality, and renormalisation Hendrik Weber Mathematics Institute University of Warwick Potsdam, 06.11.2013 An interesting model from Physics I Ising model Spin configurations: Energy: Inverse
More informationLogFeller et Ray Knight
LogFeller et Ray Knight Etienne Pardoux joint work with V. Le and A. Wakolbinger Etienne Pardoux (Marseille) MANEGE, 18/1/1 1 / 16 Feller s branching diffusion with logistic growth We consider the diffusion
More informationA stochastic formulation of a dynamical singly constrained spatial interaction model
A stochastic formulation of a dynamical singly constrained spatial interaction model Mark Girolami Department of Mathematics, Imperial College London The Alan Turing Institute, British Library Lloyds Register
More informationStochastic (intermittent) Spikes and Strong Noise Limit of SDEs.
Stochastic (intermittent) Spikes and Strong Noise Limit of SDEs. D. Bernard in collaboration with M. Bauer and (in part) A. Tilloy. IPAM-UCLA, Jan 2017. 1 Strong Noise Limit of (some) SDEs Stochastic differential
More informationSome Basic Statistical Modeling Issues in Molecular and Ocean Dynamics
Some Basic Statistical Modeling Issues in Molecular and Ocean Dynamics p. 1/2 Some Basic Statistical Modeling Issues in Molecular and Ocean Dynamics Peter R. Kramer Department of Mathematical Sciences
More informationMultiple Random Variables
Multiple Random Variables Joint Probability Density Let X and Y be two random variables. Their joint distribution function is F ( XY x, y) P X x Y y. F XY ( ) 1, < x
More informationA quantum heat equation 5th Spring School on Evolution Equations, TU Berlin
A quantum heat equation 5th Spring School on Evolution Equations, TU Berlin Mario Bukal A. Jüngel and D. Matthes ACROSS - Centre for Advanced Cooperative Systems Faculty of Electrical Engineering and Computing
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 informationFinite temperature analysis of stochastic networks
Finite temperature analysis of stochastic networks Maria Cameron University of Maryland USA Stochastic networks with detailed balance * The generator matrix L = ( L ij ), i, j S L ij 0, i, j S, i j j S
More informationTheory and applications of random Poincaré maps
Nils Berglund nils.berglund@univ-orleans.fr http://www.univ-orleans.fr/mapmo/membres/berglund/ Random Dynamical Systems Theory and applications of random Poincaré maps Nils Berglund MAPMO, Université d'orléans,
More information3.1 Derivative Formulas for Powers and Polynomials
3.1 Derivative Formulas for Powers and Polynomials First, recall that a derivative is a function. We worked very hard in 2.2 to interpret the derivative of a function visually. We made the link, in Ex.
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 informationNumerical explorations of a forward-backward diffusion equation
Numerical explorations of a forward-backward diffusion equation Newton Institute KIT Programme Pauline Lafitte 1 C. Mascia 2 1 SIMPAF - INRIA & U. Lille 1, France 2 Univ. La Sapienza, Roma, Italy September
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 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 informationHeight fluctuations for the stationary KPZ equation
Firenze, June 22-26, 2015 Height fluctuations for the stationary KPZ equation P.L. Ferrari with A. Borodin, I. Corwin and B. Vető arxiv:1407.6977; To appear in MPAG http://wt.iam.uni-bonn.de/ ferrari Introduction
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 informationAnomalous transport of particles in Plasma physics
Anomalous transport of particles in Plasma physics L. Cesbron a, A. Mellet b,1, K. Trivisa b, a École Normale Supérieure de Cachan Campus de Ker Lann 35170 Bruz rance. b Department of Mathematics, University
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 informationThe dynamics of small particles whose size is roughly 1 µmt or. smaller, in a fluid at room temperature, is extremely erratic, and is
1 I. BROWNIAN MOTION The dynamics of small particles whose size is roughly 1 µmt or smaller, in a fluid at room temperature, is extremely erratic, and is called Brownian motion. The velocity of such particles
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 informationCapacitary inequalities in discrete setting and application to metastable Markov chains. André Schlichting
Capacitary inequalities in discrete setting and application to metastable Markov chains André Schlichting Institute for Applied Mathematics, University of Bonn joint work with M. Slowik (TU Berlin) 12
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 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 informationA MODEL FOR THE LONG-TERM OPTIMAL CAPACITY LEVEL OF AN INVESTMENT PROJECT
A MODEL FOR HE LONG-ERM OPIMAL CAPACIY LEVEL OF AN INVESMEN PROJEC ARNE LØKKA AND MIHAIL ZERVOS Abstract. We consider an investment project that produces a single commodity. he project s operation yields
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 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 informationmacroscopic view (phenomenological) microscopic view (atomistic) computing reaction rate rate of reactions experiments thermodynamics
Rate Theory (overview) macroscopic view (phenomenological) rate of reactions experiments thermodynamics Van t Hoff & Arrhenius equation microscopic view (atomistic) statistical mechanics transition state
More informationBridging the Gap between Center and Tail for Multiscale Processes
Bridging the Gap between Center and Tail for Multiscale Processes Matthew R. Morse Department of Mathematics and Statistics Boston University BU-Keio 2016, August 16 Matthew R. Morse (BU) Moderate Deviations
More informationGeometric projection of stochastic differential equations
Geometric projection of stochastic differential equations John Armstrong (King s College London) Damiano Brigo (Imperial) August 9, 2018 Idea: Projection Idea: Projection Projection gives a method of systematically
More informationFree energy estimates for the two-dimensional Keller-Segel model
Free energy estimates for the two-dimensional Keller-Segel model dolbeaul@ceremade.dauphine.fr CEREMADE CNRS & Université Paris-Dauphine in collaboration with A. Blanchet (CERMICS, ENPC & Ceremade) & B.
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 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 informationHomogenization of micro-resonances and localization of waves.
Homogenization of micro-resonances and localization of waves. Valery Smyshlyaev University College London, UK July 13, 2012 (joint work with Ilia Kamotski UCL, and Shane Cooper Bath/ Cardiff) Valery Smyshlyaev
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 informationThe Central Limit Theorem: More of the Story
The Central Limit Theorem: More of the Story Steven Janke November 2015 Steven Janke (Seminar) The Central Limit Theorem:More of the Story November 2015 1 / 33 Central Limit Theorem Theorem (Central Limit
More informationGlobal well-posedness of the primitive equations of oceanic and atmospheric dynamics
Global well-posedness of the primitive equations of oceanic and atmospheric dynamics Jinkai Li Department of Mathematics The Chinese University of Hong Kong Dynamics of Small Scales in Fluids ICERM, Feb
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 informationFokker-Planck Equation with Detailed Balance
Appendix E Fokker-Planck Equation with Detailed Balance A stochastic process is simply a function of two variables, one is the time, the other is a stochastic variable X, defined by specifying: a: the
More informationEstimation of arrival and service rates for M/M/c queue system
Estimation of arrival and service rates for M/M/c queue system Katarína Starinská starinskak@gmail.com Charles University Faculty of Mathematics and Physics Department of Probability and Mathematical Statistics
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 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 informationMetastability for the Ginzburg Landau equation with space time white noise
Barbara Gentz gentz@math.uni-bielefeld.de http://www.math.uni-bielefeld.de/ gentz Metastability for the Ginzburg Landau equation with space time white noise Barbara Gentz University of Bielefeld, Germany
More informationEDP with strong anisotropy : transport, heat, waves equations
EDP with strong anisotropy : transport, heat, waves equations Mihaï BOSTAN University of Aix-Marseille, FRANCE mihai.bostan@univ-amu.fr Nachos team INRIA Sophia Antipolis, 3/07/2017 Main goals Effective
More informationMean field games via probability manifold II
Mean field games via probability manifold II Wuchen Li IPAM summer school, 2018. Introduction Consider a nonlinear Schrödinger equation hi t Ψ = h2 1 2 Ψ + ΨV (x) + Ψ R d W (x, y) Ψ(y) 2 dy. The unknown
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 informationGravitational allocation to Poisson points
Gravitational allocation to Poisson points Sourav Chatterjee joint work with Ron Peled Yuval Peres Dan Romik Allocation rules Let Ξ be a discrete subset of R d. An allocation (of Lebesgue measure to Ξ)
More informationReaction Dynamics (2) Can we predict the rate of reactions?
Reaction Dynamics (2) Can we predict the rate of reactions? Reactions in Liquid Solutions Solvent is NOT a reactant Reactive encounters in solution A reaction occurs if 1. The reactant molecules (A, B)
More informationQuantum Hydrodynamics models derived from the entropy principle
1 Quantum Hydrodynamics models derived from the entropy principle P. Degond (1), Ch. Ringhofer (2) (1) MIP, CNRS and Université Paul Sabatier, 118 route de Narbonne, 31062 Toulouse cedex, France degond@mip.ups-tlse.fr
More informationLAW OF LARGE NUMBERS FOR THE SIRS EPIDEMIC
LAW OF LARGE NUMBERS FOR THE SIRS EPIDEMIC R. G. DOLGOARSHINNYKH Abstract. We establish law of large numbers for SIRS stochastic epidemic processes: as the population size increases the paths of SIRS epidemic
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 nonlinear Schrödinger equations and modulation of solitary waves
Stochastic nonlinear Schrödinger equations and modulation of solitary waves A. de Bouard CMAP, Ecole Polytechnique, France joint work with R. Fukuizumi (Sendai, Japan) Deterministic and stochastic front
More informationSome results on interspike interval statistics in conductance-based models for neuron action potentials
Some results on interspike interval statistics in conductance-based models for neuron action potentials Nils Berglund MAPMO, Université d Orléans CNRS, UMR 7349 et Fédération Denis Poisson www.univ-orleans.fr/mapmo/membres/berglund
More informationThermalization of Random Motion in Weakly Confining Potentials
Open Systems & Information Dynamics Vol. 17, No. 3 (2010) 1 10 DOI: xxx c World Scientific Publishing Company Thermalization of Random Motion in Weakly Confining Potentials Piotr Garbaczewski and Vladimir
More informationLAN property for ergodic jump-diffusion processes with discrete observations
LAN property for ergodic jump-diffusion processes with discrete observations Eulalia Nualart (Universitat Pompeu Fabra, Barcelona) joint work with Arturo Kohatsu-Higa (Ritsumeikan University, Japan) &
More informationTRANSPORT IN POROUS MEDIA
1 TRANSPORT IN POROUS MEDIA G. ALLAIRE CMAP, Ecole Polytechnique 1. Introduction 2. Main result in an unbounded domain 3. Asymptotic expansions with drift 4. Two-scale convergence with drift 5. The case
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 informationIrreversible Langevin samplers and variance reduction: a large deviations approach
Irreversible Langevin samplers and variance reduction: a large deviations approach Konstantinos Spiliopoulos Department of Mathematics & Statistics, Boston University Partially supported by NSF-DMS 1312124
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 information