First-Principles Theoretical Study on the Structure Behaviors of NaPO3 under Compression
-
摘要: 探索PO6配位八面体的高压晶体化学行为是理解磷的高压化学性质、了解磷在下地幔中可能的赋存方式及磷的地球化学循环的重要基础。在0~80 GPa压力范围内,对MgSiO3等电子体的NaPO3开展第一性原理密度泛函理论研究,通过对其常压β相(P21/n)、透辉石相(C2/c)、钛铁矿相(R
¯3 )、斜方钙钛矿相(Pnma)和立方钙钛矿相(Pm3m) 的几何优化和总能对比,获得了NaPO3的结构相变序列及相变压力:P21/n → C2/c (2 GPa) → R¯3 (20 GPa)→ Pnma(50 GPa),相变导致的体积变化分别为7.1%、11.5%和9.0%。Pm3m-NaPO3的声子色散曲线在R点和M点呈现出显著且相似的虚频,而Pnma-NaPO3在整个布里渊区均表现为实频,表明Pnma-NaPO3动力学稳定。Pnma-NaPO3的晶格常数、P―O键长、P―O―P键角、NaO12 和PO6 多面体体积比VNaO12 /VPO6 与压力的关系表明,PO6八面体在计算的整个压力范围内都较规则,且NaO12多面体的压缩性比PO6八面体的压缩性更大。电子结构计算表明,在Pnma-NaPO3的PO6八面体中,P的3p和3s轨道与O的2p轨道强烈混合,P―O键表现出的强共价性对稳定其斜方钙钛矿结构发挥了关键作用。Abstract: Exploring the high-pressure crystal chemical behaviors of the PO6 coordinated octahedron is an important basis for understanding the high-pressure chemistry, the possible occurrence in the lower mantle, and the geochemical cycle of the phosphorus element. In this study, NaPO3, which is isoelectronic with the major component of the lower mantle MgSiO3, was studied with the first-principle density functional theory in the pressure range of 0–80 GPa. By ways of geometric optimization and total energy comparison of its ambient pressure β phase (P21/n), diopside phase (C2/c), ilmenite phase (R¯3 ), orthorhombic (Pnma) and cubic (Pm3m) perovskite phases, the structural phase transformation sequence and phase transformation pressures were obtained: P21/n→C2/c (2 GPa)→R¯3 (20 GPa)→Pnma (50 GPa), with the unit-cell volume changes of 7.1%, 11.5% and 9.0%, respectively. The phonon dispersion curves of Pm3m-NaPO3 show remarkable and similar imaginary frequencies at R and M points, while the orthorhombic perovskite structure shows real frequencies throughout the whole Brillouin zone reflecting its dynamic stability. The pressure dependence of lattice constants, P―O bond lengths, P―O―P bond angles andVNaO12 /VPO6 polyhedron volume ratio of Pnma-NaPO3 shows that the PO6 octahedron is regular in the whole calculated pressure range, and the compressibility of NaO12 polyhedron is greater than that of PO6 octahedron. The electronic structure calculation shows that the 3p and 3s orbitals of P are strongly mixed with 2p orbitals of O in the PO6 octahedron of Pnma-NaPO3, and the P―O bond exhibits strong covalency, which plays a key role in stabilizing the orthorhombic perovskite structure. -
磷(P)的价电子结构为3s23p3,电负性为2.19(Pauling标度)。在地球化学分类中,磷是一种高度不相容元素,在地壳氧化还原条件下,磷主要以+5价态与氧形成PO4配位四面体赋存于副矿物磷灰石中[1]。在玄武岩部分熔融实验中,磷在主要矿物/熔体间的分配系数
DMineral/meltP <0.15,这些硅酸盐矿物/溶体间磷的分配系数从大到小依次为石榴子石、橄榄石、单斜辉石、斜方辉石[2]。考虑到橄榄石、石榴子石和斜方辉石等硅酸盐矿物是地幔岩石的主要成分以及磷在地幔中的估计丰度值约为9.0×10−5,并且石榴子石与磷的高温高压实验显示,磷会以类质同象替代的方式进入石榴子石地幔中,在地幔中不太可能出现磷灰石独立矿物[3–5]。因此,硅酸盐矿物是磷在地幔中的主要寄主矿物。由于磷在橄榄石等矿物中扩散非常慢,玄武岩中的橄榄石晶体常出现磷的环带现象,可为岩浆演化过程和矿物结晶动力学研究提供重要依据[6–7]。石榴子石与磷的溶解度高温高压实验结果表明,当压力大于3 GPa时,磷在石榴子石中的溶解度显著增大,可以解释高压变质岩折返后出现的磷灰石出溶现象[5, 8–10]。已有实验结果和野外观察表明,在地球上地幔的高温高压环境中,磷主要以类质同象替代的方式替代硅进入石榴石、橄榄石等主要硅酸盐矿物晶格中,且以PO4四面体形式存在。但磷在地球下地幔中的可能存在方式研究几乎空白。地球下地幔主要由2种钙钛矿结构硅酸盐矿物(布里奇曼石bridgmanite和毛硅钙石davemaoite)和方镁铁矿组成,探讨下地幔环境中磷能否进入钙钛矿结构硅酸盐矿物的晶格中,关键是了解下地幔温压条件下磷能否与氧形成PO6配位八面体,并了解PO6配位八面体是否具有与SiO6八面体相似的晶体化学特征。迄今为止,高压下PO6八面体的性质仅有少量的实验和理论研究报道,主要结果总结如下。常压下柏林石(berlinite,AlPO4)是一种类似石英(SiO2)结构的矿物,在约46 GPa下经高温处理后开始转变成类似斯石英的畸变结构(P2/m),在68.4 GPa压力下,其结构中PO6八面体的P―O平均键长约为1.65 Å[11]。SiO2和AlPO4混合物经
1675 ℃、18 GPa的温压条件处理后,质量分数大约为1%的P2O5进入斯石英中,X射线吸收近边结构和核磁共振测试均显示磷以PO6配位八面体存在[12–13]。常压下为单斜结构的磷钙铍石(hurlbutite,CaBe2P2O8,P21/c)随压力的增加会经历一系列结构相变,在约90 GPa下转变成三斜结构Ⅳ相(P¯1 ),该结构中,磷占据4种不同的结晶学位置,且均与氧形成PO6八面体配位,P―O键长的变化范围为1.49~1.77 Å,PO6八面体的P―O平均键长约在1.62~1.65 Å[14]范围内。P2O5的高压结构相变行为的理论研究表明,其在约20~30 GPa区间转变成具有PO6八面体配位的单斜结构(C2/c)[15]。TbPO4的理论研究表明,随着压力的增加,TbPO4发生由锆石结构(I41/amd)向独居石结构(P21/n)、白钨矿结构(I41/a)转变,并在约64 GPa下转变成斜方SrUO4型结构(Pbcm),该结构中PO6八面体共角顶连接成层状,在约77 GPa下,P―O键的键长范围为1.55~1.72 Å[16]。InPO4和TiPO4的理论研究显示,它们在约50 GPa下转变成黑钨矿结构,结构中PO6八面体共棱连接,在约60 GPa下P―O键的平均键长约为1.65 Å[17]。下地幔矿物模型以斜方钙钛矿结构(Pnma)的布里奇曼石占主导,其主要组成为MgSiO3。从元素周期表上可以看出,NaPO3和MgSiO3为等电子体,因而,它们可能具有相似的结构和物理化学性质。因此,NaPO3是研究高压下PO6八面体的化学性质和磷在下地幔中的赋存方式的良好对象。
常压下,NaPO3可形成单斜(P21/n)、斜方(Pmcn)和四方(I41/a)等多种结构,如图1所示。在这些结构中,磷均与氧形成PO4四面体,其中,在单斜和四方结构中,PO4四面体形成类似硅酸盐辉石类矿物的单链结构,而斜方相为PO4四面体的三元环结构[18–21]。迄今已发现的磷酸盐矿物有400多种,基本上都以孤立的PO4四面体为基本结构单元[1]。借助化学和材料科学,大量以PO4四面体相互连接的链、环等为基本结构单元的材料被合成出来[22–23]。此外,PO4四面体链状结构也是生物体的基本组成单元,是携带能量的基本结构[24–25]。然而,迄今为止,鲜见碱金属磷酸盐化合物的高压行为研究报道。
因此,本研究将对NaPO3开展高压下结构和相变的第一性原理理论研究,特别关注其在结构和物性方面与MgSiO3的相似性,为探索下地幔中磷的赋存形式提供参考依据。
1. 理论研究方法
应用Materials Studio(2017)软件中的CASTEP模块对NaPO3开展密度泛函理论研究。本研究采用GGA-PBEsol交换关联泛函、OTFG norm conserving赝势,截断能取990 eV,布里渊区k点的取样密度为0.04 Å−1。几何优化收敛精度设定如下:能量取5×10−6 eV,最大力取0.01 eV/Å,最大应力取0.01 GPa,最大位移取5.0×10−4 Å。计算的价电子包括Na 2s22p63s1、P 3s23p3、O 2s22p4。
基于NaPO3常压结构的前期能量计算及其实验合成条件,发现NaPO3的几种常压结构具有相近的能量。考虑到常压结构的低对称性和大晶胞对计算量的要求,选择常压相中的高密度β型P21/n(Z=8,Z为单位晶胞中NaPO3化学式个数)进行相应的理论计算。同时,参照MgSiO3的高压相变序列[26],构建了NaPO3的常压相β型P21/n结构(Z=8)[19]、透辉石型结构(C2/c,Z=8)[27]、钛铁矿型结构(R
¯3 ,Z=6)[28]和钙钛矿型结构(立方相Pm3m(Z=1)和斜方相Pnma(Z=4)[29])开展几何优化和总能计算,以了解NaPO3的相变压力条件和可能的热力学稳定相。为了解钙钛矿结构NaPO3的动力学稳定性,采用有限位移方法(对立方相Pm3m,计算采用80个原子超胞,主要是因为立方相计算的能带结构为金属态)或线性响应方法(对Pnma斜方相)进行晶格动力学计算。2. 研究结果
2.1 NaPO3高压相变序列
NaPO3的常压β相(P21/n)、透辉石相(C2/c)、钛铁矿相(R
¯3 )以及钙钛矿相(斜方相Pnma和立方相Pm3m)分别在0~10 GPa、0~30 GPa、0~80 GPa压力范围内进行几何优化和总能计算。根据热力学定义,G = U − TS+pV = H – TS(G为吉布斯自由能,U为体系内能,T为绝对温度,S为熵,p为压强,V为体积,H为焓)。由于在密度泛函计算中T = 0 K,这种情况下吉布斯自由能G与焓H相同,因此,可以根据不同结构间的焓差(ΔH)来判断不同压力条件下热力学稳定相和相变压力。图2为计算获得的NaPO3的5种不同结构的焓差与压力的关系(焓差以钛铁矿结构为基准)。可以看到,在0 GPa压力下,常压结构P21/n的能量最低,与实验观察一致。P21/n相在约2 GPa下转变成C2/c相。随着压力的增加,C2/c相在约20 GPa下转变成R¯3 相,PO4四面体转变成PO6八面体。随后,R¯3 相在约50 GPa下转变成斜方钙钛矿Pnma相。在整个计算的压力范围内(0~80 GPa),立方钙钛矿Pm3m相的能量总是高于斜方钙钛矿Pnma相,因此,Pm3m相不可能是热力学稳定相。表1为几何优化后得到的0 GPa下NaPO3的5种结构的晶格参数、晶胞体积和原子坐标,其中:V为晶胞体积,V/Z为单位化学式体积。可以看到,计算获得的P21/n相的晶格常数比实验值略大,计算与实验的原子分数坐标基本相当,计算获得的体积比实验值大约2.26%[19],表明本研究的计算方案能较好地模拟实验结果。由于NaPO3的其他几种结构没有实验数据和研究报道,因此,本研究将其与MgSiO3的对应结构进行比较。C2/c-MgSiO3只能在高压下稳定,卸压过程中发生相变[30];在7.93 GPa下,MgSiO3的V = 381.50 Å3 [27]。C2/c-NaPO3在8 GPa下的V = 431.53 Å3,显著大于C2/c-MgSiO3的晶胞体积。R
¯3 -NaPO3在常压下的晶胞体积V0 = 293.54 Å3 ,c/a=3.446,显著大于R¯3 -MgSiO3在常压下的实验值(V0 = 262.54 Å3,c/a = 2.868)[28]。然而,当NaPO3和MgSiO3均为斜方钙钛矿Pnma结构时,计算获得的NaPO3的常压体积V0 = 162.42 Å3,较MgSiO3在常压下的实验值(V0 = 162.80 Å3)[29]略小。由于Na+离子半径比Mg2+离子半径大,而P5+离子半径比Si4+小[31],通过上述3种结构的晶胞体积对比可知,C2/c-NaPO3和R¯3 -NaPO3中Na+离子对其晶胞体积的贡献起主导作用,而Pnma-NaPO3中P5+离子起主导作用。表 1 0 GPa下NaPO3的5种结构在几何优化后的结构参数Table 1. Structural parameters of five phases for NaPO3 at 0 GPaPhase VZ−1/(Å3·formula−1) a/Å b/Å c/Å β/(°) Atom Wyckoff site x y z P21/n 66.460 11.445 6.074 7.669 85.91 Na1 4e 0.1949 0.9249 0.6067 Na2 4e − 0.0001 0.3295 0.8093 P1 4e 0.1918 0.4206 0.4687 P2 4e 0.0963 0.1538 0.1823 O1 4e 0.1914 0.3465 0.6561 O2 4e 0.1559 0.1988 0.3661 O3 4e 0.1243 0.6207 0.4222 O4 4e 0.1723 0.9435 0.1085 O5 4e 0.1136 0.3513 0.0643 O6 4e 0.0236 0.9388 0.7763 C2/c 60.831 10.193 10.040 4.804 98.15 Na1 4e 0.5000 0.3967 0.2500 Na2 4e 0.5000 0.2112 0.7500 P 8f 0.2825 0.0880 0.2236 O1 8f 0.1338 0.0912 0.1854 O2 8f 0.3631 0.2083 0.3163 O3 8f 0.3358 0.0361 − 0.0639 R¯3 48.990 4.618 15.915 Na 6c 0 0 0.3556 P 6c 0 0 0.1646 O 18f 0.3308 0.0527 0.2278 Pm3m 40.691 3.440 Na 1a 0 0 0 P 1b 0.5000 0.5000 0.5000 O 3c 0.5000 0.5000 0 Pnma 40.604 4.848 6.885 4.866 Na 4c 0.0240 0.2500 0.0042 P 4b 0.5000 0 0 O1 4c − 0.0053 0.2500 0.5567 O2 8d 0.2810 0.0285 0.2848 图3显示了计算获得的NaPO3的5种结构的晶胞体积-压力关系。为便于比较,纵坐标为单位化学式体积。根据相变压力下两相的晶胞体积来计算相变的体积变化,可以得到,P21/n相转变成C2/c相时(2 GPa),体积缩小约7.1%;C2/c相转变成R
¯3 相时(20 GPa),体积缩小约11.5%;R¯3 相转变成Pnma钙钛矿相时(50 GPa),体积缩小约9.0%。根据各相的p-V数据,利用三阶Birch-Murnaghan方程[30]拟合其状态方程参数,结果列于表2,其中:K0为常压下的体积模量,K′0 为体积模量K0对压力的导数。表 2 拟合得到的NaPO3的5种结构的状态方程参数Table 2. Fitted parameters of equation of state for the five structures of NaPO3Phase V0/Å3 K0/GPa K′0 P21/n 531.84(32) 31.7(5) 5.16(13) C2/c 486.07(32) 47.8(5) 5.78(7) R¯3 293.35(12) 116.7(7) 4.81(3) Pnma 162.34(4) 215.7(1.1) 4.55(4) Pm3m 40.67(1) 235.6(1.2) 4.58(4) 2.2 钙钛矿结构NaPO3的晶格动力学
ABX3型钙钛矿的各种畸变结构可通过理想的立方钙钛矿结构中BX6八面体的旋转、畸变,以及A、B阳离子的位移得到[32]。立方钙钛矿结构晶胞中,A阳离子占据立方体的角顶,B阳离子占据立方体的中心,X阴离子占据立方体6个面的中心位置,BX6八面体通过共角顶的X阴离子连接形成三维结构,如图4(a)所示。Glazer[33]基于群论对立方钙钛矿结构的23种刚性BX6八面体进行旋转组合,得到15种钙钛矿结构,图4(b)为斜方钙钛矿Pnma的晶体结构示意图。这些旋转畸变可以与立方钙钛矿布里渊区边界R点和M点的声子不稳联系起来[34]。其他一些电子不稳引起的畸变钙钛矿结构可以在Glazer工作的基础上通过BX6八面体畸变(如由Jahn-Teller效应引起)和A、B阳离子位移(如由铁电相变引起的)衍生出来[31],这些变化也可通过立方钙钛矿结构的声子不稳定表现出来。因此,钙钛矿结构的晶格动力学是研究钙钛矿结构动力学稳定性的有效手段。
图5(a)是采用有限位移方法计算得到的0 GPa下的Pm3m-NaPO3(a0 =
3.4462 Å)的声子色散曲线。从图5(a)中可以看到,布里渊区中心Γ点的所有光学振动模均为实频,表明阳离子不发生位移,不出现铁电畸变。布里渊区边界R点和M点出现显著且相近的虚频,分别为313i cm−1 (R+4 )和310i cm−1 (M+3 )。一般地,R+4 不稳指示PO6八面体反向旋转,而M+3 不稳指示PO6八面体同向旋转,两者组合可使Pm3m相畸变成Pnma结构(Glazer表示为a−b+a−)。ABX3型钙钛矿结构的这种八面体旋转模式使A离子配位数从12减小到8,以适应A离子半径较小的情况[32–33]。图5(b)为Pnma-NaPO3在常压时的声子色散曲线,可见,其整个布里渊区的声子模都为实频,说明其动力学稳定。基于离子晶体球体堆积几何模型的Goldschmidt容忍因子t(t=(rA+rX)/
√2 (rB+rX),其中,rA、rB、rX分别为A、B、X的离子半径)被广泛用于预测ABX3型钙钛矿结构的稳定性和可能的畸变结构[35]。通常来说,t = 1时,钙钛矿为立方相;t >1时,钙钛矿中出现BX6八面体共面的六方钙钛矿结构;t < 1时,钙钛矿为低对称结构,如常见的Pnma结构和R¯3 /c结构等。NaPO3的容忍因子t = 1.025(Na+离子半径为1.18 Å,配位数为 8;P5+离子半径为0.38 Å,配位数为6)[31]),大于MgSiO3的容忍因子(t = 0.900)和CaSiO3的容忍因子(t = 0.990)。钙钛矿MgSiO3和CaSiO3的基态结构分别为Pnma和I4/mcm[29, 36]。Na―O键是典型的离子键,因此,NaPO3采取Pnma钙钛矿结构,而非立方钙钛矿结构Pm3m,很可能与P―O键的强共价键特性有关。2.3 NaPO3斜方钙钛矿的结构特征
表3列出了NaPO3斜方钙钛矿在 0 GPa下的键长和键角数据。PO6八面体的键长范围较窄(1.729~1.757 Å),O―P―O键角略小于90°,表明结构中PO6八面体的畸变较小。Na―O的键长范围为2.181~2.396 Å(当Na的配位多面体为NaO8时)或2.181~2.800 Å(当Na的配位多面体为NaO12时),这些键长比Na+(配位数为8时,Na+的半径为1.18 Å;配位数为12时,Na+的半径为1.39 Å)与O2−(半径为1.40 Å)的半径之和(2.58 Å或2.79 Å[31])小,反映了NaPO3斜方结构中Na―O键处于压缩状态。根据键长可以计算阳离子的键价总数,并用于评估其键合状态[37],计算获得P的键价总数为4.27,小于其形式价(+5),表明P―O键处于欠键(拉伸)状态;Na的键价总数为2.15(对于NaO8)或2.49(对于NaO12),显著大于其形式价(+1),表明Na―O键处于过键(压缩)状态。在NaPO3斜方钙钛矿中,Na―O键和P―O键的这种状态与通常认为的A阳离子半径偏小进而通过BX6八面体倾斜形成斜方钙钛矿结构以匹配A―X键长和B―X键长的普遍规则不一致,表明斜方结构中存在其他因素抵消了其应变能,导致斜方结构比立方结构具有更低的能量。
表 3 0 GPa下Pnma-NaPO3的键长和键角Table 3. Bond distances and bond angles of NaPO3 with orthorhombic perovskite structure (Pnma) at 0 GPaStructure Type Bond distance/Å Bond angle/(º) PO6 P―O(1) [×2] 1.729 P―O(2) [×2] 1.744 P―O(2) [×2] 1.757 P―O1―P 161.62 P―O2―P [×2] 160.23 O1―P―O2 89.76 89.65 O2―P―O2 89.18 NaO12 Na―O(1) 2.181 Na―O(1) 2.301 Na―O(2) [×2] 2.184 Na―O(2) [×2] 2.390 Na―O(2) [×2] 2.396 Na―O(1) 2.583 Na―O(1) 2.693 Na―O(2) [×2] 2.800 磷因其高的价态、小的离子半径和较高的电负性易与其他阴离子形成较强的共价键,如常压下磷通过sp3杂化与氧形成PO4四面体,其P―O的平均键长约为1.54 Å[1],小于SiO4四面体中Si―O平均键长(约为1.63 Å)。常压下,NaPO3斜方钙钛矿中PO6八面体的P―O的平均键长约为1.743 Å,小于MgSiO3斜方钙钛矿中SiO6八面体的Si―O的平均键长(1.793 Å[29])。主族元素在BX6八面体中的B―X成键机制与过渡族元素不同,参与成键主要是s和p轨道,而过渡族元素参与成键主要是d轨道。这种成键机制的差异导致主族元素形成的B―X―B键偏离180°,从而形成多重σ键,并呈现出可以获得更大键能的趋势[38–39]。图6为计算获得的常压下NaPO3斜方钙钛矿的总态密度和分波态密度,其直接带隙Eg = 1.59 eV。从图6中可以看出,在−15~−5 eV范围内,P的3s、3p轨道与O的2p轨道强烈混合,反映出P―O键的强共价性。由此看来,NaPO3采取斜方钙钛矿结构而非立方结构是PO6八面体中P―O键的强共价性所致。
图7显示了NaPO3斜方钙钛矿结构参数随压力的变化。为更好地比较,晶格常数(a、b、c)转换为赝立方结构的格子常数(ap、bp、cp),其转换关系为:ap =
a/√2 ,bp = b/2,cp =c/√2 。可以看出,在计算的压力范围内,ap、bp、cp相差很小,且随压力升高,相对大小发生改变。对于大多数斜方钙钛矿结构化合物,如MgSiO3、CaTiO3、GdFeO3等,晶格参数受八面体倾斜模式(a−b+a−)所控,一般地表现为ap > bp > cp[29, 40–43]。NaPO3斜方钙钛矿的晶格参数的相对大小与普遍规律不同,显示出“异常”的结构特征,反映了P―O键的强共价性对其结构的重要影响。从图7(b)所示的P―O的键长随压力的变化关系可以看出,PO6八面体中,3种P―O键长相差较小,且随压力增加,键长差异越来越小。在60 GPa下,PO6八面体的P―O键的平均键长为1.648 Å,与文献[11–17]报道的实验观测和理论计算结果相当。图7(c)显示了结构中P―O1―P的键角和P―O2―P的键角随压力的变化情况。常压下P―O1―P的键角为161.62°,P―O2―P的键角为160.23°,大于MgSiO3斜方钙钛矿中相应的键角(146.54°和147.02°),表明NaPO3斜方钙钛矿中PO6八面体的倾斜较小,这一特征是Na+的大离子半径和P―O键的强共价性共同作用的结果。随着压力增加,P―O1―P的键角和P―O2―P的键角逐渐变小,且两者差异变大。图7(d)显示NaO12 与PO6 多面体体积比(VNaO12 /VPO6 )[44–45]随压力的增加而逐渐减小,表明NaO12多面体比PO6八面体更易压缩。与之对照,常压下,MgSiO3斜方钙钛矿MgO12 和SiO6 [29]多面体体积比VMgO12 /VSiO6 = 4.28[29],显著小于NaPO3斜方钙钛矿多面体体积比(4.75)。3. 结 论
第一性原理密度泛函理论计算发现,NaPO3与等电子体的MgSiO3有相似的高压相变序列。随着压力的增加,NaPO3发生如下相变:P21/n → C2/c (2 GPa) → R
¯3 (20 GPa)→ Pnma(50 GPa),相变时体积收缩分别为7.1%、11.5%和9.0%。NaPO3从透辉石型结构(C2/c)转变成钛铁矿型(R¯3 )结构时,磷从PO4四面体转变成PO6八面体,其转变压力(约20 GPa)与P2O5的相应转变压力相当。与其他结构不同,斜方钙钛矿结构NaPO3的晶胞体积比斜方钙钛矿结构MgSiO3的晶胞体积略小。晶格动力学研究显示,NaPO3的立方相在布里渊区边界R点和M点出现显著的虚频,表明立方相的动力学不稳定,而斜方相则在整个布里渊区显示实频,表明斜方相的动力学稳定。NaPO3斜方钙钛矿结构显示出的“异常”晶体化学特征与PO6八面体中P―O键的强共价性相关。 -
表 1 0 GPa下NaPO3的5种结构在几何优化后的结构参数
Table 1. Structural parameters of five phases for NaPO3 at 0 GPa
Phase VZ−1/(Å3·formula−1) a/Å b/Å c/Å β/(°) Atom Wyckoff site x y z P21/n 66.460 11.445 6.074 7.669 85.91 Na1 4e 0.1949 0.9249 0.6067 Na2 4e − 0.0001 0.3295 0.8093 P1 4e 0.1918 0.4206 0.4687 P2 4e 0.0963 0.1538 0.1823 O1 4e 0.1914 0.3465 0.6561 O2 4e 0.1559 0.1988 0.3661 O3 4e 0.1243 0.6207 0.4222 O4 4e 0.1723 0.9435 0.1085 O5 4e 0.1136 0.3513 0.0643 O6 4e 0.0236 0.9388 0.7763 C2/c 60.831 10.193 10.040 4.804 98.15 Na1 4e 0.5000 0.3967 0.2500 Na2 4e 0.5000 0.2112 0.7500 P 8f 0.2825 0.0880 0.2236 O1 8f 0.1338 0.0912 0.1854 O2 8f 0.3631 0.2083 0.3163 O3 8f 0.3358 0.0361 − 0.0639 R¯3 48.990 4.618 15.915 Na 6c 0 0 0.3556 P 6c 0 0 0.1646 O 18f 0.3308 0.0527 0.2278 Pm3m 40.691 3.440 Na 1a 0 0 0 P 1b 0.5000 0.5000 0.5000 O 3c 0.5000 0.5000 0 Pnma 40.604 4.848 6.885 4.866 Na 4c 0.0240 0.2500 0.0042 P 4b 0.5000 0 0 O1 4c − 0.0053 0.2500 0.5567 O2 8d 0.2810 0.0285 0.2848 表 2 拟合得到的NaPO3的5种结构的状态方程参数
Table 2. Fitted parameters of equation of state for the five structures of NaPO3
Phase V0/Å3 K0/GPa K′0 P21/n 531.84(32) 31.7(5) 5.16(13) C2/c 486.07(32) 47.8(5) 5.78(7) R¯3 293.35(12) 116.7(7) 4.81(3) Pnma 162.34(4) 215.7(1.1) 4.55(4) Pm3m 40.67(1) 235.6(1.2) 4.58(4) 表 3 0 GPa下Pnma-NaPO3的键长和键角
Table 3. Bond distances and bond angles of NaPO3 with orthorhombic perovskite structure (Pnma) at 0 GPa
Structure Type Bond distance/Å Bond angle/(º) PO6 P―O(1) [×2] 1.729 P―O(2) [×2] 1.744 P―O(2) [×2] 1.757 P―O1―P 161.62 P―O2―P [×2] 160.23 O1―P―O2 89.76 89.65 O2―P―O2 89.18 NaO12 Na―O(1) 2.181 Na―O(1) 2.301 Na―O(2) [×2] 2.184 Na―O(2) [×2] 2.390 Na―O(2) [×2] 2.396 Na―O(1) 2.583 Na―O(1) 2.693 Na―O(2) [×2] 2.800 -
[1] HUMINICKI D M C, HAWTHORNE F C. The crystal chemistry of the phosphate minerals [J]. Reviews in Mineralogy and Geochemistry, 2002, 48(1): 123–253. doi: 10.2138/rmg.2002.48.5 [2] BRUNET F, CHAZOT G. Partitioning of phosphorus between olivine, clinopyroxene and silicate glass in a spinel lherzolite xenolith from Yemen [J]. Chemical Geology, 2001, 176(1/2/3/4): 51–72. [3] WALTON C R, SHORTTLE O, JENNER F E, et al. Phosphorus mineral evolution and prebiotic chemistry: from minerals to microbes [J]. Earth-Science Reviews, 2021, 221: 103806. doi: 10.1016/j.earscirev.2021.103806 [4] MCDONOUGH W F, SUN S S. The composition of the Earth [J]. Chemical Geology, 1995, 120(3/4): 223–253. [5] BRUNET F, BONNEAU V, IRIFUNE T. Complete solid-solution between Na3Al2(PO4)3 and Mg3Al2(SiO4)3 garnets at high pressure [J]. American Mineralogist, 2006, 91(1): 211–215. doi: 10.2138/am.2006.2053 [6] WATSON E B, CHERNIAK D J, HOLYCROSS M E. Diffusion of phosphorus in olivine and molten basalt [J]. American Mineralogist, 2015, 100(10): 2053–2065. doi: 10.2138/am-2015-5416 [7] XING C M, WANG C Y, TAN W. Disequilibrium growth of olivine in mafic magmas revealed by phosphorus zoning patterns of olivine from mafic-ultramafic intrusions [J]. Earth and Planetary Science Letters, 2017, 479: 108–119. doi: 10.1016/j.jpgl.2017.09.005 [8] KONZETT J, FROST D J. The high P-T stability of hydroxyl-apatite in natural and simplified MORB—an experimental study to 15 GPa with implications for transport and storage of phosphorus and halogens in subduction zones [J]. Journal of Petrology, 2009, 50(11): 2043–2062. doi: 10.1093/petrology/egp068 [9] KONZETT J. From phosphates to silicates and back: an experimental study on the transport and storage of phosphorus in eclogites during uplift and exhumation [J]. American Mineralogist, 2016, 101(8): 1756–1768. doi: 10.2138/am-2016-5521 [10] YE K, CONG B L, YE D N. The possible subduction of continental material to depths greater than 200 km [J]. Nature, 2000, 407(6805): 734–736. doi: 10.1038/35037566 [11] PELLICER-PORRES J, SAITTA A M, POLIAN A, et al. Six-fold-coordinated phosphorus by oxygen in AlPO4 quartz homeotype under high pressure [J]. Nature Materials, 2007, 6(9): 698–702. doi: 10.1038/nmat1966 [12] BRUNET F, FLANK A M, ITIÉ J P, et al. Experimental evidence of sixfold oxygen coordination for phosphorus [J]. American Mineralogist, 2007, 92(7): 989–993. doi: 10.2138/am.2007.2570 [13] STEBBINS J F, KIM N, BRUNET F, et al. Confirmation of octahedrally coordinated phosphorus in AlPO4-containing stishovite by 31P NMR [J]. European Journal of Mineralogy, 2009, 21(4): 667–671. doi: 10.1127/0935-1221/2009/0021-1953 [14] PAKHOMOVA A, APRILIS G, BYKOV M, et al. Penta-and hexa-coordinated beryllium and phosphorus in high-pressure modifications of CaBe2P2O8 [J]. Nature Communications, 2019, 10(1): 2800. doi: 10.1038/s41467-019-10589-z [15] SALVADÓ M A, PERTIERRA P. Theoretical study of P2O5 polymorphs at high pressure: hexacoordinated phosphorus [J]. Inorganic Chemistry, 2008, 47(11): 4884–4890. doi: 10.1021/ic8001543 [16] LÓPEZ-SOLANO J, RODRÍGUEZ-HERNÁNDEZ P, MUÑOZ A, et al. Theoretical and experimental study of the structural stability of TbPO4 at high pressures [J]. Physical Review B, 2010, 81(14): 144126. doi: 10.1103/PhysRevB.81.144126 [17] LÓPEZ-MORENO S, ERRANDONEA D. Ab initio prediction of pressure-induced structural phase transitions of CrVO4-type orthophosphates [J]. Physical Review B, 2012, 86(10): 104112. doi: 10.1103/PhysRevB.86.104112 [18] JOST K H. Die struktur des Kurrol’schen Na-salzes (NaPO3) x typ A [J]. Acta Crystallographica, 1961, 14(8): 844–847. doi: 10.1107/S0365110X6100245X [19] JOST K H. Die struktur des Kurrol’schen Na-salzes (NaPO3) x, typ B [J]. Acta Crystallographica, 1963, 16(7): 640–642. doi: 10.1107/S0365110X63001687 [20] ONDIK H M. The structure of anhydrous sodium trimetaphosphate Na3P3O9, and the monohydrate, Na3P3O9·H2O [J]. Acta Crystallographica, 1965, 18(2): 226–232. doi: 10.1107/S0365110X65000518 [21] IMMIRZI A, PORZIO W. A new form of sodium Kurrol salt studied by the Rietveld method from X-ray diffraction data [J]. Acta Crystallographica Section B, 1982, 38(11): 2788–2792. doi: 10.1107/S0567740882009960 [22] THILO E. The structural chemistry of condensed inorganic phosphates [J]. Angewandte Chemie International Edition, 1965, 4(12): 1061–1071. doi: 10.1002/anie.196510611 [23] DURIF A. Crystal chemistry of condensed phosphates [M]. New York: Springer, 1995. [24] KORNBERG A, RAO N N, AULT-RICHÉ D. Inorganic polyphosphate: a molecule of many functions [J]. Annual Review of Biochemistry, 1999, 68: 89–125. doi: 10.1146/annurev.biochem.68.1.89 [25] AKBARI A, WANG Z J, HE P S, et al. Unrevealed roles of polyphosphate-accumulating microorganisms [J]. Microbial Biotechnology, 2021, 14(1): 82–87. doi: 10.1111/1751-7915.13730 [26] GASPARIK T. Phase diagrams for geoscientists: an atlas of the Earth’s interior [M]. 2nd ed. New York: Springer, 2014. [27] ANGEL R J, CHOPELAS A, ROSS N L. Stability of high-density clinoenstatite at upper-mantle pressures [J]. Nature, 1992, 358(6384): 322–324. doi: 10.1038/358322a0 [28] HORIUCHI H, HIRANO M, ITO E, et al. MgSiO3 (ilmenite-type): single crystal X-ray diffraction study [J]. American Mineralogist, 1982, 67(7/8): 788–793. [29] HORIUCHI H, ITO E, WEIDNER D J. Perovskite-type MgSiO3: single-crystal X-ray diffraction study [J]. American Mineralogist, 1987, 72(3/4): 357–360. [30] ANGEL R J, HUGH-JONES D A. Equations of state and thermodynamic properties of enstatite pyroxenes [J]. Journal of Geophysical Research: Solid Earth, 1994, 99(B10): 19777–19783. doi: 10.1029/94JB01750 [31] SHANNON R D. Revised effective ionic radii and systematic studies of interatomic distances in halides and chalcogenides [J]. Acta Crystallographica Section A, 1976, 32(5): 751–767. doi: 10.1107/S0567739476001551 [32] WOODWARD P M. Octahedral tilting in perovskites. Ⅰ. geometrical considerations [J]. Acta Crystallographica Section B, 1997, 53(1): 32–43. doi: 10.1107/S0108768196010713 [33] GLAZER A M. The classification of tilted octahedra in perovskites [J]. Acta Crystallographica Section B, 1972, 28(11): 3384–3392. doi: 10.1107/S0567740872007976 [34] ALEKSANDROV K S. The sequences of structural phase transitions in perovskites [J]. Ferroelectrics, 1976, 14(1): 801–805. doi: 10.1080/00150197608237799 [35] WOODWARD P M. Octahedral tilting in perovskites. Ⅱ. structure stabilizing forces [J]. Acta Crystallographica Section B, 1997, 53(1): 44–66. doi: 10.1107/S0108768196012050 [36] STIXRUDE L, COHEN R E, YU R C, et al. Prediction of phase transition in CaSiO3 perovskite and implications for lower mantle structure [J]. American Mineralogist, 1996, 81(9/10): 1293–1296. [37] BROWN I D, ALTERMATT D. Bond-valence parameters obtained from a systematic analysis of the inorganic crystal structure database [J]. Acta Crystallographica Section B, 1985, 41(4): 244–247. doi: 10.1107/S0108768185002063 [38] MIZOGUCHI H, ENG H W, WOODWARD P M. Probing the electronic structures of ternary perovskite and pyrochlore oxides containing Sn4+ or Sb5+ [J]. Inorganic Chemistry, 2004, 43(5): 1667–1680. doi: 10.1021/ic034551c [39] MIZOGUCHI H, WOODWARD P M, BYEON S H, et al. Polymorphism in NaSbO3: structure and bonding in metal oxides [J]. Journal of the American Chemical Society, 2004, 126(10): 3175–3184. doi: 10.1021/ja038365h [40] SASAKI S, PREWITT C T, BASS J D, et al. Orthorhombic perovskite CaTiO3 and CdTiO3: structure and space group [J]. Acta Crystallographica Section C, 1987, 43(9): 1668–1674. [41] MAREZIO M, REMEIKA J P, DERNIER P D. The crystal chemistry of the rare earth orthoferrites [J]. Acta Crystallographica Section B, 1970, 26(12): 2008–2022. doi: 10.1107/S0567740870005319 [42] ZHOU J S, GOODENOUGH J B. Universal octahedral-site distortion in orthorhombic perovskite oxides [J]. Physical Review Letters, 2005, 94(6): 065501. doi: 10.1103/PhysRevLett.94.065501 [43] MARTÍNEZ-LOPE M J, ALONSO J A, RETUERTO M, et al. Evolution of the crystal structure of RVO3 (R=La, Ce, Pr, Nd, Tb, Ho, Er, Tm, Yb, Lu, Y) perovskites from neutron powder diffraction data [J]. Inorganic Chemistry, 2008, 47(7): 2634–2640. doi: 10.1021/ic701969q [44] THOMAS N W. A new global parameterization of perovskite structures [J]. Acta Crystallographica Section B, 1998, 54(5): 585–599. doi: 10.1107/S0108768198001979 [45] AVDEEV M, CASPI E N, YAKOVLEV S. On the polyhedral volume ratios VA/VB in perovskites ABX3 [J]. Acta Crystallographica Section B, 2007, 63(3): 363–372. doi: 10.1107/S0108768107001140 -