搜索

x

留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

基于路径积分分子动力学与热力学积分方法的高压氢自由能计算

陈基 冯页新 李新征 王恩哥

引用本文:
Citation:

基于路径积分分子动力学与热力学积分方法的高压氢自由能计算

陈基, 冯页新, 李新征, 王恩哥

A fully quantum description of the free-energy in high pressure hydrogen

Chen Ji, Feng Ye-Xin, Li Xin-Zheng, Wang En-Ge
PDF
导出引用
  • 在相图研究中, 严格计算一个真实系统在特定温度、压强下的自由能是近年来该领域理论方法发展的前沿. 自Mermin提出有限温度密度泛函理论后, 在电子结构层面, 弱关联系统中人们就其在对自由能贡献的描述已相对完善, 但在原子核运动的描述上, 热运动与量子运动的非简谐项却总被忽视. 本文将路径积分分子动力学与热力学积分结合, 对300 GPa下氢晶体Cmca 结构中原子核热涨落与量子涨落对自由能的影响进行了分析. 发现在100 K核量子涨落非简谐项的贡献约为15 meV每原子, 远大于不同结构间静态焓的差别. 该研究提醒人们简谐近似在核量子效应描述中可能存在的不准确性(即使在低温下). 同时, 我们采取的方法 也为人们进行自由能的准确计算提供了一个简单有效的手段.
    Hydrogen is the lightest and most abundant element in the universe. Ever since Wigner and Huntington's prediction that pressure induced metallization might happen in solid hydrogen, understanding the hydrogen phase diagram has become one of the greatest challenges in condensed matter and high pressure physics. The light mass of hydrogen means that the nuclear quantum effects could be important in describing this phase diagram under high pressures. Numerical evaluations of their contributions to the structural, vibrational, and energetic properties, however, are difficult and up to now most of the theoretical simulations still remain classical. This is particularly true for the energetic properties. When the free-energies of different phases are compared in determining the ground state structure of the system at a given pressure and temperature, most of the theoretical simulations remain classical. When nuclear quantum effects must be taken into account, one often resorts to the harmonic approximation. In the very rare case, the anharmonic contributions from the nuclear statistical effects are considered by using a combination of the thermodynamic integration and the at initio molecular dynamics methods, which helps to include the classical nuclear anharmonic effects. Quantum nuclear anharmonic effects, however, are completely untouched. Here, using a self-developed combination of the thermodynamic integration and the at initio path-integral molecular dynamics methods, we calculated the free-energies of the high pressure hydrogen at 100 K from 200 GPa to 300 GPa. The harmonic lattice was taken as the reference and the Cmca phase of the solid hydrogen was chosen. When the bead number of the path-integral (P) equals one, our approach reaches the so-called classical limit. Upon increasing P until the results are converged, our approach reaches the limit when both classical and quantum nuclear anharmonic effects are included. Therefore, by comparing the free-energy of the harmonic lattice and the thermodynamic integration results at P equals one, we isolate the classical nuclear anharmonic effects. By comparing the thermodynamic integration results at P equals one and with those when they are converged with respect to P, we isolate the quantum nuclear anharmonic effects in a very clean manner. Our calculations show that the classical nuclear anharmonic contributions to the free-energy are negligible at this low temperature. Those contributions from the quantum nuclear anharmonic effects, however, are as large as ~15 meV per atom. This value also increases with pressure. This study presents an algorithm to quantitatively calculate the quantum contribution of the nuclear motion to free-energy beyond the often used harmonic approximation. The large numbers we got obtained also indicate that such quantum nuclear anharmonic effects are important in describing the phase diagram of hydrogen, at/above the pressures studied.
      通信作者: 李新征, xzli@pku.edu.cn
    • 基金项目: 国家自然科学基金(批准号: 11422431, 11275008, 11274012, 91021007)和中国博士后科学基金(批准号: 2014M550005)资助的课题.
      Corresponding author: Li Xin-Zheng, xzli@pku.edu.cn
    • Funds: Project supported by the National Natural Science Foundation of China (Grant Nos. 11422431, 11275008, 11274012, 91021007), and the National Postdoc Research Foundation of China (Grant No. 2014M550005).
    [1]

    Mermin N D 1985 Phys. Rev. A 137 1441

    [2]

    Gillan M J 1989 J. Phys.: Condens. Matter 1 689

    [3]

    Wentzcovitch R M, Martins J L, Allen P B 1992 Phys. Rev. B 45 11372

    [4]

    Xu G, Ming W, Yao Y, Dai X, Zhang S C, Fang Z 2008 EPL 82 67002

    [5]

    Lee P A, Nagaosa N, Wen X G 2006 Rev. Mod. Phys. 78 17

    [6]

    Pickard C J, Needs R J 2007 Nat. Phys. 3 473

    [7]

    Li X Z, Wang E G 2014 Computer Simulations of Molecules and Condensed Matters: From Electronic Structures to Molecular Dynamics (Beijing: Peking University Press) pp134-140

    [8]

    Frenkel D, Lekkerkerker H N W, Stroobants A 1988 Nature 332 822

    [9]

    Meijer E J, Frenkel D 1991 J. Chem. Phys. 94 2269

    [10]

    Alfé D, Gillan M J, Price G D 1999 Nature 401 462

    [11]

    Alfé D, Price G D, Gillan M J 2001 Phys. Rev. B 64 045123

    [12]

    Wigner E, Huntington H B 1935 J. Chem. Phys. 3 764

    [13]

    Babaev E, Sudbo A, Ashcroft N W 2004 Nature 431 666

    [14]

    Bonev S A, Schwegler E, Ogitsu T, Galli G 2004 Nature 431 669

    [15]

    Deemyad S, Silvera I F 2008 Phys. Rev. Lett. 100 155701

    [16]

    Li X Z, Walker B, Probert M I J, Pickard C J, Needs R J, Michaelides A 2013 J. Phys.: Condens. Matter 25 085402

    [17]

    Chen J, Li X Z, Zhang Q F, Probert M I J, Pickard C J, Needs R J, Michaelides A, Wang E G 2013 Nat. Commun. 4 2064

    [18]

    Mao H K, Hemley R J 1994 Rev. Mod. Phys. 66 671

    [19]

    McMahon J M, Morales M A, Pierleoni C, Ceperley D M 2012 Rev. Mod. Phys. 84 1607

    [20]

    Zha C S, Liu Z X, Hemley R J 2012 Phys. Rev. Lett. 108 146402

    [21]

    Liu H Y, Zhu L, Cui W W, Ma Y M 2012 J. Chem. Phys. 137 074501

    [22]

    Perez A, von Lilienfeld O A 2011 J. Chem. Theory Comput. 7 2358

    [23]

    Habershon S, Manolopoulos D E 2011 J. Chem. Phys. 135 224111

    [24]

    Kresse G, Furthmller J 1996 Phys. Rev. B 54 11169

    [25]

    Feng Y X, Chen J, Alfè D, Li X Z, Wang E G 2015 J. Chem. Phys. 142 064506

    [26]

    Alfé D 2009 Comput. Phys. Commun. 180 2622

    [27]

    Perdew J P, Burke K, Ernzerhof M 1996 Phys. Rev. Lett. 77 3865

  • [1]

    Mermin N D 1985 Phys. Rev. A 137 1441

    [2]

    Gillan M J 1989 J. Phys.: Condens. Matter 1 689

    [3]

    Wentzcovitch R M, Martins J L, Allen P B 1992 Phys. Rev. B 45 11372

    [4]

    Xu G, Ming W, Yao Y, Dai X, Zhang S C, Fang Z 2008 EPL 82 67002

    [5]

    Lee P A, Nagaosa N, Wen X G 2006 Rev. Mod. Phys. 78 17

    [6]

    Pickard C J, Needs R J 2007 Nat. Phys. 3 473

    [7]

    Li X Z, Wang E G 2014 Computer Simulations of Molecules and Condensed Matters: From Electronic Structures to Molecular Dynamics (Beijing: Peking University Press) pp134-140

    [8]

    Frenkel D, Lekkerkerker H N W, Stroobants A 1988 Nature 332 822

    [9]

    Meijer E J, Frenkel D 1991 J. Chem. Phys. 94 2269

    [10]

    Alfé D, Gillan M J, Price G D 1999 Nature 401 462

    [11]

    Alfé D, Price G D, Gillan M J 2001 Phys. Rev. B 64 045123

    [12]

    Wigner E, Huntington H B 1935 J. Chem. Phys. 3 764

    [13]

    Babaev E, Sudbo A, Ashcroft N W 2004 Nature 431 666

    [14]

    Bonev S A, Schwegler E, Ogitsu T, Galli G 2004 Nature 431 669

    [15]

    Deemyad S, Silvera I F 2008 Phys. Rev. Lett. 100 155701

    [16]

    Li X Z, Walker B, Probert M I J, Pickard C J, Needs R J, Michaelides A 2013 J. Phys.: Condens. Matter 25 085402

    [17]

    Chen J, Li X Z, Zhang Q F, Probert M I J, Pickard C J, Needs R J, Michaelides A, Wang E G 2013 Nat. Commun. 4 2064

    [18]

    Mao H K, Hemley R J 1994 Rev. Mod. Phys. 66 671

    [19]

    McMahon J M, Morales M A, Pierleoni C, Ceperley D M 2012 Rev. Mod. Phys. 84 1607

    [20]

    Zha C S, Liu Z X, Hemley R J 2012 Phys. Rev. Lett. 108 146402

    [21]

    Liu H Y, Zhu L, Cui W W, Ma Y M 2012 J. Chem. Phys. 137 074501

    [22]

    Perez A, von Lilienfeld O A 2011 J. Chem. Theory Comput. 7 2358

    [23]

    Habershon S, Manolopoulos D E 2011 J. Chem. Phys. 135 224111

    [24]

    Kresse G, Furthmller J 1996 Phys. Rev. B 54 11169

    [25]

    Feng Y X, Chen J, Alfè D, Li X Z, Wang E G 2015 J. Chem. Phys. 142 064506

    [26]

    Alfé D 2009 Comput. Phys. Commun. 180 2622

    [27]

    Perdew J P, Burke K, Ernzerhof M 1996 Phys. Rev. Lett. 77 3865

  • [1] 胡庭赫, 李直昊, 张千帆. 元素掺杂对储氢容器用高强钢性能影响的第一性原理和分子动力学模拟. 物理学报, 2024, 73(6): 067101. doi: 10.7498/aps.73.20231735
    [2] 张宇航, 李孝宝, 詹春晓, 王美芹, 浦玉学. 单层MoSSe力学性质的分子动力学模拟研究. 物理学报, 2023, 72(4): 046201. doi: 10.7498/aps.72.20221815
    [3] 田礼漫, 温永立, 王云飞, 张善超, 李建锋, 杜镜松, 颜辉, 朱诗亮. 路径积分传播子测量的研究进展. 物理学报, 2023, 72(20): 200305. doi: 10.7498/aps.72.20230902
    [4] 第伍旻杰, 胡晓棉. 单晶Ce冲击相变的分子动力学模拟. 物理学报, 2020, 69(11): 116202. doi: 10.7498/aps.69.20200323
    [5] 付宝勤, 侯氢, 汪俊, 丘明杰, 崔节超. 钨空位捕获氢及其解离过程的分子动力学. 物理学报, 2019, 68(24): 240201. doi: 10.7498/aps.68.20190701
    [6] 李杰杰, 鲁斌斌, 线跃辉, 胡国明, 夏热. 纳米多孔银力学性能表征分子动力学模拟. 物理学报, 2018, 67(5): 056101. doi: 10.7498/aps.67.20172193
    [7] 董琪琪, 胡海豹, 陈少强, 何强, 鲍路瑶. 水滴撞击结冰过程的分子动力学模拟. 物理学报, 2018, 67(5): 054702. doi: 10.7498/aps.67.20172174
    [8] 范航, 聂福德, 龙瑶, 陈军. 钝感高能炸药三氨基三硝基苯高温高压下热力学性质的分子动力学模拟研究. 物理学报, 2016, 65(6): 066201. doi: 10.7498/aps.65.066201
    [9] 张宝玲, 宋小勇, 侯氢, 汪俊. 高密度氦相变的分子动力学研究. 物理学报, 2015, 64(1): 016202. doi: 10.7498/aps.64.016202
    [10] 王成龙, 王庆宇, 张跃, 李忠宇, 洪兵, 苏折, 董良. SiC/C界面辐照性能的分子动力学研究. 物理学报, 2014, 63(15): 153402. doi: 10.7498/aps.63.153402
    [11] 常旭. 多层石墨烯的表面起伏的分子动力学模拟. 物理学报, 2014, 63(8): 086102. doi: 10.7498/aps.63.086102
    [12] 周化光, 林鑫, 王猛, 黄卫东. Cu固液界面能的分子动力学计算. 物理学报, 2013, 62(5): 056803. doi: 10.7498/aps.62.056803
    [13] 马颖. 非晶态石英的变电荷分子动力学模拟. 物理学报, 2011, 60(2): 026101. doi: 10.7498/aps.60.026101
    [14] 邵建立, 王 裴, 秦承森, 周洪强. 铁冲击相变的分子动力学研究. 物理学报, 2007, 56(9): 5389-5393. doi: 10.7498/aps.56.5389
    [15] 罗 晋, 祝文军, 林理彬, 贺红亮, 经福谦. 单晶铜在动态加载下空洞增长的分子动力学研究. 物理学报, 2005, 54(6): 2791-2798. doi: 10.7498/aps.54.2791
    [16] 王海龙, 王秀喜, 梁海弋. 应变效应对金属Cu表面熔化影响的分子动力学模拟. 物理学报, 2005, 54(10): 4836-4841. doi: 10.7498/aps.54.4836
    [17] 谢文贤, 徐 伟, 雷佑铭, 蔡 力. 随机参激和外激联合作用下非线性动力系统的路径积分解. 物理学报, 2005, 54(3): 1105-1112. doi: 10.7498/aps.54.1105
    [18] 王 平, 杨新娥, 宋小会. 具有含时平方反比项的谐振子的路径积分求解. 物理学报, 2003, 52(12): 2957-2960. doi: 10.7498/aps.52.2957
    [19] 梁海弋, 王秀喜, 吴恒安, 王宇. 纳米多晶铜微观结构的分子动力学模拟. 物理学报, 2002, 51(10): 2308-2314. doi: 10.7498/aps.51.2308
    [20] 吴恒安, 倪向贵, 王宇, 王秀喜. 金属纳米棒弯曲力学行为的分子动力学模拟. 物理学报, 2002, 51(7): 1412-1415. doi: 10.7498/aps.51.1412
计量
  • 文章访问数:  6760
  • PDF下载量:  673
  • 被引次数: 0
出版历程
  • 收稿日期:  2015-05-15
  • 修回日期:  2015-06-29
  • 刊出日期:  2015-09-05

/

返回文章
返回