Molecular Dynamics Lecture 2
|
|
- Coral Charles
- 5 years ago
- Views:
Transcription
1 Molecular Dynamics Lecture 2 Ben Leimkuhler the problem of the timestep in md splitting methods in Langevin dynamics constraints multiple timestepping Spring School, WIAS, Berlin 2017
2 fruit fly alanine dipeptide (the fruit fly of biomolecules)
3 Measure of diffusion: exploration rate in alanine dipeptide this means frequent recrossings of a particular free energy barrier recrossings
4 Harmonic Oscillator q = p Stability Threshold ṗ = 2 q leapfrog/verlet eigenvalues: i t =2 increasing t stable for t apple i symplectic for all but only useful for small t t
5 Oscillatory modes of bio-md stability threshold 3fs 6fs 12fs calibrates the stepsize in molecular dynamics
6 Ways to improve efficiency More accurate (less-biased) approximation in equilibrium Faster convergence to equilibrium larger timesteps (less force evaluations for the same rate of phase space exploration) Today s lecture: 1. Langevin dynamics and splitting methods for better accuracy and more rapid convergence of averages 2. constraints and multiple timestepping to increase the timestep in MD
7 Generalized Langevin Dynamics (Q, P ) (q i,p i ) An instance of Mori-Zwanzig reduction deterministic model for noise process friction kernel
8 Large Bath Limit Memory kernels in various scenarios: Typically approach delta-function in the limit of large auxiliary system Noise process typically models, in the weak sense, the effect of a Wiener process [See e.g. works of A. Stuart] dq = P dt dp = ru(q)dt P dt + dw
9 Fokker Planck Equation for LD For a system Forward Kolmogorov or Fokker-Planck operator: Invariant distribution: (q, p) =e H(q,p)
10 Geometric Ergodicity of Langevin Dynamics Geometric convergence relies on establishing: (1) minorization condition (2) stability condition Key part of (1): existence of a smooth transition density, usually shown by demonstrating that the operator is hypoelliptic by use of a parabolic Hörmander condition. (2) usually shown by obtaining a Lyapunov function
11 Hypoelliptic Property dx = b 0 (X)dt + kx b i (X)dW i 2 R n i=1 L is hypoelliptic in if the SDE satisfies the parabolic Hörmander condition in Span b k i=1, [b i,b j ] i=1,2,...,k,j=0,1,...,k, [b i,b j,b k ],... = R n A mixing condition for SDEs. Langevin qi = e pi = e i+n [@ pi,p qi qi
12 Lyapunov Function A Lyapunov function for the SDE is defined as a function satisfying for some positive constants, For Langevin dynamics, as shown by [Mattingley, Higham and Stuart, 2002], we may use We can choose l to bound a sufficiently large class of observables, giving exponential convergence for that class.
13 Stronger assumption: bounded q, e.g. Periodic BCs Lyapunov function 1+ p 2s Decay ( geometric ergodicity ) Regularity
14 Splitting Methods (M=I) Time stepsize p := e t t p + q (1 e 2 t ) 1 N (0, 1) p := tf(q) q := q + tp
15 Splitting Methods L = A + B + O A = p T M 1 r q B = ru(q) T r p O = p T M 1 r p + 1 p Propagator: Splitting Method: P t = e tl P t e ta e tb e to Corresponds to ABO! (Vertauschungssatz) OBA OAB ABOBA OBABO
16 In L., Matthews & Stoltz 2015 after Talay 2002, Mattingly 2002, Bou-Rabee & Owhadi 2010, Hairer and Mattingley 2010 kp n tk apple Ke n t discrete propagator Z (P n tf)(q, p) fdµ, t apple K(1 + p 2s ) e n t kfk L 1 1+ p 2s Uniform in stepsize exponential decay i.e., geometric ergodicity of the numerical method
17 Invariant Measure of Numerical Method The stable equilibrium distribution can be understood as a perturbation of the Gibbs distribution, with density ˆ,, = e H The error in long-term averages is thus directly related to the quality of this approximation. How to calculate the error?
18 Ex: P t e ta e tb e to BCH: P t = e t[l + tl 1+ t 2 L ] Invariant density: [L + tl 1 + t2 L ]ˆ =0 For symmetric splittings: [L + t 2 L ]ˆ =0 Proposal: ˆ = (1 + t 2 d )
19 Operator expansions from Baker-Campbell-Hausdorff: e.g. for BAOAB: h L 2 2 = 4 qu(q) p T U 00 (q)p + 4 p T U 00 (q)ru(q) 12 p r qp T U 00 (q)p.
20 Expansion of the invariant distribution (Talay-Tubaro expansion 1 in the ergodic limit) Leading order: L. & Matthews, AMRX, 2013 L., Matthews, & Stoltz, IMA J. Num. Anal detailed treatment of all 1st and 2nd order splittings estimates for the operator inverse discrete inv. measure treatment of nonequilibrium (e.g. transport coefficients) 1 Denis Talay and Luciano Tubaro. Expansion of the global error for numerical schemes solving stochastic differential equations. Stochastic Analysis and Applications 4 (1990)
21 Two level expansion For each of ABOBA and BAOAB, we find the first terms of the two-level expansion " =1/
22 Configurational Sampling Integrate out with respect to momenta...and discover a surprise: Proposition: The marginal (configurational) distribution of the BAOAB method has an expansion of the form = e U (1 + O( t 2 " 2 )+O( t 4 )) In the high friction limit: 4th order, and with just one force evaluation per timestep.
23 Perturbed Quartic (1D) Langevin Configurational Accuracy vs timestep Weak Friction Strong Friction
24 Testing Methods for Biological Molecules Biomolecular models include many force terms complicated landscape Stiff bonds NH, OH.. dominant harmonic terms Presence of H2O alternate diffusion mechanism basins + barriers long sampling times/rare events Accuracy requirements not always transparent goals of simulation may include dynamics or sampling Methods need to balance diffusion rate vs accuracy of basin approximation
25 Diffusion vs Accuracy With BAOAB, we would imagine it to be better to use large friction coefficient γ 4th order (BAOAB Limit Method) But... γ strongly affects the diffusion rate Friction (1/ps) BAOAB ABOBA SPV VGB LI BP BBK Step size (fs) recrossings Friction (1/ps) Step size (fs) self-diffusion (m 2 /s)
26 Measuring Accuracy Canonical density for solvated alanine dipeptide is a function in a 3200 dimensional space. Convergence and order of accuracy depend on how we collapse the space. Examples: hp 2 i i = 1 kinetic temperature hu(q)i DX qi ru(q i )E = 1 average potential energy a configurational temp
27 Solvated Systems for alanine dipeptide in flexible water Relative Error (a) Configurational Temperature Error BBK, BP SPV, ABOBA LI VGB BAOAB Relative Error (b) Average Total Potential Energy Error BBK, BP LI SPV, ABOBA VGB BAOAB ABOBA VGB SPV BP LI BBK EM EB 10 3 BAOAB Stepsize (fs) Stepsize (fs) BAOAB much better than alternatives for relevant (i.e. configurational) quantities up to the stability threshold and At moderate friction!
28 If we are not in the high friction limit, why is BAOAB so much better? Harmonic Oscillator: ABOBA and BAOAB are exact for average potential
29 Harmonic Oscillator Configurational Sampling Error in average q st order line BBK, BP VGB 2nd order line EB EM LI SPV BAOAB, ABOBA BAOAB ABOBA VGB SPV BP LI BBK EM EB Stepsize BAOAB and ABOBA both are exact for PE But...this is only part of the story since BAOAB is much better than ABOBA for real molecules
30 Anharmonic model problem U(q) = q2 2 + q4 Asymptotic analysis of sampling error of hq 2 i typical, e.g. OBABO: O( t 2 ) ABOBA: BAOAB: O( t 2 ) 9 2 t O( 2 t 4 )
31 The Timestep Problem in MD Find timestepping methods that allow t t Verlet Three possibilities: Use an implicit method Eliminate the motion of fast components constraints Isolate the stiff terms for efficient treatment multiple timestepping
32 Implicit Methods e.g. implicit midpoint (unconditionally stable) Requires solution of a nonlinear system of equations by some iterative procedure, e.g. fixed point iteration, z (k+1) n+1 = z n + h 2 (f(z n)+f(z (k) n+1 )) need to iterate to convergence requires multiple f-evaluations
33 Implicit Methods Even if stable, a method can alter model frequencies At large stepsize, high frequencies become low frequencies! Leads to nonphysical resonances of dynamical modes, e.g. bond stretch directly impacting a dihedral bend
34 Constraints
35 The Timestep Problem in MD Find timestepping methods that allow t t Verlet Three possibilities for deterministic methods Use an implicit method Eliminate the motion of fast components constraints Isolate the stiff terms for efficient treatment multiple timestepping
36 Pendulum stiff harmonic spring rigid rod (holonomic constraint) q = m 1 p tension ṗ = apple 0 g q grav 0=q q 1
37 Constrained Hamiltonian Systems d dt q = M 1 p d dt p = F X m i=1 irg i (q ) 0=g j (q ), j =1, 2,...,m 0=rg j (q ) T M 1 p, j =1, 2,...,m cotangent bundle hidden constraints T M = {(q, p) g j (q) =0, rg j (q ) T M 1 p =0,j =1, 2,...,m} symplectic: dynamics preserves [dq ^ dp] T M
38 SHAKE [Berendsen, Ciccotti, Ryckaert 1977] RATTLE [Andersen 1983] L. & Skeel, J.Comput. Phys, 1994 SHAKE and RATTLE are symplectic methods and are actually the same method (conjugate methods).
39 Microcanonical Sampling Verlet limited to about 2.7fs with stability 1fs with 1 kcal/mol energy accuracy SHAKE (RATTLE) OK up to as much as 4fs with energy accuracy (typical ~2-3fs)
40 No Solution Solution Not Unique Large stepsize relative to curvature large errors (or instability)
41 Geodesic Integrator
42 Geodesic Integrator B.L. and G. Patrick, J. Nonlin. Sci (deterministic) B.L. and C. Matthews, Proc Roy Soc A 2016 (stochastic) An alternative to SHAKE discretization Idea: preserve the configuration manifold during position moves and the cotangent space during impulse. The proper analogue of the Verlet (BAB) g-verlet: B A B h,h T M h/2,ut M h,t T M h/2,u T M Combines: geodesic flow projected kicks
43 Constrained Langevin Dynamics d dt q = M 1 p d p dt p = F p + 2k B T M 1/2 (t) 0=g j (q ), j =1, 2,...,m mx i=1 irg i (q ) 0=rg j (q ) T M 1 p, j =1, 2,...,m T M = {(q, p) g j (q) =0, rg j (q ) T M 1 p =0,j =1, 2,...,m} Some precisions: Lelievre, Rousset and Stoltz, Math. Comp., 2010
44 Geodesic Integrator Numerical Analysis Constrained Hamiltonian system in coordinates Splittings BCH (nastier commutators) BAOAB 2nd order superconvergence
45 Implementation To implement the geodesic integrator, we use a sequence of SHAKE/RATTLE steps for the A step These steps do not require re-evaluation of the force field, so each iteration is relatively cheap. The B and O (for Langevin) steps are simple RATTLE projections, so there is no significant added cost for these. No need for Hessians/normal modes, or tuning/parameterization
46 Box of H2O For rigid body water, implement the geodesic integrator using SETTLE: exact geodesic steps Standard Methods Our Schemes
47 Self-Diffusion Coefficient of TIP3P Water Self-diffusion coeff of TIP3P water in water ~ 4.51x10^-5 Appx 2x the experimental value! 2.34±0.05cm 2 /s
48 Solute-Solvent Splitting The object of bio-md simulation is virtually always a protein or nucleic acid fragment + solvent (water) bath. OnceOn Once the bonds and selected (H-X-H, X-Y-H) angles of the solute and waters are constrained, the next fastest modes are due to flexible angle bonds of the solute. These would limit the stepsize to around 5fs, even with the geodesic integrator.
49 Solute-Solvent Splitting The obvious solution is to break the interaction forces into three pieces, denoted PP (protein-protein), PS (protein-solvent), and SS (solvent-solvent) exp SS dominates the computational cost. PP and PS determine the stepsize. Therefore, consider multiple timestepping: h 2 U SS(q) exp h 2 [T (p)+u PP(q)+U PS (q)] m=2 or m=3 iterations of a geodesic Langevin integrator using k=5 RATTLE substeps exp 1. PP+PS is viewed as a many-body, stochastic system 2. Water motion can be implemented using SETTLE. h 2 U SS(q)
50 Comparisons RATTLE: Standard RATTLE + Langevin TINKER: Tinker s internal scheme [OBABO with constraint projections] MEVME: Scheme of Lelievre, Rousset, Stoltz 2010 MTS-BAOAB: geodesic-mts method, using 20 steps to compute the geodesics and m=3 (PP+PS) steps per step
51 Alanine Dipeptide, Solvated, 300K, Tinker Code illustrative FE surface
52
53 Alanine Dipeptide Effective Free Energy Barrier Height Existing methods
54 Multiple timestepping
55 The Timestep Problem in MD Find timestepping methods that allow t t Verlet Three possibilities for deterministic methods Use an implicit method Eliminate the motion of fast components constraints [~ 8 or 9fs] Isolate the stiff terms for efficient treatment multiple timestepping
56 Multiple Timestepping (RESPA) Tuckerman, Martyna, Berne 90 also Grubmüller, Heller, Windemuth, Schulten 91 Some forces are fast-changing (e.g. short-ranged) but cheap to compute Some forces are slow-changing (e.g. long-ranged) and costly to compute t exp 2 L U slow apple t t exp 2r L U exp fast r L K exp t exp 2 L U slow r t 2r L U fast Many fast cheap evaluations Few slow costly ones
57 Linear Model Problem model problem H = p2 2 + (1 + 2 )q idealized multiple timestepping (fast solve is exact) U S (q) =q 2 /2 U F (q) = 2 q 2 /2 Kick with U S Solve H F = p 2 /2 + U F (q) Kick with U S t/2 t t/2
58 Resonance Eigenvalues of timestep map: 1 2 = 1; = 2 cos( t ) " / t / < 0 sin( t ) instabilities t
59 Resonance For the harmonic model, Verlet introduces a stability restriction t < 2/ Multiple timestepping, in the idealized form described here has a stability restriction of about t< / i.e., not very dramatic improvement. In practice, we are limited to around 3.5fs
60 Mollified Impulse Method Garcia-Archilla, Sanz-Serna and Skeel 1998 To stabilize multiple timestepping, one idea is to average out over an associated fast dynamics to ``mollify the impulse in RESPA starting from e.g (q,0), solve to produce a sequence then define: 0 = q, 1, 2,..., K Advance the step using a mollified impulse: Ũ(q) :=U s (A(q))
61 Pushing the limits SHAKE-MOLLY [Izaguirre, Reich, Skeel 1999] constrains the slow force evaluations to lie exactly at the bond stretch minima (multiple timestepping + constraints). This allows fully flexible MD with larger stepsizes With the very best deterministic schemes, however, we find for a biological molecule: t. 5fs Bob Skeel s license plate ca. 2000: FEMTO5
62 Stochastic Methods Introducing random perturbations might seem to complicate the numerical integrator. In fact, if done correctly, it is possible to gain in two ways: 1) a significant stability improvement is possible 2) substantially higher accuracy is possible by shifting emphasis to the invariant distribution We have already seen this in the behavior of Langevin vis a vis Hamiltonian dynamics. These benefits carry over to multiple timestepping and constraints.
63 Langevin-RESPA L. & Mathews, Molecular Dynamics, Springer 2015 It is natural to think that stochastic perturbations may have a stabilizing influence on multiple timestepping. Consider the combination of Langevin (stochastic) dynamics with RESPA,
64 Fast Dynamics: Combine with kicks by the slow force
65 Stability if (X) apple 1
66 sampling error noise strength Verlet stability threshold timestep Langevin can stabilize RESPA but sampling accuracy is severely compromised except at high friction
67 Configurational Sampling H(q, p) =p T M 1 p/2+u(q) e H = e T (p) e U(q) In MD we typically are interested in q-dependent quantities. The Hamiltonian is used to enhance exploration but ultimately we don t much care about momenta. Question: can we use this freedom, together with stochastic dynamics, to control resonance in multiple timestepping while taking advantage of a known force field decomposition?
68 Stochastic Isokinetic Nosé-Hoover (SIN) L., Margul and Tuckerman, Mol. Phys q = p ṗ = F (q) p 1 = = 2 1 kt 2 + p 2 kt (t) = 2p F K(p, 1 ) to make K(p, 1 ):=p p =kt isokinetic constraint
69 Questions What is the invariant measure of the SIN system? Is the dynamics ergodic? How to numerically integrate the equations? What is the performance in practice?
70 Compressible statistical mechanics M. E. Tuckerman, Y. Liu, G. Ciccotti, and G. J. Martyna, J. Chem. Phys. 115, 1678 (2001). ẋ = f(x) conservation laws: C j (x) =c j,j =1, 2,...,k phase space compressibility: apple(x) =divf(x) and find: apple = L f w Then, if ergodic, the partition function is: Z (N,V, 1 )= e w(x) k j=1 [C j (x) c j ]dx
71 q = p ṗ = F (q) f(x) p 1 = OU 2 = 2 1 kt 2 + p 2 kt (t) apple(x) =divf(x) = 2 p@ /@p /@ 1 2 = L f [U(q)+ 2 2/2] an invariant distribution is given by (x) =e U(x) e 2 2 /2 [K(p, 1 ) kt ]
72 (x) =e U(x) e 2 2 /2 [K(p, 1 ) kt ] Dynamics evolves on the generalized semi-cylinder defined by periodic boundary conditions (q), and the isokinetic constraint on (p, 1 ), with the restriction 1 > 0 We need to demonstrate the Hörmander condition on the manifold, and find a Lyapunov function. Lyapunov function: ( 2 )=1+ 2s 2
73 Hörmander Condition SDE mx ẋ = b 0 (x)+ b i (x) i (t) i=1 b j 2 T M Existence of a unique invariant measure on M is assured if the collection of iterated commutators of the vector fields spans the tangent space at ever point. span{b m i=0, [b i,b j ] m i,j=0,...,} = T M 1 constraint, 4 variables, thus 3-dimensional tangent space
74 Ergodicity Property of SIN Span{ } = T M Thus conclude (subject to a minorization condition, left as an exercise) that SIN is ergodic on M with a unique invariant distribution which, after marginalization reduces to: (x) =e The Alexander Davie commutator U(x)
75 SIN(R) Stochastic Isokinetic Nosé-Hoover (RESPA)
76
77 periodic box, 25A sides, 512 fully flexible SPC water molecules (all atom) periodic boundary conditions three level splitting U = U bond + U s.r. + U l.r. Coulombic forces: Smooth Particle Mesh Ewald Long ranged forces includes reciprocal space part + screened coulomb interactions within the simulation cell Short-ranged Regime: e.g. Lennard-Jones cutoffs to 6A, real space part of Ewald. Bond part: Intramolecular interactions, e.g. angle and length bonds.
78 Flexible H2O t = 99fs t inner =0.5fs, t= 3fs mid radial distributions O-O H-O H-H
79 Comparisons for biomolecules/detailed H2O timestep Speedup Notes typical MD 1fs 1 fully flexible typical Langevin MD 2fs 2 fully flexible BAOAB Langevin 2.7fs 2.7 optimized for config sampling MTS (RESPA) 3.5fs ~3 resonances limit stability SHAKE MD 4fs 3.5 symmetric Newton SHAKE Langevin 4fs 3.5 MEVME, g-obabo, etc. Equilibrium-MOLLY Izaguirre, Reich &Skeel (JCP, 99) 6fs 4+ deterministic MTS-based force decomposition g-baoab-mts 8-9fs 6-7 Langevin, geodesics (RATTLE), PP+PS+SS decomposition Colored noise thermostats to control resonance in MTS Ceriotti et al LN, Langevin normal mode analysis, Barth&Schlick Stochastic-Isokinetic Nose (L., Margul, Tuckerman) 10-12fs?? 48fs 10 ~100fs >10 pre-analysis of fundamental frequencies (Hessians/normal modes) and tuning of parameters requires Hessians+ force decomposition force decomposition (Ewald) q-sampling *only*
80 Next Up In Lecture 3, we say goodbye to biomolecules and look at how thermostats can be extended and enhanced to address other types of problems:
Molecular Dynamics Lecture 3
Molecular Dynamics Lecture 3 Ben Leimkuhler the problem of the timestep in MD constraints - SHAKE/RATTLE multiple timestepping stochastic methods for holonomic constraints stochastic multiple timestepping
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 informationFrom Physical Modelling to Big Data Analytics: Examples and Challenges
From Physical Modelling to Big Data Analytics: Examples and Challenges Ben Leimkuhler University of Edinburgh EPSRC/NSF ERC ATI Summit on Big Data in the Physical Sciences Outline Part I: Challenges Ideas
More informationarxiv: v2 [physics.comp-ph] 24 Apr 2013
Robust and efficient configurational molecular sampling via Langevin Dynamics Benedict Leimkuhler and Charles Matthews School of Mathematics and Maxwell Institute of Mathematical Sciences, James Clerk
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 informationEnhanced sampling via molecular dynamics II: Unbiased approaches
Enhanced sampling via molecular dynamics II: Unbiased approaches Mark E. Tuckerman Dept. of Chemistry and Courant Institute of Mathematical Sciences New York University, 100 Washington Square East, NY
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 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 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 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 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 informationGeometric Numerical Integration
Geometric Numerical Integration (Ernst Hairer, TU München, winter 2009/10) Development of numerical ordinary differential equations Nonstiff differential equations (since about 1850), see [4, 2, 1] Adams
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 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 informationApproach to Thermal Equilibrium in Biomolecular
Approach to Thermal Equilibrium in Biomolecular Simulation Eric Barth 1, Ben Leimkuhler 2, and Chris Sweet 2 1 Department of Mathematics Kalamazoo College Kalamazoo, Michigan, USA 49006 2 Centre for Mathematical
More informationWORD SERIES FOR THE ANALYSIS OF SPLITTING SDE INTEGRATORS. Alfonso Álamo/J. M. Sanz-Serna Universidad de Valladolid/Universidad Carlos III de Madrid
WORD SERIES FOR THE ANALYSIS OF SPLITTING SDE INTEGRATORS Alfonso Álamo/J. M. Sanz-Serna Universidad de Valladolid/Universidad Carlos III de Madrid 1 I. OVERVIEW 2 The importance of splitting integrators
More informationMolecular dynamics simulation of Aquaporin-1. 4 nm
Molecular dynamics simulation of Aquaporin-1 4 nm Molecular Dynamics Simulations Schrödinger equation i~@ t (r, R) =H (r, R) Born-Oppenheimer approximation H e e(r; R) =E e (R) e(r; R) Nucleic motion described
More informationarxiv: v1 [cond-mat.stat-mech] 7 Mar 2019
Langevin thermostat for robust configurational and kinetic sampling Oded Farago, Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB EW, United Kingdom Department of Biomedical
More informationOvercoming stability limitations in biomolecular dynamics. I. Combining force splitting via extrapolation with Langevin dynamics in LN
JOURNAL OF CHEMICAL PHYSICS VOLUME 109, NUMBER 5 1 AUGUST 1998 Overcoming stability limitations in biomolecular dynamics. I. Combining force splitting via extrapolation with Langevin dynamics in LN Eric
More informationFree energy calculations and the potential of mean force
Free energy calculations and the potential of mean force IMA Workshop on Classical and Quantum Approaches in Molecular Modeling Mark Tuckerman Dept. of Chemistry and Courant Institute of Mathematical Science
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 informationAccepted Manuscript. Correcting Mesh-Based Force Calculations to Conserve Both Energy and Momentum in Molecular Dynamics Simulations
Accepted Manuscript Correcting Mesh-Based Force Calculations to Conserve Both Energy and Momentum in Molecular Dynamics Simulations Robert D. Seel, David J. Hardy, James C. Phillips PII: S001-1(0)001-1
More informationGear methods I + 1/18
Gear methods I + 1/18 Predictor-corrector type: knowledge of history is used to predict an approximate solution, which is made more accurate in the following step we do not want (otherwise good) methods
More informationTime-Dependent Statistical Mechanics 5. The classical atomic fluid, classical mechanics, and classical equilibrium statistical mechanics
Time-Dependent Statistical Mechanics 5. The classical atomic fluid, classical mechanics, and classical equilibrium statistical mechanics c Hans C. Andersen October 1, 2009 While we know that in principle
More informationStochastic Variational Integrators
Stochastic Variational Integrators Jeremy Schmitt Variational integrators are numerical geometric integrator s derived from discretizing Hamilton s principle. They are symplectic integrators that exhibit
More 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 informationChaotic motion. Phys 750 Lecture 9
Chaotic motion Phys 750 Lecture 9 Finite-difference equations Finite difference equation approximates a differential equation as an iterative map (x n+1,v n+1 )=M[(x n,v n )] Evolution from time t =0to
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 informationPhysics 106b: Lecture 7 25 January, 2018
Physics 106b: Lecture 7 25 January, 2018 Hamiltonian Chaos: Introduction Integrable Systems We start with systems that do not exhibit chaos, but instead have simple periodic motion (like the SHO) with
More informationAb initio molecular dynamics and nuclear quantum effects
Ab initio molecular dynamics and nuclear quantum effects Luca M. Ghiringhelli Fritz Haber Institute Hands on workshop density functional theory and beyond: First principles simulations of molecules and
More informationA GENTLE STOCHASTIC THERMOSTAT FOR MOLECULAR DYNAMICS
A GENTLE STOCHASTIC THERMOSTAT FOR MOLECULAR DYNAMICS BEN LEIMKUHLER, EMAD NOORIZADEH, AND FLORIAN THEIL Abstract. We discuss a dynamical technique for sampling the canonical measure in molecular dynamics.
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 informationGEOMETRIC INTEGRATION OF ORDINARY DIFFERENTIAL EQUATIONS ON MANIFOLDS
BIT 0006-3835/01/4105-0996 $16.00 2001, Vol. 41, No. 5, pp. 996 1007 c Swets & Zeitlinger GEOMETRIC INTEGRATION OF ORDINARY DIFFERENTIAL EQUATIONS ON MANIFOLDS E. HAIRER Section de mathématiques, Université
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 informationMultiple time step Monte Carlo
JOURNAL OF CHEMICAL PHYSICS VOLUME 117, NUMBER 18 8 NOVEMBER 2002 Multiple time step Monte Carlo Balázs Hetényi a) Department of Chemistry, Princeton University, Princeton, NJ 08544 and Department of Chemistry
More informationMolecular Dynamics Simulations. Dr. Noelia Faginas Lago Dipartimento di Chimica,Biologia e Biotecnologie Università di Perugia
Molecular Dynamics Simulations Dr. Noelia Faginas Lago Dipartimento di Chimica,Biologia e Biotecnologie Università di Perugia 1 An Introduction to Molecular Dynamics Simulations Macroscopic properties
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 informationAdvanced Molecular Molecular Dynamics
Advanced Molecular Molecular Dynamics Technical details May 12, 2014 Integration of harmonic oscillator r m period = 2 k k and the temperature T determine the sampling of x (here T is related with v 0
More informationStatistical methods in atomistic computer simulations
Statistical methods in atomistic computer simulations Prof. Michele Ceriotti, michele.ceriotti@epfl.ch This course gives an overview of simulation techniques that are useful for the computational modeling
More informationEnergy Barriers and Rates - Transition State Theory for Physicists
Energy Barriers and Rates - Transition State Theory for Physicists Daniel C. Elton October 12, 2013 Useful relations 1 cal = 4.184 J 1 kcal mole 1 = 0.0434 ev per particle 1 kj mole 1 = 0.0104 ev per particle
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 information7 To solve numerically the equation of motion, we use the velocity Verlet or leap frog algorithm. _ V i n = F i n m i (F.5) For time step, we approxim
69 Appendix F Molecular Dynamics F. Introduction In this chapter, we deal with the theories and techniques used in molecular dynamics simulation. The fundamental dynamics equations of any system is the
More informationThe Molecular Dynamics Method
H-bond energy (kcal/mol) - 4.0 The Molecular Dynamics Method Fibronectin III_1, a mechanical protein that glues cells together in wound healing and in preventing tumor metastasis 0 ATPase, a molecular
More informationA PROJECTIVE THERMOSTATTING DYNAMICS TECHNIQUE
A PROJECTIVE THERMOSTATTING DYNAMICS TECHNIQUE ZHIDONG JIA AND BENEDICT J. LEIMKUHLER Abstract. A dynamical framework is developed with several variations for modeling multiple timescale molecular dynamics
More informationBACKGROUND IN SYMPLECTIC GEOMETRY
BACKGROUND IN SYMPLECTIC GEOMETRY NILAY KUMAR Today I want to introduce some of the symplectic structure underlying classical mechanics. The key idea is actually quite old and in its various formulations
More informationUsing Molecular Dynamics to Compute Properties CHEM 430
Using Molecular Dynamics to Compute Properties CHEM 43 Heat Capacity and Energy Fluctuations Running an MD Simulation Equilibration Phase Before data-collection and results can be analyzed the system
More informationIntroduction to molecular dynamics
1 Introduction to molecular dynamics Yves Lansac Université François Rabelais, Tours, France Visiting MSE, GIST for the summer Molecular Simulation 2 Molecular simulation is a computational experiment.
More informationChaotic motion. Phys 420/580 Lecture 10
Chaotic motion Phys 420/580 Lecture 10 Finite-difference equations Finite difference equation approximates a differential equation as an iterative map (x n+1,v n+1 )=M[(x n,v n )] Evolution from time t
More informationQUANTUM MARKOVIAN KINETIC EQUATION FOR HARMONIC OSCILLATOR. Boris V. Bondarev
QUANTUM MARKOVIAN KINETIC EQUATION FOR HARMONIC OSCILLATOR Boris V. Bondarev Moscow Aviation Institute, Volokolamsk road, 4, 15871, Moscow, Russia E-mail: bondarev.b@mail.ru arxiv:130.0303v1 [physics.gen-ph]
More informationIntroduction to model potential Molecular Dynamics A3hourcourseatICTP
Introduction to model potential Molecular Dynamics A3hourcourseatICTP Alessandro Mattoni 1 1 Istituto Officina dei Materiali CNR-IOM Unità di Cagliari SLACS ICTP School on numerical methods for energy,
More informationSYMMETRIC PROJECTION METHODS FOR DIFFERENTIAL EQUATIONS ON MANIFOLDS
BIT 0006-3835/00/4004-0726 $15.00 2000, Vol. 40, No. 4, pp. 726 734 c Swets & Zeitlinger SYMMETRIC PROJECTION METHODS FOR DIFFERENTIAL EQUATIONS ON MANIFOLDS E. HAIRER Section de mathématiques, Université
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 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 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 informationGSHMC: An efficient Markov chain Monte Carlo sampling method. Sebastian Reich in collaboration with Elena Akhmatskaya (Fujitsu Laboratories Europe)
GSHMC: An efficient Markov chain Monte Carlo sampling method Sebastian Reich in collaboration with Elena Akhmatskaya (Fujitsu Laboratories Europe) 1. Motivation In the first lecture, we started from a
More informationStochastic-Dynamical Thermostats for Constraints and Stiff Restraints
EPJ manuscript No. (will be inserted by the editor) Stochastic-Dynamical Thermostats for Constraints and Stiff Restraints Janis Bajars 1,a, Jason Frank 1,b, and Ben Leimkuhler 1,c 1 Centrum Wiskunde &
More informationMolecular Dynamics Simulations
Molecular Dynamics Simulations Dr. Kasra Momeni www.knanosys.com Outline Long-range Interactions Ewald Sum Fast Multipole Method Spherically Truncated Coulombic Potential Speeding up Calculations SPaSM
More informationCanonical transformations (Lecture 4)
Canonical transformations (Lecture 4) January 26, 2016 61/441 Lecture outline We will introduce and discuss canonical transformations that conserve the Hamiltonian structure of equations of motion. Poisson
More informationDynamic force matching: Construction of dynamic coarse-grained models with realistic short time dynamics and accurate long time dynamics
for resubmission Dynamic force matching: Construction of dynamic coarse-grained models with realistic short time dynamics and accurate long time dynamics Aram Davtyan, 1 Gregory A. Voth, 1 2, a) and Hans
More informationProblem reduction, memory and renormalization
Problem reduction, memory and renormalization Panos Stinis Northwest Institute for Advanced Computing Pacific Northwest National Laboratory A simple example The linear differential system for x(t) and
More informationEfficiency of Runge-Kutta Methods in Solving Simple Harmonic Oscillators
MATEMATIKA, 8, Volume 3, Number, c Penerbit UTM Press. All rights reserved Efficiency of Runge-Kutta Methods in Solving Simple Harmonic Oscillators Annie Gorgey and Nor Azian Aini Mat Department of Mathematics,
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 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 informationarxiv: v1 [physics.chem-ph] 24 Apr 2018
arxiv:1804.08913v1 [physics.chem-ph] 24 Apr 2018 Fast-Forward Langevin Dynamics with Momentum Flips Mahdi Hijazi, 1 David M. Wilkins, 1, a) and Michele Ceriotti 1 Laboratory of Computational Science and
More information4. The Green Kubo Relations
4. The Green Kubo Relations 4.1 The Langevin Equation In 1828 the botanist Robert Brown observed the motion of pollen grains suspended in a fluid. Although the system was allowed to come to equilibrium,
More 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 informationThe Molecular Dynamics Method
The Molecular Dynamics Method Thermal motion of a lipid bilayer Water permeation through channels Selective sugar transport Potential Energy (hyper)surface What is Force? Energy U(x) F = d dx U(x) Conformation
More informationTHE Hamilton equations of motion constitute a system
Proceedings of the World Congress on Engineering 0 Vol I WCE 0, July 4-6, 0, London, U.K. Systematic Improvement of Splitting Methods for the Hamilton Equations Asif Mushtaq, Anne Kværnø, and Kåre Olaussen
More informationMolecular Dynamics Simulation of a Nanoconfined Water Film
Molecular Dynamics Simulation of a Nanoconfined Water Film Kyle Lindquist, Shu-Han Chao May 7, 2013 1 Introduction The behavior of water confined in nano-scale environment is of interest in many applications.
More informationNumerical Integration of Equations of Motion
GraSMech course 2009-2010 Computer-aided analysis of rigid and flexible multibody systems Numerical Integration of Equations of Motion Prof. Olivier Verlinden (FPMs) Olivier.Verlinden@fpms.ac.be Prof.
More informationSampling the free energy surfaces of collective variables
Sampling the free energy surfaces of collective variables Jérôme Hénin Enhanced Sampling and Free-Energy Calculations Urbana, 12 September 2018 Please interrupt! struct bioinform phys chem theoretical
More informationin order to insure that the Liouville equation for f(?; t) is still valid. These equations of motion will give rise to a distribution function f(?; t)
G25.2651: Statistical Mechanics Notes for Lecture 21 Consider Hamilton's equations in the form I. CLASSICAL LINEAR RESPONSE THEORY _q i = @H @p i _p i =? @H @q i We noted early in the course that an ensemble
More informationChapter 1 Symplectic Integrator and Beam Dynamics Simulations
Chapter 1 and Beam Accelerator Physics Group, Journal Club November 9, 2010 and Beam NSLS-II Brookhaven National Laboratory 1.1 (70) Big Picture for Numerical Accelerator Physics and Beam For single particle
More informationWhy study protein dynamics?
Why study protein dynamics? Protein flexibility is crucial for function. One average structure is not enough. Proteins constantly sample configurational space. Transport - binding and moving molecules
More informationIMPLICIT MASS-MATRIX PENALIZATION OF HAMILTONIAN DYNAMICS WITH APPLICATION TO EXACT SAMPLING OF STIFF SYSTEMS
IMPLICIT MASS-MATRIX PENALIZATION OF HAMILTONIAN DYNAMICS WITH APPLICATION TO EXACT SAMPLING OF STIFF SYSTEMS PETR PLECHÁČ AND MATHIAS ROUSSET Résumé. An implicit mass-matrix penalization (IMMP) of Hamiltonian
More informationCitation for published version (APA): Hess, B. (2002). Stochastic concepts in molecular simulation Groningen: s.n.
University of Groningen Stochastic concepts in molecular simulation Hess, Berk IMPORTANT NOTE: You are advised to consult the publisher's version (publisher's PDF) if you wish to cite from it. Please check
More informationSupporting Information Soft Nanoparticles: Nano Ionic Networks of Associated Ionic Polymers
Electronic Supplementary Material (ESI) for Nanoscale. This journal is The Royal Society of Chemistry 2016 Supporting Information Soft Nanoparticles: Nano Ionic Networks of Associated Ionic Polymers Dipak
More informationAll-atom Molecular Mechanics. Trent E. Balius AMS 535 / CHE /27/2010
All-atom Molecular Mechanics Trent E. Balius AMS 535 / CHE 535 09/27/2010 Outline Molecular models Molecular mechanics Force Fields Potential energy function functional form parameters and parameterization
More informationExercise 2: Solvating the Structure Before you continue, follow these steps: Setting up Periodic Boundary Conditions
Exercise 2: Solvating the Structure HyperChem lets you place a molecular system in a periodic box of water molecules to simulate behavior in aqueous solution, as in a biological system. In this exercise,
More informationOn the controllability of conservative systems
On the controllability of conservative systems Martin Hairer Department of Mathematics, The University of Warwick, Coventry CV4 7AL, United Kingdom. Email: hairer@maths.warwick.co.uk Abstract We consider
More informationNosé-Hoover chain method for nonequilibrium molecular dynamics simulation
PHYSICAL REVIEW E VOLUME 61, NUMBER 5 MAY 000 Nosé-Hoover chain method for nonequilibrium molecular dynamics simulation A. C. Brańka Institute of Molecular Physics, Polish Academy of Sciences, Smoluchowskiego
More informationReducing round-off errors in symmetric multistep methods
Reducing round-off errors in symmetric multistep methods Paola Console a, Ernst Hairer a a Section de Mathématiques, Université de Genève, 2-4 rue du Lièvre, CH-1211 Genève 4, Switzerland. (Paola.Console@unige.ch,
More informationImplicit-explicit variational integration of highly oscillatory problems
Implicit-explicit variational integration of igly oscillatory problems Ari Stern Structured Integrators Worksop April 9, 9 Stern, A., and E. Grinspun. Multiscale Model. Simul., to appear. arxiv:88.39 [mat.na].
More informationPotential Energy (hyper)surface
The Molecular Dynamics Method Thermal motion of a lipid bilayer Water permeation through channels Selective sugar transport Potential Energy (hyper)surface What is Force? Energy U(x) F = " d dx U(x) Conformation
More informationThe Effect of Model Internal Flexibility Upon NEMD Simulations of Viscosity
Draft: September 29, 1999 The Effect of Model Internal Flexibility Upon NEMD Simulations of Viscosity N. G. Fuller 1 and R. L. Rowley 1,2 Abstract The influence of model flexibility upon simulated viscosity
More informationControlling Mechanical Systems by Active Constraints. Alberto Bressan. Department of Mathematics, Penn State University
Controlling Mechanical Systems by Active Constraints Alberto Bressan Department of Mathematics, Penn State University 1 Control of Mechanical Systems: Two approaches by applying external forces by directly
More informationNUMERICAL METHODS FOR SECOND-ORDER STOCHASTIC DIFFERENTIAL EQUATIONS
SIAM J SCI COMPUT Vol 9, No, pp 45 64 c 7 Society for Industrial and Applied Mathematics NUMERICAL METHODS FOR SECOND-ORDER STOCHASTIC DIFFERENTIAL EQUATIONS KEVIN BURRAGE, IAN LENANE, AND GRANT LYTHE
More informationDeclaration. I, Derrick Beckedahl, declare that
Abstract In this dissertation we reformulate the configurational temperature Nosè-Hoover thermostat proposed by Braga and Travis (C. Braga and K.P. Travis, 005), using the antisymmetric matrix formalism
More informationNumerical Methods for ODEs. Lectures for PSU Summer Programs Xiantao Li
Numerical Methods for ODEs Lectures for PSU Summer Programs Xiantao Li Outline Introduction Some Challenges Numerical methods for ODEs Stiff ODEs Accuracy Constrained dynamics Stability Coarse-graining
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 informationfor changing independent variables. Most simply for a function f(x) the Legendre transformation f(x) B(s) takes the form B(s) = xs f(x) with s = df
Physics 106a, Caltech 1 November, 2018 Lecture 10: Hamiltonian Mechanics I The Hamiltonian In the Hamiltonian formulation of dynamics each second order ODE given by the Euler- Lagrange equation in terms
More informationarxiv: v1 [cond-mat.stat-mech] 28 Mar 2008
Canonical sampling through velocity-rescaling Giovanni Bussi, Davide Donadio, and Michele Parrinello Computational Science, Department of Chemistry and Applied Biosciences, ETH Zürich, USI Campus, Via
More informationOptimum protocol for fast-switching free-energy calculations
PHYSICAL REVIEW E 8, 7 Optimum protocol for fast-switching free-energy calculations Philipp Geiger and Christoph Dellago* Faculty of Physics, University of Vienna, Boltzmanngasse 5, 9 Vienna, Austria Received
More informationPhysics 115/242 Comparison of methods for integrating the simple harmonic oscillator.
Physics 115/4 Comparison of methods for integrating the simple harmonic oscillator. Peter Young I. THE SIMPLE HARMONIC OSCILLATOR The energy (sometimes called the Hamiltonian ) of the simple harmonic oscillator
More informationExploring the energy landscape
Exploring the energy landscape ChE210D Today's lecture: what are general features of the potential energy surface and how can we locate and characterize minima on it Derivatives of the potential energy
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 informationHamiltonian. March 30, 2013
Hamiltonian March 3, 213 Contents 1 Variational problem as a constrained problem 1 1.1 Differential constaint......................... 1 1.2 Canonic form............................. 2 1.3 Hamiltonian..............................
More informationExponential integrators for oscillatory second-order differential equations
Exponential integrators for oscillatory second-order differential equations Marlis Hochbruck and Volker Grimm Mathematisches Institut Heinrich Heine Universität Düsseldorf Cambridge, March 2007 Outline
More informationA Conjecture about Molecular Dynamics
A Conjecture about Molecular Dynamics P. F. Tupper June 12, 2007 Abstract An open problem in numerical analysis is to explain why molecular dynamics works. The difficulty is that numerical trajectories
More informationLecture 11: Potential Energy Functions
Lecture 11: Potential Energy Functions Dr. Ronald M. Levy ronlevy@temple.edu Originally contributed by Lauren Wickstrom (2011) Microscopic/Macroscopic Connection The connection between microscopic interactions
More information