A Nonequilibrium Molecular Dynamics Study of. the Rheology of Alkanes. S.A. Gupta, S. T. Cui, P. T. Cummings and H. D. Cochran
|
|
- Isabel Holt
- 6 years ago
- Views:
Transcription
1 A Nonequilibrium Molecular Dynamics Study of the Rheology of Alkanes S.A. Gupta, S. T. Cui, P. T. Cummings and H. D. Cochran Department of Chemical Engineering University of Tennessee Knoxville, TN and Chemical Technology Division Oak Ridge National Laboratory Oak Ridge, TN Abstract We examine the rheological properties of four dierent alkanes - n-decane (C 10 H 22 ), n-hexadecane (C 16 H 34 ), n-tetracosane (C 24 H 50 ) and squalane (C 30 H 62 ). The simulations of Couette ow are performed for a range of shear rates with 100 molecules in each case using a replicated data version of our code. The number of interaction sites ranges from 1000? We have performed extremely long simulations required to obtain acceptable statistics at low shear rates. The alkanes show a transition from non-newtonian to Newtonian behavior as the shear rate decreases to low values. KEYWORDS: Nonequilibrium Molecular Dynamics, rheology, alkanes 1 Introduction The use of molecular-based approaches to address the problem of lubricant performance has gained signicant momentum in the past few years. Much of the work has focused on relatively simple systems, the primary limitation being the high computational cost involved in characterizing these lubricants. The use of simple systems cannot give useful quantitative predictions for actual lubricants. The calculations that we have performed to determine the rheological properties of liquid alkanes in the intermediate molecular ptc@utk.edu 1
2 size range (C 20? C 40 ), can be eciently done with codes running on massively parallel supercomputers. Such calculations allow systematic studies of these systems with the aim of identifying the most important properties aecting their performance as lubricants. They also aid the design of synthetic lubricants with desired properties. Planar Couette ow has been the focus of active research using non-equilibrium molecular dynamics, in part due to its applicability in lubrication science. The theoretical foundations of this non-equilibrium system and the non-equilibrium molecular dynamics (NEMD) techniques developed for its study have been addressed [1]. In most studies of planar Couette ow, a reversible deterministic thermostat such as a Nose thermostat is incorporated into the equations of motion to permit a steady state to develop. The thermostatted SLLOD equations of motion are solved under the so-called Lees-Edwards periodic boundary conditions. The interested reader is referred to [1] for technical details. A common approach in simulations of chain molecule systems, like alkanes, has been to use the method of constraints for the internal degrees of freedom, e.g. bond bending and stretching modes. This method is computationally advantageous since it allows using a larger time step than that permissible for a fully exible model. Another approach to decreasing cpu time is to utilize the short-ranged purely repulsive potential model for the intermolecular site-site interaction. This potential was used in the NEMD simulations of eicosane (C 20 H 42 ) [2]. The disadvantage of the repulsive potential model is that it fails to yield quantitative comparison with experimental data. A more realistic simulation using the Ryckaert{Bellemans model [3] was used to determine the rheological properties of n-hexadecane [4] at a relatively high temperature and density. The calculated Newtonian viscosity was lower than the experimental value by a factor of almost three. To obtain improved quantitative predictions of the equilibrium and non-equilibrium properties of a molecular uid, it becomes essential to use a realistic and exible molecular model, at the expense of increased computational requirements. Such a model was used to predict accurately the vapor-liquid phase envelope of linear alkanes in the range C 6? C 48 [5]. The model was subsequently used by Mundy et al. [6] and Cui et al. [7] to predict quantitatively the viscosity of n-decane at a liquid state condition similar to those in the phase envelope calculation (high temperature, moderately high density). Mondello and Grest [8] used the same model to determine the equilibrium properties of linear and branched alkanes at near ambient temperature conditions, obtaining reasonable agreement with experimental results. There is a signicant practical interest in extending these studies to medium size alkanes (C 20? C 40 ) as the main constituents of lubricant basestocks. However, given the enormous computational requirements of such systems, it is essential to develop ecient parallel NEMD algorithms. Our approach uses a fully exible (including stretching, bending, and torsion) united atom model to describe the molecular sites. To aid in computational eciency we use a multiple timestep method based on a phase variable Liouvillean derived from the SLLOD equations of motion [1]. The work of Cui et al [9] implements the algorithm on the Intel Paragon in order to enable the NEMD simulation of long-chains. In this paper, we present our computational approach and the results of such an NEMD study to determine 2
3 the rheological properties of four dierent alkanes - n-decane (C 10 H 22 ), n-hexadecane (C 16 H 34 ), n-tetracosane (C 24 H 50 ) and squalane (C 30 H 62 ). The rst three alkanes are linear while squalane is a molecule with six symmetrically placed, methyl side-groups. Parallel computers make it possible to run the extremely long simulations required to obtain improved statistics at low shear rates. 2 Model, Theory and Simulation method For a system subject to planar Couette ow, the motion of the particles can be described by the SLLOD algorithm which is known to be valid in both the linear and nonlinear response regimes [1]. The SLLOD equations of motion (incorporating a Nose thermostat) for a system undergoing planar Couette are given by _r ia = p ia m ia + y ia^x; _p ia = F ia? p y;ia^x? p _ ia _ = p (1) Q _p = F where r ia and p ia are the vector coordinates and peculiar momentum of atom a in molecule i, y ia and p y;ia are its y components, m ia is the mass, ^x is a unit vector in the x direction,, p, and Q are the variables related to the Nose thermostat, F = X i;a p 2 ia m ia? 3Nk B T; Q = 3Nk B T 2 is the Nose thermostat time constant, and N is the total number of atoms in the system. The model used for the alkanes is the united atom model previously used by Siepmann et al. [5], [10], Mundy et al. [6], and by Cui et al. [7]. In this model the methyl and methylene groups are treated as spherical interaction sites with the interaction centers located at the center of the carbon atoms. Details of the model can be found in the literature [5], [10]. For completeness, we briey describe the model. The interaction between atoms on dierent molecules and atoms separated by more than three bonds in the same molecule is described by a Lennard-Jones (LJ) potential, u ij (r) = 4 ij " ij r 12 ij? r where ij is the well depth and ij is the zero point of the potential between a site of type i and a site of type j. The Lennard-Jones size parameters are CH3 = 3:93 A, CH2 = 3:93 A and CH = 3:81 A. The corresponding well depth parameters h =kb are 114 o K, 47 o K and 40 o K respectively. The Lorenz-Berthelot combining rules ij = ( i j ) 1=2 ; ij = i+ j i 2 are used for the cross interactions. A cuto distance of 9:825 A for the LJ interaction was 3 6 # (2)
4 used in this study. The bonded (intramolecular) interactions include a vibration term, a bending term and a dihedral term. The bond stretching and bond angle bending are described by harmonic potentials while the torsional potential is that by Jorgensen et al. [11]. For the simulations of squalane, an ad hoc harmonic potential similar to the bending term, is introduced to prevent umbrella inversion of the sp 3 bond conguration at tertiary carbon branch points. This potential was used by Mondello and Grest [8] and it forces the normal to the plane dened by the two backbone bonds joining at the CH group to oscillate around its equilibrium position, which forms an angle of 27:25 o with the normal to the plane dened by the three CH n groups connected to the tertiary carbon. The state points studied correspond to the densities of the liquid alkanes at atmospheric pressure and at the respective temperatures [12], [13]. These are listed in Table 1. In an NEMD calculation, the strain-rate dependent viscosity is determined from the constitutive relation = hp xyi ; (3) where hp xy i is the average of the xy component of the pressure tensor P and is the strain rate characterizing the shear eld. We have chosen the x direction to be the ow direction and the y direction to be the ow gradient direction, so that x =@y, where u x is the streaming velocity in the x direction. The pressure tensor P is calculated using the atomic formalism, which is given by P (a) V = X i;a m ia v ia v ia X X i;a i;b {z } i6=j (r ia? r jb )f ia;jb + X i;a r ia f (intra) ia ; (4) where the indices i and j refer to molecules i and j, indices a and b refer to interaction sites a and b in molecules i and j, respectively. The shear alignment is calculated using the order tensor dened by the second rank tensor [14], S = 3 * 1 NX 1 + (e i e i? 2 N 3 1) (5) i=1 Table 1: State points simulated for the liquid alkanes. Alkane T o K (gm=cm 3 ) n-decane 298 0:7247 a n-hexadecane (A) 300 0:770 b n-hexadecane (B) 323 0:753 b n-tetracosane 333 0:7728 b squalane 333 0:7837 b a Ref. [12], b Ref. [13] 4
5 where e i is the unit vector along the end-to-end direction of the molecule i, 1 is a unit second rank tensor and the summation is over all N molecules in the system. The algorithm used for the NEMD simulations is an extension to sheared systems [7] of the application to alkanes of the multi-timestep Nose dynamics algorithm developed in the Tuckermann-Berne reversible reference system propagator algorithm method [15]. Details can be found in [7]. We use two timestep sizes, with all intramolecular interactions being treated as fast motions and the intermolecular interactions as slow motions. The simulations of Couette ow are performed for a range of shear rates with 100 molecules in each case using a replicated data version of our code. The number of interaction sites ranges from 1000? A timestep of 2:35 fs is used for the slow mode of motion. For decane, hexadecane and tetracosane, a time step of 0:235 fs is used for the fast mode of motion, while the corresponding value for squalane is 0:5875 f s. The strain rate dependent properties of the systems are calculated based on the simulation runs after the systems have reached steady state under the inuence of the shear eld. 3 Replicated data code For rapid parallelization, we have used a replicated data strategy. This method is essentially what Plimpton [16] termed an \atomic decomposition" method, and was chosen for its simplicity and eectiveness in parallelizing an existing serial code. Since the action of the serial code has already been proven, the ability to rapidly parallelize an existing, scientically-complex code aids productivity enormously. The general strategy relies upon the fact that the most compute-intensive portion of the code is the calculation of the system of forces at work. There are basically two classes of force interactions considered, i.e., intermolecular and intramolecular. For the intermolecular portion of the force calculation, we require data characterizing the entire system state be present upon every node; hence the term \replicated data" refers to the replication of all of the particle positions. The force is evaluated within a double loop that compares uniquely every particle with every other particle. The use of a neighbor list optimizes the calculation time, by considering only those particles that are within a certain cuto distance from a given particle. The computational work is divided up among nodes by assigning a range of outer loop iterations to each node. These can be worked upon in parallel, and the simplicity of the division permits some straightforward load balancing to take place. At least some form of load balancing is required due to the fact that the quantity of computational work in the double loop structure is not uniformly distributed. In fact, due to the neighbor list being used as an optimization strategy, the earlier iterations will carry much more work than later iterations. To counter this imbalance, we have simply divided the work for a set of N nodes, as if we had 2N available nodes. We then assign iterations from opposite ends of the outer loop, in eect, giving each node a large and a small portion to work on. This method has the advantage that we can decompose the work according to the number of particles divided by 2, in the case of maximum parallelism. 5
6 Further parallelism can be realized if we pursue inner loop distribution as well, but this is currently not advantageous. In addition to the eectiveness of this technique, it is very straightforward to implement. Although currently a static load balancing technique is used, a dynamic action could be easily added for little extra cost, if it were required. For the intramolecular portion of the force calculation, there is a dierent action at work. Due to the multi-timestep approach adopted in this simulation, the intramolecular arrangement changes prior to the calculation of the intermolecular forces. These calculations require only local information regarding the particle positions of a single molecule, but that also xes the maximum decomposition to the molecular level. A ner level of molecular decomposition might be possible for larger molecules, but would probably require repeating calculations to avoid additional communication between nodes. Hence, in the strategy we have implemented, two communications are required to ensure a correct simulation action. One is after the integration step following the intramolecular force evaluation, and the other follows the intermolecular force evaluation. The rst communication collects and distributes all of the particle positions onto all of the processors. The intermolecular forces are then evaluated and these atomic forces are then summed globally, so that all nodes have sucient information to update the system via integration, and commence the next round of intramolecular calculations. 4 Simulation Results The longest simulations h reported in this study i are for about 9 million timesteps. The dimensionless viscosity = 2 =(m) 1=2 is plotted in Fig. 1 for the four alkanes, as a function of the dimensionless strain rate h = (m 2 =) 1=2 i. These state points correspond to the equilibrium density at atmospheric pressure and at the respective temperatures. Also plotted in Fig. 1 for comparative purposes, is the viscosity for methane (CH 4 ) at the Lennard-Jones triple point calculated using a short-ranged purely repulsive potential [17]. As can be seen from the gure, the viscosities for the alkanes exhibit a shear thinning behavior over the range of strain rates studied. This is typical of chain molecule uid systems. At high strain rates, the shear thinning follows a power law. The slopes of the log-log plots vary from to as compared to the experimentally observed slopes of -0.4 to -0.9 for polymeric uids. This likeness suggests that even though the alkane chains studied here are short in comparison to polymer systems, they nevertheless exhibit some of the generality of long-chain systems. Another important observation is that the shear viscosities for decane, hexadecane and tetracosane nearly overlap each other at high strain rates. At high strain rates these fairly short and sti alkane chains are well aligned with each other so they can slide past each other relatively easily. Furthermore, as shown by the calculated alignment angles, the longer chains align themselves better resulting in the above mentioned insensitivity of viscosity with the variation of chain length at high strain rate. For squalane, which has the same backbone chain length as tetracosane, the 6 additional methyl side groups 6
7 100 C1 (90 K) C10 (298 K) C16 (300 K) C16 (323 K) C24 (333 K) C30 (333 K) η γ Figure 1: Dimensionless viscosity ( ) for linear alkanes as a function of the dimensionless shear rate ( ) result in viscosities at high strain rate that are about twice those of tetracosane. The alignment angles for squalane at these strain rates are identical to those for tetracosane within statistical uncertainty. 5 Summary We have developed a replicated data nonequilibrium molecular dynamics code to determine the rheological properties of decane, hexadecane, tetracosane, and squalane. Due to the long simulations required to obtain good statistics for these chain molecules, these simulations can be eciently performed only on a massively parallel supercomputer. Work is currently in progress to study other branched alkanes. Acknowledgments The work of HDC has been supported by the Division of Chemical Sciences of the U. S. Department of Energy (DOE) at Oak Ridge National Laboratory (ORNL). This work was sponsored by the Laboratory Directed Research and Development Program of ORNL. The authors acknowledge the use of the Intel Paragon supercomputers in the Center 7
8 for Computational Sciences(CCS) at ORNL, funded by the DOE's Mathematical, Information, and Computational Sciences Division. ORNL is managed by Lockheed Martin Energy Research Corp. for the DOE under Contract No. DE-AC05-96OR The authors also thank Phil LoCascio of CCS for his participation in this work. References [1] D. J. Evans and G. P. Morriss, Statistical Mechanics of Nonequilibrium Liquids (Academic Press, New York, 1990). [2] G. P. Morriss, P. T. Daivis, and D. J. Evans, J. Chem. Phys. 94, 7420 (1991). [3] J. P. Ryckaert and A. Bellemans, Chem. Phys. Lett., 30, 123 (1975); Faraday Discuss. Chem. Soc., 66, 95 (1978). [4] A. Berker, S. Chinoweth, U. C. Klomp, and Y. Michopoulos, J. Chem. Soc. Faraday Trans. 88, 1719 (1992). [5] J. I. Siepmann, S. Karaboni, and B. Smit, Nature 365, 330 (1993). [6] C. J. Mundy, J. I. Siepmann, and M. L. Klein, J. Chem. Phys. 102, 3376 (1995). [7] S. T. Cui, P. T. Cummings, and H. D. Cochran, J. Chem. Phys. 104, 255 (1996). [8] M. Mondello and G. S. Grest, J. Chem. Phys. 103, 7156 (1995). [9] S. T. Cui, S. A. Gupta, P. T. Cummings, and H. D. Cochran, Molecular Dynamics Simulations of the Rheology of Normal Decane, Hexadecane and Tetracosane, J. Chem. Phys. (accepted for publication). [10] B. Smit, S. Karaboni, and J. I. Siepmann, J. Chem. Phys. 102, 2126 (1995). [11] W. L. Jorgensen, J. D. Madura, and C. J. Swenson, J. Am. Chem. Soc. 106, 6638 (1984). [12] M. Gehrig and H. Lentz, J. Chem. Thermodynamics 15, 1159 (1983). [13] API 42, Properties of hydrocarbons of high molecular weight (American Petroleum Institute, Research Project 42, Washington, DC, 1966). [14] S. Chandrasekhar, Liquid Crystals, 2nd ed. (Cambridge University Press, Cambridge, 1992); P. G. degennes and J. Prost, The Physics of Liquid Crystals, 2nd ed. (Clarendon Press, Oxford, 1993). [15] M. E. Tuckerman, B. J. Berne, and G. J. Martyna, J. Chem. Phys. 97, 1990, (1992). 8
9 [16] S. Plimpton, Fast Parallel Algorithms for Short-Range Molecular Dynamics, Technical Report SAND , Sandia National Laboratories, May [17] R. K. Bhupathiraju, P. T. Cummings, and H. D. Cochran, An ecient parallel algorithm for nonequilibrium molecular dynamics simulations of very large systems in planar Couette ow, Mol. Phys. (accepted for publication). 9
Nanorheology of liquid alkanes
Ž. Fluid Phase Equilibria 150 151 1998 125 131 Nanorheology of liquid alkanes S.A. Gupta a,b, H.D. Cochran a,b, P.T. Cummings a,b,) a Department of Chemical Engineering, UniÕersity of Tennessee, KnoxÕille,
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 informationOleg A. Mazyar, Guoai Pan and Clare McCabe*
Molecular Physics Vol. 107, No. 14, 20 July 2009, 1423 1429 RESEARCH ARTICLE Transient time correlation function calculation of the viscosity of a molecular fluid at low shear rates: a comparison of stress
More informationNonequilibrium Molecular Dynamics Simulations of Shear Viscosity: Isoamyl Alcohol, n-butyl Acetate and Their Mixtures
Draft Date: October 5, 1999 Nonequilibrium Molecular Dynamics Simulations of Shear Viscosity: Isoamyl Alcohol, n-butyl Acetate and Their Mixtures Y. Yang 1, T. A. Pakkanen, 2 and R. L. Rowley 1,3 Draft
More informationPredicting the viscosity of alkanes using nonequilibrium molecular dynamics: Evaluation of intermolecular potential models
Brigham Young University BYU ScholarsArchive All Faculty Publications 1997-06-22 Predicting the viscosity of alkanes using nonequilibrium molecular dynamics: Evaluation of intermolecular potential models
More informationEfficient viscosity estimation from molecular dynamics simulation via momentum impulse relaxation
JOURNAL OF CHEMICAL PHYSICS VOLUME 113, NUMBER 6 8 AUGUST 2000 Efficient viscosity estimation from molecular dynamics simulation via momentum impulse relaxation Gaurav Arya, Edward J. Maginn, a) and Hsueh-Chia
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 informationNosé-Hoover chain method for nonequilibrium molecular dynamics simulation
PHYSICAL REVIEW E VOLUME 61, NUMBER 5 MAY 000 Nosé-Hoover chain method for nonequilibrium molecular dynamics simulation A. C. Brańka Institute of Molecular Physics, Polish Academy of Sciences, Smoluchowskiego
More informationSupporting Information
Projection of atomistic simulation data for the dynamics of entangled polymers onto the tube theory: Calculation of the segment survival probability function and comparison with modern tube models Pavlos
More informationCode Timings of a Bulk Liquid Simulation
1 Parallel Molecular Dynamic Code for Large Simulations using Truncated Octahedron Periodics M.M. Micci, L.N. Long and J.K. Little a a Aerospace Department, The Pennsylvania State University, 233 Hammond
More informationInstitute for Condensed Matter Theory and Beckman Institute, University of Illinois at Urbana-Champaign, Urbana, IL 61801, USA
Shear-thinning of molecular fluids in Couette flow Bharath V. Raghavan 1 and Martin Ostoja-Starzewski 2 1) Department of Mechanical Science and Engineering,University of Illinois at Urbana-Champaign, Urbana,
More information4. The Green Kubo Relations
4. The Green Kubo Relations 4.1 The Langevin Equation In 1828 the botanist Robert Brown observed the motion of pollen grains suspended in a fluid. Although the system was allowed to come to equilibrium,
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 informationNon-equilibrium molecular dynamics simulation study of the behavior of hydrocarbon-isomers in silicalite
Fluid Phase Equilibria 194 197 (2002) 309 317 Non-equilibrium molecular dynamics simulation study of the behavior of hydrocarbon-isomers in silicalite S. Furukawa a,b,, C. McCabe a,c, T. Nitta b, P.T.
More informationNON-EQUILIBRIUM MOLECULAR DYNAMICS CALCULATION OF THE
FluidPhaseEquilibria, 53 (1989) 191-198 Elsevier Science Publishers B.V., Amsterdam - Printed in The Netherlands 191 NON-EQUILIBRIUM MOLECULAR DYNAMICS CALCULATION OF THE TRANSPORT PROPERTIES OF CARBON
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 informationin order to insure that the Liouville equation for f(?; t) is still valid. These equations of motion will give rise to a distribution function f(?; t)
G25.2651: Statistical Mechanics Notes for Lecture 21 Consider Hamilton's equations in the form I. CLASSICAL LINEAR RESPONSE THEORY _q i = @H @p i _p i =? @H @q i We noted early in the course that an ensemble
More informationSimulation of Nematic-Isotropic Phase Coexistence in Liquid Crystals under Shear
Simulation of Nematic-Isotropic Phase Coexistence in Liquid Crystals under Shear Guido Germano 1,2 and Friederike Schmid 2 1 Physical Chemistry Philipps-University 3532 Marburg, Germany E-mail: germano@staff.uni-marburg.de
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 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 informationDynamics of Individual Molecules of Linear Polyethylene Liquids under Shear: Atomistic Simulation and Comparison with a Free-draining Bead-rod Chain
University of Tennessee, Knoxville Trace: Tennessee Research and Creative Exchange Chemical and Biomolecular Engineering Publications and Other Works Chemical and Biomolecular Engineering 2-19-2010 Dynamics
More informationShear viscosity of liquid rubidium at the triple point
J. Phys. F: Met. Phys. 18 (1988) 1439-1447. Printed in the UK Shear viscosity of liquid rubidium at the triple point P T Cummingsi- and G P Morriss-l: i Department of Chemical Engineering, Thornton Hall,
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 informationMOLECULAR DYNAMICS SIMULATIONS OF FRICTION FORCE VERSUS LOAD
MOLECULAR DYNAMICS SIMULATIONS OF FRICTION FORCE VERSUS LOAD Ana-Camelia Pirghie Stefan cel Mare University of Suceava, Department of Mechanics and Technology, camelia.pirghie@fim.usv.ro Abstract: The
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 informationImproved model of nonaffine strain measure
Improved model of nonaffine strain measure S. T. Milner a) ExxonMobil Research & Engineering, Route 22 East, Annandale, New Jersey 08801 (Received 27 December 2000; final revision received 3 April 2001)
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 informationNon-Equilibrium Molecular Dynamics Approach to the Rheology of Model Polymer Fluids
Fluid Phase Equilibria, 88 (1993) 99-113 Elsevier Science Publishers B.V., msterdam 99 Non-Equilibrium Molecular Dynamics pproach to the Rheology of Model Polymer Fluids J. W. Rudisill and P. T. Cummings
More informationIntroduction to Nonequilibrium Molecular Dynamics
Introduction to Nonequilibrium Molecular Dynamics Dr Billy D. Todd Centre for Molecular Simulation Swinburne University of Technology PO Box 218, Hawthorn Victoria 3122, Australia Email: btodd@swin.edu.au
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 informationUniversities of Leeds, Sheffield and York
promoting access to White Rose research papers Universities of Leeds, Sheffield and York http://eprints.whiterose.ac.uk/ This is an author produced version of a paper published in The Journal of Chemical
More informationNon-equilibrium molecular dynamics simulations of organic friction modifiers
Non-equilibrium molecular dynamics simulations of organic friction modifiers James Ewen PhD Student Tribology Group (Shell UTC) Imperial College London Session 8B Bronze 2 Lubrication Fundamentals STLE
More informationShear viscosity of molten sodium chloride
Shear viscosity of molten sodium chloride Jerome Delhommelle and Janka Petravic Citation: The Journal of Chemical Physics 118, 2783 (2003); doi: 10.1063/1.1535213 View online: http://dx.doi.org/10.1063/1.1535213
More informationarxiv: v1 [cond-mat.mtrl-sci] 24 Jan 2012
Effective Viscosity of Confined Hydrocarbons I.M. Sivebaek,2,3, V.N. Samoilov,4 and B.N.J. Persson IFF, FZ-Jülich, 52425 Jülich, Germany 2 Novo Nordisk A/S, Research and Development, DK4 Hillerod, Denmark
More informationViscosity and Diffusion Constants Calculation of n-alkanes by Molecular Dynamics Simulations
Viscosity and Diffusion Constants Calculation of n-alkanes by Molecular Dynamics Simulations Song Hi Lee * and Taihyun Chang Department of Chemistry, Kyungsung University, Busan 608-736, Korea Department
More information2 Lyapunov Exponent exponent number. Lyapunov spectrum for colour conductivity of 8 WCA disks.
3 2 Lyapunov Exponent 1-1 -2-3 5 1 15 exponent number Lyapunov spectrum for colour conductivity of 8 WCA disks. 2 3 Lyapunov exponent 2 1-1 -2-3 5 1 15 exponent number Lyapunov spectrum for shear flow
More informationField indeced pattern simulation and spinodal. point in nematic liquid crystals. Chun Zheng Frank Lonberg Robert B. Meyer. June 18, 1995.
Field indeced pattern simulation and spinodal point in nematic liquid crystals Chun Zheng Frank Lonberg Robert B. Meyer June 18, 1995 Abstract We explore the novel periodic Freeredicksz Transition found
More informationMolecular Dynamics Simulation of a Nanoconfined Water Film
Molecular Dynamics Simulation of a Nanoconfined Water Film Kyle Lindquist, Shu-Han Chao May 7, 2013 1 Introduction The behavior of water confined in nano-scale environment is of interest in many applications.
More informationEffect of the Inner-Zone Vibrations on the Dynamics of Collision-Induced Intramolecular Energy Flow in Highly Excited Toluene
Notes Bull. Korean Chem. Soc. 2005, Vol. 26, No. 8 1269 Effect of the Inner-Zone Vibrations on the Dynamics of Collision-Induced Intramolecular Energy Flow in Highly Excited Toluene Jongbaik Ree, * Yoo
More informationlike carbon, has fewer than an octet. It is simply less likely but still imperative to draw.
Andrew Rosen Chapter 1: The Basics - Bonding and Molecular Structure 1.1 - We Are Stardust - Organic chemistry is simply the study of carbon-based compounds, hydrocarbons, and their derivatives, which
More informationOn the Eective Magnetic Properties of. Magnetorheological Fluids. November 9, Abstract
On the Eective Magnetic Properties of Magnetorheological Fluids November 9, 998 Tammy M. Simon, F. Reitich, M.R. Jolly, K. Ito, H.T. Banks Abstract Magnetorheological (MR) uids represent a class of smart
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 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 informationQuestion ( ) Solution. Approximate length of an atomic bond for solid aluminum.
Question (1450001) Approximate length of an atomic bond for solid aluminum. Spring A of stiness 10 N/m is hung vertically from a ringstand. From the lower end of this spring, a dierent spring (Spring B)
More informationIntroduction. Monday, January 6, 14
Introduction 1 Introduction Why to use a simulation Some examples of questions we can address 2 Molecular Simulations Molecular dynamics: solve equations of motion Monte Carlo: importance sampling Calculate
More informationAlkanes. Introduction
Introduction Alkanes Recall that alkanes are aliphatic hydrocarbons having C C and C H bonds. They can be categorized as acyclic or cyclic. Acyclic alkanes have the molecular formula C n H 2n+2 (where
More informationRheology and Constitutive Equations. Rheology = Greek verb to flow. Rheology is the study of the flow and deformation of materials.
Rheology and Constitutive Equations Rheology = Greek verb to flow Rheology is the study of the flow and deformation of materials. The focus of rheology is primarily on the study of fundamental, or constitutive,
More informationU U B U P x
Smooth Quantum Hydrodynamic Model Simulation of the Resonant Tunneling Diode Carl L. Gardner and Christian Ringhofer y Department of Mathematics Arizona State University Tempe, AZ 8587-184 Abstract Smooth
More informationTowards a computational chemical potential for nonequilibrium steady-state systems
PHYSICAL REVIEW E VOLUME 60, NUMBER 5 NOVEMBER 1999 Towards a computational chemical potential for nonequilibrium steady-state systems András Baranyai Department of Theoretical Chemistry, Eötvös University,
More informationJASS Modeling and visualization of molecular dynamic processes
JASS 2009 Konstantin Shefov Modeling and visualization of molecular dynamic processes St Petersburg State University, Physics faculty, Department of Computational Physics Supervisor PhD Stepanova Margarita
More informationvs. Chapter 4: Standard Flows Chapter 4: Standard Flows for Rheology shear elongation 2/1/2016 CM4650 Lectures 1-3: Intro, Mathematical Review
CM465 Lectures -3: Intro, Mathematical //6 Chapter 4: Standard Flows CM465 Polymer Rheology Michigan Tech Newtonian fluids: vs. non-newtonian fluids: How can we investigate non-newtonian behavior? CONSTANT
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 informationPharmaceutics I صيدالنيات 1. Unit 6
Pharmaceutics I صيدالنيات 1 Unit 6 1 Rheology of suspensions Rheology, the study of flow, addresses the viscosity characteristics of powders, fluids, and semisolids. Materials are divided into two general
More informationA. Bers, A. K. Ram, and S. D. Schultz. Plasma Science and Fusion Center,
COUPLING TO ELECTRON BERNSTEIN WAVES IN TOKAMAKS* A. Bers, A. K. Ram, and S. D. Schultz Plasma Science and Fusion Center, Massachusetts Institute of Technology, Cambridge, MA 02139. U.S.A. Abstract The
More information1 Introduction Duality transformations have provided a useful tool for investigating many theories both in the continuum and on the lattice. The term
SWAT/102 U(1) Lattice Gauge theory and its Dual P. K. Coyle a, I. G. Halliday b and P. Suranyi c a Racah Institute of Physics, Hebrew University of Jerusalem, Jerusalem 91904, Israel. b Department ofphysics,
More informationAnumerical and analytical study of the free convection thermal boundary layer or wall jet at
REVERSED FLOW CALCULATIONS OF HIGH PRANDTL NUMBER THERMAL BOUNDARY LAYER SEPARATION 1 J. T. Ratnanather and P. G. Daniels Department of Mathematics City University London, England, EC1V HB ABSTRACT Anumerical
More informationStatistical Mechanics of Nonequilibrium Liquids
1. Introduction Mechanics provides a complete microscopic description of the state of a system. When the equations of motion are combined with initial conditions and boundary conditions, the subsequent
More informationShear Properties and Wrinkling Behaviors of Finite Sized Graphene
Shear Properties and Wrinkling Behaviors of Finite Sized Graphene Kyoungmin Min, Namjung Kim and Ravi Bhadauria May 10, 2010 Abstract In this project, we investigate the shear properties of finite sized
More informationQuiz 1 Introduction to Polymers
090109 Quiz 1 Introduction to Polymers In class we discussed the definition of a polymer first by comparing polymers with metals and ceramics and then by noting certain properties of polymers that distinguish
More informationMultiscale simulations of complex fluid rheology
Multiscale simulations of complex fluid rheology Michael P. Howard, Athanassios Z. Panagiotopoulos Department of Chemical and Biological Engineering, Princeton University Arash Nikoubashman Institute of
More informationEvan Mitsoulis. School of Mining Engineering and Metallurgy. National Technical University of Athens. R.R. Huilgol
Cessation of Couette and Poiseuille ows of a Bingham plastic and nite stopping times Maria Chatzimina, Georgios C. Georgiou Department of Mathematics and Statistics, University of Cyprus P.O. Box 2537,
More informationTransferable Potentials for Phase Equilibria. 1. United-Atom Description of n-alkanes
J. Phys. Chem. B 1998, 102, 2569-2577 2569 Transferable Potentials for Phase Equilibria. 1. United-Atom Description of n-alkanes Marcus G. Martin and J. Ilja Siepmann* Department of Chemistry, UniVersity
More informationA comparative molecular dynamics study of diffusion of n-decane and 3-methyl pentane in Y zeolite
J. Chem. Sci., Vol. 11, No. 5, September 009, pp. 91 97. Indian Academy of Sciences. A comparative molecular dynamics study of diffusion of n-decane and 3-methyl pentane in Y zeolite F G PAZZONA 1, *,
More informationMolecular dynamics study of isobaric and isochoric glass transitions in a model amorphous polymer
JOURNAL OF CHEMICAL PHYSICS VOLUME 110, NUMBER 14 8 APRIL 1999 Molecular dynamics study of isobaric and isochoric glass transitions in a model amorphous polymer Liu Yang, a) David J. Srolovitz, and Albert
More informationarxiv: v1 [cond-mat.stat-mech] 6 Jan 2014
arxiv:1401.1181v1 [cond-mat.stat-mech] 6 Jan 2014 Determination of Forces from a Potential in Molecular Dynamics (note) Bernard Monasse Mines-ParisTech, Cemef bernard.monasse@mines-paristech.fr January
More informationPotential energy surfaces for polymer, biological,
Computation of Internal Coordinates, Derivatives, and Gradient Expressions: Torsion and Improper Torsion ROBERT E. TUZUN, DONALD W. NOID, 2 BOBBY G. SUMPTER 2 Computational Science Program, SUNY Brockport,
More informationShear viscosity of model mixtures by nonequilibrium molecular dynamics. I. Argon-krypton mixtures
Shear viscosity of model mixtures by nonequilibrium molecular dynamics. I. Argon-krypton mixtures Song Hi Leea) and Peter T. Cummingsb) Department of ChemicaI Engineering, Thornton Hall, University of
More informationReactive Empirical Force Fields
Reactive Empirical Force Fields Jason Quenneville jasonq@lanl.gov X-1: Solid Mechanics, EOS and Materials Properties Applied Physics Division Los Alamos National Laboratory Timothy C. Germann, Los Alamos
More informationA Corresponding State Theory for the Viscosity of Liquids Bull. Korean Chem. Soc. 2008, Vol. 29, No Articles
A Corresponding State Theory for the Viscosity of Liquids Bull. Korean Chem. Soc. 2008, Vol. 29, No. 1 33 Articles A Corresponding State Theory for the Viscosity of Liquids Wonsoo Kim * and Sukbae Lee
More informationChain length effects on aqueous alkane solubility near the solvent s critical point
Fluid Phase Equilibria 183 184 (2001) 289 294 Chain length effects on aqueous alkane solubility near the solvent s critical point Eric M. Yezdimer a,b,, Ariel A. Chialvo c,d, Peter T. Cummings b,e a Department
More informationSimulation of Alkane Adsorption in the Aluminophosphate Molecular Sieve AlPO 4-5
J. Phys. Chem. B 1998, 102, 7183-7189 7183 Simulation of Alkane Adsorption in the Aluminophosphate Molecular Sieve AlPO 4-5 Thierry Maris, Thijs J. H. Vlugt, and Berend Smit* Department of Chemical Engineering,
More informationEffect of oscillatory shear on the fluid solid phase transition of supercooled water
Ž. Fluid Phase Equilibria 150 151 1998 Effect of oscillatory shear on the fluid solid phase transition of supercooled water Istvan Borzsak a,b,), Peter T. Cummings a,b a Department of Chemical Engineering,
More informationForce Field for Water Based on Neural Network
Force Field for Water Based on Neural Network Hao Wang Department of Chemistry, Duke University, Durham, NC 27708, USA Weitao Yang* Department of Chemistry, Duke University, Durham, NC 27708, USA Department
More information6 VORTICITY DYNAMICS 41
6 VORTICITY DYNAMICS 41 6 VORTICITY DYNAMICS As mentioned in the introduction, turbulence is rotational and characterized by large uctuations in vorticity. In this section we would like to identify some
More informationA New Intermolecular Potential Model for the n- Alkane Homologous Series
A New Intermolecular Potential Model for the n- Alkane Homologous Series Jeffrey R. Errington and Athanassios Z. Panagiotopoulos School of Chemical Engineering Cornell University Ithaca, NY 14853 and Institute
More informationInstitut fur Theoretische Physik. Universitat Heidelberg. Philosophenweg 19. Germany. and. Interdisziplinares Zentrum. fur Wissenschaftliches Rechnen
A new general model with non-spherical interactions for dense polymer systems and a potential parametrization for Bisphenol-A-Polycarbonate Klaus M. Zimmer, Andreas Linke and Dieter W. Heermann Institut
More information5 Classical Molecular Dynamics (MD) Simulations with Tinker
5 Classical Molecular Dynamics (MD) Simulations with Tinker In this set of exercises, you will be briey introduced to the concept of classical force elds and molecular dynamics via propagation of atomic
More informationNote on the KaplanYorke Dimension and Linear Transport Coefficients
Journal of Statistical Physics, Vol. 101, Nos. 12, 2000 Note on the KaplanYorke Dimension and Linear Transport Coefficients Denis J. Evans, 1 E. G. D. Cohen, 2 Debra J. Searles, 3 and F. Bonetto 4 Received
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 informationNon-equilibrium molecular dynamics calculation of the shear viscosity of liquid rubidium
J. Phys. F: Met. Phys. 17 (1987) 593-604. Printed in the UK Non-equilibrium molecular dynamics calculation of the shear viscosity of liquid rubidium P T Cummingsf and G P Morrisst + Department of Chemical
More informationExample questions for Molecular modelling (Level 4) Dr. Adrian Mulholland
Example questions for Molecular modelling (Level 4) Dr. Adrian Mulholland 1) Question. Two methods which are widely used for the optimization of molecular geometies are the Steepest descents and Newton-Raphson
More informationIntroduction. Theoretical Model and Calculation of Thermal Conductivity
[P2.23] Molecular dynamic study on thermal conductivity of methyl-chemisorption carbon nanotubes H. He, S. Song*, T. Tang, L. Liu Qingdao University of Science & Technology, China Abstract. The thermal
More informationRheology, Adhesion, and Debonding of Lightly Cross-linked Polymer Gels
Rheology, Adhesion, and Debonding of Lightly Cross-linked Polymer Gels Nicholas B. Wyatt, and Anne M. Grillet 2 Materials Science and Engineering Division 2 Engineering Sciences Division Sandia National
More informationAbstract. Because mantle viscosity is temperature-dependent, cold subducting
Chapter 6 Mantle Convection with Strong Subduction Zones Abstract. Because mantle viscosity is temperature-dependent, cold subducting lithosphere should be strong, which implies that the rapid, localized
More informationLes Houches School of Foam: Rheology of Complex Fluids
Les Houches School of Foam: Rheology of Complex Fluids Andrew Belmonte The W. G. Pritchard Laboratories Department of Mathematics, Penn State University 1 Fluid Dynamics (tossing a coin) Les Houches Winter
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 informationCoupling of Organic and Inorganic Vibrational States and. Their Thermal Transport in Nanocrystal Arrays
Supporting Information: Coupling of Organic and Inorganic Vibrational States and Their Thermal Transport in Nanocrystal Arrays Wee-Liat Ong, Shubhaditya Majumdar, Jonathan A. Malen, Alan J. H. McGaughey
More informationOrigin of Particle Clustering in a Simulated Polymer Nanocomposite and its Impact on Rheology
Wesleyan University From the SelectedWorks of Francis Starr 2003 Origin of Particle Clustering in a Simulated Polymer Nanocomposite and its Impact on Rheology Francis W. Starr, Wesleyan University J. F.
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 informationHeat and Mass Transfer over Cooled Horizontal Tubes 333 x-component of the velocity: y2 u = g sin x y : (4) r 2 The y-component of the velocity eld is
Scientia Iranica, Vol., No. 4, pp 332{338 c Sharif University of Technology, October 2004 Vapor Absorption into Liquid Films Flowing over a Column of Cooled Horizontal Tubes B. Farhanieh and F. Babadi
More informationA Study of the Thermal Properties of a One. Dimensional Lennard-Jones System
A Study of the Thermal Properties of a One Dimensional Lennard-Jones System Abstract In this study, the behavior of a one dimensional (1D) Lennard-Jones (LJ) system is simulated. As part of this research,
More informationJournal of Computational Physics
Journal of Computational Physics 231 (2012) 5306 5316 Contents lists available at SciVerse ScienceDirect Journal of Computational Physics journal homepage: www.elsevier.com/locate/jcp Nonequilibrium molecular
More informationEquation of state for hard n-alkane models: Long chains
Equation of state for hard n-alkane models: Long chains C. Vega and L. G. MacDowell Departamento de Quimica Fisica, Facultad de Ciencias Quimicas, Universidad Complutense, 28040 Madrid, Spain P. Padilla
More informationMathematics Research Report No. MRR 003{96, HIGH RESOLUTION POTENTIAL FLOW METHODS IN OIL EXPLORATION Stephen Roberts 1 and Stephan Matthai 2 3rd Febr
HIGH RESOLUTION POTENTIAL FLOW METHODS IN OIL EXPLORATION Stephen Roberts and Stephan Matthai Mathematics Research Report No. MRR 003{96, Mathematics Research Report No. MRR 003{96, HIGH RESOLUTION POTENTIAL
More information2 Dissipative particle dynamics with energy conservation: Heat conduction that is dissipated due to the velocity dependent forces is transformed into
International Journal of Modern Physics C, fc World Scientic Publishing Company DISSIPATIVE PARTICLE DYNAMICS WITH ENERGY CONSERVATION: HEAT CONDUCTION M. RIPOLL, P. ESPA ~NOL Departamento de Fsica Fundamental,
More informationRHEOLOGY Principles, Measurements, and Applications. Christopher W. Macosko
RHEOLOGY Principles, Measurements, and Applications I -56081-5'79~5 1994 VCH Publishers. Inc. New York Part I. CONSTITUTIVE RELATIONS 1 1 l Elastic Solid 5 1.1 Introduction 5 1.2 The Stress Tensor 8 1.2.1
More informationORGANIC - BRUICE 8E CH.3 - AN INTRODUCTION TO ORGANIC COMPOUNDS
!! www.clutchprep.com CONCEPT: INDEX OF HYDROGEN DEFICIENCY (STRUCTURAL) A saturated molecule is any molecule that has the maximum number of hydrogens possible for its chemical structure. The rule that
More informationWATER PERMEATION THROUGH GRAPHENE NANOSLIT BY MOLECULAR DYNAMICS SIMULATION
WATER PERMEATION THROUGH GRAPHENE NANOSLIT BY MOLECULAR DYNAMICS SIMULATION Taro Yamada 1 and Ryosuke Matsuzaki 2 1 Department of Mechanical Engineering, Tokyo University of Science, 2641 Yamazaki, Noda,
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 information