Hydrogels in charged solvents

Size: px
Start display at page:

Download "Hydrogels in charged solvents"

Transcription

1 Hydrogels in charged solvents Peter Košovan, Christian Holm, Tobias Richter 27. May 2013 Institute for Computational Physics Pfaffenwaldring 27 D Stuttgart Germany Tobias Richter (ICP, Stuttgart) Hydrogels 27. May / 21

2 Outline Outline 1 Introduction 2 Motivation 3 Simulation Setup 4 Results 5 Conclusions 6 References Tobias Richter (ICP, Stuttgart) Hydrogels 27. May / 21

3 Introduction Hydrogels are charged polymer-networks Hydrogels can swell tremendously in aqueous solutions Swelling behavior can be influenced by external parameters Pressure ph Concentration of electrolytes Tobias Richter (ICP, Stuttgart) Hydrogels 27. May / 21

4 Motivation Experiment [1] Johannes Höpfner at KIT was investigating the desalination capacity of Hydrogels Poly-acrylic-acid Idea: Swelling Hydrogel in salty water solution extract Hydrogel from salt solution Deswell/Compress Hydrogel and squeeze out the incorporated water Figure 1: Experimental setup, taken from [1] Tobias Richter (ICP, Stuttgart) Hydrogels 27. May / 21

5 Motivation Experimental Results The process is reversible Desalination capacity per cycle is up to 35% [1] Tobias Richter (ICP, Stuttgart) Hydrogels 27. May / 21

6 Motivation Experimental Results The process is reversible Desalination capacity per cycle is up to 35% [1] Figure 2: Salt concentration of the extracted water, taken from [1]. Tobias Richter (ICP, Stuttgart) Hydrogels 27. May / 21

7 Motivation Experimental Results The process is reversible Desalination capacity per cycle is up to 35% [1] Open Questions How does the desalination work on the microscopic level What influences the desalination capacity Crosslinking density Charge density Figure 2: Salt concentration of the extracted water, taken from [1]. Tobias Richter (ICP, Stuttgart) Hydrogels 27. May / 21

8 Simulation Setup Simulation setup Tetrafunctional polymer-network with 16 chains and 8 nodes and equally spaced charged on the chains Periodic boundary conditions Explicit counterions Langevin thermostat Brownian dynamics Simulation in (NVT) and (µvt) ensemble as well as testing those together with the npt-isotropic barostat in ESPResSo Tobias Richter (ICP, Stuttgart) Hydrogels 27. May / 21

9 Simulation Setup Simulation setup Tetrafunctional polymer-network with 16 chains and 8 nodes and equally spaced charged on the chains Periodic boundary conditions Explicit counterions Langevin thermostat Brownian dynamics Simulation in (NVT) and (µvt) ensemble as well as testing those together with the npt-isotropic barostat in ESPResSo Translation of experiment to computer simulation 1 Determine chemical potential and pressure in salt solution 2 Determine equilibrium swelling and internal salt concentration at a certain external salt concentration c s fixed external chemical potential equilibrium conditions p = p salt (c s ) and µ i = µ i Tobias Richter (ICP, Stuttgart) Hydrogels 27. May / 21

10 Simulation Setup Figure 3: Snapshot of the simulated system Tobias Richter (ICP, Stuttgart) Hydrogels 27. May / 21

11 Simulation Setup Simulation in the grand canonical ensemble Using the (NVT) ensemble simulation ESPResSo already provides Adding particle exchange with external reservoir Monte Carlo steps Acceptance probability according to [2] ( ) V acc(n N + 1) = min 1, Λ 3 (N + 1) exp ( β ( µ out + U)) (1) ( ) acc(n N 1) = min 1, Λ3 N V exp ( β (µ out + U)) (2) Tobias Richter (ICP, Stuttgart) Hydrogels 27. May / 21

12 Simulation Setup Simulation in the grand canonical ensemble ( ) V acc(n N + 1) = min 1, Λ 3 (N + 1) exp ( β ( µ out + U)) ( ) acc(n N 1) = min 1, Λ3 N V exp ( β (µ out + U)) (1) (2) Elimination of Λ with µ out = µ ex out + k B T ln(ϱ out Λ 3 ) Tobias Richter (ICP, Stuttgart) Hydrogels 27. May / 21

13 Simulation Setup Simulation in the grand canonical ensemble ( ) V acc(n N + 1) = min 1, Λ 3 (N + 1) exp ( β ( µ out + U)) ( ) acc(n N 1) = min 1, Λ3 N V exp ( β (µ out + U)) (1) (2) Rewritten acceptance rates ( ) 1 acc(n N + 1) = min 1, exp ( β ( µ ex out 1/β ln ϱ out + U)) ϱ in acc(n N 1) = min (1, ϱ in exp ( β (µ ex out + 1/β ln ϱ out + U))) (4) (3) Tobias Richter (ICP, Stuttgart) Hydrogels 27. May / 21

14 Simulation Setup µ ex [k B T] Widom method extended Debye Hückel Experiment c s [mol/l] Figure 4: Chemical potential obtained with for Bjerrum-length equal 2 together with a fit of the extended Debye-Hückel model and experimental activity coefficients measurements done by Truesdell [3]. Tobias Richter (ICP, Stuttgart) Hydrogels 27. May / 21

15 Simulation Setup µ ex [k B T] Widom method extended Debye Hückel Experiment c s [mol/l] Figure 4: Chemical potential obtained with for Bjerrum-length equal 2 together with a fit of the extended Debye-Hückel model and experimental activity coefficients measurements done by Truesdell [3]. 3 2πλ 2 b ϱ µ DB = 1 + 8πa ϱ Tobias Richter (ICP, Stuttgart) Hydrogels 27. May / 21 (5)

16 Simulation Setup Check of the Grand canonical Algorithm: Running the simulation with an empty box, no hydrogel and measure the salt density outside density: 10 3 µ ex obtained from Widom particle insertion Density a Simulation time Figure 5: Empty box test. Tobias Richter (ICP, Stuttgart) Hydrogels 27. May / 21

17 Simulation Setup Check of the Grand canonical Algorithm: Running the simulation with an empty box, no hydrogel and measure the salt density outside density: 2.7e 3 µ ex obtained from Widom particle insertion c in 3.2e e e e e e e e e e e e e t 79 1/4 Figure 5: Development of the salt concentration inside a hydrogel (N=79 f=1/4 c s = 0.1mol/L ˆ=2.7e 3) Tobias Richter (ICP, Stuttgart) Hydrogels 27. May / 21

18 Simulation Setup Determining the pressure dependence on Volume for a given external salt solution N=59 f=1/4 external pressure Fit Internal pressure Box length [σ] Figure 6: Total pressure dependence on volume for a salt solution with concentration 0.2mol/L. Tobias Richter (ICP, Stuttgart) Hydrogels 27. May / 21

19 Simulation Setup Donnan Theory Philipse et. al derived an equation for colloid system in contact with an infinite reservoir of salt [4]. ϱ = ϱ ( ) 2 c ϱc + 1 +, (6) ϱ ext 2ϱ ext 2ϱ ext ϱ is the salt concentration in the hydrogel ϱ c = mnf V hg, m is the number of polymer chains and N the number of monomers per chain, f is the charge fraction of the hydrogel Tobias Richter (ICP, Stuttgart) Hydrogels 27. May / 21

20 Simulation Setup Osmotic Donnan Model Combining Donnan theory for the salt partitioning with an equation for the pressure balance Equation of state Tobias Richter (ICP, Stuttgart) Hydrogels 27. May / 21

21 Simulation Setup Osmotic Donnan Model Combining Donnan theory for the salt partitioning with an equation for the pressure balance Equation of state Minimal model: the chains in the hydrogel behave as ideal gaussian chains p osm in p in = p out + p elastic + p extern = p osm salt Tobias Richter (ICP, Stuttgart) Hydrogels 27. May / 21

22 Simulation Setup Osmotic Donnan Model Combining Donnan theory for the salt partitioning with an equation for the pressure balance Equation of state Minimal model: the chains in the hydrogel behave as ideal gaussian chains (i) : mnf V hg + ϱ = ϱ + (7) (ii) : ϱ 2 ext = ϱ ϱ + (8) ( (iii) : ϱ + + ϱ 1 ) Nb 2 k B T + p ext = 2ϱ ext k B T (9) R e Tobias Richter (ICP, Stuttgart) Hydrogels 27. May / 21

23 Results Equilibrium swelling Swelling degree N=39 f=1/8 theory N=39 f=1/8 N=39 f=1/4 theory N=39 f=1/4 N=39 f=1/2 theory N=39 f=1/ c s [mol/l] Swelling degree N=59 f=1/4 theory N=59 f=1/4 N=59 f=1/2 theory N=59 f=1/ c s [mol/l] Figure 7: Dependence of the equilibrium swelling degree on the salt concentration of the external reservoir. The dashed lines are the predictions of the osmotic Donnan model. Tobias Richter (ICP, Stuttgart) Hydrogels 27. May / 21

24 Results Donnan theory c int [mol/l] Donnan 0.01 mol/l N=39 N=59 N=79 f=1/8 f=1/4 f=1/ f c p [mol/l] c int [mol/l] Donnan 0.1 mol/l N=39 N=59 N=79 f=1/8 f=1/4 f=1/ f c p [mol/l] Figure 8: Testing the Donnan prediction for the salt partitioning against the simulation for two intermediately high salt concentrations: left: 0.01 mol/l, right: 0.1 mol/l. ϱ ϱ ext = ϱ c 2ϱ ext ( ϱc 2ϱ ext ) 2, (10) Tobias Richter (ICP, Stuttgart) Hydrogels 27. May / 21

25 Results Donnan 0.37 mol/l N=59 f=1/ Donnan 0.1 mol/l N=59 f=1/ c int [mol/l] c int [mol/l] f c p [mol/l] f c p [mol/l] Figure 9: Testing the Donnan prediction for salt partitioning for c s = 0.37mol/L (left) and 0.2 mol/l with randomly placed charges on the hydrogel with N=59 f=1/4. Tobias Richter (ICP, Stuttgart) Hydrogels 27. May / 21

26 Results Donnan theory in the limit of no electrostatics (λ B = 0) c int [mol/l] Donnan 0.2 mol/l N=39 N=59 N=79 f=1/8 f=1/4 f=1/ f c p [mol/l] Figure 10: Comparison of the partitioning prediction of Donnan and simulation data for c = 0.2mol/L. Tobias Richter (ICP, Stuttgart) Hydrogels 27. May / 21

27 Results Simulation in µpt Figure 11: Pressure versus volume for the N=59 f=1/4 hydrogel. The circles indicate the start configurations for testing simulations in (µpt ). Tobias Richter (ICP, Stuttgart) Hydrogels 27. May / 21

28 Results Simulation in µpt external pressure Above Below pressure c in Box length [σ] t Figure 11: Left: The final volumes that were obtained using the (µpt ) ensemble simulation for the two different start configurations. Right: The salt density inside the hydrogels for the simulations in (µpt ). Tobias Richter (ICP, Stuttgart) Hydrogels 27. May / 21

29 Conclusions Conclusions The equilibrium swelling for hydrogels in various external salt concentrations could be determined The Osmotic Donnan model and the simulations qualitatively agree The Donnan theory for the salt partitioning show increasing deviations for increasing charge fractions of the hydrogel The linear charge density on the hydrogel chains breaks the assumption of homogeneity This seems to remain true even for very high salt concentrations (strong screening of the electrostatics) and randomly placed charges on the hydrogel backbone In the case of no electrostatics (the salt behaves almost as ideal gas) the donnan theory and the simulations fully agree The barostat in combination with the grand canonical ensemble does not work to predict the equilibrium swelling volume Tobias Richter (ICP, Stuttgart) Hydrogels 27. May / 21

30 Conclusions Outlook Improvement of the Donnan model Poisson-Boltzmann Cell model to treat the elastic properties in an electrolyte solution Track down the origin of the deviations Having a closer look at the distribution of ions around the chains in the hydrogel Tobias Richter (ICP, Stuttgart) Hydrogels 27. May / 21

31 References [1] Johannes Höpfner, Christopher Klein, and Manfred Wilhelm. A novel approach for the desalination of seawater by means of reusable poly(acrylic acid) hydrogels and mechanical force. Macromol. Rapid Commun., 31:1337, [2] Daan Frenkel and Berend Smit. Understanding Molecular Simulation. Academic Press, San Diego, second edition, [3] Alfred H. Truesdell. Activity coefficients of aqueous sodium chloride from 15 to 50 c measured with a glass electrode. Science, 161(3844): , [4] A Philipse and A Vrij. The donnan equilibrium: I. on the thermodynamic foundation of the donnan equation of state. Journal of Physics: Condensed Matter, 23(19):194106, Tobias Richter (ICP, Stuttgart) Hydrogels 27. May / 21

Understanding Molecular Simulation 2009 Monte Carlo and Molecular Dynamics in different ensembles. Srikanth Sastry

Understanding 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 information

Monte Carlo Methods. Ensembles (Chapter 5) Biased Sampling (Chapter 14) Practical Aspects

Monte Carlo Methods. Ensembles (Chapter 5) Biased Sampling (Chapter 14) Practical Aspects Monte Carlo Methods Ensembles (Chapter 5) Biased Sampling (Chapter 14) Practical Aspects Lecture 1 2 Lecture 1&2 3 Lecture 1&3 4 Different Ensembles Ensemble ame Constant (Imposed) VT Canonical,V,T P PT

More information

Gas-liquid phase separation in oppositely charged colloids: stability and interfacial tension

Gas-liquid phase separation in oppositely charged colloids: stability and interfacial tension 7 Gas-liquid phase separation in oppositely charged colloids: stability and interfacial tension We study the phase behaviour and the interfacial tension of the screened Coulomb (Yukawa) restricted primitive

More information

Introduction Statistical Thermodynamics. Monday, January 6, 14

Introduction 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 information

MONTE CARLO METHOD. Reference1: Smit Frenkel, Understanding molecular simulation, second edition, Academic press, 2002.

MONTE CARLO METHOD. Reference1: Smit Frenkel, Understanding molecular simulation, second edition, Academic press, 2002. MONTE CARLO METHOD Reference1: Smit Frenkel, Understanding molecular simulation, second edition, Academic press, 2002. Reference 2: David P. Landau., Kurt Binder., A Guide to Monte Carlo Simulations in

More information

Introduction 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 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 information

An electrokinetic LB based model for ion transport and macromolecular electrophoresis

An electrokinetic LB based model for ion transport and macromolecular electrophoresis An electrokinetic LB based model for ion transport and macromolecular electrophoresis Raffael Pecoroni Supervisor: Michael Kuron July 8, 2016 1 Introduction & Motivation So far an mesoscopic coarse-grained

More information

Phase 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 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 information

Molecular Driving Forces

Molecular Driving Forces Molecular Driving Forces Statistical Thermodynamics in Chemistry and Biology SUBGfittingen 7 At 216 513 073 / / Ken A. Dill Sarina Bromberg With the assistance of Dirk Stigter on the Electrostatics chapters

More information

Introduction to molecular dynamics

Introduction 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 information

Colloidal Suspension Rheology Chapter 1 Study Questions

Colloidal 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 information

2. Thermodynamics. Introduction. Understanding Molecular Simulation

2. Thermodynamics. Introduction. Understanding Molecular Simulation 2. Thermodynamics Introduction Molecular Simulations Molecular dynamics: solve equations of motion r 1 r 2 r n Monte Carlo: importance sampling r 1 r 2 r n How do we know our simulation is correct? Molecular

More information

Contents. 1 Introduction and guide for this text 1. 2 Equilibrium and entropy 6. 3 Energy and how the microscopic world works 21

Contents. 1 Introduction and guide for this text 1. 2 Equilibrium and entropy 6. 3 Energy and how the microscopic world works 21 Preface Reference tables Table A Counting and combinatorics formulae Table B Useful integrals, expansions, and approximations Table C Extensive thermodynamic potentials Table D Intensive per-particle thermodynamic

More information

Molecular Simulation Background

Molecular Simulation Background Molecular Simulation Background Why Simulation? 1. Predicting properties of (new) materials 2. Understanding phenomena on a molecular scale 3. Simulating known phenomena? Example: computing the melting

More information

8.592J HST.452J: Statistical Physics in Biology

8.592J HST.452J: Statistical Physics in Biology Assignment # 4 8.592J HST.452J: Statistical Physics in Biology Coulomb Interactions 1. Flory Theory: The Coulomb energy of a ball of charge Q and dimension R in d spacial dimensions scales as Q 2 E c.

More information

Multiscale Materials Modeling

Multiscale Materials Modeling Multiscale Materials Modeling Lecture 02 Capabilities of Classical Molecular Simulation These notes created by David Keffer, University of Tennessee, Knoxville, 2009. Outline Capabilities of Classical

More information

V = 2ze 2 n. . a. i=1

V = 2ze 2 n. . a. i=1 IITS: Statistical Physics in Biology Assignment # 3 KU Leuven 5/29/2013 Coulomb Interactions & Polymers 1. Flory Theory: The Coulomb energy of a ball of charge Q and dimension R in d spacial dimensions

More information

Hyeyoung Shin a, Tod A. Pascal ab, William A. Goddard III abc*, and Hyungjun Kim a* Korea

Hyeyoung Shin a, Tod A. Pascal ab, William A. Goddard III abc*, and Hyungjun Kim a* Korea The Scaled Effective Solvent Method for Predicting the Equilibrium Ensemble of Structures with Analysis of Thermodynamic Properties of Amorphous Polyethylene Glycol-Water Mixtures Hyeyoung Shin a, Tod

More information

Statistical Mechanics

Statistical Mechanics Franz Schwabl Statistical Mechanics Translated by William Brewer Second Edition With 202 Figures, 26 Tables, and 195 Problems 4u Springer Table of Contents 1. Basic Principles 1 1.1 Introduction 1 1.2

More information

Benchmarking the Hartree-Fock and Hartree-Fock-Bogoliubov approximations to level densities. G.F. Bertsch, Y. Alhassid, C.N. Gilbreth, and H.

Benchmarking the Hartree-Fock and Hartree-Fock-Bogoliubov approximations to level densities. G.F. Bertsch, Y. Alhassid, C.N. Gilbreth, and H. Benchmarking the Hartree-Fock and Hartree-Fock-Bogoliubov approximations to level densities G.F. Bertsch, Y. Alhassid, C.N. Gilbreth, and H. Nakada 5th Workshop on Nuclear Level Density and Gamma Strength,

More information

Supplemental Material for Temperature-sensitive colloidal phase behavior induced by critical Casimir forces

Supplemental Material for Temperature-sensitive colloidal phase behavior induced by critical Casimir forces Supplemental Material for Temperature-sensitive colloidal phase behavior induced by critical Casimir forces Minh Triet Dang, 1 Ana Vila Verde, 2 Van Duc Nguyen, 1 Peter G. Bolhuis, 3 and Peter Schall 1

More information

INTRODUCTION TO о JLXJLA Из А lv-/xvj_y JrJrl Y üv_>l3 Second Edition

INTRODUCTION TO о JLXJLA Из А lv-/xvj_y JrJrl Y üv_>l3 Second Edition INTRODUCTION TO о JLXJLA Из А lv-/xvj_y JrJrl Y üv_>l3 Second Edition Kerson Huang CRC Press Taylor & Francis Croup Boca Raton London New York CRC Press is an imprint of the Taylor & Francis Group an Informa

More information

Proteins 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 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 information

Basics of Statistical Mechanics

Basics 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 information

Monte Carlo Methods in Statistical Mechanics

Monte Carlo Methods in Statistical Mechanics Monte Carlo Methods in Statistical Mechanics Mario G. Del Pópolo Atomistic Simulation Centre School of Mathematics and Physics Queen s University Belfast Belfast Mario G. Del Pópolo Statistical Mechanics

More information

INTERMOLECULAR AND SURFACE FORCES

INTERMOLECULAR 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 information

fiziks Institute for NET/JRF, GATE, IIT-JAM, JEST, TIFR and GRE in PHYSICAL SCIENCES

fiziks Institute for NET/JRF, GATE, IIT-JAM, JEST, TIFR and GRE in PHYSICAL SCIENCES Content-Thermodynamics & Statistical Mechanics 1. Kinetic theory of gases..(1-13) 1.1 Basic assumption of kinetic theory 1.1.1 Pressure exerted by a gas 1.2 Gas Law for Ideal gases: 1.2.1 Boyle s Law 1.2.2

More information

Non-spherical Wigner-Seitz cells in linear Poisson-Boltzmann theory

Non-spherical Wigner-Seitz cells in linear Poisson-Boltzmann theory Universiteit Utrecht Department of Physics and Astronomy Institute for Theoretical Physics Master of Science Thesis Non-spherical Wigner-Seitz cells in linear Poisson-Boltzmann theory by Kasimir van Rijn

More information

The Effect of Electrostatic Surface Charges on Photoresist Dissolution

The Effect of Electrostatic Surface Charges on Photoresist Dissolution Sean Burns 4/24/00 Dr. Bonnecaze ChE 385M The Effect of Electrostatic Surface Charges on Photoresist Dissolution Introduction/Motivation It would be very useful and economical to have a fundamental model

More information

(# = %(& )(* +,(- Closed system, well-defined energy (or e.g. E± E/2): Microcanonical ensemble

(# = %(& )(* +,(- Closed system, well-defined energy (or e.g. E± E/2): Microcanonical ensemble Recall from before: Internal energy (or Entropy): &, *, - (# = %(& )(* +,(- Closed system, well-defined energy (or e.g. E± E/2): Microcanonical ensemble & = /01Ω maximized Ω: fundamental statistical quantity

More information

EFFECTS OF ADDED ELECTROLYTES ON THE STRUCTURE OF CHARGED POLYMERIC MICELLES

EFFECTS OF ADDED ELECTROLYTES ON THE STRUCTURE OF CHARGED POLYMERIC MICELLES Soft Materials, 3(2&3): 89 120, (2006) Copyright # Taylor & Francis Group, LLC ISSN 1539-445X print/1539-4468 online DOI: 10.1080/15394450600801228 EFFECTS OF DDED ELECTROLYTES ON THE STRUCTURE OF CHRGED

More information

3.320 Lecture 18 (4/12/05)

3.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 information

Lecture 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 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 information

arxiv: v1 [cond-mat.soft] 17 Sep 2016

arxiv: v1 [cond-mat.soft] 17 Sep 2016 Polyelectrolyte Polypeptide Scaling Laws Via Mechanical and Dielectric Relaxation Measurements Jorge Monreal University of South Florida arxiv:1609.05358v1 [cond-mat.soft] 17 Sep 2016 Dated: March 1, 2018

More information

Molecular Dynamic Simulation Study of the Volume Transition of PNIPAAm Hydrogels

Molecular Dynamic Simulation Study of the Volume Transition of PNIPAAm Hydrogels Molecular Dynamic Simulation Study of the Volume Transition of PNIPAAm Hydrogels Jonathan Walter 1, Jadran Vrabec 2, Hans Hasse 1 1 Laboratory of Engineering, University of Kaiserslautern, Germany 2 and

More information

Author's personal copy

Author's personal copy Adsorption (2011) 17: 265 271 DOI 10.1007/s10450-011-9325-7 Monte Carlo simulation of polymer adsorption Christopher J. Rasmussen Aleksey Vishnyakov Alexander V. Neimark Received: 15 May 2010 / Accepted:

More information

PHASE TRANSITIONS IN SOFT MATTER SYSTEMS

PHASE 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 information

(name) Electrochemical Energy Systems, Spring 2014, M. Z. Bazant. Final Exam

(name) Electrochemical Energy Systems, Spring 2014, M. Z. Bazant. Final Exam 10.626 Electrochemical Energy Systems, Spring 2014, M. Z. Bazant Final Exam Instructions. This is a three-hour closed book exam. You are allowed to have five doublesided pages of personal notes during

More information

424 Index. Eigenvalue in quantum mechanics, 174 eigenvector in quantum mechanics, 174 Einstein equation, 334, 342, 393

424 Index. Eigenvalue in quantum mechanics, 174 eigenvector in quantum mechanics, 174 Einstein equation, 334, 342, 393 Index After-effect function, 368, 369 anthropic principle, 232 assumptions nature of, 242 autocorrelation function, 292 average, 18 definition of, 17 ensemble, see ensemble average ideal,23 operational,

More information

Electrolyte Solutions

Electrolyte Solutions Chapter 8 Electrolyte Solutions In the last few chapters of this book, we will deal with several specific types of chemical systems. The first one is solutions and equilibria involving electrolytes, which

More information

Peter Debye and Electrochemistry

Peter Debye and Electrochemistry Peter Debye and Electrochemistry A K Shukla and T Prem Kumar A K Shukla is Professor and Chairman, Solid State and Structural Chemistry Unit, Indian Institute of Science, Bangalore. His research interests

More information

AR-7781 (Physical Chemistry)

AR-7781 (Physical Chemistry) Model Answer: B.Sc-VI th Semester-CBT-602 AR-7781 (Physical Chemistry) One Mark Questions: 1. Write a nuclear reaction for following Bethe s notation? 35 Cl(n, p) 35 S Answer: 35 17Cl + 1 1H + 35 16S 2.

More information

Principles of Equilibrium Statistical Mechanics

Principles 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 information

Metropolis Monte Carlo simulation of the Ising Model

Metropolis Monte Carlo simulation of the Ising Model Metropolis Monte Carlo simulation of the Ising Model Krishna Shrinivas (CH10B026) Swaroop Ramaswamy (CH10B068) May 10, 2013 Modelling and Simulation of Particulate Processes (CH5012) Introduction The Ising

More information

2 Structure. 2.1 Coulomb interactions

2 Structure. 2.1 Coulomb interactions 2 Structure 2.1 Coulomb interactions While the information needed for reproduction of living systems is chiefly maintained in the sequence of macromolecules, any practical use of this information must

More information

Elementary Lectures in Statistical Mechanics

Elementary 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 information

On the Dynamics and Disentanglement in Thin and Two-Dimensional Polymer Films

On 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

Effect of Polyelectrolyte Adsorption on Intercolloidal Forces

Effect of Polyelectrolyte Adsorption on Intercolloidal Forces 5042 J. Phys. Chem. B 1999, 103, 5042-5057 Effect of Polyelectrolyte Adsorption on Intercolloidal Forces Itamar Borukhov, David Andelman,*, and Henri Orland School of Physics and Astronomy, Raymond and

More information

Electrical double layer

Electrical double layer Electrical double layer Márta Berka és István Bányai, University of Debrecen Dept of Colloid and Environmental Chemistry http://dragon.unideb.hu/~kolloid/ 7. lecture Adsorption of strong electrolytes from

More information

*blood and bones contain colloids. *milk is a good example of a colloidal dispersion.

*blood and bones contain colloids. *milk is a good example of a colloidal dispersion. Chap. 3. Colloids 3.1. Introduction - Simple definition of a colloid: a macroscopically heterogeneous system where one component has dimensions in between molecules and macroscopic particles like sand

More information

Table of Contents [ttc]

Table of Contents [ttc] Table of Contents [ttc] 1. Equilibrium Thermodynamics I: Introduction Thermodynamics overview. [tln2] Preliminary list of state variables. [tln1] Physical constants. [tsl47] Equations of state. [tln78]

More information

Bchem 675 Lecture 9 Electrostatics-Lecture 2 Debye-Hückel: Continued Counter ion condensation

Bchem 675 Lecture 9 Electrostatics-Lecture 2 Debye-Hückel: Continued Counter ion condensation Bchem 675 Lecture 9 Electrostatics-Lecture 2 Debye-Hückel: Continued Counter ion condensation Ion:ion interactions What is the free energy of ion:ion interactions ΔG i-i? Consider an ion in a solution

More information

Computer simulation methods (1) Dr. Vania Calandrini

Computer 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 information

Polymer Gels. Boulder Lectures in Soft Matter Physics July 2012 Yitzhak Rabin

Polymer Gels. Boulder Lectures in Soft Matter Physics July 2012 Yitzhak Rabin Polymer Gels Boulder ectures in Soft Matter Physics July Yitzhak abin M. ubinstein and.h. Colby, Polymer Physics (Oxford, ), Chapters 6 and 7 P.-G. de Gennes, Scaling Concepts in Polymer Physics (Cornell,

More information

A Brief Introduction to Statistical Mechanics

A Brief Introduction to Statistical Mechanics A Brief Introduction to Statistical Mechanics E. J. Maginn, J. K. Shah Department of Chemical and Biomolecular Engineering University of Notre Dame Notre Dame, IN 46556 USA Monte Carlo Workshop Universidade

More information

Activities and Activity Coefficients

Activities and Activity Coefficients CHEM 331 Physical Chemistry Fall 017 Activities and Activity Coefficients We now finish answering the question we asked during our last lecture, what is the form of the chemical potential i (T,P,x i )

More information

Density Functional Theory for the Nonspecific Binding of Salt to Polyelectrolytes: Thermodynamic Properties

Density Functional Theory for the Nonspecific Binding of Salt to Polyelectrolytes: Thermodynamic Properties Biophysical Journal Volume 78 February 2000 699 706 699 Density Functional Theory for the Nonspecific Binding of Salt to Polyelectrolytes: Thermodynamic Properties Chandra N. Patra and Arun Yethiraj Department

More information

Suggestions for Further Reading

Suggestions for Further Reading Contents Preface viii 1 From Microscopic to Macroscopic Behavior 1 1.1 Introduction........................................ 1 1.2 Some Qualitative Observations............................. 2 1.3 Doing

More information

Nanoscale simulation lectures Statistical Mechanics

Nanoscale simulation lectures Statistical Mechanics Nanoscale simulation lectures 2008 Lectures: Thursdays 4 to 6 PM Course contents: - Thermodynamics and statistical mechanics - Structure and scattering - Mean-field approaches - Inhomogeneous systems -

More information

Atomistic molecular simulations for engineering applications: methods, tools and results. Jadran Vrabec

Atomistic molecular simulations for engineering applications: methods, tools and results. Jadran Vrabec Atomistic molecular simulations for engineering applications: methods, tools and results Jadran Vrabec Motivation Simulation methods vary in their level of detail The more detail, the more predictive power

More information

Effective charge saturation in colloidal suspensions

Effective charge saturation in colloidal suspensions JOURNAL OF CHEMICAL PHYSICS VOLUME 117, NUMBER 17 1 NOVEMBER 2002 Effective charge saturation in colloidal suspensions Lydéric Bocquet Laboratoire de Physique de l E.N.S. de Lyon, UMR CNRS 5672, 46 Allée

More information

The XY-Model. David-Alexander Robinson Sch th January 2012

The XY-Model. David-Alexander Robinson Sch th January 2012 The XY-Model David-Alexander Robinson Sch. 08332461 17th January 2012 Contents 1 Introduction & Theory 2 1.1 The XY-Model............................... 2 1.2 Markov Chains...............................

More information

VIII. Rubber Elasticity [B.Erman, J.E.Mark, Structure and properties of rubberlike networks]

VIII. Rubber Elasticity [B.Erman, J.E.Mark, Structure and properties of rubberlike networks] VIII. Rubber Elasticity [B.Erman, J.E.Mark, Structure and properties of rubberlike networks] Using various chemistry, one can chemically crosslink polymer chains. With sufficient cross-linking, the polymer

More information

Effect of surfactant structure on interfacial properties

Effect 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 information

Ab initio molecular dynamics and nuclear quantum effects

Ab initio molecular dynamics and nuclear quantum effects Ab initio molecular dynamics and nuclear quantum effects Luca M. Ghiringhelli Fritz Haber Institute Hands on workshop density functional theory and beyond: First principles simulations of molecules and

More information

Statistical 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 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 information

Phase Transition of Acrylamide-Based Polyampholyte Gels in Water

Phase Transition of Acrylamide-Based Polyampholyte Gels in Water Journal of Macromolecular Science w, Part A: Pure and Applied Chemistry, 43:1635 1649, 2006 Copyright # Taylor & Francis Group, LLC ISSN 1060-1325 print/1520-5738 online DOI: 10.1080/10601320600897056

More information

An Introduction to namic Light Scattering by Macromole cules

An Introduction to namic Light Scattering by Macromole cules An Introduction to namic Light Scattering by Macromole cules Kenneth S. Schmitz Department of Chemistry University of Missouri-Kansas Kansas City, Missouri City ACADEMIC PRESS, INC. Harcourt Brace Jovanovich,

More information

Polymer Solution Thermodynamics:

Polymer Solution Thermodynamics: Polymer Solution Thermodynamics: 3. Dilute Solutions with Volume Interactions Brownian particle Polymer coil Self-Avoiding Walk Models While the Gaussian coil model is useful for describing polymer solutions

More information

Investigation of Surfactant Efficiency Using Dissipative Particle Dynamics

Investigation of Surfactant Efficiency Using Dissipative Particle Dynamics Langmuir 2003, 19, 8195-8205 8195 Investigation of Surfactant Efficiency Using Dissipative Particle Dynamics Live Rekvig,*, Marieke Kranenburg, Jocelyne Vreede, Bjørn Hafskjold, and Berend Smit Department

More information

Experimental Soft Matter (M. Durand, G. Foffi)

Experimental 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 information

Chapter 5. Chemical potential and Gibbs distribution

Chapter 5. Chemical potential and Gibbs distribution Chapter 5 Chemical potential and Gibbs distribution 1 Chemical potential So far we have only considered systems in contact that are allowed to exchange heat, ie systems in thermal contact with one another

More information

Multimedia : Boundary Lubrication Podcast, Briscoe, et al. Nature , ( )

Multimedia : Boundary Lubrication Podcast, Briscoe, et al. Nature , ( ) 3.05 Nanomechanics of Materials and Biomaterials Thursday 04/05/07 Prof. C. Ortiz, MITDMSE I LECTURE 14: TE ELECTRICAL DOUBLE LAYER (EDL) Outline : REVIEW LECTURE #11 : INTRODUCTION TO TE ELECTRICAL DOUBLE

More information

Statistical Mechanics in a Nutshell

Statistical Mechanics in a Nutshell Chapter 2 Statistical Mechanics in a Nutshell Adapted from: Understanding Molecular Simulation Daan Frenkel and Berend Smit Academic Press (2001) pp. 9-22 11 2.1 Introduction In this course, we will treat

More information

Copyright 2001 University of Cambridge. Not to be quoted or copied without permission.

Copyright 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 information

Set the initial conditions r i. Update neighborlist. r i. Get new forces F i

Set the initial conditions r i. Update neighborlist. r i. Get new forces F i Set the initial conditions r i t 0, v i t 0 Update neighborlist Get new forces F i r i Solve the equations of motion numerically over time step t : r i t n r i t n + v i t n v i t n + Perform T, P scaling

More information

Thus, the volume element remains the same as required. With this transformation, the amiltonian becomes = p i m i + U(r 1 ; :::; r N ) = and the canon

Thus, the volume element remains the same as required. With this transformation, the amiltonian becomes = p i m i + U(r 1 ; :::; r N ) = and the canon G5.651: Statistical Mechanics Notes for Lecture 5 From the classical virial theorem I. TEMPERATURE AND PRESSURE ESTIMATORS hx i x j i = kt ij we arrived at the equipartition theorem: * + p i = m i NkT

More information

Current Opinion in Colloid & Interface Science

Current Opinion in Colloid & Interface Science Current Opinion in Colloid & Interface Science 13 (2008) 376 388 Contents lists available at ScienceDirect Current Opinion in Colloid & Interface Science journal homepage: www.elsevier.com/locate/cocis

More information

Opening the Gate: Framework Flexibility in ZIF-8 Explored by Experiments and Simulations

Opening the Gate: Framework Flexibility in ZIF-8 Explored by Experiments and Simulations Supporting Information for: Opening the Gate: Framework Flexibility in ZIF-8 Explored by Experiments and Simulations D. Fairen-Jimenez *, S. A. Moggach, M. T. Wharmby, P. A. Wright, S. Parsons and T. Düren

More information

IV. Transport Phenomena. Lecture 23: Ion Concentration Polarization

IV. Transport Phenomena. Lecture 23: Ion Concentration Polarization IV. Transport Phenomena Lecture 23: Ion Concentration Polarization MIT Student (and MZB) Ion concentration polarization in electrolytes refers to the additional voltage drop (or internal resistance ) across

More information

Collapse of Acrylamide-Based Polyampholyte Hydrogels in Water

Collapse of Acrylamide-Based Polyampholyte Hydrogels in Water Collapse of Acrylamide-Based Polyampholyte Hydrogels in Water Saadet Dogu, Melek Kilic, Oguz Okay Department of Chemistry, Istanbul Technical University, Istanbul, Turkey Received 5 February 2008; accepted

More information

My path: BS, Rose-Hulman PhD, UIUC Postdoc, Sandia. Map from Wikimedia Commons

My path: BS, Rose-Hulman PhD, UIUC Postdoc, Sandia. Map from Wikimedia Commons CoarseGrained Modeling of Ionomers and SaltDoped Block Copolymers Lisa Hall H.C. Slip Slider Assistant Professor William G. Lowrie Department of Chemical and Biomolecular Engineering The Ohio State University

More information

MODULE No. 24: Solution Kinetics Part - III

MODULE No. 24: Solution Kinetics Part - III Subject Paper No and Title Module No and Title Module Tag 6 and PHYSICAL CHEMISTRY-II (Statistical 24 and Solution Kinetics - III CHE_P6_M24 TABLE OF CONTENTS 1. Learning outcomes 2. Introduction 3. Primary

More information

Time-Dependent Statistical Mechanics 1. Introduction

Time-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 information

Supporting Information for: Complexation of β-lactoglobulin Fibrils and Sulfated Polysaccharides

Supporting Information for: Complexation of β-lactoglobulin Fibrils and Sulfated Polysaccharides Supporting Information for: Complexation of β-lactoglobulin Fibrils and Sulfated Polysaccharides Owen G Jones 1, Stephaandschin 1, Jozef Adamcik 1, Ludger Harnau 2, Sreenath Bolisetty 1, and Raffaele Mezzenga

More information

On the Chemical Free Energy of the Electrical Double Layer

On the Chemical Free Energy of the Electrical Double Layer 1114 Langmuir 23, 19, 1114-112 On the Chemical Free Energy of the Electrical Double Layer Marian Manciu and Eli Ruckenstein* Department of Chemical Engineering, State University of New York at Buffalo,

More information

Physics of disordered materials. Gunnar A. Niklasson Solid State Physics Department of Engineering Sciences Uppsala University

Physics of disordered materials. Gunnar A. Niklasson Solid State Physics Department of Engineering Sciences Uppsala University Physics of disordered materials Gunnar A. Niklasson Solid State Physics Department of Engineering Sciences Uppsala University Course plan Familiarity with the basic description of disordered structures

More information

CHEM-UA 652: Thermodynamics and Kinetics

CHEM-UA 652: Thermodynamics and Kinetics 1 CHEM-UA 652: Thermodynamics and Kinetics Notes for Lecture 4 I. THE ISOTHERMAL-ISOBARIC ENSEMBLE The isothermal-isobaric ensemble is the closest mimic to the conditions under which most experiments are

More information

Superparamagnetic nanoparticle arrays for magnetically tunable photonics. Josh Kurzman Materials 265

Superparamagnetic nanoparticle arrays for magnetically tunable photonics. Josh Kurzman Materials 265 Superparamagnetic nanoparticle arrays for magnetically tunable photonics Josh Kurzman Materials 265 Superparamagnetism In SPM regime, thermal energy sufficient to overcome spin reversal barrier T B Below

More information

Deformation behaviour of homogeneous and heterogeneous bimodal networks

Deformation behaviour of homogeneous and heterogeneous bimodal networks Deformation behaviour of homogeneous and heterogeneous bimodal networks Natasha Kamerlin 1,2 and Christer Elvingson 1, 1 Department of Chemistry - Ångström Laboratory, Uppsala University Box 523, S-751

More information

Principles and Applications of Electrochemistry

Principles and Applications of Electrochemistry Principles and Applications of Electrochemistry Fourth edition D. R. CROW Professor of Electrochemistry and Dean of Research University of Wolverhampton BLACKIE ACADEMIC & PROFESSIONAL An Imprint of Chapman

More information

Field Theoretic Simulations of Complex Coacervates

Field Theoretic Simulations of Complex Coacervates FieldTheoretic Simulations of Complex Coacervates Debra J. Audus Glenn H. Fredrickson CFDC 2009 Annual Meeting Contents Introduction/Motivation Theory and the Model Polyelectrolyte Phase Diagrams Theoretical

More information

On the calculation of solvation free energy from Kirkwood- Buff integrals: A large scale molecular dynamics study

On the calculation of solvation free energy from Kirkwood- Buff integrals: A large scale molecular dynamics study On the calculation of solvation free energy from Kirkwood- Buff integrals: A large scale molecular dynamics study Wynand Dednam and André E. Botha Department of Physics, University of South Africa, P.O.

More information

6, Physical Chemistry -II (Statistical Thermodynamics, Chemical Dynamics, Electrochemistry and Macromolecules)

6, Physical Chemistry -II (Statistical Thermodynamics, Chemical Dynamics, Electrochemistry and Macromolecules) Subject Paper No and Title Module No and Title Module Tag 6, Physical -II (Statistical Thermodynamics, Chemical Dynamics, Electrochemistry and Macromolecules) 25, Activity and Mean Activity coefficient

More information

Foundations of. Colloid Science SECOND EDITION. Robert J. Hunter. School of Chemistry University of Sydney OXPORD UNIVERSITY PRESS

Foundations of. Colloid Science SECOND EDITION. Robert J. Hunter. School of Chemistry University of Sydney OXPORD UNIVERSITY PRESS Foundations of Colloid Science SECOND EDITION Robert J. Hunter School of Chemistry University of Sydney OXPORD UNIVERSITY PRESS CONTENTS 1 NATURE OF COLLOIDAL DISPERSIONS 1.1 Introduction 1 1.2 Technological

More information

Supporting Information for Lysozyme Adsorption in ph-responsive Hydrogel Thin-Films: The non-trivial Role of Acid-Base Equilibrium

Supporting Information for Lysozyme Adsorption in ph-responsive Hydrogel Thin-Films: The non-trivial Role of Acid-Base Equilibrium Electronic Supplementary Material (ESI) for Soft Matter. This journal is The Royal Society of Chemistry 215 Supporting Information for Lysozyme Adsorption in ph-responsive Hydrogel Thin-Films: The non-trivial

More information

Explanation of Dramatic ph-dependence of Hydrogen Binding on Noble Metal Electrode: Greatly Weakened Water Adsorption at High ph.

Explanation of Dramatic ph-dependence of Hydrogen Binding on Noble Metal Electrode: Greatly Weakened Water Adsorption at High ph. Supplementary Materials Explanation of Dramatic ph-dependence of Hydrogen Binding on Noble Metal Electrode: Greatly Weakened Water Adsorption at High ph. Tao Cheng,, Lu Wang, Boris V Merinov, and William

More information

Exercise 2: Solvating the Structure Before you continue, follow these steps: Setting up Periodic Boundary Conditions

Exercise 2: Solvating the Structure Before you continue, follow these steps: Setting up Periodic Boundary Conditions Exercise 2: Solvating the Structure HyperChem lets you place a molecular system in a periodic box of water molecules to simulate behavior in aqueous solution, as in a biological system. In this exercise,

More information

Soft Matter and Biological Physics

Soft Matter and Biological Physics Dr. Ulrich F. Keyser - ufk20 (at) cam.ac.uk Soft Matter and Biological Physics Question Sheet Michaelmas 2011 Version: November 2, 2011 Question 0: Sedimentation Initially consider identical small particles

More information

NEGATIVE SPECIFIC HEAT IN A THERMODYNAMIC MODEL OF MULTIFRAGMENTATION

NEGATIVE SPECIFIC HEAT IN A THERMODYNAMIC MODEL OF MULTIFRAGMENTATION NEGATIVE SPECIFIC HEAT IN A THERMODYNAMIC MODEL OF MULTIFRAGMENTATION C. B. Das McGill University, Montréal, Canada Collaborators: S. Das Gupta and A. Z. Mekjian Plateau in caloric curve Singularity in

More information