Dynamic modelling of morphology development in multiphase latex particle Elena Akhmatskaya (BCAM) and Jose M. Asua (POLYMAT) June
|
|
- Bertina Kelley
- 6 years ago
- Views:
Transcription
1 Dynamic modelling of morphology development in multiphase latex particle Elena Akhmatskaya (BCAM) and Jose M. Asua (POLYMAT) June
2 Publications (background and output) J.M. Asua, E. Akhmatskaya Dynamic modelling of morphology development in multiphase latex particles Akhmatskaya E., Asua J.M., Dynamic modeling of the morphology of latex particles with in situ formation of graft copolymer, Journal of Polymer Science - A - Polymer Chemistry Edition, Vol. 50, 7, pp (2012) Akhmatskaya E., Asua J.M., Dynamic modeling of the morphology of multiphase waterborne polymer particles, preprint Akhmatskaya E., Reich S., New Hybrid Monte Carlo Methods for Efficient Sampling: from Physics to Biology and Statistics, Progress in Nuclear Science and Technology, Vol. 2, pp (2012). June 2012
3 Outline Motivation Challenges Novel model Size reduction Dynamics Effective sampling Mimicking the technological conditions Prediction of the dynamic development of the morphology for waterborne systems Model testing Multiphase waterborne polymer Graft copolymer Conclusions
4 Explanation of Terms Phase domain: Region of a material that is uniform in chemical composition and physical state. Multiphase polymer: polymer comprising phaseseparated domains. Morphology: Shape or form of phase domains observed, often by microscopy or scattering techniques in polymers, polymer blends, composites, and crystals. Multiphase Equilibrium Morphologies (clays: Laponite) Multiphase Non-equilibrium Morphologies Jönsson et al. Macromolecules,1991, 24, 126 Negrete-Herrera et al. Macromol. Rapid Commun. 2007, 28, 1567
5 Motivation Waterborne polymers applications: synthetic rubber, paints, adhesives, cosmetics, leather treatments, diagnostic tests, drug delivery Multiphase particles provide advantages over particles with uniform compositions Performance of multiphase particles depends on particle morphology Assembling multiphase particles is of great practical interest Current status: The synthesis of new morphologies is time and resources consuming: largely relies on heuristic knowledge Little known: theoretical work is mainly restricted to 2-phase systems number of morphologies of 3-phase polymers is unknown no general methodology predicting the dynamic development of morphologies for multiphase waterborne systems is available. Objective: to develop the first ever computationally feasible modeling approach for prediction of dynamics of the particle morphology development in the multiphase waterborne systems Sundberg- Sundberg, J Appl Polym Sci, 1993, 47, 1277
6 Challenges Simulation scale: atomistic / coarsegrained to reproduce the properties of interest Large size systems: Practice: 200 nm particle (1500 polymer chains) surrounded by water Simulation: up to µm Slow (rare event) processes : Practice: hours Simulation: up to few seconds Dynamics is important: Practice: non-equilibrium morphologies are governed by kinetics Simulation: major enhanced samplers do not reproduce dynamics (MCMC type) Mimicking technological process: Practice: the complex multivariate process Simulation:?
7 Technological Polymerization Process Production of Composite Latex Particles Seeded emulsion polymerization Monomer(s) 1 SEED (Polymer 2) Inorganic Initial Charge Water Composite particles Emulsifier Final Product
8 Particle Morphology Development Phase I Monomer swollen particle Monomer converts to polymer Phase II Non-equilibrium morphology Equilibrium morphology process time Monomer + Polymer 1 process + Water time process time or Monomer + Polymer 1 + Inorganic + Water or More Phases process Monomer time Polymer 2 or N% Monomer Polymer 2 + (100% - N%) Monomer Copolymer 3 M% Polymer Copolymer 3 or More Phases
9 Simulation Model: Size Typical size of polymer particle: D = 200 nm Estimated minimal number of molecules:~ 10 7 not feasible Solution: particles instead of atoms Phase / inorganic particle: a polymer chain Water particle: water molecules packed into a particle of the diameter of the phase Number of particles for each species (phase / inorganic), N i : a number of polymer chains, ~ 1.5x10 3 Number of water particles, N waters : Water Polymer 1 N species N waters =! N i ~ N species " 1.5 " 10 3 Monomer i=1 Inorganic Simulation sphere: Wall Diameter = 2xD Surrounded by impenetrable structureless wall to keep phases & inorganic within polymer particle and to reproduce continuous aqueous phase out of the sphere
10 Simulation Model: Dynamics Langevin equation: Fluctuation-dissipation contributions (i) maintain temperature; (ii) mimic the impact of non- resolved finer details of an all-atom model on the coarse grained length and time scales. m d 2 r dr =!U "!m dt 2 dt + 2!k TmR( t ) B r is position, m is mass, γ is friction factor, k B is Boltzmann constant, T is temperature, U potential energy. t-time! controls the rate at which the phases move in the particle; increases with the internal viscosity of the particle:! =! 0 ( X ) 5! 0 = k T B mdf X is the conversion of monomer, Df is an effective diffusion coefficient that includes both Brownian and the interaction terms!c i = "# $ Df #C Material balance equation for polymerizing system!t i + R i suggests the way to mimic longer processes by choosing! ~ K!! exp,r i ~ K! R exp simulation _length i,k = real _ process _length C i is concentration and R i is rate of generation of phase i!r( t )" = 0;!R( t )R( t ' )" =!( t # t ' ).
11 Simulation Model: Sampler GSHMC: Generalized Shadow Hybrid Monte Carlo by Akmatskaya, Reich ( ) UK patent (2009), US patent (2011) Thermodynamically consistent implementation of constant-temperature molecular dynamics (MD) m d 2 r Based on MD and Monte Carlo (MC) Numerically integrates Newton equations dt =!U 2 and accepts new trajectories according to modified Metropolis criterion Improves sampling efficiency and retains dynamical information due to use of Modified Hamiltonians - asymptotic expansions in powers of the discretization parameter Δt, that are conserved by symplectic integrators to higher accuracy than true Hamiltonians; Introducing partial momentum update the noise vector is normally i.i.d. distributed; 0< φ π/2; φ = (2 Δt) ½ rigorous implementation of stochastic Langevin dynamics GSHMC is applied in current study after full conversion of monomer / polymers to speed up development of equilibrium morphologies (a factor of 3 in the presented simulations)
12 Interactions Lennard Jones (LJ): similar phases P1-P1; P2-P2; I-I; G-G; M2-M2; W-W; P1-M2; P1-G; P2-M2; P2-G; G-M2; W-W; W-wall U water!wall ( r ij ) = U LJ ( R c! r ij ) for r ij > R c! 2.5! 0 for r ij " R c! 2.5! Repulsive generalized soft sphere: dissimilar phases P1-P2; P1-I; P1-W; P1-wall; P2-I; P2-W; P2-wall; I-W, I-wall; M2-W; M2-wall; G-I; G-W; G-wall )! U r ( r ij ) =! " # r ij " $ & % 6 U poly!wall ( r ij ) = + * + # % $ & % " " %! $ # R c! r ' ij & 6 for r ij > R c! 2 1/ 6 " + 0 for r ij ( R c! 2 1/ 6 ", ε ij is depth potential well, σ is distance at which LJ = 0; R c is distance from the centre of the simulation sphere to the wall; ε is estimated from surface and interfacial tensions P-polymer; M-monomer; surface _tension = a!! bk B T I-Inorganic; G-graft; W-water;
13 Mimicking the Technological Conditions (Algorithm) Specify: Composition of the synthetized polymer: N p - number of monomer particles to be converted to polymer and (G) N g - number of polymer particles to be converted to graft copolymer; Conversion rate: number P p (P g ) of stages of length tp i (tg i ), i=1, P p/g required for polymerization and (G) graft formation; M p and (G) M g numbers of convertible particles per polymerization, graft formation (G) stage Phase I: Forming the monomer swollen initial particle 1. Randomly distribute simulated particles 2. Minimize initial structure using Steepest descent method 3. Equilibrate the system using Langevin Dynamics (LD) simulation Phase II: Polymerization / graft formation (G) 1. Move the simulation walls towards the centre of the simulation sphere to imitate pressure increase, applied between two phases 2. Choose randomly M p monomer particles and (G) M g polymer particles 3. Run LD for tp/g i, gradually changing LJ interaction of each of M p/g particle to the characteristic values of the new polymer / copolymer 4. Repeat Step 3 for P p / P g (G) times 5. Continue simulation using GSHMC with fixed values of LJ parameters and! value corresponding to! if non-equilibrium morphologies are of interest. Use large values of otherwise.!,!! " 2
14 Parameters Parameters of the Model Composition of the particle Composition of the formed polymer Rate of conversion Friction coefficients Parameters of the potentials Current Status: heuristic parameters based on the experimental data or obtained from expensive simulation Future Work: parameter optimization using Particle Swarm Optimization based approaches(with Universidad de Oviedo, Spain) Parameters of the Sampler Order of Time-step, Δt Number of MD steps, N MD Angle, ϕ Current Status: prior tuning of parameters through short simulation runs Future Work: development of adaptive GSHMC schemes
15 Results: 2-phase polymer systems All calculations are performed in parallel on i2basque, the computing cluster comprising 29 computing elements composed by 2 Intel QuadCore Xeon processors (GROMACS / GSHMC codes) Tested the model on 2-phase polymer systems All possible equilibrium morphologies are accurately reproduced Core-shell Inverted core-shell Hemi-spheres Sandwich
16 Results: Dynamic Modelling of 3-phase polymer particles Simulated dynamics of the development of particle morphology of 3-phase polymer systems Investigated effect of potential parameters, friction coefficients, conversion rates, composition of the particle on the morphology development X=0 ; T = 0 X=30; T=6000 X=50; T=10000 X=70; T=14000 X=90; T=18000 X=100; T=20000
17 Results: Effect of Polymerization Rate and Friction Factors on Morphology Development Faster polymerization rate X=0.3; T=1800 X=0.5; T=3000 X=0.7; T=4800 X=0.9; T=5400 X=1; T=6000 aging; T=15000 Faster polymerization rate. Increasing friction factor No se puede mostrar la imagen. Puede que su equipo no tenga suficiente memoria para abrir la imagen o que ésta esté dañada. Reinicie el equipo y, a continuación, abra el archivo de nuevo. Si sigue apareciendo la x roja, puede que tenga que borrar la imagen e insertarla de nuevo. X=0.3; T=1800 X=0.5; T=3000 X=0.7; T=4800 X=0.9; T=5400 X=1; T=6000 aging; T=25000
18 Results: Effect of Polymerization Rate and Friction Factors on Morphology Development (II) Run1: fast polymerization rate Run2: slow polymerization rate Run3: fast polymerization rate with increasing friction constants Faster polymerization rate leads to faster convergence to equilibrium morphologies Large friction constants prevent development of equilibrium morphologies
19 Results: In Situ Formation of Graft Copolymer Simulated the dynamics of of the development of particle morphology of composite waterborne systems in which a graft copolymer is produced in situ during the process Investigated effect of potential parameters, friction coefficients, conversion rates, composition of the particle, composition of copolymer on the morphology development No se puede mostrar la imagen. Puede que su equipo no tenga suficiente memoria para abrir la imagen o que ésta esté dañada. Reinicie el equipo y, a continuación, abra el archivo de nuevo. Si sigue apareciendo la x roja, puede que tenga que borrar la imagen e insertarla de nuevo. X=0; T=0 X=0.3; T=1800 X=0.5; T=3000 X=0.9; T=5400 X=1; T=6000 X=1; T=50000
20 Conclusions A novel modeling approach for prediction of dynamics of the particle morphology development in the composite waterborne systems including the systems with in situ formation of graft copolymer is presented The proposed model is based on stochastic dynamics (SD) accounts for the effect of phase compatibility and internal viscosity of the particles for the first time, is able to predict the morphologies of interesting new materials, such as polymer-polymer, polymer-polymer-inorganic complex hybrids The future work will include developing the parameter optimization technique; improving sampling and increasing accuracy through more advanced force fields and water representation
Mathematical Modeling of Chemical Reactions at Basque Center for Applied Mathematics. Simone Rusconi. December 11, 2015
Mathematical Modeling of Chemical Reactions at Basque Center for Applied Mathematics Simone Rusconi December 11, 2015 Basque Center for Applied Mathematics (BCAM) BCAM - Basque Center for Applied Mathematics,
More informationINTRODUCTION TO. Adapted from CS474/674 Prof. George Bebis Department of Computer Science & Engineering University of Nevada (UNR)
INTRODUCTION TO WAVELETS Adapted from CS474/674 Prof. George Bebis Department of Computer Science & Engineering University of Nevada (UNR) CRITICISM OF FOURIER SPECTRUM It gives us the spectrum of the
More informationTHERMODYNAMICS AND DYNAMICS OF CONFINED NANO-PHASES. Keith E. Gubbins
THERMODYNAMICS AND DYNAMICS OF CONFINED NANO-PHASES Keith E. Gubbins Institute for Computational Science & Engineering Department of Chemical and Biomolecular Engineering, North Carolina State University,
More informationNo se puede mostrar la imagen. Puede que su equipo no tenga suficiente memoria para abrir la imagen o que ésta esté dañada. Reinicie el equipo y, a
No se puede mostrar la imagen. Puede que su equipo no tenga suficiente memoria para abrir la imagen o que ésta esté dañada. Reinicie el equipo y, a continuación, abra el archivo de nuevo. Si sigue apareciendo
More informationThomas Nilsson. Status of FAIR-NuSTAR and R 3 B. MAGISOL-meeting, IEM-CSIC, Madrid, January
+ Status of FAIR-NuSTAR and R 3 B Thomas Nilsson MAGISOL-meeting, IEM-CSIC, Madrid, January 19 2009 Chalmers University of Technology/The Royal Swedish Academy of Sciences + FAIR Facility for Antiproton
More informationEnabling constant pressure hybrid Monte Carlo simulations using the GROMACS molecular simulation package
Enabling constant pressure hybrid Monte Carlo simulations using the GROMACS molecular simulation package Mario Fernández Pendás MSBMS Group Supervised by Bruno Escribano and Elena Akhmatskaya BCAM 18 October
More informationDynamics of the particle morphology during. the synthesis of waterborne polymer-inorganic. hybrids
Supporting Information for: Dynamics of the particle morphology during the synthesis of waterborne polymer-inorganic hybrids Shaghayegh Hamzehlou; Miren Aguirre; Jose R. Leiza; José M. Asua* POLYMAT, Kimika
More informationMix & Match Hamiltonian Monte Carlo
Mix & Match Hamiltonian Monte Carlo Elena Akhmatskaya,2 and Tijana Radivojević BCAM - Basque Center for Applied Mathematics, Bilbao, Spain 2 IKERBASQUE, Basque Foundation for Science, Bilbao, Spain MCMSki
More informationEngineering aspect of emulsion polymerization
Engineering aspect of emulsion polymerization Joong-In Kim Bayer Corp., Plastics, Technology Yonsei University Contents Free radical polymerization kinetics Emulsion polymerization Reactor configuration
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 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 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 informationADAPTIVE TWO-STAGE INTEGRATORS FOR SAMPLING ALGORITHMS BASED ON HAMILTONIAN DYNAMICS
ADAPTIVE TWO-STAGE INTEGRATORS FOR SAMPLING ALGORITHMS BASED ON HAMILTONIAN DYNAMICS E. Akhmatskaya a,c, M. Fernández-Pendás a, T. Radivojević a, J. M. Sanz-Serna b a Basque Center for Applied Mathematics
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 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 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 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 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 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 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 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 informationCosmic Ray (Stochastic) Acceleration from a. Background Plasma
Cosmic Ray (Stochastic) Acceleration rom a V. A. Dogiel Background Plasma P.N.Lebedev Institute o Physics,Moscow, Russia Thanks to my collaborators: D. O. Chernyshov (LPI), K. S. Cheng (HKU), A. D. Erlykin
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 informationAdvanced Molecular Dynamics
Advanced Molecular Dynamics Introduction May 2, 2017 Who am I? I am an associate professor at Theoretical Physics Topics I work on: Algorithms for (parallel) molecular simulations including GPU acceleration
More informationNew Hybrid Monte Carlo Methods for Efficient Sampling: from Physics to Biology and Statistics
TECHNICAL MATERIAL New Hybrid Monte Carlo Methods for Efficient Sampling: from Physics to Biology and Statistics Elena AKHMATSKAYA 1-3* and Sebastian REICH 4 1 Fujitsu Laboratories of Europe Ltd, Hayes
More informationpolymerization of n-butyl acrylate
SUPPORTING INFORMATION Intermolecular transfer to polymer in the radical polymerization of n-butyl acrylate Nicholas Ballard, 1 Shaghayegh Hamzehlou, 1 José M. Asua 1* 1 POLYMAT and Kimika Aplikatua Saila,
More informationPrinciples of Equilibrium Statistical Mechanics
Debashish Chowdhury, Dietrich Stauffer Principles of Equilibrium Statistical Mechanics WILEY-VCH Weinheim New York Chichester Brisbane Singapore Toronto Table of Contents Part I: THERMOSTATICS 1 1 BASIC
More informationMolecular dynamics simulation. CS/CME/BioE/Biophys/BMI 279 Oct. 5 and 10, 2017 Ron Dror
Molecular dynamics simulation CS/CME/BioE/Biophys/BMI 279 Oct. 5 and 10, 2017 Ron Dror 1 Outline Molecular dynamics (MD): The basic idea Equations of motion Key properties of MD simulations Sample applications
More informationNanoparticle Synthesis and Assembly from Atomistic Simulation Studies
Nanoparticle Synthesis and Assembly from Atomistic Simulation Studies By Takumi Hawa School of Aerospace & Mechanical Engineering Oklahoma Supercomputing Symposium 2008 October 7 Supported by NSF, ARL,
More informationGSHMC: An efficient Markov chain Monte Carlo sampling method. Sebastian Reich in collaboration with Elena Akhmatskaya (Fujitsu Laboratories Europe)
GSHMC: An efficient Markov chain Monte Carlo sampling method Sebastian Reich in collaboration with Elena Akhmatskaya (Fujitsu Laboratories Europe) 1. Motivation In the first lecture, we started from a
More informationCrystal Structure Prediction using CRYSTALG program
Crystal Structure Prediction using CRYSTALG program Yelena Arnautova Baker Laboratory of Chemistry and Chemical Biology, Cornell University Problem of crystal structure prediction: - theoretical importance
More informationStudy of Block Copolymer Lithography using SCFT: New Patterns and Methodology
Study of Block Copolymer Lithography using SCFT: New Patterns and Methodology Su-Mi Hur Glenn Fredrickson Complex Fluids Design Consortium Annual Meeting Monday, February 2, 2009 Materials Research Laboratory
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 informationAndré Schleife Department of Materials Science and Engineering
André Schleife Department of Materials Science and Engineering Length Scales (c) ICAMS: http://www.icams.de/cms/upload/01_home/01_research_at_icams/length_scales_1024x780.png Goals for today: Background
More informationBasics of Statistical Mechanics
Basics of Statistical Mechanics Review of ensembles Microcanonical, canonical, Maxwell-Boltzmann Constant pressure, temperature, volume, Thermodynamic limit Ergodicity (see online notes also) Reading assignment:
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 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 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 informationWhat is the role of simulation in nanoscience research?
ChE/MSE 557 Intro part 2 What is the role of simulation in nanoscience research? 1 Opportunities for Simulation Simulation Simulation complements both experiment and theory. Extends window of observation
More informationFile ISM02. Dynamics of Soft Matter
File ISM02 Dynamics of Soft Matter 1 Modes of dynamics Quantum Dynamics t: fs-ps, x: 0.1 nm (seldom of importance for soft matter) Molecular Dynamics t: ps µs, x: 1 10 nm Brownian Dynamics t: ns>ps, x:
More informationUnderstanding Molecular Simulation 2009 Monte Carlo and Molecular Dynamics in different ensembles. Srikanth Sastry
JNCASR August 20, 21 2009 Understanding Molecular Simulation 2009 Monte Carlo and Molecular Dynamics in different ensembles Srikanth Sastry Jawaharlal Nehru Centre for Advanced Scientific Research, Bangalore
More 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 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 informationDirected Assembly of Functionalized Nanoparticles with Amphiphilic Diblock Copolymers. Contents
Electronic Supplementary Material (ESI) for Physical Chemistry Chemical Physics. This journal is the Owner Societies 2017 Electronic Supplementary Information for Directed Assembly of Functionalized Nanoparticles
More informationMultiple time step Monte Carlo
JOURNAL OF CHEMICAL PHYSICS VOLUME 117, NUMBER 18 8 NOVEMBER 2002 Multiple time step Monte Carlo Balázs Hetényi a) Department of Chemistry, Princeton University, Princeton, NJ 08544 and Department of Chemistry
More informationMonte Carlo. Lecture 15 4/9/18. Harvard SEAS AP 275 Atomistic Modeling of Materials Boris Kozinsky
Monte Carlo Lecture 15 4/9/18 1 Sampling with dynamics In Molecular Dynamics we simulate evolution of a system over time according to Newton s equations, conserving energy Averages (thermodynamic properties)
More informationModeling Biomolecular Systems II. BME 540 David Sept
Modeling Biomolecular Systems II BME 540 David Sept Introduction Why do we perform simulations? What makes simulations possible? How do we perform simulations? What types of things/systems do we simulate?
More informationControlled healing of graphene nanopore
Controlled healing of graphene nanopore Konstantin Zakharchenko Alexander Balatsky Zakharchenko K.V., Balatsky A.V. Controlled healing of graphene nanopore. Carbon (80), December 2014, pp. 12 18. http://dx.doi.org/10.1016/j.carbon.2014.07.085
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 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 informationIntroduction to Granular Physics and Modeling Methods
Introduction to Granular Physics and Modeling Methods Stefan Luding MSM, TS, CTW, UTwente, NL Stefan Luding, s.luding@utwente.nl MSM, TS, CTW, UTwente, NL Granular Materials Real: sand, soil, rock, grain,
More informationRenormalization Group: non perturbative aspects and applications in statistical and solid state physics.
Renormalization Group: non perturbative aspects and applications in statistical and solid state physics. Bertrand Delamotte Saclay, march 3, 2009 Introduction Field theory: - infinitely many degrees of
More informationMATHEMATICAL SIMULATION OF THE STYRENE- BUTADIENE RUBBER S PRODUCTION IN THE CASCADE OF REACTORS BY THE MONTE CARLO METHOD
Int. J. Chem. Sci.: 4(4), 206, 865-876 ISSN 0972-768X www.sadgurupublications.com MATHEMATICAL SIMULATION OF THE STYRENE- BUTADIENE RUBBER S PRODUCTION IN THE CASCADE OF REACTORS BY THE MONTE CARLO METHOD
More information3.320: Lecture 19 (4/14/05) Free Energies and physical Coarse-graining. ,T) + < σ > dµ
3.320: Lecture 19 (4/14/05) F(µ,T) = F(µ ref,t) + < σ > dµ µ µ ref Free Energies and physical Coarse-graining T S(T) = S(T ref ) + T T ref C V T dt Non-Boltzmann sampling and Umbrella sampling Simple
More informationA Study of the Thermal Properties of a One. Dimensional Lennard-Jones System
A Study of the Thermal Properties of a One Dimensional Lennard-Jones System Abstract In this study, the behavior of a one dimensional (1D) Lennard-Jones (LJ) system is simulated. As part of this research,
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 informationDiffuse Interface Field Approach (DIFA) to Modeling and Simulation of Particle-based Materials Processes
Diffuse Interface Field Approach (DIFA) to Modeling and Simulation of Particle-based Materials Processes Yu U. Wang Department Michigan Technological University Motivation Extend phase field method to
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 informationPhysical Chemistry of Polymers (4)
Physical Chemistry of Polymers (4) Dr. Z. Maghsoud CONCENTRATED SOLUTIONS, PHASE SEPARATION BEHAVIOR, AND DIFFUSION A wide range of modern research as well as a variety of engineering applications exist
More informationIntroduction to Polymerization Processes
Introduction to Polymerization Processes Reference: Aspen Polymers: Unit Operations and Reaction Models, Aspen Technology, Inc., 2013. 1- Polymer Definition A polymer is a macromolecule made up of many
More informationMonte Carlo (MC) Simulation Methods. Elisa Fadda
Monte Carlo (MC) Simulation Methods Elisa Fadda 1011-CH328, Molecular Modelling & Drug Design 2011 Experimental Observables A system observable is a property of the system state. The system state i is
More 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 informationCopyright 2001 University of Cambridge. Not to be quoted or copied without permission.
Course MP3 Lecture 4 13/11/2006 Monte Carlo method I An introduction to the use of the Monte Carlo method in materials modelling Dr James Elliott 4.1 Why Monte Carlo? The name derives from the association
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 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 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 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 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 informationA Nobel Prize for Molecular Dynamics and QM/MM What is Classical Molecular Dynamics? Simulation of explicit particles (atoms, ions,... ) Particles interact via relatively simple analytical potential
More information3.320 Lecture 18 (4/12/05)
3.320 Lecture 18 (4/12/05) Monte Carlo Simulation II and free energies Figure by MIT OCW. General Statistical Mechanics References D. Chandler, Introduction to Modern Statistical Mechanics D.A. McQuarrie,
More informationEarly stages of dewetting of microscopically thin polymer films: A molecular dynamics study
JOURNAL OF CHEMICAL PHYSICS VOLUME 109, NUMBER 19 15 NOVEMBER 1998 Early stages of dewetting of microscopically thin polymer films: A molecular dynamics study Hong Liu Department of Physics, Kansas State
More informationBioengineering 215. An Introduction to Molecular Dynamics for Biomolecules
Bioengineering 215 An Introduction to Molecular Dynamics for Biomolecules David Parker May 18, 2007 ntroduction A principal tool to study biological molecules is molecular dynamics simulations (MD). MD
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 informationFor info and ordering all the 4 versions / languages of this book please visit: http://trl.lab.uic.edu/pon Contents Preface vii Chapter 1 Advances in Atomic and Molecular Nanotechnology Introduction 1
More informationAdvanced sampling. fluids of strongly orientation-dependent interactions (e.g., dipoles, hydrogen bonds)
Advanced sampling ChE210D Today's lecture: methods for facilitating equilibration and sampling in complex, frustrated, or slow-evolving systems Difficult-to-simulate systems Practically speaking, one is
More informationMOLECULAR DYNAMIC SIMULATION OF WATER VAPOR INTERACTION WITH VARIOUS TYPES OF PORES USING HYBRID COMPUTING STRUCTURES
MOLECULAR DYNAMIC SIMULATION OF WATER VAPOR INTERACTION WITH VARIOUS TYPES OF PORES USING HYBRID COMPUTING STRUCTURES V.V. Korenkov 1,3, a, E.G. Nikonov 1, b, M. Popovičová 2, с 1 Joint Institute for Nuclear
More informationHybrid Latex Systems- PU, Alkyd and Inorganic/Acrylic Latices
Emulsion Polymers Consulting and Education, LLC present: Hybrid Latex Systems- PU, Alkyd and Inorganic/Acrylic Latices WORKSHOP OBJECTIVES: This workshop is designed to provide industrial scientists and
More informationExploring the energy landscape
Exploring the energy landscape ChE210D Today's lecture: what are general features of the potential energy surface and how can we locate and characterize minima on it Derivatives of the potential energy
More informationPredici 11 Quick Overview
Predici 11 Quick Overview PREDICI is the leading simulation package for kinetic, process and property modeling with a major emphasis on macromolecular systems. It has been successfully utilized to model
More informationOverview of Accelerated Simulation Methods for Plasma Kinetics
Overview of Accelerated Simulation Methods for Plasma Kinetics R.E. Caflisch 1 In collaboration with: J.L. Cambier 2, B.I. Cohen 3, A.M. Dimits 3, L.F. Ricketson 1,4, M.S. Rosin 1,5, B. Yann 1 1 UCLA Math
More informationHybrid Latex Systems- PU, Alkyd and Inorganic/Acrylic Latices
Emulsion Polymers Consulting and Education, LLC presents Hybrid Latex Systems- PU, Alkyd and Inorganic/Acrylic Latices WORKSHOP OBJECTIVES: This workshop is designed to provide industrial scientists and
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 informationRe-entrant transition in a one-dimensional growth model
Re-entrant transition in a one-dimensional growth model Ostap Hryniv Dept. of Mathematical Sciences Durham University Ostap.Hryniv@durham.ac.uk X International Conference Stochastic and Analytic Methods
More informationAn introduction to Molecular Dynamics. EMBO, June 2016
An introduction to Molecular Dynamics EMBO, June 2016 What is MD? everything that living things do can be understood in terms of the jiggling and wiggling of atoms. The Feynman Lectures in Physics vol.
More informationResearch Statement. Shenggao Zhou. November 3, 2014
Shenggao Zhou November 3, My research focuses on: () Scientific computing and numerical analysis (numerical PDEs, numerical optimization, computational fluid dynamics, and level-set method for interface
More informationHybrid Latex Systems- PU, Alkyd and Inorganic/Acrylic Latices
Emulsion Polymers Consulting and Education, LLC present: Hybrid Latex Systems- PU, Alkyd and Inorganic/Acrylic Latices WORKSHOP OBJECTIVES: This workshop is designed to provide industrial scientists and
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 informationElectrostatic Self-assembly : A New Route Towards Nanostructures
1 Electrostatic Self-assembly : A New Route Towards Nanostructures J.-F. Berret, P. Hervé, M. Morvan Complex Fluids Laboratory, UMR CNRS - Rhodia n 166, Cranbury Research Center Rhodia 259 Prospect Plains
More informationMonte Carlo Study of the Crystalline and Amorphous NaK Alloy
Available online at www.sciencedirect.com ScienceDirect Procedia Computer Science 08C (207) 25 22 International Conference on Computational Science, ICCS 207, 2-4 June 207, Zurich, Switzerland Monte Carlo
More informationBrief introduction to Markov Chain Monte Carlo
Brief introduction to Department of Probability and Mathematical Statistics seminar Stochastic modeling in economics and finance November 7, 2011 Brief introduction to Content 1 and motivation Classical
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 informationViscous Flow Computations with Molecular Dynamics
FAKULTÄT FÜR INFORMATIK DER TECHNISCHEN UNIVERSITÄT MÜNCHEN Bachelorarbeit in Informatik Viscous Flow Computations with Molecular Dynamics Ching-Yu Kao FAKULTÄT FÜR INFORMATIK DER TECHNISCHEN UNIVERSITÄT
More informationIntroduction to Simulation - Lectures 17, 18. Molecular Dynamics. Nicolas Hadjiconstantinou
Introduction to Simulation - Lectures 17, 18 Molecular Dynamics Nicolas Hadjiconstantinou Molecular Dynamics Molecular dynamics is a technique for computing the equilibrium and non-equilibrium properties
More informationEncapsulation. Battelle Technology. Introduction
Encapsulation Introduction The encapsulation methods reported in the literature 1-7 for the production of microcapsules are generally achieved using one of the following techniques: 1. Phase separation
More informationPattern Recognition and Machine Learning. Bishop Chapter 11: Sampling Methods
Pattern Recognition and Machine Learning Chapter 11: Sampling Methods Elise Arnaud Jakob Verbeek May 22, 2008 Outline of the chapter 11.1 Basic Sampling Algorithms 11.2 Markov Chain Monte Carlo 11.3 Gibbs
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 informationJavier Junquera. Statistical mechanics
Javier Junquera Statistical mechanics From the microscopic to the macroscopic level: the realm of statistical mechanics Computer simulations Thermodynamic state Generates information at the microscopic
More informationCHAPTER V. Brownian motion. V.1 Langevin dynamics
CHAPTER V Brownian motion In this chapter, we study the very general paradigm provided by Brownian motion. Originally, this motion is that a heavy particle, called Brownian particle, immersed in a fluid
More information