Molecular dynamics study of the mechanical characteristics of Ni/Cu bilayer using nanoindentation

<正>In the present work,a three-dimensional molecular dynamics simulation is carried out to perform the nanoindentation experiment on Ni single crystal.The substrate indenter system is modeled using hybrid interatomic potentials including the many-body potential embedded atom method(EAM),and two-body morse potential.To simulate the indentation process,a spherical indenter(diameter = 80 A,1 A=0.1 nm) is chosen.The results show that the mechanical behaviour of a monolithic Ni is not affected by crystalline orientation.To elucidate the effect of a heterogeneous interface, three bilayer interface systems are constructed,namely Ni(100)/Cu(111),Ni(110)/Cu(111),and Ni(111)/Cu(111).The simulations along these systems clearly describe that mechanical behaviour directly depends on the lattice mismatch. The interface with the smaller mismatch between the specified crystal planes is proved to be harder and vice versa.To describe the relationship between film thickness and interface effect,we choose various values of film thickness ranging from 20 A to 50 A to perform the nanoindentation experiment.It is observed that the interface is significant only for the relatively small thickness of film and the separation between interface and the indenter tip.It is shown that with the increase in film thickness,the mechanical behaviour of the film shifts more toward that of monolithic material.     

纳米多晶金属样本构建的分子动力学模拟研究

研究了分子动力学模拟中纳米多晶金属样本的构建过程.首先采用Voronoi几何方法生成初始的纳米多晶铝和铜样本,然后用快速冷凝(或共轭梯度)法得到样本的局域最低能态,最后在恒温零应力周围环境下(常温常压NPT系综)退火得到最低能态样本.使用样本的残余内应力来衡量纳米多晶样本是否与实验制备的一致.通过监测这两步弛豫过程中晶界结构的变化形态、体系平均内应力和能量下降过程及具体的局域分布和不同弛豫条件下最终样本的弹性常数,发现样本的能量和残余内应力都接近实验制备的纳米多晶金属.对Voronoi几何法生成的晶界而言     

金属Ni纳米线凝固行为的分子动力学模拟

采用EAM镶嵌原子作用势，通过经典的分子动力学模拟方法研究了不同冷却速度下的金属Ni纳米线的凝固行为，并给出了纳米线在凝固区域的结构演变过程.利用键对分析技术研究了在不同冷却速度下体系中的原子团簇在降温过程中的变化情况.研究表明，纳米线的凝固起始于表面原子，并且随着冷却速度的降低，Ni纳米线的微观结构从非晶态过渡到多壳螺旋结构，最终达到稳定的面心立方结构.多壳螺旋结构同时具有确定的结晶温度和长程无序、短程有序的非晶结构的特征. 关键词： 纳米线 凝固行为 分子动力学 键对分析     

镍基单晶超合金Ni/Ni3Al晶界的分子动力学模拟

在镍基单晶超合金中,由于单晶Ni的晶格常数比单晶Ni3Al的稍小,在Ni/Ni3Al晶界面上必然要出现错配.采用分子动力学模拟了镍基单晶超合金的Ni/Ni3Al晶界的结构,考虑了两个不同的初始模型,并进行了分子动力学弛豫.弛豫的结果均表明:由于晶格的差异形成的错配能不是通过长程晶格错配的方式来释放,而是通过在局部区域形成位错的方式释放的.由于Ni3Al相周围Ni相环境的不同,形成的位错也有所不同.     

Brittle-ductile behavior of a nanocrack in nanocrystalline Ni: A quasicontinuum study

The effects of stacking fault energy, unstable stacking fault energy, and unstable twinning fault energy on the fracture behavior of nanocrystalline Ni are studied via quasicontinuum simulations. Two semi-empirical potentials for Ni are used to vary the values of these generalized planar fault energies. When the above three energies are reduced, a brittle-to-ductile transition of the fracture behavior is observed. In the model with higher generalized planar fault energies, a nanocrack proceeds along a grain boundary, while in the model with lower energies, the tip of the nanocrack becomes blunt. A greater twinning tendency is also observed in the more ductile model. These results indicate that the fracture toughness of nanocrystalline face-centered-cubic metals and alloys might be efficiently improved by controlling the generalized planar fault energies.     

深过冷液态Ni2 TiAl合金热物理性质的分子动力学模拟

采用嵌入原子法对深过冷条件下Ni2TiAl合金的焓和密度进行了分子动力学模拟,模拟获得了比热和密度在1200K至2000K范围内的温度依赖关系.结果表明过冷条件下Ni2TiAl合金的比热、密度和温度之间存在着线性关系,它们都随着温度的降低而降低.比热的大小与依据Neumann-Kopp法则估算的结果相近,同时模拟过程中的尺度效应经验证可以忽略.     

Dislocation structure and dynamics govern pop-in modes of nanoindentation on single-crystal metals

ABSTRACT

There are two types of pop-in mode that have been widely observed in nanoindentation experiments: the single pop-in, and the successive pop-in modes. Here we employ the molecular dynamics (MD) modelling to simulate nanoindentation for three face-centred cubic (FCC) metals, including Al, Cu and Ni, and two body-centred cubic (BCC) metals, such as Fe and Ta. We aim to examine the deformation mechanisms underlying these pop-in modes. Our simulation results indicate that the dislocation structures formed in single crystals during nanoindentation are mainly composed of half prismatic dislocation loops. These half prismatic dislocation loops in FCC metals are primarily constituted of extended dislocations. Lomer–Cottrell locks that result from the interactions between these extended dislocations can resist the slipping of half dislocation loops. These locks can build up the elastic energy that is needed to activate the nucleation of new half dislocation loops. A repetition of this sequence results in successive pop-in events in Al and other FCC metals. Conversely, the half prismatic dislocation loops that form in BCC metals after first pop-in are prone to slip into the bulk, which sustains plastic indentation process after first pop-in and prevents subsequent pop-ins. We thus conclude that pop-in modes are correlated with lattice structures during nanoindentation, regardless of their crystal orientations.     

Structural stability of nano-sized crystals of HMX: A molecular dynamics simulation study

The interaction potential energy and heat of sublimation of nanoparticles of HMX crystal polymorphs are studied by using molecular dynamics methods with a previously developed force field [Bedrov, et al., J. Comput.-Aided Mol. Des. 8 (2001) 77]. Molecular dynamics simulations of nanoparticles with 10, 20, 30, 40, 50, 60, 70, 80, 90, and 100 molecules of HMX are carried out at 300 K. The intermolecular, intramolecular and total interaction energies per mole for the nanoparticles are calculated at 300 K. Then, we have calculated sublimation enthalpy of HMX crystal polymorphs with different sizes. For the all sizes, the β-HMX is found to be the most stable phase, due to having the least total interaction energy. Also, α-HMX is more stable than δ-HMX. An increase in the sublimation enthalpy with the size of the nanoparticle can be seen.     

多环芳烃在超临界环己烷中的溶解：分子动力学模拟研究

重油中大量存在的多环芳烃会经过缩合反应会形成稠环芳烃,而溶解扩散的限制会导致稠环芳烃之间的聚合从形成残焦,使得对重油的利用率大大降低,所以多环芳烃热解过程的溶解扩散是重油改质的关键因素之一.因此,基于环己烷对烃类良好的溶解性,本文采用分子动力学模拟的方法研究了多环芳烃及其混合物在超临界环己烷中的溶解行为,结果表明在不同温度与密度下超临界环己烷对于具有NAP均具有良好的溶解性,而温度对于Bghip溶解的影响较小.在多环芳烃混合物的油滴溶解过程中,NAP优先溶解到环己烷相中,而Bghip则集中在粒径减小的油滴中.通过计算体系中多环芳烃的径向分配函数与溶剂化自由能发现环己烷与多环芳烃间之间能相互吸引,环己烷在多环芳烃分子周围形成的溶剂壳能抑制多环芳烃分子间的聚合.温度升高以及密度的将降低会导致多环芳烃的内聚能密度降低,增加了多环芳烃与超临界环己烷之间的相互作用.环己烷密度越小以及温度的升高可以促进多环芳烃或多环芳烃混合物在超临界环己烷中的溶解.     

Conformational properties of cyclooctane: a molecular dynamics simulation study

Atomistic molecular dynamics simulations have been used to elucidate the conformational properties of cyclooctane in the gas and bulk liquid phases. Accurate reproduction of the gas phase structure, and of the liquid phase densities and solubility parameters have been used as prerequisites to the prediction of conformational properties. The gas phase results clearly indicate the presence of a conformational mixture consisting of the crown, boat-chair, twist-boat-chair and boat-boat conformers at all temperatures (161, 313 and 400 K) studied. The fraction of the crown family of conformers was found to be relatively insensitive to temperature. However, the relative concentrations of the twist-boat-chair and boat-chair conformations was found to be highly temperature dependent with the boat-chair being favoured at low temperatures. Bulk packing was found to have a profound effect on the conformational properties in the liquid phase. At the temperatures studied (313 and 400 K) the boat-chair family was predominant, with the crown and boat families being essentially absent. The twist-boat-chair conformation was detected in the liquid phase at both temperatures. The pseudorotation pathway for the twist-boat-chair to boat-chair interconversion was prevalent in both gas and liquid phases establishing the conformational flexibility and the relative importance of the twist-boat-chair conformer in comparison to the crown family. The study successfully explains the separate experimental findings in both the gas and liquid phases of cyclooctane.     

外延生长薄膜中失配位错形成条件的分子动力学模拟研究

运用分子动力学方法对负失配条件下的外延铝簿膜中失配位错的形成进行了模拟研究.所采 用的原子间相互作用势为嵌入原子法(EAM)多体势.模拟结果显示：在500K下长时间静态弛豫 ，表面和内部结构完整的外延膜在9—80原子层厚度范围内(约为其热力学临界厚度的3—40 倍)均不形成失配位错，而在薄膜表面预置一个单原子层厚、三个原子直径大小的凸台或凹 坑时，失配位错则能够在15个原子层厚的外延膜上迅速形成：在动态沉积生长条件下，表面 自然形成凹凸，初始厚度为9个原子层厚的外延膜在沉积生长中迅速形成失配位错.在三种条 件下，所形成的位错均为伯格斯矢量与失配方向平行的全刃位错.分析发现：在压应力作用 下，表面微凸台诱发了其侧薄膜内部原子的挤出，造成位错形核；而表面微凹坑则直接因压 应力作用形成了一个表面半位错环核. 关键词： 外延薄膜 失配位错 分子动力学 铝     

A study of cavitation nucleation in pure water using molecular dynamics simulation

To discover the microscopic mechanism responsible for cavitation nucleation in pure water, nucleation processes in pure water are simulated using the molecular dynamics method. Cavitation nucleation is generated by uniformly stretching the system under isothermal conditions, and the formation and development of cavitation nuclei are simulated and discussed at the molecular level. The processes of energy, pressure, and density are analyzed, and the tensile strength of the pure water and the critical volume of the bubble nuclei are investigated. The results show that critical states exist in the process of cavitation nucleation. In the critical state, the energy, density, and pressure of the system change abruptly, and a stable cavitation nucleus is produced if the energy barrier is broken and the critical volume is exceeded. System pressure and water density are the key factors in the generation of cavitation nuclei. When the critical state is surpassed, the liquid is completely ruptured, and the volume of the cavitation nucleus rapidly increases to larger than 100 nm³; at this point, the surface tension of the bubble dominates the cavitation nucleus, instead of intermolecular forces. The negative critical pressure for bubble nucleation is -198.6 MPa, the corresponding critical volume is 13.84 nm³, and the nucleation rate is 2.42×10³² m^-3·^-1 in pure water at 300 K. Temperature has a significant effect on nucleation: as the temperature rises, nucleation thresholds decrease, and cavitation nucleation occurs earlier.     

纳米多晶铜微观结构的分子动力学模拟

基于EAM镶嵌原子势函数,采用分子动力学模拟了零温下纳米多晶铜的微观结构.首先用Voronoi几何方法构造了5个纳米多晶铜数值模型,在300K弛豫50ps并退火至0K.然后分析零温下弛豫模型的径向分布函数、原子能量、配位数、原子Voronoi体积、以及本征应力分布 关键词： 纳米多晶铜 微观结构 分子动力学     

Femtosecond laser ablation of metals:a molecular dynamics simulation study

Using molecular dynamics (MD) methods combining with two-step radiation heating model, the mechanisms of ablation and the thermodynamic states at Ni surface under femtosecond laser irradiation are investigated. Simulation results show that the main mechanisms of ablation are evaporation and tensile stresses generated inside the target. The velocity of stress wave is predicted to be nearly equal to sound velocity. The rates of ablation at different fluences obtained from simulations are in good agreement with experimental data. Superheating phenomenon is also discovered.     

Relaxation dynamics and power spectra of liquid water: a molecular dynamics simulation study

We present molecular dynamics simulations of liquid water at normal and supercooled conditions. Autocorrelation functions (ACFs) of several structural quantities and their fourier transforms are obtained and analysed. Structural correlations and relaxation times increase linearly with degree of supercooling. Power spectra of ACFs show increase in librational motion of liquid water with cooling. These modes intensify with supercooling because of structuring and ordering of water molecules. Overall, liquid water structure is homogenous over the temperatures and pressures studied and undergoes fluctuation–dissipation in its local-density variations [English and Tse, Phys. Rev. Lett. 106, 037801 (2011)].     

Anisotropic plasticity of nanocrystalline Ti:A molecular dynamics simulation

Using molecular dynamics simulations,the plastic deformation behavior of nanocrytalline Ti has been investigated under tension and compression normal to the{0001},{1010},and{1210}planes.The results indicate that the plastic deformation strongly depends on crystal orientation and loading directions.Under tension normal to basal plane,the deformation mechanism is mainly the grain reorientation and the subsequent deformation twinning.Under compression,the transformation of hexagonal-close packed(HCP)-Ti to face-centered cubic(FCC)-Ti dominates the deformation.When loading is normal to the prismatic planes(both{1010}and{1210}),the deformation mechanism is primarily the phase transformation among HCP,body-centered cubic(BCC),and FCC structures,regardless of loading mode.The orientation relations(OR)of{0001}HCP||{111}FCC and<1210>HCP||<110>FCC,and{1010}HCP||{110}FCC and<0001>HCP||<010>FCC between the HCP and FCC phases have been observed in the present work.For the transformation of HCP→BCC→HCP,the OR is{0001}α1||{110}β||{1010}α2(HCP phase before the critical strain is defined as α1-Ti,BCC phase is defined as β-Ti,and the HCP phase after the critical strain is defined as α2-Ti).Energy evolution during the various loading processes further shows the plastic anisotropy of nanocrystalline Ti is determined by the stacking order of the atoms.The results in the present work will promote the in-depth study of the plastic deformation mechanism of HCP materials.     

A technique to study doped ablation in polymethyl methacrylate using molecular dynamics simulation

A doped polymer substrate is studied using molecular dynamics simulations. Photons are preferentially absorbed as thermal energy by a spherical portion of a polymethyl methacrylate sample. Once sufficient energy is absorbed to fragment the polymer-labeled cluster and cleave the residual polymer bonds, ejection of material occurs indicating ablation. The mechanism of ejection is analogous to our previous work where a dopant cluster was represented using a cluster of carbon atoms [P.F. Conforti, M. Prasad, B.J. Garrison, Physical Chemistry Chemical Physics, 10 (2008) 6002]. This labeling procedure can serve as an approximate method to qualitatively explore the effects of different parameters without becoming consumed in the details of the simulation setup procedure.     

Generation of laser-driven flyer dominated by shock-induced shear bands: A molecular dynamics simulation study

High-power pulsed lasers provide an ingenious method for launching metal foils to generate high-speed flyers for high-pressure loading in material science or aerospace engineering.At high-temperature and high-pressure laser-induced conditions,the dynamic response of the metals and the mechanism of flyer formation remain unclear.In this study,the overall process of the laser-driven aluminum flyer,including laser ablation,rupture of metal foil,and the generation of the flyer was investigated by molecular dynamics combined with the two-temperature model.It was found that under high laser fluence(over 1.3 J/cm;with 200-fs laser pulse duration),the laser induced a shock wave with a peak pressure higher than25 GPa,which led to shear bands expanding from the edge of the laser ablation zone in the foil.Compared with the cases of low laser fluence less than 0.5 J/cm^-1,the shear band induced by high laser fluence promotes the rupture of the foil and results in a high-speed flyer(>1 km/s)with better flatness and integrity.In addition,the shock wavefront was found to be accompanied by aluminum crystal phase transformation from face-centered cubic(FCC)to body-centered cubic structure.The crystal structure reverts with the decrease of pressure,therefore the internal structure of the generated flyer is pure of FCC.The results of this study provide a better understanding of the laser-induced shock effect on the foil rupture and flyer quality and forward the development of the laser-driven flyer.     

立方铂纳米粒子的形状变化与熔化特性的分子动力学研究

采用分子动力学方法结合嵌入原子多体势,对立方铂纳米粒子的热稳定性进行了模拟研究.计算结果表明,立方纳米粒子在升温过程中首先转变为由{111}和{100}面所构成的十四面体,然后再转变为球形,最后熔化为液态.通过计算立方铂钠米粒子的统计半径,发现形状转变温度在1250 K左右.尽管形状不同,立方纳米粒子和球形纳米粒子的熔点是相同的. 关键词： 纳米粒子 热稳定性 分子动力学