Size effects on the melting of nickel nanowires: a molecular dynamics study

The melting process of nickel nanowires are simulated by using molecular dynamics with the quantum Sutten-Chen many-body force field. The wires studied were approximately cylindrical in cross-section and periodic boundary conditions were applied along their length; the atoms were arranged initially in a face-centred cubic structure with the [0 0 1] direction parallel to the long axis of the wire. The size effects of the nanowires on the melting temperatures are investigated. We find that for the nanoscale regime, the melting temperatures of Ni nanowires are much lower than that of the bulk and are linear with the reciprocal of the diameter of the nanowire. When a nanowire is heated up above the melting temperature, the neck of the nanowire begins to arise and the diameter of neck decreases rapidly with the equilibrated running time. Finally, the breaking of nanowire arises, which leads to the formation of the spherical clusters. (C) 2004 Elsevier B.V. All rights reserved.

Molecular dynamics simulation of plastic deformation of nanotwinned copper

The plastic deformation of polycrystalline Cu with ultrathin lamella twins has been studied using molecular dynamics simulations. The results of uniaxial tensile deformation simulation show that the abundance of twin boundaries provides obstacles to dislocation motion, which in consequence leads to a high strain hardening rate in the nanotwinned Cu. We also show that the twin lamellar spacing plays a vital role in controlling the strengthening effects, i.e., the thinner the thickness of the twin lamella, the harder the material. Additionally, twin boundaries can act as dislocation nucleation sites as they gradually lose coherency at large strain. These results indicate that controlled introduction of nanosized twins into metals can be an effective way of improving strength without suppression tensile ductility. (C) 2007 American Institute of Physics.

Adsorption of His-tagged peptide to Ni, Cu and Au (100) surfaces: Molecular dynamics simulation

Molecular dynamics (MD) simulations are performed to study the interaction of His-tagged peptide with three different metal surfaces in explicit water. The equilibrium properties are analyzed by using pair correlation functions (PCF) to give an insight into the behavior of the peptide adsorption to metal surfaces in water solvent. The intermolecular interactions between peptide residues and the metal surfaces are evaluated. By pulling the peptide away from the peptide in the presence of solvent water, peeling forces are obtained and reveal the binding strength of peptide adsorption on nickel, copper and gold. From the analysis of the dynamics properties of the peptide interaction with the metal surfaces, it is shown that the affinity of peptide to Ni surface is the strongest, while on Cu and An the affinity is a little weaker. In MD simulations including metals, the His-tagged region interacts with the substrate to an extent greater than the other regions. The work presented here reveals various interactions between His-tagged peptide and Ni/Cu/Au surfaces. The interesting affinities and dynamical properties of the peptide are also derived. The results give predictions for the structure of His-tagged peptide adsorbing on three different metal surfaces and show the different affinities between them, which assist the understanding of how peptides behave on metal surfaces and of how designers select amino sequences in molecule devices design. (c) 2007 Elsevier Ltd. All rights reserved.

Dislocations in the crystallites of nanocrystalline alpha-Fe - A molecular dynamics study

Molecular dynamics simulations are carried out in order to study the atomic structure of crystalline component, of nanocrystalline alpha-Fe when it is consolidated from small grains. A two-dimensional computational block is used to simulate the consolidation process. All the preset dislocations in the original grains glide out of them in the consolidation process, but new dislocations can generate when the grain size is large enough. It shows that dislocations exist in the consolidated material rather than in the original grains. Whether dislocations exist in the crystalline component of the resultant model nana-material depends upon grain size. The critical value of grain size for dislocation generation appears to be about 9 nm. This result agrees with experiments qualitatively.

Molecular dynamics simulation of interaction of a dislocation array from a crack tip with grain boundaries

The interaction of a dislocation array emitted from a crack tip under mode II loading with asymmetric tilt grain boundaries (GBs) is analysed by the molecular dynamics method. The GBs can generally be described by planar and linear matching zones and unmatching zones. All GBs are observed to emit dislocations. The GBs migrated easily due to their planar and linear matching structure and asymmetrical type. The diffusion induced by stress concentration is found to promote the GB migration. The transmissions of dislocations are either along the matched plane or along another plane depending on tilt angle theta. Alternate processes of stress concentration and stress relaxation take place ahead of the pileup. The stress concentration can be released either by transmission of dislocations, by atom diffusion along GBs, or by migration of GBs by formation of twinning bands. The simulated results also unequivocally demonstrate two processes, i.e. asymmetrical GBs evolving into symmetrical ones and unmatching zones evolving into matching ones during the loading process.

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.

The structure and the dynamics of poly(vinyl methyl ether) PVME and in PVME concentrated water solution : a study by neutron scattering and fully atomistic molecular dynamics simulations

175 p. : il.

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.