束缚条件下H2的量子蒙特卡罗研究

束缚条件下H2的量子蒙特卡罗研究

论文摘要

超高压下固态氢的行为,一直是热点研究课题。分子束缚模型能够很好地描述压力对H2体系的影响,但是以往的研究都是在固定质子位置,忽略了非常重要的质子量子运动下进行的。在本文中,我们推广了质子固定的分子束缚模型,用量子理论同时描述电子和质子,完全顾及到质子的零点运动。用量子蒙特卡罗方法精确地求解了束缚下的H2体系四体问题,避免了Born-Oppenheimer绝热近似的使用,研究了压力作用下固态分子氢的基态性质。本文中使用两种不同的方法计算固态氢在高压下的状态方程,得到的结果与LeSar等人把质子固定在特殊点上的绝热近似的结果相接近。与以往的质子固定模型相比,由于质子的零点运动,体系的基态能量变大,平衡键长变长。随着束缚的增强,两质子间的距离逐渐减小,说明氢的键长逐渐变短。束缚椭球半长轴a与半短轴b的比ratio = ab随压力增加有变大的趋势,在密度为0.28 mol/cm3和0.35mol/cm3附近存在两个突变点,对应着固态氢在低温时110GPa和150GPa压力点处的两个相变。质子固定模型(该忽略了质子的零点运动)并没有出现类似的突变现象,说明质子的零点运动引发了a/b的突变。通过计算各向异性率(δ),发现随着压力的增加δ在逐渐的降低,在两个a/b突变点处,δ也发生突变。与实际晶体中的W-S元胞计算得到的δ值对比可以看出,第二相的变化与类Pa3结构的变化趋势相同,对确认长期争论的第二相结构问题的解决有重要的意义。束缚椭球模型计算中<θ2>1/2出现两次大的变化,可以推测与分子取向的变化密切相关,进一步说明了a/b的两个突变对应着固体氢中的两次压致取向相变。质子的零点运动(ΔRp/Rp)在这两个密度附近出现很强的涨落,更加确定a/b随密度变化出现的两个突变是由于质子的零点运动引起的结论。两次势能趋势的变化是引起a/b出现突变的驱动力。

论文目录

  • 提要
  • 第一章 绪论
  • 1.1 氢的基本性质及应用
  • 1.2 高压下氢的研究现状
  • 1.3 研究目的及意义
  • 第二章 量子蒙特卡罗方法(QMC)
  • 2.1 基本思想
  • 2.2 统计基础
  • 2.3 变分原理
  • 2.4 Markov 过程
  • 2.5 Metropolis 算法
  • 2.6 变分蒙特卡罗方法
  • 2.7 扩散蒙特卡罗方法
  • 2体系'>第三章 束缚条件下的H2体系
  • 3.1 模型
  • 3.2 试探波函数的选取
  • 2的结果和讨论'>第四章 量子蒙特卡罗方法研究压力下H2的结果和讨论
  • 4.1 状态方程的计算
  • 4.2 基态能量和平衡键长的变化
  • 4.3 压力下束缚形状的改变与压致结构相变
  • 4.3.1 压力下束缚形状的改变
  • 4.3.2 束缚椭球与Wigner-Seitz元胞之间的联系
  • 4.3.3 分子轴与椭球主轴夹角的变化
  • 4.4 压缩过程中质子量子效应的变化
  • 结论
  • 参考文献
  • 中文摘要
  • 英文摘要(Abstract)
  • 致谢
  • 发表的文章
  • 相关论文文献

    标签:;  ;  ;  ;  

    束缚条件下H2的量子蒙特卡罗研究
    下载Doc文档

    猜你喜欢