Grand canonical molecular dynamics for TIP4P water systems

An algorithm was developed enabling implementation of a Nose—Hoover thermostat within the framework of grand canonical molecular dynamics [Lynch, C. G. and Pettitt, B. M., 1997, J. chem. Phys., 107, 8594]. The proposed algorithm could readily be extended to mixtures of molecular species with different chemical potentials as shown in the paper. This algorithm was first applied to simulate a μVT ensemble of TIP4P water molecules at 298 K by means of a system comprising a number of full particles and a single scaled (fractional) particle, with the scaling factor considered as a dynamic variable in its own right and chemical potential a pre-set parameter. Our finding showed that the scheme with a single fractional particle tended to freeze in metastable states as well as failed to reproduce either the real-life (−24.05 kJmol−1) or the model-specific chemical potential of water (−23.0kJ mol−1). In order to overcome the inadequacy of a single fractional particle as a chemical potential ‘probe’ the treatment of ...

[1]  W. L. Jorgensen.  Revised TIPS for simulations of liquid water and aqueous solutions , 1982 .

[2]  Grant S. Heffelfinger,et al.  Diffusion in Lennard-Jones Fluids Using Dual Control Volume Grand Canonical Molecular Dynamics Simulation (DCV-GCMD) , 1994 .

[3]  Søren Toxvaerd,et al.  Algorithms for canonical molecular dynamics simulations , 1991 .

[4]  B. Montgomery Pettitt,et al.  Grand molecular dynamics: A method for open systems , 1991 .

[5]  Shuichi Nosé,et al.  Constant Temperature Molecular Dynamics Methods , 1991 .

[6]  Bernard Pettitt,et al.  Grand canonical ensemble molecular dynamics simulations: Reformulation of extended system dynamics approaches , 1997 .

[7]  Sanat K. Kumar.  A modified real particle method for the calculation of the chemical potentials of molecular systems , 1992 .

[8]  B. Montgomery Pettitt,et al.  Ideal chemical potential contribution in molecular dynamics simulations of the grand canonical ensemble , 1994 .

[9]  J. Banavar,et al.  Computer Simulation of Liquids , 1988 .

[10]  B. Widom,et al.  Some Topics in the Theory of Fluids , 1963 .

[11]  Hideki Tanaka,et al.  Molecular design of polymer membranes using molecular simulation technique , 1995 .

[12]  David Fincham,et al.  Leapfrog Rotational Algorithms , 1992 .

[13]  A. Pohorille,et al.  Excess chemical potential of small solutes across water--membrane and water--hexane interfaces. , 1996, The Journal of chemical physics.

[14]  D. Heyes,et al.  Molecular dynamics simulations of stretched water: Local structure and spectral signatures , 1998 .

[15]  S. Furukawa,et al.  Non-Equilibrium Molecular Dynamics for Simulating Permeation of Gas Mixtures through Nanoporous Carbon Membranes , 1996 .

[16]  Cracknell,et al.  Direct molecular dynamics simulation of flow down a chemical potential gradient in a slit-shaped micropore. , 1995, Physical review letters.

[17]  Molecular dynamics implementation of the Gibbs ensemble calculation , 1994 .

[18]  Hoover,et al.  Canonical dynamics: Equilibrium phase-space distributions. , 1985, Physical review. A, General physics.

[19]  B. Kvamme,et al.  Ergodicity Range of Nosé-Hoover Thermostat Parameters and Entropy-Related Properties of Model Water Systems , 1999 .

[20]  B. Widom,et al.  Potential-distribution theory and the statistical mechanics of fluids , 1982 .

[21]  J. Hermans,et al.  Excess free energy of liquids from molecular dynamics simulations. Application to water models. , 1988, Journal of the American Chemical Society.