Load Characteristics of Underwater Explosion Shock Wave near Seabed Charge Projectile
-
摘要: 近海底爆炸冲击波会对海底光缆、海底管道等设施造成严重的破坏,不同底质条件的冲击阻抗会影响冲击波的时空演化规律,因此,研究不同底质条件下近海底爆炸冲击波载荷具有重要的意义。基于耦合欧拉-拉格朗日方法建立近海底爆炸模型,探究海底底质对近海底爆炸冲击波载荷的影响,结果显示:测点角度在20°~50°范围内时,海底底质材料会显著影响冲击波峰值压力,近海底反射的影响随爆距比的增加而增强,当测点角度超出该范围时,该现象逐渐消失;海底底质影响范围不随底质的变化而变化,但是在不同底质条件下,受影响区域的反射系数截然不同,当海底底质较软时,海底底质影响区域内的冲击波反射系数在0.81~1.05之间,当海底底质为刚壁时,海底底质影响区域内的冲击波反射系数在0.98~1.33之间;水深不会导致冲击波峰值压力发生显著变化。Abstract: The near-seabed explosion shock wave can cause serious damage to such facilities as submarine optical cables and pipelines. The spatiotemporal evolution of the underwater explosion shock wave can be affected by the impedance of different types of substrate. Therefore, it is of great significance to study the near-seabed explosion shock wave load under different substrate conditions. Based on the coupled Eulerian-Lagrangian (CEL) method, a numerical model of near seabed explosion is established to investigate the effect of seabed material on the shock wave load of near seabed explosions. The results show that the seabed material significantly affects the peak pressure of the shock wave within a certain range of measuring point angles of 20°–50°. The reflection effect of the near seabed increases with the increase of the explosion distance ratio within a certain range of measurement point angles of 20°–50°. When the measurement point angle exceeds this range, this phenomenon gradually disappears. The effect regions of these two types of seabed sediment are similar while their reflection coefficients are significantly different. When the seafloor sediment is soft, the reflection coefficient near seafloor ranges from 0.81 to 1.05. However, when the seafloor sediment is hard, the reflection coefficient near seafloor ranges from 0.98 to 1.33. The water depth has little effect on the peak pressure of the shock wave.
-
1911年,荷兰物理学家Onnes[1]首次发现汞在温度降至4.2 K以下时电阻突然降为零,并将这种现象称为超导,将电阻突然降低的温度称为超导转变温度。1933年,Meissner和Ochsenfeld发现,超导体除了具有零电阻的性质外,还会对磁场产生排斥现象,即迈斯纳效应。自超导发现以来,其优秀的性质吸引着国内外众多科研工作者,但是极低的转变温度限制了它的应用,因此寻找具有更高转变温度的超导体成为物理和材料科学的重要科学课题之一。纵观超导材料发展的历史,主要有几个大类:金属单质、合金、铜基超导体、铁基超导体、重费米子超导体及有机超导体等。金属单质与合金作为传统超导体,一般只有在极低温度下才会表现出超导特性,具有较好的力学性能,可以完全排出磁场、零电阻等性质,在电输运、粒子加速器等方面均有广泛的应用。与传统超导体相比,铜基超导体具有更高的超导转变温度,在高压条件下,HgBa2Cam–1CumO2m+2+δ具有高达164 K的超导转变温度[2]。铁基超导体在电子结构、晶格构型等方面与铜基超导体十分接近,并且具有多带超导体的特征,这也是其具有较高超导转变温度的原因之一。重费米子体系是指载流子在低温下表现出极大有效质量的一类固体材料体系。作为典型的强关联电子系统,其超导通常在铁磁、反铁磁、电四极矩等低温有序态边界附近产生,与量子临界涨落有着密切的关系。同样作为非传统超导体的有机超导体大多为低维样品,其超导相图与铜基超导体、铁基超导体、重费米子超导体均十分相似,其独特的超导性质也为我们提供了一条研究低维超导材料中电子与电子以及电子与声子相互作用的途径。对新型超导材料的不断发现和研究,不仅促进了角分辨光电子能谱等实验技术的蓬勃发展,还促进了超导理论和机制的持续探索。
超导现象是宏观量子效应。20世纪初期,爱因斯坦、玻尔、海森堡、费曼等顶级理论物理学家都试图从量子力学的基础上理解超导电性,但最终并没有获得成功。直到1957年,Bardeen、Cooper和Schrieffer提出了解释传统超导体的微观机制,即BCS理论[3]。BCS理论的基本核心是,自旋和动量相反的2个电子通过交换虚声子进而形成库珀对,由于2个电子的动量相反,故当一个电子发生能量损失时,另一个电子得到能量,库珀对总能量不变,即不会产生能量损失。根据BCS理论,超导转变温度Tc=0.85ΘDe−1/N(0)V,其中:ΘD为德拜温度,N(0)为费米面附近的电子态密度,V为有效电声子吸引相互作用。
氢作为质量最轻的物质,具有高德拜温度,因此预测金属氢是一种高温超导体[4]。理论计算表明,氢在接近500 GPa压力下具有356 K的超导转变温度[5]。在常压条件下,氢以气体分子的形态存在,在大约14 K的低温条件下固化为分子相的固体氢,是一种宽带隙绝缘体。高压是实现金属氢的重要途径之一,但是经过80多年的不懈努力,人们仍然无法突破针对固体氢的静高压加载实验极限(500 GPa以上),至今没有观察到原子相金属氢的确切证据,这被称为物质科学领域的“圣杯”。2004年,Ashcroft[6]提出了化学预压缩的概念,通过在氢体系中引入另一元素,降低体系出现金属化和超导转变的压力,富氢化物也因此被认为具有成为高温超导体的潜力。后续对富氢化物的实验研究也进一步验证了这一观点。2014年,共价富氢化物H3S首先被预测在高压下超导转变温度可以达到200 K[7],随后经过众多实验物理学家的探索,发现在155 GPa的压力下H3S的Tc=203 K[8–9]。共价H3S的成功合成极大鼓舞了科研工作者探索高温超导富氢化物的信心。此后,YH6、YH9、YH10、LaH10、CeH10等稀土笼型富氢化物也被预测为潜在的室温超导体[10–11]。Drozdov等[12]通过实验成功合成了Fm
¯3 m-LaH10,并通过原位高压电输运测量确定Fm¯3 m-LaH10在170 GPa压力下具有250 K的超导转变温度,再度刷新了超导转变温度的纪录。在随后对同为镧系元素的铈氢化物的探索中发现,在低于100 GPa条件下,Fm¯3 m-CeH10具有115 K的超导转变温度,与相同对称性的LaH10相比,合成压力降低了超过50 GPa,率先在100 GPa以下发现了超导转变温度超过百开尔文、类金属氢子晶格的超氢化物[13]。随着高压富氢化物的超导转变温度逐渐向室温逼近,国内外众多课题组对其展开了更为深入的探索,一系列高温超导富氢化物被成功合成,如CaH6、YH6、YH9和ThH10等[14–17],如图1所示[1, 2, 8, 12, 18–20]。综上所述,新型富氢化合物特别是以H3S为代表的共价型氢化物和以LaH10为代表的笼型结构都呈现了高温超导的巨大潜力,但是这些新型富氢化物离室温超导还有一步之遥。最近,美国Dias教授团队在实验中分别发现了C-S-H体系和Lu-N-H体系在高压下具有室温超导电性,遗憾的是,这两项工作均受到质疑。提高超导转变温度一直以来都是超导体研究的首要任务之一,在富氢化合物中发现近室温乃至室温超导体仍然是亟待实现的科学目标。1. 主要实验研究方法
1.1 高压实验样品装载
在高压条件下,原子间距离的改变会导致电子能带结构发生变化,致使凝聚态物质形成在常压下无法稳定的亚稳相,从而展现出一系列新奇的物理特性。目前,一系列实现高温超导的氢化物均是在高压条件下合成的。如图2所示[21],研究人员常使用的高压实验装置为金刚石对顶砧(diamond anvil cell,DAC),最高可以承受6000 K高温和超过300 GPa压力,后续改进的金刚石压砧甚至可以承受超过1 TPa的压力。采用激光打孔技术在经过金刚石预压的绝缘金属垫片上面打孔作为样品腔,填入氢源与非氢样品,且保证样品与电极充分接触,以进行后续的样品合成和测量。
1.2 X射线衍射实验:确定晶体结构及化学配比
同步辐射X射线衍射(X-ray diffraction,XRD)是目前唯一用于确定100 GPa以上富氢化合物晶体结构的方法。但是,在超高压下,受金刚石压砧的强烈散射及样品尺寸极小等因素影响,轻质元素原子,特别是具有最小原子量的氢原子的探测非常困难。因此,在富氢化合物晶体结构和化学配比的研究中,氢原子位置的确定是极为重要的一环。目前,氢原子的确定是通过对比实验与理论状态方程确定的。首先,通过同步辐射XRD数据给出非氢元素晶格的相变及对称性,并与理论预测结果进行匹配;然后,根据已有XRD数据所提供的晶胞体积变化与非氢元素的状态方程,对氢原子的含量进行估算,并与理论预测结果进行对比,从而确认氢原子位置;最后,通过理论计算的焓差、晶胞体积随压力的变化以及相稳定性等结果,确定高压富氢化物的晶体结构和化学计量[22]。
1.3 原位高压电学测量
研究氢化物的超导电性时,原位高压电输运测量是必不可少的。在高压条件下,样品一般具有较小的电阻,因此,接触电阻会对样品的电输运测量造成极大的影响。此外,为获得电输运测量信号,需要确保合成的样品与电极充分接触。高温高压合成样品的均匀性以及不规则形状均会对实验结果产生影响。为解决以上困难,通常采用四电极法进行原位高压电学测量,而电极的铺设是其中非常关键的环节。目前,有3种常见的电极铺设方法。第1种是聚焦离子束沉积法,利用镓离子束携带铂离子在金刚石表面沉积电极[23]。该方法的优点是可以控制沉积电极的厚度,附着力比较好,适合小砧面,缺点是价格昂贵。第2种是光刻电极法,如图3所示,首先采用磁控溅射将金属镀在金刚石上,然后利用掩模版在金刚石砧面形成欧姆接触的电极。其优点是成本低,缺点是光刻流程复杂。第3种方法是手工铺设电极,直接将金属电极与样品接触。该方法的优点是操作流程简单,电极连接更牢固,缺点是随着金刚石砧面的减小,操作难度逐渐增大。
1.4 原位高压迈斯纳效应
除了零电阻之外,超导体还具有迈斯纳效应特征,如图4所示。利用DAC装置在超过100 GPa压力下制备的富氢化合物样品尺寸通常小于0.05 mm×0.05 mm×0.01 mm。在这种情况下开展原位高压磁测量,样品产生的磁信号非常微弱,加之来自压机及垫片等方面的干扰,实验上获得的磁信号的信噪比极低。因此,开展超高压下富氢化合物的原位磁化率实验研究极具挑战性,是目前国际公认的难题。目前,已报道的关于富氢化合物迈斯纳效应实验研究的工作屈指可数。例如,Drozdov等[8]使用超导量子干涉仪(SQUID)在155 GPa压力下成功探测到了H3S样品的迈斯纳效应。Huang等[9]利用交流磁化率技术在高压下验证了H3S的超导性。近年来,Eremets课题组[24]对之前H3S的测量技术进行了改进,并对2个高温超导氢化物H3S和LaH10进行了磁化测量,结果表明,这2种氢化物均为第二类超导体。最近,该课题组使用SQUID探测了H3S和LaH10在场冷却和零场冷却条件下样品产生的俘获磁通量行为,对其超导电性进行了进一步验证[25]。
根据法拉第定律,当一个样品处于交变磁场中时,交变磁场只能穿透样品表层,并且穿透深度随着样品电阻的减小而减小。超导体的电阻为零,其穿透深度也为零,即为一个完全抗磁体,相对磁导率为−1。所以,当温度下降至超导体的转变温度时,会观察到交流磁化率的实部将会出现0到−1的突变。目前,高压下样品交流磁化率实部的测量方式主要有2种,分别是互感测量法和低频磁场调制法[26]。其中,互感测量法是通过互感线圈直接对样品的交流磁化率随温度的变化进行测量。然而由于样品体积很小,其信号极易淹没于背底噪声中,因此,通常会在激励线圈中增加一个与感应线圈串联反接的补偿线圈,从而抵消感应线圈中产生的感应电动势,降低噪声带来的影响。低频磁场调制法则是利用当外加磁场大于临界磁场时样品的超导态会被破坏这一原理,使用低频调制磁场来降低背景噪声。
2. 间隙型氢化物
间隙型氢化物中的氢元素一般以原子形式占据金属晶格的间隙位置(如八面体中心位点或四面体中心位点)。Mg、Pd以及LaNi5、Mg2Ni等金属单质及合金可以在一定外界条件下使H2分子解离成原子并插入晶格缝隙中,发生化学反应形成间隙氢化物。由于其储氢与释放氢气的过程比较直观,所以间隙型氢化物也可以作为储氢材料或催化材料。其中,钯作为一种高效储氢材料,其氢化物在极低温度下表现出超导电性。为了探索PdHx的结构,Worsham等[27]通过中子散射技术对钯氢化物进行探索,发现H原子位于钯晶格的八面体位置。目前,钯氢化物的合成方式主要有3种[28],即高压气体加载法、电化学方法以及离子注入法。Schirber等[29]将H2和D2引至Pd样品处,并在室温下加压至0.5 GPa,之后将样品缓慢冷却至77 K,以提高Pd金属中H的含量。通过电阻测量发现,PdH和PdD分别具有8和10 K的超导转变温度。Hemmes等[30]使用DAC将样品加压至4 GPa,合成了类似的产物。Skoskiewiczk[31]则通过电解H2SO4将产生的H2加载至Pd金属中,其使用的电流密度为50~150 mA/cm2。然而,在对样品进行电学测量的过程中,难以避免H2从Pd中溢出,导致样品中H2的分布不均匀,超导转变变宽。Stritzker等[32]在4 K温度下将
H+2 注入Pd中,穿透深度约为1 µm。离子注射方式存在一定的缺点:首先,样品中注入的离子浓度并不均匀;其次,离子注入会导致样品损坏,离子在与样品碰撞的过程中会导致原子位移,从而产生空隙,这些都会导致较宽的超导转变。随着实验技术的发展,Tripodi等[33]发现了钯氢化物新的超导相,具有51.6 K的超导转变温度,在77 K的温度下具有大于6.1×104 A/cm2的临界电流密度。PdHx还具有奇异的负同位素效应,当将H替换成更重的D后,样品转变温度提升。这种奇异的性质吸引了众多研究者的注意[34–38],对这一现象的解释主要归因于氢原子的非谐性。金属Pt和Pd具有相同的电子能带结构,理论预测发现,Pt和H2在高于21.5 GPa的压力下会合成PtH,并在70~80 GPa时转变为六方密排(hexagonal close-packed,hcp)或面心立方(face-centered cubic,fcc)结构,其中面心立方结构的PtH具有最高25 K的转变温度[39]。后续的铂氢化物实验研究中,Scheler等[40]利用高压技术将Pt单质和氢气加压至27 GPa,反应合成PtH-Ⅰ(P321)和PtH-Ⅱ(P63/mmc),并且PtH-Ⅰ和PtH-Ⅱ在42 GPa下可以共存。当压力进一步升高时,合成的样品中只存在PtH-Ⅱ。PtH-Ⅱ相属于六角密堆积结构,并且理论预测PtH-Ⅱ在90 GPa时的超导转变温度为12 K。为继续验证PtH体系的超导电性,Matsuoka等[41]利用激光加热技术,将Pt箔与氨硼烷在30 GPa下加热至1000 K合成了PtH-Ⅱ,在35 GPa下测量了合成样品的电阻与温度的关系,发现样品在7 K时发生了超导。当外加磁场为0.1 T时,电阻下降斜率突变,超导转变温度变为3.5 K;外加磁场增加到0.5 T时,电阻的下降趋势在2.8 K消失。由于超导转变温度Tc会随着外部磁场的增加而向低温转移,说明PtH-Ⅱ中存在超导性质。
3. 离子型氢化物
Ashcroft认为在纯氢体系中引入非氢元素可以实现化学预压作用,从而降低实现超导的压力。由于引入非氢元素后化学性质不同,且非氢元素与氢形成的价键也有所区别,一般来说根据氢与非氢元素的成键可以将这些氢化物分为共价型氢化物和离子型氢化物。离子型氢化物中,非氢元素与氢之间成离子键,存在电荷转移。与共价型氢化物相比,离子型氢化物的种类更加繁多,为探索高温超导富氢化物提供了更多选择。对于离子型富氢化物的超导电性,部分研究者将目光聚集在第Ⅳ主族元素(C、Si、Ge、Sn),并在实验上成功观测到了Si-H和Sn-H体系的超导现象。首先,通过第一性原理计算预测了SiH4在高压下可能存在的结构和稳定性。在0~40 GPa压力范围内,SiH4会分解为单质硅和氢;在40 GPa以上,SiH4具有I41/a和C2/c结构。并且计算了SiH4由绝缘态转变到金属态的压力范围为20~90 GPa,其超导转变温度在40~260 K之间[42–45]。Feng等[46]预测SiH4有13种不同结构,他们认为:在低压下,SiH4中的Si原子具有四面体配位时结构最稳定;在25 GPa 以上,Si原子具有六重配位和八重配位时的结构是稳定的。此外,Hanfland等[47]将单质硅和氢气直接加压到124 GPa,合成了2种结构不同的SiH4,在90~124 GPa时,SiH4的空间群为I41/a和I
¯4 2d。此外,在60~90 GPa的压力区间内,SiH4出现了非晶态现象,大于90 GPa时,SiH4重结晶,由非晶体结构转变为聚合物晶体结构。2008年,Eremets等[48]首次实验研究了SiH4的金属化和超导特性。他们发现:SiH4在50 GPa下能够发生绝缘态向金属态的转变,并且在压力为50 GPa时样品表面突然变暗;进一步加压到76 GPa时,样品表面完全不透明。如图5(c)所示,通过超导转变温度随压力的变化关系发现:在96 和120 GPa下,SiH4的超导转变温度约为17 K;在165 GPa下,超导转变温度降低到8.8 K;然而加压到192 GPa时,超导转变温度增加到11.3 K。
在研究了SiH4的电学性质后,Eremets等[48]进一步研究了SiH4的结构。SiH4的X射线衍射峰与六角密堆积结构匹配,并且SiH4的晶格参数c/a为1.679,理想六角密堆积结构的晶格参数c/a值为1.633,两者较为接近。进一步研究发现,在50~110 GPa压力范围内,SiH4的空间群为P63。当压力超过120 GPa时,SiH4的结构由P63转变为透明相。实验进一步确定了该透明相的结构,该透明相的结构与体心立方结构吻合,对应的晶格参数为:a=b=3.03 Å,c=6.87 Å,空间群为I41/a。理论预测表明,I41/a结构在热力学和动力学上可以稳定存在的压力范围为50~250 GPa,并且该结构是绝缘的[44]。图5(c)中,随着压力增加到100 GPa,超导转变温度逐渐升高;随着压力进一步增加,超导转变温度降低,说明SiH4由初始分子结构转变为金属相P63;当压力超过120 GPa时,SiH4中既存在金属相P63也有绝缘相I41/a,这也很好地解释了SiH4的超导转变温度随压力的变化趋势。
锡氢化物的超导研究最早可以追溯到2007年,Tse等[49]利用第一性原理预测了SnH4在高压下具有金属相,并且稳定的压力范围为70~160 GPa,在120 GPa下的超导转变温度为80 K。最近,Hong等[50]在实验上成功检测到了锡氢化物的超导态。他们将锡箔和氨硼烷加压到200 GPa,然后将两者加热至1700 K,合成了锡氢化物。研究发现,在200 GPa下,70 K时锡氢化物发生转变,说明该样品可能存在超导电性。为了进一步确认合成的样品具有超导性,在不同的外加磁场下进行了电阻-温度(R-T)测量。随着外加磁场增加,电阻转变位置也逐渐向低温区移动,符合超导态的性质。根据GL方程拟合,确认临界磁场μ0Hc2(0)=11.2 T,超导转变温度为67.3 K。最近,Troyan等[51]通过原位高压电输运测量,确认fcc-SnH4在180 GPa下的超导转变温度为72 K,并通过0~16 T外加磁场下的R-T曲线确认其临界磁场为16 T。
4. 共价型氢化物
如前文所述,如果起预压作用的非氢元素与氢之间成共价键,那么此类氢化物就称为共价型氢化物。金属化的共价键会导致高电声耦合,同时低原子质量的非氢元素会引入高频声子模,对电声耦合有一定的帮助。作为典型共价氢化物的硫氢体系便是其中广为人知的氢化物体系,其全氢原子的三维硫氢共价网络结构是具有高温超导电性的关键原因。
硫氢化物的高压研究最早始于1991年。通过拉曼光谱、红外光谱和同步辐射实验揭示了H2S在46 GPa下发生分解,在96 GPa下具有金属性质[52]。然而,固体H2S的高压相图还需要进一步完善。因此,Li等[53]基于粒子群优化算法研究了H2S在10~200 GPa内可能存在的结构。此外,通过电声耦合计算发现,H2S有可能是一种高温超导体。基于该理论计算,Drozdov等[8]将气态H2S在200 K低温的恒温器中液化为液态H2S。当加压到50 GPa时,样品开始导电,具有半导体性质。进一步加压到90~100 GPa,样品由半导体向金属转变。此外,如图6(a)所示,当H2S在100 K低温下压缩时,测到的超导转变温度为170 K;而在200 K下压缩后退火,可以测到一个相对较高的超导转变温度,即195 K。如图6(b)所示,通过硫化氢和硫化氘的电阻随温度的变化发现,硫化氢和硫化氘的Tc分别为195和150 K,进一步确认硫化氢在高压下存在超导态。
图 6 (a) 恒压退火过程中电阻与温度的关系,(b) 高压下硫化氢和硫化氘电阻与温度的关系,(c) H3S的超导转变温度随压力的变化[8]Figure 6. (a) Dependence of resistance to temperature in constant pressure annealing process; (b) critical temperature of sulfur hydride and sulfur deuteride at high pressure; (c) dependence of superconducting temperature with pressure for H3S[8]超导体除了具有零电阻特性,另一个重要的特征是完全抗磁性,即迈斯纳效应[24]。在高压条件下测量富氢化合物的抗磁性是十分困难的,主要原因是合成的样品尺寸太小,导致样品产生的磁性信号非常弱,很难获取有用的信号。Huang等[9]通过抑制信号噪声源、屏蔽传输噪声等方法成功测得硫化氢的迈斯纳效应,确定了硫化氢的超导性。如图7所示,通过测定的超导相图发现,硫化氢在117 GPa下开始出现超导性质,超导转变温度仅为38 K,在149 GPa下超导转变温度可以达到183 K。当压力为171 GPa时,超导转变温度降低到140 K。硫化氢的超导转变温度呈现随压力先升高后降低的趋势,理论计算分析表明,该变化趋势主要是由于硫化氢中化学计量变化引起的,并且与Drozdov等[8]的超导转变温度的变化趋势一致,从而进一步证实硫化氢确实属于高温超导体。
虽然实验证实了硫氢化物确实存在超导现象,但是具有高温超导的硫氢化物结构并没有确定。Li等[53]预测了H2S在160 GPa下的超导转变温度为80 K,但是Drozdov等[8]实验中的Tc约为200 K。这说明实验中具有高温超导的硫氢化物可能不是H2S。Fujihisa等[54]通过实验确认了H2S在43 GPa下会分解,其分解产物中包含单质硫。但是,硫的超导转变温度低于25 K,因此,可以排除单质硫对硫氢化物高温超导的影响。在此基础上,推测H2S在高压下可能产生了其他硫化物,即H3S、H4S和H6S。Duan等[7]预测了(H2S)2H2在180 GPa下会形成Im
¯3 m相H3S,Tc接近190 K。因此,H2S可能在高压下分解生成H3S,这才是样品高温超导的来源。随后的实验进一步证实了H3S的高温超导电性。Einaga等[55]按照Drozdov实验中的合成方法进行实验,除了测出了200 K的高温超导转变温度,还对合成的样品进行了同步辐射XRD分析。他们认为,实现高温超导的样品是Im
¯3 m相H3S和单质S。此外,Guigue等[56]通过激光加热固体硫和氢气,在75 GPa下合成了H3S,并且证实了H3S的结构Cccm相,Cccm-H3S可以稳定到160 GPa。与Guigue等[56]的研究结果不同,Goncharov等[57]也成功合成了H3S,根据XRD图谱,在50、70和140 GPa下对应的空间群分别为Cccm、R3m和Im¯3 m。Capitani等[58]根据原位高压电输运测量与红外光谱实验判断150 GPa时H3S的超导转变温度为200 K。此外,Troyan等[59]将放置在H2S中的Sn薄膜作为磁场传感器,通过同步辐射核共振散射检测Sn传感器的磁场,并确认在153 GPa压力下具有140 K的超导转变温度。因此,利用低温加载H2S实现高温超导的产物极有可能是H3S。表1总结了不同实验组确认H3S为高温超导体的研究情况。Groups Sample synthesis Tc/K Corresponding
pressure/GPaStructure of
superconducting phasesDrozdov et al.[8] Loading H2S gas at low temperature 203 155 Undetermined Huang et al.[9] Loading H2S gas at low temperature 183 149 Calculation: Im¯3m-H3S Einaga et al.[55] Loading H2S gas at low temperature 200 150 XRD: Im¯3m-H3S Goncharov et al.[57] Laser heated pure S and hydrogen Unmeasured Unmeasured XRD: Cccm-H3S@50 GPa,
R3m-H3S@70 GPa,
Im¯3m-H3S@140 GPaCapitani et al.[58] Loading H2S gas at low temperature 200 150 Undetermined Troyan et al.[59] Loading H2S gas at low temperature 140 153 Undetermined 钍氢化物Th4H15中的钍原子与氢原子之间成弱共价键,也可以认为是共价氢化物。常压下稳定的Th4H15具有7.5~8.5 K的超导转变温度,是具有I
¯4 3d空间群对称性的立方结构[60]。Shein等[61]通过第一性原理计算得到Th原子的5f、6d电子对Th4H15中费米面附近电子态密度起主要作用,其中5f电子贡献最大,估计其平均电声耦合常数λ为0.38,Th4H15可以视为弱耦合超导体。但后来Talantsev等[62]进一步拟合Wang等[63]的实验数据后,得出当压力从0 GPa升高至8 GPa时电声耦合常数λ在0.81~0.88之间波动,近似于恒定,并得出Th4H15为中等耦合超导体的结论。Wang等[63]将NH3BH3作为氢源,在5 GPa、1100 K的条件下合成了Th4H15,电阻测量得到其Tc为6 K,比之前的报道低2 K,同时经过磁学测量与理论计算得出Th4H15为Ⅱ型超导体。磁学和热容测量结果表明,Th4H15缺少同位素效应[64–65],这对传统电声耦合超导机理提出了质疑。部分研究认为,对应于氢的光学声子振动模式是超导出现的可能原因[63, 66]。很明显,如果与氢晶格有关的光学声子模的作用体现在超导上,那么同位素效应是不可避免的,所以关于Th4H15的超导机理尚待深入研究。高压下Th4H15的超导转变温度变化趋势为先升高后降低:Wang等[63]测得的变化趋势为当压力从0.7 GPa增大到5 GPa时,Tc从6.9 K上升至7.5 K;Dietrich等[67]测得当压力从常压增大到2.8 GPa时,Tc以0.3 K/GPa的速率近线性上升。值得一提的是,完全超导温度Tczero随化学配比变化较为明显,但超导开始转变的温度Tconset基本不随化学配比变化[62];另一方面,随着压力上升,超导转变区间有拓宽现象[63]。这些都会导致超导转变温度选取的不准确性,进而影响最终结论。与Th4H15类似的结构还有Hf4H15和Zr4H15。Kuzovnikov等[68]利用HfH2与氢气反应得到Th4H15型的新型铪氢化物Hf4H15,并在23 GPa下测得其Tc为4.5 K。Xie等[69]将样品I4/mmm-ZrH2和NH3BH3在DAC中压缩到30 GPa,然后激光加热到1600 K,第一次在ⅣB族氢化物中成功合成了Pm
¯3 n-ZrH3和I¯4 3d-Zr4H15。相比于其他Zr-H化合物,ZrH3和Zr4H15的氢含量较高,质量分数分别为3.2%和4.0%。理论计算和电学实验表明,ZrH3和Zr4H15均为常规超导材料,如图8(a)所示,在40 GPa下的超导转变温度分别为6.4和4.0 K;但在理论计算中,ZrH3在40 GPa下的Tc为8.4~12.4 K,Zr4H15在40 GPa下的Tc为0.2~0.8 K。实验所得Tc异于理论计算,其原因可能是非简谐效应。此外,从理论计算可以看出,ZrH3与之前发现的高温超导体H3S[53]相比,虽然费米能级N(εf)更大,但H的贡献要小得多,这不利于电声耦合,因而H3S的Tc 是ZrH3的19~25倍。图 8 (a) 40 GPa下Zr-H样品的电阻-温度变化曲线(插图为激光加热后的样品腔),(b) 40 GPa下超导转变温度随外加磁场的变化,(c)~(d) 分别用WHH与GL方程外推拟合上临界磁场[69]Figure 8. (a) R-T curve for Zr-H sample at 40 GPa (The sample chamber is presented in inset); (b) superconductingtemperature under applied magnetic fields at 40 GPa; (c)–(d) upper critical magnetic fieldswhich extrapolated by GL and WHH equations, respectively[69]最近,Zhang等[70]在240 GPa、2000 K的温压条件下成功合成了具有更高氢含量的锆氢化物,其超导转变温度在220 GPa时达到71 K。理论预测Cmc21-ZrH6在150~275 GPa压力下是稳定的,在215 GPa时的Tc约为70 K[71],所以预测Zhang等[70]观测到的锆氢化物超导性可能来自于Cmc21-ZrH6。
5. 笼型氢化物
CaH6、LaH10、CeH9、CeH10、YH6、YH9和ThH10等氢化物[12–17, 72–73]均具有典型的氢笼结构,见表2。非氢元素提供电子帮助H2解离形成稳定的笼型结构,故笼型氢化物也是离子型氢化物的一个分支。氢笼结构在费米面附近提供了巨大的电子态密度以及强电声耦合效应,所以这些笼型氢化物均具有很高的超导转变温度,是潜在的室温超导体。
Sample Space group Pressure/GPa Tc/K Mensurement method LaH10 Fm¯3m 170 250 XRD ThH10 Fm¯3m 170 161 XRD YH6 Im¯3m 183 220 XRD YH9 P63/mmc 201 243 XRD CeH9 P63/mmc 130 100 XRD CeH10 Fm¯3m 95 115 XRD CaH6 Im¯3m 172 215 XRD 如图9(a)所示,最先实验合成的笼型高温超导氢化物为立方相LaH10,具有H32笼结构,是目前实验合成的具有最高超导转变温度的二元氢化物之一,在实验合成之前就已经有众多理论预测LaH10在高压下的超高转变温度[10–11],吸引了众多课题组对其进行探索。Somayazulu等[73]首先将氨硼烷作为氢源将样品的温压升至175 GPa、1000~1500 K,成功合成了立方相LaH10。通过电阻测量发现,样品在180~200 GPa压力下具有260 K的超导转变温度,如图9(b)所示,并且超导转变温度随着测量电流的增加而降低。但由于实验过程中部分样品延展到金刚石砧面外,并且在激光加热时不可避免地产生温度梯度,导致合成样品为混相。随后,Drozdov等[12]以H2为氢源在150 GPa压力下对样品激光加热至700 K合成了Fm
¯3 m-LaH10,同样的产物也可以通过将镧金属加压至141 GPa并暴露在H2中2周合成。通过电阻测量发现,样品在151 GPa具有249 K的超导转变温度,如图9(c)所示,外加磁场时发现,随着磁场增强,转变温度逐渐降低。通过GL方程拟合得到临界磁场为136 相干长度处于1.56~1.86 nm之间。在对样品进行加卸压电阻测量后发现,样品的转变温度随压力增加先升高后降低,并在170 GPa时达到最高,为250~252 K。随后进行同位素效应研究,采用相似的实验方法得到立方相LaD10,在152 GPa压力下具有180 K的转变温度。除此以外,Hong等[74]将La与氨硼烷加压至165 GPa,并分别加热至1700和1000 K,合成了转变温度分别为250和74 K的镧氢化物。Sun等[75]同样在高温高压条件下合成了LaH10,其136 GPa下的超导转变温度为246 K,当卸压至120 GPa时,LaH10的结构从Fm¯3 m转变为C2/m相,转变温度降低至191 K。图 9 (a)立方相LaH10的氢笼结构[78],(b) Somayazulu等[73]合成的LaH10的电阻随温度的变化,(c) Drozdov等[12]合成的LaH10的电阻随压力的变化,(d) Huang等[76]外加磁场作用下LaH10的转变温度Figure 9. (a) Clathrate structure of LaH10[78]; (b) the superconducting critical temperature of LaH10 reported by Somayazulu et al.[73];(c) the dependence of Tc with pressure for LaH10 synthesized by Drozdov et al.[12]; (d) the critical temperature forLaH10 under applied magnetic fields reported by Huang et al.[76]近年来,本课题组同样对镧超氢化物的超导电性进行了探索[76]。将氨硼烷作为氢源,与La样品一起加压至171 GPa,并通过激光加热至1400 K。使用四电极法对样品电阻进行测量,结果如图9(d)所示,在238.4 K下观察到超导转变。为了进一步验证样品的超导电性,对其施加0~6 T的外加磁场,并观察到随着磁场强度的增大,超导转变温度减小。通过GL方程和WHH方程拟合得到其临界磁场分别为114.9和158.1 T,与Drozdov等[12]之前报道的实验数据吻合。随后,使用公式
μ 0Hc2=Φ0/(2πξ 2)计算相干长度ξ WHH(0)和ξ GL(0),分别为1.44和1.69 nm,与Shukor[77]报道的结果15.6~18.6 Å相符。Kvashnin等[79]预测锕系元素中的ThH10具有与LaH10类似的氢笼结构,并在100 GPa压力下具有241 K的转变温度。Semenok等[17]对钍氢化物进行了实验研究,通过加热钍单质和氨硼烷,在170~175 GPa压力下合成了Tc分别为146和161 K的P63/mmc-ThH9和Fm
¯3 m-ThH10。后续对样品施加0~16 T的磁场,拟合得出ThH9和ThH10的上临界磁场分别为38和45 T。对ThH10进行卸压研究发现,ThH10可以稳定至85 GPa,与理论预测结果十分吻合。对于钇氢化物,理论预测了一系列笼型钇氢化物YHx(x=4, 6, 9, 10)。其中,YH6、YH9以及YH10的超导转变温度超过200 K,是潜在的高温超导体[10–11, 80],吸引了众多实验研究者的兴趣,Shao等[81]以氨硼烷为氢源,在155 GPa压力下通过加热YH2样品至1500 K成功合成了I4/mmm-YH4。其超导转变温度在155 GPa压力下达到88 K,与之前Li等[80]理论预测的结果相符。在卸压过程中发现,I4/mmm-YH4可以稳定至113 GPa,但是在更低的压力下金刚石发生损坏。Troyan等[14]则将Y单质和NH3BH3样品加压至160 GPa以上,通过激光加热至2400 K,成功合成了具有氢笼结构的Im
¯3 m-YH6。在166 GPa压力下其Tc为224 K(见表3),通过GL方程和WHH方程拟合得到其临界磁场分别为116和158 T。之后,Griessen等[82]通过公式Jc=Jc0(1−T/Tc)5/2(1+T/Tc)−1/2对样品的临界电流密度进行计算,发现在4 K时样品的临界电流和临界电流密度分别可以达到1.75 A和3500 A/mm2。此外,还有2个课题组分别独立合成了更高H配比的P63/mmc-YH9[15, 72]。其中,Kong等[15]通过将YH3和H2样品加压至160 GPa以上并激光加热至1500 K,成功合成了体心四方(bct)YH4和体心立方(bcc)YH6,YH6在183 GPa具有220 K的超导转变温度。这2种产物除了可以通过激光加热制备外,也可以通过将YH3置于氢气中加压至200~244 GPa并放置2周得到。具有更高H配比的六方相YH9则是在185 GPa压力下对之前得到的YH4和YH6样品激光加热至2 000 K获得。YH9在201 GPa具有243 K的超导转变温度(见表3)。对样品进行加卸压电阻测量发现,当压力减小至159 GPa时,YH9分解为YH6,而YH6则稳定至135 GPa。Snider等[72]则使用Y单质作为反应物,以氨硼烷为氢源,在超过130 GPa的压力下将样品加热至1800 K,得到最高超导转变温度为262 K的钇氢化物。值得注意的是,Snider等[72]利用金属钯促进氢吸收,通过在样品钇表面镀约10 nm厚的钯层合成更高H配比的钇氢化物。与Kong等[15]的实验结果对比发现,Snider等[72]得到的样品在更低的压力下具有更高的超导转变温度。但是,由于缺乏XRD数据,无法确定Snider等[72]合成的样品为二元钇氢化物,钯金属的引入有可能使其合成了一种新型三元钯钇氢化物。表 3 钇超氢化物合成条件及超导转变温度Table 3. Superconducting critical temperature of yttrium superhydridesSamples Ref. Space
groupsPressures/
GPaTc/K Reactants Synthetic methods YH6 Troyan
et al.[14]Im¯3m 166 224 Y+NH3BH3 Laser heated to 2400 K at high pressure Kong
et al.[15]Im¯3m 237 208.5 YH3+H2 Kept the sample for three weeks under high pressure 183 220 Y+H2 Laser heated to 1500 K at high pressure 159 220 YH3+H2 Increased the sample to 201 GPa for one month
and laser heated the sample to 2000 KYH9 Snider
et al.[72]P63/mmc 182±8 262 Y+H2 Pressured the Y metal wrapped by Pd film and hydrogen to
over 130 GPa and laser heated the sample to 1800 KKong
et al.[15]P63/mmc 237 227 YH3+H2 Maintained the sample for three weeks 237 237 YH3+H2 Maintained the sample for three weeks and
laser heated sample to 700 K201 243 YH3+H2 Maintained the sample for one month
and laser heated sample to 2000 KYH4 Shao
et al.[81]I4/mmm 155 88 YH2+NH3BH3 Increased the pressure to about 150 GPa
and laser heated sample to 1500 K铈氢化物CeH9、CeH10同为镧系元素氢化物,其理论预测超导转变温度远低于LaH10,均小于56 K[11]。在另一篇理论预测报道中,F
¯4 3m-CeH9和Fm¯3 m-CeH10在94 GPa压力下分别具有142和168 K的超导转变温度[83]。本课题组[22, 84]在不同实验条件下成功合成了一系列铈氢化物,通过Ce金属和H2冷压缩至159 GPa直接反应合成了P63/mmc-CeH9,具有与立方相LaH10类似的氢笼结构,并且CeH9中氢原子间距是目前所合成的氢化物中最接近原子金属氢的。电子局域函数(electron local function,ELF)表明,铈原子与氢原子之间成离子键,电子态密度(density of electronic states,DOS)计算显示了费米面附近H的明显贡献。Salke等[85]则在80~100 GPa压力下,通过将CeH3样品与H2加热至2000 K成功合成了CeH9。通过对样品进行加卸压测量发现,CeH9可以稳定至93 GPa,随着继续卸压至50 GPa,Pm¯3 n-CeH3重新出现,完全卸压后的产物为CeH3和Ce2H5。通过理论计算得到P63/mmc-CeH9在200 GPa下具有105~117 K的超导转变温度。为了继续探索铈氢化物的超导电性,如图10(a)和图10(b)所示,Chen等[13]在高温高压条件下合成了P63/mmc-CeH9和Fm¯3 m-CeH10,并使用四电极法对其电阻进行测量。其中,CeH10在95 GPa下具有115 K的转变温度(图10(c)),而CeH9在88 GPa下具有57 K的超导转变温度,并随着加压至130 GPa,Tc提高到100 K,而当压力继续增大后,Tc逐渐减小。通过理论计算发现,这可能是由于CeH9由P63/mmc相转变到C2/c相所致。实验观察到的Tc随着外加磁场增强而减小以及同位素效应均验证了样品的超导电性。Wang等[86]在2012年最先通过理论计算预测得到立方相CaH6的笼型氢化物。在200 GPa压力下,CaH6为立方Im
¯3 m相,其结构表现为由6个H4单元相互连接组成的方钠石框架包裹位于中心的Ca原子的氢笼结构。ELF结果表明,Ca原子与H原子之间不成键,而H原子之间成弱共价键。150 GPa时其Tc为235 K,随着压力的增加,Tc逐渐减小,在200 GPa时降低至201 K,250 GPa时Tc为187 K。但是,在CaH6提出后的接近10年里,并没有成功通过实验将其合成出来。最近,2个课题组在超过150 GPa高压下通过激光加热方式成功合成出了具有超过200 K超导转变温度的钙氢化物。Li等[87]使用氨硼烷(AB)作为氢源,在160~180 GPa压力下将样品加热至1300 K合成了超导转变温度为210 K的钙氢化物,其Tc与理论预测CaH6的Tc相符。并且在电阻-温度曲线上可以观察到,除了第1个超导转变外,在160和180 K等温度处也出现了几个台阶(图11(a)),这可以归因于其他不同配比的钙氢化物的相。通过增大样品的外加磁场并进行电阻测量,观察到样品的超导转变温度逐渐降低,通过GL方程拟合得出样品的临界磁场为166 T。与此同时,另一课题组在160~190 GPa压力下同样以氨硼烷为氢源将样品加热至2000 K,成功合成出转变温度为215 K的钙氢化物(图11(b))[16],通过同步辐射XRD测量确定该样品为Im
¯3 m-CaH6。为了验证之前理论预测CaH6在150 GPa时的超导转变温度为220~235 K,对样品进行加卸压电阻测量,并在172 GPa压力下观察到样品的最高Tc为215 K,与之前理论预测结果相符。在样品Tc对压力依赖性的实验中发现,不同卸压实验得到的结果存在差异,其原因可能是卸压过程中存在不同程度的各向异性应力导致晶格变形。在172 GPa压力下对样品施加0~9 T的外加磁场发现,随着磁场强度的增大,超导转变温度逐渐降低。通过GL方程和WHH方程拟合得到临界磁场分别为142和203 T,相干长度分别为15.2和12.7 Å,远低于Ⅰ类超导体数百埃的相干长度,说明CaH6为一种Ⅱ类超导体。6. 分子型氢化物
有些氢化物中存在H2或H3单元,称为分子型氢化物。但一般来说,H2或H3单元往往对超导有负面作用,电子局域在氢分子单元里导致氢原子在费米面附近的电子态密度上很难有较大贡献,因此,即使其具有较高的氢含量,但仍然不具备良好的超导性质。
在元素周期表中,Ba和Ca同主族,和La同周期。目前,已经实验合成了Tc分别为215 和260 K的CaH6和LaH10[16, 73],因此,钡氢化物也有望成为高温超导体。Chen等[88]将单质钡和氨硼烷在高压下激光加热合成了BaH12。实验发现:在160 GPa下,BaH12的结构为Cmc21相;在84 GPa下,BaH12的结构为Fm
¯3 m相;BaH12在75~173 GPa压力范围内稳定存在。如图12所示,电学测试表明,合成的BaH12在140 GPa压力下的超导转变温度可以达到20 K。此外,也有理论和实验研究了其他钡氢化物的结构和性质[89–92]。Luo等[89]预测4 GPa下BaH2由Pnma结构相变为简单六角结构。Hooper等[90]预测P4/mmm-BaH6在100 GPa下的超导转变温度为38 K。Tse等[91]利用同步辐射XRD确定了BaH2在45 GPa压力下由hcp结构相变为简单六角结构,预测BaH2的Tc接近0 K。Zhang等[92]则测量了不同相时电输运性质的变化。Peña-Alvarez等[93]利用Ba或BaH2与氢气在50 GPa、1200 K下合成了新型Ba8H46,并通过同步辐射XRD确定了Ba8H46的结构为Pm
¯3 n,其结构特点是一个Ba原子由20或24个H原子包围形成笼型结构。但是,尚不清楚实验合成的Ba8H46的超导性质。最近,Semenok等[94]通过DAC将氨硼烷合成了一系列锶氢化物:C2/m相Sr3H13、伪立方P1-SrH6、伪四方P1-SrH22等。其中,SrH22是目前含氢量最高的氢化物之一,可以稳定到100 GPa。理论研究表明,非简谐效应使SrH22在100 GPa稳定存在,且150 GPa时SrH22为不含有序氢亚晶格的超离子氢化物,约200 GPa时可以通过能带重叠实现金属化。7. 待解决的科学问题
7.1 样品同一性和均匀性的问题
在制备氢化物样品的过程中,除让其自行反应外,更多时候需要进行各种操作来制备样品,这些操作往往是在原始样品的某个局部进行的,这就会带来样品不均匀的问题。比如,激光加热时,激光斑点正中心处可能产生高含氢量的氢化物,边缘处可能产生低含氢量的氢化物,而远离激光斑点的其他部分可能由于温度不够高或者氢含量过少而不能生成氢化物。混相的存在会给各种高压原位测试带来一定的干扰,在电学测量中混相可能会带来阶梯状的超导转变温度,XRD实验中会让峰形变得复杂难以分析,所以如何制备出均匀的优良样品是一个难题。
由于实验条件不同,往往需要用到不同种类的压机进行原位高压测试,即便制备条件相同,这些压机中所制备的样品是否会是同一种样品也很难完全确定,如果出现误差,那晶体结构与超导转变温度很可能无法反映同一个样品的性质。
受样品纯度以及实验手段的限制,很难避免在测量样品中引入杂质、晶体结构缺陷等,这些问题会在超导体中束缚住一定的磁通线,就像钉子一样把磁通线固定在晶体内部,这种现象叫做钉扎效应。考虑到束缚的磁通与外加磁场方向相同或相反,钉扎效应会导致临界磁场的测量结果偏低或偏高,这也是一个难以解决的问题。
7.2 超导类型的确认
常规超导体是属于各向同性的S波超导体。1950年,Maxwell发现了同位素效应,即随着同位素质量M的增加,Tc下降,这在H3S、CeH9和LaH10体系中得到了实验验证。目前,实验上合成的氢化物超导体基本上是符合BCS理论的常规超导体,并且根据BCS理论预测的H3S、LaH10和CaH6等高温超导体的超导转变温度和超导相可以很好地与实验结果自洽。但是,对于常规超导体而言,其超导转变宽度会随着磁场强度的增大而显著增大,如图13所示,但在这些典型高温超导氢化物并没有明显表现出这种现象。
7.3 H原子位置的直接测量
氢对于形成高压氢化物的高超导转变温度起到至关重要的作用,然而对于氢化物中氢原子的含量、原子位置等信息的研究十分困难。由于样品在高压下发生金属化,导致无法使用拉曼光谱和红外光谱对样品进行表征,并且氢原子只有一个电子,导致其XRD谱基本是不可见的。而绝大多数高温超导氢化物都是在超过百吉帕压力下合成的,这就导致中子衍射受到限制。因此,只能通过计算晶胞体积与更重的母金属原子的体积差来估算氢化物中的氢含量。但是,这种近似的计算方法存在一定的限制。首先,大多数合成的样品均为混相,这就对XRD数据的分析造成了极大的困难;其次,对样品的XRD粉末衍射图案进行分析时,需要与理论计算结果进行对比,但是很多时候只能从类似的体系中得到信息。因此,寻找一种可以准确表征高压氢化物中氢原子信息的实验手段也是一个亟待解决的问题。
7.4 迈斯纳效应的证据
对于超导体而言,除了具有零电阻的特征以外还具有迈斯纳效应,即当物体处于超导态时,其磁通量被完全排出体外。但是由于高压下样品体积极小,对高压超导氢化物迈斯纳效应的测量十分困难。到目前为止,只对H3S、LaH10等寥寥几种超氢化物进行了迈斯纳效应的测量[8–9, 24],更多超导氢化物的迈斯纳效应等待我们验证。
近日,Bhattacharyya等[95]在高温高压条件下合成了铈氢化物高温超导体,并通过NV色心量子传感技术对其局部磁通变化以及迈斯纳效应进行表征。NV色心是金刚石内部的发光点缺陷,对磁场变化具有极高的灵敏度。在本次实验中,研究团队通过测量场冷和零场冷条件下样品不同位置的抗磁性,得到与电输运性质相匹配的结果。这项研究成功验证了高达百吉帕压力条件下极微量超导富氢化物样品的迈斯纳效应,破除了对富氢化物超导体存在的质疑,极大鼓舞了科研工作者对探索富氢化物高温超导体的信心。
7.5 降压制备
氢化物超导体一般只能在高压下保持稳定,在常压下处于亚稳态,所以如何在高压制备后使其保留到常温一直是困扰氢化物向应用方向更进一步的难题。目前,通过引入另一种新元素可以使氢化物超导体合成压力有一定的降低,但其超导转变温度却同时受到影响。我们期待可以有更多的物理或化学方面的实验手段使氢化物可以在更低压力下维持超导性质。
8. 结论与展望
8.1 高温超导体的一些普适性规律
独特的结构和成键类型影响着富氢化合物超导转变温度的高低。氢的排列与非氢元素和氢之间的电荷转移有关。共价型氢化物和笼型氢化物都需要在非氢元素与氢之间转移更多的电子,以形成有利超导的结构。但是,层状化合物与共价氢化物及笼型氢化物不同,它们不会破坏固态氢的高压相,并且允许其在较低压力下出现。与上述3种氢化物不同,如BaH12 氢化物中含有孤立的氢原子,H2和H3分子单元氢化物的超导转变温度一般不会超过200 K,甚至低于100 K。在这类氢化物中,H与H之间是共价作用,并且在相同压力下,其H―H键的距离比纯氢中H―H的键长更长一些。这类氢化物超导转变温度低的原因是H2和H3分子单元对电声耦合强度贡献很小。
对于具有高超导转变温度的氢化物,总结有以下4个特点:(1) 高超导转变温度的氢化物倾向于形成对称性更高的结构,如现在广为人知的笼型结构就很好地证明了这点;(2) 在高超导转变温度的氢化物中,氢元素倾向以原子形式存在,而自然界中比较常见的H2、H3分子单元则一般对于超导转变温度的提高有不利的影响;(3) 费米面附近的电子态密度一般由氢原子提供或氢原子在其中占有较大的贡献;(4) 在晶体结构中存在较强的电声耦合[96]。因此,为了获得高温超导氢化物,首先要消除氢化物的H2分子单元,可以引入额外的电子占据H2分子的反键轨道,削弱分子内的H-H键,导致H2解离。
8.2 理论预言的高温超导富氢化合物
在实验探索的同时,理论计算也为实验探索指明了一定的方向。关于二元富氢化合物体系已经有非常多的理论预测,其中不乏超导转变温度为高温乃至室温的氢化物。另外,在二元氢化物研究逐渐趋于饱和的同时,科研工作者也将关注点放到了三元氢化物超导体上。相比于二元体系,三元体系因合成路径不同、掺入元素不同,其最终结果也会各有差异。表4为目前理论计算的部分超导转变温度为高温乃至超过室温的二元、三元富氢化合物。
表 4 理论预测的近室温或超室温二元、三元富氢化合物超导体Table 4. Theoretical prediction of binary and ternary polyhydrides with high-temperature superconductivityType Samples Space
groupsPressure/
GPaTc/K Features Binary hydride
superconductorsAcH10[97] R¯3m 200 204–251 Both are predicted to be phonon-mediated superconductors
with an almost empty layer of d atomsAcH16[97] P¯6m2 150 199–241 TbH9[98] C2/c 230 220 With a typical hydrogen cage structure, the coupling of
electrons and hydrogen phonons in the Tb-4f layer
plays a key role in superconductivityTbH10[98] R¯3m 270 270 TbH10[98] Fm¯3m 230 270 ZrH10[96] P63/mmc 250 220 Lamellar alkene H10 junction H3S0.925P0.075[99] Fm¯3m 250 280 Doping calculation based on H3S H3S0.96Si0.04[99] Fm¯3m 250 274 SrH6[100] R¯3m 150 220–235 Hydrogen atoms are clathratelike and form twisted chains SrH10[101] Cmca 300 259 Hydrogen atoms are distributed in staggered
two-dimensional honeycomb layersMgH6[102] Im¯3m 300 420 YH10[10–11, 103] Im¯3m 250 305–326 Unique H32 cagelike structure YH10[10–11, 103] Im¯3m 400 303 YH10[10–11, 103] Im¯3m 300 310 Ternary hydride
superconductorsLi2MgH16[104] Fd¯3m 250 473 Wang et al.[105] pointed that the diffusion of protons between
interstitial spaces may play a key roleCaHfH12[104] Pm¯3m 300 360 The metal skeleton material should be composed of a metal
element with an effective optimal valency of 3, and the metal
should occupy a volume of about 0.4 in the hydrideCaZrH12[106] Pm¯3m 200 290 Mg0.5Ca0.5H6[107] Im¯3m 200 288 Y3CaH24[108] Fm¯3m 150 250 Y3LuH24[109] Fm¯3m 120 283 YLu3H24[109] Fm¯3m 110 288 YLuH12[109] Fd¯3m 140 275 8.3 未来展望
室温超导体是世界各科技强国追逐的目标,富氢化合物是室温超导体的最佳候选体系之一,一旦实现,将成为国家未来战略性材料。目前,领域内主要存在2个亟待解决的问题:(1)真正通过实验获得室温超导体,给出确切的实验证据;(2)在近常压或常压下制备出富氢化合物高温超导体。与此同时,需要积极发展新的实验技术,有望带来富氢化合物的新突破,获得原创性成果。
-
ρ/(g·cm−3) D/(km·s−1) A/GPa B/GPa R1 R2 ω e/(kJ·g−1) 1.630 6.93 371.2 3.21 4.15 0.95 0.3 4.29 ρ/(g·cm−3) C0/(km·s−1) S Γ0 1.024 1.483 1.75 0.28 ρ/(kg·m−3) γ pa/MPa cV/(J·g−1·K−1) 1.17 1.4 0.10 1.012 Seafloor sediment ρ/(kg·m−3) E/MPa ν φ/(°) c/MPa Model 1 1.4 50 0.3 24 0.1 表 5 数值模拟与经验公式结果对比
Table 5. Comparison of numerical simulation and empirical formula results
R/r Peak pressure/MPa Error/% Simulation Empirical formula 7 205.03 193.72 5.84 12 96.32 86.90 10.84 17 60.65 58.62 3.45 22 43.04 43.81 1.75 27 32.58 37.91 14.07 表 6 数值模拟工况设置
Table 6. Settings of simulation cases
Case H/m Seafloor sediment Explosive environment 1 50 Nothing Free field 2 50 Model 1 Near the seabed 3 50 Model 2 Near the seabed 4 100 Model 1 Near the seabed 5 150 Model 1 Near the seabed 表 7 自由场与近海底测得的冲击波峰值压力对比
Table 7. Comparison of peak pressures of free-field and near-seabed underwater explosion shock wave
R/r α/(°) pmax/MPa Reflection
coefficientR/r α/(°) pmax/MPa Reflection
coefficientNear seabed Free field Near seabed Free field 7 20 179.29 201.44 0.89 12 60 92.73 95.67 0.97 7 30 186.21 203.80 0.91 12 70 95.00 97.86 0.97 7 40 193.43 204.80 0.94 12 80 98.33 93.98 1.05 7 50 194.98 204.80 0.95 12 90 92.42 96.30 0.96 7 60 192.78 203.80 0.95 17 20 47.78 58.76 0.81 7 70 195.10 201.43 0.97 17 30 53.24 59.25 0.90 7 80 202.46 201.74 1.00 17 40 56.52 58.58 0.96 7 90 196.41 195.79 1.00 17 50 58.40 58.58 1.00 12 20 82.82 97.87 0.85 17 60 59.44 59.27 1.00 12 30 86.62 95.67 0.91 17 70 61.07 60.10 1.02 12 40 93.96 94.66 0.99 17 80 60.00 61.17 0.98 12 50 96.23 94.65 1.02 17 90 57.77 59.91 0.96 -
[1] BLAKE J R, GIBSON D C. Cavitation bubbles near boundaries [J]. Annual Review of Fluid Mechanics, 1987, 19: 99–123. doi: 10.1146/annurev.fl.19.010187.000531 [2] 张永坤. 水下爆炸对沉底目标毁伤效应试验验证研究 [J]. 舰船电子工程, 2022, 42(5): 154–159. doi: 10.3969/j.issn.1672-9730.2022.05.034ZHANG Y K. Experimental research on the damage effects of bottom torpedoes shocked by underwater explosions [J]. Ship Electronic Engineering, 2022, 42(5): 154–159. doi: 10.3969/j.issn.1672-9730.2022.05.034 [3] 邵建军, 张永坤, 赵红光, 等. 基于相似理论的炸药海中沉底爆炸能量计算 [J]. 四川兵工学报, 2015, 36(6): 124–127. doi: 10.11809/scbgxb2015.06.031SHAO J J, ZHANG Y K, ZHAO H G, et al. Energy calculation of underwater explosion over seabed charge based on similarity law [J]. Journal of Sichuan Ordnance, 2015, 36(6): 124–127. doi: 10.11809/scbgxb2015.06.031 [4] 杨莉, 汪玉, 黄超, 等. 不同水底介质对有限域中装药沉底爆炸特性的影响 [J]. 高压物理学报, 2012, 26(5): 545–550. doi: 10.11858/gywlxb.2012.05.010YANG L, WANG Y, HUANG C, et al. Effects of different grounds on the loading characteristics of limited underwater explosion from a bottom charge [J]. Chinese Journal of High Pressure Physics, 2012, 26(5): 545–550. doi: 10.11858/gywlxb.2012.05.010 [5] 杨莉, 汪玉, 汪斌, 等. 沉底装药水中爆炸现象的实验研究 [J]. 爆炸与冲击, 2013, 33(2): 175–180. doi: 10.11883/1001-1455(2013)02-0175-06YANG L, WANG Y, WANG B, et al. Experimental investigation on loading characteristics of underwater explosion from a bottom charge [J]. Explosion and Shock Waves, 2013, 33(2): 175–180. doi: 10.11883/1001-1455(2013)02-0175-06 [6] 杨莉, 巨圆圆, 武堃, 等. 装药沉底爆炸峰值压力试验及数值模拟研究 [J]. 兵工学报, 2014, 35(Suppl 2): 368–372.YANG L, JU Y Y, WU K, et al. Experiment and simulation on peak pressure of underwater explosion of a bottom charge [J]. Acta Armamentarii, 2014, 35(Suppl 2): 368–372. [7] 杨莉, 汪玉, 杜志鹏, 等. 沉底装药水下爆炸冲击波传播规律 [J]. 兵工学报, 2013, 34(1): 100–104. doi: 10.3969/j.issn.1000-1093.2013.01.018YANG L, WANG Y, DU Z P, et al. Research on shock wave propagation of underwater explosion of bottom charge [J]. Acta Armamentarii, 2013, 34(1): 100–104. doi: 10.3969/j.issn.1000-1093.2013.01.018 [8] 黄潇, 张玮, 张阿漫. 近海底爆炸气泡对潜艇影响的研究 [J]. 中国造船, 2014, 55(3): 25–35. doi: 10.3969/j.issn.1000-4882.2014.03.003HUANG X, ZHANG W, ZHANG A M. Study on the effect of a grounding explosion bubble on submarine [J]. Shipbuilding of China, 2014, 55(3): 25–35. doi: 10.3969/j.issn.1000-4882.2014.03.003 [9] 姚熊亮, 杨文山, 陈娟, 等. 沉底水雷爆炸威力的数值计算 [J]. 爆炸与冲击, 2011, 31(2): 173–178. doi: 10.11883/1001-1455(2011)02-0173-06YAO X L, YANG W S, CHEN J, et al. Numerical calculation of explosion power of mines lying on seabed [J]. Explosion and Shock Waves, 2011, 31(2): 173–178. doi: 10.11883/1001-1455(2011)02-0173-06 [10] 邵宗战, 宋敬利, 贾则. 沉底水雷爆炸威力测量与评估方法研究 [J]. 水雷战与舰船防护, 2012, 20(1): 15–17, 46.SHAO Z Z, SONG J L, JIA Z. Research on the measurement and evaluation methods of the bottom mine’s explosive power [J]. Mine Warfare & Ship Self-Defence, 2012, 20(1): 15–17, 46. [11] 叶林征, 祝锡晶, 王建青, 等. 基于CEL不同角度超声空化微射流冲击的仿真分析 [J]. 振动与冲击, 2016, 35(16): 130–134, 157. doi: 10.13465/j.cnki.jvs.2016.16.021YE L Z, ZHU X J, WANG J Q, et al. Simulations of ultrasonic cavitation micro-jet impact with different angles based on CEL [J]. Journal of Vibration and Shock, 2016, 35(16): 130–134, 157. doi: 10.13465/j.cnki.jvs.2016.16.021 [12] 谭皓洋. 基于CEL方法的舰船近场水下爆炸全过程数值模拟研究 [D]. 哈尔滨: 哈尔滨工程大学, 2017.TAN H Y. Numerical simulation of the whole process of underwater explosion near warship based on CEL method [D]. Harbin: Harbin Engineering University, 2017. [13] AMBROSINI R D, LUCCIONI B M, DANESI R F, et al. Size of craters produced by explosive charges on or above the ground surface [J]. Shock Waves, 2002, 12(1): 69–78. doi: 10.1007/s00193-002-0136-3 [14] LUCCIONI B, AMBROSINI D, NURICK G, et al. Craters produced by underground explosions [J]. Computers & Structures, 2009, 87(21/22): 1366–1373. doi: 10.1016/j.compstruc.2009.06.002 -