Size-dependent elastic properties of Ni nanofilms by molecular dynamics simulation

Size-dependent elastic properties of Ni nanofilms are investigated by molecular dynamics ( MD) simulations with embedded atom method (EAM). The surface effects are considered by calculating the surface relaxation, surface energy, and surface stress. The Young's modulus and yield stress are obtained as functions of thickness and crystallographic orientation. It is shown that the surface relaxation has important effects on the the elastic properties at nanoscale. When the surface relaxation is outward, the Young's modulus decreases with the film thickness decreasing, and vice versa. The results also show that the yield stresses of the films increase with the films becoming thinner. With the thickness of the nanofilms decreasing, the surface effects on the elastic properties become dominant.

Molecular dynamics simulation of adhesion forces in a dipalmitoylphosphatidylcholine membrane

Adhesion forces of Dipalmitoylphosphatidylcholine ( DPPC) membrane in the gel phase are investigated by molecular dynamics ( MD) simulation. In the simulations, individual DPPC molecules are pulled out of DPPC membranes with different rates and we get the maximum adhesion forces of DPPC membrane. We find that the maximum adhesion forces increase with pull rate, from about 400 to 700 pN when pull rates are from 0.001 to 0.03 nm/ps. We analyze the relationship between pull rate and adhesion forces of different origins using Brownian dynamics and notice that viscosity of solvent plays an important role in adhesion forces. Then we simulate the motion of a single DPPC molecule in solvent and it elucidates that the maximum drag force is almost linear with respect to the pull rate. We use Stokes' relation to describe the motion of a single DPPC molecule and deduce the effective length of a DPPC molecule. Conformational analyses indicate that the free energy variation of a DPPC molecule inside and outside of the DPPC membrane is an essential part of adhesion energy.

Shape effects on the yield stress and deformation of silicon nanowires: A molecular dynamics simulation

The tension and compression of single-crystalline silicon nanowires (SiNWs) with different cross-sectional shapes are studied systematically using molecular dynamics simulation. The shape effects on the yield stresses are characterized. For the same surface to volume ratio, the circular cross-sectional SiNWs are stronger than the square cross-sectional ones under tensile loading, but reverse happens in compressive loading. With the atoms colored by least-squares atomic local shear strain, the deformation processes reveal that the failure modes of incipient yielding are dependent on the loading directions. The SiNWs under tensile loading slip in {111} surfaces, while the compressive loading leads the SiNWs to slip in the {110} surfaces. The present results are expected to contribute to the design of the silicon devices in nanosystems.

Dissipative particle dynamics simulation of wettability alternation phenomena in the chemical flooding process

Wettability alternation phenomena is considered one of the most important enhanced oil recovery (EOR) mechanisms in the chemical flooding process and induced by the adsorption of surfactant on the rock surface. These phenomena are studied by a mesoscopic method named as dissipative particle dynamics (DPD). Both the alteration phenomena of water-wet to oil-wet and that of oil-wet to water-wet are simulated based on reasonable definition of interaction parameters between beads. The wetting hysteresis phenomenon and the process of oil-drops detachment from rock surfaces with different wettability are simulated by adding long-range external forces on the fluid particles. The simulation results show that, the oil drop is liable to spread on the oil-wetting surface and move in the form of liquid film flow, whereas it is likely to move as a whole on the water-wetting surface. There are the same phenomena occuring in wettability-alternated cases. The results also show that DPD method provides a feasible approach to the problems of seepage flow with physicochemical phenomena and can be used to study the mechanism of EOR of chemical flooding.

Structural analysis of human CCR2b and primate CCR2b by molecular modeling and molecular dynamics simulation

CCR2b, a chemokine receptor for MCP-1, -2, -3, -4, plays an important role in a variety of diseases involving infection, inflammation, and/or injury, as well as being a coreceptor for HIV-1 infection. Two models of human CCR2b (hCCR2b) were generated by h

Molecular dynamics simulation on the complex of the tobramycin and 16S rRNA a site

A 3. 6 ns molecular dynamics simulation was carried out on the complex system of tobramycin and 16S rRNA in order to understand the speciality recognition mechanism between tobramycin and 16S rRNA at the molecular level. The results demonstrate that two l

Characterize dynamic conformational space of human CCR5 extracellular domain by molecular modeling and molecular dynamics simulation

The chemokine receptor CCR5 is the receptor for several chemokines and major coreceptor for R5 human immunodeficiency virus type-1 strains entry into cell. Three-dimensional models of CCR5 were built by using homology modeling approach and 1 ns molecular dynamics (MD) simulation, because studies of site-directed mutagenesis and chimeric receptors have indicated that the N-terminus (Nt) and extracellular loops (ECLs) of CCR5 are important for ligands binding and viral fusion and entry, special attention was focused on disulfide bond function, conformational flexibility, hydrogen bonding, electrostatic interactions, and solvent-accessible surface area of Nt and ECLs of this protein part. We found that the extracellular segments of CCR5 formed a well-packet globular domain with complex interactions occurred between them in a majority of time of MID simulation, but Nt region could protrude from this domain sometimes. The disulfide bond Cys20-Cys269 is essential in controlling specific orientation of Nt region and maintaining conformational integrity of extracellular domain. RMS comparison analysis between conformers revealed the ECL1 of CCR5 stays relative rigid, whereas the ECL2 and Nt are rather flexible. Solvent-accessible surface area calculations indicated that the charged residues within Nt and ECL2 are often exposed to solvent. Integrating these results with available experimental data, a two-step gp120-CCR5 binding mechanism was proposed. The dynamic interaction of CCR5 extracellular domain with gp120 was emphasized. (C) 2004 Elsevier B.V. All rights reserved.

Molecular Dynamics Simulation on Characteristics of Decomposition of Methane Hydrate in the Presence of Hydrogen Peroxide Solution

In this work, the characteristics of the decomposition of methane hydrate Structure I (SI) in the presence of hydrogen peroxide solution is investigated using the molecular dynamics simulation. The mechanism of the transformation process from the solid hydrate to the liquid is analyzed with the effect of hydrogen peroxide (HP) solution. In addition, the effect of ethylene glycol (EG) with the same molar concentration with HP on the methane hydrate dissociation is also studied. The results illustrate that both HP and EG promote well the hydrate dissociation. The work provides the important reference value for the experimental investigation into the promotion effect of HP on the hydrate dissociation.

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.

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.