We have performed long molecular dynamics simulations of water using four popular water models, namely simple point charge ~SPC!, extended simple point charge ~SPC/E!, and the three point ~TIP3P! and four point ~TIP4P! transferable intermolecular potentials. System sizes of 216 and 820 molecules were used to study the dependence of properties on the system size. All systems were simulated at 30...

A method is proposed to compute the interfacial free energy of a Lennard-Jones system in contact with a structured wall by molecular dynamics simulation. Both the bulk liquid and bulk face-centered-cubic crystal phase along the (111) orientation are considered. Our approach is based on a thermodynamic integration scheme where first the bulk Lennard-Jones system is reversibly transformed to a st...

In molecular mechanics, current generation potential energy functions provide a reasonably good compromise between accuracy and effectiveness. This paper firstly reviewed several most commonly used classical potential energy functions and their optimization methods used for energy minimization. To minimize a potential energy function, about 95% efforts are spent on the Lennard-Jones potential o...

In this study, the behavior of a one dimensional (1D) Lennard-Jones (LJ) system is simulated. As part of this research, a special purpose computer program was designed, coded, and tested. In both the model and the program calculations, energy and momentum are conserved. The modeling assumes constant pressure throughout each simulation, although several different pressures are modeled. The simul...

using the lennard-jones (12-6) potential, a new equation of state is obtained that can predict properties of both gases and liquids relatively well. this equation of state is given as (z-a)v2=(a/v2)-b, where z is the compressibility factor, a and b are constants, and a is an adjustable parameter that depends on the temperature, volume and the nature of the fluid, and its expressions depends upo...

In this study, a temperature-dependent of the dispersion coefficients is calculated from equation state.The Lennard-Jones LJ (12-6-3) effective pair potential function and simple thermodynamic argumentwith the input PVT data of liquid metals are used to calculate the dispersion coefficients. The dispersioncoefficients ( , , ) 3 6 12 C C C are found to be a linear function of 1/T1+α , where T is...

In this paper, two finite difference numerical schemes are proposed and analyzed for the droplet liquid film model with a singular Lennard--Jones energy potential involved. Both first second or...