A mathematical introduction to coarse-grained stochastic dynamics
|
|
- Maximilian Page
- 5 years ago
- Views:
Transcription
1 A mathematical introduction to coarse-grained stochastic dynamics Gabriel STOLTZ (CERMICS, Ecole des Ponts & MATHERIALS team, INRIA Paris) Work also supported by ANR Funding ANR-14-CE ( COSMOS ) Gabriel Stoltz (ENPC/INRIA) Aberdeen P.G., 14 Feb / 33
2 Outline Dissipative Particle Dynamics (Isothermal) Dissipative Particle Dynamics with energy conservation Parametrizing DPDE Some numerical results: shock waves in molecular systems Coarse-graining even further: SDPD 1 1 depending on the remaining time... Gabriel Stoltz (ENPC/INRIA) Aberdeen P.G., 14 Feb / 33
3 Dissipative Particle Dynamics Gabriel Stoltz (ENPC/INRIA) Aberdeen P.G., 14 Feb / 33
4 Stochastic differential equations (1) State of the system X R d, m-dimensional Brownian motion, diffusion matrix σ R d m dx t = b(x t ) dt + σ(x t ) dw t to be thought of as the limit as t 0 of (X n approximation of X n t ) X n+1 = X n + t b (X n ) + t σ(x n )G n, G n N (0, Id m ) Analytical study of the process: law ψ(t, x) of the process at time t distribution of all possible realizations of X t for a given initial distribution ψ(0, x), e.g. δ X 0 and all realizations of the Brownian motion Gabriel Stoltz (ENPC/INRIA) Aberdeen P.G., 14 Feb / 33
5 Stochastic differential equations (2) Averages at time t ( ) E ϕ(x t ) = ϕ(x) ψ(t, x) dx R d Generator L = b(x) σσt (x) : 2 = d b i (x) xi i=1 d i,j=1 [ ] σσ T (x) x i xj i,j Evolution of averages at time t ) t (E [ϕ(x t )] [ ] = E (Lϕ) (X t ) Proof: Taylor expansions of E[ϕ(X n+1 )] = E[ϕ(X n )] + t... Gabriel Stoltz (ENPC/INRIA) Aberdeen P.G., 14 Feb / 33
6 Stochastic differential equations (3) Rephrase as ϕ(x) t ψ(t, x) dx = R d (Lϕ)(x) ψ(t, x) dx R d Fokker-Planck equation t ψ = L ψ where L is the adjoint of L (integration by parts) ( ) (LA) (x) B(x) dx = A(x) L B (x) dx X Invariant measures are stationary solutions of the Fokker-Planck equation Invariant probability measure ψ (x) dx L ψ = 0, ψ (x) dx = 1, ψ 0 X X Gabriel Stoltz (ENPC/INRIA) Aberdeen P.G., 14 Feb / 33
7 Langevin dynamics Positions q = (q 1,..., q N ) D and momenta p = (p 1,..., p N ) R dn Hamiltonian H(q, p) = V (q) pt M 1 p Stochastic perturbation of the Hamiltonian dynamics dq t = M 1 p t dt dp t = V (q t ) dt γm 1 2γ p t dt + β dw t Friction γ > 0 (could be a position-dependent matrix) Existence and uniqueness of the invariant measure (canonical measure) µ(dq dp) = Z 1 e βh(q,p) dq dp Various ergodicity results (including exponential convergence of the law) 2 2 T. Lelièvre and G. Stoltz, Acta Numerica (2016) Gabriel Stoltz (ENPC/INRIA) Aberdeen P.G., 14 Feb / 33
8 Dissipative Particle Dynamics (1) Langevin dynamics not Galilean invariant, hence not consistent with hydrodynamics friction forces depending on relative velocities 3 Dissipative Particle Dynamics dq t = M 1 p t dt ( dp i,t = qi V (q t ) dt + i j with r ij = q i q j and v ij = p i m i p j m j γχ 2 (r ij,t )v ij,t dt + Antisymmetric stochastic forcing: W ij = W ji Cut-off function χ for the interaction, e.g. χ(r) = ) 2γ β χ(r ij,t) dw ij ( 1 r ) 1 r rcut r cut 3 P. J. Hoogerbrugge and J. M. V. A. Koelman, Europhys. Lett. (1992) Gabriel Stoltz (ENPC/INRIA) Aberdeen P.G., 14 Feb / 33
9 Dissipative Particle Dynamics (2) Can be generalized to anisotropic frictions, e.g. only along e ij = r ij r ij Generator L = L ham + γ L FD,ij with 1 i<j N L FD,ij = χ 2 (r ij ) ( v ij ( pi pj ) + 1β ) ( p1 p2 ) 2 Thermodynamic consistency 4 Invariant measure µ DPD (dq dp) = Z 1 e βh(q,p) δ { N i=1 p i P 0} (dp) dq Ergodicity is an issue proved for d = 1 only 5 4 P. Espanol and P. Warren, Europhys. Lett. (1995) 5 T. Shardlow and Y. Yan, Stoch. Dynam. (2006) Gabriel Stoltz (ENPC/INRIA) Aberdeen P.G., 14 Feb / 33
10 Dissipative Particle Dynamics: numerical schemes Many schemes, often dubious extensions of Verlet-like/BBK schemes More appropriate alternatives: splitting schemes 67 My favorite one: Verlet step to update (q, p) For each pair of particles, analytically integrate the elementary fluctuation/dissipation dynamics L FD,ij ( 1 dv ij = γ + 1 ) ( χ 2 2γ 1 (r ij,t )v ij,t dt+ + 1 ) χ(r ij,t ) dw ij m i m j β m i m j 6 T. Shardlow, SIAM J. Sci. Comput. (2003) 7 B. Leimkuhler and X. Shang, J. Comput. Phys. (2015) Gabriel Stoltz (ENPC/INRIA) Aberdeen P.G., 14 Feb / 33
11 Dissipative Particle Dynamics with energy preservation Gabriel Stoltz (ENPC/INRIA) Aberdeen P.G., 14 Feb / 33
12 Construction of the dynamics: a simplified case (1) Enforcing energy conservation in the Langevin dynamics? Gather energy variation in some additional variable ε R + ( ) ( ) dε t = d H(q t, p t ) = γ p2 t m σ2 dt σ p t 2m m dw t The variable ε can be interpreted as some coarse-grained internal energy entropy s(ε) dq t = p t m dt dp t = V (q t ) dt γ p t m dt + dε t = ( γ p2 t m σ2 2m Generator L ham + L FD with L FD = γ 2γ β dw t, ) dt σ p t m dw t. p m A + σ2 2 A2 and A = p p m ε Gabriel Stoltz (ENPC/INRIA) Aberdeen P.G., 14 Feb / 33
13 Construction of the dynamics: a simplified case (2) Invariance of measures of the form (mind the factor e s(ε)...) ( ) ρ(dq dp dε) = f H(q, p) + ε e s(ε) dq dp dε Sufficient condition L FD ρ = 0, i.e. γ p ( ) σ 2 m ρ + A 2 ρ = 0 Possible choice: γ = γ(ε) and σ R Fluctuation/dissipation relation T (ε) = 1 s (ε) σ 2 = 2γ s (ε) interpreted as some internal temperature Gabriel Stoltz (ENPC/INRIA) Aberdeen P.G., 14 Feb / 33
14 DPDE for many-particle systems (1) Coarse-graining interpretation: a (fragment of a) molecule is replaced by a mesoparticle (q i, p i ) describes the center of mass of the ith mesoparticle missing degrees of freedom described by an internal energy ε i Evolution at constant total energy H(q, p, ε) = V (q) + N i=1 p 2 i 2m i + Microscopic state law: entropies s i = s i (ε i ), internal temperature defined from the entropy as T i (ε i ) = 1 s i (ε i) Simplest case: harmonic internal degrees of freedom, T (ε) = ε/c v J. Bonet Avalos and A. Mackie, Europhys. Lett. 40, (1997) P. Español, Europhys. Lett (1997) Gabriel Stoltz (ENPC/INRIA) Aberdeen P.G., 14 Feb / 33 N i=1 ε i
15 DPDE for many-particle systems (2) dq i = p i dt m i dp i = qi V (q)dt + γ ij χ 2 (r ij )v ij dt + σ ij χ(r ij ) dw ij, i j dε i = 1 2 ( χ 2 (r ij ) j i Invariant measures γ ij v 2 ij σ2 ij 2 ( ) ρ(dq dp dε) = f H(q, p, ε) g Fluctuation-dissipation relation ( 1 m i + 1 m j ) ) dt σ ij χ(r ij)v ij dw ij. ( N ) ( N ) p i exp s i (ε i ) dq dp dε, i=1 σ ij = σ, γ ij = σ2 β ij (ε i, ε j ), β ij (ε i, ε j ) = 1 2 2k B i=1 ( 1 T i (ε i ) + 1 ) T j (ε j ) Gabriel Stoltz (ENPC/INRIA) Aberdeen P.G., 14 Feb / 33
16 Thermodynamic properties In favorable situations (...) the dynamics should be ergodic with respect to the measure ρ E0,P 0 = Z 1 E 0,P 0 δ {H(q,p,ε) E0 }δ { N i=1 p i P 0} exp ( N ) s i (ε i ) When P 0 = 0, this measure is equivalent in the thermodynamic limit to the canonical measure N ρ β (dq dp dε) = Z 1 β e βh(q,p) e βf i (ε i ) dq dp dε where f i (ε i ) = ε i s i(ε i ) β Estimators of the temperature T kin (p) = 1 dnk B i=1 i=1 is a free energy and β is such that E ρβ (H) = E 0 N i=1 p 2 i m i, T int (ε) = ( 1 N N i=1 ) 1 1 T i (ε i ) Gabriel Stoltz (ENPC/INRIA) Aberdeen P.G., 14 Feb / 33
17 Adding thermal conduction Additional elementary pairwise fluctuation/dissipation dynamics ( 1 dε i = κχ 2 (r ij ) T i (ε i ) 1 ) dt + 2κχ(r ij ) d T j (ε j ) W ij, dε j = dε i, where κ > 0 is a thermal conductivity This dynamis preserves the elementary energy ε i + ε j by construction measures of the form f (ε i + ε j ) e s i (ε i )+s j (ε j ) dε i dε j Static thermodynamic properties therefore unchanged Gabriel Stoltz (ENPC/INRIA) Aberdeen P.G., 14 Feb / 33
18 Numerical integration of DPDE (1) Splitting strategy: Hamiltonian part and elementary dynamics dp i = γ(ε i, ε j )χ 2 (r ij )v ij dt + σχ(r ij ) dw ij, dp j = dp i, dε i = χ2 (r ij ) 2 dε j = dε i, [ ( γ(ε i, ε j )vij 2 σ )] dt σ 2 m i m j 2 χ(r ij)v ij dw ij, In fact, can be reduced to a dynamics on v ij only ( 1 dv ij = γ(ε i, ε j )χ 2 (r ij ) + 1 ) ( 1 v ij dt + σχ(r ij ) + 1 m i m j m i m j ) dw ij since, by the energy conservation (recall p i + p j is constant), ( ) ε i = ε i,0 + 1 p 2 i pi,0 2 + p2 j pj,0 2 = ε i,0 +F (v ij ), ε j = ε j,0 +F (v ij ) 2 2m i 2m j Gabriel Stoltz (ENPC/INRIA) Aberdeen P.G., 14 Feb / 33
19 Numerical integration of DPDE (2) Practical integrator: integrate the dynamics on v ij with fixed friction update the internal energies to ensure the energy conservation it is possible to superimpose a Metropolis correction 8 [w.r.t. some locally invariant measure] Pro/cons of this integrator: automatically corrects for negative internal energies (stabilization) parallelization/threadability limited dedicated schemes for that 9 8 G. Stoltz, arxiv preprint A.-A. Homman, J.-B. Maillet, J. Roussel and G. Stoltz, J. Chem. Phys (2016) Gabriel Stoltz (ENPC/INRIA) Aberdeen P.G., 14 Feb / 33
20 Parametrization of DPDE Gabriel Stoltz (ENPC/INRIA) Aberdeen P.G., 14 Feb / 33
21 Summary of the quantities to provide Static properties effective masses for mesoparticles microscopic state law ε = T (ε) 0 C v (θ) dθ effective interaction potential machine learning approaches Dynamical properties fluctuation magnitude σ thermal conductivity κ cut-off functions for the fluctuation/dissipation dynamics Derivation of DPDE? From Hamiltonian dynamics P. Español, M. Serrano, I. Pagonabarraga and I. Zúñiga, Soft Matter (2016) Gabriel Stoltz (ENPC/INRIA) Aberdeen P.G., 14 Feb / 33
22 Estimating some of the quantities Microscopic state law: energy as a function of the temperature ab-initio computations of the rovibrational modes integration using some distribution function (Bose Einstein) invert to obtain the temperature as a function of the energy Fluctuation magnitude: equilibration dynamics 11 initialize internal and external degrees of freedom at different temperatures N ρ βext,βint (dq dp dε) = Z 1 β e βexth(q,p) e β intf i (ε i ) dq dp dε compute the time evolution of internal and external (kinetic) temperatures timescale for return to equilibrium reference all-atom simulation (small system) i=1 11 M. Kroonblawd, T. Sewell and J.-B. Maillet, J. Chem. Phys. (2016) Gabriel Stoltz (ENPC/INRIA) Aberdeen P.G., 14 Feb / 33
23 Some numerical results Gabriel Stoltz (ENPC/INRIA) Aberdeen P.G., 14 Feb / 33
24 Shock waves in an effective material Temperature (K) Time (ps) G. Stoltz, Europhys. Lett. (2006) Gabriel Stoltz (ENPC/INRIA) Aberdeen P.G., 14 Feb / 33
25 Detonation waves in nitromethane (1) Potential of Exp-6 type, fitted to reproduce Hugoniot curves E E (P + P 0)(V 0 V) = 0 C v taken from thermodynamic tables σ so that equilibration time between internal and external degrees of freedom is of the order of a few ps Modeling of chemical reactions through some progress variables λ i System of length 4.5 µm Gabriel Stoltz (ENPC/INRIA) Aberdeen P.G., 14 Feb / 33
26 Detonation waves in nitromethane (2) Particle velocity, temperature, progress variable, pressure J.-B. Maillet, G. Vallverdu, N. Desbiens and G. Stoltz, Europhys. Lett. (2011) Gabriel Stoltz (ENPC/INRIA) Aberdeen P.G., 14 Feb / 33
27 Coupling with higher scales Gabriel Stoltz (ENPC/INRIA) Aberdeen P.G., 14 Feb / 33
28 Coarse-graining even further: SDPD (1) Adding thermal fluctuations to particle discretization of Navier Stokes 12 Size K of particles (masses Km ) Variables: (q, p) position of discretization nodes, entropy S i or energy ε i Conservative part of the dynamics: Hamiltonian dynamics with energy E(q, p, S) = N i=1 ε i ( Si, ρ i (q) ) + p2 i 2m, ρ i(q) = N mw (r ij ) j=1 dq = p m dt, dp = q E(q, p, S) dt, ds = P. Español and M. Revenga, Phys. Rev. E (2003) Gabriel Stoltz (ENPC/INRIA) Aberdeen P.G., 14 Feb / 33
29 Coarse-graining even further: SDPD (2) Fluctuation/dissipation: energy formulation 13 dp i = γijp θ ijv θ ij dt + σijp θ ij θ db ij, dp j = dp i, [ dε i = 1 2 dε j = dε j, γ θ ijv T ] ij Pijv θ ij (σθ ij )2 m Tr(Pθ ij) dt 1 2 σθ ijv T ij P θ ij db ij, with appropriate choices for γ ij, γ ij, σ ij, σ ij to ensure the invariance of µ(dq dp dε) = g ( E(q, p, ε), ) N N p i i=1 i=1 ( ) exp Si (ε i,q) k B T i (ε i, q) Reduces to Smoothed Particle Hydrodynamics when σ ij = 0 13 G. Faure, J.-B. Maillet and G. Stoltz, Phys. Rev. E (2016) dq dp dε Gabriel Stoltz (ENPC/INRIA) Aberdeen P.G., 14 Feb / 33
30 Parametrization of SDPD Free energy F (ρ, T ), computed from NVT simulations. Deduce... ( ) F (ρ, T ) energy E(ρ, T ) = T 2 T T E(ρ, T ) F (ρ, T ) entropy S(ρ, T ) = T pressure P(ρ, T ) = ρ 2 ρ F (ρ, T ) heat capacity C(ρ, T ) = T E(ρ, T ) Fluid viscosities, which appear in γ ij, σ ij Gabriel Stoltz (ENPC/INRIA) Aberdeen P.G., 14 Feb / 33
31 SDPD: some numerical results MD NS K = 10 K = 100 K = 500 K = 1,000 K = 5,000 K = 10,000 Density (kg.m 3 ) 1,600 1,400 1, Distance from the shock front (reduced units) Density (kg.m 3 ) 1,400 1, Distance from the shock front (nm) Gabriel Stoltz (ENPC/INRIA) Aberdeen P.G., 14 Feb / 33
32 Conclusion Gabriel Stoltz (ENPC/INRIA) Aberdeen P.G., 14 Feb / 33
33 Conclusion and perspectives Coarse-grained models for atomistic simulations Dissipative Particle Dynamics with conserved energy (DPDE) can be used in nonequilibrium situations replace a molecule or some group of atoms by a mesoparticle consistent thermodynamics input: static properties (computed ab-initio) and dynamical parameters Smoothed Dissipative Particle Dynamics (SDPD) hydrodynamics-like (e.g. solvent in bulk-like regime) but particle-based: seamless coupling with DPDE can be envisioned? in progress: add chemistry (Gérôme Faure, CEA/DAM) Gabriel Stoltz (ENPC/INRIA) Aberdeen P.G., 14 Feb / 33
Approches mathématiques pour la simulation multi-échelle des matériaux
p. 1/2 Approches mathématiques pour la simulation multi-échelle des matériaux Claude Le Bris Ecole des Ponts and Inria, FRANCE http://cermics.enpc.fr/ lebris based on a series of works by F. Legoll, T.
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 informationAPMA 2811T. By Zhen Li. Today s topic: Lecture 2: Theoretical foundation and parameterization. Sep. 15, 2016
Today s topic: APMA 2811T Dissipative Particle Dynamics Instructor: Professor George Karniadakis Location: 170 Hope Street, Room 118 Time: Thursday 12:00pm 2:00pm Dissipative Particle Dynamics: Foundation,
More informationNumerical methods for computational statistical physics
Numerical methods for computational statistical physics 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 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 informationHydrodynamics, Thermodynamics, and Mathematics
Hydrodynamics, Thermodynamics, and Mathematics Hans Christian Öttinger Department of Mat..., ETH Zürich, Switzerland Thermodynamic admissibility and mathematical well-posedness 1. structure of equations
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 informationModeling of Micro-Fluidics by a Dissipative Particle Dynamics Method. Justyna Czerwinska
Modeling of Micro-Fluidics by a Dissipative Particle Dynamics Method Justyna Czerwinska Scales and Physical Models years Time hours Engineering Design Limit Process Design minutes Continious Mechanics
More informationMD Thermodynamics. Lecture 12 3/26/18. Harvard SEAS AP 275 Atomistic Modeling of Materials Boris Kozinsky
MD Thermodynamics Lecture 1 3/6/18 1 Molecular dynamics The force depends on positions only (not velocities) Total energy is conserved (micro canonical evolution) Newton s equations of motion (second order
More information2 Dissipative particle dynamics with energy conservation: Heat conduction that is dissipated due to the velocity dependent forces is transformed into
International Journal of Modern Physics C, fc World Scientic Publishing Company DISSIPATIVE PARTICLE DYNAMICS WITH ENERGY CONSERVATION: HEAT CONDUCTION M. RIPOLL, P. ESPA ~NOL Departamento de Fsica Fundamental,
More informationThermostatic Controls for Noisy Gradient Systems and Applications to Machine Learning
Thermostatic Controls for Noisy Gradient Systems and Applications to Machine Learning Ben Leimkuhler University of Edinburgh Joint work with C. Matthews (Chicago), G. Stoltz (ENPC-Paris), M. Tretyakov
More informationBROWNIAN DYNAMICS SIMULATIONS WITH HYDRODYNAMICS. Abstract
BROWNIAN DYNAMICS SIMULATIONS WITH HYDRODYNAMICS Juan J. Cerdà 1 1 Institut für Computerphysik, Pfaffenwaldring 27, Universität Stuttgart, 70569 Stuttgart, Germany. (Dated: July 21, 2009) Abstract - 1
More informationEfficient numerical integrators for stochastic models
Efficient numerical integrators for stochastic models G. De Fabritiis a, M. Serrano b P. Español b P.V.Coveney a a Centre for Computational Science, Department of Chemistry, University College London,
More informationMesoscale Simulation Methods. Ronojoy Adhikari The Institute of Mathematical Sciences Chennai
Mesoscale Simulation Methods Ronojoy Adhikari The Institute of Mathematical Sciences Chennai Outline What is mesoscale? Mesoscale statics and dynamics through coarse-graining. Coarse-grained equations
More informationFrom Molecular Dynamics to hydrodynamics a novel Galilean invariant thermostat
To appear in: Journal of Chemical Physics, 08 APR 005 From Molecular Dynamics to hydrodynamics a novel Galilean invariant thermostat Simeon D. Stoyanov and Robert D. Groot Unilever Research Vlaardingen,
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 informationCoarse-graining limits in open and wall-bounded dissipative particle dynamics systems
THE JOURNAL OF CHEMICAL PHYSICS 124, 184101 2006 Coarse-graining limits in open and wall-bounded dissipative particle dynamics systems Igor V. Pivkin and George E. Karniadakis a Division of Applied Mathematics,
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 informationPhase Equilibria and Molecular Solutions Jan G. Korvink and Evgenii Rudnyi IMTEK Albert Ludwig University Freiburg, Germany
Phase Equilibria and Molecular Solutions Jan G. Korvink and Evgenii Rudnyi IMTEK Albert Ludwig University Freiburg, Germany Preliminaries Learning Goals Phase Equilibria Phase diagrams and classical thermodynamics
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 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 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 informationCaltech Ph106 Fall 2001
Caltech h106 Fall 2001 ath for physicists: differential forms Disclaimer: this is a first draft, so a few signs might be off. 1 Basic properties Differential forms come up in various parts of theoretical
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 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 informationHydrodynamics. Stefan Flörchinger (Heidelberg) Heidelberg, 3 May 2010
Hydrodynamics Stefan Flörchinger (Heidelberg) Heidelberg, 3 May 2010 What is Hydrodynamics? Describes the evolution of physical systems (classical or quantum particles, fluids or fields) close to thermal
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 informationSmoothed Dissipative Particle Dynamics: theory and applications to complex fluids
2015 DPD Workshop September 21-23, 2015, Shanghai University Smoothed Dissipative Particle Dynamics: Dynamics theory and applications to complex fluids Marco Ellero Zienkiewicz Centre for Computational
More informationDissipative nuclear dynamics
Dissipative nuclear dynamics Curso de Reacciones Nucleares Programa Inter universitario de Fisica Nuclear Universidad de Santiago de Compostela March 2009 Karl Heinz Schmidt Collective dynamical properties
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 informationCollective behavior, from particles to fields
978-0-51-87341-3 - Statistical Physics of Fields 1 Collective behavior, from particles to fields 1.1 Introduction One of the most successful aspects of physics in the twentieth century was revealing the
More informationIntroduction. Statistical physics: microscopic foundation of thermodynamics degrees of freedom 2 3 state variables!
Introduction Thermodynamics: phenomenological description of equilibrium bulk properties of matter in terms of only a few state variables and thermodynamical laws. Statistical physics: microscopic foundation
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 informationProtein Dynamics, Allostery and Function
Protein Dynamics, Allostery and Function Lecture 3. Protein Dynamics Xiaolin Cheng UT/ORNL Center for Molecular Biophysics SJTU Summer School 2017 1 Obtaining Dynamic Information Experimental Approaches
More informationCOARSE-GRAINING AND THERMODYNAMICS IN FAR-FROM-EQUILIBRIUM SYSTEMS
Vol. 44 (2013) ACTA PHYSICA POLONICA B No 5 COARSE-GRAINING AND THERMODYNAMICS IN FAR-FROM-EQUILIBRIUM SYSTEMS J. Miguel Rubí, A. Pérez-Madrid Departament de Física Fonamental, Facultat de Física, Universitat
More informationThe Microcanonical Approach. (a) The volume of accessible phase space for a given total energy is proportional to. dq 1 dq 2 dq N dp 1 dp 2 dp N,
8333: Statistical Mechanics I Problem Set # 6 Solutions Fall 003 Classical Harmonic Oscillators: The Microcanonical Approach a The volume of accessible phase space for a given total energy is proportional
More informationAA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS
AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS 1 / 29 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS Hierarchy of Mathematical Models 1 / 29 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS 2 / 29
More informationEstimating Accuracy in Classical Molecular Simulation
Estimating Accuracy in Classical Molecular Simulation University of Illinois Urbana-Champaign Department of Computer Science Institute for Mathematics and its Applications July 2007 Acknowledgements Ben
More 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 informationarxiv: v7 [quant-ph] 22 Aug 2017
Quantum Mechanics with a non-zero quantum correlation time Jean-Philippe Bouchaud 1 1 Capital Fund Management, rue de l Université, 75007 Paris, France. (Dated: October 8, 018) arxiv:170.00771v7 [quant-ph]
More informationStochastic Particle Dynamics for Unresolved Degrees of Freedom. Sebastian Reich in collaborative work with Colin Cotter (Imperial College)
Stochastic Particle Dynamics for Unresolved Degrees of Freedom Sebastian Reich in collaborative work with Colin Cotter (Imperial College) 1. The Motivation Classical molecular dynamics (MD) as well as
More informationIntroduction Statistical Thermodynamics. Monday, January 6, 14
Introduction Statistical Thermodynamics 1 Molecular Simulations Molecular dynamics: solve equations of motion Monte Carlo: importance sampling r 1 r 2 r n MD MC r 1 r 2 2 r n 2 3 3 4 4 Questions How can
More information2 Equations of Motion
2 Equations of Motion system. In this section, we will derive the six full equations of motion in a non-rotating, Cartesian coordinate 2.1 Six equations of motion (non-rotating, Cartesian coordinates)
More informationWhat is Classical Molecular Dynamics?
What is Classical Molecular Dynamics? Simulation of explicit particles (atoms, ions,... ) Particles interact via relatively simple analytical potential functions Newton s equations of motion are integrated
More informationChapter 1. Governing Equations of GFD. 1.1 Mass continuity
Chapter 1 Governing Equations of GFD The fluid dynamical governing equations consist of an equation for mass continuity, one for the momentum budget, and one or more additional equations to account for
More informationmacroscopic view (phenomenological) microscopic view (atomistic) computing reaction rate rate of reactions experiments thermodynamics
Rate heory (overview) macroscopic view (phenomenological) rate of reactions experiments thermodynamics Van t Hoff & Arrhenius equation microscopic view (atomistic) statistical mechanics transition state
More informationAPMA 2811T. By Zhen Li. Today s topic: Lecture 3: New Methods beyond traditional DPD. Sep. 22, Division of Applied Mathematics, Brown University
Today s topic: APMA 2811T Dissipative Particle Dynamics Instructor: Professor George Karniadakis Location: 170 Hope Street, Room 118 Time: Thursday 12:00pm 2:00pm Dissipative Particle Dynamics: Foundation,
More informationSmoothed Dissipative Particle Dynamics Model for Predicting Self-Assembled Nano-Cellulose Fibre Structures
Smoothed Dissipative Particle Dynamics Model for Predicting Self-Assembled Nano-Cellulose Fibre Structures David Vidal and Tetsu Uesaka FPInnovations, Pointe-Claire, Québec, CANADA Nano-cellulose fibres
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 informationCoarse-graining, applications and mesoscopic molecular dynamics
CS Work in progress Coarse-graining, applications and mesoscopic molecular dynamics Carsten Svaneborg Institute for Physics, Chemistry, and Pharmacy University of Southern Denmark Campusvej 55, 5320 Odense
More informationIntroduction to Computer Simulations of Soft Matter Methodologies and Applications Boulder July, 19-20, 2012
Introduction to Computer Simulations of Soft Matter Methodologies and Applications Boulder July, 19-20, 2012 K. Kremer Max Planck Institute for Polymer Research, Mainz Overview Simulations, general considerations
More informationIII. Kinetic Theory of Gases
III. Kinetic Theory of Gases III.A General Definitions Kinetic theory studies the macroscopic properties of large numbers of particles, starting from their (classical) equations of motion. Thermodynamics
More informationINTRODUCTION TO о JLXJLA Из А lv-/xvj_y JrJrl Y üv_>l3 Second Edition
INTRODUCTION TO о JLXJLA Из А lv-/xvj_y JrJrl Y üv_>l3 Second Edition Kerson Huang CRC Press Taylor & Francis Croup Boca Raton London New York CRC Press is an imprint of the Taylor & Francis Group an Informa
More informationCapSel Euler The Euler equations. conservation laws for 1D dynamics of compressible gas. = 0 m t + (m v + p) x
CapSel Euler - 01 The Euler equations keppens@rijnh.nl conservation laws for 1D dynamics of compressible gas ρ t + (ρ v) x = 0 m t + (m v + p) x = 0 e t + (e v + p v) x = 0 vector of conserved quantities
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 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 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 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 informationSet the initial conditions r i. Update neighborlist. r i. Get new forces F i
Set the initial conditions r i t 0, v i t 0 Update neighborlist Get new forces F i r i Solve the equations of motion numerically over time step t : r i t n r i t n + v i t n v i t n + Perform T, P scaling
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 informationAnomalous hydrodynamics and gravity. Dam T. Son (INT, University of Washington)
Anomalous hydrodynamics and gravity Dam T. Son (INT, University of Washington) Summary of the talk Hydrodynamics: an old theory, describing finite temperature systems The presence of anomaly modifies hydrodynamics
More informationParticle-Simulation Methods for Fluid Dynamics
Particle-Simulation Methods for Fluid Dynamics X. Y. Hu and Marco Ellero E-mail: Xiangyu.Hu and Marco.Ellero at mw.tum.de, WS 2012/2013: Lectures for Mechanical Engineering Institute of Aerodynamics Technical
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 informationLinear Response Theory, The Green-Kubo formula and Langevin Dynamics
Linear Response Theory, The Green-Kubo formula and Langevin Dynamics G.A. Pavliotis Department of Mathematics Imperial College London 07/07/2015 G.A. Pavliotis (IC London) Linear Response Theory 1 / 64
More informationAn alternative approach to dissipative particle dynamics
EUROPHYSICS LETTERS 15 July 1999 Europhys. Lett., 47 (2), pp. 145-151 (1999) An alternative approach to dissipative particle dynamics C. P. Lowe Computational Physics, Delft University of Technology Lorentzweg
More informationHydrodynamic Modes of Incoherent Black Holes
Hydrodynamic Modes of Incoherent Black Holes Vaios Ziogas Durham University Based on work in collaboration with A. Donos, J. Gauntlett [arxiv: 1707.xxxxx, 170x.xxxxx] 9th Crete Regional Meeting on String
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 informationCollective Effects. Equilibrium and Nonequilibrium Physics
Collective Effects in Equilibrium and Nonequilibrium Physics: Lecture 4, April 7, 2006 1 Collective Effects in Equilibrium and Nonequilibrium Physics Website: http://cncs.bnu.edu.cn/mccross/course/ Caltech
More informationLinear Response and Onsager Reciprocal Relations
Linear Response and Onsager Reciprocal Relations Amir Bar January 1, 013 Based on Kittel, Elementary statistical physics, chapters 33-34; Kubo,Toda and Hashitsume, Statistical Physics II, chapter 1; and
More informationComparative Study of the Water Response to External Force at Nanoscale and Mesoscale
Copyright 2013 Tech Science Press CMES, vol.95, no.4, pp.303-315, 2013 Comparative Study of the Water Response to External Force at Nanoscale and Mesoscale H.T. Liu 1,2, Z. Chen 2, S. Jiang 2, Y. Gan 3,
More informationarxiv:comp-gas/ v1 28 Apr 1993
Lattice Boltzmann Thermohydrodynamics arxiv:comp-gas/9304006v1 28 Apr 1993 F. J. Alexander, S. Chen and J. D. Sterling Center for Nonlinear Studies and Theoretical Division Los Alamos National Laboratory
More information3.320: Lecture 19 (4/14/05) Free Energies and physical Coarse-graining. ,T) + < σ > dµ
3.320: Lecture 19 (4/14/05) F(µ,T) = F(µ ref,t) + < σ > dµ µ µ ref Free Energies and physical Coarse-graining T S(T) = S(T ref ) + T T ref C V T dt Non-Boltzmann sampling and Umbrella sampling Simple
More informationA Hamiltonian Numerical Scheme for Large Scale Geophysical Fluid Systems
A Hamiltonian Numerical Scheme for Large Scale Geophysical Fluid Systems Bob Peeters Joint work with Onno Bokhove & Jason Frank TW, University of Twente, Enschede CWI, Amsterdam PhD-TW colloquium, 9th
More informationAn electrokinetic LB based model for ion transport and macromolecular electrophoresis
An electrokinetic LB based model for ion transport and macromolecular electrophoresis Raffael Pecoroni Supervisor: Michael Kuron July 8, 2016 1 Introduction & Motivation So far an mesoscopic coarse-grained
More informationMesoscopic simulation for the structural change of a surfactant solution using dissipative particle dynamics
Korean J. Chem. Eng., 26(6), 1717-1722 (2009) DOI: 10.1007/s11814-009-0235-2 RAPID COMMUNICATION Mesoscopic simulation for the structural change of a surfactant solution using dissipative particle dynamics
More informationSupplementary Information for: Controlling Cellular Uptake of Nanoparticles with ph-sensitive Polymers
Supplementary Information for: Controlling Cellular Uptake of Nanoparticles with ph-sensitive Polymers Hong-ming Ding 1 & Yu-qiang Ma 1,2, 1 National Laboratory of Solid State Microstructures and Department
More informationMajor Concepts Lecture #11 Rigoberto Hernandez. TST & Transport 1
Major Concepts Onsager s Regression Hypothesis Relaxation of a perturbation Regression of fluctuations Fluctuation-Dissipation Theorem Proof of FDT & relation to Onsager s Regression Hypothesis Response
More informationV (r,t) = i ˆ u( x, y,z,t) + ˆ j v( x, y,z,t) + k ˆ w( x, y, z,t)
IV. DIFFERENTIAL RELATIONS FOR A FLUID PARTICLE This chapter presents the development and application of the basic differential equations of fluid motion. Simplifications in the general equations and common
More informationVIII.B Equilibrium Dynamics of a Field
VIII.B Equilibrium Dynamics of a Field The next step is to generalize the Langevin formalism to a collection of degrees of freedom, most conveniently described by a continuous field. Let us consider the
More informationTransport properties controlled by a thermostat: An extended dissipative particle dynamics thermostat
PAPER www.rsc.org/softmatter Soft Matter Transport properties controlled by a thermostat: An extended dissipative particle dynamics thermostat Christoph Junghans, a Matej Praprotnik ab and Kurt Kremer
More informationLangevin Dynamics in Constant Pressure Extended Systems
Langevin Dynamics in Constant Pressure Extended Systems D. Quigley and M.I.J. Probert CMMP 2004 1 Talk Outline Molecular dynamics and ensembles. Existing methods for sampling at NPT. Langevin dynamics
More informationLANGEVIN EQUATION AND THERMODYNAMICS
LANGEVIN EQUATION AND THERMODYNAMICS RELATING STOCHASTIC DYNAMICS WITH THERMODYNAMIC LAWS November 10, 2017 1 / 20 MOTIVATION There are at least three levels of description of classical dynamics: thermodynamic,
More informationAspects of nonautonomous molecular dynamics
Aspects of nonautonomous molecular dynamics IMA, University of Minnesota, Minneapolis January 28, 2007 Michel Cuendet Swiss Institute of Bioinformatics, Lausanne, Switzerland Introduction to the Jarzynski
More informationUNIVERSITY OF OSLO FACULTY OF MATHEMATICS AND NATURAL SCIENCES
UNIVERSITY OF OSLO FCULTY OF MTHEMTICS ND NTURL SCIENCES Exam in: FYS430, Statistical Mechanics Day of exam: Jun.6. 203 Problem :. The relative fluctuations in an extensive quantity, like the energy, depends
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 informationDiffusive Transport Enhanced by Thermal Velocity Fluctuations
Diffusive Transport Enhanced by Thermal Velocity Fluctuations Aleksandar Donev 1 Courant Institute, New York University & Alejandro L. Garcia, San Jose State University John B. Bell, Lawrence Berkeley
More informationNon equilibrium thermodynamics: foundations, scope, and extension to the meso scale. Miguel Rubi
Non equilibrium thermodynamics: foundations, scope, and extension to the meso scale Miguel Rubi References S.R. de Groot and P. Mazur, Non equilibrium Thermodynamics, Dover, New York, 1984 J.M. Vilar and
More informationStochastic Particle Methods for Rarefied Gases
CCES Seminar WS 2/3 Stochastic Particle Methods for Rarefied Gases Julian Köllermeier RWTH Aachen University Supervisor: Prof. Dr. Manuel Torrilhon Center for Computational Engineering Science Mathematics
More informationFrom the Newton equation to the wave equation in some simple cases
From the ewton equation to the wave equation in some simple cases Xavier Blanc joint work with C. Le Bris (EPC) and P.-L. Lions (Collège de France) Université Paris Diderot, FRACE http://www.ann.jussieu.fr/
More informationLangevin Dynamics of a Single Particle
Overview Langevin Dynamics of a Single Particle We consider a spherical particle of radius r immersed in a viscous fluid and suppose that the dynamics of the particle depend on two forces: A drag force
More informationThe Boltzmann Equation and Its Applications
Carlo Cercignani The Boltzmann Equation and Its Applications With 42 Illustrations Springer-Verlag New York Berlin Heidelberg London Paris Tokyo CONTENTS PREFACE vii I. BASIC PRINCIPLES OF THE KINETIC
More informationThus, the volume element remains the same as required. With this transformation, the amiltonian becomes = p i m i + U(r 1 ; :::; r N ) = and the canon
G5.651: Statistical Mechanics Notes for Lecture 5 From the classical virial theorem I. TEMPERATURE AND PRESSURE ESTIMATORS hx i x j i = kt ij we arrived at the equipartition theorem: * + p i = m i NkT
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 informationSymmetry of the Dielectric Tensor
Symmetry of the Dielectric Tensor Curtis R. Menyuk June 11, 2010 In this note, I derive the symmetry of the dielectric tensor in two ways. The derivations are taken from Landau and Lifshitz s Statistical
More informationn v molecules will pass per unit time through the area from left to
3 iscosity and Heat Conduction in Gas Dynamics Equations of One-Dimensional Gas Flow The dissipative processes - viscosity (internal friction) and heat conduction - are connected with existence of molecular
More informationA comparative study between dissipative particle dynamics and molecular dynamics for simple- and complex-geometry flows
THE JOURNAL OF CHEMICAL PHYSICS 123, 104107 2005 A comparative study between dissipative particle dynamics and molecular dynamics for simple- and complex-geometry flows Eric E. Keaveny, Igor V. Pivkin,
More information1 Garrod #5.2: Entropy of Substance - Equation of State. 2
Dylan J. emples: Chapter 5 Northwestern University, Statistical Mechanics Classical Mechanics and hermodynamics - Garrod uesday, March 29, 206 Contents Garrod #5.2: Entropy of Substance - Equation of State.
More informationMathematical Structures of Statistical Mechanics: from equilibrium to nonequilibrium and beyond Hao Ge
Mathematical Structures of Statistical Mechanics: from equilibrium to nonequilibrium and beyond Hao Ge Beijing International Center for Mathematical Research and Biodynamic Optical Imaging Center Peking
More informationTemperature and Pressure Controls
Ensembles Temperature and Pressure Controls 1. (E, V, N) microcanonical (constant energy) 2. (T, V, N) canonical, constant volume 3. (T, P N) constant pressure 4. (T, V, µ) grand canonical #2, 3 or 4 are
More informationarxiv: v1 [physics.chem-ph] 11 Feb 2014
Scaling properties in the adsorption of ionic polymeric surfactants on generic nanoparticles of metallic oxides by mesoscopic simulation arxiv:1402.2661v1 [physics.chem-ph] 11 Feb 2014 E. Mayoral and E.
More information