金属钛中氦团簇生长行为的分子动力学研究

运用分子动力学方法模拟了常温下金属钛中氦团簇的生长过程.从能量的角度考察了氦团簇的生长机理.研究发现,随着氦团簇的生长,在氦团簇周围逐渐形成位错环缺陷,氦团簇与氦原子的结合能有逐渐下降的趋势,当氦团簇生长到一定尺寸时会通过发射周围缺陷以使得结合能上升,从而增强了进一步吸收氦原子的能力.模拟还发现,随着氦团簇的继续生长,氦团簇的形状由原来的不规则结构逐渐变成了较为规则的棱柱形结构,在随后的生长过程中其生长仅在（001）平面进行,沿［001］轴的厚度几乎不变. 关键词： 氦团簇 缺陷发射 分子动力学模拟     

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

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

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

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

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

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

分子动力学方法研究金属Ti中He小团簇的迁移

采用分子动力学方法模拟了不同温度下He原子及He团簇在金属Ti中的迁移特性,并计算了扩散前系数和激活能. 研究发现这种扩散的各向异性非常显著,具体表现在He原子及He团簇在不同方向上扩散系数的前系数完全不同,但它们的激活能却相同. 研究表明:在预测金属中He的扩散行为时,必须采用动态模拟方法才能得到准确的前系数,仅仅考虑势垒的静态模拟方法是不行的. 另外,还发现一个不同于直觉的现象,即较低温度下He二聚物的扩散系数比单个He原子的扩散系数大;此外,在所模拟的温度范围内Arrhenius方程能够很好地描述它们的扩散. 这说明动力学模拟对预测金属中He的扩散行为具有重要的意义. 关键词： 分子动力学 扩散系数 各向异性 氦     

压力对金属铌中氦原子凝聚的影响

 采用零温条件下的赝势-平面波方法和有限温度下的Car-Parrinello分子动力学方法,模拟了不同压力环境下氦原子在金属铌中的行为特征,研究了宿主缺陷和氦泡的形成机制。结果表明,闭电子壳层的氦原子在金属铌中具有刚球模型特征,其占据区域为金属自由电子的禁区,从而破坏铌原子之间的金属性键合。在常温条件下,局域高浓度的氦原子优先凝聚于近邻宿主空位缺陷处,从而形成氦泡;完整晶格中高浓度的氦将促使铌原子易位,形成间隙-空位模式的宿主缺陷,氦原子聚集于空位区域。完整宿主在压力（40 GPa）的作用下,晶格参数减小,铌原子之间的相互作用增强,尽管氦原子的存在削弱了铌原子之间的相互作用,位于格点上的铌原子仍难以借助热振动偏离格点形成空位,因而未能形成间隙-空位对和氦泡。     

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

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

Ge/Pb/Si(111)生长中Ge原子沿团簇边缘扩散的三维Monte-Carlo模拟

用动态Monte-Carlo方法对Ge在单层表面活性剂Pb覆盖的Si(111)表面上沿团簇边缘扩散进行三维模拟.重点讨论Ge原子是否沿团簇边缘扩散,沿边缘扩散时的最大扩散步数及最近邻原子数对三维生长的影响,并计算薄膜表面粗糙度研究三维生长模式.模拟表明Ge沿团簇边缘扩散的行为对薄膜生长模式的影响很大,同时讨论了ES势对三维生长模式的影响.     

Diffusion Dynamics of Cux Cluster on Cu(111) Surface

利用分子动力学与修正分析型嵌入原子模型研究了从200K至800 K范围内,Cux (1·x·8)小团簇在(111)面的自扩散动力学行为,计算了Cu小团簇的扩散系数与激活能. 计算结果表明,密排的Cu7团簇的激活能最高,其扩散前因子比单个吸附原子的高3个数量级,这与其他类似金属的实验结果一致. 另外,还定量讨论了薄膜生长与团簇扩散的关系.     

单轴应力场下钨中氦扩散行为的分子动力学模拟研究

本文采用分子动力学方法研究了沿<100>及<111>晶向的单轴应变对钨中单个氦原子扩散的影响。结果表明,应变会使得金属钨材料发生相变,且引起相变的临界应变随温度升高而减小。相变起始的应变在达到抗拉强度的应变附近。计算结果表明,拉应变使得单个氦原子在钨中的扩散系数发生骤降,在不同应变下扩散系数变化平缓。沿<100>晶向氦扩散系数随应变的增大而线性减小,而<111>晶向则出现了震荡变化趋势。研究结果表明,沿<100>晶向应变达到+0.15%时阿纽列斯方程不再适用,而沿<111>晶向应变大+5%阿纽列斯方程仍然适用;沿<111>晶向随应变增加氦扩散激活能减小,说明应变使得单个氦原子在钨中迁移性增强。     

Ab initio potential energy curve for the helium atom pair and thermophysical properties of the dilute helium gas. II. Thermophysical standard values for low-density helium

A helium–helium interatomic potential energy curve determined from quantum-mechanical ab initio calculations and described with an analytical representation considering relativistic retardation effects (R. Hellmann, E. Bich, and E. Vogel, Molec. Phys. (in press)) was used in the framework of the quantum-statistical mechanics and of the corresponding kinetic theory to calculate the most important thermophysical properties of helium governed by two-body and three-body interactions. The second pressure virial coefficient as well as the viscosity and thermal conductivity coefficients, the last two in the so-called limit of zero density, were calculated for ³He and ⁴He from 1 to 10 000 K and the third pressure virial coefficient for ⁴He from 20 to 10 000 K. The transport property values can be applied as standard values for the complete temperature range of the calculations characterized by an uncertainty of ±0.02% for temperatures above 15 K. This uncertainty is superior to the best experimental measurements at ambient temperature.     

四体相互作用对固氦压缩特性的贡献

运用从头计算自洽场方法和原子团簇理论计算了高压下固氦原子间的四体势分量.计算结果表明四体分量对结合能的贡献为正；随着压缩度增大四体势的贡献比例变大.采用两体、三体、四体势和Aziz吸引势计算固氦零温状态方程并与实验测量相比较，结果表明两体势对固氦的压缩性贡献了过多的正效应，加入三体分量的修正，仅在低于10GPa时理论值与实验值相符很好，但考虑了四体势修正后能将理论值与实验值符合程度提高到 27GPa. 关键词： 状态方程 固氦 四体势 从头计算     

多体相互作用对高压固氦状态方程的影响

采用量子化学自洽场方法及原子团簇理论对固氦晶格原子势能进行多体展开，定量确定了最近邻原子间距R在0.26—0.175nm之间，短程的两体到五体相互作用对晶体结合能及压缩特性的贡献，并对不同间距下多体展开式的截断位置进行了讨论. 结果表明：当最近邻原子间距R在0.26—0.175 nm之间取值，在具有hcp相结构的晶格中，氦原子势能的多体展开式是一收敛的交错型级数，该级数中两体项、四体项为正值，而三体项、五体项为负值. 当R在0.26—0.23 nm之间取值时, 仅保留级数中两体和三体项就能很好地描述晶格原子势能；当R在0.23—0.20 nm之间取值，多体展式需要展开到四体相互作用项；当R在0.20—0.175 nm之间取值必须再加入五体项贡献. 考虑到五体相互作用后，理论计算已能圆满地解释目前固氦等温压缩数据，其最大压力达 60GPa. 关键词： 固氦 状态方程 多体相互作用     

Simulation of Helium Behaviour in Titanium Crystals Using Molecular Dynamics

The behaviour of helium in Ti crystals at 300 K has been investigated by means of the molecular dynamics. The study is focused on the influences of He-Ti interaction on the aggregation of helium atoms in the substrate. When a Born-Mayer potential is used to describe the He-Ti interaction, the He atoms are unable to cluster with each other due to the weak bridge barrier that cannot trap the helium atoms, Whereas using a He-Ti potential that is constructed by fitting the ab initio pairwise He-Ti potential, the clustering of He atoms can be observed. The results indicate that suitable He-Ti potential plays an important role in the formation of He clusters in metals. Moreover, it is noted that the shape of the formed He cluster is irregular, and the produced defect prefers to congregating on one side of the He cluster rather than spreading symmetrically around it.     

Equation of state of fluid helium at high temperatures and densities

Hugoniot curves and shock temperatures of gas helium with initial temperature 293 K and three initial pressures 0.6, 1.2, and 5.0 MPa were measured up to 15000 K using a two-stage light-gas gun and transient radiation pyrometer. It was found that the calculated Hugoniot EOS of gas helium at the same initial pressure using Saha equation with Debye-Hückel correction was in good agreement with the experimental data. The curve of the calculated shock wave velocity with the particle velocity of gas helium which is shocked from the initial pressure 5 MPa and temperature 293 K, i.e., theD ≈u relation,D=C ₀+λu (u<10 km/s, λ=1.32) in a low pressure region, is approximately parallel with the fittedD ≈u (λ=1.36) of liquid helium from the experimental data of Nellis et al. Our calculations show that the Hugoniot parameter λ is independent of the initial density p{in0}. TheD≈u curves of gas helium will transfer to another one and approach a limiting value of compression when their temperature elevates to about 18000 K and the ionization degree of the shocked gas helium reaches 10⁻³.     

Equation of state of fluid helium at high temperatures and densities

Helium, hydrogen, and their isotopes are the simplest monoatomic and diatomic molecules. It is relatively easy to describe their properties using the basic principles of quantum mechanics. In condensed matter physics, hydrogen and helium serve as the models for the condensed matter properties at extreme conditions so that both experi- mental and theoretical physicists pay much attention to the study of their properties[1], especially the insulator-metal transition of hydrogen[2]. The aim to st…     

Field ionization of free helium atoms: Correlation between the kinetic energy of ionized atoms and probability of their field ionization

In this paper the correlation between the kinetic energy of helium atoms and the probability of field ionization is investigated by exploiting the narrow velocity distribution of supersonic molecular beams. Field ionization measurements were carried out on supersonic helium beams at 298 K and 95 K corresponding to energies of about 65 meV and 20 meV, respectively, for the individual atoms. The field ionization was performed with a tungsten tip, radius of curvature 12 nm, kept at room temperature. The ionization probability was found to increase by about a factor 10 when the beam was cooled from 298 K to 95 K. The results presented in this paper are of importance for improving the understanding of field ionization and for the development of a new detector for helium and other molecular beams.     

The thermodynamic properties of normal liquid helium 3

The thermodynamic properties of normal liquid helium 3 are calculated by using the lowest order constrained variational (LOCV) method. The Landau Fermi liquid model and Fermi-Dirac distribution function are considered as our statistical model for the uncorrelated quantum fluid picture and the Lennard-Jones and Aziz potentials are used in our truncated cluster expansion (LOCV) to calculate the correlated energy. The single particle energy is treated variationally through an effective mass. The free energy, pressure, entropy, chemical potential and liquid phase diagram as well as the helium 3 specific heat are evaluated, discussed and compared with the corresponding available experimental data. It is found that the critical temperature for the existence of the pure gas phase is about 4.90 K (4.45 K), which is higher than the experimental prediction of 3.3 K, and the helium 3 flashing temperature is around 0.61 K (0.50 K) for the Lennard-Jones (Aziz) potential.