共振价键波函数在高压液氢量子蒙卡模拟中的适用性研究-物理学报.PDFVIP

物理学报 Acta Phys. Sin. Vol. 62, No. 15 (2013) 156101 共振价键波函数在高压液氢量子蒙卡模拟中的 适用性研究 12† 2 2 2 2 12 李名锐 周刚 初哲 戴湘晖 吴海军 范如玉 1) ( 清华大学工程物理系, 北京 100084 ) 2) ( 西北核技术研究所, 西安 710024 ) ( 2013 年3 月10 日收到; 2013 年5 月2 日收到修改稿) 在共振价键理论基础上, 选取高压液氢电子主要占据轨道的线性组合作为基组, 构建由Jastrow 项和反对称孪生 函数乘积项(AGP) 组成的波函数. 考虑电子关联作用的共振价键(RVB) 波函数得出的能量值低于LDA 能量值; 当 满足rs 175 或T 15000 K 时引入backflow 项以改善波函数结点面, 改善后的能量值下降约1 mHa/atom, 能量方 差值变小. 将构建的RVB 波函数与电子- 离子耦合的蒙特卡罗法(CEIMC) 相结合, 计算结果与实验及其他ab-initio 结果相符合, 获得的液氘单次冲击Hugoniot 曲线基本通过所有加载类型实验误差棒, 液氘在50.3 GPa 处具有最大压 缩率4.48, 在 100—120 GPa 内未发现压缩率有急剧增大的现象. 构建的RVB 波函数能够适用于较宽密度与温度范 围内(10 rs 22, 2800 K T 60000 K) 液氢的模拟, 与CEIMC 法相结合可提高液氢冲击特性的模拟精度. 关键词: 共振价键理论, 波函数, 量子蒙卡法, 液氢 PACS: 61.20.Ja, 71.15.−m, 02.70.Ss, 05.70.Ce DOI: 10.7498/aps.62.156101 法, 如Sesame15 , Ross16 , Saumon-Chabrier17 , 流体 1 引言 变分法(FVT)18 以及活度展开法(ACTEX)19 等, 当前普遍认为Sesame72 模型与Ross 模型预测了 氢是最简单也是宇宙中含量最丰富的元素. 由 100 GPa 以下液氢压缩率的上下限, 但在更高压力 1 2 于在现代天体物理 、惯性约束核聚变(ICF) 以 范围内化学模型法的适用性较差. 随着计算机技 3 及氢金属化理论 等方面中的重要作用, 需要在较 术的发展及研究的不断深入, 人们更多采取的是 宽的压力及温度范围内准确掌握氢的相关特性, 高 基于第一性原理的ab-initio 方法. Car 和Parrinello 压液氢状态成为了实验与理论研究最需解决的关 等利用密度泛函理论(DFT) 将经典的分子动力学 键问题之一. 1995 年Nellis 等利用二级轻气炮多次 20 方法(MD) 推广到第一性原理的计算中 , 为研 冲击压缩法, 在 1