An Extended van der Waals Equation of State Based on Molecular Dynamics Simulation
|
|
- Anna Russell
- 6 years ago
- Views:
Transcription
1 J. Comput. Chem. Jpn., Vol. 8, o. 3, pp (9) c 9 Society of Computer Chemistry, Japan An Extended van der Waals Equation of State Based on Molecular Dynamics Simulation Yosuke KATAOKA* and Yuri YAMADA Department of Chemical Science and Technology, Faculty of Bioscience and Applied Chemistry, Hosei University 3-7- Kajino-cho, Koganei, Tokyo , Japan * kataoka@hosei.ac.jp (Received: February 9, 9; Accepted for publication: July 7, 9; Advance publication: August 4, 9) Molecular dynamics simulations on the Lennard-Jones system are performed to obtain the pvt and UVT relations. An extended van der Waals equation of state (EOS) is derived by statistical mechanics on the perturbation approximation. A hard sphere system is used as the reference system. The attraction energy term in the canonical ensemble partition function is extended by a cluster expansion. The new EOS includes three parameters, two of which are the interaction parameters in the Lennard-Jones interaction. The last parameter is the effective volume of the hard sphere system. The extended van der Waals EOS reproduces the pvt and UVT relations, at least qualitatively, whereas the original van der Waals EOS can explain only the pvt relation. Keywords: Extended van der Waals equation of state, Lennard-Jones system, Molecular dynamics simulation 1 Introduction The van der Waals EOS is well known in physical chemistry because it has a very simple form and explains the characteristics of a real gas [1]. The following pvt relation also explains the vapor-liquid critical point: p(,v,t ) = kt ( ) V b a. (1) V here a and b are constants that characterize the molecules. The variables, V, and T are the number of molecules, the volume, and the temperature of the system. The Boltzmann constant is denoted as k. The vapor-liquid transition point is estimated by the Maxwell construction in the van der Waals loop region in the pvt relation, where the unstable state is also found [1]. This loop appears because van der Waals theory assumes a single phase. Molecular simulation also treats the phase as a single phase in a conventional method. This loop is found in molecular simulations with 1-1, molecules [ 6]. The loop appears in molecular simulation because the molecular system has uniform structure on the macroscopic scale under the periodic boundary condition. In this sense, the van der Waals EOS is important in the understanding of the molecular simulation results. On the other hand, the UVT relation in the van der Waals EOS is too simple [6]. The internal energy U is described as follows in the van der Waals EOS [7]: U(,V,T ) = 3 a kt V. () The heat capacity under the constant volume C V of real liquid depends on temperature, as confirmed by macroscopic and molecular observations [1]. The present study proposes an extended form of the van der Waals EOS to explain the pvt and UVT relations of the gas and liquid phase of the spherical molecules by a simple treatment of the attraction part in the canonical ensemble partition function. The formulation follows the perturbation theory with the hard sphere system as a reference system [7]. Although the original van der Waals EOS includes only the averaged value of the attraction energy, the extended van der Waals EOS represents the attraction part by a cluster expansion in an intuitive manner, whereby the Lennard-Jones interaction is assumed as typical in the spherical molecules. The extended van der Waals EOS for the Lennard- Jones system has only three parameters, two of which are the Lennard-Jones parameters. The remaining parameter is the effective volume of the hard sphere of the reference system. The pvt and UVT relations will be compared with those of the molecular dynamics [MD] simulation on the Lennard-Jones system. DOI: 1.477/jccj.H1 97
2 The present paper follows the statistical mechanical theory on the original van der Waals EOS [7], which is then extended to explain, at least qualitatively, the UVT relation in a broad region in the temperature-density domain. Finally, the theoretical and simulation results are compared. Theory of van der Waals EOS The theory of van der Waals EOS [7] is reviewed in this section and will be extended in the following section. A molecular system with spherical molecules of mass m at temperature T and volume V is considered. The canonical partition function Q(,V,T) is given as follows: Q(,V,T ) = 1 ( ) πmkt 3/! h Z, (3) Z Z = e U(r 1,,r )/kt dr 1 dr. (4) V Here, the potential energy of the system is written as U(r 1, r,, r ). The quantity Z is called the configuration integral. In the case of a perfect gas, Z is given as follows, because no interaction is assumed: Z = V. (5) The thermodynamic quantities are written by the partition function as follows: A(,V,T ) = kt lnq(,v,t ), (6) ( ) lnq p = kt. (7) V,T According to this formula, in the case of a perfect gas, the pressure is given as follows: p = kt V. (8) The interaction energy is approximated by the van der Waals approximation. The statistical mechanical theory on van der Waals EOS adopts the perturbation approximation. The total potential energy is divided into two parts: U = U () +U(1), (9) where U () (1) and U are the potential energy of the reference system and that of the perturbed system. The quantity β is the inverse of temperature: β = 1/kT. The configuration integral is as follows: Z Z ( Z = e β U () ) +U(1) dr 1 dr Z Z Z = e βu() ( R β R e dr 1 dr U () ) +U(1) dr 1 dr R. (1) R e βu () dr 1 dr This equation can be rewritten as follows: ( ) Z = Z () exp βu (1). (11) The reference system is the hard sphere system, and its pair potential function u HS (r) is as follows: u HS (r) =, r σ HS u HS (r) =, r > σ HS. (1) The pair correlation function g(r) of the hard sphere system is approximated as follows: g HS (r) =, r σ HS g HS (r) = 1, r > σ HS. (13) When ( βu (1) ) is sufficiently small, the term exp βu (1) in the configuration integral Z can be rewritten as follows: ( exp ) 1 β The average βu (1) U (1) = exp U (1) βu (1) (14) is estimated by the above pair correlation function as follows: U (1) πρ V Z u (1) (r)g HS (r)r dr σ HS a = π = aρ, Z Here u (1) (r) is the pair potential: U (1) σ HS u (1) (r)r dr (15) = u (1) (r). (16) i< j In this manner, the configuration integral can be written as follows: Z = Z () eβaρ. (17) In addition, Z () Z () is approximated as follows: = (V b), b = πσ3 HS. (18) 3 Using this configuration integral, one can write the pressure as follows: ( ) lnq p = kt V,T p = kt (19) 1 bρ aρ. 98 J. Comput. Chem. Jpn., Vol. 8, o. 3 (9)
3 The internal energy U is derived from the partition function as follows: Q(,V,T ) = 1 ( ) πmkt 3/! h (V b) e a /V kt, ( ) lnq(,v,t ) U = kt, T,V () U = 3 kt aρ. (1) 3 Extension of van der Waals EOS In this section, van der Waals EOS is extended by a cluster expansion of the attractive term Z (1) in the canonical partition function in the perturbation scheme described in the previous section. There may be monomers, dimers, trimers, etc., in the system according to the temperature and density. The largest cluster in the system is assumed to be an M-mer. The structure and energy of the J-mer cluster is simplified as follows. The J-mer cluster has (J 1) bonds that start from the central molecule being considered. Each bond has energy ε. The energy of the J-mer cluster is then 1 (J 1)ε. The factor 1/ appears because the present theory is a one-body approximation, and this factor is necessary in order to obtain a reasonable total energy for the system. Therefore, in the present theory, the J-mer cluster has the Boltzmann factor exp ( 1 β(j 1)ε ) according to the temperature of the system. The geometrical probability of its being found in the system is proportional to the pre-exponential factor ( ν ) (J 1)/3 ν. Here, ν is the volume per molecule, and ν is that in the case of the closest packing. The exponent (J 1)/3 is assumed because the finding probability of the J-mer cluster is proportional to the ratio of the length of the closest packing structure to that of the characteristic length of the volume per molecule of the system. Based on this assumption the partition function is written as follows: Z = Z () Z(1) Z () = (V b) () ( Z (1) = ν ) ( ) ] (J 1)/3 1 exp ν β(j 1)ε. M J=1 [ M J=1 ext, we introduce the molecular partition function of the attractive term q a as follows: [ ] q a = Z (1) 1/ ) ( ) (J 1)/3 (3) 1 = exp β(j 1)ε. ( ν ν This molecular partition function is compared with that of the original EOS, e β aρ. For this molecular partition function, the attractive term of the internal energy U a is obtained as follows: ( ) U a = kt lnqa T U a = 1 q a M J=1,V 1 ( (J 1)ε ν ν exp, ) (J 1)/3 ( 1 β(j 1)ε ). (4) The contribution of the attractive interaction to pressure is as follows: ( ) lnqa p a = kt V p a = kt 1 q a [ M J=1,T, (J 1) ( ν ) (J 1)/3 1 3 ν ν ( ) ] 1 exp β(j 1)ε. (5) The total internal energy and pressure are as follows: U = 3 kt + U a, p = kt V b + p a. (6) The Helmholtz free energy and entropy are calculated by the standard method from the partition function. If we express the molecular interaction in the Lennard-Jones form, the parameters in the attractive part are given by the Lennard-Jones parameters. The Lennard-Jones potential is written as follows: u(r) = 4ε LJ [ (σ r ) 1 ( σ r ) 6 ]. (7) The constant ε in Eq. () is equal to ε LJ in Eq. (7). The closest packing structure in the Lennard-Jones system has the following volume: ν = σ3. (8) This value is used in Eq. (). The last parameter in the present model is the effective volume of the hard sphere b in Eq. (). We adjust this parameter to fit the calculated pressure and internal energy with the observed pressure and internal energy (macroscopic or molecular observations). DOI: 1.477/jccj.H1 99
4 4 Molecular Dynamics Simulation VT MD simulations [8] on fluid argon were performed at several temperatures to obtain the average pressure p and internal energy U as a function of volume V. The Lennard-Jones potential is adopted as shown in Table 1. The time step is 1 fs, and the standard MD length is 1,, steps. The velocity scaling method was used to control temperature. The number of molecules in the cell is 1,. A periodic boundary condition is assumed on the cubic basic cell [8]. The assumed initial configuration is the simple cubic lattice, in which only the density varies according to the volume. The Materials Explorer (v4 and v5) program was used for the present simulation [9]. The obtained pvt and UVT relations are shown in Figure 1 and Figure. The abscissa shows the density Table 1. Lennard-Jones potential parameters for argon [9]. ε/(kcal/mol) (ε/k)/k σ/m E-1 ρ = /V. The obtained pressure qualitatively resembles the van der Waals EOS (see Eq. (19)). The pressure is given in units of ε/σ 3 as well as in units of atm [1]. On the other hand, the calculated internal energy depends on temperature and density in a different manner than the van der Waals EOS shown in Eq. (1). This difference is discussed in the next section. 5 pvt and UVT relations in the extended van der Waals EOS The adjustable parameter b is chosen by trial and error as σ 3 through a comparison of the extended EOS with the results of MD simulation. The fitted pvt and UVT relations are described in Figure 3 and Figure 4 on the extended van der Waals EOS. The pvt relation in Figure 3 is globally very similar to that in Figure 1. On the whole, the internal energy in Figure 4 also resembles that of the Lennard-Jones system shown in Figure. A more direct comparison is shown in Figure 5 and Figure 6. The pvt and UVT relations in Figure 5 and Figure 6 mean that the extended EOS can reproduce, at least qualitatively, the MD simulation results. Figure 1. The pvt relation calculated by MD simulation on fluid argon using the Lennard-Jones model. The ordinate shows pressure p given in Lennard-Jones units ε/σ 3 (Table 1) and in atm [1]. The abscissa shows the number density /V and the mass density ρ in units of g/cm 3. The constant temperature curves are shown at T = 1, 5, 1, 175, 5, and 5 K. Figure. The UVT relation calculated by MD simulation on fluid argon using the Lennard-Jones model as a function of density. The ordinate shows the internal energy per molecule U/ in Lennard-Jones units ε (Table 1) and in units of kj/mol. The abscissa shows the number density /V and the mass density ρ in units of g/cm 3. The constant temperature curves are shown at T = 1, 5, 1, 175, 5, and 5 K. 1 J. Comput. Chem. Jpn., Vol. 8, o. 3 (9)
5 Figure 3. The pvt relation calculated by the extended van der Waals EOS. The ordinate is pressure p in Lennard-Jones units ε/σ 3 (Table 1). The abscissa is the number density /V. The constant temperature curves are shown at T = 1, 5, 1, 175, 5, 75, and 5 K. Figure 5. The pvt relation calculated by the extended van der Waals EOS compared with that obtained by MD simulation. The ordinate shows pressure p in Lennard- Jones units ε/σ 3 (Table 1). The abscissa shows the number density /V. The constant temperature curves are shown at T = 5,, and 5 K. Figure 4. The UVT relation calculated by the extended van der Waals EOS as a function of density. The ordinate shows the internal energy per molecule U/ in Lennard- Jones units ε (Table 1). The abscissa shows the number density /V. The constant temperature curves are shown at T = 1, 5, 1, 175, 5, 75, and 5 K. Figure 6. The UVT relation calculated by the extended van der Waals EOS compared with that obtained by MD simulation. The ordinate shows the internal energy per molecule in Lennard-Jones units ε (Table 1). The abscissa shows the number density /V. The constant temperature curves are shown at T = 5,, and 5 K. DOI: 1.477/jccj.H1 11
6 6 Extended EOS compared with MD simulation The calculated results are compared with the MD simulations in Figure 5 and Figure 6. In these figures, the present EOS provides reasonable results for pressure and internal energy, although the agreement is only semiquantitative. Since we have only three parameters in the present EOS, the agreement is satisfactory. Although the function is not so simple as the original van der Waals EOS, the physical meanings of the three parameters are clear in the present case. The construction rule is also unambiguous in the canonical partition function. Table. The critical point determined by the extended van der Waals EOS is compared with that in the fitted EOS of MD simulations on Lenard-Jones system. T c /(ε/k) p c /(ε/σ 3 ) ρ c /σ -3 icolas et al [11] Kolafa & ezbeda [1] Present Study T c /K p c /atm V c /(cm 3 /mol) icolas et al [11] Kolafa & ezbeda [1] Present Study The simulation results are usually fitted by EOS with several fitting parameters [11, 1]. The critical point is determined by such EOS [11, 1]. The critical point determined by the present EOS is compared with that shown in Table. The obtained result does not differ greatly from the MD simulation results even with only three parameters in the present extended EOS. As shown in Figure 7, the extended EOS pressure is different from the MD results at low temperatures. The effect of attractive interaction is not well incorporated in this case. This is due to the simplified cluster structure, which is determined by the short-range term in the molecular interaction. The original van der Waals EOS also yields poor results at low temperatures, as shown in Figure 7. This curve shows that the attractive force is too strong in this EOS. The van der Waals coefficients used in Figure 7 and Figure 8 are shown in Table 3 [1]. A comparison of internal energy values is shown in Figure 8 at low temperatures. Table 3. The van der Waals coefficients for fluid argon [1]. a/(atm L /mol ) b/(l/mol) E- Figure 7. Pressures are compared for the extended van der Waals EOS, the original van der Waals EOS, and the MD results on the Lennard-Jones system as a model of argon at T = 1 K =.8 ε/k. The abscissa shows the number density /V. Figure 8. Average potential energies per molecule U e / = U/ 1.5kT are compared for the extended van der Waals EOS, the original van der Waals EOS, and the MD results on the Lennard-Jones system as a model of argon at T = 1 K =.8 ε/k. The abscissa shows the number density /V. 1 J. Comput. Chem. Jpn., Vol. 8, o. 3 (9)
7 The value 6ε shown in Figure 8 of the present study is obtained from the 1 nearest neighbors in the FCC lattice. With respect to the internal energy, the extended van der Waals EOS is better than the original van der Waals EOS, as shown in Figure 8. The shortcomings in the present study (as shown in Figure 7 and Figure 8) may be overcome by the inclusion of the long-range part of the attractive interaction. 7 Conclusion An extended van der Waals EOS is obtained by a statistical mechanical treatment on the Lennard-Jones system. The new EOS reproduces the pvt and UVT relations calculated by molecular dynamics simulations, at least qualitatively. As the effective Lennard-Jones parameters are known for many molecular systems, the new EOS is useful to estimate the pvt and UVT relations on such systems. The authors would like to thank the Research Center for Computing and Multimedia Studies, at Hosei University for use of computers. References [1] P. W. Atkins, Physical Chemistry, 6th edition, Oxford Univ. Press (1998). [] R. Yamamoto, H. Tanaka, K. akanishi, and X. C. Zeng, Chem. Phys. Lett., 31, 41 (1994). [3] R. Yamamoto and K. akanishi, Phys. Rev. B, 49, (1994). [4] R. Yamamoto, O. Kitao, and K. akanishi, Mol. Phys., 84, 757 (1995). [5] R. Yamamoto and K. akanishi, Phys. Rev. B, 51, 715 (1995). [6] R. Yamamoto and K. akanishi, Mol. Simulation, 16, 119 (1996). [7] D. A. McQuarrie, Statistical Mechanics, Harper Collins (1976). [8] M. P. Allen and D. J. Tildesley, Computer Simulation of Liquids, Oxford Univ. Press, Oxford (1987). [9] [1] The relation between the units of pressure is given as follows: 1 atm = 11,35 Pa. [11] J. J. icolas, K. E. Gubbins, W. B. Streett, and D. J. Tildesley, Mol. Phys., 37, 149 (1979). [1] J. Kolafa and I. ezbeda, Fluid Phase Equil., 1, 1 (1994). DOI: 1.477/jccj.H1 13
8 14 J. Comput. Chem. Jpn., Vol. 8, o. 3 (9)
Thermodynamics of Three-phase Equilibrium in Lennard Jones System with a Simplified Equation of State
23 Bulletin of Research Center for Computing and Multimedia Studies, Hosei University, 28 (2014) Thermodynamics of Three-phase Equilibrium in Lennard Jones System with a Simplified Equation of State Yosuke
More informationPhase Equilibrium Calculations by Equation of State v2
Bulletin of Research Center for Comuting and Multimedia Studies, Hosei University, 7 (3) Published online (htt://hdl.handle.net/4/89) Phase Equilibrium Calculations by Equation of State v Yosuke KATAOKA
More informationImperfect Gases. NC State University
Chemistry 431 Lecture 3 Imperfect Gases NC State University The Compression Factor One way to represent the relationship between ideal and real gases is to plot the deviation from ideality as the gas is
More informationChapter 2 Experimental sources of intermolecular potentials
Chapter 2 Experimental sources of intermolecular potentials 2.1 Overview thermodynamical properties: heat of vaporization (Trouton s rule) crystal structures ionic crystals rare gas solids physico-chemical
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 informationG : Statistical Mechanics
G5.651: Statistical Mechanics Notes for Lecture 8 I. DISTRIBUTION FUNCTIONS AND PERTURBATION THEORY A. General formulation Recall the expression for the configurational partition function: Z N = dr 1 dr
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 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 informationMOLECULAR DYNAMICS STUDY OF THE NUCLEATION OF BUBBLE
CAV2:sessionA.5 MOLECULAR DYNAMICS STUDY OF THE NUCLEATION OF BUBBLE Takashi Tokumasu, Kenjiro Kamijo, Mamoru Oike and Yoichiro Matsumoto 2 Tohoku University, Sendai, Miyagi, 98-8577, Japan 2 The University
More informationSIMULATIONAL ANALYSIS OF GLASSES PREPARED VIA DIFFERENT INTERATOMIC POTENTIALS
Journal of Optoelectronics and Advanced Materials Vol. 7, No. 4, August 2005, p. 1915-1922 SIMULATIONAL ANALYSIS OF GLASSES PREPARED VIA DIFFERENT INTERATOMIC POTENTIALS Y. Sano, J. K ga, F. Yonezawa Department
More informationIntermolecular Potentials and The Second Virial Coefficient
Intermolecular Potentials and The Second Virial Coefficient Patrick L. Holt Department of Chemistry & Physics Bellarmine University Louisville, Kentucky 40205 pholt@bellarmine.edu Copyright 2004 by the
More informationEquation of State for Free Energy of Homogeneous Nucleation in Supersaturated Lennard-Jones Vapor Phase Derived by Monte Carlo Simulations
Equation of State for Free Energy of Homogeneous Nucleation in Supersaturated Lennard-Jones Vapor Phase Derived by Monte Carlo Simulations YAMADA Yuri March, 2003 Hosei University ABSTRACT Monte Carlo
More information5.62 Physical Chemistry II Spring 2008
MIT OpenCourseWare http://ocw.mit.edu 5.62 Physical Chemistry II Spring 2008 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms. 5.62 Lecture #20: Virial
More informationEQUATION OF STATE DEVELOPMENT
EQUATION OF STATE DEVELOPMENT I. Nieuwoudt* & M du Rand Institute for Thermal Separation Technology, Department of Chemical Engineering, University of Stellenbosch, Private bag X1, Matieland, 760, South
More informationCH.7 Fugacities in Liquid Mixtures: Models and Theories of Solutions
CH.7 Fugacities in Liquid Mixtures: Models and Theories of Solutions The aim of solution theory is to express the properties of liquid mixture in terms of intermolecular forces and liquid structure. The
More informationChE 524 A. Z. Panagiotopoulos 1
ChE 524 A. Z. Panagiotopoulos 1 VIRIAL EXPANSIONS 1 As derived previously, at the limit of low densities, all classical fluids approach ideal-gas behavior: P = k B T (1) Consider the canonical partition
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 informationMolecular Dynamics Simulation Study of Transport Properties of Diatomic Gases
MD Simulation of Diatomic Gases Bull. Korean Chem. Soc. 14, Vol. 35, No. 1 357 http://dx.doi.org/1.51/bkcs.14.35.1.357 Molecular Dynamics Simulation Study of Transport Properties of Diatomic Gases Song
More informationBinary Hard-Sphere Mixtures Within Spherical Pores
Journal of the Korean Physical Society, Vol. 35, No. 4, October 1999, pp. 350 354 Binary Hard-Sphere Mixtures Within Spherical Pores Soon-Chul Kim Department of Physics, Andong National University, Andong
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 informationPre-yield non-affine fluctuations and a hidden critical point in strained crystals
Supplementary Information for: Pre-yield non-affine fluctuations and a hidden critical point in strained crystals Tamoghna Das, a,b Saswati Ganguly, b Surajit Sengupta c and Madan Rao d a Collective Interactions
More informationCHEM-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 informationChem 4501 Introduction to Thermodynamics, 3 Credits Kinetics, and Statistical Mechanics
Chem 4501 Introduction to hermodynamics, 3 Credits Kinetics, and Statistical Mechanics Module Number 2 Active Learning Answers and Optional Problems/Solutions 1. McQuarrie and Simon, 2-6. Paraphrase: How
More informationThermodynamic study of liquid with solver ζ & suitable η max as a pole in basic two parameter Khasare s equation of state
Available online at www.pelagiaresearchlibrary.com Advances in Applied Science Research,, 3 (5):353-359 ISSN: 976-86 CODEN (USA): AASRFC Thermodynamic study of liquid with solver ζ & suitable η max as
More informationA New Uniform Phase Bridge Functional: Test and Its Application to Non-uniform Phase Fluid
Commun. Theor. Phys. (Beijing, China) 39 (2003) pp. 231 237 c International Academic Publishers Vol. 39, No. 2, February 15, 2003 A New Uniform Phase Bridge Functional: Test and Its Application to Non-uniform
More informationMOLECULAR DYNAMICS SIMULATION OF HETEROGENEOUS NUCLEATION OF LIQUID DROPLET ON SOLID SURFACE
MOLECULAR DYNAMICS SIMULATION OF HETEROGENEOUS NUCLEATION OF LIQUID DROPLET ON SOLID SURFACE Tatsuto Kimura* and Shigeo Maruyama** *Department of Mechanical Engineering, The University of Tokyo 7-- Hongo,
More informationCritical Properties of Isobaric Processes of Lennard-Jones Gases
Critical Properties of Isobaric Processes of Lennard-Jones Gases Akira Matsumoto Department of Material Sciences, College of Integrated Arts Sciences, Osaka Prefecture University, Sakai, Osaka, 599-8531,
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 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 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 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 informationFinal Exam, Chemistry 481, 77 December 2016
1 Final Exam, Chemistry 481, 77 December 216 Show all work for full credit Useful constants: h = 6.626 1 34 J s; c (speed of light) = 2.998 1 8 m s 1 k B = 1.387 1 23 J K 1 ; R (molar gas constant) = 8.314
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 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 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 informationA MOLECULAR DYNAMICS SIMULATION OF A BUBBLE NUCLEATION ON SOLID SURFACE
A MOLECULAR DYNAMICS SIMULATION OF A BUBBLE NUCLEATION ON SOLID SURFACE Shigeo Maruyama and Tatsuto Kimura Department of Mechanical Engineering The University of Tokyo 7-- Hongo, Bunkyo-ku, Tokyo -866,
More informationVariational Approach to the Equilibrium Thermodynamic Properties of Simple Liquids. I*
THE JOURNAL OF CHEMICAL PHYSICS VOLUME 51, NUMBER 11 1 DECEMBER 1969 Variational Approach to the Equilibrium Thermodynamic Properties of Simple Liquids. I* G.Ali MANSOORi (1) AND Frank B. CANFIELD (2)
More informationarxiv: v1 [cond-mat.soft] 20 Jun 2008
Accurate determination of crystal structures based on averaged local bond order parameters Wolfgang Lechner and Christoph Dellago Faculty of hysics, University of Vienna, Boltzmanngasse, 19 Vienna, Austria
More informationICCP Project 2 - Advanced Monte Carlo Methods Choose one of the three options below
ICCP Project 2 - Advanced Monte Carlo Methods Choose one of the three options below Introduction In statistical physics Monte Carlo methods are considered to have started in the Manhattan project (1940
More informationThree Semi-empirical Analytic Expressions for the Radial Distribution Function of Hard Spheres
Commun. Theor. Phys. (Beijing, China) 4 (2004) pp. 400 404 c International Academic Publishers Vol. 4, No. 3, March 5, 2004 Three Semi-empirical Analytic Expressions for the Radial Distribution Function
More informationCH 240 Chemical Engineering Thermodynamics Spring 2007
CH 240 Chemical Engineering Thermodynamics Spring 2007 Instructor: Nitash P. Balsara, nbalsara@berkeley.edu Graduate Assistant: Paul Albertus, albertus@berkeley.edu Course Description Covers classical
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 informationTheoretical comparative study on hydrogen storage of BC 3 and carbon nanotubes
J. At. Mol. Sci. doi: 10.4208/jams.121011.011412a Vol. 3, No. 4, pp. 367-374 November 2012 Theoretical comparative study on hydrogen storage of BC 3 and carbon nanotubes Xiu-Ying Liu a,, Li-Ying Zhang
More informationFig. 3.1? Hard core potential
6 Hard Sphere Gas The interactions between the atoms or molecules of a real gas comprise a strong repulsion at short distances and a weak attraction at long distances Both of these are important in determining
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 informationTheory of infinite dilution chemical potential
Fluid Phase Equilibria Journal Volume 85, 141-151, 1993 141 Esam Z. Hamada and G.Ali Mansoorib, * a Department of Chemical Engineering, King Fahd University of Petroleum and Minerals, Dhahran 31261 (Saudi
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 informationMultiple time step Monte Carlo
JOURNAL OF CHEMICAL PHYSICS VOLUME 117, NUMBER 18 8 NOVEMBER 2002 Multiple time step Monte Carlo Balázs Hetényi a) Department of Chemistry, Princeton University, Princeton, NJ 08544 and Department of Chemistry
More informationScientific Computing II
Scientific Computing II Molecular Dynamics Simulation Michael Bader SCCS Summer Term 2015 Molecular Dynamics Simulation, Summer Term 2015 1 Continuum Mechanics for Fluid Mechanics? Molecular Dynamics the
More informationMD simulation of methane in nanochannels
MD simulation of methane in nanochannels COCIM, Arica, Chile M. Horsch, M. Heitzig, and J. Vrabec University of Stuttgart November 6, 2008 Scope and structure Molecular model for graphite and the fluid-wall
More informationA new algorithm for Reverse Monte Carlo simulations
A new algorithm for Reverse Monte Carlo simulations Fernando Lus B. da Silva, Bo Svensson, Torbjörn Åkesson, and Bo Jönsson Citation: The Journal of Chemical Physics 109, 2624 (1998); doi: 10.1063/1.476861
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 informationA Nobel Prize for Molecular Dynamics and QM/MM What is Classical Molecular Dynamics? Simulation of explicit particles (atoms, ions,... ) Particles interact via relatively simple analytical potential
More 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 information2. Derive ideal mixing and the Flory-Huggins models from the van der Waals mixture partition function.
Lecture #5 1 Lecture 5 Objectives: 1. Identify athermal and residual terms from the van der Waals mixture partition function.. Derive ideal mixing and the Flory-Huggins models from the van der Waals mixture
More informationElastic constants and the effect of strain on monovacancy concentration in fcc hard-sphere crystals
PHYSICAL REVIEW B 70, 214113 (2004) Elastic constants and the effect of strain on monovacancy concentration in fcc hard-sphere crystals Sang Kyu Kwak and David A. Kofke Department of Chemical and Biological
More information1. Thermodynamics 1.1. A macroscopic view of matter
1. Thermodynamics 1.1. A macroscopic view of matter Intensive: independent of the amount of substance, e.g. temperature,pressure. Extensive: depends on the amount of substance, e.g. internal energy, enthalpy.
More informationComputer simulation study of liquid CH 2 F 2 with a new effective pair potential model
JOURNAL OF CHEMICAL PHYSICS VOLUME 110, NUMBER 6 8 FEBRUARY 1999 Computer simulation study of liquid CH 2 F 2 with a new effective pair potential model Pál Jedlovszky a),b) and Mihaly Mezei Department
More informationMOLECULAR DYNAMICS SIMULATION OF VAPOR BUBBLE NUCLEATION ON A SOLID SURFACE. Tatsuto Kimura and Shigeo Maruyama
MOLECULAR DYNAMICS SIMULATION OF VAPOR BUBBLE NUCLEATION ON A SOLID SURFACE Tatsuto Kimura and Shigeo Maruyama * Department of Mechanical Engineering, The University of Tokyo, 7-- Hongo, Bunkyo-ku, Tokyo
More informationAdvantages of a Finite Extensible Nonlinear Elastic Potential in Lattice Boltzmann Simulations
The Hilltop Review Volume 7 Issue 1 Winter 2014 Article 10 December 2014 Advantages of a Finite Extensible Nonlinear Elastic Potential in Lattice Boltzmann Simulations Tai-Hsien Wu Western Michigan University
More informationMolecular Dynamics Study on the Binary Collision of Nanometer-Sized Droplets of Liquid Argon
Molecular Dynamics Study on the Binary Collision Bull. Korean Chem. Soc. 2011, Vol. 32, No. 6 2027 DOI 10.5012/bkcs.2011.32.6.2027 Molecular Dynamics Study on the Binary Collision of Nanometer-Sized Droplets
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 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 informationDevelopment of a Water Cluster Evaporation Model using Molecular Dynamics
Development of a Water Cluster Evaporation Model using Molecular Dynamics Arnaud Borner, Zheng Li, Deborah A. Levin. Department of Aerospace Engineering, The Pennsylvania State University, University Park,
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 informationPolymer 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 informationto satisfy the large number approximations, W W sys can be small.
Chapter 12. The canonical ensemble To discuss systems at constant T, we need to embed them with a diathermal wall in a heat bath. Note that only the system and bath need to be large for W tot and W bath
More informationSunyia Hussain 06/15/2012 ChE210D final project. Hydration Dynamics at a Hydrophobic Surface. Abstract:
Hydration Dynamics at a Hydrophobic Surface Sunyia Hussain 6/15/212 ChE21D final project Abstract: Water is the universal solvent of life, crucial to the function of all biomolecules. Proteins, membranes,
More informationHow the Nonextensivity Parameter Affects Energy Fluctuations. (Received 10 October 2013, Accepted 7 February 2014)
Regular Article PHYSICAL CHEMISTRY RESEARCH Published by the Iranian Chemical Society www.physchemres.org info@physchemres.org Phys. Chem. Res., Vol. 2, No. 2, 37-45, December 204. How the Nonextensivity
More informationStatistical physics. Emmanuel Trizac LPTMS / University Paris-Sud
Statistical physics Emmanuel Trizac LPTMS / University Paris-Sud Outline I Introduction to phase transitions and critical phenomena 1- The problems raised by phase transitions, from a statistical mechanics
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 informationVelocity cross-correlations and atomic momentum transfer in simple liquids with different potential cores
PHYSICAL REVIEW E VOLUME 62, NUMBER 1 JULY 2000 Velocity cross-correlations and atomic momentum transfer in simple liquids with different potential cores A. Verdaguer and J. A. Padró Departament de Física
More informationDerivation of Van der Waal s equation of state in microcanonical ensemble formulation
arxiv:180.01963v1 [physics.gen-ph] 9 Nov 017 Derivation of an der Waal s equation of state in microcanonical ensemble formulation Aravind P. Babu, Kiran S. Kumar and M. Ponmurugan* Department of Physics,
More informationStructure of the First and Second Neighbor Shells of Water: Quantitative Relation with Translational and Orientational Order.
Structure of the First and Second Neighbor Shells of Water: Quantitative Relation with Translational and Orientational Order Zhenyu Yan, Sergey V. Buldyrev,, Pradeep Kumar, Nicolas Giovambattista 3, Pablo
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 informationAtomic and molecular interaction forces in biology
Atomic and molecular interaction forces in biology 1 Outline Types of interactions relevant to biology Van der Waals interactions H-bond interactions Some properties of water Hydrophobic effect 2 Types
More informationHyeyoung 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 informationA theory for calculating the number density distribution of. small particles on a flat wall from pressure between the two
theory for calculating the number density distribution of small particles on a flat wall from pressure between the two walls Kota Hashimoto a and Ken-ichi mano a a Department of Energy and Hydrocarbon
More informationFrom Dynamics to Thermodynamics using Molecular Simulation
From Dynamics to Thermodynamics using Molecular Simulation David van der Spoel Computational Chemistry Physical models to describe molecules Software to evaluate models and do predictions - GROMACS Model
More informationPerturbation approach for equation of state for hard-sphere and Lennard Jones pure fluids
PRAMANA c Indian Academy of Sciences Vol. 76, No. 6 journal of June 2011 physics pp. 901 908 Perturbation approach for equation of state for hard-sphere and Lennard Jones pure fluids S B KHASARE and M
More informationWhat is Classical Molecular Dynamics?
What is Classical Molecular Dynamics? Simulation of explicit particles (atoms, ions,... ) Particles interact via relatively simple analytical potential functions Newton s equations of motion are integrated
More informationStructural transition and solid-like behavior of alkane films confined in nano-spacing
Fluid Phase Equilibria 183 184 (2001) 381 387 Structural transition and solid-like behavior of alkane films confined in nano-spacing S.T. Cui a,b,, P.T. Cummings b,c, H.D. Cochran a,b a Department of Chemical
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 informationStatistical Theory and Learning from Molecular Simulations
Statistical Theory and Learning from Molecular Simulations Lawrence R. Pratt 1 and Susan B. Rempe 2 1 Department of Chemical & Biomolecular Engineering 2 Center for Biological and Material Sciences, Sandia
More informationSurface Tension of the Vapor Liquid Interface with Finite Curvature
Colloid Journal, Vol. 65, No. 4, 00, pp. 44045. Translated from Kolloidnyi Zhurnal, Vol. 65, No. 4, 00, pp. 480494. Original Russian Text Copyright 00 by Zhukhovitskii. Surface Tension of the VaporLiquid
More information2. 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 informationMultiphase Flow and Heat Transfer
Multiphase Flow and Heat Transfer Liquid-Vapor Interface Sudheer Siddapuredddy sudheer@iitp.ac.in Department of Mechanical Engineering Indian Institution of Technology Patna Multiphase Flow and Heat Transfer
More informationScattering Properties of Gas Molecules on a Water Adsorbed Surface
Scattering Properties of Gas Molecules on a Water Adsorbed Surface Hideki Takeuchi a, Koji Yamamoto b, and Toru Hakutake c a Department of Mechanical Engineering, Kochi National College of Technolog, Nankoku
More informationGuest-Guest Interactions between Small Molecules and THF in Binary Clathrate Hydrates
Guest-Guest Interactions between Small Molecules and THF in Binary Clathrate Hydrates Abstract Tim Strobel CHGN 54 Term Project April 26 th, 26 The application of the statistical thermodynamic model, originally
More informationDENSITY FUNCTIONAL THEORY FOR STUDIES OF MULTIPLE STATES OF INHOMOGENEOUS FLUIDS AT SOLID SURFACES AND IN PORES.
J. Smith, D. Stroud, MRS Symposium Proceedings Series, v.49, p.7-33, 998. DENSITY FUNCTIONAL THEORY FOR STUDIES OF MULTIPLE STATES OF INHOMOGENEOUS FLUIDS AT SOLID SURFACES AND IN PORES. A.. NEIMARK, and
More informationStatistical Physics. Solutions Sheet 11.
Statistical Physics. Solutions Sheet. Exercise. HS 0 Prof. Manfred Sigrist Condensation and crystallization in the lattice gas model. The lattice gas model is obtained by dividing the volume V into microscopic
More informationOn 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 informationMultiscale 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 informationFrom Atoms to Materials: Predictive Theory and Simulations
From Atoms to Materials: Predictive Theory and Simulations Week 3 Lecture 4 Potentials for metals and semiconductors Ale Strachan strachan@purdue.edu School of Materials Engineering & Birck anotechnology
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 informationLennard-Jones as a model for argon and test of extended renormalization group calculations
JOURNAL OF CHEMICAL PHYSICS VOLUME 111, NUMBER 2 22 NOVEMBER 1999 Lennard-Jones as a model for argon and test of extended renormalization group calculations John A. White Department of Physics, American
More informationProperties of real fluids in critical region: third virial coefficient
Indian J hys (February 2014) 88(2):185 191 DOI 10.1007/s12648-013-0402-5 ORIGINAL AER roperties of real fluids in critical region: third virial coefficient R Khordad*, B Mirhosseini and M M Mirhosseini
More informationMolecular Dynamics Simulation of Diffusion of Gases and Liquids in Conditions of Phase Transition
, October 19-21, 2016, San Francisco, USA Molecular Dynamics Simulation of Diffusion of Gases and Liquids in Conditions of Phase Transition Georgii V. Kharlamov, Sergey V. Zhilkin Abstract The diffusion
More informationPressure Dependent Study of the Solid-Solid Phase Change in 38-Atom Lennard-Jones Cluster
University of Rhode Island DigitalCommons@URI Chemistry Faculty Publications Chemistry 2005 Pressure Dependent Study of the Solid-Solid Phase Change in 38-Atom Lennard-Jones Cluster Dubravko Sabo University
More information