model reduction in stochastic dynamics and applications to...
TRANSCRIPT
![Page 1: Model reduction in stochastic dynamics and applications to ...helper.ipam.ucla.edu/publications/elws4/elws4_14925.pdf · Our aim dX t = −∇V(X t)dt + p 2β−1 dW t in Rn Given](https://reader033.vdocument.in/reader033/viewer/2022060406/5f0f6f7e7e708231d444261f/html5/thumbnails/1.jpg)
Model reduction in stochastic dynamics
and applications to Molecular Dynamics
Frederic Legoll
Ecole des Ponts & project-team MATHERIALS, INRIA Paris
Joint with T. Lelievre (ENPC and INRIA), U. Sharma (ENPC) and S. Olla
(CEREMADE)
IPAM long program on Complex High-Dimensional Energy Landscapes
Workshop IV: UQ for Stochastic Systems and Applications, Nov. 13-17, 2017
Frederic Legoll (ENPC & INRIA) IPAM W4, 13-17 Nov. 2017 1 / 31
![Page 2: Model reduction in stochastic dynamics and applications to ...helper.ipam.ucla.edu/publications/elws4/elws4_14925.pdf · Our aim dX t = −∇V(X t)dt + p 2β−1 dW t in Rn Given](https://reader033.vdocument.in/reader033/viewer/2022060406/5f0f6f7e7e708231d444261f/html5/thumbnails/2.jpg)
Molecular systems
All the physics is encoded in the potential energy function,
V : Rn → R
Frederic Legoll (ENPC & INRIA) IPAM W4, 13-17 Nov. 2017 2 / 31
![Page 3: Model reduction in stochastic dynamics and applications to ...helper.ipam.ucla.edu/publications/elws4/elws4_14925.pdf · Our aim dX t = −∇V(X t)dt + p 2β−1 dW t in Rn Given](https://reader033.vdocument.in/reader033/viewer/2022060406/5f0f6f7e7e708231d444261f/html5/thumbnails/3.jpg)
Molecular simulation
Quantities of interest in molecular dynamics:
thermodynamical averages wrt Gibbs measure:
〈Φ〉 =ˆ
Rn
Φ(X ) dµGibbs, dµGibbs = Z−1 exp(−βV (X )) dX
or dynamical quantities (diffusion coefficients, rate constants, . . . ).
Consider the dynamics (overdamped Langevin equation)
dXt = −∇V (Xt) dt +√2β−1 dWt in R
n,
prototypical of those used in MD (e.g. ergodic for the Gibbs measure).
Frederic Legoll (ENPC & INRIA) IPAM W4, 13-17 Nov. 2017 3 / 31
![Page 4: Model reduction in stochastic dynamics and applications to ...helper.ipam.ucla.edu/publications/elws4/elws4_14925.pdf · Our aim dX t = −∇V(X t)dt + p 2β−1 dW t in Rn Given](https://reader033.vdocument.in/reader033/viewer/2022060406/5f0f6f7e7e708231d444261f/html5/thumbnails/4.jpg)
Molecular simulation
Quantities of interest in molecular dynamics:
thermodynamical averages wrt Gibbs measure:
〈Φ〉 =ˆ
Rn
Φ(X ) dµGibbs, dµGibbs = Z−1 exp(−βV (X )) dX
or dynamical quantities (diffusion coefficients, rate constants, . . . ).
Consider the dynamics (overdamped Langevin equation)
dXt = −∇V (Xt) dt +√2β−1 dWt in R
n,
prototypical of those used in MD (e.g. ergodic for the Gibbs measure).
In practice, quantities of interest often depend on a few variables.
Reduced description of the system, that still includes some dynamicalinformation?
Frederic Legoll (ENPC & INRIA) IPAM W4, 13-17 Nov. 2017 3 / 31
![Page 5: Model reduction in stochastic dynamics and applications to ...helper.ipam.ucla.edu/publications/elws4/elws4_14925.pdf · Our aim dX t = −∇V(X t)dt + p 2β−1 dW t in Rn Given](https://reader033.vdocument.in/reader033/viewer/2022060406/5f0f6f7e7e708231d444261f/html5/thumbnails/5.jpg)
Metastability and reaction coordinate
dXt = −∇V (Xt) dt +√
2β−1 dWt , Xt ≡ position of all atoms
in practice, the dynamics is metastable: the system stays a long timein a well of V before jumping to another well:
25002000150010005000
1
0
-1
we assume that wells are fully described through a well-chosenreaction coordinate ξ : Rn 7→ R
ξ(x) may be e.g. the coordinate of some particular atom.
Quantity of interest: path t 7→ ξ(Xt).
Frederic Legoll (ENPC & INRIA) IPAM W4, 13-17 Nov. 2017 4 / 31
![Page 6: Model reduction in stochastic dynamics and applications to ...helper.ipam.ucla.edu/publications/elws4/elws4_14925.pdf · Our aim dX t = −∇V(X t)dt + p 2β−1 dW t in Rn Given](https://reader033.vdocument.in/reader033/viewer/2022060406/5f0f6f7e7e708231d444261f/html5/thumbnails/6.jpg)
Our aim
dXt = −∇V (Xt) dt +√
2β−1 dWt in Rn
Given a reaction coordinate ξ : Rn 7→ R,propose a dynamics zt that approximates ξ(Xt).
preservation of equilibrium properties:when X ∼ dµGibbs, then ξ(X ) is distributed according toexp(−βA(z)) dz , where A is the free energy.
The dynamics zt should be ergodic wrt exp(−βA(z)) dz .
recover in zt some dynamical information included in ξ(Xt).
what about non-reversible cases:
dXt = F(Xt) dt +√
2β−1 σ(Xt) dWt , F not gradient
Related works: Mori-Zwanzig approaches, Schuette, Pavliotis and Stuart,Hartmann, Papanicolaou, E & Vanden-Eijnden, . . .Frederic Legoll (ENPC & INRIA) IPAM W4, 13-17 Nov. 2017 5 / 31
![Page 7: Model reduction in stochastic dynamics and applications to ...helper.ipam.ucla.edu/publications/elws4/elws4_14925.pdf · Our aim dX t = −∇V(X t)dt + p 2β−1 dW t in Rn Given](https://reader033.vdocument.in/reader033/viewer/2022060406/5f0f6f7e7e708231d444261f/html5/thumbnails/7.jpg)
Construction of
an effective dynamics
Frederic Legoll (ENPC & INRIA) IPAM W4, 13-17 Nov. 2017 6 / 31
![Page 8: Model reduction in stochastic dynamics and applications to ...helper.ipam.ucla.edu/publications/elws4/elws4_14925.pdf · Our aim dX t = −∇V(X t)dt + p 2β−1 dW t in Rn Given](https://reader033.vdocument.in/reader033/viewer/2022060406/5f0f6f7e7e708231d444261f/html5/thumbnails/8.jpg)
Effective dynamics using conditional expectations - 1
dXt = −∇V (Xt) dt +√
2β−1 dWt , ξ : Rn → R
From the dynamics on Xt , we obtain (chain rule)
d [ξ(Xt)] =(−∇V · ∇ξ + β−1∆ξ
)(Xt) dt +
√2β−1 |∇ξ|(Xt) dBt
where Bt is a 1D brownian motion.
Frederic Legoll (ENPC & INRIA) IPAM W4, 13-17 Nov. 2017 7 / 31
![Page 9: Model reduction in stochastic dynamics and applications to ...helper.ipam.ucla.edu/publications/elws4/elws4_14925.pdf · Our aim dX t = −∇V(X t)dt + p 2β−1 dW t in Rn Given](https://reader033.vdocument.in/reader033/viewer/2022060406/5f0f6f7e7e708231d444261f/html5/thumbnails/9.jpg)
Effective dynamics using conditional expectations - 1
dXt = −∇V (Xt) dt +√
2β−1 dWt , ξ : Rn → R
From the dynamics on Xt , we obtain (chain rule)
d [ξ(Xt)] =(−∇V · ∇ξ + β−1∆ξ
)(Xt) dt +
√2β−1 |∇ξ|(Xt) dBt
where Bt is a 1D brownian motion.
Introduce the average of the drift term:
b(z) := −ˆ (
−∇V · ∇ξ + β−1∆ξ)(X ) ψGibbs(X ) δξ(X )−z dX
= EGibbs
[(−∇V · ∇ξ + β−1∆ξ
)(X )
∣∣∣ ξ(X ) = z]
Frederic Legoll (ENPC & INRIA) IPAM W4, 13-17 Nov. 2017 7 / 31
![Page 10: Model reduction in stochastic dynamics and applications to ...helper.ipam.ucla.edu/publications/elws4/elws4_14925.pdf · Our aim dX t = −∇V(X t)dt + p 2β−1 dW t in Rn Given](https://reader033.vdocument.in/reader033/viewer/2022060406/5f0f6f7e7e708231d444261f/html5/thumbnails/10.jpg)
Effective dynamics using conditional expectations - 1
dXt = −∇V (Xt) dt +√
2β−1 dWt , ξ : Rn → R
From the dynamics on Xt , we obtain (chain rule)
d [ξ(Xt)] =(−∇V · ∇ξ + β−1∆ξ
)(Xt) dt +
√2β−1 |∇ξ|(Xt) dBt
where Bt is a 1D brownian motion.
Introduce the average of the drift term:
b(z) := −ˆ (
−∇V · ∇ξ + β−1∆ξ)(X ) ψGibbs(X ) δξ(X )−z dX
= EGibbs
[(−∇V · ∇ξ + β−1∆ξ
)(X )
∣∣∣ ξ(X ) = z]
and likewise for the diffusion term:
σ2(z) := −ˆ
|∇ξ(X )|2 ψGibbs(X ) δξ(X )−z dX
Frederic Legoll (ENPC & INRIA) IPAM W4, 13-17 Nov. 2017 7 / 31
![Page 11: Model reduction in stochastic dynamics and applications to ...helper.ipam.ucla.edu/publications/elws4/elws4_14925.pdf · Our aim dX t = −∇V(X t)dt + p 2β−1 dW t in Rn Given](https://reader033.vdocument.in/reader033/viewer/2022060406/5f0f6f7e7e708231d444261f/html5/thumbnails/11.jpg)
Effective dynamics using conditional expectations - 2
Chain rule:
d [ξ(Xt)] =(−∇V · ∇ξ + β−1∆ξ
)(Xt) dt +
√2β−1 |∇ξ|(Xt) dBt
Average of the drift and diffusion terms:
b(z) := −ˆ (
−∇V · ∇ξ + β−1∆ξ)(X ) ψGibbs(X ) δξ(X )−z dX
σ2(z) := −ˆ
|∇ξ(X )|2 ψGibbs(X ) δξ(X )−z dX
Eff. dyn. we propose: dzt = b(zt) dt +√
2β−1 σ(zt) dBt
Frederic Legoll (ENPC & INRIA) IPAM W4, 13-17 Nov. 2017 8 / 31
![Page 12: Model reduction in stochastic dynamics and applications to ...helper.ipam.ucla.edu/publications/elws4/elws4_14925.pdf · Our aim dX t = −∇V(X t)dt + p 2β−1 dW t in Rn Given](https://reader033.vdocument.in/reader033/viewer/2022060406/5f0f6f7e7e708231d444261f/html5/thumbnails/12.jpg)
Effective dynamics using conditional expectations - 2
Chain rule:
d [ξ(Xt)] =(−∇V · ∇ξ + β−1∆ξ
)(Xt) dt +
√2β−1 |∇ξ|(Xt) dBt
Average of the drift and diffusion terms:
b(z) := −ˆ (
−∇V · ∇ξ + β−1∆ξ)(X ) ψGibbs(X ) δξ(X )−z dX
σ2(z) := −ˆ
|∇ξ(X )|2 ψGibbs(X ) δξ(X )−z dX
Eff. dyn. we propose: dzt = b(zt) dt +√
2β−1 σ(zt) dBt
The approximation makes sense if, in the manifold
Σz = {X ∈ Rn, ξ(X ) = z} ,
Xt quickly reaches equilibrium: no metastability in Σz .
Frederic Legoll (ENPC & INRIA) IPAM W4, 13-17 Nov. 2017 8 / 31
![Page 13: Model reduction in stochastic dynamics and applications to ...helper.ipam.ucla.edu/publications/elws4/elws4_14925.pdf · Our aim dX t = −∇V(X t)dt + p 2β−1 dW t in Rn Given](https://reader033.vdocument.in/reader033/viewer/2022060406/5f0f6f7e7e708231d444261f/html5/thumbnails/13.jpg)
Remarks
dzt = b(zt) dt +√
2β−1 σ(zt) dBt
OK from the statistical viewpoint: the dynamics is ergodic wrtexp(−βA(z))dz .
Reformulation:
dzt =[−σ2(zt)A′(zt) + 2β−1σ′(zt)σ(zt)
]dt +
√2β−1 σ(zt) dBt
with A the free energy associated to ξ, and σ2(z) = 〈|∇ξ|2〉Σz.
Interesting particular case: ξ(X ) = X 1
Then b = −A′ and the effective dynamics reads
dzt = −A′(zt) dt +√
2β−1 dBt
Frederic Legoll (ENPC & INRIA) IPAM W4, 13-17 Nov. 2017 9 / 31
![Page 14: Model reduction in stochastic dynamics and applications to ...helper.ipam.ucla.edu/publications/elws4/elws4_14925.pdf · Our aim dX t = −∇V(X t)dt + p 2β−1 dW t in Rn Given](https://reader033.vdocument.in/reader033/viewer/2022060406/5f0f6f7e7e708231d444261f/html5/thumbnails/14.jpg)
An example without clear time-scale separation
���������������������
���������������������
������������������������
������������������������
���������������������
���������������������
������������������
������������������
���������������������
���������������������
����������������������������
����������������������������
����������������������������
����������������������������
������������������
������������������
������������������������
������������������������
���������������������
���������������������
����������������������������
����������������������������
���������������������
���������������������
������������������������
������������������������
������������������
������������������
����������������������������
����������������������������
���������������������
���������������������
������������������
������������������
���������������������
���������������������
���������������������
���������������������
������������������
������������������
������������������
������������������ ���
������������������
���������������������
������������������
������������������
������������������
������������������
����������������������������
����������������������������
������������������������
������������������������
solvent-solvent, solvent-monomer: truncated LJ on r = ‖xi − xj‖:VWCA(r) = 4ε
(σ12
r12− 2
σ6
r6
)if r ≤ σ, 0 otherwise (repulsive potential)
monomer-monomer: double well on r = ‖x1 − x2‖
Reaction coordinate: the distance between the two monomers
Frederic Legoll (ENPC & INRIA) IPAM W4, 13-17 Nov. 2017 10 / 31
![Page 15: Model reduction in stochastic dynamics and applications to ...helper.ipam.ucla.edu/publications/elws4/elws4_14925.pdf · Our aim dX t = −∇V(X t)dt + p 2β−1 dW t in Rn Given](https://reader033.vdocument.in/reader033/viewer/2022060406/5f0f6f7e7e708231d444261f/html5/thumbnails/15.jpg)
An example without clear time-scale separation
���������������������
���������������������
������������������������
������������������������
���������������������
���������������������
������������������
������������������
���������������������
���������������������
����������������������������
����������������������������
����������������������������
����������������������������
������������������
������������������
������������������������
������������������������
���������������������
���������������������
����������������������������
����������������������������
���������������������
���������������������
������������������������
������������������������
������������������
������������������
����������������������������
����������������������������
���������������������
���������������������
������������������
������������������
���������������������
���������������������
���������������������
���������������������
������������������
������������������
������������������
������������������ ���
������������������
���������������������
������������������
������������������
������������������
������������������
����������������������������
����������������������������
������������������������
������������������������
solvent-solvent, solvent-monomer: truncated LJ on r = ‖xi − xj‖:VWCA(r) = 4ε
(σ12
r12− 2
σ6
r6
)if r ≤ σ, 0 otherwise (repulsive potential)
monomer-monomer: double well on r = ‖x1 − x2‖
Reaction coordinate: the distance between the two monomers
β Reference Eff. dyn. Dyn. based on A
0.5 262 ± 6 245 ± 5 504 ± 11
0.25 1.81 ± 0.04 1.68 ± 0.04 3.47 ± 0.08
Frederic Legoll (ENPC & INRIA) IPAM W4, 13-17 Nov. 2017 10 / 31
![Page 16: Model reduction in stochastic dynamics and applications to ...helper.ipam.ucla.edu/publications/elws4/elws4_14925.pdf · Our aim dX t = −∇V(X t)dt + p 2β−1 dW t in Rn Given](https://reader033.vdocument.in/reader033/viewer/2022060406/5f0f6f7e7e708231d444261f/html5/thumbnails/16.jpg)
Effective drift and diffusion
-100
-80
-60
-40
-20
0
20
40
60
80
100
1 1.2 1.4 1.6 1.8 2 2.2
b(z)−A′(z)
z 1.98
1.985
1.99
1.995
2
2.005
2.01
2.015
2.02
1 1.2 1.4 1.6 1.8 2 2.2
σ2(z)
z
b(z) = EGibbs
[−∇V · ∇ξ + β−1∆ξ
∣∣∣ ξ(X ) = z]
−A′(z) = EGibbs
[− ∇V · ∇ξ
|∇ξ|2 + β−1div
( ∇ξ|∇ξ|2
) ∣∣∣ ξ(X ) = z]
σ2(z) = EGibbs
[|∇ξ|2
∣∣∣ ξ(X ) = z]
Frederic Legoll (ENPC & INRIA) IPAM W4, 13-17 Nov. 2017 11 / 31
![Page 17: Model reduction in stochastic dynamics and applications to ...helper.ipam.ucla.edu/publications/elws4/elws4_14925.pdf · Our aim dX t = −∇V(X t)dt + p 2β−1 dW t in Rn Given](https://reader033.vdocument.in/reader033/viewer/2022060406/5f0f6f7e7e708231d444261f/html5/thumbnails/17.jpg)
Accuracy of the effective dynamics:
the entropy approach
FL, T. Lelievre, Nonlinearity 2010FL, T. Lelievre, Springer LN Comput. Sci. Eng., vol. 82, 2012
Our effective dynamics is ergodic wrt exp(−βA(z))dzAt any time t, law of zt ≈ law of ξ(Xt)?
Frederic Legoll (ENPC & INRIA) IPAM W4, 13-17 Nov. 2017 12 / 31
![Page 18: Model reduction in stochastic dynamics and applications to ...helper.ipam.ucla.edu/publications/elws4/elws4_14925.pdf · Our aim dX t = −∇V(X t)dt + p 2β−1 dW t in Rn Given](https://reader033.vdocument.in/reader033/viewer/2022060406/5f0f6f7e7e708231d444261f/html5/thumbnails/18.jpg)
A convergence result
dXt = −∇V (Xt) dt +√
2β−1dWt , consider ξ(Xt)
Let ψexact(t, z) be the probability distribution function of ξ(Xt):
P (ξ(Xt) ∈ I ) =
ˆ
I
ψexact(t, z) dz
We have introduced the effective dynamics
dzt = b(zt) dt +√
2β−1 σ(zt) dBt
Let φeff(t, z) be the probability distribution function of zt .
Introduce the error
E (t) :=
ˆ
R
ψexact(t, ·) lnψexact(t, ·)φeff(t, ·)
We would like ψexact ≈ φeff , i.e. E small . . .Frederic Legoll (ENPC & INRIA) IPAM W4, 13-17 Nov. 2017 13 / 31
![Page 19: Model reduction in stochastic dynamics and applications to ...helper.ipam.ucla.edu/publications/elws4/elws4_14925.pdf · Our aim dX t = −∇V(X t)dt + p 2β−1 dW t in Rn Given](https://reader033.vdocument.in/reader033/viewer/2022060406/5f0f6f7e7e708231d444261f/html5/thumbnails/19.jpg)
Error estimate (FL, T. Lelievre, Nonlinearity 2010)
E (t) = error =
ˆ
R
ψexact(t, ·) lnψexact(t, ·)φeff(t, ·)
Under some assumptions that formalize the fact that
fast ergodicity in Σz (quantified by some ρ≫ 1),small coupling between dynamics in Σz and on zt (quantified by someκ≪ 1),
we have, for all t ≥ 0,
E (t) ≤ Cκ2
ρ2
The effective dynamics is accurate in the sense that
at any time t, law of ξ(Xt) ≈ law of zt .
Remark 1: this is not an asymptotic result, and this holds for any ξ.Remark 2: this estimate does not contain any information aboutcorrelations in time . . .Frederic Legoll (ENPC & INRIA) IPAM W4, 13-17 Nov. 2017 14 / 31
![Page 20: Model reduction in stochastic dynamics and applications to ...helper.ipam.ucla.edu/publications/elws4/elws4_14925.pdf · Our aim dX t = −∇V(X t)dt + p 2β−1 dW t in Rn Given](https://reader033.vdocument.in/reader033/viewer/2022060406/5f0f6f7e7e708231d444261f/html5/thumbnails/20.jpg)
Accuracy of the effective dynamics:
a pathwise approach
FL, T. Lelievre, S. Olla, Stoch. Processes and their Applications 127, 2017
Frederic Legoll (ENPC & INRIA) IPAM W4, 13-17 Nov. 2017 15 / 31
![Page 21: Model reduction in stochastic dynamics and applications to ...helper.ipam.ucla.edu/publications/elws4/elws4_14925.pdf · Our aim dX t = −∇V(X t)dt + p 2β−1 dW t in Rn Given](https://reader033.vdocument.in/reader033/viewer/2022060406/5f0f6f7e7e708231d444261f/html5/thumbnails/21.jpg)
Setting
For the sake of simplicity, we restrict ourselves to the case ξ(X ) = X 1.
The reference dynamics is
dXt = −∇V (Xt) dt +√
2β−1dWt
while the effective dynamics approximating t 7→ ξ(Xt) = X 1t is
dzt = b(zt) dt +√
2β−1 dW 1t
We aim at controlling E
[sup
0≤t≤T
|ξ(Xt)− zt |]= E
[sup
0≤t≤T
∣∣X 1t − zt
∣∣].
Frederic Legoll (ENPC & INRIA) IPAM W4, 13-17 Nov. 2017 16 / 31
![Page 22: Model reduction in stochastic dynamics and applications to ...helper.ipam.ucla.edu/publications/elws4/elws4_14925.pdf · Our aim dX t = −∇V(X t)dt + p 2β−1 dW t in Rn Given](https://reader033.vdocument.in/reader033/viewer/2022060406/5f0f6f7e7e708231d444261f/html5/thumbnails/22.jpg)
Exact and approximate dynamics
Non-closed dynamics:
dX 1t = −∇1V (Xt) dt +
√2β−1dBt , Bt = W 1
t .
Effective dynamics:dzt = b(zt) dt +
√2β−1dBt
with
−b(z) = EGibbs
(∂1V (X )
∣∣∣X 1 = z)=
ˆ
Rn−1
∂1V (z , xn2 ) ψzGibbs(x
n2 ) dx
n2
where
ψzGibbs(x
n2 ) =
ψGibbs(z , xn2 )
ˆ
Rn−1
ψGibbs(z , xn2 ) dx
n2
, xn2 = (x2, . . . , xn)
The above non-closed dynamics can be written
dX 1t = b(X 1
t ) + f (Xt) dt +√
2β−1dBt , mean of f vanishes
Frederic Legoll (ENPC & INRIA) IPAM W4, 13-17 Nov. 2017 17 / 31
![Page 23: Model reduction in stochastic dynamics and applications to ...helper.ipam.ucla.edu/publications/elws4/elws4_14925.pdf · Our aim dX t = −∇V(X t)dt + p 2β−1 dW t in Rn Given](https://reader033.vdocument.in/reader033/viewer/2022060406/5f0f6f7e7e708231d444261f/html5/thumbnails/23.jpg)
Assumptions – 1
For any z , the conditional probability measures ψzGibbs(x
n2 ) satisfy a
Poincare inequality for a constant ρ independent of z : for anyv ∈ H1(ψz
Gibbs),
ˆ
Rn−1
(v −ˆ
Rn−1
v ψzGibbs
)2
ψzGibbs ≤
1
ρ
ˆ
Rn−1
∣∣∣∇v∣∣∣2ψzGibbs
where ∇v = (∂2v , . . . , ∂nv).
Recall that a Poincare inequality holds on a probability measureexp(−βW (x)) dx under relatively mild assumption on W .
Frederic Legoll (ENPC & INRIA) IPAM W4, 13-17 Nov. 2017 18 / 31
![Page 24: Model reduction in stochastic dynamics and applications to ...helper.ipam.ucla.edu/publications/elws4/elws4_14925.pdf · Our aim dX t = −∇V(X t)dt + p 2β−1 dW t in Rn Given](https://reader033.vdocument.in/reader033/viewer/2022060406/5f0f6f7e7e708231d444261f/html5/thumbnails/24.jpg)
Assumptions – 2
The cross derivative ∇∂1V is in L2(ψGibbs):
κ2 :=
ˆ
Rn
∣∣∣∇∂1V (x)∣∣∣2ψGibbs(x)dx <∞.
Frederic Legoll (ENPC & INRIA) IPAM W4, 13-17 Nov. 2017 19 / 31
![Page 25: Model reduction in stochastic dynamics and applications to ...helper.ipam.ucla.edu/publications/elws4/elws4_14925.pdf · Our aim dX t = −∇V(X t)dt + p 2β−1 dW t in Rn Given](https://reader033.vdocument.in/reader033/viewer/2022060406/5f0f6f7e7e708231d444261f/html5/thumbnails/25.jpg)
Assumptions – 2
The cross derivative ∇∂1V is in L2(ψGibbs):
κ2 :=
ˆ
Rn
∣∣∣∇∂1V (x)∣∣∣2ψGibbs(x)dx <∞.
The effective drift b is one-sided Lipschitz on R: there exists Lb > 0 suchthat ∀x ∈ R, b′(x) ≤ Lb
This assumption is satisfied if b is Lipschitz on bounded domains anddecreasing at infinity, which corresponds to a case when the associatedfree energy is smooth and convex at infinity.
Frederic Legoll (ENPC & INRIA) IPAM W4, 13-17 Nov. 2017 19 / 31
![Page 26: Model reduction in stochastic dynamics and applications to ...helper.ipam.ucla.edu/publications/elws4/elws4_14925.pdf · Our aim dX t = −∇V(X t)dt + p 2β−1 dW t in Rn Given](https://reader033.vdocument.in/reader033/viewer/2022060406/5f0f6f7e7e708231d444261f/html5/thumbnails/26.jpg)
Assumptions – 2
The cross derivative ∇∂1V is in L2(ψGibbs):
κ2 :=
ˆ
Rn
∣∣∣∇∂1V (x)∣∣∣2ψGibbs(x)dx <∞.
The effective drift b is one-sided Lipschitz on R: there exists Lb > 0 suchthat ∀x ∈ R, b′(x) ≤ Lb
This assumption is satisfied if b is Lipschitz on bounded domains anddecreasing at infinity, which corresponds to a case when the associatedfree energy is smooth and convex at infinity.
For any x > 0, introduce α(x) = sups∈[−x ,x ] |b′(s)|. We assume that
E
[(α(∣∣X 1
∣∣))2]
<∞.
This assumption is satisfied e.g. if V has polynomial growth.Frederic Legoll (ENPC & INRIA) IPAM W4, 13-17 Nov. 2017 19 / 31
![Page 27: Model reduction in stochastic dynamics and applications to ...helper.ipam.ucla.edu/publications/elws4/elws4_14925.pdf · Our aim dX t = −∇V(X t)dt + p 2β−1 dW t in Rn Given](https://reader033.vdocument.in/reader033/viewer/2022060406/5f0f6f7e7e708231d444261f/html5/thumbnails/27.jpg)
Error estimate
We work under the above assumptions, and we assume that the systemstarts under equilibrium:
X0 ∼ ψGibbs. (⋆)
Consider (Xt)0≤t≤T solution to the reference dynamics and (zt)0≤t≤T
solution to the effective dynamics over a bounded time interval [0,T ].Then, there exists a constant C , independent of ρ and κ, such that
E
[sup
0≤t≤T
∣∣X 1t − zt
∣∣]≤ C
κ
ρ.
Remark: assumption (⋆) can be relaxed if assume
∥∥∥∥ψ0
ψGibbs
∥∥∥∥L∞(Rn)
<∞.
Frederic Legoll (ENPC & INRIA) IPAM W4, 13-17 Nov. 2017 20 / 31
![Page 28: Model reduction in stochastic dynamics and applications to ...helper.ipam.ucla.edu/publications/elws4/elws4_14925.pdf · Our aim dX t = −∇V(X t)dt + p 2β−1 dW t in Rn Given](https://reader033.vdocument.in/reader033/viewer/2022060406/5f0f6f7e7e708231d444261f/html5/thumbnails/28.jpg)
Ingredients of the proof
dX 1t = b(X 1
t ) + f (Xt) dt +√
2β−1dBt , mean of f vanishes
dzt = b(zt) dt +√
2β−1dBt
which implies that
X 1t − zt =
ˆ t
0
(b(X 1
s )− b(zs))ds +
ˆ t
0f (Xs) ds
Frederic Legoll (ENPC & INRIA) IPAM W4, 13-17 Nov. 2017 21 / 31
![Page 29: Model reduction in stochastic dynamics and applications to ...helper.ipam.ucla.edu/publications/elws4/elws4_14925.pdf · Our aim dX t = −∇V(X t)dt + p 2β−1 dW t in Rn Given](https://reader033.vdocument.in/reader033/viewer/2022060406/5f0f6f7e7e708231d444261f/html5/thumbnails/29.jpg)
Ingredients of the proof
dX 1t = b(X 1
t ) + f (Xt) dt +√
2β−1dBt , mean of f vanishes
dzt = b(zt) dt +√
2β−1dBt
which implies that
X 1t − zt =
ˆ t
0
(b(X 1
s )− b(zs))ds +
ˆ t
0f (Xs) ds
as the mean of f vanishes, we use the Poincare inequality assumptionto directly obtain a bound on E
[f 2(Xt)
];
Frederic Legoll (ENPC & INRIA) IPAM W4, 13-17 Nov. 2017 21 / 31
![Page 30: Model reduction in stochastic dynamics and applications to ...helper.ipam.ucla.edu/publications/elws4/elws4_14925.pdf · Our aim dX t = −∇V(X t)dt + p 2β−1 dW t in Rn Given](https://reader033.vdocument.in/reader033/viewer/2022060406/5f0f6f7e7e708231d444261f/html5/thumbnails/30.jpg)
Ingredients of the proof
dX 1t = b(X 1
t ) + f (Xt) dt +√
2β−1dBt , mean of f vanishes
dzt = b(zt) dt +√
2β−1dBt
which implies that
X 1t − zt =
ˆ t
0
(b(X 1
s )− b(zs))ds +
ˆ t
0f (Xs) ds
as the mean of f vanishes, we use the Poincare inequality assumptionto directly obtain a bound on E
[f 2(Xt)
];
through the introduction of a Poisson equation (with f as right-handside), and using an argument due to T. Lyons and T. Zhang, we get
an estimate on E
(sup
0≤t≤T
∣∣∣∣ˆ t
0f (Xs)ds
∣∣∣∣2);
Frederic Legoll (ENPC & INRIA) IPAM W4, 13-17 Nov. 2017 21 / 31
![Page 31: Model reduction in stochastic dynamics and applications to ...helper.ipam.ucla.edu/publications/elws4/elws4_14925.pdf · Our aim dX t = −∇V(X t)dt + p 2β−1 dW t in Rn Given](https://reader033.vdocument.in/reader033/viewer/2022060406/5f0f6f7e7e708231d444261f/html5/thumbnails/31.jpg)
Ingredients of the proof
dX 1t = b(X 1
t ) + f (Xt) dt +√
2β−1dBt , mean of f vanishes
dzt = b(zt) dt +√
2β−1dBt
which implies that
X 1t − zt =
ˆ t
0
(b(X 1
s )− b(zs))ds +
ˆ t
0f (Xs) ds
as the mean of f vanishes, we use the Poincare inequality assumptionto directly obtain a bound on E
[f 2(Xt)
];
through the introduction of a Poisson equation (with f as right-handside), and using an argument due to T. Lyons and T. Zhang, we get
an estimate on E
(sup
0≤t≤T
∣∣∣∣ˆ t
0f (Xs)ds
∣∣∣∣2);
Gronwall type argument to deduce a bound on E
(sup
0≤t≤T
∣∣X 1t − zt
∣∣).
Frederic Legoll (ENPC & INRIA) IPAM W4, 13-17 Nov. 2017 21 / 31
![Page 32: Model reduction in stochastic dynamics and applications to ...helper.ipam.ucla.edu/publications/elws4/elws4_14925.pdf · Our aim dX t = −∇V(X t)dt + p 2β−1 dW t in Rn Given](https://reader033.vdocument.in/reader033/viewer/2022060406/5f0f6f7e7e708231d444261f/html5/thumbnails/32.jpg)
Gronwall type argument
dXt = b(Xt)dt + σ dBt , dYt = b(Yt)dt + σ dBt + ftdt
The standard Gronwall Lemma gives an estimate of the difference(|Xt − Yt |)t∈[0,T ] in terms of (|ft |)t∈[0,T ].
We have an estimate on (|ft |)t∈[0,T ], but we have a much better one
on
(∣∣∣∣ˆ t
0fsds
∣∣∣∣)
t∈[0,T ]
that we wish to use.
This is why, in addition to the one-sided Lipschitz assumption on b,we assume that
Cα(β) = E
[(α(∣∣X 1
∣∣))2]
=
ˆ
R
(sup
s∈[−|x |,|x |]|b′(s)|
)2
ϕGibbs(x) dx <∞
The Gronwall argument is of interest by its own.
Frederic Legoll (ENPC & INRIA) IPAM W4, 13-17 Nov. 2017 22 / 31
![Page 33: Model reduction in stochastic dynamics and applications to ...helper.ipam.ucla.edu/publications/elws4/elws4_14925.pdf · Our aim dX t = −∇V(X t)dt + p 2β−1 dW t in Rn Given](https://reader033.vdocument.in/reader033/viewer/2022060406/5f0f6f7e7e708231d444261f/html5/thumbnails/33.jpg)
Conclusions on the reversible case
We proposed a “natural” way to obtain a closed equation on ξ(Xt)
Encouraging numerical results and rigorous error bounds
Our approach can also be applied to the standard problem
{dX ε,1
t = −∂1V (X εt ) dt +
√2β−1 dW 1
t
dX ε,it = −ε−1 ∂iV (X ε
t ) dt +√2β−1ε−1 dW i
t for i = 2, . . . , n
without Lipschitz assumptions on ∇V .
FL, T. Lelievre, Nonlinearity 23, 2010FL, T. Lelievre, Springer LN Comput. Sci. Eng., vol. 82, 2012FL, T. Lelievre, S. Olla, Stoch. Processes and their Applications 127, 2017
See also works by Duong, Lamacz, Pelletier, Schlichting and Sharma.
Frederic Legoll (ENPC & INRIA) IPAM W4, 13-17 Nov. 2017 23 / 31
![Page 34: Model reduction in stochastic dynamics and applications to ...helper.ipam.ucla.edu/publications/elws4/elws4_14925.pdf · Our aim dX t = −∇V(X t)dt + p 2β−1 dW t in Rn Given](https://reader033.vdocument.in/reader033/viewer/2022060406/5f0f6f7e7e708231d444261f/html5/thumbnails/34.jpg)
Extension to some
non-reversible cases
dXt = F(Xt) dt +√2 dWt , F is not a gradient
FL, T. Lelievre, U. Sharma, ongoing works
Frederic Legoll (ENPC & INRIA) IPAM W4, 13-17 Nov. 2017 24 / 31
![Page 35: Model reduction in stochastic dynamics and applications to ...helper.ipam.ucla.edu/publications/elws4/elws4_14925.pdf · Our aim dX t = −∇V(X t)dt + p 2β−1 dW t in Rn Given](https://reader033.vdocument.in/reader033/viewer/2022060406/5f0f6f7e7e708231d444261f/html5/thumbnails/35.jpg)
A motivating example – 1
dXt = F(Xt) dt +√2 dWt on T
2
for a Z2 periodic function F of the form
F(x , y) = (−∂yψ(x , y), ∂xψ(x , y)).
Unique stationary measure:
µ(x , y) = 1
Conditional stationary measure:
µx(y) =µ(x , y)
´
µ(x , y) dy= 1
Free energy associated to the reaction coordinate ξ(x , y) = x :in non-reversible cases, there is no obvious definitionone possible definition (reversible setting spirit):
exp(−A(x)) =
ˆ
µ(x , y) dy
and hence A(x) = 1 here.
Frederic Legoll (ENPC & INRIA) IPAM W4, 13-17 Nov. 2017 25 / 31
![Page 36: Model reduction in stochastic dynamics and applications to ...helper.ipam.ucla.edu/publications/elws4/elws4_14925.pdf · Our aim dX t = −∇V(X t)dt + p 2β−1 dW t in Rn Given](https://reader033.vdocument.in/reader033/viewer/2022060406/5f0f6f7e7e708231d444261f/html5/thumbnails/36.jpg)
A motivating example – 2
dXt = F(Xt) dt +√2 dWt on T
2, F = (−∂yψ, ∂xψ)and with ξ(X ) = ξ(x , y) = x , we have
dxt = −∂yψ(xt , yt) dt +√2 dW 1
t
Closure procedure:
dzt = b(zt) dt +√2 dW 1
t
with
b(z) = Eµ
[drift
∣∣∣ ξ(X ) = z]=
ˆ (− ∂yψ(z , y)
)µz(y) dy
In the simple case
ψ(X ) = ψper(X ) + L · X , L constant,
we get b(z) = −L2.
Frederic Legoll (ENPC & INRIA) IPAM W4, 13-17 Nov. 2017 26 / 31
![Page 37: Model reduction in stochastic dynamics and applications to ...helper.ipam.ucla.edu/publications/elws4/elws4_14925.pdf · Our aim dX t = −∇V(X t)dt + p 2β−1 dW t in Rn Given](https://reader033.vdocument.in/reader033/viewer/2022060406/5f0f6f7e7e708231d444261f/html5/thumbnails/37.jpg)
A motivating example – 3
dXt = F(Xt) dt +√2 dWt on T
2, F = (−∂yψ, ∂xψ)and
ξ(x , y) = x , ψ(X ) = ψper(X ) + L · X
Following our general procedure, we propose the effective dynamics
dzt = −L2 dt +√2 dW 1
t
On the other hand:
In the reversible case with ξ(x , y) = x , the effective dynamics is
dzt = −A′(zt) dt +√2 dW 1
t
A natural definition of A here leads to A′ = 0.
The two propositions differ! Which one is accurate?
Frederic Legoll (ENPC & INRIA) IPAM W4, 13-17 Nov. 2017 27 / 31
![Page 38: Model reduction in stochastic dynamics and applications to ...helper.ipam.ucla.edu/publications/elws4/elws4_14925.pdf · Our aim dX t = −∇V(X t)dt + p 2β−1 dW t in Rn Given](https://reader033.vdocument.in/reader033/viewer/2022060406/5f0f6f7e7e708231d444261f/html5/thumbnails/38.jpg)
Setting (non-reversible case)
For the sake of simplicity, we restrict ourselves to the case ξ(X ) = X 1.
Reference dynamics:
dXt = F(Xt) dt +√2dWt , unique stat. measure µ
Non-closed dynamics:
dX 1t = F1(Xt) dt +
√2dBt , Bt = W 1
t
Frederic Legoll (ENPC & INRIA) IPAM W4, 13-17 Nov. 2017 28 / 31
![Page 39: Model reduction in stochastic dynamics and applications to ...helper.ipam.ucla.edu/publications/elws4/elws4_14925.pdf · Our aim dX t = −∇V(X t)dt + p 2β−1 dW t in Rn Given](https://reader033.vdocument.in/reader033/viewer/2022060406/5f0f6f7e7e708231d444261f/html5/thumbnails/39.jpg)
Setting (non-reversible case)
For the sake of simplicity, we restrict ourselves to the case ξ(X ) = X 1.
Reference dynamics:
dXt = F(Xt) dt +√2dWt , unique stat. measure µ
Non-closed dynamics:
dX 1t = F1(Xt) dt +
√2dBt , Bt = W 1
t
Effective dynamics:
dzt = b(zt) dt +√2dBt
with
b(z) = Eµ
(F1(X )
∣∣∣X 1 = z)=
ˆ
Rn−1
F1(z , xn2 ) µ
z(xn2 ) dxn2
where
µz(xn2 ) =µ(z , xn2 )
ˆ
Rn−1
µ(z , xn2 ) dxn2
, xn2 = (x2, . . . , xn)
Frederic Legoll (ENPC & INRIA) IPAM W4, 13-17 Nov. 2017 28 / 31
![Page 40: Model reduction in stochastic dynamics and applications to ...helper.ipam.ucla.edu/publications/elws4/elws4_14925.pdf · Our aim dX t = −∇V(X t)dt + p 2β−1 dW t in Rn Given](https://reader033.vdocument.in/reader033/viewer/2022060406/5f0f6f7e7e708231d444261f/html5/thumbnails/40.jpg)
Assumptions (non-reversible case)
Similar assumptions as in the reversible case:
For any z , the conditional probability measures µz(xn2 ) satisfy aPoincare inequality for a constant ρ independent of z .
The cross derivative ∇F1 is in L2(µ):
κ2 :=
ˆ
Rn
∣∣∣∇F1
∣∣∣2µ <∞.
The effective drift b is one-sided Lipschitz on R.
For any x > 0, introduce α(x) = sups∈[−x ,x ] |b′(s)|. We assumethat
E
[(α(∣∣X 1
∣∣))2]
<∞.
Frederic Legoll (ENPC & INRIA) IPAM W4, 13-17 Nov. 2017 29 / 31
![Page 41: Model reduction in stochastic dynamics and applications to ...helper.ipam.ucla.edu/publications/elws4/elws4_14925.pdf · Our aim dX t = −∇V(X t)dt + p 2β−1 dW t in Rn Given](https://reader033.vdocument.in/reader033/viewer/2022060406/5f0f6f7e7e708231d444261f/html5/thumbnails/41.jpg)
Error estimate (non-reversible case)
We work under the above assumptions, and we assume that the systemstarts under equilibrium:
X0 ∼ µ ≡ stationary measure.
Consider (Xt)0≤t≤T solution to the reference dynamics and (zt)0≤t≤T
solution to the effective dynamics over a bounded time interval [0,T ].Then, there exists a constant C , independent of ρ and κ, such that
E
[sup
0≤t≤T
∣∣X 1t − zt
∣∣]≤ C
κ
ρ.
The proof essentially follows the same steps as in the reversible case.
Frederic Legoll (ENPC & INRIA) IPAM W4, 13-17 Nov. 2017 30 / 31
![Page 42: Model reduction in stochastic dynamics and applications to ...helper.ipam.ucla.edu/publications/elws4/elws4_14925.pdf · Our aim dX t = −∇V(X t)dt + p 2β−1 dW t in Rn Given](https://reader033.vdocument.in/reader033/viewer/2022060406/5f0f6f7e7e708231d444261f/html5/thumbnails/42.jpg)
On-going works and open questions
Ongoing work: general non-reversible SDE of the form
dXt = F(Xt) dt +√2σ(Xt) dWt
for a non-degenerate diffusion coefficient σ.
Open question: SDEs with degenerate noise, such as the Langevinequation:
dXt = Pt dt, dPt = F(Xt) dt − Pt dt +√2σ(Xt) dWt
https://team.inria.fr/matherials/
Frederic Legoll (ENPC & INRIA) IPAM W4, 13-17 Nov. 2017 31 / 31