Dynamics of Polymers and Membranes P. B. Sunil Kumar
|
|
- Edward Lewis
- 6 years ago
- Views:
Transcription
1 Dynamics of Polymers and Membranes P. B. Sunil Kumar Department of Physics, Indian Institute of Technology Madras, Chennai,
2 Dissipative Particle Dynamics (DPD) - Basic model P. Espagnol and P. Warren, Europhys. Lett. 30, 191 (1995) N DPD" particles interacting via conservative, dissipative and random pair wise forces S'10"%"(+)53#)%* &1#3("1@"7#+:"CCC DPD (mesoscopic) Soft-repulsive +@@+)(#7+"conservative 313,-13.+."8F[F,?"-+%.$ Molecular Dynamics (MD) Special (DPD-)thermostat, #0&*+0+3(+."-2"0+%3$"1@! dissipative stochastic DPD beadsr '+&'+$+3(")*4$(+'$"1@" +*+0+3($ MD with special thermostat 8:1'A$"43.+'"+64#*#-'#40"B 313,+64#*#-'#40")13.#(#13$/ +C=C/"$5+%'#3=? DC"E1..+0%33/"9C"FG3:+=/"%3."HC"H'+0+'/ Phys. Rev. E 68 8IJJK?"JLMNJI" OP2.'1.23%0#) #3(+'%)(#13$Q 8@*1:"+@@+)($?"%'+"@4**2 (%A+3"#3(1"%))143(>
3 Conservative non-bonded pair forces C F i j r C ij Fi j :!(# ij &( rij ) eij, with &( rij ) :! (1" ) rc ' ' conservative force strength "weight function" r i e i j r i j r j F C ji!" F C ij r :! r " r ri j :! ri j e :! r / r i j i j i j i j i j F C ij 0$ # i j % e ij cut-off radius r c r ij 3
4 Conservative pair forces for for bonded interactions B$C$'2/'6%-'%(#$%2$6)$$%-*%,-+)3$I6)+/'/'6%/EC5/$2%/'%(#$%BGB%E-2$5J FAH%%&'%+22/(/-'J%G-($'(/+53%(#+(%$EC#+3/M$%(#$%entropy-dominated character -'$%,+'%,#--3$J!!#+)E-'/, C-($'(/+5!! OPQP%FO/'/($%PR($'3/N5$%Q-'5/'$+)%P5+3(/,H C-($'(/+5!! $'()-C/,%3C)/'6 C-($'(/+5 S%%N$'2/'6%C-($'(/+539%(-)3/-'%C-($'(/+53 4
5 Thermostat forces: Dissipative pair forces vi D F i j r i e i j r ij r j 2 2 D F :#! $! %( r ) ( e " v ) e #! $! %( r ) ( e & e v j r :# r! r ri j :# ri j e :# r / r i j i j i j i j i j!v j!"#$%&'($)*+,$%-*%./*$0%1 &'2/+'%&'3(/(4($%-*%"$,#'-5-67%8+2)+39%:#$''+/9%&'2/+9%;<1=>%?+'@%A;;> =; ) v ij ij ij ij ij ij ij ij ij " dissipative force strengh t F D ji #! F D ij %! ( r ij )% is another weight function (see below) v :# v! v i j i j v e i j i j vi i j 5 ij v j ( e "v ) e i j
6 Thermostat forces: Random pair forces (continuum limit) R!" ( r Fi j : $ " "!( r ij )# i j e ij )% is the same weight ij function as in the case! random force strength of the dissipative force (see below) # % % % ij are random variables with the following properties : symmetric for each interacting particle pair at all times : # ( t) $ # ( t) ij zero mean for all particle pairs at all times : ji &# ( t) '$ 0 ij unit variance for all particle pairs at all times, and uncorrelated for different pairs of particles and different times, : &# ( t) # ( t( )' $ () ) * ) ) ) )+ t, t( ) ij kl ik jl il jk 2 2 This can be fixed by the substitution!! /, % I% Gaussian%+'2%+% uniform %% C)-J+J/5/(7%2$'3/( 7%*-)%(#$ %%%%)+'2-E%D+)/+J5$3%K-)L% %% $M4+557%K$ 55@ % "#$%2$*/'/(/-'%*-)%(#$% %%%%)+'2-E%*-),$%+3%3#-K'% %%%%#$)$%/3%*-)%(#$% continuum limit FNO/$' $) C)-, $ 33NH@ %%%% P-)% discrete (/E$%3($C3% %%%% 0-.t # 1 (#/3%#+3%(-%J$% %%%%E-2/*/$2%F3$$%'$Q(%C+ 6$HR ( # t i.e., by! R R!! (, and F "! ( F " # t # t! 6
7 Thermostat forces: Random pair forces (discrete case) t 2 # $ 2 R (2) x( t) ) 5 F ( t% ) dt% * + 0, & ' ( & '!!!! Discretization: Divide [0, t]! into! n! steps of size! 0t " t/ n. 2 n n n 2 2 # R - t. $ R R - t. ( ) ) 6 Fi 1 2 Fi Fj. i" 1 n * ' " 66 & ' i" 1 j" 1 3 n4 / & x t ' ( & +, R R 2 Assume the discrete analogon t o & F ( t) F ( t% )' "! 7 ( t8 t% ) is & F F ' "! 7 R R 2 i j ij n t. 2 t 2 - t. 2 ( ). 6 i" 1. / & x t ' (! 1 2 "! "! t 1 2 "! t 0t 3 n 4 n 3 n 4 But x t of t; for an ordinary random walk & x t ' ( t Assu ' " 7 2 & ( ) ' must be independent 0 2 R R 2 ( ) must hold / mption & Fi Fj! ij was wrong. 7
8 Why is the DPD thermostat able to establish a chosen (set) temperature kbt* in the system? Dissipative forces always take energy out of the system. This mimics energy transfer to atomistic degrees of freedom that are not explicitly taken into account in the mesoscopic DPD model. Random kick forces on average transfer energy(from hidden atomistic degrees of freedom) into the system. A given temperature kbt* can only be established if the to and from energy transfers on average balance each other over time. 8
9 Lett. 30 G=QQRI%=Q=%% # S7%F$+'3%-*%(#$%CHC%(#$)F-3(+(%+%373($F%-*%D+)(/,5$3%G!CHCTK$+230I /3%2$3,)/K$2%/'%(#$%canonical ensemble9%/@$9%at constant temperature kt B * /**! () r G=I "#$%same!j$/6#(%*4',(/-'0 +DD$+)3% linear /'%(#$%)+'2-F%*-),$%F R 9%+'2%quadratic /'%(#$%2/33/D+(/E$%*-),$%F () r &'%D)/',/D5$%(#$%*-)F%-*%%%%%%%%%/3%+)K/()+)7%G+3%5-'6%+3%/(3%*/)3(%(J-%F-F$'(3 $L/3(I9%+'2%/'%(#$%standard formulation -'$%,#--3$3%(#$%same weight function as in the conservative forcem!! ( r) "! ( r). GAI "#$%2/33/D+(/E$%*-),$%3()$'6(#%$ +'2%(#$%)+'2-F%*-),$%3()$'6(#%% 3+(/3*7%(#$ fluctuation-dissipation relation kt* & B 2 % 2$ 9
10 Discretized equations of motion in DPD: (Modified) velocity-verlet integration scheme!"#$%&'()*"#+$#(,!"#$%&'()*+!,+-,./ &0+12+,&)-"!'&$0!-.-/"#$%&'() "#+$#(,!"#$%&'() 345%.!/+6+5-''.%78 9:;+<.,$"/*+=;+>!''-"!&0.0*+ Phys. Rev. E 01 3BCCC8+DEFGGH 10
11 How is a Lipid Simulated? Model Hydrophilic Head Group Choline Phosphate Glycerol Hydrophobic Tail Group: Fatty Acids e.g. Phosphatidylcholine
12 Lipids - Interaction parameters for conservative forces In addition the particles in a lipid are connected by springs We simulate 16,000 lipids in a box of 86x86x40, with a total of approximately one million particles.
13 13
14 DPD simulations of membranes in constant tension ensemble The method is similar to the Nos - Hoover method for NVT simulations. é We will introduce a variable ɛ = ln(v/v 0 ) p ɛ and associated mass W. The DPD integration scheme is now modified as, v ɛ v ɛ v ɛ v ɛ ɛ ɛ v i v i p ɛ p ɛ F ɛ t Update ɛ ɛ + p ɛ W t V,F C i, F D i, F R i and P = 1 dv [ i r i exp(ɛ ɛ ) {r i + v i t} p 2 i m i + i F C i r i ] Iterate for v i and p ɛ such that P is equal to P 0 ṽ i v i p ɛ p ɛ Initial guess : v i exp(ɛ ɛ )v i p ɛ p ɛ F ɛ t A. F. Jakobsen, JCP 122, ( 2005) ; 125, ( 2006) {( F C i m i + FD i m i 2 p ɛ W v i {( F C i m i 1+p ɛ t/w ) + FD i m i ) t + FR i m i t } Compute F D i and construct a piston with momentum } t + FR i t m i p ɛ p ɛ F ɛ t F e = dv [P P 0 ]+ d dn d 14 i p 2 i m i
15 Tension in a membrane is given by γ s =<L z [P z 1 2 (P x + P y )] > Lz is the dimension of the simulation box normal to the membrane. To obtain the desired values of the surface tension ϒs and the normal pressure PN, we will replace the parameter pε by a tensor. d The first term in the equation F e = dv [P P 0 ]+ dn d should then be, for α =(x, y), and for the normal direction F α = da [ γ 0 s L z (P 0 N P α ) ] + F α = dv [ P z P 0 N ] + d dn d d dn d i p 2 i m i i i p 2 i m i p 2 i m i It is also possible to control the tension by a Langevin piston method in which the force on the piston is given by d p 2 i F α = dv [P α P 0 ]+ γ p p αα + σ p ξ α dn d m i i with σ 2 p =2γ p Wk b T 15 A. F. Jakobsen, JCP 122, ( 2005) ; 125, ( 2006)
16 Hydrophilic Hydrophobic (a) Spherical Micelle, (b) Cylindrical Micelle, (c) Reverse Micelle, (d) Bilayer, (e) Hexagonal phase, (f) Vesicles Cylindrical micelles also called living polymers Can undergo scission/recombination under thermal fluctuations 2
17 DPD model for living polymers ( LP) For, W W, W I W A, I I, I A For, A-A r 1 r 3 r 5 9
18 A trimer based model for living polymers The three-body bending potential I 1 A 1 A 2 & A 1 A 2 I 2 Linear spring connecting A 1 -I 1, A 1 -I 1 and A 1 -A 1 For LP For solvent Simulation box size L = 40, Periodic Boundary Condition density ρ = 3, number of particles Resting length A-I = 1.0, A-A = 2.0 Spring constant k = 200 Time step = 0.01 τ, with τ of the order of 10-6 sec We use LAMMPS ( modified) 10
19 Fluid-gel transition C -- % ratio of number of monomers to the total number of particles Phase behavior was characterized by diffusivity of trimers Fluid phase at low C Gel phase at higher C Transition point C = Cp = 3.5 Branched point coordination number of particle is > 2 11
20 Organization of polymers 12
21 Organization of polymers C = 2 12
22 Organization of polymers C = 2 C = 3 12
23 Organization of polymers C = 2 Fluid C = 3 12
24 Organization of polymers C = 2 Fluid C = 3 Increasing C 12
25 Organization of polymers C = 2 Fluid C = 3 Increasing C Gel 12
26 Organization of polymers C = 2 Fluid C = 3 Increasing C C = 4 Gel 12
27 Organization of polymers C = 2 Fluid C = 3 Increasing C C = 4 Gel C = 5 12
28 Segment length distribution Cluster size distribution branching of polymers lead to random percolation l Φ Average polymer length Monomer concentration Branching reduces the average segment length M. E. Cates and S. J. Candau, J. Phys 2,6869 (1990) Exponential distribution at low C < Cp Power law distribution at C = Cp: random percolation Gel phase contains spanning clusters and small segments 13
29 Self-intermediate scattering function N: total number of monomers Relaxation time Zimm dynamics at large q: Diffusive at low q: 14
30 First recombination time τ R = scission time (t 2 ) recombination time (t 1 ) Two classes of recombination kinetic: Mean field (MF) and Diffusion controlled (DC) Ben O Shaughnessy and Jane Yu, PRL, 74, 4329 (1995) t d t d t s t d 15
31 First recombination time τ R = scission time (t 2 ) recombination time (t 1 ) Two classes of recombination kinetic: Mean field (MF) and Diffusion controlled (DC) Ben O Shaughnessy and Jane Yu, PRL, 74, 4329 (1995) t d t d t s t d 15
32 stress relaxation Residual stress is due to the spanning clusters with stretched bonds Early time oscillation is the result of bond potential within a trimer G(t) decays as: G. Faivre and J. L. Gardissat, Macromolecule, 19, 1988 (1986) 16
33 Shear V* z y -V* x Density of wall = 20* Density of fluids Periodic boundary in x & y directions Results were verified using Lees-Edward method to rule out wall effect 17
34 Shear induced structures x z x y Lamellar phase in the system with C = 5 y 18
35 Shear induced Lamallar phase After 2000 DPD time step Initiated by migration of polymers towards the higher velocity regions First layer forms near the shear boundaries After DPD time step All the polymers within a range of r 5 is pulled towards the boundary Next layer forms exactly at a distance r 5 from the first 19
36 Shear induced Lamallar phase After 2000 DPD time step Initiated by migration of polymers towards the higher velocity regions First layer forms near the shear boundaries After DPD time step All the polymers within a range of r 5 is pulled towards the boundary Next layer forms exactly at a distance r 5 from the first 19
37 Shear induced Lamallar phase After 2000 DPD time step Initiated by migration of polymers towards the higher velocity regions First layer forms near the shear boundaries After DPD time step All the polymers within a range of r 5 is pulled towards the boundary Next layer forms exactly at a distance r 5 from the first 19
38 Shear induced Lamallar phase After 2000 DPD time step Initiated by migration of polymers towards the higher velocity regions First layer forms near the shear boundaries After DPD time step All the polymers within a range of r 5 is pulled towards the boundary Next layer forms exactly at a distance r 5 from the first 19
39 Shear induced structures x z x y Columnar phase in the system with C = 5 Layer spacing is decided by r 5 y 20
40 Stress and viscosity Alignment of polymers leads to decrease in viscosity 21
41 Polyelctrolytes Macromolecules which dissociates into charged polymer and counterions (CI) in polar solvents are known as polyelectrolytes (PE). Anchored polyelectrolytes are of importance for many biologically and industrially relevant processes We will look at some of the equilibrium and dynamical properties of a single grafted polyelectrolyte using Brownian dynamics and Monte Carlo simulations.
42 Polyelctrolytes Macromolecules which dissociates into charged polymer and counterions (CI) in polar solvents are known as polyelectrolytes (PE). Anchored polyelectrolytes are of importance for many biologically and industrially relevant processes We will look at some of the equilibrium and dynamical properties of a single grafted polyelectrolyte using Brownian dynamics and Monte Carlo simulations.
43 The Model Total energy of the system consists of three parts, (i) Lennard-Jones potential Diameter of the bead a Length of the box L Number of monomers N Impenetrable, non-polarizable walls at x=0 and x=l Periodic minimum image boundary conditions along the y and z directions. (ii) harmonic spring potential between connected beads (iii) Electrostatic potential between any pair of beads,
44 The Model Total energy of the system consists of three parts, (i) Lennard-Jones potential Diameter of the bead a Length of the box L Number of monomers N Impenetrable, non-polarizable walls at x=0 and x=l Periodic minimum image boundary conditions along the y and z directions. (ii) harmonic spring potential between connected beads (iii) Electrostatic potential between any pair of beads,
45 Polyelectrolyte Phase Diagram K. Jayasree, P. Ranjith, Madan Rao and PBSK J. Chem. Phys. 130, (2009) In our simulation, we model the nonelectrostatic interactions by a Lennard-Jones potential monomer-monomer interaction ε= ε 1 = ε 2 CI-CI and monomer-ci interactions ε 1 =1, ε 2 =0 ε=3,a=3.0 ε=3,a=0.1
46 The free energy from umbrella sampling Monte Carlo
47 The free energy from umbrella sampling Monte Carlo
48 Fluctuations The asphericity Where, are the eigenvalues of the gyration tensor
49 Fluctuations The asphericity Where, are the eigenvalues of the gyration tensor
50 Brownian Dynamics Effect of counter-ion condensation on polyelectrolytes
51 Brownian Dynamics Effect of counter-ion condensation on polyelectrolytes
52 Brownian Dynamics Effect of counter-ion condensation on polyelectrolytes
53 Brownian dynamics with external shear The equation of motion for the monomers is dr i dt = j M ij rj U(t)+v(r i )+ 2k B TB ij f j (t) v(r) = γyˆx M ij = 1 6πηa I [ 1 3 6πηr ij is the solvent velocity profile M = B B T 4 (I + r ijr ij r 2 ij ] )+ a2 (I 3 r ijr ij ) 2rij 2 rij 2 is the Rotne-Prager tensor i=j i j, r ij 2a 1 6πηa [ (1 9r ij 32a )I +3r ijr ij 32ar ij ] i j, r ij 2a The components of B is obtained recursively by the Cholesky decomposition of the 36 mobility tensor M
54 Lees Edward periodic boundary condition Our simple periodic boundary condition should now be modified to the Lees-Edward periodic boundary condition. If the flow is in the X-direction and the velocity gradient in the Y-direction then, a particle leaving the top boundary at y=l with velocity (vx,vy,vz) will appear at the bottom at y=0 with new x-coordinate and new x-component of velocity given by, x new = x n t U t int[n t U t/l]l 37
Supplementary 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 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 informationDissipative Particle Dynamics: Foundation, Evolution and Applications
Dissipative Particle Dynamics: Foundation, Evolution and Applications Lecture 4: DPD in soft matter and polymeric applications George Em Karniadakis Division of Applied Mathematics, Brown University &
More informationKobe-Brown Simulation Summer School 2015 Project: DPD Simulation of a Membrane
Kobe-Brown Simulation Summer School 2015 Project: DPD Simulation of a Membrane Clark Bowman Karen Larson Yuji Funaki Ross Parker Tae Woo Kim Week 1: Introduction to Molecular Dynamics (MD) Computer simulation
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 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 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 informationCEE 618 Scientific Parallel Computing (Lecture 12)
1 / 26 CEE 618 Scientific Parallel Computing (Lecture 12) Dissipative Hydrodynamics (DHD) Albert S. Kim Department of Civil and Environmental Engineering University of Hawai i at Manoa 2540 Dole Street,
More informationFORCE ENERGY. only if F = F(r). 1 Nano-scale (10 9 m) 2 Nano to Micro-scale (10 6 m) 3 Nano to Meso-scale (10 3 m)
What is the force? CEE 618 Scientific Parallel Computing (Lecture 12) Dissipative Hydrodynamics (DHD) Albert S. Kim Department of Civil and Environmental Engineering University of Hawai i at Manoa 2540
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 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 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 informationModeling and simulation of nanoparticle aggregation in colloidal systems
Retrospective Theses and Dissertations Iowa State University Capstones, Theses and Dissertations 2008 Modeling and simulation of nanoparticle aggregation in colloidal systems Sergiy Markutsya Iowa State
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 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 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 informationLecture 5: Macromolecules, polymers and DNA
1, polymers and DNA Introduction In this lecture, we focus on a subfield of soft matter: macromolecules and more particularly on polymers. As for the previous chapter about surfactants and electro kinetics,
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 informationDissipative particle dynamics: Bridging the gap between atomistic and mesoscopic simulation
Dissipative particle dynamics: Bridging the gap between atomistic and mesoscopic simulation Robert D. Groot and Patrick B. Warren Unilever Research Port Sunlight Laboratory, Quarry Road East, Bebington,
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 informationSimulation of Coarse-Grained Equilibrium Polymers
Simulation of Coarse-Grained Equilibrium Polymers J. P. Wittmer, Institut Charles Sadron, CNRS, Strasbourg, France Collaboration with: M.E. Cates (Edinburgh), P. van der Schoot (Eindhoven) A. Milchev (Sofia),
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 informationMicroscopic Picture of Aging in SiO 2 : A Computer Simulation
Microscopic Picture of Aging in SiO 2 : A Computer Simulation Katharina Vollmayr-Lee, Robin Bjorkquist, Landon M. Chambers Bucknell University & Göttingen 7 td tb r n (t) 6 R 5 4 3 2 ti t [ns] waiting
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 informationOrganization of NAMD Tutorial Files
Organization of NAMD Tutorial Files .1.1. RMSD for individual residues Objective: Find the average RMSD over time of each residue in the protein using VMD. Display the protein with the residues colored
More informationModeling the combined effect of surface roughness and shear rate on slip flow of simple fluids
Modeling the combined effect of surface roughness and shear rate on slip flow of simple fluids Anoosheh Niavarani and Nikolai Priezjev www.egr.msu.edu/~niavaran November 2009 A. Niavarani and N.V. Priezjev,
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 informationDissipative Particle Dynamics
CHAPTER 2 Dissipative Particle Dynamics Igor V. Pivkin, a,b Bruce Caswell, a and George Em Karniadakis a a Division of Applied Mathematics, Brown University, Providence, Rhode Island b Department of Materials
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 informationOn the Dynamics and Disentanglement in Thin and Two-Dimensional Polymer Films
J. Phys. IV France 1 (006) Pr1-1 c EDP Sciences, Les Ulis On the Dynamics and Disentanglement in Thin and Two-Dimensional Polymer Films H. Meyer, T. Kreer, A. Cavallo, J. P. Wittmer and J. Baschnagel 1
More informationφ(z) Application of SCF to Surfaces and Interfaces (abridged from notes by D.J. Irvine)
Application of SCF to Surfaces and Interfaces (abridged from notes by D.J. Irvine) Edwards continuum field theory reviewed above is just one flavor of selfconsistent mean field theory, but all mean field
More informationGeneral-purpose Coarse-grained. Molecular Dynamics Program COGNAC
General-purpose Coarse-grained Molecular Dynamics Program COGNAC (COarse-Grained molecular dynamics program by NAgoya Cooperation) Takeshi Aoyagi JCII, Doi project 0 sec -3 msec -6 sec -9 nsec -12 psec
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 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 informationFluctuating Hydrodynamics Approaches for Mesoscopic Modeling and Simulation Applications in Soft Materials and Fluidics
Fluctuating Hydrodynamics Approaches for Mesoscopic Modeling and Simulation Applications in Soft Materials and Fluidics Theoretical Background and Applications Summer School on Multiscale Modeling of Materials
More informationSupporting Information for: Physics Behind the Water Transport through. Nanoporous Graphene and Boron Nitride
Supporting Information for: Physics Behind the Water Transport through Nanoporous Graphene and Boron Nitride Ludovic Garnier, Anthony Szymczyk, Patrice Malfreyt, and Aziz Ghoufi, Institut de Physique de
More informationThe dynamics of small particles whose size is roughly 1 µmt or. smaller, in a fluid at room temperature, is extremely erratic, and is
1 I. BROWNIAN MOTION The dynamics of small particles whose size is roughly 1 µmt or smaller, in a fluid at room temperature, is extremely erratic, and is called Brownian motion. The velocity of such particles
More informationMesoscale simulations of stimuli-sensitive polymer networks
Mesoscale simulations of stimuli-sensitive polymer networks Alexander Alexeev Complex Fluids Modeling and Simulation Group The George W. Woodruff School of Mechanical Engineering Georgia Institute of Technology
More informationRelaxation time of non-linear polymers in dilute solution via computer simulation
Journal of Non-Crystalline Solids 352 (2006) 5081 5086 www.elsevier.com/locate/jnoncrysol Relaxation time of non-linear polymers in dilute solution via computer simulation J.G. Hernández Cifre *, R. Pamies,
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 information1 Bulk Simulations in a HCP lattice
1 Bulk Simulations in a HCP lattice 1.1 Introduction This project is a continuation of two previous projects that studied the mechanism of melting in a FCC and BCC lattice. The current project studies
More informationStructure of a tethered polymer under flow using molecular dynamics and hybrid molecular-continuum simulations
Physica A 362 (2006) 30 35 www.elsevier.com/locate/physa Structure of a tethered polymer under flow using molecular dynamics and hybrid molecular-continuum simulations Rafael Delgado-Buscalioni a, Peter
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 informationCoarse-Grained Models!
Coarse-Grained Models! Large and complex molecules (e.g. long polymers) can not be simulated on the all-atom level! Requires coarse-graining of the model! Coarse-grained models are usually also particles
More informationGlass-Transition and Side-Chain Dynamics in Thin Films: Explaining. Dissimilar Free Surface Effects for Polystyrene and Poly(methyl methacrylate)
Supporting Information for Glass-Transition and Side-Chain Dynamics in Thin Films: Explaining Dissimilar Free Surface Effects for Polystyrene and Poly(methyl methacrylate) David D. Hsu, Wenjie Xia, Jake
More informationModelling across different time and length scales
18/1/007 Course MP5 Lecture 1 18/1/007 Modelling across different time and length scales An introduction to multiscale modelling and mesoscale methods Dr James Elliott 0.1 Introduction 6 lectures by JAE,
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 informationSupporting Information
Projection of atomistic simulation data for the dynamics of entangled polymers onto the tube theory: Calculation of the segment survival probability function and comparison with modern tube models Pavlos
More informationAutocorrelation Study for a Coarse-Grained Polymer Model
Autocorrelation Study for a Coarse-Grained Polymer Model KaiK.Qi and andmichael M. Bachmann Bachmann he Center for CSP,UGA Simulational Physics he University of Georgia alk given at the 14th Leipzig Workshop
More informationTime-Dependent Statistical Mechanics 1. Introduction
Time-Dependent Statistical Mechanics 1. Introduction c Hans C. Andersen Announcements September 24, 2009 Lecture 1 9/22/09 1 Topics of concern in the course We shall be concerned with the time dependent
More informationEXAM I COURSE TFY4310 MOLECULAR BIOPHYSICS December Suggested resolution
page 1 of 7 EXAM I COURSE TFY4310 MOLECULAR BIOPHYSICS December 2013 Suggested resolution Exercise 1. [total: 25 p] a) [t: 5 p] Describe the bonding [1.5 p] and the molecular orbitals [1.5 p] of the ethylene
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 informationFor one such collision, a new particle, "k", is formed of volume v k = v i + v j. The rate of formation of "k" particles is, N ij
Coagulation Theory and The Smoluchowski Equation: Coagulation is defined as growth of particles by collisions among particles. It is usually associated with dense, 3-d growth, in contrast to aggregation
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 informationDistribution of chains in polymer brushes produced by a grafting from mechanism
SUPPLEMENTARY INFORMATION Distribution of chains in polymer brushes produced by a grafting from mechanism Andre Martinez, Jan-Michael Y. Carrillo, Andrey V. Dobrynin,, * and Douglas H. Adamson, * Department
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 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 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 informationExplaining and modelling the rheology of polymeric fluids with the kinetic theory
Explaining and modelling the rheology of polymeric fluids with the kinetic theory Dmitry Shogin University of Stavanger The National IOR Centre of Norway IOR Norway 2016 Workshop April 25, 2016 Overview
More informationHydrodynamics in Espresso. Kai Grass 1
Hydrodynamics in Espresso Kai Grass 1 Hydrodynamics in Espresso Explicit solvent: Lennard-Jones liquid at constant energy Dissipative Particle Dynamics (DPD): thermostated particle fluid Lattice Boltzmann
More informationFast Brownian Dynamics for Colloidal Suspensions
Fast Brownian Dynamics for Colloidal Suspensions Aleksandar Donev, CIMS and collaborators: Florencio Balboa CIMS) Andrew Fiore and James Swan MIT) Courant Institute, New York University Modeling Complex
More informationPre-yield non-affine fluctuations and a hidden critical point in strained crystals
Supplementary Information for: Pre-yield non-affine fluctuations and a hidden critical point in strained crystals Tamoghna Das, a,b Saswati Ganguly, b Surajit Sengupta c and Madan Rao d a Collective Interactions
More informationPolymers Dynamics by Dielectric Spectroscopy
Polymers Dynamics by Dielectric Spectroscopy What s a polymer bulk? A condensed matter system where the structural units are macromolecules Polymers Shape of a Macromolecule in the Bulk Flory's prediction
More informationAnalytical and Numerical Study of Coupled Atomistic-Continuum Methods for Fluids
Analytical and Numerical Study of Coupled Atomistic-Continuum Methods for Fluids Weiqing Ren Courant Institute, NYU Joint work with Weinan E Multiscale modeling for two types of problems: Complex fluids
More informationCOMPUTER SIMULATION OF TRACER DIFFUSION IN GEL NETWORK
COMPUTER SIMULATION OF TRACER DIFFUSION IN GEL NETWORK ZHOU HUAI A THESIS SUBMITED FOR THE DEGREE OF DOCTOR OF PHILOSOPHY DEPARTMENT OF CHEMICAL AND BIOMOLECULAR ENGINEERING NATIONAL UNIVERSITY OF SINGAPORE
More informationDynamic modelling of morphology development in multiphase latex particle Elena Akhmatskaya (BCAM) and Jose M. Asua (POLYMAT) June
Dynamic modelling of morphology development in multiphase latex particle Elena Akhmatskaya (BCAM) and Jose M. Asua (POLYMAT) June 7 2012 Publications (background and output) J.M. Asua, E. Akhmatskaya Dynamic
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 informationSupplementary Information for Atomistic Simulation of Spinodal Phase Separation Preceding Polymer Crystallization
Supplementary Information for Atomistic Simulation of Spinodal Phase Separation Preceding Polymer Crystallization Richard H. Gee * Naida Lacevic and Laurence E. Fried University of California Lawrence
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 information6 Similarities between principal components of protein dynamics and random diffusion
6 Similarities between principal components of protein dynamics and random diffusion B. Hess Phys. Rev. E 6:8438-8448 () Principal component analysis, also called essential dynamics, is a powerful tool
More informationThermodynamic behaviour of mixtures containing CO 2. A molecular simulation study
Thermodynamic behaviour of mixtures containing. A molecular simulation study V. Lachet, C. Nieto-Draghi, B. Creton (IFPEN) Å. Ervik, G. Skaugen, Ø. Wilhelmsen, M. Hammer (SINTEF) Introduction quality issues
More informationminimal models for lipid membranes: local modifications around fusion objects
minimal models for lipid membranes: local modifications around fusion objects Giovanni Marelli Georg August Universität, Göttingen January 21, 2013 PhD defense collective phenomena structures and morphologies
More informationChapter 1 Introduction
Chapter 1 Introduction This thesis is concerned with the behaviour of polymers in flow. Both polymers in solutions and polymer melts will be discussed. The field of research that studies the flow behaviour
More informationINTERMOLECULAR AND SURFACE FORCES
INTERMOLECULAR AND SURFACE FORCES SECOND EDITION JACOB N. ISRAELACHVILI Department of Chemical & Nuclear Engineering and Materials Department University of California, Santa Barbara California, USA ACADEMIC
More informationColloidal Suspension Rheology Chapter 1 Study Questions
Colloidal Suspension Rheology Chapter 1 Study Questions 1. What forces act on a single colloidal particle suspended in a flowing fluid? Discuss the dependence of these forces on particle radius. 2. What
More informationPolymer Dynamics and Rheology
Polymer Dynamics and Rheology 1 Polymer Dynamics and Rheology Brownian motion Harmonic Oscillator Damped harmonic oscillator Elastic dumbbell model Boltzmann superposition principle Rubber elasticity and
More informationChapter 6: The Rouse Model. The Bead (friction factor) and Spring (Gaussian entropy) Molecular Model:
G. R. Strobl, Chapter 6 "The Physics of Polymers, 2'nd Ed." Springer, NY, (1997). R. B. Bird, R. C. Armstrong, O. Hassager, "Dynamics of Polymeric Liquids", Vol. 2, John Wiley and Sons (1977). M. Doi,
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 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 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 informationProteins polymer molecules, folded in complex structures. Konstantin Popov Department of Biochemistry and Biophysics
Proteins polymer molecules, folded in complex structures Konstantin Popov Department of Biochemistry and Biophysics Outline General aspects of polymer theory Size and persistent length of ideal linear
More informationGEM4 Summer School OpenCourseWare
GEM4 Summer School OpenCourseWare http://gem4.educommons.net/ http://www.gem4.org/ Lecture: Polymer Chains by Ju Li. Given August 16, 2006 during the GEM4 session at MIT in Cambridge, MA. Please use the
More informationDynamics of the constrained polymer collapse
EUROPHYSICS LETTERS 15 June 2003 Europhys. Lett., 62 (6), pp. 848 854 (2003) Dynamics of the constrained polymer collapse I. S. Aranson 1 and L. S. Tsimring 2 1 Materials Science Division, Argonne National
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 informationAn efficient hybrid method for modeling lipid membranes with molecular resolution
An efficient hybrid method for modeling lipid membranes with molecular resolution G.J.A. (Agur) Sevink, M. Charlaganov & J.G.E.M Fraaije Soft Matter Chemistry Group, Leiden University, The Netherlands
More informationNew from Ulf in Berzerkeley
New from Ulf in Berzerkeley from crystallization to statistics of density fluctuations Ulf Rørbæk Pedersen Department of Chemistry, University of California, Berkeley, USA Roskilde, December 16th, 21 Ulf
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 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 informationA comparison of lattice-boltzmann and Brownian-dynamics simulations of polymer migration in confined flows. Abstract
A comparison of lattice-boltzmann and Brownian-dynamics simulations of polymer migration in confined flows Rahul Kekre, Jason E. Butler, and Anthony J. C. Ladd Department of Chemical Engineering, University
More informationTutorial School on Fluid Dynamics: Aspects of Turbulence Session I: Refresher Material Instructor: James Wallace
Tutorial School on Fluid Dynamics: Aspects of Turbulence Session I: Refresher Material Instructor: James Wallace Adapted from Publisher: John S. Wiley & Sons 2002 Center for Scientific Computation and
More informationMigration Behavior of Bead-spring Dumbbell Models under Microchannel Flow from Dissipative Particle Dynamics Simulations
2426 Bull. Korean Chem. Soc. 2007, Vol. 28, No. 12 Kwang Jin Oh Migration Behavior of Bead-spring Dumbbell Models under Microchannel Flow from Dissipative Particle Dynamics Simulations Kwang Jin Oh Supercomputing
More informationStatistical and hydrodynamic properties of topological polymers for various graphs showing enhanced short-range correlation
arxiv:1607.08726v1 [cond-mat.soft] 29 Jul 2016 Statistical and hydrodynamic properties of topological polymers for various graphs showing enhanced short-range correlation Erica Uehara and Tetsuo Deguchi
More information1.3 Molecular Level Presentation
1.3.1 Introduction A molecule is the smallest chemical unit of a substance that is capable of stable, independent existence. Not all substances are composed of molecules. Some substances are composed of
More informationLines of Renormalization Group Fixed Points for Fluid and Crystalline Membranes.
EUROPHYSICS LETTERS 1 October 1988 Europhys. Lett., 7 (3), pp. 255-261 (1988) Lines of Renormalization Group Fixed Points for Fluid and Crystalline Membranes. R. LIPOWSKY Institut für Festkörperforschung
More informationICCP Project 2 - Advanced Monte Carlo Methods Choose one of the three options below
ICCP Project 2 - Advanced Monte Carlo Methods Choose one of the three options below Introduction In statistical physics Monte Carlo methods are considered to have started in the Manhattan project (1940
More informationSurveying Free Energy Landscapes: Applications to Continuum Soft Matter Systems
Surveying Free Energy Landscapes: Applications to Continuum Soft Matter Systems Halim Kusumaatmaja Department of Physics, University of Durham Acknowledgement David Wales (Cambridge) Discussions on various
More informationEffect of surfactant structure on interfacial properties
EUROPHYSICS LETTERS 15 September 2003 Europhys. Lett., 63 (6), pp. 902 907 (2003) Effect of surfactant structure on interfacial properties L. Rekvig 1 ( ), M. Kranenburg 2, B. Hafskjold 1 and B. Smit 2
More information2. Molecules in Motion
2. Molecules in Motion Kinetic Theory of Gases (microscopic viewpoint) assumptions (1) particles of mass m and diameter d; ceaseless random motion (2) dilute gas: d λ, λ = mean free path = average distance
More informationPart III. Polymer Dynamics molecular models
Part III. Polymer Dynamics molecular models I. Unentangled polymer dynamics I.1 Diffusion of a small colloidal particle I.2 Diffusion of an unentangled polymer chain II. Entangled polymer dynamics II.1.
More information