A Molecular Dynamics Simulation: the Effect of Finite Size on the Thermal Conductivity in a Single Crystal Silicon

Non-equilibrium molecular dynamics (NEMD) simulations are performed to calculate thermal conductivity. The environment-dependent interatomic potential (EDIP) potential on crystal silicon is adopted as a model system. The issues are related to nonlinear response, local thermal equilibrium and statistical averaging. The simulation results by non-equilibrium molecular dynamics show that the calculated thermal conductivity decreases almost linearly as the film thickness reduced at the nanometre scale. The effect of size on the thermal conductivity is also obtained by a theoretic analysis of the kinetic theory and formulas of the heat capacity. The analysis reveals that the contributions of phonon mean free path (MFP) and phonon number in a finite cell to thermal conductivity are very important.

Molecular Dynamics And Density Functional Studies Of A Body-Centered-Tetragonal Polymorph Of Zno

We report a previously unknown body-centered-tetragonal structure for ZnO. This structure results from a phase transformation from wurtzite in [0001]-oriented nanorods during uniaxial tensile loading and is the most stable phase for ZnO when stress is above 7 GPa. The stress-induced phase transformation has important implications for the electronic, piezoelectric, mechanical, and thermal responses of ZnO. The discovery of this polymorph brings about a more complete understanding of the extent and nature of polymorphism in ZnO. A crystalline structure-load triaxiality map is developed to summarize the relationship between structure and loading.

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.

Molecular Dynamics Study Of Mechanical Behaviour Of Screw Dislocation During Cutting With Diamond Tip On Silicon

By means of Tersoff and Morse potentials, a three-dimensional molecular dynamics simulation is performed to study atomic force microscopy cutting on silicon monocrystal surface. The interatomic forces between the workpiece and the pin tool and the atoms of workpiece themselves are calculated. A screw dislocation is introduced into workpiece Si. It is found that motion of dislocations does not occur during the atomic force microscopy cutting processing. Simulation results show that the shear stress acting on dislocation is far below the yield strength of Si.

Grain Boundary Effects On Plastic Deformation And Fracture Mechanisms In Cu Nanowires: Molecular Dynamics Simulations

Metallic nanowires have many attractive properties such as ultra-high yield strength and large tensile elongation. However, recent experiments show that metallic nanowires often contain grain boundaries, which are expected to significantly affect mechanical properties. By using molecular dynamics simulations, here, we demonstrate that polycrystalline Cu nanowires exhibit tensile deformation behavior distinctly different from their single-crystal counterparts. A significantly lowered yield strength was observed as a result of dislocation emission from grain boundaries rather than from free surfaces, despite of the very high surface to volume ratio. Necking starts from the grain boundary followed by fracture, resulting in reduced tensile ductility. The high stresses found in the grain boundary region clearly play a dominant role in controlling both inelastic deformation and fracture processes in nanoscale objects. These findings have implications for designing stronger and more ductile structures and devices on nanoscale.

Molecular Dynamics Study Of The Mechanical Behavior Of Nickel Nanowire: Strain Rate Effects

We present the analysis of uniaxial deformation of nickel nanowires using molecular dynamics simulations, and address the strain rate effects on mechanical responses and deformation behavior. The applied strain rate is ranging from 1 x 10(8) s(-1) to 1.4 x 10(11) s(-1). The results show that two critical strain rates, i.e., 5 x 10(9) s(-1) and 8 x 10(10) s(-1), are observed to play a pivotal role in switching between plastic deformation modes. At strain rate below 5 x 10(9) s(-1), Ni nanowire maintains its crystalline structure with neck occurring at the end of loading, and the plastic deformation is characterized by {111} slippages associated with Shockley partial dislocations and rearrangements of atoms close to necking region. At strain rate above 8x10(10) s(-1), Ni nanowire transforms from a fcc crystal into a completely amorphous state once beyond the yield point, and hereafter it deforms uniformly without obvious necking until the end of simulation. For strain rate between 5 x 10(9) s(-1) and 8 x 10(10) s(-1), only part of the nanowire exhibits amorphous state after yielding while the other part remains crystalline state. Both the {111} slippages in ordered region and homogenous deformation in amorphous region contribute to the plastic deformation. (C) 2007 Published by Elsevier B.V.

Study Of Mechanical Properties Of Amorphous Copper With Molecular Dynamics Simulation

The formation and mechanical properties of amorphous copper are studied using molecular dynamics simulation. The simulations of tension and shearing show that more pronounced plasticity is found under shearing, compared to tension. Apparent strain hardening and strain rate effect are observed. Interestingly, the variations of number density of atoms during deformation indicate free volume creation, especially under higher strain rate. In particular, it is found that shear induced dilatation does appear in the amorphous metal.

Molecular dynamics simulations of DFZ

Dislocation emission from the crack tip in copper under mode II loading is simulated with molecular dynamics method. After 26 partial dislocations are emitted and then relaxed to reach the equilibrium under the constant displacement, the double pile-ups (including an inverse pile-up and a pile-up) are formed. i.e., the first dislocation is piled up before the obstruction, and the last dislocation is piled up ahead of the crack tip. These results conform to the TEM observations.

Molecular dynamics studies on dislocations in crystallites of nanocrystalline alpha-iron

Molecular dynamics simulations have been carried our to study the atomic structure of the crystalline component of nanocrystalline alpha-iron. A two-dimensional computational block is used to simulate the consolidation process. It is found that dislocations are generated in the crystallites during consolidation when the grain size is large enough. The critical value of the grain size for dislocation generation appears to be about 9 nm. This result agrees with experiment qualitatively. AN dislocations that are preset in the original grains glide out during consolidation. It shows that dislocations in the crystallites we generated in consolidation process, but not in the original grains. Higher consolidation pressure results in more dislocations. Furthermore, new interfaces are found within crystallites. These interfaces might result from the special environment of nanomaterial. (C) 1998 Acta Metallurgica Inc.

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.