Sharif Digital Repository

We have investigated the size dependence of a nano-cavity properties produced in a Xe fluid using molecular dynamics simulations. We have created a nano-cavity of different sizes at 170 and 200 K (cavities diameters are within 1-10 nm). Liquid pressure, vapor pressure and surface tension of the nano-cavity for some given values of diameter are calculated. Within 1-10 nm cavity diameter, we have observed two opposite behavior for dependency of surface tension on the cavity diameter: for the range of 1-5 nm, it increases with the diameter, while, for the range of 5-10 nm remains constant. Also, the value of liquid pressure becomes less negative, when the size of cavity increases. Vapor... 

The concept of effective pair potential (EPP), in which medium effects and contribution of long-range interactions are added to near-neighbor pair interactions, is a useful tool to derive and calculate the thermodynamic quantities. It has been shown that the extended Lennard-Jones (12, 6, 3) potential is an accurate EPP for many systems. However, the parameters of this EPP are not known for any systems. We introduce a new approach to use experimental pvT data of a fluid to calculate its EPP parameters. Such calculation has been carried out for some fluids at some temperatures. Unlike the isolated pair potential, it has been shown that the parameters significantly vary with temperature. We... 

We have studied the topology of the energy landscape of a spin-glass model and the effect of frustration on it by looking at the connectivity and disconnectivity graphs of the inherent structure. The connectivity network shows the adjacency of energy minima whereas the disconnectivity network tells us about the heights of the energy barriers. Both graphs are constructed by the exact enumeration of a two-dimensional square lattice of a frustrated spin glass with nearest-neighbor interactions up to the size of 27 spins. The enumeration of the energy-landscape minima as well as the analytical mean-field approximation show that these minima have a Gaussian distribution, and the connectivity... 

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... 

We derive exact expressions for the forces and torques between biaxial molecules interacting via the RE-squared potential, a recent variant of the Gay-Berne potential. Moreover, efficient routines have been provided for rigid body MD simulations, resulting in 1.6 times speedup compared to the two-point finite difference approach. It has also been shown that the time cost of a MD simulation will be almost equal to a similar MC simulation, making use of the provided routines. © 2006 Elsevier Inc. All rights reserved  

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... 

Collective behavior of active matter is observed for self-propelled particles, such as vibrated disks and active Brownian particles, as well as for cytoskeletal filaments in motile cells. Here, a system of quasi two-dimensional penetrable self-propelled rods inside rigid rings is used to construct a complex self-propelled particle. The rods interact sterically with each other and with a stationary or mobile ring via a separation-shifted Lennard-Jones potential. They either have a sliding attachment to the inside of the ring at one of their ends, or can move freely within the ring confinement. We study the inner structure and dynamics of the mobile self-propelled rings. We find that these... 

Classical molecular dynamics simulations are used to compute the solvation free energy of two pharmaceutical solids, namely ibuprofen and acetaminophen in carbon dioxide (CO2), over the density range of interest in supercritical processes. In order to examine the influence of the solvent model on the resulting free energies, three popular CO2 models (Zhang, EPM2, and TraPPE) are studied. Relatively large discrepancies for the solvation free energy exist between these CO2 models, suggesting that the former is sensitive to the different balances between dispersive and electrostatic forces used in these models. In particular, for the solvation of the highly polar (dipole moment of ∼5.2 Debye)... 

Macroscale experiments show that a train of two immiscible liquid drops, a bislug, can spontaneously move in a capillary tube because of surface tension asymmetries. We use molecular dynamics simulation of Lennard-Jones fluids to demonstrate this phenomenon for NVT ensembles in submicron tubes. We deliberately tune the strength of intermolecular forces and control the velocity of bislug in different wetting and viscosity conditions. We compute the velocity profile of particles across the tube and explain the origin of deviations from the classical parabolae. We show that the self-generated molecular flow resembles the Poiseuille law when the ratio of the tube radius to its length is less... 

Natural nanochannels and those used for applied purposes are often several nanometers in diameter and have lengths up to micrometer scales. Since normally no pressure gradient, mechanical force or electrical field is exerted on fluids in these situations, the fluid flow within them is not single-file and can be bidirectional. For this reason, studying the behavior of bidirectional flows in the channels is of paramount importance. In this study, a comprehensive investigation on the straight and non-straight bidirectional water transfer through carbon nanotubes is conducted via extensive molecular dynamics simulations. The results indicate that by changing the length, the diameter and the... 

Too much effort has been done for manipulating individual atoms, using nano-manipulators and Scanning Tunneling Microscopes (STM). On the other hand, characterization and manipulation of nano-flows is of great concern. In the current work a molecular valve has been considered, which is made up of six atoms placed on the circumstance of a circle. A fuzzy controller has been designed for controlling the diameter of this molecular valve. The designed fuzzy controller used singleton fuzzifier, Mamdani inference engine, center average defuzzifier and exponential membership functions. A model based on the classical Molecular Dynamics (MD) is used for modeling the nano-system and passing the states... 

CNT-reinforced polymer composites have attracted attention due to their exceptional high strength. The high strength can be affected by the presence of defects in the nanotubes used as reinforcements in the practical nanocomposites. In this paper, a Molecular Structural Mechanics / Finite Element (MSM/FE) multiscale modeling is used to study the effect of carbon nanotube geometrical defects on the stability of SWCNT-polymer composites. Here, two types of representative volume elements (RVEs) for these nanocomposites are considered with perfect and defected CNT. These RVEs have the same dimensions. The nanotube is modeled at the atomistic scale using molecular structural mechanics whereas the... 

Shape memory alloys are referred to as a group of alloys that can retrieve the permanent deformation and strain applied to them and eventually return to their original form. So far, various studies have been done to determine the behavior of these alloys under cyclic loading. Most of the studies have mainly been conducted by using the foundations of Continuum Mechanics in order to examine the properties of memory alloys. In this study, instead of using the Continuum Mechanics, a Molecular Dynamics simulation method using Lennard-Jones potential is utilized. The changes in the behavior and properties of memory alloy under cyclic loading are being examined. First, the functional form... 

Shape memory alloys are referred to as a group of alloys that can retrieve the permanent deformation and strain applied to them and eventually return to their original form. So far, various studies have been done to determine the behavior of these alloys under cyclic loading. Most of the studies have mainly been conducted by using the foundations of Continuum Mechanics in order to examine the properties of memory alloys. In this study, instead of using the Continuum Mechanics, a Molecular Dynamics simulation method using Lennard-Jones potential is utilized. The changes in the behavior and properties of memory alloy under cyclic loading are being examined. First, the functional form... 

In this paper, using the continuum approximation together with Lennard-Jones potential, a new semi-analytical expression is given to evaluate the van der Waals interaction between two single-walled carbon nanotubes. Based on this expression, two new formulations are also proposed to model multi-walled carbon nanotubes. In the first one, the interactions between each pair of shells from the inner and outer tubes are summed up over all of the pairs, whereas in the second formulation, a set of correction factors are applied to convert the results of double-walled carbon nanotubes to the correlated multi-walled ones. With respect to the present formulations, extensive studies on the variations... 

In this article, a Molecular Structural Mechanics/Finite Element (MSM/FE) multiscale modeling of carbon nanotube/polymer composites with viscoelastic (VE) polymer matrix is introduced. The nanotube is modeled at the atomistic scale using structural molecular mechanics. The matrix deformation is analyzed by nonlinear finite element method considering VE behavior. The nanotube and matrix are assumed to be bonded by van der Waals interactions based on the Lennard-Jones potential at the interface. Using the MSM/FE multiscale model, we investigate the effect of carbon nanotube (CNT) on the improvement of mechanical stability of the nanocomposite. Also, the buckling behavior of these... 

A simple functional form for a general equation of state based on an effective near-neighbor pair interaction of an extended Lennard-Jones (12,6,3) type is given and tested against experimental data for a wide variety of fluids and solids. Computer simulation results for ionic liquids are used for further evaluation. For fluids, there appears to be no upper density limitation on the equation of state. The lower density limit for isotherms near the critical temperature is the critical density. The equation of state gives a good description of all types of fluids, nonpolar (including long-chain hydrocarbons), polar, hydrogen-bonded, and metallic, at temperatures ranging from the triple point... 

We investigate heat transfer between parallel plates separated by liquid argon using two-dimensional molecular dynamics (MD) simulations incorporating with 6-12 Lennard-Jones potential between molecule pairs. In molecular dynamics simulation of nanoscale flows through nanochannels, it is customary to fix the wall molecules. However, this approach cannot suitably model the heat transfer between the fluid molecules and wall molecules. Alternatively, we use thermal walls constructed from the oscillating molecules, which are connected to their original positions using linear spring forces. This approach is much more effective than the one which uses a fixed lattice wall modeling to simulate the... 

In this article, a multiscale modeling procedure is implemented to study the effect of interphase on the Young's modulus of CNT/polymer composites. For this purpose, a three-phase RVE is introduced which consists of three components, i.e., a carbon nanotube, an interphase layer, and an outer polymer matrix. The nanotube is modeled at the atomistic scale using molecular structural mechanics. Moreover, three-dimensional elements are employed to model the interphase layer and polymer matrix. The nanotube and polymer matrix are assumed to be bonded by van der Waals interactions based on the Lennard-Jones potential at the interface. Using this Molecular Structural Mechanics/Finite Element...