高压固氩物态方程的量子理论计算

基于第一性原理,运用ab initio超分子单、双(三)重激发耦合簇理论(CCSD(T))和aug-cc-pVQZ基矢,精确计算了fcc晶体固氩在最近邻原子间距R=0.20~0.39nm时的两体、三体和四体结合能,零点振动能及物态方程。结果表明:固氩的多体势对结合能的贡献在高压区域是一正负交叉级数;零点振动能占多体相互总能的比例较小,但不可忽略;在高压区域,只考虑两体势时对固氩的压缩特性表现过硬,加入三体势后与实验结果在60GPa内完全吻合,考虑到四体势后对整个实验区间0~114GPa内做出令人满意的描述,且在压强达到114GPa时与实验值相差约3GPa,吻合程度达到97%。最后,通过与密度泛函理论的局域密度近似和广义梯度近似方法比较发现,泛函理论(DFT)只有在50~114GPa范围内与实验值符合较好,不如本研究所采用方法适用的压强范围宽。     

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

采用量子化学自洽场方法及原子团簇理论对固氦晶格原子势能进行多体展开，定量确定了最近邻原子间距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. 关键词： 固氦 状态方程 多体相互作用     

基于原子晶体构型的高压固氩中的多体相互作用

X射线衍射实验显示固氩是面心立方(fcc)晶格结构,目前对晶体氩的研究只限于两体,三体以及四体相互作用势.本文利用多体展开方法和超分子单、双(三)重激发耦合簇理论(CCSD(T))对固氩fcc晶格结构的三体和四体的几何构型、几何参数、不同体积下所有三体和四体构型的势能以及各构型所占比例等几个方面进行了准确的量子化学计算.结果表明:所有三体构型中对总的三体势能贡献最大的是构型1、构型6、构型12和构型23;三体势及其交换部分和色散部分的计算结果与现有解析经验势在长程部分符合得非常好,但在短程部分有较小差异.所有的四体构形中对总的四体势能贡献最大的是构型1,构型2,构型4,构型5,构型7和构型8;四体势及其交换势部分和色散部分的计算结果尚无解析经验势可比较.利用这些特殊构型的相关数据并结合其它构型,可拟合出更准确的三体经验势函数及其参数,也为拟合四体经验势函数及其参数提供了重要的参考价值.     

Ne-Xe体系势能曲线和光谱研究

本论文对Ne-Xe体系的势能曲线和束缚态能级作了系统的研究.采用耦合簇CCSD(T)方法和超分子近似,以及aug-cc-PVXZ(X=T,Q,5)基组,计算了Ne-Xe体系的相互作用势,并采用三种方案外推得到基底限值,给出了不同基底和外推方法下的势能曲线平衡位置和势阱深度.计算了Ne-Xe体系振转能级和各同位素基振动态的纯转动跃迁频率,及相应的光谱常数,并与实验结果进行了比较.     

基于原子晶体构型的高压固氩中的多体相互作用

X射线衍射实验显示固氩是面心立方(fcc)晶格结构,目前对晶体氩的研究只限于两体,三体以及四体相互作用势。本文利用多体展开方法和超分子单、双 (三)重激发耦合簇理论(CCSD(T))对固氩fcc晶格结构的三体和四体的几何构型、几何参数、不同体积下所有三体和四体构型的势能以及各构型所占比例等几个方面进行了准确的量子化学计算。结果表明:所有三体构型中对总的三体势能贡献最大的是构型1、构型6、构型12和构型23;对三体势及其交换部分和色散部分的计算结果与现有解析经验势在长程部分符合得非常好,但在短程部分有较小差异。所有的四体构形中对总的四体势能贡献最大的是构型1,构型2,构型4,构型5,构型7和构型8;对四体势及其交换势部分和色散部分的计算结果尚无解析经验势可比较。利用这些特殊构型的相关数据并结合其它构型,可拟合出更准确的三体经验势函数及其参数,也为拟合四体经验势函数及其参数提供了重要的参考价值。     

结合能的新势函数对高压稠密惰性元素压缩特性的影响

基于量子理论和原子团簇理论,运用多体展开方法和第一性原理的从头算方法,提出了一种计算稠密惰性元素(氦、氖、氩和氪)原子结合能的新势函数,运用新公式研究了结合能对高压稠密惰性元素高压压缩特性的影响。此公式引入了一个物理参量β(其值为0.5),使得势函数的表达形式更加简单、准确。对比结果表明,结合能的新势函数能够准确地描述多体相互作用对结合能的贡献,且平均相对误差在5%以内。结合能的新势函数对压缩特性的影响在当前实验压强范围内(氦60 GPa、氖238 GPa、氩114 GPa、氪128 GPa)做出了令人满意的描述,且与实验值及理论计算结果基本完全吻合,平均相对误差在3%以内。最后,以固氩的压强数据为例,验证了势函数的准确性。该势函数不仅适用于更宽密度和更高压强范围,而且对所有惰性元素原子各种状态的结合能、高压压缩特性、定容比热容、熔化曲线和弹性模量的研究具有重要的指导意义。     

一种适用于原子团簇的多体展开新方案: 以氮团簇为例

多体展开方法虽然已经广泛地用于估算弱相互作用体系的能量,但是其并不适用于计算共价团簇和金属团簇的能量. 本文提出了一种适用于计算共价体系能量的相互作用多体展开(IMBE)方法. 在相互作用多体展开方法中,体系的能量表示为孤立原子的能量及该原子与其他周围原子间相互作用的和. 首先将该方法应用于计算氮团簇的能量,且多体展开截断至四体项. 结果表明：与传统的多体展开方法相比,相互作用多体展开方法可以显著地降低能量误差. 另外,以密度泛函理论计算结果为参考,相互作用多体展开方法估算能量的误差不依赖于体系的大小和结构,说明相互作用多体展开方法比较适合用于估算共价相互作用大体系的能量.     

一种适用于原子团簇的多体展开新方案:以氮团簇为例（英文）

多体展开方法虽然已经广泛地用于估算弱相互作用体系的能量,但是其并不适用于计算共价团簇和金属团簇的能量.本文提出了一种适用于计算共价体系能量的相互作用多体展开(IMBE)方法.在相互作用多体展开方法中,体系的能量表示为孤立原子的能量及该原子与其他周围原子间相互作用的和.首先将该方法应用于计算氮团簇的能量,且多体展开截断至四体项.结果表明:与传统的多体展开方法相比,相互作用多体展开方法可以显著地降低能量误差.另外,以密度泛函理论计算结果为参考,相互作用多体展开方法估算能量的误差不依赖于体系的大小和结构,说明相互作用多体展开方法比较适合用于估算共价相互作用大体系的能量.     

碱金属Li原子与CO分子体系相互作用势从头计算研究

采用三重激发项校正CCSD的CCSD(T)方法和AUG-CC-PVnZ(n=2,3,4)基组,优化了CO基态分子结构,并计算了碱金属Li原子与CO分子的相互作用势,共1010个构型势能点得到体系的势能面.结果表明:同一方法下,不同基组得到的CO基态分子的键长、能量等均与实验符合很好.Li-CO势能面体现较小的各向异性势,存在两个势阱,且都为非严格T型结构;基组重叠误差(BSSE)对相互作用势的影响比较明显,采用CP方法(Counterposie method)消除基组重叠误差,不同基组计算的相互作用势显示较好一致.体系各项异性势势阱值远高于CO基态分子转动常数,使得碰撞产生强烈的非弹性碰撞,这将阻止协同冷却制备超冷分子.     

基于Ne-HBr从头算CCSD(T)势的散射截面

利用非线性最小二乘法拟合在CCSD(T)/aug-cc-pVQZ理论水平下计算的相互作用能,得到了基态Ne-HBr复合物相互作用势的解析表达式.基于拟合的从头算CCSD(T)势,通过收敛的密耦计算得到了入射能量分别为40,60,80和100 meV时,Ne-HBr散射的微分截面和分波截面,详细讨论了态-态转动激发截面对总非弹性散射截面的影响和散射截面随能量的变化趋势.希望研究结果对该体系的散射实验和进一步的理论研究能提供参考信息.     

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

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

Efimov physics in heteronuclear four-body systems

We study three- and four-body Efimov physics in a heteronuclear atomic system with three identical heavy bosonic atoms and one light atom. We show that exchange of the light atom between the heavy atoms leads to both three- and four-body features in the low-energy inelastic rate constants that trace to the Efimov effect. Further, the effective interaction generated by this exchange can provide an additional mechanism for control in ultracold experiments. Finally, we find that there is no true four-body Efimov effect-that is, no infinite number of four-body states in the absence of two- and three-body bound states-resolving a decades-long controversy.     

Model Study of Three-Body Forces in the Three-Body Bound State

The Faddeev equation for the three-body bound state with two- and three-body forces is solved directly as three-dimensional integral equation. The numerical feasibility and stability of the algorithm, which does not employ partial wave decomposition is demonstrated. The three-body binding energy and the full wave function are calculated with Malfliet-Tjon-type two-body potentials and scalar two-meson exchange three-body forces. For two- and three- body forces of ranges and strengths typical of nuclear forces the single-particle momentum distribution and the two-body correlation function are similar to the ones found for realistic nuclear forces.     

Bound-state calculations of the three-dimensional Yakubovsky equations with the inclusion of three-body forces

The four-body Yakubovsky equations in a three-dimensional approach with the inclusion of the three-body forces are proposed. The four-body bound state with two- and three-body interactions is formulated in the three-dimensional approach for identical particles as a function of vector Jacobi momenta, specifically, the magnitudes of the momenta and the angles between them. The modified three-dimensional Yakubovsky integral equations are successfully solved with the scalar two-meson exchange three-body force, where the Malfliet-Tjon-type two-body force is implemented. The three-body force effects on the energy eigenvalue and the four-body wave function, as well as accuracy of our numerical calculations are presented.     

Three-body interaction-induced light scattering in krypton gas: a computer simulation of the spectral line shapes

The two-, three- and four-body effective collision induced scattering spectral line shapes are calculated for dense gaseous krypton using the pairwise additivity (PA) approximation and different polarizability models. These spectra and several interaction induced spectra calculated at various densities are compared with the experimental measurements of Barocchi et al. [1988, Europhys. Lett., 5, 607]. The potential effect on the spectrum is found to be weak. The results obtained with the Meinander et al. [1986, J. chem. Phys., 84, 3005] empirical polarizability model and molecular dynamics fit well the experimental two- and three-body spectral shapes. The irreducible contribution to the spectral shape is evaluated using the dipole induced dipole irreducible polarizability [buckingham, A. D., and Hands, I. D., 1991, Chem. Phys. Lett., 185, 544]. This contribution is found to be relatively weak for the anisotropic spectra in the frequency and density range studied, explaining the good agreement between the pairwise approximation calculations and the experimental data. The spectra radiated by the quasi-molecules Kr₂, Kr₃, and Kr₄ (the total spectrum within the PA approximation) are also simulated.     

Crystal Structure Transformations of Rare-Gas Solids Under Pressure

Pressure-induced structural changes in solid krypton (Kr) and xenon (Xe) have been studied using angle dispersive X-ray diffraction in a diamond-anvil cell (DAC) up to 50 GPa. The analysis of the results shows that in solid Kr (Xe) the phase transition from fcc to hcp starts below 3.2 GPa (1.5 GPa). Albeit the hcp/fcc ratio increases under pressure, both phases coexist up to the highest pressure reached in this study. Room temperature (RT) equations of state (EOS) are determined.     

Ab initio calculations of energy transfer and non-additivity in the He-Ne laser system

Ab initio calculations were performed for several suggested mechanisms of energy transfer between helium metastable particles and neon. Optimized geometries and excited-state energies were calculated for neon excited-state complexes and the convergence properties of the non-additive contributions to the interaction energies were examined. The most probable excitation-transfer mechanism was found to be based on an energy difference of 0.0674 eV between the triplet excited state of and the singlet excited state of . No theoretical evidence was found for the production of neon singlet excited-state complexes other than 20.0858 to 20.4875 eV by the considered two-, three- and four-body models of energy transfer processes. The energy curves of the reactions involving the excited-state complexes and are provided and compared with the previously reported experimental results on the reaction . The relation between the probability of energy transfer and laser activity is discussed. The non-additive contribution to the total interaction energy of the nominated intermediate complex was found to be negligible, pointing to the possibility of constructing model potentials and simulation of larger systems. Received: 15 December 1998 / Received in final form: 20 March 1999     

固氖高压物态方程的量子理论计算

 运用Hartree-Fock计算方法和原子团簇理论,研究了高压下面心立方晶体氖中的二体相互作用、三体相互作用对中心原子势能及冷能的贡献。结合零点振动能,计算得到了固氖在4~208 GPa压强范围内的等温物态方程。结果表明,在较高压强区域内,等温压缩线与高压实验数据符合得非常好;在压强超过100 GPa的区域内,计算得到的等温压缩线较目前已有的理论研究结果系统地改善并提高了3%~5%,精确地解释了固氖的高压实验数据。