First-Principles Study of the High-Pressure Phase Transition and Physical Properties of Rubidium Nitrate
-
摘要: 采用基于密度泛函理论的第一性原理计算,并结合CALYPSO晶体结构预测软件,系统地探索了零温下RbNO3的高压结构与物理性质。基于RbNO3-Ⅳ相的实验数据,比较了4种不同泛函的准确性,发现PBEsol泛函的准确性最高。基于该泛函预测到RbNO3在零温下的相变序列为R3m → Pnma → Pmmn (实验Ⅴ相),相变压强依次为1.7 和8.2 GPa,2次相变均为一级相变,体积坍塌率分别为3.73%和2.54%。研究结果表明,低压下RbNO3可能存在有别于实验(常温常压下)给出的P31结构的低温新相。3个相在各自能量稳定的压强范围内均满足“伯恩-黄昆”弹性稳定性判据,且在整个布里渊区中均不存在声子虚频现象,表明它们均具有动力学稳定性。电子性质分析表明,3个相均为半导体,相变引起的带隙变化普遍较小,加压普遍地抑制了铷离子的电荷向硝酸根离子转移。所预测的高压相变序列及各个相的弹性、晶格动力学、电子结构性质可为后续实验与理论研究提供参考。Abstract: The high-pressure structure and physical properties of RbNO3 at zero temperature was systematically explored using first-principles calculations based on density generalized theory combined with CALYPSO crystal structure predictions. The accuracy of four different functionals was compared based on experimental data of the RbNO3-Ⅳ phase, and the revised PBE for solids (PBEsol) functional was found to be the most reliable. The zero-temperature phase transition sequence of RbNO3 predicted is R3m→Pnma→Pmmn (experimental phaseⅤ) based on the PBEsol functional, and the two phase transition pressures are 1.7 and 8.2 GPa. The two phase transitions are first-order phase transitions, and the volume collapse rates reach 3.73% and 2.54%, respectively. This suggests that RbNO3 at low pressure may have new low-temperature phases those are different from P31 structure given by experiment in the room temperature and ambient pressure. In the energy-stabilized pressure interval, all three phases satisfy the “Born-Huangkun” elastic stability criteria, and there is no phonon virtual frequency phenomenon in the whole Brillouin zone, which indicates that they are dynamically stable structures. The electronic property analysis showed that the three phases are semiconductors, and the band gap changes caused by the phase transition are generally small, but the pressure generally inhibits the charge transfer from alkali metal ions to nitrate ions. The high-pressure phase transition sequence predicted in this paper and the elasticity, lattice dynamics, and electronic structure properties of the individual phases can provide a reference for subsequent experimental and theoretical studies.
-
Key words:
- high pressure /
- phase transition /
- first principles /
- RbNO3
-
盐类矿物是指以矿石或岩石形式存在于自然界的碱金属/碱土金属的卤化物、硫酸盐、碳酸盐、硝酸盐等,被广泛应用于工业和农业领域[1]。该类矿物在高温高压下的结构相变在物理与地球科学领域备受关注,例如:2013 年理论与实验研究发现,典型离子晶体 NaCl 在高压下转变为 Na3Cl和NaCl3[2];Santos等[3]发现,作为地球深部碳循环重要载体的碳酸盐在高温高压下具有10余种结构。硝酸根离子(
NO−3 )和碳酸根离子(CO2−3 )具有等电子结构,形成的盐类矿物均为离子晶体。离子晶体中的离子键易受外部温度和压强调控而发生结构相变,鉴于NO−3 和CO2−3 在结构和电子性质上具有相似之处,高温高压下碱金属硝酸盐是否具有丰富的结构相变是值得深入探讨的问题。碱金属硝酸盐被广泛用于多个领域,如作为传热、储热介质等运用于工业炸药、冶金淬火硬化、太阳能和核能领域[4–6]。碱金属硝酸盐在较低压强下已经展现出丰富的相变行为[7],引起了许多关注。研究发现,在常压下硝酸铷(RbNO3)从熔融态退火至常温的过程中历经了4个不同的结构,即Ⅰ、Ⅱ、Ⅲ、Ⅳ相。Ⅰ相在熔点至564 K的区间内稳定,为NaCl型的立方无序结构(空间群Fm
¯3 m)[8];Ⅱ相在564~492 K区间稳定,具有菱方结构(空间群R¯3 m)[7];Ⅲ相在492~437 K区间稳定,一般认为是CsCl型立方结构(空间群Pm¯3 m)[8];Ⅳ相在常温下稳定,为菱方结构(空间群P31)[7]。2001年,Liu等[9]利用分子动力学模拟方法从理论上发现了RbNO3在Ⅳ→Ⅲ相变过程中NO−3 基团从定向有序变为无序。在高压实验研究方面,1984年Dean等[10]通过原位X射线衍射研究发现了RbNO3的Ⅳ→Ⅴ结构相变,其相变压强约为1.2 GPa,高压Ⅴ相为正交结构(空间群Pmmn)。与RbNO3的高温研究相比,RbNO3的高压研究并不丰富。本研究将基于第一性原理计算,结合CALYPSO晶体结构预测软件,系统地探索零温下RbNO3在0~100 GPa区间的结构与物理性质。
1. 计算方法
在0 K和0、20、40、60、80和100 GPa的压强条件下,使用CALYPSO软件[11–12]进行RbNO3的结构搜索。结构预测过程中演化代数为20,模拟采用2倍胞和4倍胞,每个分子式的体积在32~70 Å3区间,体积随压强的升高而减小。基于第一性原理计算,交换关联泛函采用广义梯度近似(generalized gradient approximation,GGA)[13]的Perdew-Burke-Ernzerhof(PBE)泛函[13]、Original vdW-DF泛函[14]、Revised PBE for solids(PBEsol)泛函[15]以及Revised PBE Hammer et al.(RPBE)泛函[16],数值模拟在VASP[17]软件包中实现,采用投影缀加波(projected augmented wave,PAW)[18]方法,Rb、N、O的价电子分别选为4s24p65s1、2s22p3和2s22p4,平面波截断能取520 eV,对于R3m、Pnma和Pmmn相,布里渊区采样的网格尺寸分别为10×10×10、6×8×6和10×8×10,总能量收敛至小于1×10−5 eV/atom。采用有限位移法结合PHONOPY软件[19]计算声子谱,扩胞数分别为3×3×2、2×3×2和3×3×3。弹性模量计算采用Voigt-Reuss-Hill(VRH)平均法[20–21],晶体结构和声子模式的原子位移可视化通过VESTA软件[22]和V_Sim软件[23]实现。
2. 结果与讨论
RbNO3的常温常压结构为RbNO3-Ⅳ相,空间群为P31。采用4种不同交换关联泛函优化0 GPa下RbNO3-Ⅳ相的平衡体积与晶胞参数,结果列于表1。对比实验结果可以发现,4种交换关联泛函的准确性从低到高依次为RPBE、vdW-DF、 PBE、 PBEsol。因此,后续研究将主要基于PBE和PBEsol泛函展开。
表 1 RbNO3-Ⅳ相的晶格常数的理论值与实验值的比较Table 1. Comparison between experimental and calculated lattice parameters of RbNO3-Ⅳ phaseMethod Condition a/Å c/Å Volume/Å3 Simulation PBE 10.68(+1.62%) 7.60(+1.88%) 867.56(+5.29%) vdW-DF 10.72(+2.00%) 7.62(+2.14%) 876.29(+6.35%) PBEsol 10.39(−1.14%) 7.39(−0.94%) 797.69(−3.19%) RPBE 11.14(+5.99%) 7.91(+6.03%) 982.63(+19.26%) Experiment T=296 K[24] 10.47 7.44 815.58 T=298 K[25] 10.55 7.47 831.43 Average 10.51 7.46 823.96 采用PBE和PBEsol泛函在0~100 GPa压强范围内对能量较低的备选结构进行精细优化后发现,R3m、Pnma、P31、Pmmn相在不同的压强区间内具有热力学稳定性,其焓差随压强变化曲线如图1所示。在0 GPa下,2种泛函的预测结果均表明R3m相的能量最稳定;随压强的升高,R3m相均转变为Pnma相。继续加压,PBE泛函与PBEsol泛函预测的相变序列不同:PBE泛函预测结果表明,Pnma相先相变至Ⅳ相(即P31结构),再相变至Ⅴ相(即Pmmn结构);PBEsol泛函预测结果表明,Pnma相直接转变为Ⅴ相。压强-体积关系如图2所示,其中,ΔV/V为体积坍塌率。除PBE泛函预测的Pnma→Ⅳ相变中不能看到明显的体积坍缩以外,PBE泛函所预测的其余2个相变和PBEsol预测的2个相变均存在明显的体积坍塌,这种体积不连续表明它们均为一级相变。
考虑到PBEsol泛函对实验Ⅳ相的平衡体积与晶胞参数的描述更为准确,后续研究将主要采用PBEsol泛函预测零温下RbNO3的相变序列(R3m→Pnma→Pmmn)。该序列所对应的2个相变点的压强分别为1.7和8.2 GPa,相应的体积坍塌率分别为3.73%和2.54%。由于碳酸钙的结构与碱金属硝酸盐的结构相似,两者具有等电子结构,因此,对比了它们的高压相变序列。碳酸钙的相变序列为:R
¯3 c→Pnma→Pmmn→P21/c,相变压强分别为1、40和67 GPa[26]。可以发现,RbNO3的Pnma→Pmmn的相变压强远低于CaCO3发生相同相变时的压强,其可能的原因如下:一方面,Rb+与NO−3 之间的离子相互作用比Ca2+与CO2−3 之间的离子作用更弱;另一方面,Rb+的半径大于Ca2+的半径,Rb+产生的“化学预压缩”更强。上述零温理论预测出了能量比常温常压下的Ⅳ相更小的2个低压新相(R3m相和Pnma相),表明低压下RbNO3可能存在不同于实验给出的常温常压下Ⅳ相的低温新相。预测获得的0 GPa下的R3m相、4 GPa下的Pnma相和12 GPa下的Pmmn相的平衡态晶格常数、原子位置列于表2,相应的晶体结构见图3。3个结构存在明显差异,三方相的单胞包含3个RbNO3分子式,而正交相的单胞则分别包含4和2个RbNO3分子式。3个相在0、4、12 GPa压强下的弹性常数列于表3,可以发现,三方相和正交相分别有7和9个独立的弹性常数。在R3m相中,C44随压强的升高先减小后增大,其余的弹性常数均与压强呈正相关。3个相的弹性模量,即体积模量B、剪切模量G以及杨氏模量E也列于表3中,分别表示抵抗脆性变形、抵抗塑性变形和抵抗张力的能力。同一压强下Pmmn相的弹性模量最大,表明Pmmn相的抗变形能力比另2个相更强。此外,根据“伯恩-黄昆”弹性稳定性判据[27–28]分析3个相的机械稳定性,判别条件见式(1)和式(2)。力稳定性分析显示:12 GPa时,Pnma相的C55 < 0,不满足弹性稳定性判据,因而Pnma相不是弹性稳定结构;除此之外,各相在相应压强下的弹性常数均满足弹性稳定判据。
表 2 高压下R3m、Pnma和Pmmn相的晶格参数和原子位置Table 2. Lattice parameters and atomic coordinates of R3m, Pnma and Pmmn phases at different pressuresPhase Pressure/GPa Lattice parameters Wyckoff positions R3m 0 a=b=5.64 Å, c=9.48 Å,
α=β=90°, γ=120°Rb1: 3a(0.3333, 0.6667, 0.1699)
N1: 3a(0.3333, 0.6667, 0.7275)
O1: 9b(0.5373, 0.0746, 0.3953)Pnma 4 a=7.18 Å, b=5.61 Å, c=6.77 Å,
α=β=γ=90°Rb1: 4c(0.483, 0.250, 0.686)
N1: 4c(0.350, 0.250, 0.122)
O1: 4c(0.448, 0.250, 0.275)
O5: 8d(0.302, 0.444, 0.043)Pmmn 12 a=4.67 Å, b=5.34 Å, c=4.63 Å,
α=β=γ=90°Rb1: 2b(0.500, 0, 0.400)
N1: 2a(0, 0, 0.992)
O1: 2a(0, 0, 0.721)
O3: 4e(0, 0.204, 0.128)表 3 高压下R3m、Pnma和Pmmn相的弹性常数和弹性模量Table 3. Elastic constant and elastic moduli of R3m, Pnma and Pmmn phase at high pressureGPa Phase Pressure C11 C12 C13 C14 C22 C23 C33 C44 C55 C66 B G E R3m 0 37 14 11 −4 13 7 12 15 7 19 4 75 28 21 −8 39 6 23 35 11 29 12 124 46 42 −15 94 12 39 66 19 52 Pnma 0 22 7 9 41 9 32 7 10 9 16 10 24 4 50 19 30 73 19 69 12 17 19 36 17 44 12 96 43 59 129 35 119 19 −5 34 69 −57 −236 Pmmn 0 26 9 9 50 16 54 15 5 5 21 10 26 4 60 16 25 84 32 83 33 20 9 40 21 53 12 113 25 53 143 59 129 60 37 16 72 35 91 C11>|C12|,C44>0,C213<12C33(C11+C12),C214<12C44(C11−C12) (1) C11>0,C11C22>C212,C11C22C33+2C12C13C23−C11C223−C22C213−C33C212>0C44>0,C55>0,C66>0 (2) 为了确定RbNO3的3个相的动力学稳定性,采用有限位移法计算了它们的声子谱。若声子谱具有负频率,则意味着动力学不稳定[29]。在简谐近似下,分别计算了R3m、Pnma和Pmmn相在不同压强下的声子色散曲线以及声子分波态密度(phonon density of states, PHDOS),结果如图4所示。可以看出,声子谱均未出现负频率,表明3个相在相应压强下均具有动力学稳定性。
声子态密度分析研究表明,RbNO3的3个相的声子态密度分布趋势皆与常压下的CaCO3[30]十分相似,例如:当RbNO3(CaCO3)的振动频率低于10 THz(12 THz)时,声子态密度的贡献主要来自于O原子与金属阳离子Rb+(Ca2+);在振动频率32 THz附近,声子态密度的贡献主要来自于O原子;在振动频率大于40 THz的高频区,声子态密度的贡献来自于N―O键或者C―O键。采用V_Sim软件对Γ点处声子的振动模式进行可视化分析,结果表明:当频率小于10 THz时,N原子与O原子之间的距离不变,Rb和
NO−3 协同振动;在32 THz附近,R3m和Pnma相的N―O振动和C―O振动均为对称伸缩振动,与方解石和文石的特性类似;在21 THz附近,N―O振动和C―O振动均为平面内弯曲振动;在RbNO3(CaCO3)的振动频率约为24 THz(25 THz)时,N―O振动(C―O振动)均为平面外弯曲振动;在大于40 THz的范围内,N―O振动和C―O振动均为反对称弯曲振动。相对于0 GPa的R3m相和4 GPa的Pnma相而言,12 GPa的Pmmn相的声子振动频率整体升高,但不同模式的频率分布范围大致相同。为了研究R3m、Pnma和Pmmn相中原子间的成键行为,计算了0 GPa的R3m相、4 GPa的Pnma相和12 GPa的Pmmn相的电子局域函数(electron localization functionm, ELF),计算结果如图5所示。三维ELF图中等值面的取值为0.8,二维ELF图的晶格切面分别为(1.000 14.005 31.363)、(20.145 2.048 −1.000)和(4.542 1.000 −1.516)。可以看出,N原子与O原子间的ELE值趋近于1,表明N―O之间存在强共价键。
此外,计算了0 GPa的R3m相、4 GPa的Pnma相和12 GPa的Pmmn相的晶体轨道哈密顿布居(crystal orbital Hamiton population, COHP),计算结果如图6所示。COHP为负表示成键态。相对于O―O键和N―O键而言,Rb―O键对COHP的贡献非常小,几乎可以忽略不计,符合Rb+与
NO−3 之间为离子相互作用的事实。在−1.5~0 eV范围内,O―O键的COHP为正值,表明其主要来源为反键轨道,包含NO−3 从Rb+处获取的1个价电子的贡献。在0 eV以下,N―O键的COHP几乎全部为负值,表明NO−3 中的N―O键是极为稳定的强共价键,与ELF的分析结果一致。在0 GPa的R3m相、4 GPa的Pnma相和12 GPa的Pmmn相中,N原子与O原子的键长分别为1.26、1.25和1.26 Å。计算Bader电荷转移以确定Rb原子、N原子和O原子之间的电子转移情况,计算结果如表4所示。这3个结构中,硝酸根离子分别从Rb+获得了0.90e、0.87e和0.86e的电荷,表明压强抑制Rb+的电荷向NO−3 离子转移。表 4 高压下R3m、Pnma和Pmmn相的元胞内的Bader电荷转移Table 4. Calculated Bader charges within R3m, Pnma and Pmmn phases primitive cells at different pressuresPhase Pressure/GPa Atoms Number Charge/e Charge transfer/e R3m 0 Rb 1 8.10 0.90 N 1 4.16 0.84 O 3 6.58 −0.58 Pnma 4 Rb 4 8.13 0.87 N 4 4.13 0.87 O1 4 6.57 −0.57 O2 4 6.58 −0.58 O3 4 6.59 −0.59 Pmmn 12 Rb 2 8.14 0.86 N 2 4.07 0.93 O1 2 6.58 −0.58 O2 2 6.59 −0.59 O3 2 6.62 −0.62 图7给出了0 GPa的R3m相、4 GPa的Pnma相和12 GPa的Pmmn相的电子能带结构和电子态密度( density of states, DOS)。电子态密度结果表明,价带主要来自O的p轨道,导带主要来自N-p与O-p的杂化。3个结构在相应压强下的带隙分别为3.06、3.40和3.10 eV。带隙的存在说明它们均为半导体材料,然而3个结构的带隙类型不同,其中,R3m结构具有直接带隙,Pnma和Pmmn结构具有间接带隙。此外,还可以发现,在结构相变过程中,带隙随压强的增加呈现非单调变化的规律,即第1次相变带隙增大,而第2次相变带隙减小。
3. 结 论
采用第一性原理计算与CALYPSO结构预测相结合的方法系统地研究了RbNO3在0~100 GPa压强范围内的高压结构与物理性质。以RbNO3常温常压结构的实验数据为参照,测试了PBE、vdW-DF、RPBE、PBEsol泛函对结构描述的准确性,发现PBEsol泛函给出的结果最准确。基于PBEsol泛函预测出零温下RbNO3的高压相变序列为R3m→Pnma→Pmmn,相变压强分别约为1.7 和8.2 GPa;2次相变均为一级结构相变,相应的体积坍塌率分别为3.73%和2.54%。预测零温下RbNO3存在2个低压新相,即R3m相和Pnma相,它们在能量上均优于早期实验在常温常压下给出的P31结构,表明降温可能会诱导RbNO3发生结构相变。0 GPa的R3m相、4 GPa的Pnma相和12 GPa的Pmmn相均具有力学稳定性和动力学稳定性,其带隙分别为3.06、3.40和3.10 eV,对比分析发现,压强会抑制Rb+的电荷向
NO−3 转移。本研究获得的硝酸铷结构与电子性质随压强的变化规律可以为其他碱金属硝酸盐的理论与实验研究提供参考。 -
表 1 RbNO3-Ⅳ相的晶格常数的理论值与实验值的比较
Table 1. Comparison between experimental and calculated lattice parameters of RbNO3-Ⅳ phase
Method Condition a/Å c/Å Volume/Å3 Simulation PBE 10.68(+1.62%) 7.60(+1.88%) 867.56(+5.29%) vdW-DF 10.72(+2.00%) 7.62(+2.14%) 876.29(+6.35%) PBEsol 10.39(−1.14%) 7.39(−0.94%) 797.69(−3.19%) RPBE 11.14(+5.99%) 7.91(+6.03%) 982.63(+19.26%) Experiment T=296 K[24] 10.47 7.44 815.58 T=298 K[25] 10.55 7.47 831.43 Average 10.51 7.46 823.96 表 2 高压下R3m、Pnma和Pmmn相的晶格参数和原子位置
Table 2. Lattice parameters and atomic coordinates of R3m, Pnma and Pmmn phases at different pressures
Phase Pressure/GPa Lattice parameters Wyckoff positions R3m 0 a=b=5.64 Å, c=9.48 Å,
α=β=90°, γ=120°Rb1: 3a(0.3333, 0.6667, 0.1699)
N1: 3a(0.3333, 0.6667, 0.7275)
O1: 9b(0.5373, 0.0746, 0.3953)Pnma 4 a=7.18 Å, b=5.61 Å, c=6.77 Å,
α=β=γ=90°Rb1: 4c(0.483, 0.250, 0.686)
N1: 4c(0.350, 0.250, 0.122)
O1: 4c(0.448, 0.250, 0.275)
O5: 8d(0.302, 0.444, 0.043)Pmmn 12 a=4.67 Å, b=5.34 Å, c=4.63 Å,
α=β=γ=90°Rb1: 2b(0.500, 0, 0.400)
N1: 2a(0, 0, 0.992)
O1: 2a(0, 0, 0.721)
O3: 4e(0, 0.204, 0.128)表 3 高压下R3m、Pnma和Pmmn相的弹性常数和弹性模量
Table 3. Elastic constant and elastic moduli of R3m, Pnma and Pmmn phase at high pressure
GPa Phase Pressure C11 C12 C13 C14 C22 C23 C33 C44 C55 C66 B G E R3m 0 37 14 11 −4 13 7 12 15 7 19 4 75 28 21 −8 39 6 23 35 11 29 12 124 46 42 −15 94 12 39 66 19 52 Pnma 0 22 7 9 41 9 32 7 10 9 16 10 24 4 50 19 30 73 19 69 12 17 19 36 17 44 12 96 43 59 129 35 119 19 −5 34 69 −57 −236 Pmmn 0 26 9 9 50 16 54 15 5 5 21 10 26 4 60 16 25 84 32 83 33 20 9 40 21 53 12 113 25 53 143 59 129 60 37 16 72 35 91 表 4 高压下R3m、Pnma和Pmmn相的元胞内的Bader电荷转移
Table 4. Calculated Bader charges within R3m, Pnma and Pmmn phases primitive cells at different pressures
Phase Pressure/GPa Atoms Number Charge/e Charge transfer/e R3m 0 Rb 1 8.10 0.90 N 1 4.16 0.84 O 3 6.58 −0.58 Pnma 4 Rb 4 8.13 0.87 N 4 4.13 0.87 O1 4 6.57 −0.57 O2 4 6.58 −0.58 O3 4 6.59 −0.59 Pmmn 12 Rb 2 8.14 0.86 N 2 4.07 0.93 O1 2 6.58 −0.58 O2 2 6.59 −0.59 O3 2 6.62 −0.62 -
[1] LAUE W, THIEMANN M, SCHEIBLER E, et al. Nitrates and nitrites [J]. Ullmann’s Encyclopedia of Industrial Chemistry, 2000, 24: 153–157. [2] MIAO M S, SUN Y H, ZUREK E, et al. Chemistry under high pressure [J]. Nature Reviews Chemistry, 2020, 4(10): 508–527. doi: 10.1038/s41570-020-0213-0 [3] SANTOS S S M, MARCONDES M L, JUSTO J F, et al. Calcium carbonate at high pressures and high temperatures: a first-principles investigation [J]. Physics of the Earth and Planetary Interiors, 2020, 299: 106327. doi: 10.1016/j.pepi.2019.106327 [4] LEE W E, GILBERT M, MURPHY S T, et al. Opportunities for advanced ceramics and composites in the nuclear sector [J]. Journal of the American Ceramic Society, 2013, 96(7): 2005–2030. doi: 10.1111/jace.12406 [5] SOLEIMANPOUR S, SADRAMELI S M, MOUSAVI S A H S, et al. Preparation and characterization of high temperature shape stable NaNO3/diatomite phase change materials with nanoparticles for solar energy storage applications [J]. Journal of Energy Storage, 2022, 45: 103735. doi: 10.1016/j.est.2021.103735 [6] DIÉGUEZ O, VANDERBILT D. Theoretical study of ferroelectric potassium nitrate [J]. Physical Review B, 2007, 76(13): 134101. doi: 10.1103/PhysRevB.76.134101 [7] RAO C N R, PRAKASH B, NATARAJAN M. Crystal structure transformations in inorganic nitrites, nitrates, and carbonates: NSRDS-NBS 53 [R]. Washington, USA: National Bureau of Standards, 1975. [8] AHTEE M, HEWAT A W. Structures of the high temperature phases of rubidium nitrate [J]. Physica Status Solidi (A), 1980, 58(2): 525–531. doi: 10.1002/pssa.2210580224 [9] LIU J J, DUAN C G, OSSOWSKI M M, et al. Molecular dynamics simulation of structural phase transitions in RbNO3 and CsNO3 [J]. Journal of Solid State Chemistry, 2001, 160(1): 222–229. doi: 10.1006/jssc.2001.9226 [10] DEAN C, HAMBLEY T W, SNOW M R. Structures of phase Ⅳ rubidium nitrate, RbNO3, and phase Ⅱ caesium nitrate, CsNO3 [J]. Acta Crystallographica Section C, 1984, 40(9): 1512–1515. [11] WANG Y C, LV J, ZHU L, et al. Crystal structure prediction via particle-swarm optimization [J]. Physical Review B, 2010, 82(9): 094116. doi: 10.1103/PhysRevB.82.094116 [12] 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. [13] 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 [14] DION M, RYDBERG H, SCHRÖDER E, et al. Van der Waals density functional for general geometries [J]. Physical Review Letters, 2004, 92(24): 246401. doi: 10.1103/PhysRevLett.92.246401 [15] PERDEW J P, RUZSINSZKY A, CSONKA G I, et al. Restoring the density-gradient expansion for exchange in solids and surfaces [J]. Physical Review Letters, 2008, 100(13): 136406. doi: 10.1103/PhysRevLett.100.136406 [16] HAMMER B, HANSEN L B, NØRSKOV J K. Improved adsorption energetics within density-functional theory using revised Perdew-Burke-Ernzerhof functionals [J]. Physical Review B, 1999, 59(11): 7413–7421. doi: 10.1103/PhysRevB.59.7413 [17] 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–11186. doi: 10.1103/PhysRevB.54.11169 [18] BLÖCHL P E. Projector augmented-wave method [J]. Physical Review B, 1994, 50(24): 17953–17979. doi: 10.1103/PhysRevB.50.17953 [19] TOGO A, TANAKA I. First principles phonon calculations in materials science [J]. Scripta Materialia, 2015, 108: 1–5. doi: 10.1016/j.scriptamat.2015.07.021 [20] HILL R. The elastic behaviour of a crystalline aggregate [J]. Proceedings of the Physical Society Section A, 1952, 65(5): 349–354. doi: 10.1088/0370-1298/65/5/307 [21] OUADHA I, RACHED H, AZZOUZ-RACHED A, et al. Study of the structural, mechanical and thermodynamic properties of the new MAX phase compounds (Zr1−xTix)3AlC2 [J]. Computational Condensed Matter, 2020, 23: e00468. doi: 10.1016/j.cocom.2020.e00468 [22] MOMMA K, IZUMI F. VESTA 3 for three-dimensional visualization of crystal, volumetric and morphology data [J]. Journal of Applied Crystallography, 2011, 44(6): 1272–1276. doi: 10.1107/S0021889811038970 [23] MISHRA K K, CHANDRA S, SALKE N P, et al. Soft modes and anharmonicity in H3[Co(CN)6]: Raman spectroscopy and first-principles calculations [J]. Physical Review B, 2015, 92(13): 134112. doi: 10.1103/PhysRevB.92.134112 [24] POHL J, POHL D, ADIWIDJAJA G. Phase transition in rubidium nitrate at 346 K and structure at 296, 372, 413 and 437 K [J]. Acta Crystallographica Section B, 1992, 48(2): 160–166. doi: 10.1107/S0108768191013459 [25] SHAMSUZZOHA M, LUCAS B W. Structure (neutron) of phase Ⅳ rubidium nitrate at 298 and 403 K [J]. Acta Crystallographica Section B, 1982, 38(9): 2353–2357. doi: 10.1107/S0567740882008772 [26] SMITH D, LAWLER K V, MARTINEZ-CANALES M, et al. Postaragonite phases of CaCO3 at lower mantle pressures [J]. Physical Review Materials, 2018, 2(1): 013605. doi: 10.1103/PhysRevMaterials.2.013605 [27] WU Z J, ZHAO E J, XIANG H P, et al. Crystal structures and elastic properties of superhard IrN2 and IrN3 from first principles [J]. Physical Review B, 2007, 76(5): 054115. doi: 10.1103/PhysRevB.76.054115 [28] MOUHAT F, COUDERT F X. Necessary and sufficient elastic stability conditions in various crystal systems [J]. Physical Review B, 2014, 90(22): 224104. doi: 10.1103/PhysRevB.90.224104 [29] BORN M, HUANG K, LAX M. Dynamical theory of crystal lattices [J]. American Journal of Physics, 1955, 23(7): 474. [30] UKITA M, TOYOURA K, NAKAMURA A, et al. Pressure-induced phase transition of calcite and aragonite: a first principles study [J]. Journal of Applied Physics, 2016, 120(14): 142118. doi: 10.1063/1.4961723 -