Study of AFM-based nanometric cutting process using molecular dynamics

Three-dimensional molecular dynamics (MD) simulations are conducted to investigate the atomic force microscope (AFM)-based nanometric cutting process of copper using diamond tool. The effects of tool geometry, cutting depth, cutting velocity and bulk temperature are studied. It is found that the tool geometry has a significant effect on the cutting resistance. The friction coefficient (cutting resistance) on the nanoscale decreases with the increase of tool angle as predicted by the macroscale theory. However, the friction coefficients on the nanoscale are bigger than those on the macroscale. The simulation results show that a bigger cutting depth results in more material deformation and larger chip volume, thus leading to bigger cutting force and bigger normal force. It is also observed that a higher cutting velocity results in a larger chip volume in front of the tool and bigger cutting force and normal force. The chip volume in front of the tool increases while the cutting force and normal force decrease with the increase of bulk temperature.     

Molecular dynamics study of groove fabrication process using AFM-based nanometric cutting technique

Three-dimensional molecular dynamics simulation of groove fabrication using atomic force microscopy (AFM)-based nanometric cutting technique is set up, fabrication processes of grooves with two types (line, and folder line) and five folder angles (0°, 30°, 45°, 60°, 90°) are simulated to investigate the effect of groove geometry on the fabrication process. The results show that the Normal force, Lateral force, and Resultant forces are almost symmetric with respect to the critical folder angle of 45°. The best surface quality of fabricated groove can be obtained at the folder angle of 45°. It reveals that the groove geometry has a significant effect on the groove fabrication process due to the material anisotropy on the atomic scale.     

导电石墨烯孔道内双电层结构的研究

为有效开发和利用新能源,人们迫切需要高性能的超级电容器提供能量的存储和转换.在超级电容器中双电层结构扮演着关键性的角色.本文利用分子动力学方法通过建立开放的石墨烯纳米孔道(1～2 nm),研究了KCl溶液在纳米孔道内的双电层结构,同时也比较了恒电量模拟(Q)和恒电势模拟法(U)下双电层结构的异同.结果表明在恒电势模拟法考虑了导电石墨烯壁的镜像作用使结果更符合实验中的材料系统.而石墨烯壁的镜像作用能额外吸附离子从而增强孔道内部的阴阳离子,这可能有助于电极电容的提升.通过对不同孔道高度的研究,本文发现水分子作为介电材料在水基超级电容器中发挥着决定性的作用.它能在很大程度上抵消不同离子和不同孔道高度下双电层的变化,从而在不同情况下获得了相似的电容.     

Quasicontinuum analysis of the effect of tool geometry on nanometric cutting of single crystal copper

A dynamic multiscale simulation based on quasicontinuum method (QC) has been conducted to study the effect of tool geometry in nanometric cutting process of single crystal copper. In the simulation, the many-body EAM potential is used for the interactions between copper atoms in of the workpiece. The simulation captures the atomistic behaviors of material removal mechanisms from the free surface and the mobility of dislocations and their interactions with the computational cost of local atomistic simulation method. Simulations are performed on single crystal copper to study the atomistic details of material removal, chip formation, sub-surface deformation, and machining mechanism. The simulation results demonstrate that tool edge radius has significant effect on chip formation and subsurface deformation, because the effective rake angle varies with the tool edge radius. In addition, different effective rake angles result in different stress states and smoother surface can be obtained under bigger clearance angle. The variations of tangential force, normal force as well as the ratio of normal force to tangential force are obtained to analyze the effects of tool edge radius, rake angle and clearance angle in quantitative way.     

Molecular dynamics simulation of pervaporation in zeolite membranes

The pervaporation separation of liquid mixtures of water/ethanol and water/methanol using three zeolite (Silicalite, NaA and Chabazite) membranes has been examined using the method of molecular dynamics. The main goal of this study was to identify intermolecular interactions between water, methanol, ethanol and the membrane surface that play a critical role in the separations. This would then allow better membranes to be designed more efficiently and systematically than the trial-and-error procedures often being used. Our simulations correctly exhibited all the qualitative experimental observations for these systems, including the hydrophobic or hydrophilic behaviour of zeolite membranes. The simulations showed that, for Silicalite zeolite, the separation is strongly influenced by the selective adsorption of ethanol. The separation factor, as a consequence, increases almost exponentially as the ethanol composition decreases. For ethanol dehydration in NaA and Chabazite, pore size was found to play a very important role in the separation; very high separation factors were therefore possible. Simulations were also used to investigate the effect of pore structure, feed compositions and operating conditions on the pervaporation efficiency. Finally, our simulations also demonstrated that molecular simulations could serve as a useful screening tool to determine the suitability of a membrane for potential pervaporation separation applications. Simulations can cost only a small fraction of an experiment, and can therefore be used to design experiments most likely to be successful.     

基于桥域理论的Cu单晶纳米切削跨尺度仿真研究

桥域方法是一种典型的跨尺度仿真研究方法.基于桥域理论,本文分析了原子和连续介质耦合区域的处理问题,即在耦合区采用不同的权重计算系统的能量,通过Lagrange乘子法对原子和连续介质位移进行约束.采用桥域方法,建立了单晶Cu米纳切削的跨尺度仿真模型,获得了单晶Cu纳米切削的材料变形机理.同时,研究了不同切削速度对纳米切削过程和原子受力分布的影响,仿真结果表明:随着切削速度的提高,切削区原子所受的力值增大,切屑变形系数减小,已加工表面变质层厚度增加.本文基于桥域理论,实现了Cu单晶纳米切削跨尺度的建模和仿真, 关键词： 桥域法 纳米切削 单晶Cu 切削速度     

Molecular dynamics simulation study on mechanical responses of nanoindented monolayer-graphene-nanoribbon

We investigated the mechanical responses of the nanoindented graphene-nanoribbon (GNR)-resonator using classical molecular dynamics simulations. The nanoindented force in this work was applied to the GNR's local point and then, GNR-resonator's frequency could be tuned by a nanoindented depth. We found the hardening or the softening of the GNR during its nanoindented-deflections, and such properties were recognized by the shift of the resonance frequency. The linear elastic regime in low applied force is explicitly separated with the non-linear elastic regime in high applied force. In particular, at the threshold point, a very small change of the nanoindented depth can cause great change in the resonance frequency, and this property can enable the GNR to be applied to electromechanical relay switching devices and the quantum-computer in quantum-mechanical coupling as well as mass detectors, pressure sensors, accelerometers, and alarms.     

Ni-Co合金快淬过程的分子动力学模拟

采用分子动力学模拟熔体旋淬技术的合金快淬过程,在不同的冷却速度下研究Ni-Co合金在快淬后的结构特征.模拟发现:Ni-Co二元合金的凝固过程对冷却速率的变化较为敏感,体系对形成非晶对冷却速率要求较高,约在80 K/ps以上,可以采用增加添加元素的方法来降低材料对冷却速率的要求;在75 K/ps的冷却速率下合金最终完全晶...     

Molecular dynamics simulation of defects production in the vicinity of voids

We report on the results of computer simulation of point defect production near voids in crystalline Cu at primary knock-on atom (PKA) kinetic energies ranging from 5 to 1000?eV. The PKA energy dependence of numbers of created defects are revealed. The threshold energy for a stable vacancy formation is found to be much smaller than that for an interstitial atom, which results in a biased formation of vacancies in the void proximity in the whole investigated range of PKA energies. Dissolution of small voids by subthreshold irradiation is simulated. The impact of considered radiation effects on kinetics of radiation damage is discussed.     

Ni-Co合金快淬过程的分子动力学模拟

采用分子动力学模拟熔体旋淬技术的合金快淬过程,在不同的冷却速度下研究Ni-Co合金在快淬后的结构特征。模拟发现:Ni-Co二元合金的凝固过程对冷却速率的变化较为敏感,体系对形成非晶对冷却速率要求较高,约在80 K/ps以上,可以采用增加添加元素的方法来降低材料对冷却速率的要求;在75 K/ps的冷却速率下合金最终完全晶化,发生结构转变,而且在较低温度(375 K)下很明显。     

CO_2水合物导热特性的分子动力学模拟

使用分子动力学模拟方法在NVT系综下对结构完整CO_2水合物以及结构缺陷CO_2水合物进行了导热模拟计算.对于结构完整的CO_2水合物,在200-230 K温度区间内,体系导热系数由0.4684 W·m~(-1)·K~(-1)变化到0.4836 W·m~(-1)·K~(-1),温度相关性较弱;而在230-280 K温度区间内,体系导热系数由0.4836 W·m~(-1)·K~(-1)变化到0.7494 W·m~(-1)·K~(-1),温度相关性变强;另外,通过计算功率图谱发现主体分子对水合物体系的导热贡献更大.对于结构缺陷CO_2水合物,发现晶穴占有率和笼形结构缺陷对体系导热均有一定影响,空笼晶胞导热系数约为完整晶胞导热系数的86.67%,体系的导热能力主要取决于主体结构的性质.     

Molecular dynamics simulation study of deep hole drilling in iron by ultrashort laser pulses

Classical molecular dynamics simulation technique is applied for investigation of the iron ablation by ultrashort laser pulses at conditions of deep hole for the first time. Laser pulse duration of 0.1 ps at wavelength of 800 nm is considered. The evolution of the ablated material in deep hole geometry differs completely from the free expansion regime as two major mechanisms are important for the final hole shape. The first one is the deposition of the ablated material on the walls, which narrows the hole at a certain height above its bottom. The second mechanism is related to ablation of the material from the walls (secondary ablation) caused by its interaction with the primary ablated particles. Properties of the secondary ablated particles in terms of the velocity and the angular distribution are obtained. The material removal efficiency is estimated for vacuum or in Ar environment conditions. In the latter case, the existence of well-defined vapor cloud having low center of the mass velocity is found. The processes observed affect significantly the material expulsion and can explain the decrease of the drilling rate with the hole depth increase, an effect observed experimentally.     

考虑非傅立叶效应的固液相变分子动力学模拟

采用耦合双温度模型的分子动力学方法对飞秒激光照射金箔的固液相变过程进行了模拟研究,利用序参数法对固液原子进行判定从而确定了金箔发生相变时的固液界面位置和温度,对基于傅立叶定律的抛物线模型和考虑非傅立叶效应的双曲线模型模拟得到的结果进行对比分析,在此基础上采用耦合双曲线模型的分子动力学方法研究了激光能量密度和脉冲宽度对金箔相变过程的影响.结果表明,当激光作用于金箔时,金箔上表面首先熔化,固液界面随时间不断向金箔底部移动,并且在相同条件下,双曲线模型下的金箔熔化深度和固液界面温度均大于抛物线模型的结果.当考虑非傅里叶效应时,激光能量密度越大,固液界面温度越高,金箔熔化时间越短;激光脉冲宽度越小,固液界面温度越大,金箔熔化速度越快.     

C+轰击铍的分子动力学模拟

使用分子动力学模拟方法研究了入射能量对C⁺离子与Be样品表面相互作用的影响。模拟结果表明,随着C⁺离子的入射能量增大,C⁺离子注入深度也增加,Be原子的溅射产额近似线性增加,而滞留在样品中的C原子数量变化不大,在C⁺离子轰击Be样品的初始阶段,样品中Be原子的溅射产额较大,而随着C⁺离子注入剂量的增加,Be原子的溅射产额逐渐减小并趋于稳定。在此作用过程中,在样品表面形成一个富C层,减缓了样品中Be原子的溅射速率,起到了保护Be样品的作用。     

Molecular dynamics simulations of atomic assembly in the process of GaN film growth

Molecular dynamics simulations using a Coulomb–Buckingham potential have been used to investigate the process of wurtzite GaN films growth, including the appearance of films in early stage, regulation of growth, structure of the surface and the dynamic features. The simulations show that the N atoms and Ga atoms are absorbed on the lattice of substrate and take on a distinct sandwich structure. Time evolution of the mean square displacements and diffusion coefficient of the deposited atoms are observed, the results show that the clusters will become stable with the increase of time steps and the atoms reach the initial stable state after 25 ps; N atoms reach the equilibrium positions more quickly than Ga atoms. It is proved by radial distribution function and the ratio of vacancy of every deposited layer that the crystalline characters of the films will become better as the time steps increase and weaker from bottom to top.     

纳米MgO掺杂聚乙烯微观特性的分子动力学模拟研究

电力系统高压电缆的主要绝缘材料为聚乙烯,为了提升聚乙烯的热稳定性以及减弱水分对其的渗透能力,采用纳米MgO掺杂聚乙烯,利用分子动力学模拟方法建立包含低密度聚乙烯(LDPE)、不同颗粒半径的MgO纳米团簇以及相同质量分数水分的复合模拟模型.研究结果表明,水分会降低复合体系的玻璃化温度,MgO的掺杂则会提高复合体系的玻璃化温度,减弱聚乙烯分子链的运动并减小复合体系的自由体积,使得复合体系结构更加稳定,从而增强了聚乙烯材料的热稳定性能.此外发现水分子的扩散随着温度的上升而增大,纳米MgO的添加会与水分子形成氢键抑制水分子的扩散,同时自由体积的缩减使水分子的溶解度系数与扩散系数都减小,导致水分子的渗透能力减弱,更难以渗透进聚乙烯材料破坏其结构.研究结果可为聚乙烯的水树枝生长以及老化过程的抑制提供有益的参考.     

超短激光诱导金属薄膜后向层裂的分子动力学模拟

 采用结合双温模型的分子动力学方法详尽描述了应力约束区域内部金属薄膜后向层裂的动力学过程。与辐照表面在激光加热作用下机械稳定性受到强烈影响而发生的前向喷射不同，后向层裂是冷材料的断裂。分析了层裂机制，得出靶材是在卸载波及被反射的压力波的共同作用下发生层裂;探讨了激光诱导压力波的传播规律，预测了不同靶厚下的层裂厚度及其对层裂开始时间的影响。     

氯化钾水溶液微观结构的分子动力学模拟研究

氯化钾水溶液不同于氯化钠的,其在低温下只能形成无水盐和冰.现有研究表明在氯化钠水溶液中不仅存在着水分子连续分布的区域与离子团簇,而且这二者的物质组成也与该溶液所形成晶体的分别对应.为了探寻其他溶液中是否存在此类对应关系,并考察溶液的微观结构,本文采用分子动力学方法对氯化钾水溶液进行了研究,表明了该溶液的一些性质.氯化钾水溶液中K⁺-K⁺和Cl^--Cl^-径向分布函数的特征具有一致性,峰的最大值所对应的位置都分别相同,明显不同于氯化钠水溶液的.系列时刻下瞬态图像内O到其最近离子距离中最大值的统计结果表明氯化钾水溶液中存在着一定大小的水分子连续分布的区域,其平均尺寸至少为2.26 nm;瞬态图像中K⁺与其最近邻Cl^-之间的距离主要分布在0.28 nm～0.38 nm之间,占比约为97.4%;溶液中存在着较大和较小两类离子团簇,较大团簇的平均尺寸为1.73nm,平均离子数是25.0,其内部的离子与周围离子之间具有与氯化钾晶体类似的结构;这些结果表明氯化钾水溶液中也...     

单晶硅纳米切削中C–C键断裂对金刚石刀具磨损的影响

本文基于分子动力学方法模拟金刚石刀具纳米切削单晶硅, 从刀具的弹塑性变形、C–C键断裂对碳原子结构的影响以及金刚石刀具的石墨化磨损等方面对金刚石刀具的磨损进行分析, 采用配位数法和6元环法表征刀具上的磨损碳原子. 模拟结果表明: 在纳米切削过程中, 金刚石刀具表层C–C键的断裂使其两端碳原子由sp³杂化转变为sp²杂化, 同时, 表面上的杂化结构发生变化的碳原子与其第一近邻的sp²杂化碳原子所构成的区域发生平整, 由金刚石的立体网状结构转变为石墨的平面结构, 导致金刚石刀具发生磨损; 刀具表面低配位数碳原子的重构使其近邻区域产生扭曲变形, C–C键键能随之减弱, 在高温和高剪切应力的作用下, 极易发生断裂; 在切削刃的棱边上, 由于表面碳原子的配位严重不足, 断开较少的C–C键就可以使表面6 元环中碳原子的配位数都小于4, 导致金刚石刀具发生石墨化磨损.