Cu刃型扩展位错附近局部应变场的原子模拟研究

通过结合virial应变分析技术的准连续介质多尺度模拟方法研究了金属Cu刃型扩展位错的局部应变场.结果表明在距离位错核心几十纳米的区域内晶体处于小变形状态,virial应变计算结果与弹性理论预测结果符合得相当好,当距离位错核心仅几纳米时,晶格畸变加剧,virial应变极大值出现在扩展位错两端的Shockley分位错芯部.进一步分析表明Shockley分位错芯部严重畸变区大致呈长轴7b1、短轴3b1的椭圆形,其中b1为分位错柏氏矢量的长度.     

纳米接触行为的三维多尺度数值模拟

应用分子动力学与有限元耦合的桥域多尺度算法,模拟三维刚性球形压头与光滑基体表面的纳米尺度接触行为,并与全原子分子模拟结果比较.考察在一定载荷下的系统弛豫行为、两种模型桥接区位移和应力的连续性、法向力和接触面积随压头位移变化等,结果表明:一定外载荷下,桥域多尺度算法能较快达到平衡状态,且压头的振荡幅度更小,系统初始温度为0 K时该算法的相对误差最小.在准静态加载过程中,该算法能够将原子区的位移、应力等连续的过渡到连续介质区,具有较好的耦合效果;法向力-压头位移和接触半径-压头位移曲线几乎与分子模拟结果重合,表明算法具有较高的计算精度.     

界面结构对Cu/Ni多层膜纳米压痕特性影响的分子动力学模拟

金属多层膜调制周期下降到纳米级时,其力学性质会发生显著改变. Cu-Ni晶格失配度约为2.7%,可以形成共格界面和半共格界面,实验中实现沿[111]方向生长的调制周期为几纳米且具有异孪晶界面结构的Cu/Ni多层膜,其力学性质发生显著改变.本文采用分子动力学方法对共格界面、共格孪晶界面、半共格界面、半共格孪晶界面等四种不同界面结构的Cu/Ni多层膜进行纳米压痕模拟,研究压痕过程中不同界面结构类型的形变演化规律以及位错与界面的相互作用,获取Cu/Ni多层膜不同界面结构对其力学性能的影响特征.计算结果表明,不同界面结构的样品在不同压痕深度时表现出的强化或软化作用机理不同,软化机制主要是由于形成了平行于界面的分位错以及孪晶界面的迁移,强化机制主要是由于界面对位错的限定作用以及失配位错网状结构与孪晶界面迁移时所形成的弓形位错之间的相互作用.     

纳米流体多相流动的多尺度模拟方法

针对纳米流体多相流动的微观特征,提出一种基于格子Boltzmann的多尺度耦合方法,在速度和悬浮纳米粒子分布变化比较剧烈的区域采用细网格多相模型,在其它区域视纳米流体为均匀混合的单相流体,使用粗网格单相模型.为保证不同尺度区域之间的物理信息(参数)的准确传递,运用质量和动量守恒原理,建立跨区域的边界耦合模型.几个算例表明,该方法既可以反映纳米流体流动的微观特征,又能提高计算效率,与单纯使用多相模型相比,节省大量时间.     

单空位体系多尺度模拟

基于多尺度协同算法，实行了统一动力学模拟方案，研究了单空位体系的位移分布与位移振荡行为、能量随时间起伏以及差分电荷密度与外部环境的影响.计算结果发现位移分布及其相关行为反映了单空位体系的量子力学效应；同时，体系总能、差分电荷密度计算则给出了外部环境的影响，反映了经典效应与量子效应的耦合. 关键词： 多尺度杂化算法 统一动力学 电子态 原子位移     

动力学晶格蒙特卡洛方法模拟Cu薄膜生长

利用动力学晶格蒙特卡洛方法模拟了Cu薄膜在Cu(100)面上的三维生长过程。模型中考虑了四个动力学过程:原子沉积、增原子迁移、双原子迁移和台阶边缘原子迁移,各动力学过程发生的概率由多体势函数确定。讨论了基底温度、沉积速率及原子覆盖率对Cu原子迁移、成核和表面岛生长等微观生长机制的影响;获得了Cu薄膜的表面形貌图并计算了表面粗糙度。模拟结果表明,随基底温度升高或沉积速率下降,岛的平均尺寸增大,数目减少,形状更加规则。低温时,Cu薄膜表现为分形的离散生长,高温时,Cu原子迁移能力增强形成密集的岛。Cu薄膜表面粗糙度随着基底温度的升高而迅速减小;当基底温度低于某一临界温度时,表面粗糙度随原子覆盖率或沉积速率的增大而增大;当基底温度超过临界温度时,表面粗糙度随原子覆盖率或沉积速率的变化很小,基本趋于稳定。     

Au/Au(100)外延薄膜生长的计算机模拟及其微观机制研究

复旦大学现代物理研究所基于加速器的原子及原子核物理开放实验室,上海200433 关键词：     

超薄膜多中心生长过程的计算机模拟

利用计算机模拟了不同的允许扩散步数下超薄膜的多中心分形生长和团状生长现象,研究了成核及长大的动力学过程．分形生长时分形维数随团簇大小的增大而增加．分形生长和团状生长时成核率随扩散步数的增大而减小,随时间的增大而急速下降．团簇长大过程可用团簇大小S和生长时间t-t₀的幂函数(t-t₀)^κ描述．由于团簇间的分流作用,生长指数κ比经典理论值1略小,并且存在着非线性现象,即长得较大的团簇的生长指数Κ也较大. 关键词：     

论FeCr合金辐照损伤的多尺度模拟

铁素体/马氏体钢具有优异的抗辐照肿胀和抗腐蚀能力,是聚变堆、ADS和四代堆结构材料重要候选者之一.多尺度模拟在核武器自辐照损伤、压力容器辐照脆化、医学材料辐照损伤等方面有着广泛的应用.国际上围绕F/M FeCr合金辐照损伤开展了大量的多尺度模拟工作.论文概述了F/M FeCr合金辐照损伤多尺度模拟研究现状及进展,并针对目前存在的某些问题提出了解决建议. 关键词： 铁素体/马氏体FeCr合金 辐照损伤 多尺度模拟     

Atomistic simulations of fluid structure and solvation forces in atomic force microscopy

We describe results of atomistic molecular dynamics simulations modelling an atomic force microscope (AFM) tip immersed in a fluid. Both the tip and the surface are modelled by rigid arrays of atoms. The tip is pyramidal and the surface is the (100) face of a fcc crystal. The focus is on the solvation forces acting on the tip and on the surface and their relation to the structural and dynamic properties of the fluid. Fluid particles in the neighborhood of the tip-surface junction are found to be highly ordered compared to the bulk, as shown by localized variations in the average fluid density. The atomistic nature of the model gives rise to several effects related to the discrete sizes of the fluid, tip, and surface particles which are not observed in continuum-based theories. A number of simulated force-distance curves are presented, along with an analysis of the effect of changing fluid particle size, tip (lateral) position, tip shape, and the lyocompatability of the tip and surface materials. The atomic-scale distribution of fluid-surface forces is examined for various positions of the tip, and the extent to which the fluid can act as a “cushion” by increasing the effective area of the tip-surface interaction is studied. The effect of a fluid on AFM imaging is investigated by generating “fluid images”, which are shown to be comparable in magnitude to the direct tip-surface interaction in the noncontact mode. We compare images generated by defective and defect-free surfaces, and analyse the fluid-tip forces acting in a lateral direction. An image formed from fluid forces acting in the direction of the surface normal does not show the presence of a vacancy, but an image formed from lateral fluid forces does.     

Dynamical behavior of the phase transition of strained from atomistic simulations

Using an atomistic shell model we study the temperature dependence of the ferroelectric properties of BaTiO₃ under biaxial compressive strain applicable to growth on perovskites substrate. Molecular dynamics simulations show a “r→c→p” sequence of phase transitions when temperature is increased, and the absence of an “ac phase”. The first-order paraelectric-to-ferroelectric phase transition presents in bulk changes to a second-order one as a consequence of the in-plane constraint imposed by the mechanical boundary conditions. From the tetragonal ferroelectric c phase, the transition takes place in a finite range of temperature where the lattice parameter normal to the plane keeps approximately constant until T_c is reached. Analysis of the local polarization behavior reveals an order–disorder dynamics as the dominant mechanism of the transition.     

Atomistic tight-binding computations of the structural and optical properties of CdTe/CdX (X=S and Se)/ZnS core/shell/shell nanocrystals

A study of CdTe/CdX (X=S and Se)/ZnS core/shell/shell nanocrystals is carried out using atomistic tight-binding theory and the configuration interaction method to provide information for applications in bioimaging, biolabeling, display devices and near-infrared electronic instruments. The calculations yield the dependences of the internal and external passivated shells on the natural behaviours of CdTe/CdX (X=S and Se)/ZnS core/shell/shell nanocrystals. The reduction of the optical band gaps is observed with increasing numbers of monolayers in the external ZnS shell due to quantum confinement. Interestingly, the optical band gaps of CdTe/CdS/ZnS core/shell/shell nanocrystals are greater than those of CdTe/CdSe/ZnS core/shell/shell nanocrystals. In the presence of an external ZnS-coated shell, electron–hole wave function overlaps, oscillation strengths, ground-state exchange energies and Stokes shift are improved, whereas ground-state coulomb energies and fine-structure splitting are reduced. The oscillation strengths, Stokes shift and fine-structure splitting are reduced with the increase in external ZnS shell thickness. The oscillation strengths, Stokes shift and fine-structure splitting of CdTe/CdS/ZnS core/shell/shell nanocrystals are larger than those of CdTe/CdSe/ZnS core/shell/shell nanocrystals. Reduction of the atomistic electron–hole interactions is observed with increasing external ZnS shell size. The strong electron–hole interactions are more probed in CdTe/CdS/ZnS core/shell/shell nanocrystals than in CdTe/CdSe/ZnS core/shell/shell nanocrystals.     

Atomistic modeling of Zn deposition on the low index faces of Cu

The formation process of Zn/Cu surface alloys is investigated using the Bozzolo-Ferrante-Smith (BFS) method for alloys. The effects of the crystallographic orientation on the deposition process, formation of surface alloys as a function of temperature and coverage, Zn surface migration, and interdiffusion in the Cu substrate, are modeled and discussed with atom-by-atom energy analyses and large scale simulations.     

Microstructures and strengthening mechanisms of Cu/Ni/W nanolayered composites

Cu/Ni/W nanolayered composites with individual layer thickness ranging from 5?nm to 300?nm were prepared by a magnetron sputtering system. Microstructures and strength of the nanolayered composites were investigated by using the nanoindentation method combined with theoretical analysis. Microstructure characterization revealed that the Cu/Ni/W composite consists of a typical Cu/Ni coherent interface and Cu/W and Ni/W incoherent interfaces. Cu/Ni/W composites have an ultrahigh strength and a large strengthening ability compared with bi-constituent Cu–X (X?=?Ni, W, Au, Ag, Cr, Nb, etc.) nanolayered composites. Summarizing the present results and those reported in the literature, we systematically analyze the origin of the ultrahigh strength and its length scale dependence by taking into account the constituent layer properties, layer scales and heterogeneous layer/layer interface characteristics, including lattice and modulus mismatch as well as interface structure.     

Magnetic properties of Ni/Au core/shell studied by Monte Carlo simulations

The magnetic properties of ferromagnetic Ni/Au core/shell have been studied using Monte Carlo simulations within the Ising model framework. The considered Hamiltonian includes the exchange interactions between Ni–Ni, Au–Au and Ni–Au and the external magnetic field. The thermal total magnetizations and total magnetic susceptibilities of core/shell Ni/Au are computed. The critical temperature is deduced. The exchange interaction between Ni and Au atoms is obtained. In addition, the total magnetizations versus the external magnetic field and crystal filed for different temperature are also established.     

Computer simulations of the Stranski-Krastanov growth of heteroepitaxial films with elastic anisotropy

A three-dimensional finite element method is developed to simulate the surface morphological evolution during the Stranski-Krastanov heteroepitaxial growth. In the formulation, the surface evolves through surface diffusion driven by the gradient of the surface chemical potential, which includes the elastic strain energy, elastic anisotropy and surface energy. Surface condensation rate is assumed to depend on the difference between the surface chemical potential and the chemical potential of the vapor phase. Our simulations reveal that the self-assembly of quantum dots are strongly dependent on the variation of growth rate and elastic anisotropy strength. With appropriate choice of growth rate and elastic anisotropy strength, a relatively more uniform and regular quantum dot array can be obtained.     

Reversible H-induced switching of the magnetic easy axis in Ni/Cu(001) thin films

A reversible switching of the easy axis of magnetization for Ni on Cu(001) from in plane to out of plane is found by changing the partial pressure of hydrogen in the gas phase around the sample, allowing even for oscillations of the magnetization direction. A quantitative low-energy electron diffraction study of the diffracted intensity versus electron energy [I(E)] shows that the hydrogen-induced spin reorientation transition is accompanied by changes of the tetragonal distortion of the topmost Ni layer upon hydrogen adsorption. Surprisingly, the orientation switch to perpendicular to the surface comes with a relaxation, i.e., reduction of the film's tetragonal distortion rather than its amplification.     

SHG investigations of the magnetization of thin Ni and Co films on Cu(001)

Received: 22 November 1998