Modified Nosé-Hoover thermostat for solid state for constant temperature molecular dynamics simulation

Nose-Hoover (NH) thermostat methods incorporated with molecular dynamics (MD) simulation have been widely used to simulate the instantaneous system temperature and feedback energy in a canonical ensemble. The method simply relates the kinetic energy to the system temperature via the particles' momenta based on the ideal gas law. However, when used in a tightly bound system such as solids, the method may suffer from deriving a lower system temperature and potentially inducing early breaking of atomic bonds at relatively high temperature due to the neglect of the effect of the potential energy of atoms based on solid state physics. In this paper, a modified NH thermostat method is proposed for solid system. The method takes into account the contribution of phonons by virtue of the vibrational energy of lattice and the zero-point energy, derived based on the Debye theory. Proof of the equivalence of the method and the canonical ensemble is first made. The modified NH thermostat is tested on different gold nanocrystals to characterize their melting point and constant volume specific heat, and also their size and temperature dependence. Results show that the modified NH method can give much more comparable results to both the literature experimental and theoretical data than the standard NH. Most importantly, the present model is the only one, among the six thermostat algorithms under comparison, that can accurately reproduce the experimental data and also the T^3-law at temperature below the Debye temperature, where the specific heat of a solid at constant volume is proportional to the cube of temperature.

[1]  T. H. Geballe,et al.  The Heat Capacity and Entropy of Gold from 15 to 300°K.1 , 1952 .

[2]  Federico Frascoli,et al.  Molecular dynamics simulation of planar elongational flow at constant pressure and constant temperature. , 2007, The Journal of chemical physics.

[3]  Jonathan W Essex,et al.  Conformational Motions of HIV-1 Protease Identified Using Reversible Digitally Filtered Molecular Dynamics. , 2009, Journal of chemical theory and computation.

[4]  M. Tuckerman,et al.  On the classical statistical mechanics of non-Hamiltonian systems , 1999 .

[5]  Denis J. Evans,et al.  The Nose–Hoover thermostat , 1985 .

[6]  James F. Shackelford,et al.  The CRC Materials Science And Engineering Handbook , 1991 .

[7]  G. Ciccotti,et al.  Non-Hamiltonian molecular dynamics: Generalizing Hamiltonian phase space principles to non-Hamiltonian systems , 2001 .

[8]  Young Whan Cho,et al.  Thermal stability of unsupported gold nanoparticle: a molecular dynamics study , 2002 .

[9]  M. Tuckerman,et al.  Generalized Gaussian moment thermostatting: A new continuous dynamical approach to the canonical ensemble , 2000 .

[10]  Kim,et al.  Melting the fullerenes: A molecular dynamics study. , 1994, Physical review letters.

[11]  B. Berne,et al.  Efficient molecular dynamics and hybrid Monte Carlo algorithms for path integrals , 1993 .

[12]  L. Girifalco,et al.  Application of the Morse Potential Function to Cubic Metals , 1959 .

[13]  Glenn J. Martyna,et al.  Molecular dynamics simulations of a protein in the canonical ensemble , 1993 .

[14]  J. Gillis,et al.  Classical dynamics of particles and systems , 1965 .

[15]  Michele Parrinello,et al.  Integrating the Car–Parrinello equations. I. Basic integration techniques , 1994 .

[16]  Alessandro Laio,et al.  A Collective Variable for the Efficient Exploration of Protein Beta-Sheet Structures: Application to SH3 and GB1. , 2009, Journal of chemical theory and computation.

[17]  M. Klein,et al.  Nosé-Hoover chains : the canonical ensemble via continuous dynamics , 1992 .

[18]  H. C. Andersen.  Molecular dynamics simulations at constant pressure and/or temperature , 1980 .

[19]  Susan B. Sinnott,et al.  Constant temperature molecular dynamics simulations of energetic particle-solid collisions: comparison of temperature control methods , 2004 .

[20]  P. Buffat,et al.  Size effect on the melting temperature of gold particles , 1976 .

[21]  Victor Szebehely,et al.  Time transformations in the extended phase-space , 1975 .

[22]  Linda S. Schadler,et al.  Temperature dependence of radial breathing mode Raman frequency of single-walled carbon nanotubes , 2002 .

[23]  Eun-Ha Kim,et al.  Size dependency of melting point of crystalline nano particles and nano wires: A thermodynamic modeling , 2009 .

[24]  K. Travis,et al.  A configurational temperature Nosé-Hoover thermostat. , 2005, The Journal of chemical physics.

[25]  H. Berendsen,et al.  Molecular dynamics with coupling to an external bath , 1984 .

[26]  E. Lomba,et al.  Hubbard corrections in a tight-binding Hamiltonian for Se: Effects on the band structure, local order, and dynamics , 2000 .

[27]  J. D. Doll,et al.  Generalized Langevin equation approach for atom/solid-surface scattering: General formulation for classical scattering off harmonic solids , 1976 .

[28]  HengAn Wu,et al.  Molecular dynamics study on mechanics of metal nanowire , 2006 .

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

[30]  Akira Takada,et al.  Molecular dynamics simulations of vitreous silica structures , 2004 .

[31]  L. V. Woodcock.  Isothermal molecular dynamics calculations for liquid salts , 1971 .

[32]  I László.  Formation of cage-like C60 clusters in molecular-dynamics simulations , 1998 .

[33]  Thomas Gilbert.  Fluctuation theorem applied to the Nosé-Hoover thermostated Lorentz gas. , 2006, Physical review. E, Statistical, nonlinear, and soft matter physics.

[34]  John S. Rowlinson,et al.  Molecular Theory of Capillarity , 1983 .

[35]  Subramanian K. R. S. Sankaranarayanan,et al.  Molecular dynamics simulation study of the melting and structural evolution of bimetallic Pd-Pt nanowires , 2006 .

[36]  HengAn Wu,et al.  Molecular dynamics study of the mechanics of metal nanowires at finite temperature , 2006 .

[37]  Z. Kovács,et al.  Fluctuation formula in the Nosé-Hoover thermostated Lorentz gas. , 2004, Physical review. E, Statistical, nonlinear, and soft matter physics.

[38]  S. Nosé.  A molecular dynamics method for simulations in the canonical ensemble , 1984 .

[39]  J. M. Haile,et al.  Molecular dynamics simulation : elementary methods / J.M. Haile , 1992 .

[40]  Ko Higashitani,et al.  Revised Implicit Solvent Model for the Simulation of Surfactants in Aqueous Solutions. 2. Modeling of Charged Headgroups at Oil-Water Interface. , 2007, Journal of chemical theory and computation.

[41]  Koichiro Nakanishi,et al.  Constant temperature molecular dynamics calculation on Lennard‐Jones fluid and its application to watera) , 1983 .

[42]  Stephen D. Bond,et al.  The Nosé-Poincaré Method for Constant Temperature Molecular Dynamics , 1999 .

[43]  Shozaburo Saito,et al.  Molecular dynamics study on solvent clustering in supercritical fluid solutions based on particle radial kinetic energy , 1995 .

[44]  Brad A. Bauer,et al.  Incorporating Phase-Dependent Polarizability in Non-Additive Electrostatic Models for Molecular Dynamics Simulations of the Aqueous Liquid-Vapor Interface. , 2009, Journal of chemical theory and computation.

[45]  J. Haile.  Molecular Dynamics Simulation , 1992 .

[46]  C. Kittel.  Introduction to solid state physics , 1954 .

[47]  Farid F. Abraham,et al.  Computer-Simulation Dynamics of an Unstable Two-Dimensional Fluid: Time-Dependent Morphology and Scaling , 1982 .

[48]  Akira Onuki,et al.  Molecular dynamics simulation of heat conduction in near-critical fluids. , 2005, Physical review. E, Statistical, nonlinear, and soft matter physics.