非晶态合金表面的水润湿动力学

采用分子动力学模拟方法研究了改进的Simple pointcharge模型SPC/E水滴在Cu_(50)Zr_(50)非晶薄膜上的润湿行为和铺展过程.通过与CuZr(110)和(100)晶面对比研究发现,水滴在Cu_(50)Zr_(50)非晶薄膜表面上表现出较高的铺展速度.水滴在非晶合金表面的铺展过程中形成了明显的吸附层;而在晶态表面,水滴铺展前沿呈脚状形态.分析结果表明非晶表面的水分子在吸附层内呈现完全无序的单层排列方式,而在晶态表面,特别是(100)晶面,吸附层水分子呈双层有序排列.这种吸附层结构的差异导致了吸附层内水分子方向的差异:非晶表面吸附层内水分子方向倾向平行于表面,而晶态基底上吸附层内的水分子倾向于指向液滴内部.由此造成了非晶表面吸附层中的水分子与液滴内部以相对较弱的氢键相互作用,这使得上层水分子能够较容易扩散至吸附层前沿,促进液滴铺展.     

fcc金属表面能的各向异性分析及表面偏析的预测

本文将元素变量(φ^*和n_WS)和MAEAM相结合,从原子尺度上对10种fcc金属Cu,Ag,Au,Ni,Pd,Pt,Rh,Al,Ir和Pb的38个不同晶面的表面能进行模拟计算及各向异性分析. 结果表明,fcc金属的密排面(111)的表面能最小,则该晶粒取向优先生长,与实验结果和第一原理的LMTO-ASA计算结果一致;各个晶面的表面能均随着其他晶面与(111)晶面的夹角cosθ_(hkl)的增长而呈线性 关键词： FCC金属 MAEAM 表面能 表面偏析     

H2O分子在Fe（100）, Fe（110）, Fe（111）表面吸附的第一性原理研究

采用第一性原理研究了H₂O分子在Fe（100）,Fe（110）,Fe（111）三个高对称晶面上的表面吸附.结果表明,H₂O分子在三个晶面上的最稳定结构皆为平行于基底表面的顶位吸附结构.H₂O分子与三个晶面相互作用的吸附能及几何结构计算结果表明H₂O分子与三个晶面的相互作用程度不同,H₂O分子与Fe（111）晶面的相互作用最强,其次是Fe（100）,相互作用最弱的是Fe（110）表面,而这与晶面原子 关键词： 第一性原理 Fe单晶表面 2O分子')" href="#">H₂O分子 分子吸附     

UO2 晶体中低密勒指数晶面表面能的分子动力学模拟

采用基于刚性离子势的分子动力学模拟方法初步计算了UO₂晶体中(100), (110)和(111) 3种低密勒指数晶面在300–1500 K范围内的表面能大小. 结果表明, 3种晶面的表面能大小随温度的升高而降低, 与实验结果趋势一致; 原子排列最紧密的(111)晶面具有最低的表面能, 3种晶面的表面能大小从高到低依次为(100), (110)和(111)晶面; 达到平衡状态下的表面层原子相对于体内原子层在表面的法线方向上发生了明显的压缩并且表面层原子的对称性也降低了, 表面原子的弛豫效应一直影响到了第5层. 计算研究结果将有助于深入认识UO₂燃料中裂变气体气泡的聚集长大以及燃料的辐照肿胀开裂行为. 关键词： 分子动力学 2')" href="#">UO₂ 低密勒指数晶面 表面能     

H2O分子在Fe(100), Fe(110), Fe(111)表面吸附的第一性原理研究

采用第一性原理研究了H2O分子在Fe(100),Fe(110),Fe(111)三个高对称晶面上的表面吸附.结果表明,H2O分子在三个晶面上的最稳定结构皆为平行于基底表面的顶位吸附结构.H2O分子与三个晶面相互作用的吸附能及几何结构计算结果表明H2O分子与三个晶面的相互作用程度不同,H2O分子与Fe(111)晶面的相互作用最强,其次是Fe(100),相互作用最弱的是Fe(110)表面,而这与晶面原子的排列密度相关.吸附体系的电子结构计算结果也得出了相似的结论.同时电荷布居分析表明,H2O分子与Fe表面相互作用时,O原子与基底原子之间的电荷交换使基底Fe原子表面带负电,导致表面电位降低,也促使Fe表面更易于发生电化学腐蚀反应.     

分离压作用下波纹状基底上含活性剂液滴铺展稳定性

针对波纹状基底上含不溶性活性剂液滴的铺展过程,引入受活性剂浓度影响的分离压模型,应用非模态稳定性理论分析液滴铺展稳定性.研究表明:与不计及分离压情形相比,分离压作用下的液滴铺展前沿高度明显下降,液滴铺展速率加快;长波扰动有利于增强液滴的演化稳定性,且随扰动波数增大稳定性增强;然而随短波波数增加,液滴演化稳定性逐渐减弱甚至转变为不稳定;较小波数下(k=3)减小引力强度系数α₁和提高斥力强度系数α₂液滴铺展稳定性增强;而较大波数下(k=30)增大α₁和减小α₂有利于液滴铺展稳定性.     

碱基分子的原子电距矢量表达及核磁共振碳谱模拟

提出以原子电性距离矢量(VAED)描述各种碱基分子中不同等价碳原子的化学环境,结合γ效应校正与碳原子类型,建立核磁共振碳谱(¹³C NMR)化学位移(CS)的五参数线性模型.用于碱基分子中4类不同碳(其中主要是仲叔季碳)的等价碳原子化学位移的估计,复相关系数R和均方根误差RMS[ppm]分别为仲叔季碳:1)对C2n=12,R=0.9998,R²=0.9997,F=2732.2498,EV=0.9994,RMS=0.1804,SD=0.18842)对C3n=41,R=0.9334,R²=0.8712,F=38.3299,EV=0.8528,RMS=9.5268,SD=9.64523)对C4n=23,R=0.9183,R²=0.8433,F=14.3479,EV=0.7972,RMS=6.5304,SD=6.67721)对C2n=12,R=0.9998,R²=0.9997,F=2592.7072,EV=0.9994,RMS=0.1852,SD=0.19342)对C3n=41,R=0.9270,R²=0.8594,F=34.6283,EV=0.8393,RMS=9.548,SD=10.07843)对C4n=23,R=0.9071,R²=0.8229,F=12.3881,EV=0.7708,RMS=6.9425,SD=7.0985经交互校验,模拟稳定性较好.并综合几种处理方法,找到一种较好的建模方法,将它用于4个外部样本化学位移的定量预测,结果良好.     

层流圆管潜射流生成蘑菇形涡结构特性数值研究

基于不可压Navier-Stokes方程, 采用计算流体力学方法, 数值模拟与分析了层流圆管潜射流在密度均匀黏性流体中的演化机理及其表现特征, 定量研究了蘑菇形涡结构无量纲射流长度L^*、螺旋型涡环半径R^*及其包络外形长度d^*等几 何特征参数随无量纲时间t^*的变化规律. 数值结果表明, 蘑菇形涡结构的形成与演化过程可分为三个不同的阶段: 启动阶段、发展阶段和衰退阶段. 在启动阶段, L^*和d^*随t^* 线性变化, 而R^*则近似为一个常数; 在发展阶段, 蘑菇形涡结构的演化具有自相似性, L^*, R^*和d^*与t^*1/2均为同一正比关系, 而且雷诺数和无量纲射流时间不影响该正比关系; 在衰退阶段, L^* 和R^* 正比于t^*1/5, 而d^*则近似为一个常数. 此外, 还对蘑菇形涡结构二次回流点、 动量源作用中心及其几何中心的速度变化规律、垂向涡量分布特征和 涡量-流函数关系进行了分析. 关键词： 圆管潜射流 蘑菇形涡结构 演化机理     

波动方程辛算法的迭代求解

对(∂²)/(∂x²)利用中心差商算子,对expt作对角Padé逼近,由波动偏微分方程可得到两类具有O(Δx²+Δt^2l)和O(Δx⁴+Δt^2l)精度的辛格式.对由此类辛格式产生的线性方程组构造了两种迭代解法,并对l=1,2,3,4给出了它们的收敛条件.并进行了数值实验.     

不同晶面银纳米晶高温熔化的各向异性

采用嵌入原子势, 使用分子动力学方法, 模拟研究了银纳米晶高温弛豫过程中的热稳定性和熔化机制, 并引入均方位移和稳定寿命来分析它的结构和形状的演化过程. 结果表明: 对于沿相互垂直{110}, {211}和{111}面切割形成的近正方体截面纳米晶, 高温弛豫熔化存在明显的各向异性行为; (112) 面热稳定性最低, 最易熔化, 其次是(110) 面, 热稳定性最高的是(111) 面, 最难熔化; 三个不同晶面的最外层和次外层原子的稳定寿命极短, 且三个不同晶面之间相差很小, 没有明显差异; 对于具有相同晶面指数的晶面, 第三层及其以内的稳定寿命较长, 且依次微量增长, 但不同晶面第三层及其以内的寿命相差明显.     

Precursor film controlled wetting of Pb on Cu

Wetting in a system where the kinetics of drop spreading are controlled by the rate of formation of a precursor film is modeled for the first time at the atomistic scale. Molecular dynamics simulations of Pb(l) wetting Cu(111) and Cu(100) show that a precursor film of atomic thickness evolves and spreads diffusively. This precursor film spreads significantly faster on Cu(111) than on Cu(100). For Cu(100), the kinetics of drop spreading are dramatically decreased by slow advancement of the precursor film. Slow precursor film kinetics on Cu(100) are partly due to the formation of a surface alloy at the solid-liquid interface which does not occur on Cu(111).     

The definition and calculation of interfacial energies for thin films

We provide a new definition of the interfacial energy which eliminates three physically extraneous contributions from the conventional definition: (1) the strain or stress energy due to lattice mismatch between film and substrate; (2) the surface energy of the film-vacuum interface; and, (3) the substrate surface energy contribution from substrate layers below the film layers. This new interface energy then quantifies the variation in interactions among film/substrate, film/film and substrate/substrate bonding. Using this new definition, we derive the equations for evaluation of the interfacial energy in terms of the interaction energy for any atom in each layer of the film/substrate, film/film and substrate/substrate systems. With this formulation, it is simple to determine the dependence of the interfacial energy on the film thickness using virtually any interaction potential. Using a corrected effective medium theory, we present results for a few pseudomorphic film systems containing Ni/Cu, Ni/Ag, Cu/Ag and Rh/Ag on (111) and (100) surfaces. These systems cover a wide range of lattice mismatch and alloy formation energies. The results demonstrate that the new definition of interfacial energy converges after only 3–4 film layers, regardless of the degree of lattice mismatch. We also show that the interfacial energies at (100) and (111) interfaces differ and that the interfacial energy for a given pair of materials depends on which of the materials is the film.     

Surface wetting of liquid nanodroplets: droplet-size effects

The spreading of liquid nanodroplets of different initial radii R0 is studied using molecular dynamics simulation. Results for two distinct systems, Pb on Cu(111), which is nonwetting, and a coarse-grained polymer model, which wets the surface, are presented for Pb droplets ranging in size from approximately 55,000 to 220,000 atoms and polymer droplets ranging in size from approximately 200,000 to 780 000 monomers. In both cases, a precursor foot precedes the spreading of the main droplet. This precursor foot spreads as r(2)(f)(t) = 2D(eff)t with an effective diffusion constant that exhibits a droplet-size dependence D(eff) approximately R(1/2)(0). The radius of the main droplet r(b)(t) approximately R(4/5)(0) is in agreement with kinetic models for the cylindrical geometry studied.     

固体表面液滴铺展与润湿接触线的移动分析

液滴在固体表面上的铺展行为与润湿特性对许多工业生产过程的研究具有重要意义.根据液滴在光滑表面上的受力情况,建立了液滴平壁铺展的动力学模型.应用润滑近似方法和二维Navier-Stokes方程,建立了液滴沿理想表面铺展的动量和连续性方程.根据建立的方程,应用数值解法求解并详细分析了液滴在铺展过程中膜厚、接触线铺展半径以及铺展速度随时间的变化关系.研究结果表明:液滴的铺展过程可分为扩展和收缩两个阶段,铺展过程伴随着表面能、动能以及各种势能的相互转化,液滴最终的铺展半径大小由固体基面固有的润湿特性所决定;液滴在铺展过程中出现的"坍塌效应"与弯曲液面处的Laplace压力差有关;铺展半径随时间变化的标定律近似满足"1/7"次方标度律.     

LEED structure determination of the Ni(1 1 1)(√3×√3)R30°-Sn surface

Quantitative low energy electron diffraction has been used to determine the structure of the Ni(1 1 1)(√3×√3)R30°-Sn surface phase. The results confirm that the surface layer comprises a substitutional alloy of composition Ni₂Sn as previously found by low energy ion scattering (LEIS), and also shows that there is no stacking fault at the substrate/alloy interface as has been found in (√3×√3)R30°-Sb surface alloys on Ag and Cu(1 1 1). The surface alloy layer is rumpled with the Sn atoms 0.45 ± 0.03 Å higher above the substrate than the surrounding Ni atoms. This rumpling amplitude is almost identical to that previously reported on the basis of the LEIS study. Comparison with similar results for Sn-induced surface alloy phases on Ni(1 0 0) and Ni(1 1 0) shows a clear trend to reduced rumpling with reduced surface atomic layer density, an effect which can be rationalised in terms of the different effects of valence electron charge smoothing at the surface.     

Sonoprocessing of wetting of SiC by liquid Al: A thermodynamic and kinetic study

The sonoprocessing of droplet spreading during the wetting process of molten aluminum droplets on SiC ceramic substrates at 700 °C is investigated in this paper. When wetting is assisted by a 20 kHz frequency ultrasonic field, the wettability of liquid metal gets enhanced, which has been determined by the variations in thermodynamic energy and wetting kinetics. Wetting kinetic characteristics are divided into two stages according to pinning and depinning states of substrate/droplet contact lines. The droplet is static when the contact line is pinning, while it is forced to move when the contact line is depinning. When analyzing the pinning stage, high-speed photography reveals the evidence of oxide films being rapidly crushed outside the aluminum droplet. In this work, atomic models of spherical Al core being wrapped by alumina shell are tentatively built, whose dioxide microstructures are being transformed from face-centered cubic into liquid at the atomic scale. At the same time, the wetting experiment reveals that the oxide films show changes in the period of sonoprocessing from 3rd to 5th second.During the ultrasonic spreading behavior in the late stage, there is a trend of evident expansion of the base contact area. The entire ultrasonic process lasts for no longer than 10 s. With the aid of ultrasonic sinusoidal waves, the wettability of metal Al gets a rapid improvement. Both molecular dynamic (MD) investigations and the experiments results reveal that the precursor film phenomenon is never found unless wetting is assisted by ultrasonic treatments. However, the precursor film appears near the triple line after using ultrasonics in the droplet wetting process, whose formation is driven by ultrasonic oscillations. Due to the precursor film, the ultrasonic wetting contact angle is lower than the non-ultrasonic contact angle. In addition, the time-variant effective ultrasonic energy has been quantitatively evaluated. The numerical expressions of thermodynamic variables are well verified by former ultrasonic spreading test results, which altogether provide an intrinsic explanation of the fast-decreasing contact angle of Al/SiC.     

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.     

Coverage-dependent collective diffusion of a dense Pb wetting layer on Si(111)

The dynamics of a dense Pb wetting layer on the Si(111) surface is studied in the framework of a generalized Frenkel-Kontorova model. Instead of the typical diffusion by random hopping processes, a liquidlike collective motion of the Pb atoms within the dense wetting layer is revealed to give rise to ultrafast kinetics of the wetting layer even at low temperatures. A kinetic Monte Carlo simulation including this collective spreading mechanism of the dense wetting layer quantitatively reproduces the experimental observations.     

Light-enhanced reactivation of passivated boron in hydrogen treated silicon

The thermal reactivation of the boron acceptor in hydrogenated silicon is strongly enhanced if the anneals are performed under illumination. The changes of the inactive boron concentration R versus annealing time t satisfy the equation dR/dt = -rR²/(N_A − R)², where N_A is the total boron concentration. Under strong illumination the annealing parameter r is proportional to the light intensity.