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

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

固氩高压物态方程和弹性性质的密度泛函理论计算

 用平面波赝势方法结合局域密度近似密度泛函理论（DFT-LDA）计算了零温下固态氩晶体在压力0~82 GPa的p-V关系和弹性性质,计算结果与静高压实验数据符合较好,计算结果表明局域密度近似方法能较好地描述固氩晶体高压下的性质,采取合理的方法和计算参数,惰性气体固态晶体高压下的力学性质可以比较准确地计算出来,这可为一些还不能通过实验进行研究的物态分析提供借鉴。     

固氪物态方程的关联量子化学计算

运用多体展开理论和量子化学方法—–超分子单、双(三重)激发微扰处理耦合簇CCSD(T)方法,首次系统地计算了面心立方固氪在较宽(从晶格平衡位置到体积压缩率超过3倍)区间的两体、三体和四体相互作用对结合能和物态方程的贡献大小,包括Hartree-Fock自洽场项和范德瓦耳斯长程关联作用项;并与实验数据进行比较.结果表明,在考虑到两体、三体、四体相互作用能后,多体展开理论以及CCSD(T)方法对平衡位置结合能测量数据0—130 GPa整个研究区间的实验物态方程数据都做出令人满意的描述.     

一种计算氩等离子物态方程的简单模型

将Thomas-Feimi统计模型电离势的数值结果进行函数逼近，给出一个便于近似数值求解的解析表达式和计算电离度的近似计算方法，计算了Ar元素LTE情况下的电离度和物态方程，结果与Saha模型的计算结果和实验结果符合较好.所提出的简单模型也适用于计算混合物物态方程，可以在电磁发射技术领域中的强电离等离子体中有广阔的应用前景. 关键词： 等离子体 物态方程 电离势 电离度     

应用Saha方程计算氩等离子体的Hugoniot物态方程

应用Saha模型加Debye-Huckel修正计算氩等离子体的Hugoniot物态方程,温度在10000-30000K之间,压力在0.0133-0.166GPa范围内,非理想参数Γ≈0.38.实验值和理论计算结果符合较好,说明Saha模型加Debye-Huckel修正能较好描述该热力学条件下氩等离子体的热力学行为.     

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

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

氩等离子体物态方程的TF简化模型

对于高温高压下氩等离子体的电离度和物态方程,本文给出了一种基于Thomas-Feimi(TF)统计模型的简化计算新方法:首先将TF模型电离势的数值结果进行函数逼近,给出一个便于数值求解的计算电离度的近似计算方法,并由此计算了局部热动平衡下的氩等离子体在10～1000 eV高温范围内的物态方程.计算结果与国外报道的其他几种理论模型的计算结果均符合很好,与实验值也吻合较好.本文所提出的简单模型也适用于计算混合物物态方程,可以在电磁发射技术领域中的强电离等离子体中有更为广阔的应用前景.     

混合物物态方程的计算

在采用体积相加原理计算混合物物态方程的基础上,建立了一种物理模型确定混合物温度。根据混合物中各组分温度和压强平衡条件,采用压强-密度迭代方法计算给出混合物物态方程,编制了两种组分的混合物物态方程计算程序。为检验建立的温度模型的合理性及程序的有效性,分析了不同密度、温度状态的氢（H2）和钨（W）组成的混合物状态参量,计算了以下情形及其组合情形的混合物物态方程：H2和W以不同质量比混合;质量比固定,单组分状态不同;温度区间和密度区间不同。研究表明:实际应用中在建立的混合物温度模型基础上确定的混合物物态方程是合理的。     

混合物物态方程的计算

 在采用体积相加原理计算混合物物态方程的基础上,建立了一种物理模型确定混合物温度。根据混合物中各组分温度和压强平衡条件,采用压强-密度迭代方法计算给出混合物物态方程,编制了两种组分的混合物物态方程计算程序。为检验建立的温度模型的合理性及程序的有效性,分析了不同密度、温度状态的氢（H2）和钨（W）组成的混合物状态参量,计算了以下情形及其组合情形的混合物物态方程：H₂和W以不同质量比混合;质量比固定,单组分状态不同;温度区间和密度区间不同。研究表明:实际应用中在建立的混合物温度模型基础上确定的混合物物态方程是合理的。     

基于第一性原理的固氦零点振动能的量子理论计算

基于第一性原理,运用原子团簇理论和Hartree-Fock Self-Consistent-Field从头算法计算了高压固氦hcp结构在平衡位置r=1.75-2.6Å附近的原子相互作用能,并使用1stOpt软件拟合出固氦的零点振动频率,得到零点振动能、Gruneisen系数及零点振动压强,最后与实验值做了比较.结果表明：固氦晶体在其各自平衡位置附近的运动可近似为简谐振动;固氦晶体的零点振动频率随着摩尔体积的增大而减小;零点振动能占多体相互总能的比例超过5%,已不可忽略,其比例随摩尔体积的增加而增大;零点振动能随原子间距和摩尔体积的增大而减小;Gruneisen系数随着摩尔体积的增大而增大,并始终处于1和2之间;零点振动引起的压强与实验值曲线趋势完全一致,但比实验值稍小,零点振动压强对总压强的贡献随着摩尔体积的增大由8%增大到50%左右,由此引起的零点振动压强值的贡献越来越大.     

基于第一性原理的固氦零点振动能的量子理论计算

基于第一性原理,运用原子团簇理论和Hartree-Fock Self-Consistent-Field从头算法计算了高压固氦hcp结构在平衡位置r=1.75-2.6附近的原子相互作用能,并使用1stOpt软件拟合出固氦的零点振动频率,得到零点振动能、Gruneisen系数及零点振动压强,最后与实验值做了比较.结果表明:固氦晶体在其各自平衡位置附近的运动可近似为简谐振动;固氦晶体的零点振动频率随着摩尔体积的增大而减小;零点振动能占多体相互总能的比例超过5%,已不可忽略,其比例随摩尔体积的增加而增大;零点振动能随原子间距和摩尔体积的增大而减小;Gruneisen系数随着摩尔体积的增大而增大,并始终处于1和2之间;零点振动引起的压强与实验值曲线趋势完全一致,但比实验值稍小,零点振动压强对总压强的贡献随着摩尔体积的增大由8%增大到50%左右,由此引起的零点振动压强值的贡献越来越大.     

绝热量子计算理论引介

本文引介绝热量子计算理论评述量子绝热定理最新的应用.基于量子绝热方法在最新的前沿领域量子计算中建立量子绝热算法.我们引介"局域量子绝热",并考虑了这一局域绝热概念对量子算法的可能应用.     

绝热量子计算理论引介

本文引介绝热量子计算理论评述量子绝热定理最新的应用。基于量子绝热方法在最新的前沿领域量子计算中建立量子绝热算法。我们引介"局域量子绝热",并考虑了这一局域绝热概念对量子算法的可能应用。     

第一性原理对金的高压相变和零温物态方程的计算

在密度泛函理论下，用缀加平面波加局域轨道方法，分别采用广义梯度近似(GGA)和局域密度近似(LDA)对金的面心立方晶格结构(fcc)、体心立方晶格结构(bcc)和六角密堆积结构(hcp)的结构能量进行了计算.在GGA下，计算得出fcc向hcp和hcp向bcc的相变分别发生在380GPa 和1250GPa；而LDA下相变分别发生在490GPa和790GPa.当计算压强达到2TPa时，bcc在这两种近似下仍然保持稳定的结构.根据不同体积下不同结构的电子态密度的特征，对发生相变的物理原因进行了定性的分析，在此基础上得到了金的零温状态方程. 关键词： 缀加平面波方法 固态相变 电子态密度 物态方程     

金属铝固液气完全物态方程研究

基于GRAY模型建立了金属铝的固液气三相完全物态方程,并与等温压缩线、Hugoniot线、熔化线以及零压热力学函数的实验结果进行对比,表明本物态方程可合理描述金属铝在宽广热力学空间的热力学状态.     

液氩冲击压缩特性的理论计算

用修正的WCA状态方程和Ross变分微扰理论计算了液氩冲击压缩曲线，计算中体系分子间相互作用势选择EXP-6有效两体势模型，计算结果与Nellis等人的实验数据符合较好。     

液氩冲击压缩特性的理论计算

用修正的WCA状态方程和Ross变分微扰理论计算了液氩冲击压缩曲线,计算中体系分子间相互作用势选择EXP-6有效两体势模型,计算结果与Nellis等人的实验数据符合较好.     

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

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

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

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