氢是宇宙中最简单的元素,由一个质子和一颗电子组成。然而,当科学家将这种气体挤压进极小的空间时,它便转化成了物理学中最难解的谜题之一。在巨大的压力下,氢分子发生分解,原子失去了各自的身份,材料开始像金属一样导电。几十年来,研究人员一直试图绘制出这一过程在三维世界(即我们生活的世界)中究竟是如何发生的。他们曾争论原子是会排列成整齐、重复的晶体模式,还是即使在极低温度下也会熔化成混沌的类液体状态。这个答案至关重要,因为理解这种转变有助于解释像木星这样的巨行星是如何构成的,甚至可能为寻找无电阻导电材料提供指导。尽管三维版本的氢已被研究了近一个世纪,但一个更简单的版本——当氢被挤压成一层扁平的二维薄片时会发生什么——在很大程度上仍未被探索。
一支研究团队现在将注意力转向了这种扁平的二维氢。他们并没有使用物理样本(因为在没有支撑面的情况下,不可能保持一层扁平的氢片),而是利用强大的计算机模拟来构建并测试数千种不同的氢原子排列方式。他们的目标是找到这些原子在高压下能够共同存在的最稳定方式。在物理学领域,“最稳定”的状态是指维持该状态所需能量最少的状态,科学家称之为基态。他们首先测试了许多不同的构想,包括每个原子都有特定、重复邻居的完美晶格,以及随机、混乱的原子堆。他们每次运行包含多达512个原子的模拟,这个数量足以观察材料是自然倾向于形成晶体,还是更倾向于保持无序状态。
结果令人惊讶。在压力中间范围内——即氢分子分解点与形成简单原子晶体点之间——计算机模拟显示,最稳定的状态并非晶体。相反,原子沉淀为一种无序、杂乱的排列,没有表现出长程重复模式。即便研究人员将原子核视为简单的经典粒子,忽略了通常使微观粒子表现得像波一样的奇异量子效应,这种情况依然发生了。这种无序状态如此稳定,以至于其能量低于团队测试过的任何完美晶体结构。这表明,在扁平的二维世界中,氢变成金属的过程采取了与三维世界不同的路径,跳过了整齐的晶体阶段,首先沉淀为一种杂乱且稳定的状态。
研究人员非常谨慎地确保这并非其计算机模型的误差产物。他们通过测试较小规模的原子组进行了验证,发现无序状态仅在他们使用较大的原子组(具体为70个原子以上)时才会出现。这证实了该结果是材料的真实属性,而非模拟规模造成的偏差。他们还通过测量结构中模式的强度来观察原子的排列情况。在完美的晶体中,随着原子的增加,这些模式会变得越来越强;但在无序氢中,这些模式增长得非常缓慢,证实了原子并未形成长程有序排列。即使团队加入了使原子振动和抖动的量子力学效应,无序状态依然保持最稳定,并变得更加受到青睐。
这一发现挑战了长期以来的假设,即压力下的物质最终必须形成晶体。虽然研究人员承认,不能排除存在一种模式如此庞大且复杂的晶体从而被其模拟所遗漏的可能性,但证据有力地指向了这种特定扁平形态氢的无序基态。这一发现表明,将材料限制在二维空间可以从根本上改变其行为,抑制了其自然结晶的倾向。它为理解受限环境如何影响物质开启了一扇新门,表明通过将氢“扁平化”,科学家可以创造出一种在三维世界中并不存在的稳定无序相。这项工作并未解决我们自身世界中高压氢的谜团,但它提供了一个清晰、具体的例子,展示了改变系统的形状如何能重写其自我组织的方式。
技术摘要:高密度二维氢的低对称性基态
问题陈述
高压下氢的行为仍然是凝聚态物理学中一个核心且悬而未决的问题,特别是关于分子解离、金属化以及异常熔融的机制。尽管三维(3D)相图已研究了近一个世纪,但关于分子到原子转变的性质,以及在解离附近是否存在中间液体或非晶态相的争论仍然存在。相比之下,二维(2D)氢的对应物在很大程度上仍未被探索,现有的研究大多局限于非弛豫的高对称性候选结构或极低压力的情形。鉴于维度降低可以从本质上改变键合、堆积与量子涨落之间的平衡,确定二维限制是否能稳定独特的基态,或者是否会改变金属化机制,是一个关键的开放性问题。
方法论
作者对二维氢(2DH)的零温相图进行了全面的从头算(ab initio)研究,涵盖了从分子相到原子相的转变过程。
- 计算框架: 研究采用了基于 PBE 能级的平面波密度泛函理论(DFT)。电子结构在三维空间中求解,而质子被限制在二维平面(x−y)内运动,并在 z 方向设置较大的层间距以模拟单层结构。
- 结构搜索: 作者采用了一种广泛的搜索策略,包括:
- 用于寻找候选结构的随机结构搜索(AIRSS)。
- 受候选 3D 晶体(如 $Cmca$-12, C2/c)启发的超胞构型。
- 具有各种晶格对称性(正方形、三角形、六角形)和取向的分子晶体。
- 原子晶体(正方形、三角形、蜂窝状、无序态)。
- 系统尺寸与有限尺寸效应: 模拟使用了包含多达 512 个原子的超胞。通过改变系统尺寸 N(从 48 到 512 个原子)进行了有限尺寸标度分析,以确保所观察到的相在热力学极限下的稳定性。
- 热力学: 相稳定性通过计算焓(H=E+pA)来确定,其中 A 是二维胞面积,p 是二维压力。研究考虑了经典原子核和核量子效应(NQE),后者通过路径积分分子动力学(PIMD)进行评估。
- 验证: 结果使用 CASTEP 软件包进行了交叉检查,并通过结构因子分析来区分晶体态与无序态。
主要贡献与结果
主要发现是:在高密度的二维氢中,存在一个无序基态,在一定的中间压力区间内,其焓值低于所有考虑过的晶体候选结构。
相图:
- 低压阶段: 基态为分子“人字形”(herringbone)相(分子位于三角形晶格上且取向平行),该相在 ≈0.02 Ha/Bohr2 以下保持稳定。
- 中间压力阶段(0.029(1)−0.054(4) Ha/Bohr2): 在此压力范围内,对应的维格纳-赛茨半径(rs)约为 1.17 至 0.99 Bohr,无序超胞结构成为热力学基态。这些从随机原子构型弛豫而来的结构并未表现出长程晶体有序性。在此窗口期内,它们的焓值低于晶体候选结构(包括 P62m、$Amm2$ 和蜂窝状晶格)。
- 高压阶段: 在 ≈0.054 Ha/Bohr2 以上,系统转变为晶体原子三角形相。
- 金属化: 在平均场 PBE 描述下,金属化(以费米能级处非零的态密度为标志)发生在从分子相向第一个原子结构(P62m)的转变过程中,这意味着中间无序相也具有金属性。
有限尺寸标度:
- 只有在原子数大于 ∼70 的超胞中,无序相的稳定性才得到确认。较小的系统(如 48 个原子)会产生定性不同的相图,其中晶体相更占优势,这凸显了利用大尺寸超胞来克服有限尺寸人工效应的必要性。
- 结构因子分析(Smax)表明,虽然晶体相的 Smax 随 N 线性增长,但无序相的增长显著缓慢,这为不存在长程有序或准晶相提供了严谨的证据。
核量子效应 (NQE):
- PIMD 模拟表明,质子量子涨落进一步稳定了无序相。当在无序相为经典基态的密度下模拟晶体 P62m 结构时,一旦引入 NQE,该结构会迅速熔变为无序相。
意义与主张
本文声称,高密度二维氢实现了一种独特的结构无序化路径,这种无序化即使在经典的 Born-Oppenheimer 能级下也是稳定的。这一发现表明,维度限制可以作为一种控制参数,用以抑制氢的结晶过程。
- 理论意义: 无序基态的存在挑战了“简单元素在零温下必然结晶”的假设。作者指出,虽然对于两体势存在结晶的数学证明,但 ab initio 能面本质上是多体性质的,因此允许此类例外。
- 与其他系统的区别: 该无序相不同于由零点运动克服微弱范德华力而产生的典型零温量子液体(如 3He 或 4He),也不同于在有限温度下终止的普通低熔点单质液体。它也不同于质子无序冰,后者通常是动力学受阻态而非真正的基态。
- 局限性与未来工作: 作者谦逊地承认,他们无法严格排除存在比其模拟超胞更大原胞的晶体基态的可能性(类似于碱金属中的不共度结构)。他们强调,其结论是基于 DFT-PBE 能量学的,而更精确的方法(如量子蒙特卡洛方法)是下一步的必要工作。此外,他们指出实验实现需要基底,而基底会改变能量学特性,并且包括量子质子效应在内的全面表征仍是未来的课题。
最终,本研究将高密度二维氢视为研究结晶与结构无序竞争关系的理想环境,并提出通过维度限制来设计低温液体的途径。
每周获取最佳 materials science 论文。
受到斯坦福、剑桥和法国科学院研究人员的信赖。
请查收邮箱确认订阅。
出了点问题,再试一次?
无垃圾邮件,随时退订。