Brownian dynamics simulation of the crystallization dynamics of charged colloidal particles

Crystal formation process of charged colloidal particles is investigated using Brownian dynamics (BD) simulations. The particles are assumed to interact with the pair-additive repulsive Yukawa potential. The time evolution of crystallization process and the crystal structure during the simulation are characterized by means of the radial distribution functions (RDF) and mean square displacement (MSD). The simulations show that when the interaction is featured with long-range, particles can spontaneously assemble into body-centered-cubic (BCC) arrays at relatively low particle number density. When the interaction is short-ranged, with increasing the number density particles become trapped into a stagnant disordered configuration before the crystallization could be actualized. The simulations further show that as long as the trapped configurations are bypassed, the face-centered-cubic (FCC) structures can be achieved and are actually more stable than BCC structures. (C) 2010 Elsevier Inc. All rights reserved.

Molecular dynamics simulation of sub-transition for polyethersulfone

Molecular dynamics (MD) simulations of a polyethersulfone (PES) chain are carried out in the amorphous state by using the Dreiding 2.21 force field at four temperatures. Two types of molecular motion, i.e, rotations of phenylene rings and torsions of large segments containing two oxygen atoms, two sulfur atoms, and five phenylene rings on the backbone, are simulated. The modeling results show that the successive phenylene rings should be in-phase cooperative rotations, whereas the successive large segments should be out-of-phase cooperative torsions. By calculating the diffusion coefficient for the phenylene ring rotations, it is found that this rotation contributes to the beta -transition of PES.

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.