IMPROVED METHOD FOR CALCULATING SURFACE TENSION AND APPLICATION TO WATER
|
|
- Elwin Willis
- 6 years ago
- Views:
Transcription
1 IMPROVED METHOD FOR CALCULATING SURFACE TENSION AND APPLICATION TO WATER ABSTRACT Hong Peng 1, Anh V Nguyen 2 and Greg R Birkett* School of Chemical Engineering, The University of Queensland Brisbane, Queensland 472, Australia 1 Author s h.peng2@uq.edu.au 2 Author s a.nguyen1@uq.edu.au *Corresponding author: g.birkett@uq.edu.au We have proposed an improved method for calculating surface tension of liquid-vapour interface based on the thermodynamic route. This method combines the central difference scheme used in the test-area method and the Bennett weighting factor. This improved method shows faster convergence than the test area method as the weighting factor reduces the asymmetric bias of two perturbations. In addition, this method is used to calculate the surface tension of water-methanol mixtures. The values obtained agree with experimental data and other simulation data. This method provides the new surface tension calculation procedure which is independent with on evaluation of complex pressure tensor when using the classical Kirkwood and Buff method. INTRODUCTION Surface tension is central to many daily and industrial activities, ranging from washing with soap, foaming in kitchen sinks to bubble coalescence in chemical and biological reactors. The measurement of surface tension is now well established. However, the determination of surface tension by molecular simulation still remains a challenge to scientists. There are two general routes employed to determine surface tension by molecular simulation. The first and most widely used method is based on the mechanical definition of surface tension. The mechanical route evaluates surface tension, γ, using the normal, P N, and tangential, P T, pressure tensors and is given in equation(1). ( ) PN PT dz (1) γ = This equation was first introduced by Bakker in 1928 (Bakker, 1928) but is often referred to as the Kirkwood and Buff method (Kirkwood and Buff, 1949). The mechanical route has proven useful for a number of applications but the evaluation of the pressure tensors can be cumbersome and is a large overhead in Monte Carlo simulations. The second route is to use the thermodynamic definition of surface tension. In the canonical ensemble (where the temperature, T, volume, V, and number of particles, N, are held constant), the thermodynamic definition is given in equation (2) along with the finite difference approximation. Page 1
2 F F γ = = lim A A A NVT NVT (2) where F is the Helmholtz free energy and A is the interfacial area of the system. The accuracy of the finite difference approximation in equation (2) comes down to two factors: 1. the accuracy of the finite difference approximation for a given system and step size and 2. The accuracy of the method used to evaluate the free energy difference between the states (which differ in this case by their interfacial surface area). The difference in the free energy between any two systems can be calculated using the partition functions, Q and Q 1, for systems and 1, respectively, as per equation(3). 1 Q ln F 1 = F1 F = β Q1 (3) where β is given by 1/ ( kbt ). The canonical partition function is dependent on the configurational energy, U, and the configurational space, r N, and is described as 1 Q = β N! ( ) N exp U d 3N Λ r (4) where Λ is the Thermal de Broglie wavelength. Various methods can be used for calculating the free energy difference in equation(4). For example, the Bennett method (BM) (Bennett, 1976) the free-energy difference is obtained by considering two equilibrated systems, denoted as and 1, with the same values for N, V and T, but different interfacial areas, A and A 1, respectively. A perturbation in interfacial area for system from A to A 1, whilst maintaining all values for N, V and T, is denoted as 1 and referred to as a forward perturbation. And for system 1, a change from A 1 to A, is denoted as 1 and referred to as the reverse perturbation. The difference in configurational energy between system, U, and the perturbed system 1, U 1, is the change in the potential energy due to the perturbation, U 1. In a canonical ensemble, Bennett (Bennett, 1976) showed that the ratio of the partition functions can be calculated using equation (5). N ( r ) exp( βu' ) N ( r ) exp( β ) Q W = Q W U 1 1 1' (5) where the angle brackets with the subscripts and 1 describe the canonical ensemble averages of the perturbation steps for equilibrated systems and 1, respectively, and ( ) N W r is a weighting function of some description. A prudent choice in the weighting function can improve the accuracy of the calculation in equation (5) by maximizing the weighting of the calculation in the area of configuration space where the forward and reverse perturbations overlap. Bennett investigated what form of the weighting function was optimal and found that using the Fermi-Dirac function, f(x) = (1+exp(x)) -1 together Page 2
3 with a shifted potential was optimal. Implementing this in equation(5), gives equation(6) which is the basis of the Bennett method. ( β 1 ' + ) ( β ) 1 Q 1 ln ln f U C C β β β 1 F 1 = = + Q1 f U 1' C (6) where the and C is an adjustable parameter used to maximize the overlap of the forward and revere perturbations. The optimal value for C requires knowledge of the overlap between the two simulations and is given by equation(7). C ln Q = (7) Q1 We cannot know the result of equation (7) before we start a simulation, however, C can be estimated with increasing accuracy during the simulation using equation(7). It may seem inconsistent to use an equation that requires C, to calculate C, it must be remembered that a property of equation(7) is that it will converge for any value of C. A suitable value for C decreases convergence time and better estimates of C will be calculated as the simulation progresses. The Bennett method has been rarely used to evaluate the surface tension due to the need of two more simulations (Errington and Kofke, 27) although this procedure has been widely used to calculate other properties in simulations (Lu et al., 23, Harismiadis et al., 1996, Lu and Kofke, 21). Another method for estimating the free energy difference is the Zwanzig theory for high-temperature perturbation, as given in equation(8) (Zwanzig, 1954). This equation is for a single system which is perturbed to system 1. 1 F 1 = ln exp ( β U ) (8) β The test-area method (Gloor et al., 25), uses equation (8) to calculate the free energy difference for perturbations that increase (1) and decrease (-1) the interfacial area from a specific reference system by A and A, respectively. The resulting free-energy differences, F 1 and F -1, are averaged and subbed into 2 as the estimate for F to give equation(9). F = F 2 A 1 1 γ (9) The test-area method has been used to calculate the surface tension of simple Lenard-Jones (LJ) fluids (Gloor et al., 25, Errington and Kofke, 27), water (Vega and Miguel, 27, Ghoufi et al., 28) and related mixtures (Bahadur et al., 27). It is noted that the test-area scheme only requires a single simulation to determine surface tension. However, under some special conditions the test-area method can provide less accurate surface tension than the Bennett method. Examples include the cases where asymmetric bias is inherent in the simulation of free-energy differences (e.g., the square-well fluid) (Errington and Kofke, 27). This problem with averaging forward and reverse free energy differences has also been discussed by Lu et al (Errington and Kofke, 27, Lu et al., 23) who recommend that this method should be avoided due to the unpredictable consequences of combing the free energy changes this way. In this paper, we propose an improved method by using the advantages of the Bennett method (improved accuracy and precision) and the Test Area method (requires Page 3
4 a single simulation). The proposed scheme is shown schematically in Table 1. The key assumption to the proposed method is that F 1 = F -1 such that the two perturbations have equivalent weights in the Bennett method. Substituting this into equation (6) and using the slightly neater (but equivalent) weighting factor used by Lu et al. (Lu et al., 23), gives equation (1), which is used in this study to calculate surface tension. Table 1 The Bennett method, the test-area method and the improved method proposed in this paper. The state (-1) presents the negative (compression) perturbation, while the state (1) is the positive (expansion) perturbation, all relative to the reference state (). Dash line square reflects the virtual perturbation state. W represents the weighting function. Method Perturbation Scheme Test area method -1 1 Bennett method w 1 w 1 A single simulation is performed in which a reference system equilibrates and then samples configurations relevant to an interfacial area A. An area change, A, is then performed to generate a virtual perturbed state (System 1, expansion of interfacial area: A 1 = A (1+ A), or System -1, compression of interfacial area: A- 1 = A (1- A). The perturbation in area is made along x and y axis by equal ratio, L = L 1±.5 A. To keep the volume constant, the corresponding change ( x, y) 1, 1 ( x, y) ( ) ( ) along z axis is, L ( ) L ( ) ( A) = 1±.5.The corresponding energy differences U ( z) 1,1 ( z) 1 and U -1 are obtained by perturbation steps are repeated many times through the simulation to collect the ensemble averages in equation (1) to calculated the surface tension. ( ) exp ( β / 2) ( ) ( β ) 1 W U U -1 ln γ = (1) β A W U exp U / 2 1 The weighting function W ( U) in the equation (1) is described by a Gaussianlike hyperbolic secant function, i.e., ( ) = ( β ) Improved method w -1 1 W U cosh U C 2 (11) Applying the Bennett equation (6) for determining the free-energy differences of the forward (expansion) and reverse (compression) perturbations, F 1 and F 1, and employing the central difference equation, i.e. equation (9) of the test-area method give. Page 4 w
5 The required Bennett factor can be determined at the crossing point of the probability distribution functions of free-energy perturbations for the forward and reverse directions in the Bennett method. The crossing point is selected to gain the greatest overlap of samplings and maximize the weighting on the configurational energy differences, giving ( ) exp ( β / 2) ( ) exp ( β / 2) W U U -1 C = ln W U U 1 (12) Equation (1) in conjunction with equations (11) and (12) was used to calculate surface tension. The improved method and the test-area method are used to evaluate the surface tension of water-vapour and water-methanol mixtures using an in-house Monte Carlo (MC) simulation code. METHODS Interaction Energy Models Isotropic site-site pair-wise potentials were used in this study. The interaction potential between two molecules is: 12 6 Α Β a b cd cd q C D i q j σ cd ij σ ij φij = + 4 ab ε ij cd cd a= 1 b= 1 4πε rij c= 1 d = 1 r ij r ij (13) where ϕ ij is the interaction energy between fluid molecules i and j, Α and Β are the number of charges on the molecules i and j, respectively, C and D are the number of LJ sites on the molecules i and j, respectively, ε is the permittivity of a vacuum, r ij ab is the separation between the charge a on molecule i and the charge b on molecule j having charges q i a and q j b, respectively and r ij cd is the separation between the LJ site c on molecule i and the LJ site d on molecule j with combined LJ well depth of ε ij cd and a combined LJ collision diameter of σ ij cd for the two sites. The combined LJ well-depth and collision diameter of two sites is calculated based on the Lorentz-Berthelot mixing rules (Potoff et al., 1999). For water models, the TIP4P/25 (Abascal and Vega, 25) model was mainly used. It consists of one LJ site located at the oxygen and three coulomb sites; two positive charges located at the two hydrogen sites and a negative charge placed along the bisector of the H-O-H angle (the M site). The distance of OH bond is.9572å and the angle of H-O-H is Both LJ and electrostatic interactions were truncated at 1.3 nm. No long-range corrections were used for either LJ or electrostatic interactions. The OPLS model described by Jorgensen (Jorgensen et al., 1993) is used to simulate the methanol molecules. The mole of fraction of methanol in the liquid water varied from % to 1% with interval value of 25%. The ratio of number of molecules is listed in Table 2 (Chang and Dang, 25). The temperature used in simulation is 3K. Page 5
6 Table 2. Simulation Mole fraction of water-methanol mixture Fraction (%) N water N methanol Simulation Procedure The simulations were performed in the canonical NVT ensemble by using an in-house Monte Carlo code. For each simulation, a two-stage run was used to obtain the properties of the system. Firstly, a cubic box of the bulk liquid was run and equilibrated by at least 2, cycles (each cycle consists of number of molecules moves). The second stage was specially designed to determine the interface properties; surface tension and the local density distribution. The equilibrium bulk box of liquid from the last run of the first stage was used for the initial configuration. To make the interface of vapor-liquid, two empty boxes with the same geometry, were added on each side of bulk box of liquid along z-axis (Figure 1). This new system with two interfaces is equilibrated for 2, cycles following by 4, cycles to produce constant (reliable) values for surface tension. The last 1, cycles are divided 5 blocks to obtain the mean values and standard deviation for each simulation condition. Figure 1. A schematic of the simulation box at the second computation stage with the liquid cubic slab sandwiched between two vapor phase cubic boxes. RESULTS y x L L L Convergence of thermodynamic route method There is the issue in simulations of inhomogeneous system is whether the runs have been performed over a sufficiently long period time to obtain reliable value of surface tension. From our simulation in Figure 2, it is shown that the value at the first 1, cycles fluctuated followed a slow convergence trend which ranged from 2, cycles to 4, cycles. The surface tension profile by improved method proposed shows the much earlier convergence just at near 2, cycles than the value by test-area method which converged near the end of simulation around 4, cycles. z Page 6
7 Surface tension[nn/ m] Improved method Test area x 1 4 Monte Carlo Cycles Figure 2. Surface tension versus Monte Carlo simulation cycles as obtained by the testarea method and the improved method developed in this paper. Effect of box size on surface tension In the molecular simulation, box size is the critical parameter on the accuracy of surface tension. We vary the simulation box size from 2.7 nm to 3.5 nm with the number of particles ranging from 656 to 1428 at 3K. As shown in Figure 3, for the box size more than 3.nm, the average values of surface tension are similar and close with the reported value (Vega and Miguel, 27). When box size is larger than 3nm, the effect of system size is limited Surface tension [mn/ m] Box size L[nm] Figure 3. Surface tension versus box sizes as obtained by the improved method developed in this paper. Page 7
8 Surface tension of water-vapour phase To verify our approach could be applied for water, improved method is used to calculate its surface tension at two temperatures. The simulated values of surface tension are summarized in Table 3 in comparison with test-area method and the reported values (Vega and Miguel, 27). In our simulation, no long range correction for both LJ and CL interactions is applied but the simulation data is reasonable consistent with the values with long range correction. Maybe, the large cutoff distance in the simulation ensures that the effect of tail correction is small. The relative errors of surface tension by improved method for both temperatures are smaller than the value from test area method. This could be attributed to the use of weight function optimizing the two perturbation ways. Table 3. Simulation value of surface tension Surface Tension (mn/m) T (K) Improved method Test-area method Experiments ± ± ± ± Surface tension of water-methanol mixtures We also use the improved method to calculate the surface tension of water methanol mixtures. The mole of fraction of methanol in the liquid water varied from % to 1% with interval value of 25%. The calculated surface tension is shown in Figure 4 compared with experimental value (Vazquez et al., 1995) and the simulation values by Chang and Dang (Chang and Dang, 25). The calculated value is expected higher than experimental value and close to the MD model except that pure water value which is lower than both experimental and MD simulation data. As we discussed before, the rigid non-polarization water model only could provide the value lower than experimental data. In general, the overall dependence of the surface tension on mixture composition agrees well with the experimental trend. When Liquid water is mixed with the amount of methanol molecular, the surface tension of mixture decreased with the elevated the mole fraction of methanol. It may be caused by the strong tendency for methanol to concentrate at the interface. 1 Surface tension [mn/ m] Experiment This study MD Methanol concentration in water [%] Page 8
9 Figure 4. Surface tension of water-methanol mixtures at 3 K. Experimental data are shown with continuous line (Vazquez et al., 1995). MD results are taken from (Chang and Dang, 25). CONCLUSION Here, we proposed an improved method on evaluation of surface tension of liquidvapour interface based on thermodynamic route. This method combined the central difference scheme used the test-area method and the weighting factor used in the Bennett technique. This improved method shows faster convergence than the test area method as the weighting factor reduces the asymmetric bias of two perturbations. The value of surface tension of water-vapour obtained by this improved method agrees with the expected values within the small uncertainty. In addition, this method is also used to evaluate the surface tension of water-methanol mixtures and obtained the close value with experimental data and other simulation data. REFERENCES ABASCAL, J. L. F. & VEGA, C. (25) A general purpose model for the condensed phases of water: TIP4P/25. The Journal of Chemical Physics, 123, BAHADUR, R., RUSSELL, L. M. & ALAVI, S. (27) Surface Tensions in NaCl- Water-Air Systems from MD Simulations. The Journal of Physical Chemistry B, 111, BAKKER, G. (1928) Kapillarität und Oberflächenspannung, Leipzig, Akad. Verlagsgesellschaft. BENNETT, C. H. (1976) Efficient estimation of free-energy differences from montecarlo data. Journal of Computational Physics, 22, CHANG, T.-M. & DANG, L. X. (25) Liquid-Vapor Interface of Methanol-Water Mixtures: A Molecular Dynamics Study. The Journal of Physical Chemistry B, 19, ERRINGTON, J. R. & KOFKE, D. A. (27) Calculation of surface tension via area sampling. The Journal of Chemical Physics, 127, GHOUFI, A., GOUJON, F., LACHET, V. & MALFREYT, P. (28) Surface tension of water and acid gases from Monte Carlo simulations. The Journal of Chemical Physics, 128, GLOOR, G. J., JACKSON, G., BLAS, F. J. & MIGUEL, E. D. (25) Test-area simulation method for the direct determination of the interfacial tension of systems with continuous or discontinuous potentials. The Journal of Chemical Physics, 123, HARISMIADIS, V., VORHOLZ, J. & PANAGIOTOPOULOS, A. (1996) Efficient pressure estimation in molecular simulations without evaluating the virial. Journal of Chemical Physics, 15, JORGENSEN, W., MAXWELL, D. & TIRADO-RIVES, J. (1993) Development and testing of the OPLS all-atom force field on conformational energetics and properties of organic liquids. Q. ReV. Biophys, 26, 49. KIRKWOOD, J. G. & BUFF, F. P. (1949) The Statistical Mechanical Theory of Surface Tension. The Journal of Chemical Physics, 17, LU, N. & KOFKE, D. A. (21) Accuracy of free-energy perturbation calculations in molecular simulation. I. Modeling. Journal of Chemical Physics, 114, Page 9
10 LU, N., SINGH, J. & KOFKE, D. (23) Appropriate methods to combine forward and reverse free-energy perturbation averages. The Journal of Chemical Physics, 118, POTOFF, J. J., ERRINGTON, J. R. & PANAGIOTOPOULOS, A. Z. (1999) MOLECULAR PHYSICS LECTURE Molecular simulation of phase equilibria for mixtures of polar and non-polar components. Molecular Physics: An International Journal at the Interface Between Chemistry and Physics, 97, VAZQUEZ, G., ALVAREZ, E. & NAVAZA, J. M. (1995) Surface Tension of Alcohol Water + Water from 2 to 5.degree.C. Journal of Chemical & Engineering Data, 4, VEGA, C. & MIGUEL, E. D. (27) Surface tension of the most popular models of water by using the test-area simulation method. The Journal of Chemical Physics, 126, ZWANZIG, R. W. (1954) High-Temperature Equation of State by a Perturbation Method. I. Nonpolar Gases. The Journal of Chemical Physics, 22, Brief Biography of Presenter Hong Peng, is the Ph.D candidate in school of chemical Engineering, The university of Queensland. His main project is about molecular simulation on interface phenomena of water and its mixtures with ions and surfactants cross liquid-vapour-solid phases. Before taken this Ph.D scholarship in Uni, he has worked in BHPBilliton technology and Olympic Dam for three years on various hydrometallurgy process development of titanium, nickel and uranium immediately after obtained his master degree in bioleaching area. Page 1
Supporting Information for: Physics Behind the Water Transport through. Nanoporous Graphene and Boron Nitride
Supporting Information for: Physics Behind the Water Transport through Nanoporous Graphene and Boron Nitride Ludovic Garnier, Anthony Szymczyk, Patrice Malfreyt, and Aziz Ghoufi, Institut de Physique de
More informationCE 530 Molecular Simulation
1 CE 530 Molecular Simulation Lecture 14 Molecular Models David A. Kofke Department of Chemical Engineering SUNY Buffalo kofke@eng.buffalo.edu 2 Review Monte Carlo ensemble averaging, no dynamics easy
More informationAn Introduction to Two Phase Molecular Dynamics Simulation
An Introduction to Two Phase Molecular Dynamics Simulation David Keffer Department of Materials Science & Engineering University of Tennessee, Knoxville date begun: April 19, 2016 date last updated: April
More informationCE 530 Molecular Simulation
CE 530 Molecular Simulation Lecture 20 Phase Equilibria David A. Kofke Department of Chemical Engineering SUNY Buffalo kofke@eng.buffalo.edu 2 Thermodynamic Phase Equilibria Certain thermodynamic states
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 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 information510 Subject Index. Hamiltonian 33, 86, 88, 89 Hamilton operator 34, 164, 166
Subject Index Ab-initio calculation 24, 122, 161. 165 Acentric factor 279, 338 Activity absolute 258, 295 coefficient 7 definition 7 Atom 23 Atomic units 93 Avogadro number 5, 92 Axilrod-Teller-forces
More informationMolecular Simulation Study of the Vapor Liquid Interfacial Behavior of a Dimer-forming Associating Fluid
Molecular Simulation, Vol. 30 (6), 15 May 2004, pp. 343 351 Molecular Simulation Study of the Vapor Liquid Interfacial Behavior of a Dimer-forming Associating Fluid JAYANT K. SINGH and DAVID A. KOFKE*
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 informationSurface analysis algorithms in the mardyn program and the ls1 project
Surface analysis algorithms in the mardyn program and the ls1 project Stuttgart, 15 th December 1 M. T. Horsch Surface tension The virial route Bakker-Buff equation: γ R 2 out in dz z Normal pressure decays
More informationSupplemental 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 informationControlling the Long-Range Corrections in Atomistic Monte Carlo Simulations of Two-Phase Systems
pubs.acs.org/jctc Controlling the Long-Range Corrections in Atomistic Monte Carlo Simulations of Two-Phase Systems Florent Goujon, Aziz Ghoufi, Patrice Malfreyt,*, and Dominic J. Tildesley Institut de
More informationSurface thermodynamics of planar, cylindrical, and spherical vapour-liquid interfaces of water
Surface thermodynamics of planar, cylindrical, and spherical vapour-liquid interfaces of water Gabriel V. Lau, Ian J. Ford, Patricia A. Hunt, Erich A. Müller, and George Jackson Citation: The Journal of
More informationEvaluation of the pressure tensor and surface tension for molecular fluids with discontinuous potentials using the volume perturbation method
Evaluation of the pressure tensor and surface tension for molecular fluids with discontinuous potentials using the volume perturbation method Guadalupe Jiménez-Serratos, Carlos Vega, and Alejandro Gil-Villegas
More informationUB association bias algorithm applied to the simulation of hydrogen fluoride
Fluid Phase Equilibria 194 197 (2002) 249 256 UB association bias algorithm applied to the simulation of hydrogen fluoride Scott Wierzchowski, David A. Kofke Department of Chemical Engineering, University
More informationclassical molecular dynamics method
Numerical study on local pressure c Titletensions of liquid film in the vici classical molecular dynamics method Author(s) 藤原, 邦夫 Citation Issue Date Text Version ETD URL http://hdl.handle.net/11094/52207
More informationSurface tension of the most popular models of water by using the test-area simulation method
THE JOURNAL OF CHEMICAL PHYSICS 126, 154707 2007 Surface tension of the most popular models of water by using the test-area simulation method C. Vega a Departamento de Química Física, Facultad de Ciencias
More informationThermodynamic behaviour of mixtures containing CO 2. A molecular simulation study
Thermodynamic behaviour of mixtures containing. A molecular simulation study V. Lachet, C. Nieto-Draghi, B. Creton (IFPEN) Å. Ervik, G. Skaugen, Ø. Wilhelmsen, M. Hammer (SINTEF) Introduction quality issues
More 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 informationMolecular dynamics simulation of the liquid vapor interface: The Lennard-Jones fluid
Molecular dynamics simulation of the liquid vapor interface: The Lennard-Jones fluid Matthias Mecke and Jochen Winkelmann a) Institut für Physikalische Chemie, Universität Halle-Wittenberg, Geusaer Str.,
More informationMonte Carlo Calculations of Effective Surface Tension for Small Clusters
Monte Carlo Calculations of Effective Surface Tension for Small Clusters Barbara N. Hale Physics Department and Cloud and Aerosol Science Laboratory, University of Missouri- Rolla, Rolla, MO 65401, USA
More informationChapter 11. Intermolecular Forces and Liquids & Solids
Chapter 11 Intermolecular Forces and Liquids & Solids The Kinetic Molecular Theory of Liquids & Solids Gases vs. Liquids & Solids difference is distance between molecules Liquids Molecules close together;
More informationDirect determination of the Tolman length from the bulk pressures of liquid drops via molecular dynamics simulations
Direct determination of the Tolman length from the bulk pressures of liquid drops via molecular dynamics simulations Alan E. van Giessen, and Edgar M. Blokhuis Citation: The Journal of Chemical Physics
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 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 informationCommunication: Slab thickness dependence of the surface tension: Toward a criterion of liquid sheets stability
THE JOURNAL OF CHEMICAL PHYSICS 141, 081103 (2014) Communication: Slab thickness dependence of the surface tension: Toward a criterion of liquid sheets stability G. Filippini, 1 E. Bourasseau, 1,a) A.
More informationk θ (θ θ 0 ) 2 angles r i j r i j
1 Force fields 1.1 Introduction The term force field is slightly misleading, since it refers to the parameters of the potential used to calculate the forces (via gradient) in molecular dynamics simulations.
More informationMonte 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 informationMolecular modeling of hydrogen bonding fluids: vapor-liquid coexistence and interfacial properties
12 th HLRS Results and Review Workshop Molecular modeling of hydrogen bonding fluids: vapor-liquid coexistence and interfacial properties High Performance Computing Center Stuttgart (HLRS), October 8,
More informationMolecular Dynamics Simulation of Methanol-Water Mixture
Molecular Dynamics Simulation of Methanol-Water Mixture Palazzo Mancini, Mara Cantoni University of Urbino Carlo Bo Abstract In this study some properties of the methanol-water mixture such as diffusivity,
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 informationPhase Equilibria of binary mixtures by Molecular Simulation and PR-EOS: Methane + Xenon and Xenon + Ethane
International Journal of ChemTech Research CODEN( USA): IJCRGG ISSN : 0974-4290 Vol.5, No.6, pp 2975-2979, Oct-Dec 2013 Phase Equilibria of binary mixtures by Molecular Simulation and PR-EOS: Methane +
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 informationMolecular simulation of adsorption from dilute solutions
Vol. 52 No. 3/2005 685 689 on-line at: www.actabp.pl Molecular simulation of adsorption from dilute solutions Werner Billes Rupert Tscheliessnig and Johann Fischer Institut für Verfahrens- und Energietechnik
More informationChapter 12. Liquids: Condensation, Evaporation, and Dynamic Equilibrium
Chapter 12 Liquids: Condensation, Evaporation, and Dynamic Equilibrium Chapter Map Condensation (Gas to Liquid) Evaporation For a particle to escape from the surface of the liquid, it must meet the following
More informationGibbs ensemble simulation of phase equilibrium in the hard core two-yukawa fluid model for the Lennard-Jones fluid
MOLECULAR PHYSICS, 1989, VOL. 68, No. 3, 629-635 Gibbs ensemble simulation of phase equilibrium in the hard core two-yukawa fluid model for the Lennard-Jones fluid by E. N. RUDISILL and P. T. CUMMINGS
More informationThe gas-liquid surface tension of argon: A reconciliation between experiment and simulation
THE JOURNAL OF CHEMICAL PHYSICS 140, 44710 (014) The gas-liquid surface tension of argon: A reconciliation between experiment and simulation Florent Goujon, 1,a) Patrice Malfreyt, 1 and Dominic J. Tildesley
More informationChapter 8: Physical Equilibria
Chapter 8: Physical Equilibria Our first foray into equilibria is to examine phenomena associated with two phases of matter achieving equilibrium in which the free energy in each phase is the same and
More informationAtomistic 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 informationMolecular 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 informationComparison of different mixing rules for prediction of density and residual internal energy of binary and ternary Lennard Jones mixtures
Fluid Phase Equilibria 178 (2001) 87 95 Comparison of different mixing rules for prediction of density and residual internal energy of binary and ternary Lennard Jones mixtures Jian Chen a,, Jian-Guo Mi
More informationMolecular dynamics study on the evaporation part of the kinetic boundary condition at the interface between water and water vapor
Molecular dynamics study on the evaporation part of the kinetic boundary condition at the interface between water and water vapor Tatsuya Ishiyama, Takeru Yano and Shigeo Fujikawa Division of Mechanical
More informationMolecular Interactions F14NMI. Lecture 4: worked answers to practice questions
Molecular Interactions F14NMI Lecture 4: worked answers to practice questions http://comp.chem.nottingham.ac.uk/teaching/f14nmi jonathan.hirst@nottingham.ac.uk (1) (a) Describe the Monte Carlo algorithm
More informationName: Class: Date: SHORT ANSWER Answer the following questions in the space provided.
CHAPTER 10 REVIEW States of Matter SECTION 1 SHORT ANSWER Answer the following questions in the space provided. 1. Identify whether the descriptions below describe an ideal gas or a real gas. a. The gas
More informationOn the Calculation of the Chemical Potential. Using the Particle Deletion Scheme
On the Calculation of the Chemical Potential Using the Particle Deletion Scheme Georgios C. Boulougouris,2, Ioannis G. Economou and Doros. Theodorou,3,* Molecular Modelling of Materials Laboratory, Institute
More informationUnusual Entropy of Adsorbed Methane on Zeolite Templated Carbon. Supporting Information. Part 2: Statistical Mechanical Model
Unusual Entropy of Adsorbed Methane on Zeolite Templated Carbon Supporting Information Part 2: Statistical Mechanical Model Nicholas P. Stadie*, Maxwell Murialdo, Channing C. Ahn, and Brent Fultz W. M.
More informationDirect calculation of interfacial tensions from computer simulation: Results for freely jointed tangent hard sphere chains
Direct calculation of interfacial tensions from computer simulation: Results for freely jointed tangent hard sphere chains Luis G. MacDowell Departamento de Química Física, Facultad de Ciencias Químicas,
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 informationChapter 13 States of Matter Forces of Attraction 13.3 Liquids and Solids 13.4 Phase Changes
Chapter 13 States of Matter 13.2 Forces of Attraction 13.3 Liquids and Solids 13.4 Phase Changes I. Forces of Attraction (13.2) Intramolecular forces? (forces within) Covalent Bonds, Ionic Bonds, and metallic
More informationContents. 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 informationThe Second Virial Coefficient & van der Waals Equation
V.C The Second Virial Coefficient & van der Waals Equation Let us study the second virial coefficient B, for a typical gas using eq.v.33). As discussed before, the two-body potential is characterized by
More informationIntermolecular Model Potentials and Virial Coefficients from Acoustic Data
JASEM ISSN 1119-8362 All rights reserved Full-text Available Online at https://www.ajol.info/index.php/jasem http://www.bioline.org.br/ja J. Appl. Sci. Environ. Manage. Vol.22 (2) 246-251. February 2018
More informationPhase transitions of quadrupolar fluids
Phase transitions of quadrupolar fluids Seamus F. O Shea Department of Chemistry, University of Lethbridge, Lethbridge, Alberta, Canada, T1K 3M4 Girija S. Dubey Brookhaven National Laboratory, Upton, New
More informationChapter 11. Intermolecular forces. Chapter 11 1
Chapter 11 Intermolecular Attractions and the Properties of Liquids and Solids 1 2 Intermolecular forces Forces of attraction between molecules Directly dependent on the distance between the molecules
More informationChapter Intermolecular attractions
Chapter 11 11.2 Intermolecular attractions Intermolecular Attractions and the Properties of Liquids and Solids Intermolecular forces control the physical properties of the substance. Intramolecular forces
More informationChemical Potential of Benzene Fluid from Monte Carlo Simulation with Anisotropic United Atom Model
Chemical Potential of Benzene Fluid from Monte Carlo Simulation with Anisotropic United Atom Model Mahfuzh Huda, 1 Siti Mariyah Ulfa, 1 Lukman Hakim 1 * 1 Department of Chemistry, Faculty of Mathematic
More informationEquations of State. Equations of State (EoS)
Equations of State (EoS) Equations of State From molecular considerations, identify which intermolecular interactions are significant (including estimating relative strengths of dipole moments, polarizability,
More informationExperimental measurements and numerical modelling of interfacial tension in water-methane systems
Proceedings of the 4th WSEAS International Conference on Heat and Mass Transfer, Gold Coast, Queensland, Australia, January 17-19, 07 Experimental measurements and numerical modelling of interfacial tension
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 information2D Clausius Clapeyron equation
2D Clausius Clapeyron equation 1/14 Aim: Verify the Clausius Clapeyron equation by simulations of a 2D model of matter Model: 8-4 type potential ( Lennard-Jones in 2D) u(r) = 1 r 8 1 r 4 Hard attractive
More informationIdeal Gas Behavior. NC State University
Chemistry 331 Lecture 6 Ideal Gas Behavior NC State University Macroscopic variables P, T Pressure is a force per unit area (P= F/A) The force arises from the change in momentum as particles hit an object
More informationUniversity of Groningen. Characterization of oil/water interfaces van Buuren, Aldert Roelf
University of Groningen Characterization of oil/water interfaces van Buuren, Aldert Roelf IMPORTANT NOTE: You are advised to consult the publisher's version (publisher's PDF) if you wish to cite from it.
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 informationCARBON 2004 Providence, Rhode Island. Adsorption of Flexible n-butane and n-hexane on Graphitized Thermal Carbon Black and in Slit Pores
CARBON Providence, Rhode Island Adsorption of Flexible n-butane and n-hexane on Graphitized Thermal Carbon Black and in Slit Pores D. D. Do* and H. D. Do, University of Queensland, St. Lucia, Qld 7, Australia
More informationMODELING OF PHASE EQUILIBRIA FOR BINARY AND TERNARY MIXTURES OF CARBON DIOXIDE, HYDROGEN AND METHANOL
MODELING OF PHASE EQUILIBRIA FOR BINARY AND TERNARY MIXTURES OF CARBON DIOXIDE, HYDROGEN AND METHANOL Neil R. Foster *, Keivan Bezanehtak, Fariba Dehghani School of Chemical Engineering and Industrial
More informationMolecular Modeling of Matter
Molecular Modeling of Matter Keith E. Gubbins Lecture 1: Introduction to Statistical Mechanics and Molecular Simulation Common Assumptions Can treat kinetic energy of molecular motion and potential energy
More informationPORE SIZE DISTRIBUTION OF CARBON WITH DIFFERENT PROBE MOLECULES
PORE SIZE DISTRIBUTION OF CARBON WITH DIFFERENT PROBE MOLECULES Atichat Wongkoblap*, Worapot Intomya, Warangkhana Somrup, Sorod Charoensuk, Supunnee Junpirom and Chaiyot Tangsathitkulchai School of Chemical
More informationJoy Billups Woller A THESIS. Presented to the Faculty of. The Graduate College in the University of Nebraska
DETERMINATION OF THE THERMODYNAMIC PROPERTIES OF FLUIDS BY GIBBS ENSEMBLE MONTE CARLO COMPUTER SIMULATIONS by Joy Billups Woller A THESIS Presented to the Faculty of The Graduate College in the University
More informationUsing Fundamental Measure Theory to Treat the Correlation Function of the Inhomogeneous Hard-Sphere Fluid
Using Fundamental Measure Theory to Treat the Correlation Function of the Inhomogeneous Hard-Sphere Fluid Jeff B. Schulte, Patrick A. Kreitzberg, Chris V. Haglund, and David Roundy Department of Physics,
More informationIon-Gated Gas Separation through Porous Graphene
Online Supporting Information for: Ion-Gated Gas Separation through Porous Graphene Ziqi Tian, Shannon M. Mahurin, Sheng Dai,*,, and De-en Jiang *, Department of Chemistry, University of California, Riverside,
More informationDensity dependence of dielectric permittivity of water and estimation of the electric field for the breakdown inception
Journal of Physics: Conference Series PAPER OPEN ACCESS Density dependence of dielectric permittivity of water and estimation of the electric field for the breakdown inception To cite this article: D I
More informationChemical Potential, Helmholtz Free Energy and Entropy of Argon with Kinetic Monte Carlo Simulation
Chemical Potential, Helmholtz Free Energy and Entropy of Argon with Kinetic Monte Carlo Simulation Chunyan Fan a, D. D. Do a *, D. icholson a and E. Ustinov b a School of Chemical Engineering, University
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 information424 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 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 informationGases and the Virial Expansion
Gases and the irial Expansion February 7, 3 First task is to examine what ensemble theory tells us about simple systems via the thermodynamic connection Calculate thermodynamic quantities: average energy,
More informationIdeal Gas Laws Empirical Gas Laws The Mole Equations of State Dalton's Law The Mole Fraction Extensive and Intensive Variables Graham's Law of
Ideal Gas Laws Empirical Gas Laws The Mole Equations of State Dalton's Law The Mole Fraction Extensive and Intensive Variables Graham's Law of Effusion The Maxwell-Boltzmann Distribution A Digression on
More informationTHE DETAILED BALANCE ENERGY-SCALED DISPLACEMENT MONTE CARLO ALGORITHM
Molecular Simulation, 1987, Vol. 1, pp. 87-93 c Gordon and Breach Science Publishers S.A. THE DETAILED BALANCE ENERGY-SCALED DISPLACEMENT MONTE CARLO ALGORITHM M. MEZEI Department of Chemistry, Hunter
More informationof Nebraska - Lincoln
University of Nebraska - Lincoln DigitalCommons@University of Nebraska - Lincoln Xiao Cheng Zeng Publications Published Research - Department of Chemistry 10-1-2006 Homogeneous nucleation at high supersaturation
More informationAs we ended the lectures on gases, we were introduced to an idea that serves as foundation for the material in this lecture:
LECTURE 16. INTRODUCTION TO INTERMOLECULAR FORCES As we ended the lectures on gases, we were introduced to an idea that serves as foundation for the material in this lecture: As we were introduced to ideal
More information2D Clausius Clapeyron equation
2D Clausius Clapeyron equation 1/14 Aim: Verify the Clausius Clapeyron equation by simulations of a 2D model of matter Model: 8-4 type potential ( Lennard-Jones in 2D) u(r) = 1 r 8 1 r 4 Hard attractive
More informationThe Effect of Model Internal Flexibility Upon NEMD Simulations of Viscosity
Draft: September 29, 1999 The Effect of Model Internal Flexibility Upon NEMD Simulations of Viscosity N. G. Fuller 1 and R. L. Rowley 1,2 Abstract The influence of model flexibility upon simulated viscosity
More informationMolecular dynamics investigation of thickness effect on liquid films
JOURNAL OF CHEMICAL PHYSICS VOLUME 113, NUMBER 14 8 OCTOBER 2000 Molecular dynamics investigation of thicness effect on liquid films Jian-Gang Weng, Seungho Par, a) Jennifer R. Lues, and Chang-Lin Tien
More informationElectrochemical Properties of Materials for Electrical Energy Storage Applications
Electrochemical Properties of Materials for Electrical Energy Storage Applications Lecture Note 3 October 11, 2013 Kwang Kim Yonsei Univ., KOREA kbkim@yonsei.ac.kr 39 Y 88.91 8 O 16.00 7 N 14.01 34 Se
More informationPressure Volume Temperature Relationship of Pure Fluids
Pressure Volume Temperature Relationship of Pure Fluids Volumetric data of substances are needed to calculate the thermodynamic properties such as internal energy and work, from which the heat requirements
More informationMayer-sampling Monte Carlo calculations of methanol virial. Katherine R. S. Shaul, Andrew J. Schultz, David A. Kofke*
Mayer-sampling Monte Carlo calculations of methanol virial coefficients Katherine R. S. Shaul, Andrew J. Schultz, David A. Kofke* Department of Chemical & Biological Engineering, University at Buffalo,
More informationV.C The Second Virial Coefficient & van der Waals Equation
V.C The Second Virial Coefficient & van der Waals Equation Let us study the second virial coefficient B, for a typical gas using eq.(v.33). As discussed before, the two-body potential is characterized
More informationChapter 10 Practice. Name: Multiple Choice Identify the choice that best completes the statement or answers the question.
Name: Score: 0 / 18 points (0%) [3 open ended questions not graded] Chapter 10 Practice Multiple Choice Identify the choice that best completes the statement or answers the question. 1. A sample of gas
More informationModule 5: "Adsoption" Lecture 25: The Lecture Contains: Definition. Applications. How does Adsorption occur? Physisorption Chemisorption.
The Lecture Contains: Definition Applications How does Adsorption occur? Physisorption Chemisorption Energetics Adsorption Isotherms Different Adsorption Isotherms Langmuir Adsorption Isotherm file:///e
More informationEffect of polarizability of halide anions on the ionic salvation in water clusters
University of Nebraska - Lincoln DigitalCommons@University of Nebraska - Lincoln Xiao Cheng Zeng Publications Published Research - Department of Chemistry 9-22-2003 Effect of polarizability of halide anions
More informationLecture 1: Physical Equilibria The Temperature Dependence of Vapor Pressure
Lecture 1: Physical Equilibria The Temperature Dependence of Vapor Pressure Our first foray into equilibria is to examine phenomena associated with two phases of matter achieving equilibrium in which the
More informationInterfacial Thermodynamics of Water and Six Other Liquid Solvents
pubs.acs.org/jpcb Interfacial Thermodynamics of Water and Six Other Liquid Solvents Tod A. Pascal*, and William A. Goddard, III* Materials and Process Simulation Center, California Institute of Technology,
More informationSurface tension of the two center Lennard-Jones plus quadrupole model fluid
Surface tension of the two center Lennard-Jones plus quadrupole model fluid Stephan Werth, Martin Horsch 1, Hans Hasse Laboratory of Engineering Thermodynamics, Department of Mechanical and Process Engineering,
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 informationChapter #16 Liquids and Solids
Chapter #16 Liquids and Solids 16.1 Intermolecular Forces 16.2 The Liquid State 16.3 An Introduction to Structures and Types of Solids 16.4 Structure and Bonding of Metals 16.5 Carbon and Silicon: Network
More informationStatistical 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 informationDetermination of Phase Equilibria and the Critical Point Using Two-Phase Molecular Dynamics Simulations with Monte Carlo Sampling
Brigham Young University BYU ScholarsArchive All Theses and Dissertations 2012-06-15 Determination of Phase Equilibria and the Critical Point Using Two-Phase Molecular Dynamics Simulations with Monte Carlo
More informationModule 4: "Surface Thermodynamics" Lecture 21: "" The Lecture Contains: Effect of surfactant on interfacial tension. Objectives_template
The Lecture Contains: Effect of surfactant on interfacial tension file:///e /courses/colloid_interface_science/lecture21/21_1.htm[6/16/2012 1:10:36 PM] Surface Thermodynamics: Roles of Surfactants and
More informationImproved association in a classical density functional theory for water. Abstract
Improved association in a classical density functional theory for water Eric J. Krebs, Jeff B. Schulte, and David Roundy Department of Physics, Oregon State University, Corvallis, OR 97331, USA Abstract
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 information