Sharif Digital Repository

Water inside carbon nanotubes is an interesting confined system, and its theoretical and applied aspects have been studied extensively. The confinement in nanometer- and sub-nanometer-sized nanotubes gives rise to new molecular arrangements of water and affects its physical properties drastically. In order to study these new arrangements, Monte Carlo simulations of water inside carbon nanotubes have been performed. Simulations are carried out for water with a wide range of density inside carbon nanotubes with different diameters. It is observed that at constant temperature, the density of water dictates the presence of water clusters, filled states, and different ordered phases. The... 

A general regularity was found based on an effective pair potential of Lennard-Jones LJ (12, 6), for both dense, nonmetallic and nonionic fluids and solids according to which (Z-1) v2 linearly varies with respect to ρ2 for each isotherm, and this equation of state (EoS I) is known as LIR. However, despite the fact that Ne is a simple spherical species, unexpectedly, its solid and liquid phases both show a significant deviation from EoS I. In this work, we have investigated the accuracy of the EoS I for other systems, including quantum light molecules, such as D 2, H 2 and He, in both fluid and solid states at different temperatures. Like Ne, we have noticed that these systems do not well... 

Water inside carbon nanotubes as an example of nanoconfined water has gained noticeable attention, in both theoretical and applied aspects. Molecular simulation has played a major role in the studies in this field. Yet, there is a need for systematic study of simulation results and compilation of scientifically reliable predictive relations. Here we present Monte Carlo simulations of water inside carbon nanotubes with different radii. An equation of state which was derived on the basis of the extended Lennard-Jones (12,6,3) as the effective pair potential is chosen for the system of water inside the carbon nanotubes. The equation of state is modified to take the effects of anisotropic... 

Metallic nanoclusters are interesting because of their utility in catalysis and sensors. The thermal and physical characteristics of metallic Pt nanoclusters with different sizes were investigated via molecular-dynamics simulations using Quantum Sutton-Chen (QSC) potential. This force field accurately predicts solid and liquid states properties as well as melting of the bulk platinum. Molecular dynamic simulations of Pt nanoclusters with 256, 456, 500, 864, 1372, 2048, 2916, 4000, 5324, 6912, 8788 atoms have been carried out at various temperatures. The Pt-Pt radial distribution function, internal energy, heat capacity, enthalpy, entropy of the nanoclusters were calculated at some... 

We have derived an analytical equation of state (EOS) based on the soft-core statistical mechanical perturbation theory for fluids, using the Weeks-Chandler-Andersen (WCA) theory recently developed by Ben-Amotz-Stell (BAS) for the choice of the hard-sphere diameter, but with a new algorithm for calculation of the pair and many-body interactions. We have used Carnahan-Starling expression with the Boltzmann factor criterion (BFC) as an effective hard-sphere diameter for the reference system, and also decomposed the perturbed pair potential to symmetric and asymmetric terms. The former term is due to the many-body interactions at high densities as was used in the linear isotherm regularity... 

The effect of repulsive steepness of the soft-core square well (SCSW) potential model on the second virial coefficient, critical behaviour (two- phase region and the position of critical point), and coordination number are investigated. The soft-core thermodynamic perturbation theory (TPT) presented by Weeks-Chandler-Anderson (WCA) recently developed by Ben-Amotz and Stell (BAS) has been used for the reference system, and the Barker-Henderson TPT for the perturbed system. The Barker-Henderson macroscopic compressibility approximation has been used for all order perturbation terms in which the second-order one is improved by assuming that the molecules in every two neighbouring shells are... 

The hydrogen abstraction reaction of the OH radical with CH 2FCH22F (HFC-152) is studied theoretically over the 150-3000 K temperature range. In this study, the two most recently developed hybrid density functional theories, namely, BBlK and MPWBlK, are applied, and their efficiency in reaction dynamics calculation is discussed. The BBlK/6-31+G(d,p) method gives the best result for the potential energy surface (PES) calculations, including barrier heights, reaction path information (the first and second derivatives of PES), geometry of transition state structures, and even weak hydrogen bond orientations. The rate constants were obtained by the dual-level direct dynamics with the... 

We present in this paper a simple method of obtaining various equations of state for hard sphere fluid in a simple unifying way. Using the first several virial coefficients of hard sphere fluid, we will guess equations of state by using the asymptotic expansion method. Among the equations of state obtained in this way are Percus-Yevick, Scaled Particle Theory, and Carnahan-Starling equations of state. Also by combining the Monte Carlo results on hard sphere fluid with the asymptotic expansion method many other equations of state for hard sphere fluid can be found where all of them give essentially similar results in the region of isotropic hard sphere liquid, i.e., up to η < 0.5, in which η... 

Using the newly introduced asymptotic expansion method to obtain equations of state for hard sphere fluid, new simple equations of state for hard disks based on known virial coefficients are derived. Comparison of the obtained equations of state with computer simulation data shows that they are accurate in the whole fluid region. © 2007 Elsevier B.V. All rights reserved  

The interaction potential energy and heat of sublimation of nanoparticles of HMX crystal polymorphs are studied by using molecular dynamics methods with a previously developed force field [Bedrov, et al., J. Comput.-Aided Mol. Des. 8 (2001) 77]. Molecular dynamics simulations of nanoparticles with 10, 20, 30, 40, 50, 60, 70, 80, 90, and 100 molecules of HMX are carried out at 300 K. The intermolecular, intramolecular and total interaction energies per mole for the nanoparticles are calculated at 300 K. Then, we have calculated sublimation enthalpy of HMX crystal polymorphs with different sizes. For the all sizes, the β-HMX is found to be the most stable phase, due to having the least total... 

In this paper the density and temperature dependencies of surface tension are investigated. Using the Lennard-Jones (12,6), as an effective pair interaction potential, a linear expression is derived for isotherms of γ /ρ 2 versus ρ 2 for some normal and ChloroFluoroCarbons (CFCs) fluids, where is surface tension and ρ = 1/v is molar density. The linearity behavior of the derived equation is well fitted onto the experimental data of surface tension for monatomic, diatomic, nonpolar, polar, hydrogen-bonded and chlorofluorocarbons. In addition, the temperature dependence of surface tension for 20 different fluids is examined, in which the contributions of both terms of the average effective... 

Using a general expression for the second term in the density expansion of the hard core potentials, we obtained an analytical expression for the direct correlation function of the hard core double Yukawa potential. The results of this calculation using the proposed direct correlation function and its related radial distribution function through the Ornstein-Zernike equation show good agreements with those obtained by the computer simulation of the same model  

As the delay time and hence nuclei formation play a crucial role in the pathophysiology of sickle cell disease, MD simulation and molecular mechanics-Poisson-Boltzmann surface area (MM-PBSA) calculations have been performed on three systems of hemoglobin; namely dimer of hemoglobin with valine (Hb S), tryptophan (Hbβ6W), and phenylalanine (Hbβ6F) at β6 position. The structural changes due to these aromatic substitutions are investigated. It is shown that β subunits have significant impact on the differences between a dimer and other crystal structures. Transition from a dimer to polymer for Hb S system affects the donor molecule more than that of the acceptor. In the case of donor and... 

Metallic nanoparticles are interesting because of their use in catalysis and sensors. The surface energy of the FCC platinum nanoparticles are investigated via molecular dynamics simulation using Quantum Sutton-Chen (QSC) potential. We have calculated the Gibbs free energy for the FCC platinum bulk and also for its nanoparticle. All calculations have been carried out at zero pressure. We have used the thermodynamic integration method to obtain the Gibbs free energy. The total Gibbs free energy is taken as the sum of its central bulk and its surface free energy. We have calculated the free energy of a platinum nanoparticle as a function of temperature  

In the present work, the thermodynamic of dense fluids, both compressed liquids and dense supercritical fluids, has been modeled, solely, based on the contribution of attraction of effective pair potential. The intermolecular interaction is modeled by the hard-core Yukawa potential (HCY) as an effective pair potential (EPP) with temperature dependent hard-core diameter. Using this EPP in the exact thermodynamic relations, an equation of state (EoS) for the compressibility factor of dense fluid has been derived. This EoS shows that (Z - ZCS) as function of ρ1/3 must be linear for each isotherm of fluid where ZCS is the compressibility factor of the reference fluid (Carnahan-Starling (CS) EoS)... 

In the present work, several new thermodynamic linear isotherm regularities for the dense fluids have been derived for the first time. For this purpose, the thermodynamic perturbation theory (TPT) employing only the attractive effective pair potential (AEPP) as u(r) = - ϵeff (σeff / r)m was used, where σeff which is the effective hard core diameter, is temperature dependent and m > 0. Based on the derived regularities, the isotherm (Z - Z(0))v2 is a linear function of ρ2, ρ and 1/ρ, depends on the values of m = 12, 9 and 6, respectively where Z - Z(0) is the difference between the experimental compressibility factor of the real fluid (Z) and that of the reference fluid (Z(0)). Also, Z - Z(0)... 

It is of interest to explore transfer of fluid through nanopores because of widespread applications for such systems. Carbon Nanotubes (CNTs) with their exceptional properties are the best candidates as building blocks for nanostructures. Water transfer in lots of biological systems acts as an important role for keeping the tissue working properly. Peptide nanotube is one of the best biological channels which was proposed recently. While the mechanism of water permeation through channels is very complex, however, investigations such as effect of charge distributions and temperature on water permeation could shed light on the determinants of water and proton conduction rates in biol ogical... 

Molecular dynamics, MD, simulation of calcite (CaCO3) is selected to compare the p-v-T behaviour of some universal equations of state, UEOS, for the temperature range 100 K ≤ T 800 ≤ K, and pressures up to 3000 kbar. The isothermal sets of p-v-T data generated by simulation were each fitted onto some three- and two-parameter EOSs including Parsafar and Mason (PM), Linear Isotherm Regularity (LIR), Birch-Murnaghan (BM), Shanker, Vinet, Baonza and Modified generalized Lennard-Jones (MGLJ) EOSs. It is found that the MD data satisfactorily fit these UEOS with reasonable precision. Some features for a good UEOS criteria such as temperature dependencies of coefficients, pressure deviation,... 

In this paper, we have derived two equations of state, one for the metallic and ionic solids and the other for the remaining solids on the basis of the concept of the average effective pair potential (AEPP). According to the former EOS, (Z-1)υ2 is linear with respect to 1/ρ, where Z is the compressibility factor and υ = 1/ρ is the molar volume for each isotherm. On the basis of the latter EOS, (Z-1)υ2 is a linear function in terms of ρ2 for each isotherm. As these EOSs suggest, the temperature dependencies of the internal energy is separable from its density dependencies. Hence, the heat capacity of a solid is independent of its density, interaction potential parameters and non-ideal thermal... 

Based on a new derived radial distribution function (RDF) for potentials with a hard-core we have presented in this paper a method to apply the derived RDF for calculating thermodynamic properties of real fluids up to moderate densities. In order to use the derived RDF for real fluids, one of the potential parameters is chosen in such a way that the RDF behaves more like that for a real fluid. Hence we have been able to calculate all thermodynamic properties of a simple fluid analytically. We have then applied our procedure to a Lennard-Jones fluid and compared the results with simulation data. The agreement is good up to moderate densities, i.e. ρ* ≤ 0.6, which lies in the liquid range of...