Rare Event Simulations
|
|
- Charleen Goodwin
- 5 years ago
- Views:
Transcription
1 Rare Event Simulations
2 Homogeneous nucleation is a rare event (e.g. Liquid Solid) Crystallization requires supercooling (µ solid < µ liquid ) Crystal nucleus 2r Free - energy gain r 3 4! 3! GBulk = ñr! ì 3 Free - energy loss :! G Surface = 4!r 2 ã s,l s,l > < 0 0 r 2
3 !G tot!g * ~ r 2 ~ -r 3 0 r r c
4 Also, in many chemical reactions, a free-energy barrier separates reactants from products.
5 Classical chemical kinetics: A B
6 because
7 Statistical mechanics approach.
8 reactant product 1 A B 0 q q * g A (q-q * )
9 Common (not best) choice for g A (q-q * ): If we normalize the total concentration to 1, then <c A > = <g A >, and hence
10 Linear response theory (next slide) Because:
11 Intermezzo: linear response theory (ultra short).
12 Consider the response of an observable A due to an external field f B that couples to an observable B: For simplicity, assume that For small f B we can linearize:
13 Hence Now consider a weak field that is switched off at t=0. f B ΔA t 0
14 Using exactly the same reasoning as in the static case, we find: In the present case,
15 Time dependent linear response theory: We eliminate the external field, using then
16 The macroscopic equation was: And therefore
17 More convenient: consider at time derivative
18 Because in equilibrium: Hence: And so, for t=0:
19
20 Specific choice for g A : g q " q * = 1 " g q " q* =! q " q * B ( ) ( ) ( ) " g q # q* B ( ) " q A g q " q * = 1 "! q " q* =! q *" q A ( ) ( ) ( ) =! q # q * ( )
21 We get With
22 This can be rewritten as: Now the second term on the right is simply equal to For the first term on the right we have:
23 Our final result for k A B is Kinetic prefactor
24 In words: Rate= equilibrium probability to find system at the top of the barrier, multiplied by The cross-correlation of the velocity in the direction of the reaction coordinate at t=0 and the probability that the system has crossed the barrier at time t.
25 Transition state theory k A! B ( t) = = ( ) ( 0 *) " gb q # q q& ( 0) gb t " q c ( )! ( ) " " ( q* # q) A ( ) ( # ) ( ( )# ) q& 0 q 0 q* q t q * Let us consider the limit: t 0 + k TST A# B ( ( ) ) ( ) ( 0)! ( ( 0 ) $ *)" (&) " ( q * $ q) lim = # * = & " ( t) ( ) +! q t q! q t t 0 = q& q q q
26 The motion of the particle is ballistic The product or reactant states can equilibrate k TST A# B ( t) q& ( 0)! ( q( 0) $ q1) " ( q& )! ( q( 0) $ q1) = %! ( q( 0) $ q1) " ( q1 $ q) ( q( 0) & q1) ) dq ( q( 0) & q1) exp & U ( q) = ( q1 q) ) dq ( q( 0) & q1) exp & U ( q)!! $ ' " % ( # &! $ ' " % ( ( ) ( )! q( 0) = q 1 ( "# ( ) ( ) ( q1) ( ) exp $ &"! U q1 %' ( ) dq exp $ &"! U q %' ) ( ) ( ) ( ) Transition state theory ( ) ( ) ( ) ( & 1) $ & ( )% dq! ( q( 0) % q1) exp # & %" U ( q) $ ' dq! ( q( 0) % q1) exp #%" U ( q) $ q& 0! q 0 & q1 " q& dqq& 0! q 0 & q1 " q& exp $ ' &# U q % ( = & dq! q 0 q exp ' # U q ( ) ( ) = 0.5 q& 0 ( ) = 0.5 q& 0 ( ( & '
27 ( ) = 0.5 &( 0) k t q TST A" B q kbt = 2" m 1 ) #$ q 1 ) #$ ( ) exp %#! U q & ' 1 ( ( ) dq exp %' #! U q &( ( ) exp %#! U q & ' 1 ( ( ) dq exp %' #! U q &( depends on the choice for q 1.
28 What happens if you choose the wrong value for q *? Bennett-Chandler: less efficient, but still OK TST: wrong
29 Example: Ballistic barrier crossing
30
31 Transition state theory Requires accurate knowledge of free energy Yields an upper bound to the reaction rate Assumption: no re-crossings
32
33 k TST A# B k A# B ( t) ( t) ( $ ) (& 1 ) ( ( ) $ 1) ( $ q1) " ( q1 $ q) ( $ 1) ( ( ) $ 1) ( ( ) $ 1)! ( q( 0) $ q1) " ( q1 $ q) ( ) ( )! q( 0) q& 0! q 0 q " q! q 0 q = % ( ) ( ) q& 0! q 0 q " q t q! q 0 q = % Transmission coefficient! ( t) # k A" B TST ka" B ( t) ( 0) ( 0) ( # ) ( ( )# ) q! q q " q t q = & 0.5 q& ( ) Bennett Chandler Approach
34
35
36 Bennett-Chandler approach Results are independent of the precise location of the estimate of the transition state, but the accuracy does. If the transmission coefficient is very low Poor estimate of the reaction coordinate Diffuse barrier crossing
37 Simulating rare events: the diffusive limit reactant product ω A B q * q
38 How do concentration fluctuations decay? If the concentration increase on the left side of the barrier is ΔC A (0), then the concentration increase on the top of the barrier is, on the left side: And, on the right side: 0 ( absorbing boundary ) The diffusive flux across the barrier is:
39 Therefore The solution is with
40 In practice, the expression is only slightly more complex. To compute the rate of an activated process, we must compute: 1. The free energy barrier (ΔG) 2. The diffusion constant at the top of the barrier (D) 3. The effective width of the barrier (for a square barrier, this is W in other cases it is a bit different)
41 Diffusive barrier crossing Consider the following membrane TST: Flux = P(0) <v>/ω Bennett Chandler Flux= P(0) <v>/ω κ Diffusive barrier crossing Flux= P(0) D eff /ω 2
42 k A! B ( t) =! 0 ( ( 0 ) *) " gb q # q q& ( 0) gb t " q d = $ dt d t! = $ dt 0 q& c A ( ) ( ( 0 ) *) " gb q # q q& ( 0) gb t " q ( 0) c A ( ) ( ( 0 ) *) q& ( t ) " gb q # q " g " q " q " 0 reactant: q<0.5! # gb ( q) = $ q! on the barrier # % 1 product: q>0.5! k " 1 A# B ( t) = 2! $ 0 dt q& c A ( 0) q& ( t) * c A B Diffusive barrier crossing " 1 # g B $ barrier = %! # q $ & 0 elsewhere
43
44 Reaction coordinate F(q) r ( ) ( ) t : θ=1 q& 0 and! q& 0 For both ( ) ( ) ( # ) ( ( )# ) q& 0! q 0 q * " q t q * Low value of κ q
45 Computing free-energy landscapes Suppose we have an order parameter Q(r N ). The probability that Q(r N ) has a value Q is given by From this probability we can derive the variation of the free energy with Q:
46 This indicates a route to determine F(Q): just make a histogram of the spontaneous fluctuations of Q. P(Q) F(Q) Q Problem: F(Q) is very noisy, except near its minimum.
47 Application of biased MC simulation to determine P w (Q) But this we can rewrite as
48 Clearly, or We can choose the bias w(q) such that any desired range of Q-values is sampled. And we can correct for the bias.
49 The only remaining problem is that the different parts of F(Q) are shifted with respect to each other.
50 Solutions: 1. Fit to a single curve (not very elegant, but effective) 2. Use Ferrenberg-Swendsen scheme to combine different parts of the histogram (more elegant, but more sensitive to noise) Result:
51
52 SOLUTION 2: Use path sampling techniques The advantage of the latest path-sampling techniques (R.J. Allen et al. Phys. Rev. Lett. 94, 18104(2005)) is that they can be used for NON-EQUILIBRIUM systems.
53 Basic idea: 1. Define a reaction coordinate that distinguishes START from FINISH. 2. Compute the rate Γ 0 (low, but not very low) to reach surface Generate many trajectories from first crossing points and compute the fraction P(12) that reach 2 (without first returning to START) Repeat for all subsequent planes - and for many starting points. 5. Rate= Γ 0 P(12)P(23)..P(5F)
54 Illustration: protein crystallization.
55 !G tot!g * ~ r 2 ~ -r 3 0 r r c
56 Study nucleation events: 1. By experiment DIFFICULT BECAUSE NUCLEATION EVENTS ARE INFREQUENT, BUT FAST 2. By simulation IMPOSSIBLE WHAT TAKES LONG IN EXPERIMENTS, TAKES FOREVER [O(10 29 ) MD STEPS] IN A SIMULATION
57 SOLUTION 1: 1. Compute height of the free-energy barrier ΔG * (MC/MD) using biased sampling 2. Compute the transmission coefficient Γ ( using MD) Rate = " exp(#$g % /kt) Kinetic Prefactor (usually weak function of T) Probability of critical fluctuation (strong function of T)
58 PROTEIN STRUCTURE Problem: HUMAN GENOME PROJECT (?) protein sequences BUT WHAT IS THEIR 3D STRUCTURE? X-RAY CRYSTALLOGRAPHY REQUIRES GOOD CRYSTALS Only a small fraction of proteins have been successfully crystallized.
59 Crystallization of globular proteins QUOTE:... mainly trial and error... much like prospecting for gold... (McPherson Preparation and Analysis of Protein Crystals ) N
60 T Fluid F + S Solid Typical phase diagram of globular proteins T c ρ concentration Meta-stable fluid-fluid
61 M. Broide et al., PNAS 88,5660(1991) Phase diagram of GLOBULAR PROTEINS (γ-crystallin) T F F + S T c ρ
62 Protein crystallization conditions T F T c F + S F+F Protein concentration
63 # of crystalline particles Crystallization Condensation # of particles in a dense cluster (e.g. a droplet)
64 Crystal nucleation barriers
65 Nucleation speedup: 10 10
66 Experimental Evidence?? Experiments by Galkin and Vekilov PNAS 97, (2000) Lyzosyme crystallization
67 Liquid-liquid phase boundary Nucleation Rate (s -1 cm -3 ) Temperature (ºC)
68 Crystallization on structured templates If the crystal fits well on the template, it is easy to nucleate
69 SURPRISE:
70 RANDOM (!) SUBSTRATES ACT AS UNIVERSAL NUCLEATION SEEDS (Experiments: Prof. Naomi Chayen, Imperial College, London) Figure 1: Electron micrograph of micro-porous silicon (left) shows that the sample contains a network of pores with sizes between 5 and 10 nanometers. The figure on the right hand shows a number of protein crystals (phycobili protein) that have formed on the micro-porous silicon. (N.E. Chayen et al., Ref.. [1])
71 Explanation: For every Cinderella there is a shoe that fits
72 Theory: Richard Sear D. Frenkel, NATURE, 443, 641, 2006
73 Pore-facilitated vapour-liquid nucleation Off-lattice simulations (J.A. van Meel )
74 Pore-facilitated ``protein crystal nucleation. Almost epitaxial: crystallises too easily.
75 but with defects.
76 Pore-facilitated ``protein crystal nucleation. Far from epitaxial: no crystal at all.
77 Pore-facilitated ``protein crystal nucleation. Disordered substrate yields perfect crystal that nucleates in wetting layer.
78
(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 informationModeling the Free Energy Landscape for Janus Particle Self-Assembly in the Gas Phase. Andy Long Kridsanaphong Limtragool
Modeling the Free Energy Landscape for Janus Particle Self-Assembly in the Gas Phase Andy Long Kridsanaphong Limtragool Motivation We want to study the spontaneous formation of micelles and vesicles Applications
More informationRare Events. Transition state theory Bennett-Chandler Approach 16.2 Transition path sampling16.4
Rare Events Transition state theory 16.1-16.2 Bennett-Chandler Approach 16.2 Transition path sampling16.4 Outline Part 1 Rare event and reaction kinetics Linear Response theory Transition state theory
More informationNucleation rate (m -3 s -1 ) Radius of water nano droplet (Å) 1e+00 1e-64 1e-128 1e-192 1e-256
Supplementary Figures Nucleation rate (m -3 s -1 ) 1e+00 1e-64 1e-128 1e-192 1e-256 Calculated R in bulk water Calculated R in droplet Modified CNT 20 30 40 50 60 70 Radius of water nano droplet (Å) Supplementary
More informationCrystal nucleation for a model of globular proteins
JOURNAL OF CHEMICAL PHYSICS VOLUME 120, NUMBER 17 1 MAY 2004 Crystal nucleation for a model of globular proteins Andrey Shiryayev and James D. Gunton Department of Physics, Lehigh University, Bethlehem,
More informationNucleation: Theory and Applications to Protein Solutions and Colloidal Suspensions
Nucleation: Theory and Applications to Protein Solutions and Colloidal Suspensions Richard P. Sear Department of Physics, University of Surrey Guildford, Surrey GU2 7XH, United Kingdom r.sear@surrey.ac.uk
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 informationTemperature and Pressure Controls
Ensembles Temperature and Pressure Controls 1. (E, V, N) microcanonical (constant energy) 2. (T, V, N) canonical, constant volume 3. (T, P N) constant pressure 4. (T, V, µ) grand canonical #2, 3 or 4 are
More informationSUPPLEMENTARY INFORMATION
S1 Transition from WZ to ZB by addition of Zinc. By p-doping with zinc, the WZ fraction in the wires changes dramatically at a partial pressure of DEZn around 4.610-5 mbar. At high dopant levels fewer
More informationPhysical Biochemistry. Kwan Hee Lee, Ph.D. Handong Global University
Physical Biochemistry Kwan Hee Lee, Ph.D. Handong Global University Week 3 CHAPTER 2 The Second Law: Entropy of the Universe increases What is entropy Definition: measure of disorder The greater the disorder,
More informationOlivier Bourgeois Institut Néel
Olivier Bourgeois Institut Néel Outline Introduction: necessary concepts: phonons in low dimension, characteristic length Part 1: Transport and heat storage via phonons Specific heat and kinetic equation
More informationField Method of Simulation of Phase Transformations in Materials. Alex Umantsev Fayetteville State University, Fayetteville, NC
Field Method of Simulation of Phase Transformations in Materials Alex Umantsev Fayetteville State University, Fayetteville, NC What do we need to account for? Multi-phase states: thermodynamic systems
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 informationAnatoly B. Kolomeisky Department of Chemistry Center for Theoretical Biological Physics How to Understand Molecular Transport through Channels: The
Anatoly B. Kolomeisy Department of Chemistry Center for Theoretical Biological Physics How to Understand Molecular Transport through Channels: The Role of Interactions Transport Through Channels Oil pumping
More informationLecture 30: Kinetics of Epitaxial Growth: Surface Diffusion and
Lecture 30: Kinetics of Epitaxial Growth: Surface Diffusion and Nucleation Today s topics Understanding the basics of epitaxial techniques used for surface growth of crystalline structures (films, or layers).
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 informationTemperature and Pressure Controls
Ensembles Temperature and Pressure Controls 1. (E, V, N) microcanonical (constant energy) 2. (T, V, N) canonical, constant volume 3. (T, P N) constant pressure 4. (T, V, µ) grand canonical #2, 3 or 4 are
More information4/18/2011. Titus Beu University Babes-Bolyai Department of Theoretical and Computational Physics Cluj-Napoca, Romania
1. Introduction Titus Beu University Babes-Bolyai Department of Theoretical and Computational Physics Cluj-Napoca, Romania Bibliography Computer experiments Ensemble averages and time averages Molecular
More informationEquilibria in Materials
009 fall Advanced Physical Metallurgy Phase Equilibria in Materials 10. 13. 009 Eun Soo Park Office: 33-316 Telephone: 880-71 Email: espark@snu.ac.kr Office hours: by an appointment 1 Contents for previous
More informationSliding Friction in the Frenkel-Kontorova Model
arxiv:cond-mat/9510058v1 12 Oct 1995 to appear in "The Physics of Sliding Friction", Ed. B.N.J. Persson (Kluwer Academic Publishers), 1995 Sliding Friction in the Frenkel-Kontorova Model E. Granato I.N.P.E.,
More informationImportance Sampling in Monte Carlo Simulation of Rare Transition Events
Importance Sampling in Monte Carlo Simulation of Rare Transition Events Wei Cai Lecture 1. August 1, 25 1 Motivation: time scale limit and rare events Atomistic simulations such as Molecular Dynamics (MD)
More informationPrecipitation. Size! Shape! Size distribution! Agglomeration!
Precipitation Size! Shape! Size distribution! Agglomeration! Precipitation Four major questions: 1. Why do molecules/ions precipitate? 2. What determines the size? 3. What determines the size distribution?
More informationThermodynamics of phase transitions
Thermodynamics of phase transitions Katarzyna Sznajd-Weron Institute of Physics Wroc law University of Technology, Poland 7 Oct 2013, SF-MTPT Katarzyna Sznajd-Weron (WUT) Thermodynamics of phase transitions
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 informationFree energy simulations
Free energy simulations Marcus Elstner and Tomáš Kubař January 14, 2013 Motivation a physical quantity that is of most interest in chemistry? free energies Helmholtz F or Gibbs G holy grail of computational
More informationLocal Hyperdynamics. Arthur F. Voter Theoretical Division Los Alamos National Laboratory Los Alamos, NM, USA
Local Hyperdynamics Arthur F. Voter Theoretical Division National Laboratory, NM, USA Summer School in Monte Carlo Methods for Rare Events Division of Applied Mathematics Brown University Providence, RI
More informationDr.Abel MORENO CARCAMO Instituto de Química, UNAM. address:
Dr.Abel MORENO CARCAMO Instituto de Química, UNAM. E-mail address: carcamo@sunam.mx When? Understanding the macromolecular scale in time for crystal growth phenomena It is typical for crystal growth that
More informationMatSci 331 Homework 4 Molecular Dynamics and Monte Carlo: Stress, heat capacity, quantum nuclear effects, and simulated annealing
MatSci 331 Homework 4 Molecular Dynamics and Monte Carlo: Stress, heat capacity, quantum nuclear effects, and simulated annealing Due Thursday Feb. 21 at 5pm in Durand 110. Evan Reed In this homework,
More informationSwelling and Collapse of Single Polymer Molecules and Gels.
Swelling and Collapse of Single Polymer Molecules and Gels. Coil-Globule Transition in Single Polymer Molecules. the coil-globule transition If polymer chains are not ideal, interactions of non-neighboring
More informationHeterogenous Nucleation in Hard Spheres Systems
University of Luxembourg, Softmatter Theory Group May, 2012 Table of contents 1 2 3 Motivation Event Driven Molecular Dynamics Time Driven MD simulation vs. Event Driven MD simulation V(r) = { if r < σ
More informationAnalysis of the simulation
Analysis of the simulation Marcus Elstner and Tomáš Kubař January 7, 2014 Thermodynamic properties time averages of thermodynamic quantites correspond to ensemble averages (ergodic theorem) some quantities
More informationMolecular Dynamics Simulations of Fusion Materials: Challenges and Opportunities (Recent Developments)
Molecular Dynamics Simulations of Fusion Materials: Challenges and Opportunities (Recent Developments) Fei Gao gaofeium@umich.edu Limitations of MD Time scales Length scales (PBC help a lot) Accuracy of
More informationMolecular Dynamics and Accelerated Molecular Dynamics
Molecular Dynamics and Accelerated Molecular Dynamics Arthur F. Voter Theoretical Division National Laboratory Lecture 3 Tutorial Lecture Series Institute for Pure and Applied Mathematics (IPAM) UCLA September
More informationQuantum confined materials
Quantum confined materials Quantum Wells Quantum wells structures can be fabricated by seeral methods in semiconductors alloys. A quantum well is fabricated by depositing ery thin layers of different materials
More informationPhysics is time symmetric Nature is not
Fundamental theories of physics don t depend on the direction of time Newtonian Physics Electromagnetism Relativity Quantum Mechanics Physics is time symmetric Nature is not II law of thermodynamics -
More informationKinetic Monte Carlo: from transition probabilities to transition rates
Kinetic Monte Carlo: from transition probabilities to transition rates With MD we can only reproduce the dynamics of the system for 100 ns. Slow thermallyactivated processes, such as diffusion, cannot
More informationImportance Sampling in Monte Carlo Simulation of Rare Transition Events
Importance Sampling in Monte Carlo Simulation of Rare Transition Events Wei Cai Lecture 1. August 1, 25 1 Motivation: time scale limit and rare events Atomistic simulations such as Molecular Dynamics (MD)
More informationOptics of complex micro structures
Optics of complex micro structures dielectric materials λ L disordered partially ordered ordered random multiple scattering liquid crystals quasi crystals (Fibonacci) photonic crystals Assembly of photonic
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 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 informationEnhanced sampling via molecular dynamics II: Unbiased approaches
Enhanced sampling via molecular dynamics II: Unbiased approaches Mark E. Tuckerman Dept. of Chemistry and Courant Institute of Mathematical Sciences New York University, 100 Washington Square East, NY
More informationModeling the sputter deposition of thin film photovoltaics using long time scale dynamics techniques
Loughborough University Institutional Repository Modeling the sputter deposition of thin film photovoltaics using long time scale dynamics techniques This item was submitted to Loughborough University's
More information3.320 Lecture 23 (5/3/05)
3.320 Lecture 23 (5/3/05) Faster, faster,faster Bigger, Bigger, Bigger Accelerated Molecular Dynamics Kinetic Monte Carlo Inhomogeneous Spatial Coarse Graining 5/3/05 3.320 Atomistic Modeling of Materials
More 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 informationPhase transitions and finite-size scaling
Phase transitions and finite-size scaling Critical slowing down and cluster methods. Theory of phase transitions/ RNG Finite-size scaling Detailed treatment: Lectures on Phase Transitions and the Renormalization
More informationKinetics. Rate of change in response to thermodynamic forces
Kinetics Rate of change in response to thermodynamic forces Deviation from local equilibrium continuous change T heat flow temperature changes µ atom flow composition changes Deviation from global equilibrium
More informationMajor Concepts Lecture #11 Rigoberto Hernandez. TST & Transport 1
Major Concepts Onsager s Regression Hypothesis Relaxation of a perturbation Regression of fluctuations Fluctuation-Dissipation Theorem Proof of FDT & relation to Onsager s Regression Hypothesis Response
More informationBasics of Thermodynamics: Easy learning by Dr. Anjana Sen
Basics of Thermodynamics: Easy learning by Dr. Anjana Sen Part 1: Theory and concept Part 2: Definitions and equations Part 3: Laws of Thermodynamics Part 1: theory and concept Thermodynamics means conversion
More informationStep 1. Step 2. g l = g v. dg = 0 We have shown that over a plane surface of water. g v g l = ρ v R v T ln e/e sat. this can be rewritten
The basic question is what makes the existence of a droplet thermodynamically preferable to the existence only of water vapor. We have already derived an expression for the saturation vapor pressure over
More informationsubmitted to Physical Review E Numerical Tests of Nucleation Theories for the Ising Models Seunghwa Ryu 1 and Wei Cai 2
submitted to Physical Review E Numerical Tests of Nucleation Theories for the Ising Models Seunghwa Ryu 1 and Wei Cai 2 1 Department of Physics, Stanford University, Stanford, California 94305 2 Department
More informationQuantum Monte Carlo simulation of spin-polarized tritium
Higher-order actions and their applications in many-body, few-body, classical problems Quantum Monte Carlo simulation of spin-polarized tritium I. Bešlić, L. Vranješ Markić, University of Split, Croatia
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 informationMorphological evolution of single-crystal ultrathin solid films
Western Kentucky University From the SelectedWorks of Mikhail Khenner March 29, 2010 Morphological evolution of single-crystal ultrathin solid films Mikhail Khenner, Western Kentucky University Available
More informationVIII. Phase Transformations. Lecture 38: Nucleation and Spinodal Decomposition
VIII. Phase Transformations Lecture 38: Nucleation and Spinodal Decomposition MIT Student In this lecture we will study the onset of phase transformation for phases that differ only in their equilibrium
More informationModel Solutions Spring 2003
Exam 2 BE.462J/3.962J Model Solutions Spring 2003 (80 points total possible) 1. (10 points) Explain the phenomenon of phsensitive swelling in polyelectrolyte hydrogels. Why does the swelling depend on
More informationLesmahagow High School AHChemistry Inorganic and Physical Chemistry Lesmahagow High School CfE Advanced Higher Chemistry
Lesmahagow High School CfE Advanced Higher Chemistry Unit 1 Inorganic and Physical Chemistry Reaction Feasibility 1 Thermochemistry Thermochemistry is the study of energy changes in reactions. The First
More informationDependence of the surface tension on curvature
Dependence of the surface tension on curvature rigorously determined from the density profiles of nanodroplets Athens, st September M. T. Horsch,,, S. Eckelsbach, H. Hasse, G. Jackson, E. A. Müller, G.
More informationCollective Effects. Equilibrium and Nonequilibrium Physics
Collective Effects in Equilibrium and Nonequilibrium Physics: Lecture 3, 3 March 2006 Collective Effects in Equilibrium and Nonequilibrium Physics Website: http://cncs.bnu.edu.cn/mccross/course/ Caltech
More informationDensity Functional Theory: from theory to Applications
Density Functional Theory: from theory to Applications Uni Mainz May 27, 2012 Large barrier-activated processes time-dependent bias potential extended Lagrangian formalism Basic idea: during the MD dynamics
More informationOn the Reaction Rate Calculation Methods using Path Sampling
On the Reaction Rate Calculation Methods using Path Sampling Paris, France France 28 November, 2006 Titus S. van Erp Centre for Surface Chemistry and Catalysis, KULeuven Leuven, Belgium Introduction Straightforward
More informationAppendix 4. Appendix 4A Heat Capacity of Ideal Gases
Appendix 4 W-143 Appendix 4A Heat Capacity of Ideal Gases We can determine the heat capacity from the energy content of materials as a function of temperature. The simplest material to model is an ideal
More informationWhite, Marc A. (2012) A dynamic light scattering study of small organic molecule and ionic salt solutions through cooling crystallisation.
White, Marc A. (2012) A dynamic light scattering study of small organic molecule and ionic salt solutions through cooling crystallisation. MRes thesis http://theses.gla.ac.uk/3780/ Copyright and moral
More informationSuperfluidity and Superconductivity
Superfluidity and Superconductivity These are related phenomena of flow without resistance, but in very different systems Superfluidity: flow of helium IV atoms in a liquid Superconductivity: flow of electron
More informationSupplementary Information for Observation of dynamic atom-atom correlation in liquid helium in real space
3 4 5 6 7 8 9 0 3 4 5 6 7 8 9 0 Supplementary Information for Observation of dynamic atom-atom correlation in liquid helium in real space Supplementary Note : Total PDF The total (snap-shot) PDF is obtained
More informationDynamics of materials with X-ray Photon Correlation Spectroscopy - Opportunities and detector requirements
Dynamics of materials with X-ray Photon Correlation Spectroscopy - Opportunities and detector requirements Quasi-static speckles from colloidal suspension near random compact packing volume fraction Speckles
More informationDynamics of Nucleation in the Ising Model
J. Phys. Chem. B 2004, 108, 19681-19686 19681 Dynamics of Nucleation in the Ising Model Albert C. Pan and David Chandler* Department of Chemistry, UniVersity of California, Berkeley, California 94720-1460
More informationOrganization of NAMD Tutorial Files
Organization of NAMD Tutorial Files .1.1. RMSD for individual residues Objective: Find the average RMSD over time of each residue in the protein using VMD. Display the protein with the residues colored
More informationMulti-Ensemble Markov Models and TRAM. Fabian Paul 21-Feb-2018
Multi-Ensemble Markov Models and TRAM Fabian Paul 21-Feb-2018 Outline Free energies Simulation types Boltzmann reweighting Umbrella sampling multi-temperature simulation accelerated MD Analysis methods
More informationQuantum Corrections for Monte Carlo Simulation
Quantum Corrections for Monte Carlo Simulation Brian Winstead and Umberto Ravaioli Beckman Institute University of Illinois at Urbana-Champaign Outline Quantum corrections for quantization effects Effective
More informationarxiv:cond-mat/ v2 8 Jan 2004
On the Coalescence of Nanoscale Metal Clusters S. Hendy Applied Mathematics, Industrial Research Ltd, Lower Hutt, New Zealand S. A. Brown and M. Hyslop Nanostructure Engineering Science and Technology
More informationMultiscale Modeling of Epitaxial Growth Processes: Level Sets and Atomistic Models
Multiscale Modeling of Epitaxial Growth Processes: Level Sets and Atomistic Models Russel Caflisch 1, Mark Gyure 2, Bo Li 4, Stan Osher 1, Christian Ratsch 1,2, David Shao 1 and Dimitri Vvedensky 3 1 UCLA,
More informationEvaporation/Condensation of Ising Droplets
, Elmar Bittner and Wolfhard Janke Institut für Theoretische Physik, Universität Leipzig, Augustusplatz 10/11, D-04109 Leipzig, Germany E-mail: andreas.nussbaumer@itp.uni-leipzig.de Recently Biskup et
More informationTheoretical physics. Deterministic chaos in classical physics. Martin Scholtz
Theoretical physics Deterministic chaos in classical physics Martin Scholtz scholtzzz@gmail.com Fundamental physical theories and role of classical mechanics. Intuitive characteristics of chaos. Newton
More informationThe Phases of Hard Sphere Systems
The Phases of Hard Sphere Systems Tutorial for Spring 2015 ICERM workshop Crystals, Quasicrystals and Random Networks Veit Elser Cornell Department of Physics outline: order by disorder constant temperature
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 information16. Working with the Langevin and Fokker-Planck equations
16. Working with the Langevin and Fokker-Planck equations In the preceding Lecture, we have shown that given a Langevin equation (LE), it is possible to write down an equivalent Fokker-Planck equation
More informationSupporting information for: Anomalous Stability of Two-Dimensional Ice. Confined in Hydrophobic Nanopore
Supporting information for: Anomalous Stability of Two-Dimensional Ice Confined in Hydrophobic Nanopore Boxiao Cao, Enshi Xu, and Tianshu Li Department of Civil and Environmental Engineering, George Washington
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 informationGood Diffraction Practice Webinar Series
Good Diffraction Practice Webinar Series High Resolution X-ray Diffractometry (1) Mar 24, 2011 www.bruker-webinars.com Welcome Heiko Ress Global Marketing Manager Bruker AXS Inc. Madison, Wisconsin, USA
More informationDetermination of redox reaction rates and orders by in-situ liquid cell. electron microscopy of Pd and Au solution growth
BNL-107596-2014-JA Determination of redox reaction rates and orders by in-situ liquid cell electron microscopy of Pd and Au solution growth Eli A. Sutter * and Peter W. Sutter Center for Functional Nanomaterials,
More informationKinetic Monte Carlo: Coarsegraining
Kinetic Monte Carlo: Coarsegraining Time and Space Peter Kratzer Faculty of Physics, University Duisburg-Essen, Germany morphology Time and length scales Ga As 2D islands surface reconstruction Methods
More informationWang-Landau Monte Carlo simulation. Aleš Vítek IT4I, VP3
Wang-Landau Monte Carlo simulation Aleš Vítek IT4I, VP3 PART 1 Classical Monte Carlo Methods Outline 1. Statistical thermodynamics, ensembles 2. Numerical evaluation of integrals, crude Monte Carlo (MC)
More informationEntropic Crystal-Crystal Transitions of Brownian Squares K. Zhao, R. Bruinsma, and T.G. Mason
Entropic Crystal-Crystal Transitions of Brownian Squares K. Zhao, R. Bruinsma, and T.G. Mason This supplementary material contains the following sections: image processing methods, measurements of Brownian
More informationhydrate systems Gránásy Research Institute for Solid State Physics & Optics H-1525 Budapest, POB 49, Hungary László
Phase field theory of crystal nucleation: Application to the hard-sphere and CO 2 Bjørn Kvamme Department of Physics, University of Bergen Allégaten 55, N-5007 N Bergen, Norway László Gránásy Research
More information1.3 Molecular Level Presentation
1.3.1 Introduction A molecule is the smallest chemical unit of a substance that is capable of stable, independent existence. Not all substances are composed of molecules. Some substances are composed of
More informationThe first three categories are considered a bottom-up approach while lithography is a topdown
Nanowires and Nanorods One-dimensional structures have been called in different ways: nanowires, nanorod, fibers of fibrils, whiskers, etc. The common characteristic of these structures is that all they
More informationPower Law of Molecular Weight of the Nucleation Rate of Folded Chain Crystals of Polyethylene
Macromolecules 2002, 35, 6985-6991 6985 Power Law of Molecular Weight of the Nucleation Rate of Folded Chain Crystals of Polyethylene Swapan K. Ghosh, Masamichi Hikosaka,*, Akihiko Toda, Shinichi Yamazaki,
More informationModified Maxwell-Garnett equation for the effective transport coefficients in inhomogeneous media
J. Phys. A: Math. Gen. 3 (998) 7227 7234. Printed in the UK PII: S0305-4470(98)93976-2 Modified Maxwell-Garnett equation for the effective transport coefficients in inhomogeneous media Juris Robert Kalnin
More informationAGuideto Monte Carlo Simulations in Statistical Physics
AGuideto Monte Carlo Simulations in Statistical Physics Second Edition David P. Landau Center for Simulational Physics, The University of Georgia Kurt Binder Institut für Physik, Johannes-Gutenberg-Universität
More informationA droplet of colloidal solution is left to evaporate on a superhydrophobic surface. Avijit Baidya
A droplet of colloidal solution is left to evaporate on a superhydrophobic surface. Avijit Baidya 14.03.15 In this paper Evaporation-driven particle self-assembly can be used to generate three-dimensional
More informationLiquid-vapour oscillations of water in hydrophobic nanopores
Alicante 2003 4th European Biophysics Congress Liquid-vapour oscillations of water in hydrophobic nanopores 7 th July 2003 Oliver Beckstein and Mark S. P. Sansom Department of Biochemistry, Laboratory
More informationChE 210B: Advanced Topics in Equilibrium Statistical Mechanics
ChE 210B: Advanced Topics in Equilibrium Statistical Mechanics Glenn Fredrickson Lecture 1 Reading: 3.1-3.5 Chandler, Chapters 1 and 2 McQuarrie This course builds on the elementary concepts of statistical
More informationGraupel and Hail Growth
Graupel and Hail Growth I. Growth of large ice particles In this section we look at some basics of graupeln and hail growth. Important components of graupeln and hail growth models include production of
More informationBiomolecular modeling. Theoretical Chemistry, TU Braunschweig (Dated: December 10, 2010)
Biomolecular modeling Marcus Elstner and Tomáš Kubař Theoretical Chemistry, TU Braunschweig (Dated: December 10, 2010) IX. FREE ENERGY SIMULATIONS When searching for a physical quantity that is of most
More informationThe Coulomb Barrier not Static in QED A correction to the Theory by Preparata on the Phenomenon of Cold Fusion and Theoretical hypothesis
Frisone, F. "The Coulomb Barrier not Static in QED," A correction to the Theory by Preparata on the Phenomenon of Cold Fusion and Theoretical hypothesis. in ICCF-14 International Conference on Condensed
More informationLECTURE 11: Monte Carlo Methods III
1 LECTURE 11: Monte Carlo Methods III December 3, 2012 In this last chapter, we discuss non-equilibrium Monte Carlo methods. We concentrate on lattice systems and discuss ways of simulating phenomena such
More informationTopological defects and its role in the phase transition of a dense defect system
Topological defects and its role in the phase transition of a dense defect system Suman Sinha * and Soumen Kumar Roy Depatrment of Physics, Jadavpur University Kolkata- 70003, India Abstract Monte Carlo
More informationIn situ studies on dynamic properties of carbon nanotubes with metal clusters
In situ studies on dynamic properties of carbon nanotubes with metal clusters Jason Chang, Yuan-Chih Chang, Der-Hsien Lien, Shaw-Chieh Wang*, Tung Hsu*, and Tien T. Tsong Institute of Physics, Academia
More informationMolecular Machines and Enzymes
Molecular Machines and Enzymes Principles of functioning of molecular machines Enzymes and catalysis Molecular motors: kinesin 1 NB Queste diapositive sono state preparate per il corso di Biofisica tenuto
More informationSupporting Information: On Localized Vapor Pressure Gradients Governing Condensation and Frost Phenomena
Supporting Information: On Localized Vapor Pressure Gradients Governing Condensation and Frost Phenomena Saurabh Nath and Jonathan B. Boreyko Department of Biomedical Engineering and Mechanics, Virginia
More information