Molecular dynamics simulation of plastic effects upon spalling

The molecular dynamics method is applied to simulate spalling during the plane shock interaction between plates. The effect of lattice defects in a material on the propagation of a shock wave and the process of spalling is studied. The plastic effects are described using a model of imperfect particle packing with defects (vacancies). The model proposed can describe the separation of the shock-wave front into an elastic precursor and a plastic front and give velocity profiles for the free target surface close to the experimental profiles.     

金属Al表面熔化各向异性的分子动力学模拟

采用嵌入原子势,使用分子动力学方法对金属Al不同低指数晶面的表面熔化现象分别进行了模拟.分析了熔化过程中样品结构组态的变化.模拟结果表明对于不同的自由表面,表面熔化呈现出明显的各向异性行为.Al（110）面在低于熔点的温度之下发生预熔化;（111）与（001）面都出现过热现象.与（111）面不同,（001）面发生过热现象时表面原子层为类液层,而（111）面仍然保持很好的晶格结构,即预熔化的Al（001）面在高于熔点的温度下,仍可以在很长的时间内处于相对稳定的亚稳态.由模拟得到Al的热力学熔点为950 K左右,与实验值基本符合. 关键词： 分子动力学 表面熔化 过热     

Molecular dynamics simulation of deposition of nickel nanocluster on copper surface

The burrowing of nickel nanocluster deposited onto the copper surface is investigated by molecular dynamics (MD) simulation. The simulation is carried out at different temperatures for 40.1 ns with three different lattice orientations Cu(100), Cu(110), and Cu(111) that have 20 × 20 × 15, 14 × 14 × 15, and 12 × 14 × 9 lattice units, respectively. The Ni(100) nanocluster consists of 249 atoms (five lattices diameter) and the initial kinetic energy is assigned to be 0 eV. The results show that the burrowing process goes extremely slow as temperature is at or under 900 K. There is virtually no burrowing observed when the system temperature is below 500 K. The burrowing processes at different temperatures are discussed in terms of kinetic energy of the cluster exerted by the strong capillary force. It is found that the kinetic energy will play a key role in the acceleration of the burrowing process.     

金属钨中氦行为的分子动力学模拟

采用分子动力学方法模拟了氦在金属钨中的扩散聚集行为. 首先,建立了氦与钨原子间相互作用势,短程部分采用ZBL势形式,长程部分采用从头算法数据,实现了两者之间的平滑连接. 通过计算氦在钨中不同间隙位的形成能发现,单个氦原子更易存在于金属钨中的四面体间隙位,这与最新的研究成果是一致的. 在400-1200 K的温度范围内,考察了氦原子在金属钨中的扩散行为,获得了扩散迁移能,其值介于实验值和从头算法结果之间. 最后,研究了氦的聚集行为,从能量的角度考察了氦团簇形成初期的生长机理. 研究发现,在氦团簇形成初期,氦团簇对氦的结合能随着氦团簇的生长有逐渐增大的趋势,说明氦团簇吸收氦的能力逐渐增强. 关键词： 氦扩散 氦团簇 辐照损伤 分子动力学模拟     

Ni熔体凝固过程中临界晶核和亚临界晶核的分子动力学模拟

本文用分子动力学模拟研究了Ni熔体以不同冷速凝固后微观结构的演变规律, 并通过理论计算确定出了Ni熔体凝固后获得理想非晶的临界条件. 模拟结果发现冷速小于10¹¹ K/s时, Ni 熔体凝固后形成晶态组织; 冷速在10¹¹ K/s到10^14.5 K/s之间时, Ni熔体凝固后形成由晶态结构与非晶态结构组成的混合组织. 冷速小于10¹⁰ K/s, Ni 熔体凝固后形成的晶态组织具有fcc结构; 冷速在10¹⁰ K/s到10^14.5 K/s之间时, Ni熔体凝固后组织中的晶态由fcc和hcp结构层状镶嵌排列构成. 通过分析模拟结果和计算结果, 确定出了Ni熔体凝固后形成理想非晶的临界冷速为10^14.5 K/s. 并发现Ni熔体中临界晶核(冷速等于10^14.5 K/s)和亚临界晶核(冷速大于10^14.5 K/s) 均由fcc和hcp组成层状偏聚结构, 这表明Ni熔体中生长的晶体、临界晶核和晶胚的结构是相同的. 关键词： 分子动力学模拟 晶体团簇 临界冷速 结构     

Pb液滴在Ni基底润湿铺展行为的分子动力学模拟

采用分子动力学方法研究了Pb液滴在Ni(100)晶面、(110)晶面和(111)晶面的铺展润湿行为. 结果表明: Pb液滴在Ni(100)及(111)基底上的润湿铺展现象呈各向同性, 而在Ni(110)基底上的润湿铺展现象呈明显的各向异性, 且这种各向异性源于Ni(110)晶面点阵结构上Pb原子沿不同晶向的扩散机制及速度的明显差异; Pb液滴在Ni(111)晶面上铺展时, 未发生表面合金化, 液滴铺展动力学描述近似满足 R² ∝ t, 而液滴在(100)晶面和(110)晶面上铺展时表面产生合金化现象, 铺展动力学关系近似满足 R⁴ ∝ t, 且液滴在(100)晶面上的铺展速度高于(110)晶面上的铺展速度. 关键词： 分子动力学 润湿各向异性 铺展膜 扩散机制     

单晶Cu冲击加载及卸载下塑性行为的微观模拟

使用分子动力学方法对室温下单晶铜沿［001］和［111］方向冲击加载及卸载下的塑性行为进行了模拟,得到了Hugoniot关系以及冲击熔化压力,与实验基本符合. 加载过程中,较高的初始温度有利于位错的形核与发展. 通过对冲击波在自由表面卸载过程的模拟和分析发现：卸载过程呈现“准弹性卸载行为”;沿［001］方向卸载后大量不全位错环与堆积层错消失,而沿［111］方向卸载后只有少量层错消失,部分层错甚至会发展扩大. 关键词： 分子动力学 冲击波 塑性     

Molecular dynamics simulation of thermal accommodation coefficients for laser-induced incandescence sizing of nickel particles

Extending time-resolved laser-induced incandescence (TiRe-LII), a diagnostic traditionally used to characterize soot and other carbonaceous particles, into a tool for measuring metal nanoparticles requires knowledge of the thermal accommodation coefficient for those systems. This parameter can be calculated using molecular dynamics (MD) simulations provided the interatomic potential is known between the gas molecule and surface atoms, but this is not often the case for many gas/surface combinations. In this instance, researchers often resort to the Lorentz–Berthelot combination rules to estimate the gas/surface potential using parameters derived for homogeneous systems. This paper compares this methodology with a more accurate approach based on ab initio derived potentials to estimate the thermal accommodation coefficient for laser-energized nickel nanoparticles in argon. Results show that the Lorentz–Berthelot combining rules overestimate the true potential well depth by an order of magnitude, resulting in perfect thermal accommodation, whereas the more accurate ab initio derived potential predicts an accommodation coefficient in excellent agreement with experimentally-determined values for other metal nanoparticle aerosols. This result highlights the importance of accurately characterizing the gas/surface potential when using MD to estimate thermal accommodation coefficients for TiRe-LII.     

铜-铝扩散焊及拉伸的分子动力学模拟

采用分子动力学方法模拟了铜-铝扩散焊过程，分析了理想平面铜-铝试件(001)晶面间扩散焊的过渡层厚度，并利用径向分布、键对分析方法分析了在不同的降温速率下过渡层的结构变化.降温速率大时，过渡层保持原有无序结构，降温速率小时，过渡层从无序结构向面心立方结构转变.还对扩散焊后的铜-铝试件进行了拉伸模拟，并与尺寸大小相近的单晶铜和单晶铝的拉伸模拟结果进行比较.结果发现焊接后的强度比单晶铝和单晶铜的强度都要小，最大应变值也小.     

纳米团簇熔化过程的分子动力学模拟

本文采用分子动力学结合嵌入原子多体势,模拟了不同半径的Ni纳米团簇的升温熔化过程,研究团簇尺寸对熔点和表面能的影响.模拟结果表明:团簇的熔点显著低于体材料的熔点.团簇熔化的过程首先是在团簇的表面出现预熔,然后向团簇内部扩展,直到整个团簇完全熔为液态.在模拟的纳米尺度范围内,团簇的熔点与团簇尺寸基本成线性关系.团簇的表面能随着团簇尺寸的增大而减小,而且表面能均高于体材料的表面能.     

Molecular dynamics study of effects of nickel coating on torsional behavior of single-walled carbon nanotube

The effects of nickel coating on the torsional behaviors of single-walled carbon nanotubes (SWCNTs) subject to torsion motion are investigated using the molecular dynamics (MD) simulation method. The simulation results show that regardless of chirality, defect or radius, nickel coating can considerably enhance the critical torque of SWCNTs. However, by comparing the critical torsion angle between nickel-coated SWCNTs and corresponding pristine SWCNTs, it is found that nickel coating in small-radius nanotubes does induce a reduction in the critical torsion angle. The results also show that the structural failure of nickel coated imperfect (9,0) SWCNT occurs at an obviously higher critical torque in comparison with uncoated (9,0) SWCNT with a vacancy defect. Furthermore, we also find that the critical torque of a short nickel coated SWCNT is bigger than that of a long one, while the critical torsion angle for a short tube is smaller than that for a long one.     

气体分子与纳米粒子碰撞的分子动力学仿真分析

当气体分子与纳米粒子碰撞的时候,纳米粒子传输理论预测到当纳米粒子的直径由小变大时,碰撞会由镜面反射转化为漫反射.文章利用分子动力学仿真研究了气体分子与纳米粒子碰撞的过程.在验证了这种转化存在同时,又探讨了碰撞转化的机理,即漫反射的起因.仿真结果揭示了漫反射的起因是由于纳米粒子对气体分子的吸附作用.这种吸附作用是由于纳米粒子对能量的容纳特性而产生的.     

气体分子与纳米粒子碰撞的分子动力学仿真分析

当气体分子与纳米粒子碰撞的时候,纳米粒子传输理论预测到当纳米粒子的直径由小变大时,碰撞会由镜面反射转化为漫反射.文章利用分子动力学仿真研究了气体分子与纳米粒子碰撞的过程.在验证了这种转化存在同时,又探讨了碰撞转化的机理,即漫反射的起因.仿真结果揭示了漫反射的起因是由于纳米粒子对气体分子的吸附作用.这种吸附作用是由于纳米粒子对能量的容纳特性而产生的.     

Molecular dynamics simulation of zirconia melting

The melting point for the tetragonal and cubic phases of zirconia (ZrO₂) was computed using Z-method microcanonical molecular dynamics simulations for two different interaction models: the empirical Lewis-Catlow potential versus the relatively new reactive force field (ReaxFF) model. While both models reproduce the stability of the cubic phase over the tetragonal phase at high temperatures, ReaxFF also gives approximately the correct melting point, around 2900 K, whereas the Lewis-Catlow estimate is above 6000 K.     

Molecular dynamics simulation of the sintering of metallic nanoparticles

The sintering of two different-sized nickel nanoparticles is simulated by a molecular dynamics method in this work. The particles are partitioned into different regimes where tracing atoms are arranged to investigate the sintering kinetics. The detailed sintering process of two nanoparticles, 3.52 and 1.76 nm in diameter, respectively, is subsequently examined by the shrinkage ratio, gyration radius, mean square displacement, sintering diffusivity, and activation energy. A three-stage sintering scenario is established, and the layered structure shows a regime dependent behavior of diffusivity during the sintering process. Besides the surface diffusion, sintering of different-sized nanoparticles is found to be affected by a few other mechanisms.     

Molecular dynamics simulation of the formation of metal nanocontacts

The process of the formation of nanocontacts has been studied by the molecular dynamics methods for a group of metals (Cu, Rh, Pd, Ag, Pt, Au). It has been shown that the disruption forces of nanocontacts substantially depend on the orientation ((100), (110), or (111)) of the contact-surface interface. The possibility of forming linear atomic chains as a result of the disruption of nanocontacts has been analyzed for different orientations of the electrode surfaces. The possibility of forming quasi-one-dimensional nanostructures from the Co/Au alloy, which represent a periodic alternation of gold atoms and cobalt trimers, has been predicted.     

Molecular dynamics simulation of cylindrical Richtmyer-Meshkov instability

The microscopic-scale Richtmyer-Meshkov(RM) instability of a single-mode Cu-He interface subjected to a cylindrically converging shock is studied through the classical molecular dynamics simulation. An unperturbed interface is first considered to examine the flow features in the convergent geometry, and notable distortions at the circular inhomogeneity are observed due to the atomic fluctuation. Detailed processes of the shock propagation and interface deformation for the single-mode interface impacted by a converging shock are clearly captured. Different from the macroscopic-scale situation, the intense molecular thermal motions in the present microscale flow introduce massive small wavelength perturbations at the single-mode interface, which later significantly impede the formation of the roll-up structure. Influences of the initial conditions including the initial amplitude,wave number and density ratio on the instability growth are carefully analyzed. It is found that the late-stage instability development for interfaces with a large perturbation does not depend on its initial amplitude any more. Surprisingly, as the wave number increases from 8 to 12, the growth rate after the reshock drops gradually. The distinct behaviors induced by the amplitude and wave number increments indicate that the present microscopic RM instability cannot be simply characterized by the amplitude over wavelength ratio(η). The pressure history at the convergence center shows that the first pressure peak caused by the shock focusing is insensitive to η, while the second one depends heavily on it.     

Molecular dynamics simulation of gas diffusion behavior in polyethylene terephthalate/aluminium/polyethylene interface

Molecular dynamics (MD) simulations were performed to estimate the diffusion coefficients of O₂ and H₂O molecules in polyethylene terephthalate/aluminum/polyethylene interface at the temperature of 298 K. It came out that the diffusion coefficient of gasses in the interface is smaller than that of a single polymer, and the diffusion coefficients compare well with experimental data as well as previously published work. Furthermore, the diffusion coefficients of H₂O molecules in the interface are preferable to that of O₂ molecules. Interestingly, the largest diffusion coefficient was detected in the polyethylene terephthalate/aluminum(1 0 0)/polyethylene interface, while the smallest value of the diffusion coefficients was found in the polyethylene terephthalate/aluminum(1 1 1)/polyethylene interface. Calculation and analysis of the interaction between aluminum and polymers indicated that the interaction of polymer/aluminum(1 1 0) has the most interface strength, and crystal density of the metal surface has a definite effect on the planar interface energy. What’s more, the figure of gas molecule concentration is further resulted that the interface make contribution to adsorption of gas molecules. Moreover, the diffusion is belonging to the Einstein diffusion in the multilayer materials, and this work provides some key clues to improve the performance of polymer materials.