高密度氦状态方程的分子动力学模拟研究

 用分子动力学模拟和与固氦实验等温线比较的方法,确定了高密度氦指数6势函数中劲度系数α的优选值为12.7。然后我们用这个新的α值计算了fcc固氦的状态方程及其径向分布函数,发现当把ρ限定为1.60 g/cm³时,其径向分布函数的第二个峰将在2 000 K到3 040 K区间消失,表明此时固氦的长程有序性降低,或者说发生了固-液相变过程。     

1.5μm处CO2与CO高温线强的实验分析与理论计算

本文在采用乘积近似方法计算二氧化碳、一氧化碳分子总的配分函数(其中分子的振动配分函数采用谐振子近似, 转动配分函数采用非刚性转子模型, 并考虑了离心扭曲修正)的基础上, 利用所得配分函数和振动跃迁矩平方的实验值以及Herman-Wallis系数, 计算了1.5 μm 附近二氧化碳30012–00001跃迁带和一氧化碳3–0跃迁带在300–6000 K温度范围内部分温度下的吸收线强; 为验证计算方法和结果的准确性, 在基于可调谐二极管激光吸收光谱技术搭建的高温测量系统中, 对300–800 K温度范围内部分谱线线强进行了测量, 并把计算结果、测量结果及HITRAN数据库中对应数据进行了对比, 发现相对偏差小于3%, 证明了本方法的有效性, 同时计算及测量所得高温线强数据可对HITRAN数据库进行有效的校正和补充.     

派热克斯玻璃的固态相变区物态方程

派热克斯玻璃的Hugoniot实验数据显示了在20GPa以下存在相变,我们认为这是一种固态-固态的一级相变,根据Gibbs相平衡条件并应用了相变前的低压物态方程,计算了这种固态-固态相变的熵增。计算表明熵增为负值,即相变为放热过程。在此基础上给出了相变区的物态方程,其中冷压为平台,内能计入了相变熵的影响。由相变区的物态方程得到的理论的与实验的Hugoniot曲线基本一致。     

高温高密度氩物态方程研究

用Exp-6有效两体势模型和液体变分微扰理论计算了液Ar冲击压缩曲线,在35 GPa以下的压力范围内计算的冲击压缩曲线与Thiel 及Nellis等人的实验数据进行了比较.液体分子间Exp-6有效两体势参数根据实验数据计算优选出最佳值,计算结果表明得到的优选参数较为准确地描述了液Ar分子间相互作用.由理论结合实验数据可以认为液氩体系在冲击压力36 GPa以下无相变现象.对较高冲击压力下理论计算的冲击曲线和实验结果之间的偏差作了分析,结合不透明度实验的结果,得到当压力超过35 GPa,温度在12 000 K     

温度对金刚石涂层膜基界面力学性能的影响

利用分子动力学方法建立了硬质合金基底金刚石涂层膜基界面模型, 并采用Morse势函数和Tersoff势函数相互耦合的方法来表征模型内原子间的相互作用关系, 在此基础上对不同温度(0–800 K)条件下硬质合金基底金刚石涂层膜基界面的力学性能进行了分子动力学仿真计算. 结果表明: 当温度由0 K上升到800 K的过程中, 金刚石涂层膜基界面拉伸强度呈下降趋势, 并且在0–300 K范围内下降趋势明显, 在300–800 K范围内下降趋势缓和; 体系能量随温度的变化具有相同的下降趋势.     

氦泡生长实验数据的反演模拟计算方法

介绍了基于“迁移-合并”机理的氦泡生长理论模型与蒙特卡罗模拟方法，还研制了以上述算法为正算方式，对实验数据进行反演的迭代方法.最后应用这些方法对钯中氦泡生长实验数据进行了反演模拟计算，得到了氦泡半径与迁移速率的关系，并对计算结果作了分析.     

相变材料热物理性质的分子动力学模拟

为从微观尺度探寻相变材料的热物性变化机理, 本文采用分子动力学的方法, 构建了由正二十二烷组成的无定形结构的相变材料体系, 采用周期性边界条件以及COMPASS力场对相变材料的比热以及导热系数进行了模拟, 并对纯正二十二烷进行了DSC测试. 结果表明, 模拟所得的相变材料热容与文献实验值的偏差是6.5%, 熔点与DSC实验值的偏差是0.98%. 当温度为288–318 K时, 相变材料的导热系数在0.1–0.4 W·m^-1·K^-1 范围内波动, 且随着压力增大略呈下降趋势. 关键词： 扩散系数 比热 导热系数 分子动力学     

α-Fe中氦泡阻碍位错移动硬化的分子动力学研究

基于周期性位错阵列模型,利用分子动力学方法,研究了α-Fe中1/2a0111{110}刃位错与3nm氦泡的相互作用导致硬化的过程,重点研究了不同温度下(300 K和623 K),氦空位比例对位错与3 nm氦泡作用机制影响的差异性.研究结果表明,当氦空位比例在0～1之间变化时,氦空位比例的变化对临界剪切应力的影响微弱,当氦空位比例在1～1.5之间变化时,氦泡内氦原子的数量的增加,临界剪切应力随之降低.氦空位比例升高至1.75时,高温623 K与室温相比,临界剪切应力反常升高,原因是高温引起的位错与氦泡的排斥机制:氦泡严重超压,与位错接触瞬间,沿位错拉伸侧踢出自间隙原子团簇,被位错吸收后产生割阶;由于位错割阶与超压氦泡同为压应力场,割阶被氦泡强烈排斥反弹.     

钨中氢同位素热脱附实验的速率理论模拟研究

本文采用基于速率理论的模拟方法研究钨材料中氢同位素氘的热脱附谱. 热脱附数据来源于520 K下受等离子体辐照的多晶钨, 入射离子能量为40 eV, 剂量为1× 10²⁶ D/m². 通过调节速率理论中的俘获能、俘获率等参数, 最终获得与实验相符合的热脱附拟合谱. 拟合结果表明, 钨中俘获的氘存在于三种俘获态, 俘获能分别为1.14 eV, 1.40 eV和1.70 eV, 相应脱附温度峰值为500 K, 600 K和730 K. 这三个俘获能分别应对应于第一原理计算得到的空位俘获第3–5个氢原子的俘获能(含零点振动能修正)、空位俘获第1–2个氢原子的俘获能, 空位团簇对氢原子的俘获能. 模拟结果表明, 在本辐照实验条件下, 钨中空位及空位团簇是氘在钨中的主要俘获态.     

金属表面下氦泡融合与释放研究

采用分子动力学方法研究了钛金属表面下不同深度处氦泡的行为,分析了氦泡融合与释放的竞争,对比了不同深度处氦泡的释放对金属的影响.结果表明:在接近金属表面处,氦泡很难通过融合无限长大,当达到临界尺寸后,氦泡将会释放而不再与邻近的氦泡发生融合;植入深度对氦泡的融合有一定的影响,深度越大,越有利于形成具有较高氦密度的大氦泡;较深处氦泡的释放会在金属表面形成较大的突起和表面针孔.实验中观察到的不同尺寸的表面孔,其部分原因来自于金属表面下不同深度处氦泡的释放.     

Equations of State and the Phase Diagram of Dielectric and Metallic Helium-4

A new form of the semiempirical equation of state for the liquid phase of helium-4 is proposed that is based on the assumption that the structure of this phase consists of a mixture of dielectric and metallic components. It is postulated that solid dielectric helium with density higher than 5.3 g/cm³ becomes a metal. The values of the parameters of the equations of state for both solid phases and the liquid phase of helium are calculated. The unknown values of the initial data for helium are taken by analogy with the parameters for deuterium. The phase diagram, shock adiabat, isentropes, isotherms, and electrical conductivity in these processes are calculated with the use of the equations of state of solid and liquid phases of helium-4. The results of calculation are compared with experimental data in the range of pressures of up to 35 GPa for an isotherm, up to 150 GPa for a shock adiabat, up to 42 GPa for the melting curve, and up to 2000 GPa for isentropes, and showed quite satisfactory agreement. Numerical extrapolation of the melting curve is performed to a range of ultrahigh pressures of up to 8000 GPa.     

Radiation defects induced by helium implantation in gold-based alloys investigated by positron annihilation spectroscopy

The formation of gas bubbles in metallic materials may result in drastic degradation of in-service properties. In order to investigate this effect in high density and medium-low melting temperature (T _M) alloys, positron annihilation spectroscopy measurements were performed on helium-implanted gold–silver solid solutions after isochronal annealing treatments. Three recovery stages are observed, attributed to the migration and elimination of defects not stabilized by helium atoms, helium bubble nucleation and bubble growth. Similarities with other metals are found for the recovery stages involving bubble nucleation and growth processes. Lifetime measurements indicate that He implantation leads to the formation of small and over-pressurized bubbles that generate internal stresses in the material. A comprehensive picture is drawn for possible mechanisms of helium bubble evolution. Two values of activation energy (0.26 and 0.53 eV) are determined below and above 0.7T _M, respectively, from the variation of the helium bubble radius during the bubble growth stage. The migration and coalescence mechanism, which accounts for these very low activation energies, controls the helium bubble growth.     

Annealing Behaviour of Helium Bubbles in Titanium Films by Thermal Desorption Spectroscopy and Positron Beam Analysis

Helium-containing Ti Glms are prepared using magnetron sputtering in the helium-argon atmosphere. Isochronal annealing at different temperatures for an hour is employed to reveal the behaviour of helium bubble growth. Ion beam analysis is used to measure the retained helium content. Helium can release largely when annealing above 970K. A thermal helium desorption spectroscopy system is constructed for assessment of the evolution of helium bubbles in the annealed samples by linear heating （OAK/s） from room temperature to 1500K. Also, Doppler broadening measurements of positron annihilation radiation spectrum are performed by using changeable energy positron beam. Bubble coarsening evolves gradually below 680K, migration and coalescence of small bubbles dominates in the range of 68-970 K, and the Ostwald ripening mechanism enlarges the bubbles with a massive release above 970K.     

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

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

Using the Ab Initio Molecular Dynamics Method for Simulating the Peculiarities in the Temperature Dependence of Liquid Bismuth Properties

The specific features pertinent to the temperature dependence of the electronic and atomic properties of liquid bismuth that have been observed in experiments are investigated according to the ab initio molecular dynamics method using the SIESTA open software package. The density of electronic states, the radial distribution function of atoms, and the self-diffusion coefficient are calculated for the temperature range from the melting point equal to 545 K to 1500 K. The calculated data are in good agreement with the experimental data. It is found that the position of the first peak in the radial distribution function of atoms and the self-diffusion coefficient are characterized by a nonmonotonic dependence under the conditions of superheating by approximately 150 K above the melting temperature. In the authors’ opinion, this dependence feature is attributed to a change in the liquid short-range order structure.     

Coexistence of superfluid and solid helium in aerogel

The results of recent neutron scattering studies of solid helium in silica aerogel are discussed. Previously I.V. Kalinin et al., Pis’ma Zh. éksp. Teor. Fiz. 87 (1), 743 (2008) [JETP Lett. 87 (1), 645 (2008)], we detected the existence of a superfluid phase in solid helium at a temperature below 0.6 K and a pressure of 51 bar, although, according to the phase diagram, helium should be in the solid state under these conditions. This work is a continuation of the above studies whose main goal was to examine the detected phenomenon and to establish basic parameters of the existence of a superfluid phase. We have determined the temperature of the superfluid transition from solid to superfluid helium, T _C = 1.3 K, by analyzing experimental data. The superfluid phase excitation parameters (lifetime, intensity, and energy) have a temperature dependence similar to that of bulk helium. The superfluid phase coexists with the solid phase in the entire measured temperature range from T = 0.05 K to T _C and is a nonequilibrium one and disappears at T _C.     

Theoretical study of a melting curve for tin

The melting curve of Sn has been calculated using the dislocation-mediated melting model with the `zone-linking method'. The results are in good agreement with the experimental data. According to our calculation, the melting temperature of γ -Sn at zero pressure is about 436~K obtained by the extrapolation of the method from the triple point of Sn. The results show that this calculation method is better than other theoretical methods for predicting the melting curve of polymorphic material Sn.     

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.     

Role of thermal conduction in single-bubble cavitation

Effect of thermal conduction on radiation from a single cavitating bubble has been studied in a hydrochemical framework including variation of heat conductivity of noble gases up to 2500 K. Results of numerical simulation show that thermal conductivity plays an important role in determining ultimate cavitation temperature. Higher thermal conductivity of lighter noble gases causes to more thermal dissipation during the bubble collapse, leading to a lower peak temperature. Moreover, at the same driving conditions, radius of light emitting region is greater for heavier noble bubbles. Therefore, sonoluminescence radiation is more intensive from heavier noble gases. Phase diagrams of single-bubble sonoluminescence have also been calculated and in comparison with available experimental data, there is a relatively good agreement between the theory and experiment for noble gases.     

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⁻³.