Transitions Between Different Contact Models

The transitions between the different contact models which include the Hertz, Bradley, Johnson-Kendall-Roberts (JKR), Derjaguin-Muller-Toporov (DMT) and Maugis-Dugdale (MD) models are revealed by analyzing their contact pressure profiles and surface interactions. Inside the contact area, surface interaction/adhesion induces tensile contact pressure around the contact edge. Outside the contact area, whether or not to consider the surface interaction has a significant influence on the contact system equilibrium. The difference in contact pressure due to the surface interaction inside the contact area and the equilibrium influenced by the surface interaction outside the contact area are physically responsible for the different results of the different models. A systematic study on the transitions between different models is shown by analyzing the contact pressure profiles and the surface interactions both inside and outside the contact area. The definitions of contact radius and the flatness of contact surfaces are also discussed. (C) Koninklijke Brill NV, Leiden, 2008.

Molecular/Cluster Statistical Thermodynamics Methods To Simulate Quasi-Static Deformations At Finite Temperature

The rapid evolution of nanotechnology appeals for the understanding of global response of nanoscale systems based on atomic interactions, hence necessitates novel, sophisticated, and physically based approaches to bridge the gaps between various length and time scales. In this paper, we propose a group of statistical thermodynamics methods for the simulations of nanoscale systems under quasi-static loading at finite temperature, that is, molecular statistical thermodynamics (MST) method, cluster statistical thermodynamics (CST) method, and the hybrid molecular/cluster statistical thermodynamics (HMCST) method. These methods, by treating atoms as oscillators and particles simultaneously, as well as clusters, comprise different spatial and temporal scales in a unified framework. One appealing feature of these methods is their "seamlessness" or consistency in the same underlying atomistic model in all regions consisting of atoms and clusters, and hence can avoid the ghost force in the simulation. On the other hand, compared with conventional MD simulations, their high computational efficiency appears very attractive, as manifested by the simulations of uniaxial compression and nanoindenation. (C) 2008 Elsevier Ltd. All rights reserved.

Effect Of Load Triaxiality On Polymorphic Transitions In Zinc Oxide

Molecular dynamics (MD) simulations and first-principles calculations are carried out to analyze the stability of both newly discovered and previously known phases of ZnO under loading of various triaxialities. The analysis focuses on a graphite-like phase (FIX) and a body-centered-tetragonal phase (BCT-4) that were observed recently in [0 1 (1) over bar 0]- and [0 0 0 1]-oriented nanowires respectively under uniaxial tensile loading as well as the natural state of wurtzite (WZ) and the rocksalt (RS) phase which exists under hydrostatic pressure loading. Equilibrium critical stresses for the transformations are obtained. The WZ -> HX transformation is found to be energetically favorable above a critical tensile stress of 10 GPa in [0 1 (1) over tilde 0] nanowires. The BCT-4 phase can be stabilized at tensile stresses above 7 GPa in [0 0 0 1] nanowires. The RS phase is stable at hydrostatic pressures above 8.2 GPa. The identification and characterization of these phase transformations reveal a more extensive polymorphism of ZnO than previously known. A crystalline structure-load triaxiality map is developed to summarize the new understanding. (c) 2007 Elsevier Ltd. All rights reserved.

Surface-induced melting of metal nanoclusters

We investigate the size effect on melting of metal nanoclusters by molecular dynamics simulation and thermo dynamic theory based on Kofman's melt model. By the minimization of the free energy of metal nanoclusters with respect to the thickness of the surface liquid layer, it has been found that the nanoclusters of the same metal have the same premelting temperature T-pre = T-0 - T-0(gamma(su) - gamma(lv) - gamma(sl))/(rhoLxi) (T-0 is the melting point of bulk metal, gamma(sv) the solid-vapour interfacial free energy, gamma(sl) the liquid-vapour interfacial free energy, gamma(sl),l the solid-liquid interfacial free energy, p the density of metal, L the latent heat of bulk metal, and xi the characteristic length of surface-interface interaction) to be independent of the size of nanoclusters, so that the characteristic length of a metal can be obtained easily by T-pre, which can be obtained by experiments or molecular dynamics (MD) simulations. The premelting temperature T-pre of Cu is obtained by AID simulations, then xi is obtained. The melting point T-cm is further predicted by free energy analysis and is in good agreement with the result of our MD simulations. We also predict the maximum premelting-liquid width of Cu nanoclusters with various sizes and the critical size, below which there is no premelting.

A dynamic coupling model for hybrid atomistic-continuum computations

A dynamic coupling model is developed for a hybrid atomistic-continuum computation in micro- and nano-fluidics. In the hybrid atomistic-continuum computation, a molecular dynamics (MD) simulation is utilized in one region where the continuum assumption breaks down and the Navier-Stokes (NS) equations are used in another region where the continuum assumption holds. In the overlapping part of these two regions, a constrained particle dynamics is needed to couple the MD simulation and the NS equations. The currently existing coupling models for the constrained particle dynamics have a coupling parameter, which has to be empirically determined. In the present work, a novel dynamic coupling model is introduced where the coupling parameter can be calculated as the computation progresses rather than inputing a priori. The dynamic coupling model is based on the momentum constraint and exhibits a correct relaxation rate. The results from the hybrid simulation on the Couette flow and the Stokes flow are in good agreement with the data from the full MD simulation and the solutions of the NS equations, respectively. (c) 2007 Elsevier Ltd. All rights reserved.

Atomistic investigation of the effects of temperature and surface roughness on diffusion bonding between Cu and Al

Molecular dynamics (MD) simulations are carried out to analyze the diffusion bonding at Cu/Al interfaces. The results indicate that the thickness of the interfacial layer is temperature-dependent, with higher temperatures yielding larger thicknesses. At temperatures below 750 K, the interface thickness is found to increase in a stepwise manner as a function of time. At temperatures above 750 K, the thickness increases rapidly and smoothly. When surface roughness is present, the bonding process consists of three stages. In the first stage, surfaces deform under stress, resulting in increased contact areas. The second stage involves significant plastic deformation at the interface as temperature increases, resulting in the disappearance of interstices and full contact of the surface pair. The last stage entails the diffusion of atoms under constant temperature. The bonded specimens show tensile strengths reaching 88% of the ideal Cu/Al contact strength. (c) 2007 Acta Materialia Inc. Published by Elsevier Ltd. All rights reserved.

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.

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).

纳米硬度测量中接触面积及压头曲率半径效应的分子动力学模拟

用分子动力学(MD)模拟了单晶铜的纳米压痕。提出“接触原子”方法来计算接触面积,并和实验上采用的方法进行比较,发现此方法无论是对“挤出”还是“沉陷”情况都能更准确地计算接触面积。同时也模拟了压头曲率半径对纳米压痕的影响,发现在同一压痕深度下材料纳米硬度的测量值随着压头半径的增大而变大。

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 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.

Investigations of cohesive zone properties of the interface between metal film and ceramic substrate at nano- and micro-scales

Cohesive zone characterizations of the interface between metal film and ceramic substrate at micro- and nano-scales are performed in the present research. At the nano-scale, a special potential for special material interface (Ag/MgO) is adopted to investigate the interface separation mechanism by using MD simulation, and stress-separation relationship will be obtained. At the micro-scale, peeling experiment is performed for the Al film/Al2O3 substrate system with an adhesive layer at the interface. Adhesive is a mixture of epoxy and polyimide with mass ratio 1:1, by which a brittle cohesive property is obtained. The relationships between energy release rate, the film thickness and the adhesive layer thickness are measured during the steady-state peeling. The experimental result has a similar trend as modeling result for a weak adhesion interface case.

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.