Molecular dynamics studies on the dislocation gliding near a tilt boundary

The gliding behavior of edge dislocation near a grain boundary(QB) in copper under pure shear stresses is simulated by using molecular dynamics(MD) method. Many-body potential incorporating the embedded atom method (EAM) is used. The critical shear stresses for a single disocation to pass across GB surface are obtained at values of sigma(c)=23MPa similar to 68 MPa and 137 MPa similar to 274 MPa for Sigma=165 small angle tilt GB at 300 K and 20 K, respectively. The first result agrees with the experimental yield stress sigma(y)(=42 MPa) quite well. It suggests that there might be one of the reasons of initial plastic yielding caused by single dislocation gliding across GB. In addition, there might be possibility to obtain yield strength from microscopic analysis. Moreover, the experimental value of sigma(y) at low temperature is generally higher than that at room temperature. So, these results are in conformity qualitatively with experimental fact. On the other hand, the Sigma=25 GB is too strong an obstacle to the dislocation. In this case, a dislocation is able to pass across GB under relatively low stress only when it is driven by other dislocations. This is taken to mean that dislocation pile-up must be built up in front of this kind of GB, if this GB may take effect on the process of plastic deformation.

Molecular dynamics studies on the local disordering of σ3 and σ11 grain boundaries in aluminium bicrystal

Molecular dynamics (MD) simulations using Morse interaction potential are performed in studies of [110] symmetrical tilt grain boundary (GB) structures with mis-orientation angles 50.5 degrees(Sigma 11), 129.5 degrees(Sigma 11), 70.5 degrees(Sigma 3) and 109.5 degrees(Sigma 3) at various tempratures. The GB structures are found to start local disordering at about 0.5T(m)(T-m is the melting point of aluminium) for 50.5 degrees(Sigma 11), 0.32T(m) for 129.5 degrees(Sigma 11) and 0.38T(m) for 70.5 degrees(Sigma 3), respectively. These results agree with conclusions deduced from the anelastic measurements. But, for twin-boundary structure 109.5 degrees(Sigma 3), this disordering has not been found even when temperature increases up to 0.9T(m).

Molecular dynamics simulation of crack-tip processes in copper

The crack tip processes in copper under mode II loading have been simulated by a molecular dynamics method. The nucleation, emission, dislocation free zone (DFZ) and pile-up of the dislocations are analyzed by using a suitable atom lattice configuration and Finnis & Sinclair potential. The simulated results show that the dislocation emitted always exhibits a dissociated fashion. The stress intensity factor for dislocation nucleation, DFZ and dissociated width of partial dislocations are strongly dependent on the loading rate. The stress distributions are in agreement with the elasticity solution before the dislocation emission, but are not in agreement after the emission. The dislocation can move at subsonic wave speed (less than the shear wave speed) or at transonic speed (greater than the shear wave speed but less than the longitudinal wave speed), but at the longitudinal wave speed the atom lattice breaks down.

QM/MM and Classical Molecular Dynamics Simulation of Histidine-Tagged Peptide Immobilization on Nickel Surface

The hybrid quantum mechanics (QM) and molecular mechanics (MM) method is employed to simulate the His-tagged peptide adsorption to ionized region of nickel surface. Based on the previous experiments, the peptide interaction with one Ni ion is considered. In the QM/MM calculation, the imidazoles on the side chain of the peptide and the metal ion with several neighboring water molecules are treated as QM part calculated by "GAMESS", and the rest atoms are treated as MM part calculated by "TINKER". The integrated molecular orbital/molecular mechanics (IMOMM) method is used to deal with the QM part with the transitional metal. By using the QM/MM method, we optimize the structure of the synthetic peptide chelating with a Ni ion. Different chelate structures are considered. The geometry parameters of the QM subsystem we obtained by QM/MM calculation are consistent with the available experimental results. We also perform a classical molecular dynamics (MD) simulation with the experimental parameters for the synthetic peptide adsorption on a neutral Ni(100) surface. We find that half of the His-tags are almost parallel with the substrate, which enhance the binding strength. Peeling of the peptide from the Ni substrate is simulated in the aqueous solvent and in vacuum, respectively. The critical peeling forces in the two environments are obtained. The results show that the in-tidazole rings are attached to the substrate more tightly than other bases in this peptide.

Molecular Dynamics Simulation of Barnacle Cement

Barnacle cement is an underwater adhesive that is used for permanent settlement. Its main components are insoluble protein complexes that have not been fully studied. In present article, we chose two proteins of barnacle cement for study, 36-KD protein and Mrcp-100K protein. In order to investigate the characteristic of above two proteins, we introduced the method of molecular modeling. And the simulation package GROMACS was used to simulate the behavior of these proteins. In this article, before the simulations, we introduce some theories to predict the time scale for polymer relaxation. During the simulation, we mainly focus on two properties of these two proteins: structural stability and adhesive force to substrate. First, we simulate the structural stability of two proteins in water, and then the stability of 36-KD protein in seawater environment is investigated. We find that the stability varies in the different environments. Next, to study adhesive ability of two proteins, we simulate the process of peeling the two proteins from the substrate (graphite). Then, we analyze the main reasons of these results. We find that hydrogen bonds in proteins play an important role in the protein stability. In the process of the peeling, we use Lennard-Jones 12-6 potential to calculate the van der Waals interactions between proteins and substrate.

Molecular Dynamics Simulation of the Bio-Adhesion in Molecular Motors

Molecular dynamics (MD) simulation is employed to study the bio-adhesion in F1 ATP molecular motor. Histidine-peptide is widely used as linkage in micro systems because of its strong binding strength to metals. This paper focuses on the adhesion between a synthetic peptide containing 6xHis-tag (Gly-Gly-Lys-Gly-Gly-Lys-Gly-Gly-His-His-His-His-His-His) and metal substrate, which is used to define the position of the F1 ATP molecular motor on the metal substrate. It is shown that the binding strength between histidine and nickel substrate is the strongest, while that of copper is smaller and that of gold substrate is the smallest. From the result of simulation, we find that the stability of adhesion between histidine and the metal substate result of the ringed structure in histidine.

Molecular dynamics studies of brittle fracture of SiO2

By applying for molecular dynamics (MD) simulation and Griffith fracture criterion, the brittle behavior of crack extension of mode I type is investigated. The critical stress intensity factor (SIF)K-Ic(MD) of crack extension is calculated, and the evolution of atoms near crack tip is observed. It is found that K-Ic(MD) is in good agreement with the Griffith ftacture criterion K-Ic(Griffith).

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.

Effect Of Water On Brittle Fracture Of Sio2 By Molecular Dynamics Study

The influence of water on the brittle behavior of beta-cristobalite is studied by means of molecular dynamics (MD) simulation With the TTAM potential. Crack extension of mode 1 type is observed as the crack opening is filled LIP With water. The critical stress intensity factor K-lc(MD) is used to characterize the crack extension of MD simulation. The surface energy of SiO2 covered with layers of water is calculated at temperature of 300 K. Based oil the Griffith fracture criterion, the critical stress intensity factor K-lc(Griffith) is calculated, and it is in good agreement with that of MD simulation. (C) 2008 Elsevier B.V. All rights reserved.

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.

Splitting the fast and slow motions in molecular dynamics simulations based on the change of cold potential well bottom

Abstract. The atomic motion is coupled by the fast and slow components due to the high frequency vibration of atoms and the low frequency deformation of atomic lattice, respectively. A two-step approximate method was presented to determine the atomic slow motion. The first step is based on the change of the location of the cold potential well bottom and the second step is based on the average of the appropriate slow velocities of the surrounding atoms. The simple tensions of one-dimensional atoms and two-dimensional atoms were performed with the full molecular dynamics simulations. The conjugate gradient method was employed to determine the corresponding location of cold potential well bottom. Results show that our two-step approximate method is appropriate to determine the atomic slow motion under the low strain rate loading. This splitting method may be helpful to develop more efficient molecular modeling methods and simulations pertinent to realistic loading conditions of materials.