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

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

氦团簇及氦间隙原子在钨中的稳定性研究

首先采用分子动力学方法研究了在钨中预存氦-空位团簇（He_nV₂₂）后氦原子结合能与氦-空位比的关系。研究发现:当氦-空位比小于4.5时,氦原子结合能随氦-空位比呈线性减小趋势;当氦-空位比大于4.5时,氦原子的结合能随氦-空位比出现剧烈振荡的现象,这种现象是由于钨中预存氦-空位团簇随机挤出位错环使体系能量骤降所导致的。与此同时,氦-空位团簇周围出现了一些处于亚稳态的fcc结构和hcp结构的钨。为了研究氦团簇周围压强对钨基体相变的影响,本文利用第一性原理对钨的三种结构进行了高压相变计算,发现静水压力不能使钨的三种结构互相转变。另外,通过对bcc钨和fcc钨中四面体间隙氦原子和八面体间隙氦原子电荷密度差的计算,发现bcc钨中四面体间隙氦原子的稳定性高于八面体间隙氦原子的稳定性,而在fcc钨中四面体间隙氦原子的稳定性弱于八面体间隙氦原子的稳定性。     

银团簇结构与特性的分子动力学研究

用分子动力学方法模拟了银团簇的结构与力能学。计算模拟中使用了一种基于第一性原理的原子间互作用多体势函数。通过分子动力学模拟确定了银徽团簇（原子个数3－13）的稳态结构;模拟了原子个数为13－141的银FCC晶体结构理想球形团簇的力能学,发现球形银簇形成三维紧密结构;计算了平均结合能,给出了结合能随团簇原子数N的变化图,发现随N增大团簇结合能逐渐接近块材的数值。     

嵌入原子法计算金属钚中点缺陷的能量

钚因放射性衰变而出现老化效应.钚中点缺陷的性质和行为是理解钚老化效应的一个基础和前提.运用分子动力学模拟技术,计算了金属钚中点缺陷和点缺陷团簇的形成能和结合能.其中钚-钚、钚-氦和氦-氦相互作用势分别采用嵌入原子多体势、Morse对势和Lennard-Jones对势.计算结果表明,单个自间隙原子易以〈100〉哑铃状形态存在;间隙氦原子在理想晶格的八面体间隙位置相对较为稳定;氦原子与空位的结合能较大,在钚的自辐照过程中两者易于结合并形成氦-空位团簇;氦-空位团簇的形成能随氦原子数的增加而增大,当氦与空位的数     

金属钛中氦团簇融合的分子动力学模拟

运用分子动力学方法研究了金属钛中氦的扩散聚集行为.在300—800K的温度范围内,模拟了钛基底中氦团簇之间的融合过程.研究发现,温度的升高会加快氦团簇的融合.在300—800K,融合后的氦团簇在所模拟的时间尺度内三维结构保持不变.模拟结果还表明,常温下氦团簇之间的吸引力是导致氦团簇融合的重要因素. 关键词： 氦团簇 团簇融合 分子动力学模拟     

低温下Cu13团簇负载于Cu(001)表面上结构变化的分子动力学研究

应用分子动力学方法研究温度为10和50 K时具有二十面体结构的Cu₁₃团簇以不同接触条件与Cu（001）表面结合后的结构变化,原子间的相互作用势采用Johnson的嵌入原子方法模型.通过基于原子密度分布函数的分析表明,负载团簇与表面的结合能主要受团簇与载体相接触的最低层原子数及这些原子所具有的不同几何构型影响,同时更高层的原子呈现出不同的几何结构.温度为10 K时,负载团簇的初始位置对团簇几何结构和结合能影响较大. 关键词： 分子动力学 团簇 表面 计算机模拟     

分子动力学方法研究钛中预存缺陷对氦融合的影响

采用分子动力学方法模拟了在有预存氦空位复合物的体系中氦的迁移和融合,研究了氦融合的各向异性以及氦空位比和氦空位簇大小对氦融合的影响.研究发现:氦迁移的各向异性与体系中是否有氦空位簇无关;氦与氦空位簇的融合主要沿[001]方向;氦二聚物与氦空位簇的融合远快于单个氦原子与氦空位簇的融合.     

低温下Cu13团簇负载于Cu(001)表面上结构变化的分子动力学研究

应用分子动力学方法研究温度为10和50 K时具有二十面体结构的Cu13团簇以不同接触条件与Cu(001)表面结合后的结构变化,原子间的相互作用势采用Johnson的嵌入原子方法模型.通过基于原子密度分布函数的分析表明,负载团簇与表面的结合能主要受团簇与载体相接触的最低层原子数及这些原子所具有的不同几何构型影响,同时更高层的原子呈现出不同的几何结构.温度为10 K时,负载团簇的初始位置对团簇几何结构和结合能影响较大.     

低能Pt原子团簇沉积过程的分子动力学模拟

利用分子动力学模拟系统研究了低能Pt38,Pt141和Pt266原子团簇与Pt(001)表面的相互作用过程,详细分析了初始原子平均动能为0.1,1.0和10eV的原子团簇的沉积演化过程及其对基体表面形貌的影响.研究表明,初始原子平均动能是描述低能原子团簇的重要参量.当团簇的平均原子动能较低时,团簇对基体表层原子点阵损伤较小,基本属于沉积团簇;随着入射团簇的原子平均动能的增加,团簇对表层原子点阵结构的破坏能力增强,当团簇的原子平均动能增加到10eV时,团簇已经显现出注入特征.低能原子团簇对基体表面形貌的影响 关键词： 分子动力学模拟 低能原子团簇 载能沉积     

AuCu249合金团簇热稳定性的原子尺度计算研究

采用基于嵌入原子法的NVT正则系综分子动力学方法,在原子尺度上计算了包含249个原子的金属间化合物AuCu₂₄₉合金团簇由固态转变为熔体的结构演化过程. 根据对分布函数、原子密度分布函数和主要原子键对数目随温度的变化,发现在温度从低温升到高温的过程中,合金团簇内伴随着原子的连续位置交换,团簇呈现由外及里的分阶段结构转变. 同时还发现在团簇内原子堆积结构转变过程中,Au原子出现由团簇内层向外层运动的趋势,而Cu原子则有由外层向内层运动的趋势. 关键词： 合金团簇 分子动力学 计算机模拟 相变     

Bound states of an electron and a macroscopic cluster at a liquid helium surface

In a strong electric field, there are bound states of an electron at the surface of liquid helium, interacting with a large cluster of atoms in the bulk of liquid. This phenomenon is related to long-range interaction between the electron and the dipole moment of the cluster. The electron, holding the cluster under the liquid surface, is localized at this surface. One electron is capable of binding a cluster of up to 10⁶ atoms. The value of the binding energy may reach up to several kelvins.     

Ab Initio Investigation of Helium Bubble Formation from Trapped He+ Ions in Solid Hydrogen

We are involved in a program aimed at a first-principles examination of the possibility of helium bubble formation in solid hydrogen. Our procedure is based on an extension of the Hartree–Fock cluster approach that we have been using for study of individual He⁺ ions in solid hydrogen. The preliminary results of our cluster investigations of a particular configuration of two He⁺ ions in neighboring tetrahedral interstitial positions are presented here. It has been found that the hydrogen molecules common to the two He⁺ ions lead to a substantial attractive force between the ions that almost overcomes their inherent repulsion, leading to a binding energy for this configuration only slightly less than the binding energy for two well separated trapped He⁺ ions. The analysis of our results have suggested future investigations of different pairs of trapping sites which appear to be more favorable for providing higher binding energies for the He⁺ pairs. This revised version was published online in August 2006 with corrections to the Cover Date.     

Molecular Dynamics Simulations of Helium Behaviour in Titanium Crystals

Molecular dynamics simulations are performed to investigate the behaviour of helium atoms in titanium at a temperature of 30OK. The nucleation and growth of helium bubble has been simulated up to 50 helium atoms. The approach to simulate the bubble growth is to add helium atoms one by one to the bubble and let the system evolve. The titanium cohesion is based on the tight binding scheme derived from the embedded atom method, and the helium-titanium interaction is characterized by fitted potential in the form of a Lennard-Jones function. The pressure in small helium bubbles is approximately calculated. The simulation results show that the pressure will decrease with the increasing bubble size, while increase with the increasing helium atoms. An analytic function about the quantitative relationship of the pressure with the bubble size and number of helium atoms is also fitted.     

The positron-helium bound state

A variational calculation with a trial wave function including up to 17 parameters has been used to investigate the existence of a positron-helium bound state. The resulting binding energy was found to be greater than the helium binding energy, and thus it is believed that the positron does not get bound to a helium atom. The calculations provide also values of ground state helium energy and the corresponding wave functions.     

Hartree-Fock interaction between a helium atom and lithium metal: A test for the effective medium theory

The interaction of a helium atom with 6 and 10 atom clusters of lithium has been calculated using the unrestricted Hartree-Fock method, Hartree-Fock method with correlation corrections, and the effective medium theory. Inside the cluster the helium embedding energy is found to be proportional to the electron density of the cluster. The proportionality constant obtained by the Hartree-Fock method is in fair agreement with that calculated in the homogeneous electron gas using the local density approximation. Outside the cluster, in the region compatible with the helium diffraction experiments, the self-consistent calculations give much larger repulsion than the effective medium theory.     

Simulation of nanostructure evolution under helium implantation in Fe–(2.5–12.5)at% Cr alloys at a temperature of 343 K

Characterization of helium-implanted Fe–(2.5–12.5)at% Cr alloys with the flux of 7?×?10^–6?dpa/s at a temperature of 343?K has been performed by means of cluster dynamics simulations. We have suggested a model for simulating an Fe–C–Cr system under helium implantation based on a selection of the latest data from atomistic studies and available experiments. Kinetics of carbon-vacancy and helium-vacancy complexes has been studied. Only one parameter is used to guarantee the best reproduction possible of experimental positron annihilation lifetime spectroscopy data for Fe–Cr alloys on dependences of vacancy cluster size on chromium content and irradiation dose via fitting. This is an effective binding energy of self-interstitial atoms to dislocation loops decorated by chromium atoms. It has a “snaky” dependence of chromium content with a minimum of about 9%Cr.