Molecular Dynamics Simulation of Methane Hydrate Dissociation Process in the Presence of Thermodynamic Inhibitor

The dissociation of methane hydrate in the presence of ethylene glycol (11.45 mol.L-1) at 277.0 K was studied using canonical ensemble (NVT) molecular dynamics simulations. Results show that hydrate dissociation starts from the surface layer of the solid hydrate and then gradually expands to the internal layer. Thus, the solid structure gradually shrinks until it disappears. A distortion of the hydrate lattice structure occurs first and then the hydrate evolves from a fractured frame to a fractional fragment. Finally, water molecules in the hydrate construction exist in the liquid state. The inner dissociating layer is, additionally, coated by a liquid film formed from outer dissociated water molecules outside. This film inhibits the mass transfer performance of the inner molecules during the hydrate dissociation process.

Molecular dynamics simulation of entropy and surface tension for grain boundary of α-Fe

The grain boundary is an interface and the surface tension is one of its important thermodynamic properties. In this paper, the surface tension of the ∑9 grain boundary for α-Fe at various temperatures and pressures is calculated by means of Computer Molecular Dynamics (CMD). The results agree satisfactorily with the experimental data. It is shown that the contribution of entropy to surface tension of grain boundary can be ignored.

Phase transition temperatures of side chain liquid crystal polymer from molecular dynamics simulation

In this paper, phase transition temperatures of side chain liquid crystal polymer were predicted by molecular dynamics simulation. We analyzed the change of energy and the degree of similarity(S) with the temperature varying. The simulated phase transition temperatures agree with the experimental values in a proper deviation.

Molecular dynamics simulation of polystyrene-block-poly(methyl methacrylate)

Molecular dynamics is applied to the system of polystyrene-block-poly(methyl methacrylate). The simulation shows that for the block copolymer system, a layered structure, which reflects microphase separation, is obtained and this structure is stable. In order to elucidate that the formation of the layered structure is reasonable, some static properties such as the radial distribution function and the dipole moment are analyzed in some detail.

Molecular dynamics simulation of observed c(4x4) and c(4x3) C60 alignments on the Si(100) reconstructed surface.

We have implemented a large-scale classical molecular dynamics simulation at constant temperature to provide a theoretical insight into the results of a recently performed experiment on the monolayer and multi-layer formations of molecular films on the Si(100) reconstructed dimerized surface. Our simulation has successfully reproduced all of the morphologies observed on the monolayer film by this experiment. We have obtained the formation of both c(4 4) and c(4 3) structures of the molecules and have also obtained phase transitions of the former into the latter.

Molecular dynamics simulation of adsorption of Ag particles on a graphite substrate

We have performed for the first time a molecular dynamics simulation of the adsorption of gas-phase Ag particles on a graphite substrate to provide an insight into the results of a comprehensive STM-based experiment on this system. Both pair-wise and many-body interatomic potentials have been employed, and a Morse-type Ag–C potential was specifically constructed to describe the interactions at the interface. Our simulation has successfully reproduced a significant portion of the experimental findings. We have also observed the intercalation of silver in graphite.

Molecular dynamics simulation of fractures using an N-body potential

Molecular dynamics has been employed to model the fracture of a twodimensional triangular atomic lattice. The N-body Sutton-Chen potential developed for fcc metals and its extended version (Rafii-Tabar and Sutton) for fcc random binary alloys were used for the interatomic interactions. It is shown that at low temperatures cleavage fractures can occur in both an elemental metal and an alloy. At elevated temperatures the nucleation of dislocations is shown to cause a brittle-to-ductile transition. For the brittle crack propagation in the elemental metal, crack propagation speeds have been computed for different stress rates, and a crack instability found to exist as the speed reaches a critical value of about 32% of the Rayleigh wave speed. For the random alloy, we find that the dislocation movement can be affected by the distorted lattice.

Molecular dynamics simulation of fractures using an N-body potential

Molecular dynamics has been employed to model the fracture of a two dimensional triangular atomic lattice. The N-body Sutton-Chen potential developed for fcc metals and its extended version (Rafii-Tabar and Sutton) for fcc random binary alloys were used for the interatomic interactions. It is shown that at low temperatures cleavage fractures can occur in both an elemental metal and an alloy. At elevated temperatures the nucleation of dislocations is shown to cause a brittle-to-ductile transition. For the brittle crack propagation in the elemental metal, crack propagation speeds have been computed for different stress rates, and a crack instability found to exist as the speed reaches a critical value of about 32% of the Rayleigh wave speed. For the random alloy, we find that the dislocation movement can be affected by the distorted lattice.

Numerical modelling of adsorption of metallic particles on graphite substrate via molecular dynamics simulation

A computer-based numerical modelling of the adsorption process of gas phase metallic particles on the surface of a graphite substrate has been performed via the application of molecular dynamics simulation method. The simulation relates to an extensive STM-based experiment performed in this field, and reproduces part of the experimental results. Both two-body and many-body inter-atomic potentials have been employed. A Morse-type potential describing the metal-carbon interactions at the interface was specifically formulated for this modelling. Intercalation of silver in graphite has been observed as well as the correct alignments of monomers, dimers and two-dimensional islands on the surface. PACS numbers: 02.60.Cb, 07.05.Tp, 68.55.-a, 81.05.Tp

Ab initio molecular dynamics simulation of a room temperature ionic liquid

Ab initio molecular dynamics simulations have been performed for the first time on the room-temperature organic ionic liquid dimethyl imidazolium chloride [DMIM][Cl] using density functional theory. The aim is to compare the local liquid structure with both that obtained from two different classical force fields and from neutron scattering experiments. The local structure around the cation shows significant differences compared to both the classical calculations and the neutron results. In particular, and unlike in the gas-phase ion pair, chloride ions tend to be located near a ring C-H proton in a position suggesting hydrogen bonding. The results are used to suggest ways in which the classical potentials may be improved.

Molecular dynamics simulation of ionic liquids: The effect of electronic polarizability

An electronically polarizable model has been developed for the ionic liquid 1-ethyl-3-methylimidazolium nitrate (EMIM+/NO3-), Molecular dynamics simulation studies were then performed on both the polarizable and nonpolarizable versions of the model. Comparisons of shear viscosity and diffusion constants at 400 K show that the effects of polarizability are quite substantial and the polarizable model results are in better agreement with the experimental values.

Binding modes of diketo-acid inhibitors of HIV-1 integrase: A comparative molecular dynamics simulation study

HIV-1 integrase (IN) has become an attractive target since drug resistance against HIV-1 reverse transcriptase (RT) and protease (PR) has appeared. Diketo acid (DKA) inhibitors are potent and selective inhibitors of HIV-1 IN: however the action mechanism is not well understood. Here, to study the inhibition mechanism of DKAs we performed 10 ns comparative molecular dynamics simulations on HIV-1 IN bound with three most representative DMA inhibitors: Shionogi inhibitor, S-1360 and two Merck inhibitors L-731,988 and L-708,906. Our simulations show that the acidic part of S-1360 formed salt bridge and cation-pi interactions with Lys159. In addition, the catalytic Glu152 in S-1360 was pushed away from the active site to form an ion-pair interaction with Arg199. The Merck inhibitors can maintain either one or both of these ion-pair interaction features. The difference in potencies of the DMA inhibitors is thus attributed to the different binding modes at the catalytic site. Such structural information at atomic level, not only demonstrates the action modes of DMA inhibitors but also provides a novel starting point for structural-based design of HIV-1 IN inhibitors.