原子尺度摩擦行为的分子动力学模拟

应用大规模分子动力学方法,模拟了具有原子级光滑和原子级粗糙形貌的刚性球形探头与弹性平面基体的干摩擦行为,研究了无/有粘附条件下的载荷与摩擦力、载荷与真实接触面积,以及摩擦力与真实接触面积之间的关系,对纳米尺度下的摩擦行为规律进行了分析.几种系统的真实接触面积-载荷关系都与相应的连续力学接触模型定性的一致,它们分别是Hertz光滑表面接触模型、Greenwood-Williamson粗糙表面接触模型和Mau-gis Dugdale粘着接触模型.无论是由光滑表面还是粗糙表面构成的摩擦系统,在无粘附条件下摩擦力与载荷成正比,而摩擦力与真实接触面积之间没有一个简单的关系;在粘附条件下摩擦力与真实接触面积成正比,而摩擦力与载荷之间表现为Maugis Dugdale模型预测的亚线性关系.研究表明,当表面作用从无粘附到粘附时,控制摩擦力的决定因素从载荷转变为接触面积,摩擦行为从载荷控制摩擦转变为粘着控制摩擦.     

Molecular dynamics simulation of the role of dislocations in microcrack healing

The molecular dynamics method is used to simulate microcrack healing during heating or/and under compressive stress. A centre microcrack in Cu crystal would be sealed under compressive stress or by heating. The role of compressive stress and heating in crack healing was additive. During microcrack healing, dislocation generation and motion occurred. When there were pre-existing dislocations around the microcrack, the critical temperature or compressive stress necessary for microcrack healing would decrease, and, the higher the number of dislocations, the lower the critical temperature or compressive stress. The critical temperature necessary for microcrack healing depended upon the orientation of the crack plane. For example, the critical temperature for the crack along the (001) plane was the lowest, i.e. 770K. The project supported by the Special Fund for the Major State Basic Research Projects (No. G19990650) and by the National Natural Science Foundation of China (No. 19891180, 59871010)     

Molecular dynamics studies on the dislocation gliding near a tilt boundary

The gliding behavior of edge dislocation near a grain boundary (GB) 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 σ _c =23MPa ∼ 68 MPa and 137MPa∼274 MPa for Ω=165 small angle tilt GB at 300K and 20K, respectively. The first result agrees with the experimental yield stress σ _y (=42MPa) 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 σ _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 Ω=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. The project supported by KM85-33 of Academia Sinica and the National Natural Science Foundation of China     

Molecular dynamics of dewetting of ultra-thin water films on solid substrate

Molecular dynamics simulation is applied to study the instability and rupture process of ultra-thin water films on a solid substrate.Results show the small disturbance of the film will develop linearly due to the spinodal instability,whereas the interaction between solid and liquid has less influences on the initial growth.Then the rupture occurs and the rim recedes with a dynamic contact angle.The radius of the rim varies with time as the square root of t,which is consistent with the macroscopic theory available. Stronger interaction between solid and liquid will postpone rupture time,decline the dynamic contact angle and raise the density of water near the interface between solid and liquid.     

固体塑性实验室时间尺度的分子动力学模拟概述

随着超级计算机软硬件的飞速提升,基于经验势函数的分子动力学模拟在解析固体塑性的微观机制方面发挥着关键作用.但是,由于传统分子动力学基于牛顿运动方程数值积分,积分时间步长通常为飞秒量级,其模拟的时间尺度通常限于纳秒量级,从而为模拟长时间尺度固体塑性机制带来了巨大的挑战.本文从分子动力学模拟的时间尺度限制切入,介绍目前国际...     

分子沉积膜纳米压痕过程的分子动力学模拟

在利用联合原子模型并考虑静电力作用的基础上，采用分子动力学模拟方法研究了金探针作用下分子沉积膜的纳米压痕过程，并对其纳米力学行为进行了理论分析．结果表明，分子沉积膜在探针压下过程中出现了明显的接触跳跃现象，探针下面的膜的分子倾角和法向载荷呈现出相同的滞后趋势，这可能是探针和分子沉积膜之间粘着力作用的结果．     

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°(Σ11), 129.5°(Σ11), 70.5°(Σ3) and 109.5°(Σ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°(Σ11), 0.32T _m for 129.5° (Σ11) and 0.38T _m for 70.5°(Σ3), respectively. These results agree with conclusions deduced from the anelastic measurements. But, for twin-boundary structure 109.5°(Σ3), this disordering has not been found even when temperature increases up to 0.9T _m . The project supported by the National Natural Science Foundation of China and Laboratory for Non-linear Mechanics of Continuous Media, Institute of Mechanics, Academia Sinica.     

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 project supported by the National Natural Science Foundation of China     

A study on the measurement of mean velocity and its convergence in molecular dynamics simulations

In many particle‐based simulations, measurement of local mean flow velocity and other continuum‐based properties are of utmost importance. Macroscopic quantities, such as mean flow velocity, temperature, and density, can be estimated by averaging the corresponding microscopic behavior of the particles. The two main subjects that should be considered in the averaging over the particles in a specific problem are spatial and temporal behaviors of them. In this paper, we study the latter. Because of the chaotic nature of the collisions among the molecules and consequently their random path, extracted macroscopic values fluctuate about their average values causing statistical errors. In this paper, an averaging method called SAM‐Modified‐CAM (SMC) will be proposed for the measurement of mean velocity that reduces statistical errors in its calculation. This proposal is based on the study conducted here on the implementations of two common averaging methods, sample‐averaged measurement (SAM) and cumulative average measurement (CAM) in molecular dynamics. In addition, convergence of mean flow velocity measurement is thoroughly discussed, and a convergence criterion is proposed for this purpose. Implementation of the proposed method in different test cases has approved its reliable performance. Copyright © 2010 John Wiley & Sons, Ltd.     

纳米晶铜单向拉伸变形的分子动力学模拟

纳米材料是由尺度在1－100nm的微小颗粒组成的体系，由于它具有独特的性能而备受关注。本文简要地回顾了分子动力学在纳米材料研究中的应用，并运用它模拟了平均晶粒尺寸从1．79－5．38nm的纳米晶体的力学性质。模拟结果显示：随着晶粒尺寸的减小，系统与晶粒内部的原子平均能量升高，而晶界上则有所下降；纳米晶体的弹性模量要小于普通多晶体，并随着晶粒尺寸的减小而减小；纳米晶铜的强度随着晶粒的减小而减小，显示了反常的Hall-Petch效应；纳米晶体的塑性变形主要是通过晶界滑移与运动，以及晶粒的转动来实现的；位错运动起着次要的、有限的作用；在较大的应变下（约大于5％），位错运动开始起作用；这种作用随着晶粒尺寸的增加而愈加明显。     

Molecular dynamics study on mechanics of metal nanowire

The new concept of using nanowires as building blocks for logic and memory circuits makes it very necessary to fully understand the mechanical behaviors of these nanowires. Embedded-atom method is employed to carry out three-dimensional molecular dynamics simulations of the mechanical properties of rectangular cross-section copper nanowire. A stable free-relaxation state and the stress–strain relation of nanowire under extension are obtained. The elastic modulus, yielding strength and deformation are studied. The surface effect, size effect, and temperature effect on the extension property of metal nanowire are discussed in detail. The simulation results from our present work show that at nanoscale surface atoms play an important role on the mechanical behaviors of nano-structures. This study of mechanical properties of metal nanowires will be helpful to the design, manufacture and manipulation of nano-devices.     

自组装膜纳米压痕的分子动力学模拟

研究了金探针对沉积在金 (111)表面的 CH3 (CH2 ) 1 5S自组装膜的纳米压痕的分子动力学模拟 .结果表明 ,自组装膜在金探针的作用下出现了跳跃接触现象 ,倾角和法向载荷都出现明显的滞后 ,并进一步表明与探针和膜之间的粘着力有关     

分子动力学模拟含不同初始缺陷单晶镍的断裂行为与应力演化

断裂是一个跨尺度复杂的物理过程,对宏观尺度的断裂行为已有深入的研究和发展,然而对微观尺度的断裂行为及断裂过程中应力场的变化缺乏深入的理解。本文通过分子动力学模拟,研究了具有不同初始缺陷(尖锐裂纹、钝裂纹和孔洞)的单晶镍的断裂行为和应力分布特征。结果表明,不同的初始缺陷导致了不同的断裂机制、断裂强度和抗断裂性能。含初始孔洞的单晶镍样品有最高的断裂强度和最强的抗断裂性能,这与孔洞扩展过程中堆积层错的形成密切相关。其次是含初始钝裂纹的样品,在裂纹扩展过程中出现由[100]超位错发射引起的裂尖钝化;含尖锐裂纹的样品表现为脆性断裂,裂尖原子没有出现微结构的变化,其强度和抗断裂性能最低。此外,不同的初始缺陷也会导致断裂过程中应力分布的变化,对含有尖锐裂纹的脆性断裂试样,高应力(拉伸应力、平均应力和米塞斯应力)总是出现在扩展裂纹的裂尖。而对于含有钝裂纹或孔洞的韧性断裂试样,高应力不仅分布在裂尖,也分布在位错发射和堆积层错形成的区域,在裂纹/孔洞扩展之前,应力随着加载时间的增加而迅速增加,而一旦裂纹或孔洞开始扩展,应力增加非常缓慢或几乎不增加,但拉伸应力值始终大于平均应力和米塞斯应力值。这表明,在I型...     

Bin size determination for the measurement of mean flow velocity in molecular dynamics simulations

Extracting macroscopic properties from molecular simulation of fluids is required in most of the molecular dynamics problems. However, methods used for this purpose, their accuracy, and their dependence on the sampling and averaging schemes are still ambiguous. Macroscopic properties at any point of a molecular domain can be extracted via sampling and averaging of molecular behavior within a control region around that point, called a bin. The size of this bin has a significant impact on the accuracy of the macroscopic properties' measurement. In this research, we will focus on the measurement of mean flow velocity using ‘binning method’. On the basis of the results of this study, the most appropriate range of the bin size in which mean velocity can be measured accurately is determined. Although this range is determined based on the measurement of mean flow velocity, we believe that it can be used to estimate the proper range for the measurement of other flow quantities as well. Copyright © 2012 John Wiley & Sons, Ltd.     

Explicit solvent molecular dynamics simulations of chaperonin-assisted rhodanese folding

Chaperonins are known to facilitate the productive folding of numerous misfolded proteins, Despite their established importance, the mechanism of chaperonin-assisted protein folding remains unknown. In the present article, all-atom explicit solvent molecular dynamics (MD) simulations have been performed for the first time on rhodanese folding in a series of cavity-size and cavity-charge chaperonin mutants. A compromise between stability and flexibility of chaperonin structure during the substrate folding has been observed and the key factors affecting this dynamic process are discussed.     

Size and temperature effects on the fracture mechanisms of silicon nanowires: Molecular dynamics simulations

We present molecular dynamics simulations of [1 1 0]-oriented Si nanowires (NWs) under a constant strain rate in tension until failure, using the modified embedded-atom-method (MEAM) potential. The fracture behavior of the NWs depends on both temperature and NW diameter. For NWs of diameter larger than 4 nm, cleavage fracture on the transverse (1 1 0) plane are predominantly observed at temperatures below 1000 K. At higher temperatures, the same NWs shear extensively on inclined {1 1 1} planes prior to fracture, analogous to the brittle-to-ductile transition (BDT) in bulk Si. Surprisingly, NWs with diameter less than 4 nm fail by shear regardless of temperature. Detailed analysis reveals that cleavage fracture is initiated by the nucleation of a crack, while shear failure is initiated by the nucleation of a dislocation, both from the surface. While dislocation mobility is believed to be the controlling factor of BDT in bulk Si, our analysis showed that the change of failure mechanism in Si NWs with decreasing diameters is nucleation controlled. Our results are compared with a recent in situ tensile experiment of Si NWs showing ductile failure at room temperature.     

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. The project supported by the National Natural Science Foundation of China and the Science Foundation of Chinese Academy of Sciences.     

Molecular dynamics simulation of peeling a DNA molecule on substrate

Molecular dynamics (MD) simulations are performed to study adhesion and peeling of a short fragment of single strand DNA (ssDNA) molecule from a graphite surface. The critical peel-off force is found to depend on both the peeling angle and the elasticity of ssDNA. For the short ssDNA strand under investigation, we show that the simulation results can be explained by a continuum model of an adhesive elastic band on substrate. The analysis suggests that it is often the peak value, rather than the mean value, of adhesion energy which determines the peeling of a nanoscale material.The project supported by the Distinguished Young Scholar Fund of NSFC (10225209) and key project from the Chinese Academy of Sciences (KJCX2-SW-L2)