非静水压强下固体分子氢的结构与电子性质演化

宋贤齐; 刘畅; 刘子恺; 王建云; 李全

doi:10.11858/gywlxb.20230720

非静水压强下固体分子氢的结构与电子性质演化

doi: 10.11858/gywlxb.20230720

宋贤齐^{1, 2,},
刘畅^{1, 2},
刘子恺^{1, 2},
王建云^{1, 2,*, ,},
李全^{1, 2,*, ,}

1.
吉林大学物理学院物质模拟方法与软件教育部重点实验室, 吉林长春　130012
2.
吉林大学物理学院超硬材料国家重点实验室, 吉林长春　130012

基金项目: 国家自然科学基金（T2325013，52288102，52090024，12074140，12202158）；国家重点研发计划（2021YFA1400503，2018YFA0703404）；中国博士后科学基金（2023M731296）

详细信息

作者简介:
宋贤齐（1994－），男，博士，主要从事极端高压下计算凝聚态物理研究.E-mail：sxq@calypso.cn

通讯作者:
王建云（1992－），女，硕士，工程师，主要从事极端高压下计算凝聚态物理研究.E-mail：wjy@calypso.cn

李　全（1980－），男，博士，教授，主要从事极端高压下计算凝聚态物理研究.E-mail：liquan777@jlu.edu.cn

中图分类号: O521.2
计量
- 文章访问数: 345
- HTML全文浏览量: 114
- PDF下载量: 127
出版历程
- 收稿日期: 2023-08-16
- 修回日期: 2023-09-01
- 网络出版日期: 2023-10-09
- 刊出日期: 2023-11-07

Structural and Electronic Properties of Solid Hydrogen at Non-Hydrostatic Pressures

SONG Xianqi^{1, 2
,},
LIU Chang^{1, 2},
LIU Zikai^{1, 2},
WANG Jianyun^{1, 2,*
, ,},
LI Quan^{1, 2,*
, ,}

1.
Key Laboratory of Material Simulation Methods & Software of Ministry of Education, College of Physics, Jilin University, Changchun 130012, Jilin, China
2.
State Key Laboratory of Superhard Materials, College of Physics, Jilin University, Changchun 130012, Jilin, China

摘要

摘要: 固体氢金属化所需的压强超过400 GPa，这种超高压条件给固体氢的实验制备和表征带来了极大的挑战。为此，利用第一性原理计算方法，系统研究了固体分子氢在非静水压强下的结构和物性演化。研究发现，在非静水高压条件下，固体分子氢具有良好的结构稳定性。非静水压条件将导致固体氢晶格的对称性破缺和电荷的重新分布，使固体分子氢在较低压强下（如压强低于300 GPa）转变为金属和超导体。据此提出了引入各向异性非静水压环境从而在较低压强下获得金属氢和高温超导氢的新思路。
- 非静水压 /
- 金属氢 /
- 超导体 /
- 第一性原理计算
Abstract: The pressure required for the metallization of solid hydrogen exceeds 400 GPa, thereby presenting a formidable challenge for its experimental preparation and characterization. Here, we systematically explore the structures and properties undergone by solid hydrogen under non-hydrostatic pressure conditions by first-principle calculations. Our findings reveal that solid molecular hydrogen can retain good structural stability under non-hydrostatic pressure conditions, which induces symmetry breaking and charge redistribution within the solid hydrogen lattice, facilitating the transformation of solid molecular hydrogen into metallic and superconducting states at lower pressures (e.g., pressures are lower than 300 GPa). This study proposes a new idea of introducing an anisotropic non-hydrostatic pressure environment for achieving metallic hydrogen at lower pressure.
- non-hydrostatic pressure /
- metallic hydrogen /
- superconductor /
- first-principles calculations

HTML全文

碳纤维复合材料（carbon-fibre reinforced composites，CFRP）具有轻质、高强、高韧、耐腐蚀、耐高温、可设计性强等优异特性，在轨道交通、新能源、火箭导弹以及航空航天（如飞机蒙皮、高空飞艇、发动机等）等领域具有重要应用价值^[1–2]。然而，在某些应用领域中，碳纤维复合材料面临极端服役环境，从而对其性能提出了极其严苛的要求。例如，作为飞机蒙皮或高空飞艇中的关键材料，碳纤维复合材料时常遭遇低温环境和结冰问题，不仅引起飞行器表面平整度、外形、气动布局的改变^[3–4]以及重量的增加，导致不可预知的飞行问题，而且由于纤维与基体的热膨胀系数差异较大，在低温环境中还可能出现收缩不协调导致复合材料内部产生裂纹，进而引发灾难性破坏^[5]。因此，调控碳纤维复合材料在低温环境下的力学性能，并实现其在低温环境下的除冰功能是非常必要的，这对于保障碳纤维复合材料结构件，尤其是应用于高空飞行器结构件中的碳纤维复合材料服役安全性具有重大意义。

在低温除冰方面，当前最常用的机翼除冰方法有热气除冰、化学除冰、机械除冰和电热除冰等^[6–8]。热气除冰主要利用发动机的热气加热进行除冰，但是可能会降低发动机的推力。化学除冰是通过喷涂防冰液等进行除冰，例如，Xue等^[9]开发出聚二甲基硅氧烷(PDMS)/TiO₂纳米复合材料，将其作为牺牲层喷涂在基材上，一旦在超疏水表面形成冰，便会导致牺牲层自动脱离，达到除冰效果，但是这种方法存在涂层寿命较短、需要重复喷涂的问题。电加热除冰是将加热垫布置在需要保护的部件表面，以代替热气除冰，在其服役期间可重复利用。然而，现在的电加热除冰技术尚不成熟，在民用飞机上的应用较少，尤其是机翼表面很难布置加热垫，因此需要进一步深入研究。

不锈钢超薄带，又名手撕钢（hand-torn steel，HTS），是一种强度高、柔韧性良好、导电率高的金属材料，在航天航空等高精尖领域具有重要应用^[10]。手撕钢是一种各向同性材料，将手撕钢与碳纤维材料复合，有望弥补各向异性碳纤维层合板在某些方向上的性能短板，从而有效地增强碳纤维复合材料的综合力学性能。此外，手撕钢作为即时发热材料，将其与碳纤维材料复合不仅能够实现其整体温度调控，改善材料在低温服役环境下的力学性能，而且通过发挥其致热功能，可能实现飞机机翼等结构件表面的快速除冰。将碳纤维与手撕钢相结合，形成碳纤维-金属复合材料，发挥纤维和金属的优点，不仅可以实现飞行器的轻量化设计，还能够有效避免飞行器因表面结构件结冰而造成的安全事故。

目前，学者们已经对纤维-金属层合板进行了相关研究^[11–13]。例如，李艳^[11]采用阳极氧化法和热压成型法制备碳纤维金属铝层合板和碳纤维金属钛层合板，结果表明碳纤维金属钛层合板相比碳纤维金属铝层合板具有更好的力学性能，这是由添加金属种类所决定的。马其华等^[12]制备了金属-CFRP混合管，发现CFRP与金属管相互协调可提升整体的承载能力，改善整个混合结构的能量吸收能力。研究表明，温度对纤维复合材料的力学性能有一定程度的影响^[14–15]。然而，已有研究主要针对纤维-金属复合材料在常温环境下的力学性能，对于纤维-金属复合材料的电流驱动温度调节方案，及其在电流致热效应调控下的力学性能和除冰功能尚缺乏深入研究。

为此，本研究通过设计和裁剪不同形状的手撕钢，采用铺层固化方法制备“碳纤维预浸料”、“碳纤维预浸料-未裁剪手撕钢”和“碳纤维预浸料-弓字形手撕钢”3种复合材料，通过通电致热试验、三点弯曲性能测试和通电致热除冰测试，探索电流通路和电流强度引起的致热效应对复合材料力学性能和除冰功能的影响，以期为低温服役环境下碳纤维复合材料的力学性能优化和除冰功能实现提供指导，研究结果有望在航天航空等领域得到重要应用。

1.   样品制备

1.1   主要实验原料

本研究采用的CFRP为T300单向碳纤维增强环氧树脂薄层预浸料，型号为USN 20000（碳纤维与环氧树脂的质量比为6∶4），购自威海光威集团有限责任公司；手撕钢超薄带购自上海宏伟有限公司。材料参数如表1所示。

表 1 材料的基本属性

Table 1. Basic material properties

Material Tensile strength/MPa Density/(g·cm⁻³) Thickness/mm

HTS 1130 7.26 0.02

CFRP 1800 1.80 0.20

下载: 导出CSV
| 显示表格

1.2   三点弯曲样品制备

为了研究电流致热效应对材料力学性能的影响，将加热层与支撑层合并，即采用2层碳纤维预浸料和1层手撕钢交替叠层铺设，共11层，如图1(a)所示，其中预浸料采用单向铺层。加热层手撕钢起支撑和通电加热作用，另外两层手撕钢则主要起支撑作用。

图 1 三点弯曲层合板铺层设计方案

Figure 1. Lay-up design scheme of three-point bending laminate

下载: 全尺寸图片幻灯片

根据三点弯曲测试标准《纤维金属层板弯曲性能试验方法》（GB/T 41049—2021），将样品尺寸设计为15 mm×80 mm，对加热层手撕钢进行裁剪处理，形成宽度约为3 mm的“弓字形”电流加热通路，如图1(b)所示。采用铺层固化方法制备样品：(1) 将碳纤维预浸料和手撕钢裁剪成合适的尺寸并解冻6 h；(2) 按照铺层方案逐层铺设样品；(3) 将样品进行模压和抽真空消泡处理，并在干燥箱中高温固化，干燥箱中的温度先以1～2 ℃/min的速率由室温升至120 ℃并保温90 min，然后以1～2 ℃/min的速率降至60 ℃以下进行脱模处理。所有样品均采用相同的制备工艺。最终制备出“碳纤维预浸料”、“碳纤维预浸料-未裁剪手撕钢”和“碳纤维预浸料-弓字形手撕钢”3组样品，如图1(c)所示。

1.3   通电致热除冰样品制备

按照图2(a)所示形状将手撕钢裁剪成“弓字形”，以研究电流通路宽度（a）对致热效应的影响。为获取加热层温度，按照图1(a)中第3～11层方案，制备尺寸为100 mm×100 mm的层合板，如图2(b)～图2(f)所示，其中a分别为5、10、15和20 mm。此外，为了验证除冰功能，按照图1(a)中第1～11层方案制备电流通路宽度为10 mm的样品。

图 2 不同电流通路宽度的手撕钢层合板

Figure 2. Hand-torn steel laminates with different current path widths

下载: 全尺寸图片幻灯片

2.   实验测试

对制备的样品进行测试，先对碳纤维-弓字形手撕钢长方形层合板通入不同的电流，测试其力学性能，然后对正方形层合板进行通电致热和除冰功能测试。

2.1   通电致热和三点弯曲测试

为研究电流强度对材料力学性能的影响，将样品接入电路中通电加热，如图3(a)所示，通过调压器（型号TDGC2-0.5 kVA）控制输出电压，调控样品中的电流大小（见表2），利用红外热成像仪（型号FLIR C2，测量误差±2 ℃，测量范围−20～180 ℃）采集样品表面温度数据，当样品表面温度达到稳态后，采用石棉对样品进行保温处理，然后断电，并立即在万能试验机（型号STD5000）上进行三点弯曲性能测试（如图3(b)所示），并与碳纤维预浸料-未裁剪手撕钢样品进行对比。测试中，跨距为40 mm，加载速率为2 mm/min，每组样品进行3次测试，最终测试结果为多次测量的平均值。

图 3 三点弯曲样品通电实验和力学性能测试

Figure 3. Three-point bending sample energization test and mechanical property test

下载: 全尺寸图片幻灯片

表 2 三点弯曲样品通电方案

Table 2. Three-point bending sample energization scheme

Samples Current/A

CFRP No power
HTS-CFRP No power
Cut HTS-CFRP 0, 0.5, 1.0, 1.5

下载: 导出CSV
| 显示表格

2.2   电流致热除冰升温特性测试

为测试碳纤维预浸料-弓字形手撕钢复合材料的电流致热除冰功能，将样品接入如图3(a)所示的电路中，调节调压器，使样品中的电流依次增大（间隔约为0.5 A），通电方案见表3，利用红外热成像仪测量样品表面温度，记录样品在不同电流下的稳态温度。此外，对电流通路宽度为10 mm的11层除冰验证样品的升温过程和除冰功能进行测试。

表 3 电流致热样品通电方案

Table 3. Energizing samples with current heating

Width of hand-torn steel/mm Current/A

20 0, 0.5, 1.0, 1.5, 2.0, 2.5
15 0, 0.5, 1.0, 1.5, 2.0, 2.5
10 0, 0.5, 1.0, 1.5, 2.0, 2.5
5 0, 0.5, 1.0, 1.5, 1.8

下载: 导出CSV
| 显示表格

3.   结果与讨论

3.1   电流通路和电流强度对复合材料力学性能的影响

对3种复合材料在未通电条件下进行三点弯曲性能测试，结果如图4所示。本研究主要考虑比强度和比吸能两方面因素，其中，比强度是材料的强度与其表观密度之比，而比吸能E_SA的计算公式为

图 4 3种样品的三点弯曲测试结果

Figure 4. Three-point bending test results of three kinds of samples

下载: 全尺寸图片幻灯片

$E_{\rm {SA}} = \int_0^{{x_0}} \frac {F}{m}{\text{d}}x$
(1)

式中：F为三点弯曲载荷；x为压头位移；x₀为破坏位移（载荷下降30%对应的压头位移）；m为样品的表观密度，纯碳纤维预浸料、碳纤维预浸料-未裁剪手撕钢、碳纤维预浸料-弓形手撕钢复合材料的表观密度分别为1.80、2.03、2.02 g/cm³。图4(a)表明，3种复合材料的应力-应变曲线变化趋势相同，达到弯曲强度后，载荷迅速下降，复合材料发生破坏。由图4(b)和图4(c)可知：与纯碳纤维预浸料复合材料相比，碳纤维预浸料-未裁剪手撕钢复合材料的比强度有所提升（由0.817 N·m/kg提升至0.876 N·m/kg），但是比吸能有所降低（由412.9 J/kg降低为359.2 J/kg），碳纤维预浸料-弓字形手撕钢复合材料的比强度虽然有所降低（由0.817 N·m/kg降为0.733 N·m/kg），但是破坏应变显著增大，延缓了复合材料的破坏，从而使其比吸能略有提升（由412.9 J/kg提升为414.6 J/kg）。对比碳纤维预浸料-未裁剪手撕钢复合材料与碳纤维预浸料-弓字形手撕钢复合材料，可以看出，后者的比强度降低，但比吸能提高。这说明加入未裁剪的手撕钢可显著提升碳纤维预浸料复合材料的强度，加入裁剪的手撕钢可在一定程度上提升碳纤维预浸料复合材料的韧性。下面将对该实验结果进行进一步的机理分析。

在纯碳纤维预浸料复合材料的基础上添加未裁剪手撕钢后，比强度得到提升，而添加裁剪手撕钢后比强度有所降低，这可能是由于：(1) 手撕钢具有良好的强度和柔性，有助于提升复合材料的整体强度；(2) 与碳纤维预浸料不同的是（沿纤维方向的承载能力显著高于垂直纤维方向的承载能力），手撕钢是各向同性材料，能够很好地承载并传递层间应力，从而使其整体强度得到提升；(3) 裁剪后的手撕钢由于相邻通路间含有间隙，间隙处无法承受和传递载荷，且材料质量增加，致使其比强度有所降低。

从吸能方面来看，碳纤维预浸料-弓字形手撕钢层合板的比吸能显著提升，而碳纤维预浸料-未裁剪手撕钢层合板的比吸能有所降低，其原因可能是：(1) 对于碳纤维预浸料-未裁剪手撕钢层合板，虽然强度得到提高，但是极限应变减低（手撕钢的极限应变低于碳纤维预浸料的极限应变），且质量增加导致其比吸能降低；(2) 对于碳纤维预浸料-弓字形手撕钢层合板，裁剪手撕钢留下间隙，此处的碳纤维预浸料将在固化过程中结合到一起，从而在材料破坏过程中形成更多的断裂面，提高材料的能量吸收特性。

上述分析表明，碳纤维预浸料-弓字形手撕钢层合板与纯碳纤维预浸料层合板相比，比强度降低，比吸能显著提升，即可以在一定程度上延缓层合板的破坏，同时吸收更多能量以发挥其承载和防护作用。下面对碳纤维预浸料-手撕钢复合材料在不同电流作用下的力学性能做进一步研究。

分别对碳纤维预浸料-弓字形手撕钢层合板在不同电流强度下的力学性能进行测试，并与碳纤维预浸料-未裁剪手撕钢层合板（不通电）的力学性能进行对比分析，结果如图5所示。可以看出，通电后碳纤维预浸料-弓字形手撕钢层合板的比强度和比吸能提高，但是随着电流强度增加，比强度和比吸能降低。通电不仅能够提升碳纤维预浸料-手撕钢复合材料的强度，而且增大破坏时的极限应变，使其在破坏过程中吸收更多的能量，例如：通入0.5 A的电流后，比强度大幅提升（由0.733 N·m/kg提高至0.957 N·m/kg），甚至超过了碳纤维预浸料-未裁剪手撕钢层合板，且变形伸长量较大，通电后碳纤维预浸料-弓字形手撕钢层合板的比吸能显著提升（由414.6 J/kg降低为579.4 J/kg），相对于纯碳纤维预浸料层合板，比强度和比吸能分别提升17.2% 和40.3%。这可能是由于通电后样品温度上升，达到一定温度时预浸料中的环氧树脂基体发生链段重排并进一步固化，导致环氧树脂分子链的交联密度增加^[14]，增强了碳纤维与树脂之间的结合强度，从而提高了结构整体的强度和能量吸收能力。但是，随着电流强度的增加，样品温度过高会引起部分分子链断裂^[15]，从而降低层合板的强度和能量吸收特性。因此，需要选择合适的电流，以满足实际应用中对层合板强度和吸能特性的要求。

图 5 碳纤维-手撕钢层合板在不同电流强度下的弯曲测试结果

Figure 5. Bending test results of carbon fibre-hand-torn steel laminates at different current strengths

下载: 全尺寸图片幻灯片

综上所述，在纯碳纤维预浸料中加入手撕钢会提升材料的力学性能，添加经过裁剪处理的手撕钢后，与由纯碳纤维预浸料制备的样品相比，虽然比强度有所下降，但是比吸能显著提升。通电对于碳纤维-手撕钢复合材料的力学性能也有重要影响，通过调控电流强度可实现强度和能量吸收特性的优化。

3.2   碳纤维-手撕钢层合板的电流致热除冰功能

为了实现碳纤维-手撕钢层合板电流致热除冰功能，测试得到具有不同电流通路宽度（“弓”字形手撕钢宽度）的碳纤维-手撕钢层合板在通电后的稳态温度，结果如图6所示。由图6(a)可知，随着电流通路宽度的减小，表面温度的上升趋势存在一定的差异。其中：碳纤维预浸料-未裁剪手撕钢样品通入电流后，温度几乎不变，可能是由于整片手撕钢的电阻很小，导致发热功率很小，样品表面升温过程和散热过程几乎处于动态平衡状态，因此表面温度基本维持在某一稳定值；碳纤维预浸料-弓字形手撕钢样品的电阻显著增大，通电后样品表面温度随着电流增加而升高。

图 6 不同电流通路宽度层合板表面稳态温度随电流变化和温度场分布

Figure 6. Steady-state temperature variation of laminate surface with current and temperature field distribution for different current path widths

下载: 全尺寸图片幻灯片

此外，还发现当电流一定时，随着电流通路宽度的减小，样品表面稳态温度逐渐升高，样品表面温度场逐渐均匀。当电流达到2.5 A时，电流通路宽度为20、15、10 mm的样品的表面温度分别达到65.7、78.5、101.0 ℃，如图6(b)～图6(d)所示。需要说明的是，当电流通路宽度为5 mm时，未能完成电流为2.5 A的实验，这是因为随着电流的增大，温度迅速上升，当电流为1.6 A时，样品表面温度达到152.0 ℃，如图6(e)所示，当电流达到1.8 A时，温度超过180 ℃，且一段时间后样品发生变形损坏。图6(b)～图6(e)的结果表明：当手撕钢宽度为15和20 mm时，温度场分布均较为分散，且温度提升缓慢，不利于大面积除冰；随着电流通路宽度减小，样品表面温度场趋于均匀；当手撕钢宽度为5和10 mm时，温度场分布较为均匀，可以实现大面积除冰。结合3.1节的力学性能分析，可以发现：碳纤维预浸料-弓字形手撕钢层合板中的裂纹或缺陷增加，导致力学性能下降；电流通路为10 mm的碳纤维-手撕钢层合板相比于其他层合板能够在保证良好力学性能的情况下，达到大面积除冰效果。

在2.5 A的电流下，电流通路宽度为10 mm的样品的加热功率为2 kW/m²，其温度变化过程如图7(a)～图7(f)所示。可见，接通电路后，样品表面温度迅速提升到49 ℃，如图7(a)所示，但初始热量分布不均匀；随着时间的延长，样品表面温度变化逐渐减缓；40 s时最高温度达到86.6 ℃，如图7(d)所示，样品表面高温区域趋于均匀；120 s时，最高温度基本稳定在94.9 ℃，样品中间区域的温度相对均匀，更有利于大范围除冰。此外，将边长约为15 mm的冰块置于样品表面，大约5 min后冰块完全融化，表明手撕钢铺层材料具有很好的除冰功能。

图 7 电流通路宽度为10 mm的层合板通入2.5 A电流后不同时刻样品表面的温度场

Figure 7. Temperature field on the surface of a laminate with a current path width of 10 mm after applying 2.5 A current at different moments

下载: 全尺寸图片幻灯片

上述针对碳纤维-手撕钢复合材料的电流致热除冰功能研究表明，不同电流通路宽度对复合材料电流致热功能具有重要影响，电流通路宽度减小，复合材料在相同电流下的致热效应提升，可用来通电致热除冰。此外，电流强度对于碳纤维-手撕钢复合材料的致热效应也有重要影响，增加电流强度会引起复合材料温度过高，可能引起材料破坏失效。通过调控手撕钢尺寸，在碳纤维-手撕钢复合材料的力学性能与除冰功能之间找到平衡点。实验表明，手撕钢宽度为10 mm的碳纤维-手撕钢复合材料可以在不降低材料性能的基础上实现通电除冰的目的。

4.   结论与展望

设计和制备了碳纤维预浸料、碳纤维预浸料-手撕钢、碳纤维预浸料-弓字形手撕钢3种复合材料，通过通电实验和三点弯曲测试，获得了手撕钢电流通路和电流强度对复合材料力学性能的影响，据此设计了具有不同电流通路宽度的碳纤维-手撕钢复合材料，研究了其电流致热除冰功能，得到如下结论。

(1) 在纯碳纤维预浸料层合板中加入未裁剪手撕钢会提升比强度，加入弓字形手撕钢会在比强度有所降低的情况下提升能量吸收特性。电流强度对碳纤维-手撕钢复合材料的力学性能也有重要影响，施加电流时，其比强度和比吸能显著提升。因此，在实际应用过程中，可通过调节电流强度优化复合材料的强度和能量吸收特性，以满足不同领域对复合材料力学性能的要求，从而为该类复合材料的应用提供积极指导。

(2) 碳纤维-手撕钢复合材料的电流致热除冰特性与其中手撕钢的电流通路宽度有密切联系。随着电流通路宽度的减小，相同电流强度下复合材料的致热温度及其均匀性提高，但是随着电流强度提升，致热温度过高会引起层合板破坏失效。通过调控手撕钢的尺寸，可以平衡其力学性能和除冰功能。相对于其他尺寸，电流通路宽度为10 mm的碳纤维-手撕钢层合板在保证材料性能和安全性的前提下，实现通电致热和快速除冰。这对于碳纤维-金属复合材料的结构设计和功能应用提供了参考依据，并有望在航天航空等领域得到应用。

本研究对于电流强度和电流通路宽度变化下复合材料的测试次数较少，实验环境为室温，且仅对比了复合材料的强度和吸能，研究结果具有一定的局限性。后续研究将细化电流通路宽度和电流强度，同时对刚度等其他重要力学性能进行分析，结合实际应用场景，增加风洞试验，研究复合材料在不同条件下的力学性能和致热效果，实现金属-碳纤维层合板力学性能和功能特性优化。

图 1 常压和高压条件下固体单质的金属化

Figure 1. Metallization of solid elementary substance at atmospheric and high pressures

下载: 全尺寸图片幻灯片

图 2 分子氢和原子氢的结构和电子轨道分布示意图

Figure 2. Structures and electron orbital distribution of molecular and atomic hydrogen

下载: 全尺寸图片幻灯片

图 3 固体氢候选相的焓-压强关系与代表性晶体结构

Figure 3. Enthalpy-pressure relations and typical crystal structures of solid hydrogen candidate phases

下载: 全尺寸图片幻灯片

图 4 高压下固体分子氢C2/c相沿(001)[010]晶向的剪切应力-应变曲线

Figure 4. Shear stress-strain curves of solid molecular hydrogen C2/c phase along the (001)[010] direction under high pressure

下载: 全尺寸图片幻灯片

图 5 250 GPa下分子氢C2/c相发生应变前后的电荷转移与超导体转变示意图

Figure 5. Electron transfer and superconducting transition of undeformed and deformed molecular hydrogen at 250 GPa

下载: 全尺寸图片幻灯片

图 6 350 GPa静水压下固体分子氢Cmca-12相的剪切应力-应变曲线

Figure 6. Shear stress-strain curves of solid molecular hydrogen Cmca-12 phase along various directions at 350 GPa

下载: 全尺寸图片幻灯片

参考文献(74)

[1]	WIGNER E, HUNTINGTON H B. On the possibility of a metallic modification of hydrogen [J]. The Journal of Chemical Physics, 1935, 3(12): 764–770. doi: 10.1063/1.1749590
[2]	ASHCROFT N W. Metallic hydrogen: a high-temperature superconductor? [J]. Physical Review Letters, 1968, 21(26): 1748–1749. doi: 10.1103/PhysRevLett.21.1748
[3]	BALL P. Metallic hydrogen in the spotlight [J]. Nature Materials, 2017, 16(3): 288. doi: 10.1038/nmat4872
[4]	DIAS R P, SILVERA I F. Observation of the Wigner-Huntington transition to metallic hydrogen [J]. Science, 2017, 355(6326): 715–718. doi: 10.1126/science.aal1579
[5]	EREMETS M I, DROZDOV A P, KONG P P, et al. Semimetallic molecular hydrogen at pressure above 350 GPa [J]. Nature Physics, 2019, 15(12): 1246–1249. doi: 10.1038/s41567-019-0646-x
[6]	LOUBEYRE P, OCCELLI F, DUMAS P. Synchrotron infrared spectroscopic evidence of the probable transition to metal hydrogen [J]. Nature, 2020, 577(7792): 631–635. doi: 10.1038/s41586-019-1927-3
[7]	MONACELLI L, ERREA I, CALANDRA M, et al. Black metal hydrogen above 360 GPa driven by proton quantum fluctuations [J]. Nature Physics, 2021, 17(1): 63–67. doi: 10.1038/s41567-020-1009-3
[8]	SONG X Q, LIU C, LI Q, et al. Stress-induced high- T_c superconductivity in solid molecular hydrogen [J]. Proceedings of the National Academy of Sciences of the United States of America, 2022, 119(26): e2122691119. doi: 10.1073/PNAS.2122691119
[9]	ZHANG L J, WANG Y C, LV J, et al. Materials discovery at high pressures [J]. Nature Reviews Materials, 2017, 2(4): 17005. doi: 10.1038/natrevmats.2017.5
[10]	李全, 马琰铭. 典型双原子分子晶体的高压解离和单原子相[J]. 高压物理学报, 2013, 27(3): 313–324. doi: 10.11858/gywlxb.2013.03.001 LI Q, MA Y M. High pressure dissociation of typical diatomic molecular solids and their atomic phases [J]. Chinese Journal of High Pressure Physics, 2013, 27(3): 313–324. doi: 10.11858/gywlxb.2013.03.001
[11]	DALLADAY-SIMPSON P, BINNS J, PEÑA-ALVAREZ M, et al. Band gap closure, incommensurability and molecular dissociation of dense chlorine [J]. Nature Communications, 2019, 10(1): 1134. doi: 10.1038/s41467-019-09108-x
[12]	EREMETS M I, GAVRILIUK A G, TROJAN I A, et al. Single-bonded cubic form of nitrogen [J]. Nature Materials, 2004, 3(8): 558–563. doi: 10.1038/nmat1146
[13]	WANG X L, WANG Y C, MIAO M S, et al. Cagelike diamondoid nitrogen at high pressures [J]. Physical Review Letters, 2012, 109(17): 175502. doi: 10.1103/PhysRevLett.109.175502
[14]	JI C, ADELEKE A A, YANG L X, et al. Nitrogen in black phosphorus structure [J]. Science Advances, 2020, 6(23): eaba9206. doi: 10.1126/sciadv.aba9206
[15]	DUAN D F, LIU Z T, LIN Z Y, et al. Multistep dissociation of fluorine molecules under extreme compression [J]. Physical Review Letters, 2021, 126(22): 225704. doi: 10.1103/PhysRevLett.126.225704
[16]	BARDEEN J, COOPER L N, SCHRIEFFER J R. Microscopic theory of superconductivity [J]. Physical Review, 1957, 106(1): 162–164. doi: 10.1103/PhysRev.106.162
[17]	BARDEEN J, COOPER L N, SCHRIEFFER J R. Theory of superconductivity [J]. Physical Review, 1957, 108(5): 1175–1204. doi: 10.1103/PhysRev.108.1175
[18]	PICKARD C J, NEEDS R J. Structure of phase Ⅲ of solid hydrogen [J]. Nature Physics, 2007, 3(7): 473–476. doi: 10.1038/nphys625
[19]	PICKARD C J, MARTINEZ-CANALES M, NEEDS R J. Density functional theory study of phase Ⅳ of solid hydrogen [J]. Physical Review B, 2012, 85(21): 214114. doi: 10.1103/PhysRevB.85.214114
[20]	MCMINIS J, CLAY III R C, LEE D, et al. Molecular to atomic phase transition in hydrogen under high pressure [J]. Physical Review Letters, 2015, 114(10): 105305. doi: 10.1103/PhysRevLett.114.105305
[21]	DALLADAY-SIMPSON P, HOWIE R T, GREGORYANZ E. Evidence for a new phase of dense hydrogen above 325 gigapascals [J]. Nature, 2016, 529(7584): 63–67. doi: 10.1038/nature16164
[22]	MONSERRAT B, DRUMMOND N D, DALLADAY-SIMPSON P, et al. Structure and metallicity of phase Ⅴ of hydrogen [J]. Physical Review Letters, 2018, 120(25): 255701. doi: 10.1103/PhysRevLett.120.255701
[23]	MAO H K, HEMLEY R J. Ultrahigh-pressure transitions in solid hydrogen [J]. Reviews of Modern Physics, 1994, 66(2): 671–692. doi: 10.1103/RevModPhys.66.671
[24]	LORENZANA H E, SILVERA I F, GOETTEL K A. Orientational phase transitions in hydrogen at megabar pressures [J]. Physical Review Letters, 1990, 64(16): 1939–1942. doi: 10.1103/PhysRevLett.64.1939
[25]	HEMLEY R J, MAO H K. Phase transition in solid molecular hydrogen at ultrahigh pressures [J]. Physical Review Letters, 1988, 61(7): 857–860. doi: 10.1103/PhysRevLett.61.857
[26]	HOWIE R T, GUILLAUME C L, SCHELER T, et al. Mixed molecular and atomic phase of dense hydrogen [J]. Physical Review Letters, 2012, 108(12): 125501. doi: 10.1103/PhysRevLett.108.125501
[27]	JI C, LI B, LIU W J, et al. Ultrahigh-pressure isostructural electronic transitions in hydrogen [J]. Nature, 2019, 573(7775): 558–562. doi: 10.1038/s41586-019-1565-9
[28]	MEIER T, LANIEL D, PENA-ALVAREZ M, et al. Nuclear spin coupling crossover in dense molecular hydrogen [J]. Nature Communications, 2020, 11(1): 6334. doi: 10.1038/s41467-020-19927-y
[29]	ASHCROFT N W. Hydrogen dominant metallic alloys: high temperature superconductors? [J]. Physical Review Letters, 2004, 92(18): 187002. doi: 10.1103/PhysRevLett.92.187002
[30]	LI Y W, HAO J, LIU H Y, et al. The metallization and superconductivity of dense hydrogen sulfide [J]. The Journal of Chemical Physics, 2014, 140(17): 174712. doi: 10.1063/1.4874158
[31]	DROZDOV A P, EREMETS M I, TROYAN I A, et al. Conventional superconductivity at 203 Kelvin at high pressures in the sulfur hydride system [J]. Nature, 2015, 525(7567): 73–76. doi: 10.1038/nature14964
[32]	DUAN D F, LIU Y X, TIAN F B, et al. Pressure-induced metallization of dense (H₂S)₂H₂ with high- T_c superconductivity [J]. Scientific Reports, 2014, 4: 6968. doi: 10.1038/srep06968
[33]	PENG F, SUN Y, PICKARD C J, et al. Hydrogen clathrate structures in rare earth hydrides at high pressures: possible route to room-temperature superconductivity [J]. Physical Review Letters, 2017, 119(10): 107001. doi: 10.1103/PhysRevLett.119.107001
[34]	LIU H Y, NAUMOV I I, HOFFMANN R, et al. Potential high- T_c superconducting lanthanum and yttrium hydrides at high pressure [J]. Proceedings of the National Academy of Sciences of the United States of America, 2017, 114(27): 6990–6995. doi: 10.1073/pnas.1704505114
[35]	SOMAYAZULU M, AHART M, MISHRA A K, et al. Evidence for superconductivity above 260 K in lanthanum superhydride at megabar pressures [J]. Physical Review Letters, 2019, 122(2): 027001. doi: 10.1103/PhysRevLett.122.027001
[36]	DROZDOV A P, KONG P P, MINKOV V S, et al. Superconductivity at 250 K in lanthanum hydride under high pressures [J]. Nature, 2019, 569(7757): 528–531. doi: 10.1038/s41586-019-1201-8
[37]	WANG H, TSE J S, TANAKA K, et al. Superconductive sodalite-like clathrate calcium hydride at high pressures [J]. Proceedings of the National Academy of Sciences of the United States of America, 2012, 109(17): 6463–6466. doi: 10.1073/pnas.1118168109
[38]	MA L, WANG K, XIE Y, et al. High-temperature superconducting phase in clathrate calcium hydride CaH₆ up to 215 K at a pressure of 172 GPa [J]. Physical Review Letters, 2022, 128(16): 167001. doi: 10.1103/PhysRevLett.128.167001
[39]	SUN Y, LV J, XIE Y, et al. Route to a superconducting phase above room temperature in electron-doped hydride compounds under high pressure [J]. Physical Review Letters, 2019, 123(9): 097001. doi: 10.1103/PhysRevLett.123.097001
[40]	LI B, JI C, YANG W G, et al. Diamond anvil cell behavior up to 4 Mbar [J]. Proceedings of the National Academy of Sciences of the United States of America, 2018, 115(8): 1713–1717. doi: 10.1073/pnas.1721425115
[41]	卢志鹏, 祝文军, 刘绍军, 等. 非静水压条件下铁从 α到 ε结构相变的第一性原理计算[J]. 物理学报, 2009, 58(3): 2083–2089. doi: 10.7498/aps.58.2083 LU Z P, ZHU W J, LIU S J, et al. Structure phase transition from α to ε in Fe under non-hydrostatic pressure: an ab initio study [J]. Acta Physica Sinica, 2009, 58(3): 2083–2089. doi: 10.7498/aps.58.2083
[42]	CHENG C. Uniaxial phase transition in Si: ab initio calculations [J]. Physical Review B, 2003, 67(13): 134109. doi: 10.1103/PhysRevB.67.134109
[43]	GAÁL-NAGY K, STRAUCH D. Transition pressures and enthalpy barriers for the cubic diamond→ β-tin transition in Si and Ge under nonhydrostatic conditions [J]. Physical Review B, 2006, 73(13): 134101. doi: 10.1103/PhysRevB.73.134101
[44]	DANG C Q, LU A L, WANG H Y, et al. Diamond semiconductor and elastic strain engineering [J]. Journal of Semiconductors, 2022, 43(2): 021801. doi: 10.1088/1674-4926/43/2/021801
[45]	DANG C Q, CHOU J P, DAI B, et al. Achieving large uniform tensile elasticity in microfabricated diamond [J]. Science, 2021, 371(6524): 76–78. doi: 10.1126/science.abc4174
[46]	LIU C, SONG X Q, LI Q, et al. Smooth flow in diamond: atomistic ductility and electronic conductivity [J]. Physical Review Letters, 2019, 123(19): 195504. doi: 10.1103/PhysRevLett.123.195504
[47]	LIU C, SONG X Q, LI Q, et al. Superconductivity in compression-shear deformed diamond [J]. Physical Review Letters, 2020, 124(14): 147001. doi: 10.1103/PhysRevLett.124.147001
[48]	LIU C, SONG X Q, LI Q, et al. Superconductivity in shear strained semiconductors [J]. Chinese Physics Letters, 2021, 38(8): 086301. doi: 10.1088/0256-307X/38/8/086301
[49]	KRESSE G, FURTHMÜLLER J. Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set [J]. Physical Review B, 1996, 54(16): 11169. doi: 10.1103/PhysRevB.54.11169
[50]	KRESSE G, JOUBERT D. From ultrasoft pseudopotentials to the projector augmented-wave method [J]. Physical Review B, 1999, 59(3): 1758–1775. doi: 10.1103/PhysRevB.59.1758
[51]	PERDEW J P, BURKE K, ERNZERHOF M. Generalized gradient approximation made simple [J]. Physical Review Letters, 1996, 77(18): 3865–3868. doi: 10.1103/PhysRevLett.77.3865
[52]	MONKHORST H J, PACK J D. Special points for Brillouin-zone integrations [J]. Physical Review B, 1976, 13(12): 5188–5192. doi: 10.1103/PhysRevB.13.5188
[53]	ZHANG Y, SUN H, CHEN C F. Superhard cubic BC₂N compared to diamond [J]. Physical Review Letters, 2004, 93(19): 195504. doi: 10.1103/PhysRevLett.93.195504
[54]	PAN Z C, SUN H, CHEN C F. Colossal shear-strength enhancement of low-density cubic BC₂N by nanoindentation [J]. Physical Review Letters, 2007, 98(13): 135505. doi: 10.1103/PhysRevLett.98.135505
[55]	PAN Z C, SUN H, CHEN C F. Indenter-angle-sensitive fracture modes and stress response at incipient plasticity [J]. Physical Review B, 2009, 79(10): 104102. doi: 10.1103/PhysRevB.79.104102
[56]	PAN Z C, SUN H, ZHANG Y, et al. Harder than diamond: superior indentation strength of wurtzite BN and lonsdaleite [J]. Physical Review Letters, 2009, 102(5): 055503. doi: 10.1103/PhysRevLett.102.055503
[57]	BARONI S, DE GIRONCOLI S, DAL CORSO A, et al. Phonons and related crystal properties from density-functional perturbation theory [J]. Reviews of Modern Physics, 2001, 73(2): 515–562. doi: 10.1103/RevModPhys.73.515
[58]	GIANNOZZI P, BARONI S, BONINI N, et al. Quantum ESPRESSO: a modular and open-source software project for quantum simulations of materials [J]. Journal of Physics: Condensed Matter, 2009, 21(39): 395502. doi: 10.1088/0953-8984/21/39/395502
[59]	TSE J S, KLUG D D, YAO Y S, et al. Structure and spectroscopic properties of dense solid hydrogen at 160 GPa [J]. Solid State Communications, 2008, 145(1/2): 5–10. doi: 10.1016/j.ssc.2007.10.018
[60]	MONSERRAT B, NEEDS R J, GREGORYANZ E, et al. Hexagonal structure of phase Ⅲ of solid hydrogen [J]. Physical Review B, 2016, 94(13): 134101. doi: 10.1103/PhysRevB.94.134101
[61]	SINGH R, AZADI S, KÜHNE T D. Anharmonicity and finite-temperature effects on the structure, stability, and vibrational spectrum of phase III of solid molecular hydrogen [J]. Physical Review B, 2014, 90(1): 014110. doi: 10.1103/PhysRevB.90.014110
[62]	MONACELLI L, CASULA M, NAKANO K, et al. Quantum phase diagram of high-pressure hydrogen [J]. Nature Physics, 2023, 19(6): 845–850. doi: 10.1038/s41567-023-01960-5
[63]	BORINAGA M, RIEGO P, LEONARDO A, et al. Anharmonic enhancement of superconductivity in metallic molecular Cmca-4 hydrogen at high pressure: a first-principles study [J]. Journal of Physics: Condensed Matter, 2016, 28(49): 494001. doi: 10.1088/0953-8984/28/49/494001
[64]	WEN L B, WU H, SUN H, et al. Profound softening and shear-induced melting of diamond under extreme conditions: an ab- initio molecular dynamics study [J]. Carbon, 2019, 155: 361–368. doi: 10.1016/j.carbon.2019.08.079
[65]	LI Z G, CHEN Q F, GU Y J, et al. Multishock compression of dense cryogenic hydrogen-helium mixtures up to 60 GPa: validating the equation of state calculated from first principles [J]. Physical Review B, 2018, 98(6): 064101. doi: 10.1103/PhysRevB.98.064101
[66]	RANIERI U, CONWAY L J, DONNELLY M E, et al. Formation and stability of dense methane-hydrogen compounds [J]. Physical Review Letters, 2022, 128(21): 215702. doi: 10.1103/PhysRevLett.128.215702
[67]	SONG X Q, YIN K T, WANG Y C, et al. Exotic hydrogen bonding in compressed ammonia hydrides [J]. The Journal of Physical Chemistry Letters, 2019, 10(11): 2761–2766. doi: 10.1021/acs.jpclett.9b00973
[68]	HEMLEY R J, MAO H K, SHEN G Y, et al. X-ray imaging of stress and strain of diamond, iron, and tungsten at megabar pressures [J]. Science, 1997, 276(5316): 1242–1245. doi: 10.1126/science.276.5316.1242
[69]	MAO H K, BADRO J, SHU J F, et al. Strength, anisotropy, and preferred orientation of solid argon at high pressures [J]. Journal of Physics: Condensed Matter, 2006, 18(25): S963–S968. doi: 10.1088/0953-8984/18/25/S04
[70]	HEMLEY R J, MAO H K. Optical studies of hydrogen above 200 gigapascals: evidence for metallization by band overlap [J]. Science, 1989, 244(4911): 1462–1465. doi: 10.1126/science.244.4911.1462
[71]	NARAYANA C, LUO H, ORLOFF J, et al. Solid hydrogen at 342 GPa: no evidence for an alkali metal [J]. Nature, 1998, 393(6680): 46–49. doi: 10.1038/29949
[72]	LOUBEYRE P, OCCELLI F, LETOULLEC R. Optical studies of solid hydrogen to 320 GPa and evidence for black hydrogen [J]. Nature, 2002, 416(6881): 613–617. doi: 10.1038/416613a
[73]	EREMETS M I, TROYAN I A. Conductive dense hydrogen [J]. Nature Materials, 2011, 10(12): 927–931. doi: 10.1038/nmat3175
[74]	ZHA C S, LIU Z X, HEMLEY R J. Synchrotron infrared measurements of dense hydrogen to 360 GPa [J]. Physical Review Letters, 2012, 108(14): 146402. doi: 10.1103/PhysRevLett.108.146402