Systematic Coarse-Graining of Molecular Models by the Inverse Monte Carlo and Newton Inversion
|
|
- Lisa Taylor
- 5 years ago
- Views:
Transcription
1 Systematic Coarse-Graining of Molecular Models by the Inverse Monte Carlo and Newton Inversion Alexander Lyubartsev Division of Physical Chemistry Department of Material and Environmental Chemistry Stockholm University NSASM10 Dreseden September 2010
2 Computer modeling from 1st principles to mesoscale Three levels of molecular modeling: First-principles (Quantum Mechanics) Classical Molecular Mechanics / Dynamics nuclei electrons nano scale tape measure 1 nm = 10-9 m atoms coarse-grained molecules atoms molecules atoms molecules 0.1 nm Meso-scale soft matter 1.0 nm small molecules 10 nm soft matter 100 nm large biomolecules nm organelles, bacteria The problem: Larger scale more approximations
3 Coares-graining: reduction degrees of freedom All-atom model 118 atoms Coarse-grained model 10 sites Large-scale simulations We need to: 1) Keep only important degrees of freedom and remove unimportant (solvent; local details of large molecules) 2) For important degrees of freedom we need interaction potential Question: what is the interaction potential for the coarse-grained model?
4 Formal solution: N-body mean force potential Hamiltonian (potential energy) of the original (FG) system: H(R1,R2,...RN, r1,r2,...,rn) important unimportant Ri can be combinations of original degrees of freedom (e.g. centres of mass) Partition function: Z = dr1... dr N dr 1... dr n exp H R1,..., R N, r 1,... r n = dr1... drn exp βh CG R1,..., R N 1 kt 1 CG H R 1,..., R N = ln dr 1... dr n exp βh R 1,..., R N,r 1,..., r n where = is the effective N-body coarse-grained potential = potential of mean force = free energy of the non-important degrees of freedom
5 HCG(R1,...RN) : provides the same structural properties for the important (CG) degrees of freedom R1,...RN as the FG system. Even thermodynamics can be restored (but: HCG(R1,...RN) is temperature and density dependent ). Problem with HCG(R1,...RN) : simulation with N-body potential is unrealistic. We need something usable e.g. pair potentials! H CG R 1, R2,..., R N V ij R ij i j Criteria for approximation: microscopical : - best fit for the energy - best fit for the force (force matching) ensemble averages (observables) - RDF (can be fit exactly) - other ensemble properties Rij = R i R j
6 Newton Inversion (Inverse Monte Carlo) Model Interaction potential direct inverse Properties radial distribution functions, E, P,... Effective potentials for coarse-grained models from "lower level" simulations (ab-initio atomistic ; atomistic coarse grained;) Reconstruct interaction potential from experimental RDF or experimental thermodynamic properties Use both data of fine grained simulations and experiment
7 The method (A.P.Lyubartsev, A.Mirzoev, L.J. Chen, A.Laaksonen, Faraday Discussions, 144, 2010) Assume: H(λ1,λ2,...λM): Hamiltonian (potential energy) depending on a set of parameters λ1,λ2,...λm. We wish to reproduce M canonical averages <S1>,<S2>,...,<SM> Set of λα, α=1,,m Space of Hamiltonians direct {<Sα>} {λα} inverse
8 In the vicinity of an arbitrary point in the space of Hamiltonians one can write: Δ S α = γ S α λγ Δλ γ O Δλ 2 where S S = λγ H H S γ S γ λα λα β = 1/kT This allows us to solve the inverse problem iteratively
9 Choose trial values λα(0) Direct MC or MD Calculate <Sα> (n) and differences <Sα>(n) = <Sα >(n) - Sα* Solve linear equations system Obtain λα(n) Repeat until convergence Algorithm: New Hamiltonian: λα(n1) =λα(n) λα(n) Newton method of solving non-linear equation: S S(λ) S λ S* λ* λ 1 λ0 λ if no convergence, a regularization procedure can be applied: <Sα>(n) = a(<sα >(n) - Sα*) with a<1
10 Reconstruction of interaction parameters from thermodynamic data LJ potential σ, ε direct inverse Pressure P, evaporation heat Ev Derivatives in the matrix: formulas like : 2 2 P 3 6 PV NkT 8 E V P P = σ σ V kt V PE P E E 3 = NkT PV σ σ kt
11 Toy example: Single site LJ water model Can we fit water molecule by a particle with LJ potential? Experimentally: at T=298K, ρ = 1g/cm3 P=1 bar E = kj/mol Start NI from LJ parameters of oxygen in SPC water... ε (kj/mol) 0.65 P(bar) E(kJ/mol) 0.2 σ (Å) Iteration Reg.param 1 But no use of this, with σ=3.121 Å and ε=5.192 kj/mol this is a solid...
12 Inverse Monte Carlo: Reconstruction of pair potentials from RDFs A.Lyubartsev, A.Laaksonen,Phys,Rev.E, 52, 1995) Consider Hamiltonian with pair interaction: H= V r ij Vα i, j Make grid approximation : H= V α S α α α=1,,m Here: Vα=V(Rcutα/M) - potential within α-interval = λ- parameters Sα - number of particle s pairs with distance between them within α-interval = 1 V g r = Sα estimator of RDF: α 2 2 4πr α Δ r N /2 Inversion matrix became: S = S S S S V So: we can reconstruct pair potential from RDF Rcut
13 Some comments Solution of the inverse problem is unique for pair potentials reconstructed from the corresponding RDFs gik(r) Vik(r)const A simpler scheme to correct the potential: V(n1)(r) = V(n)(r) kt ln(g(n)(r)/gref(r)) (A.K.Soper, Chem.Phys.Lett, 202, 295 (1996)) - also known as Iterative Boltzmann inversion May not converge in multicomponent case; May be reasonable to use in the beginning of the iteration process Intramolecular potentials (bonds, angles, torsions) can be easily included Origin: Renormalization Group method
14 Single site water model: reconstructing from atomistic RDF O-O RDF Potential Problems: P 9000 bar ; not a liquid state,...
15 Let us generate g(r) at lower density, e.g g/cm3 OO RDF at different densities; for 500 water molecules
16 Effective potentials generated from different densities Potential obtained in inhomogeneous system: keep correct RDF and provides realistic pressure RDFs at density ρ=1 g/cm3
17 Computer modeling from 1st principles to mesoscale Three levels of molecular modeling: First-principles (Quantum Mechanics) Classical Molecular Mechanics / Dynamics nuclei electrons nano scale tape measure 1 nm = 10-9 m atoms coarse-grained molecules atoms molecules atoms molecules 0.1 nm Meso-scale soft matter 1.0 nm small molecules 10 nm soft matter 100 nm large biomolecules nm organelles, bacteria
18 From ab-initio to atomistic: ab-initio derived 3-site water model Reference simulation: ab-initio (CPMD) simulation of 32 water 3.0 Box sixe 9.86 Å, BLYP functional Vanderbilt pseudopotentials 25 Ha cutoff timestep 0.15 fs 25 ps simulation at 300 C (Chem.Phys.Lett, 325, 15 (2000)) O-O CPMD Exp SPC 2.5 RDF O-H H-H O-O H-H O-H r (Å) Ab-initio eff. pot SPC model 200 Eff. potential / kt Some details: r (Å) 6
19 Some insight Average energy is -25 kj/mol (flexible model ) Experimental: -41 kj/mol. E(32 H2O) 32 E(H2O,opt) = -29 kj/mol CPMD energy Comparison with CPMD: Difference about 4 kj/mol due to many body effects (?) effective potential energy
20 Li - water ab-initio potentials Input RDF-s obtained from 520 ps Car-Parrinello MD for 1 Li ion in 32 water molecules (A.P.Lyubartsev, K.Laasonen and A.Laaksonen, J.Chem.Phys., 114,3120,2001) IMC procedure: SPC model for water; only LiO potential varied to match ab-initio LiO RDF 60 Inversion RDF Fit Dang, J.Chem.Phys.96,6970(1992) Periole et al, J.Phys.Chem.B,102,8579(1998) Heinzinger, Physica BC, 131, 196 (1985) Eff.pot (kj/m) 40 Non-electrostatic part of the LiO potential may be well fitted by: ULiO = A exp(-br) where A = kj/m b = 3.63 Å-1 r(å) Atomistic simulations: Egorov et al, J.Phys.Chem. B, 107, 3234(2003)
21 From atomistic to coarse-grained: solvent-mediated potentials. Atomistic: All-atom simulations (MD) with explicit solvent (water). State of art: < atoms - box size ~ Å Time scale problem CG: Coarse-grain solutes; use continuum solvent (McMillan-Mayer level), From 60-ties: primitive electrolyte model: ions - hard spheres interacting by Coulombic potential with suitable ε, ion radius - adjustable parameter We can now build effective solvent-mediated potential, which takes into account molecular structure of the solvent. RDF from atomistic MD Implicit solvent effective potentials
22 NaCl ion solution Ion-ion RDFs (from atomistic MD) Ion-ion effective potentials 4,0 6 3,0 RDF 2,5 2,0 1,5 1,0 0,5 0,0 NaCl, L=39Å NaCl, L=24Å NaNa, L=39Å NaNa, L=24Å ClCl, L=39Å ClCl, L=24Å Prim. model 4 Eff. Potential / kt NaCl, L=39Å NaCl, L=24Å NaNa, L=39Å NaNa, L=24Å ClCl, L=39Å ClCl, L=24Å 3, r (Å) 10 5 r (Å) 10 15
23 NaCl osmotic and activity coefficients Solvent-mediated effective potentials were applied to calculate osmotic and activity coefficients of Na and Cl- ions in the whole concentration range. MC simulations are carried out for 200 ion pairs using effective potentials Osmotic coefficient: osm. coef activity coef P osm V 1 F / c = T = kt c NkT Activity coefficient: μ ex λ=exp kt E Concentration (M) Lines are calculated values and points are experimental data Important to have water model with correct permittivity 1
24 Temperature dependence increase T U short r =U r qi q j 4 0 r ij
25 From tail asymptotic extract dielectric permitttivity
26 Short range effective potential Temperature dependence can be effectively included in ε
27 Ion-DNA solvent-mediated potentials All-atom model: Coarse-grained model Na
28 Atomistic MD RDFs for coarse-grained sites Ion - DNA effective potentials Ion - P Ion - C4 (base) 6 Li 6 Na Na * K 2 Cs Li Na Eff. potential / kt Eff. potential / kt 4 Li 4 Na * K* Cs 2 0 Na Eff. potential / kt 6 Ion - C4 (sugar) 4 Na * K* Cs r (Å) r (Å) r (Å)
29 MC simulation: a bigger DNA fragment (3 turns) in a box 100x100x102Å, ions interacting by effective solvent-mediated potentials; no explicit water. These are results for the density profile and integral charge 1.0 Na NA * Integral charge Density profile (M/l) Li K* Cs PB Li Na r (Å) Na * 0.4 K* Cs PB r (Å) 40 50
30 Relative binding affinities of ions The order of relative binding affinities of alkali counterions to DNA, defined by MC simulation with effective potentials, is: Cs > Li > Na > K The binding order was defined also in a number of experimental works: P.D.Ross, R.L.Scruggs, Biopolymers, 2, 89 (1964) ; Electrophoresis: Li>Na>K U.P.Strauss, C.Helfgott, H.Pink, J.Phys.Chem.,71,2550 (1967); Donnan equilibrium: Li>Na>K S.Hallon et al, Biochemistry, 14, 1648 (1975); Circular dichroism Cs>Li>K>Na P.Anderson, W.Bauer, Biochemistry, 17, 594 (1978), DNA supercoiling Cs>Li>K>Na M.L.Bleam, C.F.Anderson, M.T.Record, Proc. Natl.Acad.Sci USA,77,3085 (1980), NMR: I.A.Kuznetsov et al, Reactive Polymers, 3, 37 (1984), Ion exchange Cs>Li>K>Na Li>K Na Qualitative agreement with results of experiments of very different nature.
31 Proline in DMSO (A.P.Lyubartsev, A.Mirzoev, L.J. Chen, A.Laaksonen, Faraday Discussions, 144, 2010) D proline L-proline DMSO solvent Atomistic CG implicit
32 Experimentally: prolines in DMSO self-assemble into different structures, depending on ratio R- and D- forms. This can be also observed in atomistic simulations: but these are tedious slow Still, RDF can be computed and converted to effective potentials for CG model
33 Spatial distribution functions: From atomistic MD From CG simulations
34 Coarse-grained lipid model Lipid bilayer simulations: atoms Possible, but... Coarse-graining:
35 Atomistic MD: 16 DMPC lipids in 1600 H2O CHARMM 27 ff 25 ns simulations
36 R D F calculations N P 4 different groups -> 10 pairs 10 RDFs and eff. intermolecular potentials 5 bond potentials (N-P, P-CO, CO CO, CO-C, C-C) CO C 7 NP 6 6 solid - 16 lipids dashed - 64 lipids PP 3 RDF RDF 5 NN r(a) CC 2 PC NC 5 10 r(a) 15 20
37 Effective potentials: N - CO N - CH P - CO P - CH NN 10 Veff (kj/m) Veff (kj/m) 15 PP r(a) r (A) 30 CO - CO CO - CH CH - CH 10 Veff (kj/m) Veff (kj/m) Bond potentials NP N-P P - CO CO - CH CH - CH r(a) r(a) 6 7 8
38 Bilayer structure: Z Periodic box Density distribution Coarse-grained MC All-atom MD
39 Vesicle formation Start from a square plain piece of membrane, 325x325 Å, 3592 lipids: cut plane
40 Membrane self-assembly MD simulation of 392 CG lipids with Lowe-Andersen thermostat
41 Larger systems (1000 or more lipids) can form vesicles: spherical vesicle (cut in the middle) formed by self-assembly of 1000 DMPC lipids Sometimes: double bicell d-spacing 60-65Å
42 Formation of multilammelar vesicles Start from random configuration, 5000 lipids in 300Å After After 200 ns After After ns ns Phase separation on multilammelar phase with d-spacing ~ 60Å (20 Å water layer between lipids = 30 H2O per lipid) and bulk phase (here nothing) 60-62Å
43 Conclusions The multiscale approach based on the inversion of radial distribution functions provides a straightforward way to build effective potentials starting from ab-inition level to atomistic models and further to coarsegrained models, Effective potentials as good as underlying detailed model (or worse) Effective potentials derived from inhomogeneous setup seems to provide better transferability Use of Newton invertion combining fitting to RDF from detailed simulations and to experimentally known thermodynamical propeties may be an option (to do).
44 Acknowledgements Aatto Laaksonen Alexander Mirzoev Andrei Egorov LiJun Chen $$: Vetenskapsrådet CPU: SNIC resources at PDC, NSC, HPC2N
Multiscale Simulations of DNA and Nucleosome Core Particles
Multiscale Simulations of DNA and Nucleosome Core Particles Alexander Lyubartsev Department of Materials and Environmental Chemistry, Stockholm University Molecular Mobility and Order in Polymer Systems
More informationSystematic Coarse-Graining and Concurrent Multiresolution Simulation of Molecular Liquids
Systematic Coarse-Graining and Concurrent Multiresolution Simulation of Molecular Liquids Cameron F. Abrams Department of Chemical and Biological Engineering Drexel University Philadelphia, PA USA 9 June
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 informationReproducible Estimation of Osmotic Coefficients Using the Inverse Monte Carlo Method
Reproducible Estimation of Osmotic Coefficients Using the Inverse Monte Carlo Method Jaydeep P. Bardhan, 1 Gary K. Leaf, 2 and Bob Eisenberg 1, 2 1 Department of Molecular Biophysics and Physiology Rush
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 informationMichael W. Mahoney Department of Physics, Yale University, New Haven, Connecticut 06520
JOURNAL OF CHEMICAL PHYSICS VOLUME 115, NUMBER 23 15 DECEMBER 2001 Quantum, intramolecular flexibility, and polarizability effects on the reproduction of the density anomaly of liquid water by simple potential
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 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 informationScientific Computing II
Scientific Computing II Molecular Dynamics Simulation Michael Bader SCCS Summer Term 2015 Molecular Dynamics Simulation, Summer Term 2015 1 Continuum Mechanics for Fluid Mechanics? Molecular Dynamics the
More informationCE 530 Molecular Simulation
1 CE 530 Molecular Simulation Lecture 14 Molecular Models David A. Kofke Department of Chemical Engineering SUNY Buffalo kofke@eng.buffalo.edu 2 Review Monte Carlo ensemble averaging, no dynamics easy
More informationMolecular Dynamics. A very brief introduction
Molecular Dynamics A very brief introduction Sander Pronk Dept. of Theoretical Physics KTH Royal Institute of Technology & Science For Life Laboratory Stockholm, Sweden Why computer simulations? Two primary
More informationA Nobel Prize for Molecular Dynamics and QM/MM What is Classical Molecular Dynamics? Simulation of explicit particles (atoms, ions,... ) Particles interact via relatively simple analytical potential
More informationThis semester. Books
Models mostly proteins from detailed to more abstract models Some simulation methods This semester Books None necessary for my group and Prof Rarey Molecular Modelling: Principles and Applications Leach,
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 informationParameterization of Coarse-Grained Molecular Interactions through Potential of Mean Force Calculations and Cluster Expansion Techniques
entropy Article Parameterization of Coarse-Grained Molecular Interactions through Potential of Mean Force Calculations and Cluster Expansion Techniques Anastasios Tsourtis, *, Vagelis Harmandaris,, * and
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 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 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 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 informationMolecular Interactions F14NMI. Lecture 4: worked answers to practice questions
Molecular Interactions F14NMI Lecture 4: worked answers to practice questions http://comp.chem.nottingham.ac.uk/teaching/f14nmi jonathan.hirst@nottingham.ac.uk (1) (a) Describe the Monte Carlo algorithm
More informationCE 530 Molecular Simulation
1 CE 530 Molecular Simulation Lecture 1 David A. Kofke Department of Chemical Engineering SUNY Buffalo kofke@eng.buffalo.edu 2 Time/s Multi-Scale Modeling Based on SDSC Blue Horizon (SP3) 1.728 Tflops
More informationFree energy calculations
Free energy calculations Berk Hess May 5, 2017 Why do free energy calculations? The free energy G gives the population of states: ( ) P 1 G = exp, G = G 2 G 1 P 2 k B T Since we mostly simulate in the
More informationA Molecular Dynamics Simulation of a Homogeneous Organic-Inorganic Hybrid Silica Membrane
A Molecular Dynamics Simulation of a Homogeneous Organic-Inorganic Hybrid Silica Membrane Supplementary Information: Simulation Procedure and Physical Property Analysis Simulation Procedure The molecular
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 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 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 informationStatistical Theory and Learning from Molecular Simulations
Statistical Theory and Learning from Molecular Simulations Lawrence R. Pratt 1 and Susan B. Rempe 2 1 Department of Chemical & Biomolecular Engineering 2 Center for Biological and Material Sciences, Sandia
More informationSpecific Ion Solvtion in Ethylene Carbonate and Propylene Carbonate
Specific Ion Solvtion in Ethylene Carbonate and Propylene Carbonate A. Arslanargin, A. Powers, S. Rick, T. Pollard, T. Beck Univ Cincinnati Chemistry Support: NSF, OSC TSRC 2016 November 2, 2016 A. Arslanargin,
More informationAn introduction to Molecular Dynamics. EMBO, June 2016
An introduction to Molecular Dynamics EMBO, June 2016 What is MD? everything that living things do can be understood in terms of the jiggling and wiggling of atoms. The Feynman Lectures in Physics vol.
More informationLecture 2+3: Simulations of Soft Matter. 1. Why Lecture 1 was irrelevant 2. Coarse graining 3. Phase equilibria 4. Applications
Lecture 2+3: Simulations of Soft Matter 1. Why Lecture 1 was irrelevant 2. Coarse graining 3. Phase equilibria 4. Applications D. Frenkel, Boulder, July 6, 2006 What distinguishes Colloids from atoms or
More informationTheoretical models for the solvent effect
Theoretical models for the solvent effect Benedetta Mennucci Dipartimento di Chimica e Chimica Industriale Web: http://benedetta.dcci.unipi.it Email: bene@dcci.unipi.it 8. Excited electronic states in
More informationDeveloping Monovalent Ion Parameters for the Optimal Point Charge (OPC) Water Model. John Dood Hope College
Developing Monovalent Ion Parameters for the Optimal Point Charge (OPC) Water Model John Dood Hope College What are MD simulations? Model and predict the structure and dynamics of large macromolecules.
More informationMolecular modeling and simulation in chemistry
Molecular modeling and simulation in chemistry 1/20 Elements of modeling: Nuclei + electrons + photons: QED Nuclei + electrons: Schrödinger equation Auxiliary interaction centers (sites) Atoms classical
More informationHydrogels in charged solvents
Hydrogels in charged solvents Peter Košovan, Christian Holm, Tobias Richter 27. May 2013 Institute for Computational Physics Pfaffenwaldring 27 D-70569 Stuttgart Germany Tobias Richter (ICP, Stuttgart)
More informationCoarse-grained force field for simulating polymer-tethered silsesquioxane self-assembly in solution
THE JOURNAL OF CHEMICAL PHYSICS 127, 114102 2007 Coarse-grained force field for simulating polymer-tethered silsesquioxane self-assembly in solution Elaine R. Chan a Department of Chemical Engineering,
More informationWater models in classical simulations
Water models in classical simulations Maria Fyta Institut für Computerphysik, Universität Stuttgart Stuttgart, Germany Water transparent, odorless, tasteless and ubiquitous really simple: two H atoms attached
More informationAdaptive resolution molecular-dynamics simulation: Changing the degrees of freedom on the fly
THE JOURNAL OF CHEMICAL PHYSICS 123, 224106 2005 Adaptive resolution molecular-dynamics simulation: Changing the degrees of freedom on the fly Matej Praprotnik, a Luigi Delle Site, and Kurt Kremer b Max-Planck-Institut
More informationAndré Schleife Department of Materials Science and Engineering
André Schleife Department of Materials Science and Engineering Length Scales (c) ICAMS: http://www.icams.de/cms/upload/01_home/01_research_at_icams/length_scales_1024x780.png Goals for today: Background
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 informationElectrostatics of membrane adhesion
Electrostatics of membrane adhesion S. Marcelja Department of Applied Mathematics, Research School of Physical Sciences and Engineering, The Australian National University, Canberra ACT 6, Australia ABSTRACT
More informationOrigins of Mechanical and Rheological Properties of Polymer Nanocomposites. Venkat Ganesan
Department of Chemical Engineering University of Texas@Austin Origins of Mechanical and Rheological Properties of Polymer Nanocomposites Venkat Ganesan $$$: NSF DMR, Welch Foundation Megha Surve, Victor
More informationUnderstanding Molecular Simulation 2009 Monte Carlo and Molecular Dynamics in different ensembles. Srikanth Sastry
JNCASR August 20, 21 2009 Understanding Molecular Simulation 2009 Monte Carlo and Molecular Dynamics in different ensembles Srikanth Sastry Jawaharlal Nehru Centre for Advanced Scientific Research, Bangalore
More informationComputer simulation methods (1) Dr. Vania Calandrini
Computer simulation methods (1) Dr. Vania Calandrini Why computational methods To understand and predict the properties of complex systems (many degrees of freedom): liquids, solids, adsorption of molecules
More informationComputer simulation methods (2) Dr. Vania Calandrini
Computer simulation methods (2) Dr. Vania Calandrini in the previous lecture: time average versus ensemble average MC versus MD simulations equipartition theorem (=> computing T) virial theorem (=> computing
More informationBrief Review of Statistical Mechanics
Brief Review of Statistical Mechanics Introduction Statistical mechanics: a branch of physics which studies macroscopic systems from a microscopic or molecular point of view (McQuarrie,1976) Also see (Hill,1986;
More informationChapter 2 Experimental sources of intermolecular potentials
Chapter 2 Experimental sources of intermolecular potentials 2.1 Overview thermodynamical properties: heat of vaporization (Trouton s rule) crystal structures ionic crystals rare gas solids physico-chemical
More informationTitle Super- and subcritical hydration of Thermodynamics of hydration Author(s) Matubayasi, N; Nakahara, M Citation JOURNAL OF CHEMICAL PHYSICS (2000), 8109 Issue Date 2000-05-08 URL http://hdl.handle.net/2433/50350
More informationAb Ini'o Molecular Dynamics (MD) Simula?ons
Ab Ini'o Molecular Dynamics (MD) Simula?ons Rick Remsing ICMS, CCDM, Temple University, Philadelphia, PA What are Molecular Dynamics (MD) Simulations? Technique to compute statistical and transport properties
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 information3. Solutions W = N!/(N A!N B!) (3.1) Using Stirling s approximation ln(n!) = NlnN N: ΔS mix = k (N A lnn + N B lnn N A lnn A N B lnn B ) (3.
3. Solutions Many biological processes occur between molecules in aqueous solution. In addition, many protein and nucleic acid molecules adopt three-dimensional structure ( fold ) in aqueous solution.
More informationComputing free energy: Thermodynamic perturbation and beyond
Computing free energy: Thermodynamic perturbation and beyond Extending the scale Length (m) 1 10 3 Potential Energy Surface: {Ri} 10 6 (3N+1) dimensional 10 9 E Thermodynamics: p, T, V, N continuum ls
More informationStructure of Cement Phases from ab initio Modeling Crystalline C-S-HC
Structure of Cement Phases from ab initio Modeling Crystalline C-S-HC Sergey V. Churakov sergey.churakov@psi.ch Paul Scherrer Institute Switzerland Cement Phase Composition C-S-H H Solid Solution Model
More informationCOSMO-RS Theory. The Basics
Theory The Basics From µ to properties Property µ 1 µ 2 activity coefficient vapor pressure Infinite dilution Gas phase Pure compound Pure bulk compound Partition coefficient Phase 1 Phase 2 Liquid-liquid
More informationComparison of iterative inverse coarse-graining methods
EPJ manuscript No. (will be inserted by the editor) Comparison of iterative inverse coarse-graining methods David Rosenberger 1, Martin Hanke 2 and Nico F.A.van der Vegt 1,a 1 Eduard-Zintl-Institut für
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 informationHands-on : Model Potential Molecular Dynamics
Hands-on : Model Potential Molecular Dynamics OUTLINE 0. DL_POLY code introduction 0.a Input files 1. THF solvent molecule 1.a Geometry optimization 1.b NVE/NVT dynamics 2. Liquid THF 2.a Equilibration
More informationAccuracy, Transferability, and Efficiency of Coarse-Grained Models of Molecular Liquids
This is an open access article published under an ACS AuthorChoice License, which permits copying and redistribution of the article or any adaptations for non-commercial purposes. Cite This: pubs.acs.org/jpcb
More informationTitle Theory of solutions in the energy r of the molecular flexibility Author(s) Matubayasi, N; Nakahara, M Citation JOURNAL OF CHEMICAL PHYSICS (2003), 9702 Issue Date 2003-11-08 URL http://hdl.handle.net/2433/50354
More information2.3 Modeling Interatomic Interactions Pairwise Potentials Many-Body Potentials Studying Biomolecules: The Force
Contents 1 Introduction to Computational Meso-Bio-Nano (MBN) Science and MBN EXPLORER.... 1 1.1 Meso-Bio-Nano Science: A Novel Field of Interdisciplinary Research.... 1 1.1.1 Structure and Dynamics of
More informationMethods for searching and sampling configurational space
International Spring School Statistical Thermodynamics, Santiago de Chile Tuesday, ovember 28, 2017 Lecture 18 Methods for searching and sampling configurational space Prof. Dr. Wilfred F. van Gunsteren
More informationResearch Statement. Shenggao Zhou. November 3, 2014
Shenggao Zhou November 3, My research focuses on: () Scientific computing and numerical analysis (numerical PDEs, numerical optimization, computational fluid dynamics, and level-set method for interface
More informationThe change in free energy on transferring an ion from a medium of low dielectric constantε1 to one of high dielectric constant ε2:
The Born Energy of an Ion The free energy density of an electric field E arising from a charge is ½(ε 0 ε E 2 ) per unit volume Integrating the energy density of an ion over all of space = Born energy:
More informationAdvanced Topics in Equilibrium Statistical Mechanics
Advanced Topics in Equilibrium Statistical Mechanics Glenn Fredrickson 2. Classical Fluids A. Coarse-graining and the classical limit For concreteness, let s now focus on a fluid phase of a simple monatomic
More informationCoarse-grained Models for Oligomer-grafted Silica Nanoparticles
Modeling So+ Ma-er: Linking Mul3ple Length and Time Scales KITP Conference, Santa Barbara, June 4-8, 2012 Coarse-grained Models for Oligomer-grafted Silica Nanoparticles Bingbing Hong Alexandros Chremos
More informationSummary of the new Modelling Vocabulary
Summary of the new Modelling Vocabulary These two pages attempts to summarise in a concise manner the Modelling Vocabulary. What are Models? What are Simulations? Materials Models consist of Physics or
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 informationSupporting information. The role of halide ions in the anisotropic growth. of gold nanoparticles: a microscopic, atomistic.
Electronic Supplementary Material (ESI) for Physical Chemistry Chemical Physics. This journal is the Owner Societies 26 Supporting information The role of halide ions in the anisotropic growth of gold
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 informationTHE JOURNAL OF CHEMICAL PHYSICS 126,
THE JOURNAL OF CHEMICAL PHYSICS 16 44503 007 Development of an equation of state for electrolyte solutions by combining the statistical associating fluid theory and the mean spherical approximation for
More informationHyeyoung Shin a, Tod A. Pascal ab, William A. Goddard III abc*, and Hyungjun Kim a* Korea
The Scaled Effective Solvent Method for Predicting the Equilibrium Ensemble of Structures with Analysis of Thermodynamic Properties of Amorphous Polyethylene Glycol-Water Mixtures Hyeyoung Shin a, Tod
More informationNuggets on coarse- graining and mul,scale computa,onal schemes Maria Fyta Ins,tut für Computerphysik, Universität Stu<gart Stu<gart, Germany
Nuggets on coarse- graining and mul,scale computa,onal schemes Maria Fyta Ins,tut für Computerphysik, Universität Stu
More informationComputational Chemistry - MD Simulations
Computational Chemistry - MD Simulations P. Ojeda-May pedro.ojeda-may@umu.se Department of Chemistry/HPC2N, Umeå University, 901 87, Sweden. May 2, 2017 Table of contents 1 Basics on MD simulations Accelerated
More informationMonte Carlo Methods. Ensembles (Chapter 5) Biased Sampling (Chapter 14) Practical Aspects
Monte Carlo Methods Ensembles (Chapter 5) Biased Sampling (Chapter 14) Practical Aspects Lecture 1 2 Lecture 1&2 3 Lecture 1&3 4 Different Ensembles Ensemble ame Constant (Imposed) VT Canonical,V,T P PT
More informationSupporting Information for. Dynamics Study"
Supporting Information for "CO 2 Adsorption and Reactivity on Rutile TiO 2 (110) in Water: An Ab Initio Molecular Dynamics Study" Konstantin Klyukin and Vitaly Alexandrov,, Department of Chemical and Biomolecular
More informationThe effect of concentration- and temperature-dependent dielectric constant on the activity coefficient of NaCl electrolyte solutions.
The effect of concentration- and temperature-dependent dielectric constant on the activity coefficient of NaCl electrolyte solutions Mónika Valiskó, and Dezső Boda Department of Physical Chemistry, University
More informationBasic Ingredients of Free Energy Calculations: A Review
Feature Article Basic Ingredients of Free Energy Calculations: A Review CLARA D. CHRIST, 1 ALAN E. MARK, 2 WILFRED F. van GUNSTEREN 1 1 Laboratory of Physical Chemistry, Swiss Federal Institute of Technology,
More information8.333: Statistical Mechanics I Problem Set # 5 Due: 11/22/13 Interacting particles & Quantum ensembles
8.333: Statistical Mechanics I Problem Set # 5 Due: 11/22/13 Interacting particles & Quantum ensembles 1. Surfactant condensation: N surfactant molecules are added to the surface of water over an area
More informationMD simulation of methane in nanochannels
MD simulation of methane in nanochannels COCIM, Arica, Chile M. Horsch, M. Heitzig, and J. Vrabec University of Stuttgart November 6, 2008 Scope and structure Molecular model for graphite and the fluid-wall
More information3.320 Lecture 23 (5/3/05)
3.320 Lecture 23 (5/3/05) Faster, faster,faster Bigger, Bigger, Bigger Accelerated Molecular Dynamics Kinetic Monte Carlo Inhomogeneous Spatial Coarse Graining 5/3/05 3.320 Atomistic Modeling of Materials
More informationNon-bonded interactions
speeding up the number-crunching continued Marcus Elstner and Tomáš Kubař December 3, 2013 why care? number of individual pair-wise interactions bonded interactions proportional to N: O(N) non-bonded interactions
More informationHigh-order Chin actions in path integral Monte Carlo
High-order Chin actions in path integral Monte Carlo High-order actions and their applications, Barcelona 2009 Jordi Boronat Departament de Física i Enginyeria Nuclear Universitat Politècnica de Catalunya
More informationpolymers ISSN
Polymers 2014, 6, 1655-1675; doi:10.3390/polym6061655 Article OPEN ACCESS polymers ISSN 2073-4360 www.mdpi.com/journal/polymers A Coarse-Grained DNA Model Parameterized from Atomistic Simulations by Inverse
More informationStructuring of hydrophobic and hydrophilic polymers at interfaces Stephen Donaldson ChE 210D Final Project Abstract
Structuring of hydrophobic and hydrophilic polymers at interfaces Stephen Donaldson ChE 210D Final Project Abstract In this work, a simplified Lennard-Jones (LJ) sphere model is used to simulate the aggregation,
More informationNucleation rate (m -3 s -1 ) Radius of water nano droplet (Å) 1e+00 1e-64 1e-128 1e-192 1e-256
Supplementary Figures Nucleation rate (m -3 s -1 ) 1e+00 1e-64 1e-128 1e-192 1e-256 Calculated R in bulk water Calculated R in droplet Modified CNT 20 30 40 50 60 70 Radius of water nano droplet (Å) Supplementary
More informationCan a continuum solvent model reproduce the free energy landscape of a β-hairpin folding in water?
Can a continuum solvent model reproduce the free energy landscape of a β-hairpin folding in water? Ruhong Zhou 1 and Bruce J. Berne 2 1 IBM Thomas J. Watson Research Center; and 2 Department of Chemistry,
More informationLinking electronic and molecular structure: Insight into aqueous chloride solvation. Supplementary Information
Linking electronic and molecular structure: Insight into aqueous chloride solvation Ling Ge, Leonardo Bernasconi, and Patricia Hunt Department of Chemistry, Imperial College London, London SW7 2AZ, United
More informationSolvation and Macromolecular Structure. The structure and dynamics of biological macromolecules are strongly influenced by water:
Overview Solvation and Macromolecular Structure The structure and dynamics of biological macromolecules are strongly influenced by water: Electrostatic effects: charges are screened by water molecules
More informationSUPPLEMENTARY INFORMATION
DOI: 10.1038/NCHEM.2161 Coulomb explosion during the early stages of the reaction of alkali metals with water Philip E. Mason, 1 Frank Uhlig, 1 Václav Vaněk, 1 Tillmann Buttersack, 2 Sigurd Bauerecker,
More informationMONTE CARLO SIMULATION OF I -, Br -, AND Cl - IN WATER USING AB INITIO PAIR POTENSIAL FUNCTIONS
4 Indo. J. Chem., 7, 7 (), 4-9 MONTE CARLO SIMULATION OF I -, Br -, AND Cl - IN WATER USING AB INITIO PAIR POTENSIAL FUNCTIONS Austrian-Indonesian Centre for Computational Chemistry (AIC), Department of
More informationClassical and Quantum Mechanics Dr Mark R. Wormald. Bibliography
Classical and Quantum Mechanics 1 Classical and Quantum Mechanics Dr Mark R. Wormald Lecture 1. Introduction. Building blocks of matter, fundamental particles. Forces between particles. Electrostatics,
More informationMultiscale Materials Modeling
Multiscale Materials Modeling Lecture 06 Polymers These notes created by David Keffer, University of Tennessee, Knoxville, 2012. Outline Multiscale Modeling of Polymers I. Introduction II. Atomistic Simulation
More information2. Thermodynamics. Introduction. Understanding Molecular Simulation
2. Thermodynamics Introduction Molecular Simulations Molecular dynamics: solve equations of motion r 1 r 2 r n Monte Carlo: importance sampling r 1 r 2 r n How do we know our simulation is correct? Molecular
More informationJournal of Physics: Conference Series PAPER. To cite this article: Jiramate Kitjanon et al 2017 J. Phys.: Conf. Ser
Journal of Physics: Conference Series PAPER Transferability of Polymer Chain Properties between Coarse-Grained and Atomistic Models of Natural Rubber Molecule Validated by Molecular Dynamics Simulations
More informationContents. 1 Introduction and guide for this text 1. 2 Equilibrium and entropy 6. 3 Energy and how the microscopic world works 21
Preface Reference tables Table A Counting and combinatorics formulae Table B Useful integrals, expansions, and approximations Table C Extensive thermodynamic potentials Table D Intensive per-particle thermodynamic
More informationMachine learning the Born-Oppenheimer potential energy surface: from molecules to materials. Gábor Csányi Engineering Laboratory
Machine learning the Born-Oppenheimer potential energy surface: from molecules to materials Gábor Csányi Engineering Laboratory Interatomic potentials for molecular dynamics Transferability biomolecular
More informationMonte Carlo (MC) Simulation Methods. Elisa Fadda
Monte Carlo (MC) Simulation Methods Elisa Fadda 1011-CH328, Molecular Modelling & Drug Design 2011 Experimental Observables A system observable is a property of the system state. The system state i is
More informationRustam Z. Khaliullin University of Zürich
Rustam Z. Khaliullin University of Zürich Molecular dynamics (MD) MD is a computational method for simulating time evolution of a collection of interacting atoms by numerically integrating Newton s equation
More informationAb-initio simulation of liquid water by quantum Monte Carlo
Ab-initio simulation of liquid water by quantum Monte Carlo Sandro Sorella G. Mazzola & Y. Luo SISSA, IOM DEMOCRITOS, Trieste A. Zen, L. Guidoni U. of L Aquila, L Aquila 28 July 2014, Mike Towler Institute,
More informationPeptide folding in non-aqueous environments investigated with molecular dynamics simulations Soto Becerra, Patricia
University of Groningen Peptide folding in non-aqueous environments investigated with molecular dynamics simulations Soto Becerra, Patricia IMPORTANT NOTE: You are advised to consult the publisher's version
More informationGromacs Workshop Spring CSC
Gromacs Workshop Spring 2007 @ CSC Erik Lindahl Center for Biomembrane Research Stockholm University, Sweden David van der Spoel Dept. Cell & Molecular Biology Uppsala University, Sweden Berk Hess Max-Planck-Institut
More information