Adaptive Resolution Simulation of Molecular Liquids
|
|
- Rosanna Dennis
- 5 years ago
- Views:
Transcription
1 Adaptive Resolution Simulation of Molecular Liquids Matej Praprotnik Max Planck Institute for Polymer Research, Theory Group, Mainz, Germany On leave from the National Institute of Chemistry, Ljubljana, Slovenia Universität Stuttgart - SFB 76, Pforzheim-Hohenwart, September 8-, 28 p. /59
2 Collaborators Kurt Kremer, Max Planck Institute for Polymer Research, Mainz, Germany Luigi Delle Site, Max Planck Institute for Polymer Research, Mainz, Germany Christoph Junghans, Max Planck Institute for Polymer Research, Mainz, Germany Rafael Delgado Buscalioni, Universidad Autonoma de Madrid, Madrid, Spain Cecilia Clementi, Rice University, Houston, Texas Silvina Matysiak, Rice University, Houston, Texas Universität Stuttgart - SFB 76, Pforzheim-Hohenwart, September 8-, 28 p. 2/59
3 Outline Introduction: molecular liquids and soft matter, length- and time-scales, multiscale modeling Hybrid Adaptive Resolution Scheme (AdResS): linking atomic and mesoscopic length-scales Triple-scale AdResS-HybridMD method: coupling particle-based and continuum descriptions of a liquid Universität Stuttgart - SFB 76, Pforzheim-Hohenwart, September 8-, 28 p. 3/59
4 Molecular systems gases: intermolecular distances are large in comparison with molecular sizes the intermolecular interactions are negligible hard condensed matter: strong intermolecular interactions, long-range orientational and positional order molecular liquids and soft matter: energy-entropy interplay. The free-energy scale: k B T. The relevant properties of the system are determined by the interplay of the various temporal and spatial scales involved. Universität Stuttgart - SFB 76, Pforzheim-Hohenwart, September 8-, 28 p. 4/59
5 Length-scales Wide range of characteristic length-scales: bond length:.nm persistence length of polymers: nm radius of gyration: nm macroscopic domains size: µm Universität Stuttgart - SFB 76, Pforzheim-Hohenwart, September 8-, 28 p. 5/59
6 Time-scales Wide range of characteristic time-scales: covalent bond vibrations: 5 s angle bending: 2 5 s atom dynamics governed by torsional and van der Waals interactions: 4 5 s atom dynamics governed by long-range Coulomb interactions: 2 s relaxation times of macromolecular segments: 2 9 s maximal macromolecular relaxation times: s Universität Stuttgart - SFB 76, Pforzheim-Hohenwart, September 8-, 28 p. 6/59
7 Multiscale modeling Universita t Stuttgart - SFB 76, Pforzheim-Hohenwart, September 8-, 28 p. 7/59
8 Adaptive resolution simulation Motivation: to treat in a simulation only as many degrees of freedom (DOFs) as absolutely necessary for the problem considered. Method: AdResS (Adaptive Resolution Scheme) allows for an dynamical switching between the atomistic and mesoscopic levels of detail = on-the-fly changing of the number of DOFs tailor-made for molecular systems where spatially localized domains with required atomistic resolution exchange particles with the remainder of the system sufficiently described on the mesoscopic scale. Results: accurately reproduces the statistical properties of the reference all-atom system. Universität Stuttgart - SFB 76, Pforzheim-Hohenwart, September 8-, 28 p. 8/59
9 Molecular Dynamics (MD) simulation V F Rauschen Noise Langevin Thermostat t "hard core" Damping Reibung V r F = m. a Electrostatic Elektrostatik Periodic Periodische Boundary Randbedingungen Conditions Geometry r 3 P M Algorithm (Mesh Ewald) hexagonal cell Universität Stuttgart - SFB 76, Pforzheim-Hohenwart, September 8-, 28 p. 9/59
10 MD simulation All-Atom MD simulation: allows to study processes at the atomic level of detail is often incapable to bridge a gap between a wide range of length and time scales involved in molecular systems Mesoscopic MD simulation: reduces the number of DOFs by retaining only those that are relevant for the property of interest = longer length and time scales can be reached specific chemical details are usually lost in the coarse-graining procedure Combining the best from both approaches: Hybrid Adaptive MD Schemes Universität Stuttgart - SFB 76, Pforzheim-Hohenwart, September 8-, 28 p. /59
11 Hybrid atomistic/mesoscopic liquid MP, L. Delle Site, K. Kremer, J. Chem. Phys. 23, 2246, 25. MP, L. Delle Site, K. Kremer, Phys. Rev. E 73, 667, 26. Universität Stuttgart - SFB 76, Pforzheim-Hohenwart, September 8-, 28 p. /59
12 Changing number of degrees of freedom A tetrahedral molecule has a defined spatial orientation and 3N = 2 DOFs: 3 translational 3 rotational 3N 6 = 6 vibrational One particle mesoscopic molecule has no defined spatial orientation and only 3 translational DOFs. MP, L. Delle Site, K. Kremer, J. Chem. Phys. 23, 2246, 25. MP, L. Delle Site, K. Kremer, Annu. Rev. Phys. Chem. 59, 545, 28. Universität Stuttgart - SFB 76, Pforzheim-Hohenwart, September 8-, 28 p. 2/59
13 Geometrically induced phase transition In thermodynamical equilibrium,boundary conditions analogous to to two-phase coexistence must be satisfied: µ ex = µ cg, p ex = p cg, T ex = T cg. The rotational and vibrational parts of the free energy can be viewed as thelatent heat, which is supplied or taken by the thermostat, at this transition. Universität Stuttgart - SFB 76, Pforzheim-Hohenwart, September 8-, 28 p. 3/59
14 Transition region F(x) A B d +d X Molecules in A and B are physically identical but differently represented. The number of DOFs is n = n(x) with: n A = const A ; n B = const B ; and n = n(x) The system is in equilibrium which implies: lim F A(x) x d = lim F B(x) x x d + = = x lim n A(x) x d = lim n B(x) x x d + = x Universität Stuttgart - SFB 76, Pforzheim-Hohenwart, September 8-, 28 p. 4/59
15 Weighting Function = Order Parameter w.5..2 r/a The values w = and w = correspond to the atomistic and coarse-grained regions, respectively, while the values < w < correspond to the transition (hyb) regime. Universität Stuttgart - SFB 76, Pforzheim-Hohenwart, September 8-, 28 p. 5/59
16 AdResS: Linking atomic and mesoscopic length-scales AdResS consists of two main steps:. Derive the effective pair potential U cm between coarse-grained molecules on the basis of the reference all-atom system. 2. Couple the atomistic and mesoscopic scales: F αβ = w(x α )w(x β )F atom αβ +[ w(x α )w(x β )]F cm αβ, where F atom αβ = iα,jβ F atom iαjβ is the sum of all pair interactions between explicit atoms of molecules α and β and F atom iαjβ F cm αβ = Uatom r iαjβ, = Ucm R αβ. Universität Stuttgart - SFB 76, Pforzheim-Hohenwart, September 8-, 28 p. 6/59
17 May the Force be with you One must interpolate the forces and not the interaction potentials if the Newton s Third Law is to be satisfied! MP, K. Kremer, L. Delle Site, J. Phys. A: Math. Theor. 4, F28, 27. Universität Stuttgart - SFB 76, Pforzheim-Hohenwart, September 8-, 28 p. 7/59
18 Continuous change of the phase space dimensionality The switching procedure implies that in the transition regime, where < w(x) <, we deal with fractional DOFs, i.e., by switching on/off a DOF we continously change the dimensionality of the phase space. To proper describe this situation we resort here to the fractional calculus. Let us consider a given w(x) = α and use α as a variable parameter. The infinitesimal volume element of the fractional configurational space is defined as dv α = d α xγ(α/2)/2π α/2 Γ(α) = x α dx/γ(α) = dx α /αγ(α) where the positive real parameter α denotes the order of the fractional coordinate differential. MP, K. Kremer, L. Delle Site, Phys. Rev. E 75, 77, 27. Universität Stuttgart - SFB 76, Pforzheim-Hohenwart, September 8-, 28 p. 8/59
19 The extension of equipartition theorem to non-integer DOFs For the fractional quadratic DOF Θ with the weight w = α we can write the partition function as: exp( βf α ) = C exp( βαp 2 Θ /2)dV α = = 2C exp( βαp 2 Θ /2) p Θ α dp Θ Γ(α) = = 2α/2 CΓ(α/2) Γ(α) α α/2 β α/2 β α/2. K α = d(βf α) dβ = α 2β = αk BT 2. In equilibriumt A = T B = T = T and thus: n α α. MP, K. Kremer, L. Delle Site, Phys. Rev. E 75, 77, 27. MP, K. Kremer, L. Delle Site, J. Phys. A: Math. Theor. 4, F28, 27. Universität Stuttgart - SFB 76, Pforzheim-Hohenwart, September 8-, 28 p. 9/59
20 Langevin thermostat Equation of motion for thei-th particle: Fluctuation-dissipation theorem: m i d 2 r i dt 2 = F i m i Γ dr i dt +W i(t). < W i (t) W j (t ) >= δ ij δ(t t )6 m i m j k B TΓ. Universität Stuttgart - SFB 76, Pforzheim-Hohenwart, September 8-, 28 p. 2/59
21 Mapping of structural properties Effective pair potential that reproduces the structure of the all-atom system as closely as possible is determined using the RDF cm of the reference all-atom system via the potential of mean force (PMF ) as U cm (r) PMF(r) = k B T logg cm ex (r), g cm ex (r) is the all-atom RDF cm and U cm (r) is the effective potential. Effective potential is: in general temperature and density dependent softer than the interatomic potentials Universität Stuttgart - SFB 76, Pforzheim-Hohenwart, September 8-, 28 p. 2/59
22 Iterative Boltzmann inversion In the limit ρ = : U eff (r) = F(r) = k B T lng(r). For systems with ρ > the above relation is used as an initial approximation in the iteration scheme: U eff i+ (r) = Ueff i (r)+k B T ln g i (r) g target (r) Ramp pressure correction: U lin = A( r r cut ). D. Reith, M. Pütz, F. Müller-Plathe, J. Comput. Chem., 24, 624, 23. Universität Stuttgart - SFB 76, Pforzheim-Hohenwart, September 8-, 28 p. 22/59
23 Tetrahedral liquid: effective potential Utab cm U cm ex cg U RDF r r (c) Tabulated effective potential. (d) RDF cm. Universität Stuttgart - SFB 76, Pforzheim-Hohenwart, September 8-, 28 p. 23/59
24 Equation of state 6 5 ex cg ex ex cg 4.5 p 3 2 RDF ρ r (e) EOS forρ =.75 andt =. (f) RDF cm. Universität Stuttgart - SFB 76, Pforzheim-Hohenwart, September 8-, 28 p. 24/59
25 Density profile profile.4.3 ex cg.4.3 ex cg ρ/ρ ρ/ρ r r 2 4 (g) Radius of the explicit regime: (h) Radius of the explicit regime: r = 6.. r =.. Universität Stuttgart - SFB 76, Pforzheim-Hohenwart, September 8-, 28 p. 25/59
26 Artifacts of the method p * RDF ex ex-cg:w=const=/2 U * PMF ex PMF ex-cg(w=const=/2) r * w (i) Pressure r * (j) RDF cm. Universität Stuttgart - SFB 76, Pforzheim-Hohenwart, September 8-, 28 p. 26/59
27 Interface pressure correction (ic) F cm αβ = s[w(r α )w(r β )]F cm αβ o + ( s[w(r α )w(r β )] ) F cm αβ ic. s [,] is defined as s[x] = 4( x 2 )2, where s[] =, s[] = and s[/4] = Universität Stuttgart - SFB 76, Pforzheim-Hohenwart, September 8-, 28 p. 27/59
28 System of hybrid molecules with w = /2 RDF U Uic cm U cm.6.8 r RDF U Utab cm ic Utab cm r.5 2 ex ex cg(w = /2) ex cg(w = /2) ic r ex ex cg(w = /2) ex cg(w = /2) ic r (k) ρ =.. (l) ρ =.75. Universität Stuttgart - SFB 76, Pforzheim-Hohenwart, September 8-, 28 p. 28/59
29 Density profile with ic correction.4.3 ex cg ic.4.3 ex cg ic ρ/ρ ρ/ρ r r 2 4 (m)ρ =.. (n)ρ =.75. Universität Stuttgart - SFB 76, Pforzheim-Hohenwart, September 8-, 28 p. 29/59
30 Macromolecule in solvent MP, L. Delle Site, K. Kremer, J. Chem. Phys. 26, 3492, 27. Universität Stuttgart - SFB 76, Pforzheim-Hohenwart, September 8-, 28 p. 3/59
31 Maximal monomer distance r = 7. r =. r = 2. rmax t Universität Stuttgart - SFB 76, Pforzheim-Hohenwart, September 8-, 28 p. 3/59
32 Static properties of the solvent ex N= ex N=2 ex N= ex N=3 ex cg N=3 ex cg ic N=3 RDFcm.5 RDFcm r r Universität Stuttgart - SFB 76, Pforzheim-Hohenwart, September 8-, 28 p. 32/59
33 Static properties of the polymer.4.2 ex N=3 ex cg N=3 ex cg ic N=3 q 2 S(q) ρ/ρ q ex N=3 ex cg N=3 ex cg icn= r S(q) = N exp(iq (r i r j )) ij S(q) q /ν q 2 S(q) q 2 /ν,ν.57 Universität Stuttgart - SFB 76, Pforzheim-Hohenwart, September 8-, 28 p. 33/59
34 Dynamic properties of the polymer MSD MSD.. < r 2 > < r 2 exn= > < r 2 exn=2 > < R 2 exn=3 > < R 2 exn= > < R 2 exn=2 > exn=3 t.. < r 2 > < r 2 ex cg/icn= > < r 2 ex cg/icn=2 > < R 2 ex cg/icn=3 > < R 2 ex cg/icn= > < R 2 ex cg/icn=2 > ex cg/icn=3 t D = 6 lim t R i (t) R i () 2 t = 6 lim t R 2 t r 2 = (r i (t) r i ()) 2 t 2/z,2/z.67(N = 3) Universität Stuttgart - SFB 76, Pforzheim-Hohenwart, September 8-, 28 p. 34/59
35 Diffusion constant of the chain N 2 3 D(ex) D(ex-cg) D(ex-cg ic ) D bulkex =.36,D bulkcg =.57 D polymer D solvent Universität Stuttgart - SFB 76, Pforzheim-Hohenwart, September 8-, 28 p. 35/59
36 Liquid water The simulation speed-up is 7 2 compared to atomistic simulations. MP, S. Matysiak, L. Delle Site, K. Kremer, C. Clementi, J. Phys.: Condens. Matter, 9, 2922, 27. Universität Stuttgart - SFB 76, Pforzheim-Hohenwart, September 8-, 28 p. 36/59
37 Electrostatics: Reaction field method The electrostatic forces interactions are described using the Reaction field (RF) method: F atom C i αj β (r iα j β ) = e i α e jβ 4πǫ [ r 3 i α j β R 3 c 2(ǫ RF ) +2ǫ RF ] r iα j β. The RF is suitable to be used with AdResS because: it is pairwise like AdResS it must also be applied with a thermostat Universität Stuttgart - SFB 76, Pforzheim-Hohenwart, September 8-, 28 p. 37/59
38 Weighting function for flat geometry w x/a Universität Stuttgart - SFB 76, Pforzheim-Hohenwart, September 8-, 28 p. 38/59
39 Static properties q = j=k=j+ ( cosψ jk + ) 2 3 Universität Stuttgart - SFB 76, Pforzheim-Hohenwart, September 8-, 28 p. 39/59
40 Interface effect of the cg water transition all-atom transition The transition regime neutralizes the interface effect of the cg water = the structure of water in the explicit regime is the same as in the bulk. Universität Stuttgart - SFB 76, Pforzheim-Hohenwart, September 8-, 28 p. 4/59
41 Diffusion across the transition regime I. coarse-grain transition all-atom Universität Stuttgart - SFB 76, Pforzheim-Hohenwart, September 8-, 28 p. 4/59
42 Diffusion coefficient across the simulation box S. Matysiak, C. Clementi, MP, K. Kremer, L. Delle Site, J. Chem. Phys 28, 2453, 28. Universität Stuttgart - SFB 76, Pforzheim-Hohenwart, September 8-, 28 p. 42/59
43 Position dependent Langevin thermostat The Langevin equation with a position dependent coefficient Γ(x) can be written as: m i dv i /dt = F i m i Γ(x)v i +R i (x,t) () where R i (x,t) is: R i (x,t) =, (2) R i (x,t )R j (x,t 2 ) = 2Γ(x)m i ktδ(t t 2 )δ ij (3) Γ(x) = Γ cg if x.6 αx+β if.6 < x. (4) This choice provides a simple interpolation between the two limit values of Γ(.6) = Γ() = Γ cg = 5ps and Γ() = Γ all atom = 5ps. The parameters α and β are 25ps and 3ps, respectively. Universität Stuttgart - SFB 76, Pforzheim-Hohenwart, September 8-, 28 p. 43/59
44 Diffusion coefficient in the hybrid system Universität Stuttgart - SFB 76, Pforzheim-Hohenwart, September 8-, 28 p. 44/59
45 Diffusion coefficient across the simulation box II. (w) Position dependent thermostat (x) Regular thermostat Universität Stuttgart - SFB 76, Pforzheim-Hohenwart, September 8-, 28 p. 45/59
46 Transverse DPD Thermostat The variation of the dissipative particle dynamics (DPD) thermostat includes the damping of the perpendicular components of the relative velocity, yet keeping the advantages of conserving Galilei invariance and within our error bar also hydrodynamics. It allows for controlling transport properties of molecular fluids. C. Junghans, MP, K. Kremer, Soft Matter 4, 56, 28. Universität Stuttgart - SFB 76, Pforzheim-Hohenwart, September 8-, 28 p. 46/59
47 AdResS: Conclusions Adaptive Resolution MD simulation: Changing resolution is formally equivalent to a phase transition latent heat. For a smooth variation of the resolution we introduce a transition regime. The temperature in the transition region can be obtained by extending the equipartition theorem to non-integer dimensions. Hybrid method AdResS: Allows for a dynamical switching of the spatial resolution. We treat only as many DOFs as absolutely necessary for the problem considered. AdResS was so far applied to MD simulations of a simple tetrahedral liquid, a macromolecule in an explicit solvent, and liquid water at standard conditions. The simulation speed-up for liquid water is 7 2 compared to atomistic simulations. Review: MP, L. Delle Site, K. Kremer, Annu. Rev. Phys. Chem. 59, 545, 28. Universität Stuttgart - SFB 76, Pforzheim-Hohenwart, September 8-, 28 p. 47/59
48 AdResS: Future work Universität Stuttgart - SFB 76, Pforzheim-Hohenwart, September 8-, 28 p. 48/59
49 Concurrent triple-scale simulation Motivation: to simplify the model to the largest extend possible while keeping all the necessary details where this is required Method: Triple-scale AdResS-HybridMD scheme is a combination of two dual-scale models: a particle-based Adaptive Resolution Scheme (AdResS), which couples the atomic and mesoscopic scales, and a hybrid continuum-molecular dynamics scheme (HybridMD) covers the length-scales ranging from the micro- to macro-scale successfully sorts out the problem of large molecule insertion in the hybrid particle-continuum simulations of molecular liquids opens up the possibility to perform efficient grand-canonical molecular dynamics simulations of truly open molecular liquid systems Results: the structural and dynamical properties of the liquid are accurately captured Universität Stuttgart - SFB 76, Pforzheim-Hohenwart, September 8-, 28 p. 49/59
50 Triple-scale model R. Delgado Buscalioni, K. Kremer, MP, J. Chem. Phys. 28, 4, 28. Universität Stuttgart - SFB 76, Pforzheim-Hohenwart, September 8-, 28 p. 5/59
51 HybridMD: Coupling particle-based and continuum descriptions The hybrid particle-continuum scheme (HybridMD) is designed to connect the dynamics of a molecular domain with that obtained from a continuum description of the surrounding fluid flow. The method is based on flux-exchange. The system is divided in (at least) two domains, described via classical molecular dynamics (MD) and continuum fluid dynamics (CFD), i.e., solving the Navier-Stokes equations. The MD and CFD domains share one unique hybrid interface, H: Flux balance implies the conservation of mass and momentum across H. G. De Fabritiis, R. Delgado Buscalioni, P. Coveney, Phys. Rev. Lett 97, 345, 26. R. Delgado Buscalioni, G. De Fabritiis, Phys. Rev. E 76, 3679, 27. Universität Stuttgart - SFB 76, Pforzheim-Hohenwart, September 8-, 28 p. 5/59
52 AdResS-HybridMD: Combined scheme HybridMD setup AdResS setup 3σ C cg 2.5σ 2.5σ hyb ex hyb cg P P C wall e H B H H x B e H 6 wall Domain decomposition of the combined scheme. The top part of the figure shows the location of the fluid model layers (cg, hyb and ex) within the HybridMD setup. The bottom part of the figure shows the set of control cells used in the HybridMD setup. R. Delgado Buscalioni, K. Kremer, MP, J. Chem. Phys. 28, 4, 28. Universität Stuttgart - SFB 76, Pforzheim-Hohenwart, September 8-, 28 p. 52/59
53 Weighting function w x/l Universität Stuttgart - SFB 76, Pforzheim-Hohenwart, September 8-, 28 p. 53/59
54 Molecular density profile molecular density molecular density B cg cg hyb ex hyb cg x (σ) ex H hyb cg H H cg B cg B buffer (a) (b) x(σ) (a) ρ m =.σ 3. (b) ρ m =.75σ 3. Universität Stuttgart - SFB 76, Pforzheim-Hohenwart, September 8-, 28 p. 54/59
55 Radial distribution functions 2.5 ex+hyb (HybridMD) ex+hyb+cg (HybridMD) ex (PBC) RDF r (σ) RDF cm s of the liquid in the atomistic and transition domains (ex+hyb) and in the total molecular region (ex+hyb+cg) of the triple-scale model together with the reference RDF cm of the all-atom system (ex(pbc)) at ρ =.75/σ 3. Universität Stuttgart - SFB 76, Pforzheim-Hohenwart, September 8-, 28 p. 55/59
56 Couette flow.2 cg cg hyb ex hyb cg cg velocity (σ/τ) B H H x(σ) B Velocity profile at the particle region of an hybrid simulation of a Couette flow. Universität Stuttgart - SFB 76, Pforzheim-Hohenwart, September 8-, 28 p. 56/59
57 Stokes flow velocity (σ/τ) velocity (σ/τ) cell #5 cell # time (τ) cell #8 cell # time (τ) Velocity in the y-direction at some selected cells in a hybrid simulation of a Stokes flow. Universität Stuttgart - SFB 76, Pforzheim-Hohenwart, September 8-, 28 p. 57/59
58 AdResS-HybridMD: Conclusions AdResS-HybridMD Scheme: We performed a triple-scale simulation of a molecular liquid. Length scales from the micro- to macro-scale are concurrently coupled. The method allows us to perform efficient grand-canonical molecular dynamics simulations of molecular liquids. Future work: Application to realistic systems, e.g., liquid water. Universität Stuttgart - SFB 76, Pforzheim-Hohenwart, September 8-, 28 p. 58/59
59 Triple-scale liquid water URL: praprot/ Universität Stuttgart - SFB 76, Pforzheim-Hohenwart, September 8-, 28 p. 59/59
Coupling Atomistic and Continuum Hydrodynamics
Coupling Atomistic and Continuum Hydrodynamics Matej Praprotnik praprot@cmm.ki.si Laboratory for Molecular Modeling National Institute of Chemistry Ljubljana, Slovenia & Department of Physics Faculty of
More informationHybrid Atomistic-Continuum Methods for Dense Liquids
Hybrid Atomistic-Continuum Methods for Dense Liquids Matej Praprotnik Laboratory for Molecular Modeling National Institute of Chemistry Ljubljana Slovenia TUD Autumn School 2009, Darmstadt, September 24-25,
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 informationMultiscale Simulation of Soft Matter: From Scale Bridging to Adaptive Resolution
Annu Rev Phys Chem 200859:545-571 Downloaded from arjournalsannualreviewsorg by Dr Matej Praprotnik on 04/08/08 For personal use only Annu Rev Phys Chem 2008 59:545 71 First published online as a Review
More informationAtomistic, mesoscopic and continuum hydrodynamics
Atomistic, mesoscopic and continuum hydrodynamics coupling liquid models with different resolution Rafael Delgado-Buscalioni Departamento de Fisica Teorica de la Materia Condensada Universidad Autonoma
More informationAdaptive Resolution Simulations: Towards Open Systems Molecular Dynamics Simulations. K. Kremer Max Planck Institute for Polymer Research, Mainz
Adaptive Resolution Simulations: Towards Open Systems Molecular Dynamics Simulations K. Kremer Max Planck Institute for Polymer Research, Mainz ACMAC Crete * June, 30, 2011 COWORKERS: Luigi Delle Site
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 informationA Control Algorithm for Multiscale Simulations of Liquid Water
A Control Algorithm for Multiscale Simulations of Liquid Water Evangelos M. Kotsalis and Petros Koumoutsakos Computational Science, ETH Zürich, CH-8092, Switzerland petros@ethz.ch Abstract. We present
More informationA multiscale link between atomistic and continuum hydrodynamics via adaptive coarse-graining.
A multiscale link between atomistic and continuum hydrodynamics via adaptive coarse-graining. Rafael Delgado-Buscalioni Departamento de Fisica Teorica de la Materia Condensada Universidad Autonoma de Madrid
More informationTransport properties controlled by a thermostat: An extended dissipative particle dynamics thermostat
PAPER www.rsc.org/softmatter Soft Matter Transport properties controlled by a thermostat: An extended dissipative particle dynamics thermostat Christoph Junghans, a Matej Praprotnik ab and Kurt Kremer
More 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 informationSupporting information
Electronic Supplementary Material ESI for Physical Chemistry Chemical Physics. This journal is the Owner Societies 2018 Supporting information Understanding three-body contributions to coarse-grained force
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 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 informationTools for multiscale simulation of liquids using open molecular dynamics
Tools for multiscale simulation of liquids using open molecular dynamics Rafael Delgado-Buscalioni PLEASE CITE as: "Tools for..." in Numerical Analysis and Multiscale Computations, Series Lect. Notes Comput.
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 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 informationMultiscale Materials Modeling
Multiscale Materials Modeling Lecture 10 Adaptive Resolution Fundamentals of Sustainable Technology These notes created by David Keffer, University of Tennessee, Knoxville, 2012. Reminder In the previous
More informationSmoothed Dissipative Particle Dynamics: theory and applications to complex fluids
2015 DPD Workshop September 21-23, 2015, Shanghai University Smoothed Dissipative Particle Dynamics: Dynamics theory and applications to complex fluids Marco Ellero Zienkiewicz Centre for Computational
More 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 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 informationDirect Modeling for Computational Fluid Dynamics
Direct Modeling for Computational Fluid Dynamics Kun Xu February 20, 2013 Computational fluid dynamics (CFD) is new emerging scientific discipline, and targets to simulate fluid motion in different scales.
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 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 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 informationMACROSCOPIC VARIABLES, THERMAL EQUILIBRIUM. Contents AND BOLTZMANN ENTROPY. 1 Macroscopic Variables 3. 2 Local quantities and Hydrodynamics fields 4
MACROSCOPIC VARIABLES, THERMAL EQUILIBRIUM AND BOLTZMANN ENTROPY Contents 1 Macroscopic Variables 3 2 Local quantities and Hydrodynamics fields 4 3 Coarse-graining 6 4 Thermal equilibrium 9 5 Two systems
More informationParticle-Simulation Methods for Fluid Dynamics
Particle-Simulation Methods for Fluid Dynamics X. Y. Hu and Marco Ellero E-mail: Xiangyu.Hu and Marco.Ellero at mw.tum.de, WS 2012/2013: Lectures for Mechanical Engineering Institute of Aerodynamics Technical
More informationNon-equilibrium phenomena and fluctuation relations
Non-equilibrium phenomena and fluctuation relations Lamberto Rondoni Politecnico di Torino Beijing 16 March 2012 http://www.rarenoise.lnl.infn.it/ Outline 1 Background: Local Thermodyamic Equilibrium 2
More informationMassively parallel molecular-continuum simulations with the macro-micro-coupling tool Neumann, P.; Harting, J.D.R.
Massively parallel molecular-continuum simulations with the macro-micro-coupling tool Neumann, P.; Harting, J.D.R. Published in: Hybrid particle-continuum methods in computational materials physics, 4-7
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 informationA Coupling Tool for Parallel Molecular Dynamics Continuum Simulations
A Coupling Tool for Parallel Molecular Dynamics Continuum Simulations ISPDC 2012 Philipp Neumann and Nikola Tchipev 29.06.2012 ISPDC 2012, 29.06.2012 1 Contents Motivation The Macro Micro Coupling Tool
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 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 informationStructure Analysis of Bottle-Brush Polymers: Simulation and Experiment
Jacob Klein Science 323, 47, 2009 Manfred Schmidt (Mainz) Structure Analysis of Bottle-Brush Polymers: Simulation and Experiment http://www.blueberryforest.com Hsiao-Ping Hsu Institut für Physik, Johannes
More informationHydrodynamics, Thermodynamics, and Mathematics
Hydrodynamics, Thermodynamics, and Mathematics Hans Christian Öttinger Department of Mat..., ETH Zürich, Switzerland Thermodynamic admissibility and mathematical well-posedness 1. structure of equations
More 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 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 informationHydrodynamics. Stefan Flörchinger (Heidelberg) Heidelberg, 3 May 2010
Hydrodynamics Stefan Flörchinger (Heidelberg) Heidelberg, 3 May 2010 What is Hydrodynamics? Describes the evolution of physical systems (classical or quantum particles, fluids or fields) close to thermal
More 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 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 informationFluid-Particles Interaction Models Asymptotics, Theory and Numerics I
Fluid-Particles Interaction Models Asymptotics, Theory and Numerics I J. A. Carrillo collaborators: T. Goudon (Lille), P. Lafitte (Lille) and F. Vecil (UAB) (CPDE 2005), (JCP, 2008), (JSC, 2008) ICREA
More information1. Thermodynamics 1.1. A macroscopic view of matter
1. Thermodynamics 1.1. A macroscopic view of matter Intensive: independent of the amount of substance, e.g. temperature,pressure. Extensive: depends on the amount of substance, e.g. internal energy, enthalpy.
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 informationMesoscale fluid simulation of colloidal systems
Mesoscale fluid simulation of colloidal systems Mingcheng Yang Institute of Physics, CAS Outline (I) Background (II) Simulation method (III) Applications and examples (IV) Summary Background Soft matter
More informationStatistical Thermodynamics and Monte-Carlo Evgenii B. Rudnyi and Jan G. Korvink IMTEK Albert Ludwig University Freiburg, Germany
Statistical Thermodynamics and Monte-Carlo Evgenii B. Rudnyi and Jan G. Korvink IMTEK Albert Ludwig University Freiburg, Germany Preliminaries Learning Goals From Micro to Macro Statistical Mechanics (Statistical
More informationThe Truth about diffusion (in liquids)
The Truth about diffusion (in liquids) Aleksandar Donev Courant Institute, New York University & Eric Vanden-Eijnden, Courant In honor of Berni Julian Alder LLNL, August 20th 2015 A. Donev (CIMS) Diffusion
More informationNanoscale Energy Transport and Conversion A Parallel Treatment of Electrons, Molecules, Phonons, and Photons
Nanoscale Energy Transport and Conversion A Parallel Treatment of Electrons, Molecules, Phonons, and Photons Gang Chen Massachusetts Institute of Technology OXFORD UNIVERSITY PRESS 2005 Contents Foreword,
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 informationAnomalous Transport and Fluctuation Relations: From Theory to Biology
Anomalous Transport and Fluctuation Relations: From Theory to Biology Aleksei V. Chechkin 1, Peter Dieterich 2, Rainer Klages 3 1 Institute for Theoretical Physics, Kharkov, Ukraine 2 Institute for Physiology,
More informationIntroduction to Molecular Dynamics
Introduction to Molecular Dynamics Dr. Kasra Momeni www.knanosys.com Overview of the MD Classical Dynamics Outline Basics and Terminology Pairwise interacting objects Interatomic potentials (short-range
More informationV (r,t) = i ˆ u( x, y,z,t) + ˆ j v( x, y,z,t) + k ˆ w( x, y, z,t)
IV. DIFFERENTIAL RELATIONS FOR A FLUID PARTICLE This chapter presents the development and application of the basic differential equations of fluid motion. Simplifications in the general equations and common
More 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 informationPHASE TRANSITIONS IN SOFT MATTER SYSTEMS
OUTLINE: Topic D. PHASE TRANSITIONS IN SOFT MATTER SYSTEMS Definition of a phase Classification of phase transitions Thermodynamics of mixing (gases, polymers, etc.) Mean-field approaches in the spirit
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 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 informationLong-Term Atomistic Simulation of Hydrogen Diffusion in Metals
Long-Term Atomistic Simulation of Hydrogen Diffusion in Metals K.G. Wang 1, M. and M. Ortiz 2 Universidad de Sevilla 1 Virginia Polytechnic Institute and State University 2 California Institute of Technology
More informationReaction Dynamics (2) Can we predict the rate of reactions?
Reaction Dynamics (2) Can we predict the rate of reactions? Reactions in Liquid Solutions Solvent is NOT a reactant Reactive encounters in solution A reaction occurs if 1. The reactant molecules (A, B)
More informationAdvantages and challenges in coupling an ideal gas to atomistic models in adaptive resolution simulations
EPJ manuscript No. (will be inserted by the editor) Advantages and challenges in coupling an ideal gas to atomistic models in adaptive resolution simulations Karsten Kreis 1,2,a, Aoife C. Fogarty 1,b,KurtKremer
More informationExperimental Soft Matter (M. Durand, G. Foffi)
Master 2 PCS/PTSC 2016-2017 10/01/2017 Experimental Soft Matter (M. Durand, G. Foffi) Nota Bene Exam duration : 3H ecture notes are not allowed. Electronic devices (including cell phones) are prohibited,
More informationElementary Lectures in Statistical Mechanics
George DJ. Phillies Elementary Lectures in Statistical Mechanics With 51 Illustrations Springer Contents Preface References v vii I Fundamentals: Separable Classical Systems 1 Lecture 1. Introduction 3
More informationCoupling the Level-Set Method with Variational Implicit Solvent Modeling of Molecular Solvation
Coupling the Level-Set Method with Variational Implicit Solvent Modeling of Molecular Solvation Bo Li Math Dept & CTBP, UCSD Li-Tien Cheng (Math, UCSD) Zhongming Wang (Math & Biochem, UCSD) Yang Xie (MAE,
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 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 information5. Systems in contact with a thermal bath
5. Systems in contact with a thermal bath So far, isolated systems (micro-canonical methods) 5.1 Constant number of particles:kittel&kroemer Chap. 3 Boltzmann factor Partition function (canonical methods)
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 informationA MOLECULAR DYNAMICS STUDY OF POLYMER/GRAPHENE NANOCOMPOSITES
A MOLECULAR DYNAMICS STUDY OF POLYMER/GRAPHENE NANOCOMPOSITES Anastassia N. Rissanou b,c*, Vagelis Harmandaris a,b,c* a Department of Applied Mathematics, University of Crete, GR-79, Heraklion, Crete,
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 informationImproved Adaptive Resolution Molecular Dynamics Simulation
Improved Adaptive Resolution Molecular Dynamics Simulation Iuliana Marin, Virgil Tudose, Anton Hadar, Nicolae Goga, Andrei Doncescu To cite this version: Iuliana Marin, Virgil Tudose, Anton Hadar, Nicolae
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 informationCollective behavior, from particles to fields
978-0-51-87341-3 - Statistical Physics of Fields 1 Collective behavior, from particles to fields 1.1 Introduction One of the most successful aspects of physics in the twentieth century was revealing the
More informationmacroscopic view (phenomenological) microscopic view (atomistic) computing reaction rate rate of reactions experiments thermodynamics
Rate Theory (overview) macroscopic view (phenomenological) rate of reactions experiments thermodynamics Van t Hoff & Arrhenius equation microscopic view (atomistic) statistical mechanics transition state
More 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 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 informationStatistical Mechanics and Thermodynamics of Small Systems
Statistical Mechanics and Thermodynamics of Small Systems Luca Cerino Advisors: A. Puglisi and A. Vulpiani Final Seminar of PhD course in Physics Cycle XXIX Rome, October, 26 2016 Outline of the talk 1.
More informationThermodynamics of Solids: Harmonic and Quasi-harmonic Approximations
Thermodynamics of Solids: Harmonic and Quasi-harmonic Approximations, USA, July 9-14, 2017 Alessandro Erba Dipartimento di Chimica, Università di Torino (Italy) alessandro.erba@unito.it 2017 Outline -
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 informationElectronic Supplementary Information
Electronic Supplementary Material (ESI) for Physical Chemistry Chemical Physics. This journal is the Owner Societies 2014 Electronic Supplementary Information Computational investigation of structural
More informationAnnual Report for Research Work in the fiscal year 2006
JST Basic Research Programs C R E S T (Core Research for Evolutional Science and Technology) Annual Report for Research Work in the fiscal year 2006 Research Area : High Performance Computing for Multi-scale
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-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 informationMOLECULAR DYNAMICS SIMULATION OF HETEROGENEOUS NUCLEATION OF LIQUID DROPLET ON SOLID SURFACE
MOLECULAR DYNAMICS SIMULATION OF HETEROGENEOUS NUCLEATION OF LIQUID DROPLET ON SOLID SURFACE Tatsuto Kimura* and Shigeo Maruyama** *Department of Mechanical Engineering, The University of Tokyo 7-- Hongo,
More information(Crystal) Nucleation: The language
Why crystallization requires supercooling (Crystal) Nucleation: The language 2r 1. Transferring N particles from liquid to crystal yields energy. Crystal nucleus Δµ: thermodynamic driving force N is proportional
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 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 informationChapter 9: Differential Analysis of Fluid Flow
of Fluid Flow Objectives 1. Understand how the differential equations of mass and momentum conservation are derived. 2. Calculate the stream function and pressure field, and plot streamlines for a known
More informationComputer Simulation of Shock Waves in Condensed Matter. Matthew R. Farrow 2 November 2007
Computer Simulation of Shock Waves in Condensed Matter Matthew R. Farrow 2 November 2007 Outline of talk Shock wave theory Results Conclusion Computer simulation of shock waves Shock Wave Theory Shock
More informationsimulations Extracting coarse-grained elastic behavior of capsid proteins from molecular dynamics Outline Feb. 2009
Extracting coarse-grained elastic behavior of capsid proteins from molecular dynamics simulations Stephen D. Hicks Christopher L. Henley Cornell University Gordon Conf. on Physical Virology Feb. 2009 Outline
More informationContinuum Model of Avalanches in Granular Media
Continuum Model of Avalanches in Granular Media David Chen May 13, 2010 Abstract A continuum description of avalanches in granular systems is presented. The model is based on hydrodynamic equations coupled
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 informationFluctuations and continuity in particle-continuum hybrid simulations of unsteady flows based on flux-exchange
EUROPHYSICS LETTERS 15 March 25 Europhys. Lett., 69 (6), pp. 959 965 (25) DOI: 1.129/epl/i24-1431-y Fluctuations and continuity in particle-continuum hybrid simulations of unsteady flows based on flux-exchange
More informationOnce Upon A Time, There Was A Certain Ludwig
Once Upon A Time, There Was A Certain Ludwig Statistical Mechanics: Ensembles, Distributions, Entropy and Thermostatting Srinivas Mushnoori Chemical & Biochemical Engineering Rutgers, The State University
More informationI. Collective Behavior, From Particles to Fields
I. Collective Behavior, From Particles to Fields I.A Introduction The object of the first part of this course was to introduce the principles of statistical mechanics which provide a bridge between the
More informationON ALGORITHMS FOR BROWNIAN DYNAMICS COMPUTER SIMULATIONS
COMPUTATIONAL METHODS IN SCIENCE AND TECHNOLOGY 4,35-42 (1998) ON ALGORITHMS FOR BROWNIAN DYNAMICS COMPUTER SIMULATIONS ARKADIUSZ C. BRAŃKA Institute of Molecular Physics, Polish Academy of Sciences, Smoluchowskiego
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 informationDiffusive Transport Enhanced by Thermal Velocity Fluctuations
Diffusive Transport Enhanced by Thermal Velocity Fluctuations Aleksandar Donev 1 Courant Institute, New York University & Alejandro L. Garcia, San Jose State University John B. Bell, Lawrence Berkeley
More informationChapter 1 Direct Modeling for Computational Fluid Dynamics
Chapter 1 Direct Modeling for Computational Fluid Dynamics Computational fluid dynamics (CFD) is a scientific discipline, which aims to capture fluid motion in a discretized space. The description of the
More informationMultiscale simulations of complex fluid rheology
Multiscale simulations of complex fluid rheology Michael P. Howard, Athanassios Z. Panagiotopoulos Department of Chemical and Biological Engineering, Princeton University Arash Nikoubashman Institute of
More informationCECAM Tutorial Stuttgart, Oct. 2012
CECAM Tutorial Stuttgart, Oct. 2012 09.10.2012 ESPResSo++ CECAM Tutorial 1 / 40 ESPResSo++ ESPResSo++ is an open-source simulation package designed to perform molecular dynamics and monte carlo simulations
More informationIntermolecular Forces in Density Functional Theory
Intermolecular Forces in Density Functional Theory Problems of DFT Peter Pulay at WATOC2005: There are 3 problems with DFT 1. Accuracy does not converge 2. Spin states of open shell systems often incorrect
More information