-
双折射晶体是制作偏振器、Q开关、光隔离器等的重要材料[1–8]. 近年来已经发现了许多较好的双折射晶体, 如CaCO3[9], YVO4[10], α-BBO[11–13]等, 其中YVO4以其突出的特点, 如低成本、高温稳定性、在较宽波长范围(400—5000 nm)内具有较高的光学透过率和良好的物理性能. 虽然许多双折射晶体已经被广泛研究[9–16], 但仍然缺乏能够满足社会发展需求的优秀晶体, 如经济实惠、出色的温度稳定性、在宽波长范围(400—5000 nm)内的高传输速率以及良好的物理特性. 因此, 寻找新的双折射晶体越来越受到关注, 受到YVO4大双折射(Δn = 0.21@532 nm)的启发, 期望能阐述钒酸盐具有良好双折射率的内在机理.
钒酸盐晶体作为光学晶体家族中的重要成员具有许多独特的特性, 在许多对钒酸盐中掺杂稀土元素研究中[17–19], 特别是掺杂Nd3+的钒酸盐GdVO4, YVO4等晶体是设计光谱近红外范围激光源的优异介质, 同时也适用于基于受激拉曼散射(SRS)技术的多种激光器[20]. Lei等[21]利用响应电子分布各向异性(REDA)方法得出[VO4]基团对YVO4的双折射的贡献约为70%; Huang等[22]通过对VOF多面体-碘酸盐体系的探索, 成功地合成了Sr[VO2F(IO3)2]和Sr3F2(VO2F4)(IO3)两种在550 nm处的双折射率分别为0.250和0.406的材料, 显著高于市售的钒酸盐(YVO4, 在532 nm处为0.204)的双折射率; Zhang等[23]系统研究了关于硫代钒酸盐中[VOxS4–x](x = 0—4)四面体对其双折射率的增强作用; 程金磊等[24]合成了具有[VO3]∞链的Li0.8Na0.2CsV2O6·H2O晶体, 并系统研究了[VO3]∞对双折射率的增大效应.
碱金属钒酸盐AV3O8 (A = Li, Na, K, Rb)晶体也拥有独特的[VO]基团: LiV3O8晶体由3个[VO6]的八面体阴离子基团围绕在Li离子周围; NaV3O8晶体中[V3O12]基团与Na原子交叉排列形成一维无限链; KV3O8和RbV3O8晶体中[V3O12]基团通过O原子相连接无限延展. 但目前对碱金属钒酸盐AV3O8 (A = Li, Na, K, Rb)的研究和应用还主要集中在用于提高电池正极材料的电化学性能[25–33]. 在对其光学特性的研究中, 通常使用XRD (X射线衍射)来测试碱金属钒酸盐的光化学稳定性, 但在双折射领域的研究相对较少. 因此本文通过第一性原理的计算碱金属钒酸盐AV3O8 (A = Li, Na, K, Rb)晶体结构与光学性质之间的关系, 系统研究AV3O8 (A = Li, Na, K, Rb)晶体电子结构及双折射率, 从而理清碱金属钒酸盐AV3O8 (A = Li, Na, K, Rb)具有良好双折射率的理论机理.
-
通过调研The Material Project晶体数据库, 本文的晶体结构选择LiV3O8 (mp-27503)[34], NaV3O8 (mp-755861)[35], KV3O8 (mp-19172)[36]和 RbV3O8 (mp-16864)[37]作为研究对象. 基于密度泛函理论中的CASTEP软件包[38], 选取平面基截止能(cutoff energy)迭代收敛精度(SCF)均为5×10–7 eV/atom, 能量截止为830 eV, 使用Monkhorst-Pack方法设置K点密度为0.005 nm–1[39–41], 其中电子交换关联能采用广义梯度近似GGA-PBE形式的泛函[42]处理, 离子实与价电子的相互作用势使用模守恒赝势(NCP)[43]确定. 体系电子波函数通过平面波基组展开, 选取的价层电子组态分别为 Li: 2s1, Na: 2p63s1, K: 3p64s1, Rb: 4p65s1, V: 3p63d34s2, O: 2s22p4, 并对所选参数均进行收敛性测试以保证计算结果的准确性.
LiV3O8, NaV3O8, KV3O8和RbV3O8四种晶体均属于单斜晶系P21/m空间群, 其β角分别为107.830°, 101.993°, 98.1353°, 96.943°. 图1为4种晶体结构图, 可以看出NaV3O8, KV3O8和RbV3O8三种晶体结构中的一个显著特征是以[V3O12]为基本的构建模块. NaV3O8通过碱金属原子Na相互连接, KV3O8和RbV3O8通过O原子相互连接, 形成沿b轴为方向的一维无限[V3O12]∞链; LiV3O8晶体由3个[VO6]的八面体阴离子基团围绕在Li离子周围组成的[LiV3O21]为基本结构单元, 沿b轴为方向的一维无限[LiV3O21]∞链. 这种[V3O12]和[VO6]阴离子基团的原子排列方式更容易获得较大的光学各向异性[44].
-
为了保证计算结果的准确性, 基于4种晶体结构的实验数据进行几何结构优化, 在此基础上计算其电子结构与光学性质, 表1给出了几何结构优化前后AV3O8 (A = Li, Na, K, Rb)的晶格参数和晶胞体积. 通过对比, AV3O8 (A = Li, Na, K, Rb)晶胞参数的最大误差仍小于5%, 说明本工作计算参数的可行性.
-
能带结构描述了固体材料中电子能量状态分布, 其主要用来解释和预测材料的电子性质、导电性、磁性和光学性质[45]. 图2为AV3O8 (A = Li, Na, K, Rb)晶体的能带结构图, LiV3O8, NaV3O8, KV3O8和RbV3O8的禁带宽度分别为1.695, 1.898, 1.965和1.984 eV, 属于半导体材料. 结合图3可知, 在2—4 eV范围内的能带主要来源于O原子的s和p轨道电子贡献. 价带则主要分布在–5 eV到费米能级的范围内, 在这一范围内的能带以O-2p态及V-3d态杂化贡献为主. 从能带图中还可以看到LiV3O8, NaV3O8, KV3O8和RbV3O8晶体的带隙类型均为间接带隙.
对态密度的分析有助于理解从占据态到非占据态的跃迁, 并有助于研究其光学性质. 为了进一步研究不同轨道对电子态的贡献以及成键特性, 本文对4种晶体的态密度(DOS)和分波态密度(PDOS)进行分析, 使用Gaussian 09计算软件[46]计算LiV3O8晶体的阴离子基团[VO4]3–和NaV3O8, KV3O8, RbV3O8三种晶体的阴离子基团[V3O8]–的最高占据轨道(highest occupied molecular orbital, HOMO)和最低未占据分子轨道(lowest unoccupied molecular orbital, LUMO).
图3为AV3O8 (A = Li, Na, K, Rb) 4种晶体中不同原子的电子轨道对能带的贡献及HOMO&LUMO图. 对图3中4种晶体的态密度分析可知: LiV3O8晶体在–20—–15 eV的价带由O的2s轨道所占据; NaV3O8晶体在–23—–20 eV的价带由Na的2p轨道所占据; KV3O8晶体在–12—–10 eV的价带被K的3p轨道所占据; RbV3O8晶体在–9—–6 eV的价带被Rb的4p轨道所占据, NaV3O8, KV3O8, RbV3O8三种晶体在–18—–15 eV区域的价带均由O的2s轨道所占据. 4种晶体的价带顶部和导带底部的贡献均主要由V-3d轨道和O-2p轨道共同提供. 并结合图3中4种晶体的HOMO&LUMO图像进一步说明LiV3O8晶体的阴离子基团[VO4]3–和AV3O8 (A = Na, K, Rb) 3种晶体的阴离子基团[V3O8]–由O-p轨道与V-d轨道组成d-p杂化, 在HOMO中电子在O原子周围聚集, LUMO中O原子和V原子周围存在大量空轨道. 4种化合物在费米能级附近d-p能级交织更加突出, 充分表明V-3d和O-2p之间的轨道杂化得到显著增强. 在[VO]基团中V-3d轨道与O-2p轨道形成σ键和π键, 增强了[VO]基团的稳定性.
由于碱金属阳离子半径的增大会使晶格参数膨胀削弱V—O的d-p轨道杂化, 同时由于离子层间距离扩张导致碱金属阳离子的p轨道局域化增强[47], 使4种晶体的带隙宽度随碱金属阳离子的半径增大而增大, 且碱金属阳离子的p轨道电子对态密度的贡献逐渐靠近费米能级.
综上可知, 在费米能级附近, AV3O8 (A = Li, Na, K, Rb) 4种晶体的价带和导带由O的2p轨道和V的3d轨道主导, 并且表现出从–5—5 eV的d-p轨道杂化, 结合HOMO&LUMO和分波态密度进一步证实了4种碱金属钒酸盐化合物导带底与价带顶的轨道来源: LiV3O8晶体的阴离子基团[VO4]3–和AV3O8 (A = Na, K, Rb) 3种晶体的阴离子基团[V3O8]–的O-p轨道与V-d轨道组成d-p杂化. 由此可知, [V3O12]和[VO6]阴离子基团是该化合物带隙的主要贡者.
-
Mulliken布居数分析可以定量地分析所形成的化学键, 原子布居数可以反映出各轨道的占据情况和电子的电荷量. 通常情况下, 键布居数值越接近1, 呈现的共价性越强, 键布居数值越接近0, 呈现的离子性越强[48,49]. AV3O8 (A = Li, Na, K, Rb)的原子布居数与键布居数分析列于表2.
对比优化前的价电子构型, 碱金属原子与V原子为电子供体, O原子为电子受体. 在这4种晶体中V—O的键布居数为0.01和0.97, 说明V, O原子之间形成了共价键和离子键, 并不是唯一的键型, 主要以共价键为主, 而离子键形成的原因是相聚较远的V, O原子的相互作用. 碱金属原子与O原子成键的键布居数值均较小, 表现出离子性.
电子局域函数(electron localization function, ELF)常用来分析共价系统中原子之间的成键情况, 可以给出结构模型中电子的局域化程度, 并反映两原子之间所成化学键的相对强度[50]. AV3O8 (A = Li, Na, K, Rb)的电子局域函数如图4所示, 电子云主要聚集在O原子周围, 且O原子与V原子之间具有方向性, 说明O-P之间呈现共价性; 而碱金属金属阳离子与O之间无方向性, 说明金属阳离子与O之间呈现离子键特性. 所得结论与键布居数所得结果相互验证.
-
材料的双折射率被广泛运用于偏振光学器件以及偏光技术中, 同时也是非线性光学材料一个重要的参数[51–54]. 钒酸盐晶体PbVTeO5F被报道在590 nm处的双折射率为0.183, Pb2TeV2O10在532 nm处的双折射率为0.275[55], Pb2(V2O4F)(VO2)(SeO3)3在1064 nm处的双折射率为0.105[56], Sr[VO2F(IO3)2]和Sr3F2(VO2F4)(IO3)在550 nm处的双折射率分别为0.250和0.406[22]. 进一步说明了含[VO]基团化合物具有良好的双折射效应. 本文研究的4个晶体均为双轴晶体, 双轴晶体的3个主介电系数都不相等, 即ε1 ≠ ε2 ≠ ε3. 通常主介电系数按ε1 < ε2 < ε3取值有两个光轴, 当光沿该光轴方向传播时, 其相应的二特许线偏振光波的传播速度(或折射率)相等. 图5为双折射原理示意图.
图6为LiV3O8, NaV3O8, KV3O8和RbV3O8的双折射率图谱, 在1064 nm处LiV3O8, NaV3O8, KV3O8和RbV3O8的双折射率分别为0.28, 0.30, 0.28和0.27, 这种大的双折射率是由于其结构特点导致的: 每个[LiV3O21]∞和[V3O12]∞层都有其光学各向异性, 而且这些[LiV3O21]∞和[V3O12]∞层是平行或反平行排列的, 结合对HOMO和LUMO的分析可知每一个[VO]基团中的V-3d轨道和O-2p轨道均形成了d-p轨道杂化, 这种特殊的[VO]基团排列结构造成了光学各向异性和大的双折射率的叠加效应. 为了进一步探讨阴离子基团的光学各向异性响应机理, 本文计算了各晶体的响应电子分布各向异性指数(REDA)[21,57], 由表3可知NaV3O8, KV3O8和RbV3O8三种晶体的双折射随REDA指数的变化正相关, 但LiV3O8晶体却表现出独特的双折射关系, 探究其原因是阴离子基团的差异性造成的. 进一步验证了阴离子基团对双折射率的影响.
双折射率还与晶体的各向异性有关[58–60], 因此采用杨氏模量图对比, 4种钒酸盐材料的各向异性. 图7为4种晶体在三维空间中的杨氏模量分布图, 可以看出, AV3O8 (A = Li, Na, K, Rb)结构的杨氏模量E的三维表面结构各向异性特性较为明显, 计算出LiV3O8的杨氏模量E的最小值(44.29 GPa)是沿 [100] 方向, 而最大值(178.65 GP)是沿[010]方向. 得到LiV3O8弹性各向异性[61] AE
$ = (|{E}_{{\mathrm{max}}} | + |{E}_{{\mathrm{min}}} |)/(|{E}_{{\mathrm{max}}}| + |{E}_{{\mathrm{min}}}|) $ 约为60.27%. NaV3 O8, KV3O8和RbV3O8三种晶体的弹性各向异性分别为68.48%, 53.35%, 39.58%. 其中NaV3O8晶体的弹性各向异性最大, 对应的双折射率也大于其他3种晶体; RbV3O8晶体的弹性各向异性最小, 对应的双折射率也最小. 进一步证实了双折射率与晶体各向异性呈正相关[62–64].结合碱金属原子半径分析可知, NaV3O8, KV3O8和RbV3O8三种晶体的双折射率随着碱金属原子半径的增大而减小, 而LiV3O8晶体的双折射率却小于NaV3O8晶体的双折射率, 其可能原因是其独特的3个[VO6]的八面体阴离子基团围绕在Li离子周围, 不满足同类化合物中双折射率随着碱金属原子半径增大而减小的规律.
-
基于密度泛函理论的第一性原理, 对AV3O8 (A = Li, Na, K, Rb)的电子结构和光学性质进行了理论计算. 研究结果表明这4种晶体均为间接带隙, 带隙宽度分别为1.695, 1.898, 1.965, 1.984 eV. 导带底主要为V原子的最外层轨道所占据, 价带顶的主要贡献者是O-2p轨道, O原子的2p轨道在费米能级附近表现出较强的局域性, 在费米能级附近O的2p轨道和V的3d轨道形成p-d轨道杂化. 通过Mulliken布居数和电子局域函数分析四种化合物的O—V键以共价性为主.
4种晶体LiV3O8, NaV3O8, KV3O8和RbV3O8在1064 nm处的双折射率分别为0.28, 0.30, 0.28和0.27. 这与晶体阴离子基团的排列方式、d-p轨道杂化、响应电子分布各向异性指数和晶体各向异性相关. 通过研究晶体双折射率的内在机理, 为开发具有大双折射率的钒酸盐晶体材料提供了重要的理论依据和优化方向.
基于第一性原理对碱金属钒酸盐AV3O8 (A = Li, Na, K, Rb)大双折射率的机理研究
First-principles study of mechanism of high birefringence in alkali metal vanadates AV3O8 (A = Li, Na, K, Rb)
-
摘要: 双折射作为光学晶体的基本参数, 在相位调制、分光、偏振等许多光学应用中发挥着重要的作用, 是激光科学与技术中的关键材料, 而钒酸盐多面体较大双折射率为开发双折射材料提供了一条新的途径. 本文采用第一性原理研究4种碱金属钒酸盐AV3O8 (A = Li, Na, K, Rb)晶体的能带结构、态密度、电子局域函数和双折射率. 计算结果表明碱金属钒酸盐AV3O8 (A = Li, Na, K, Rb)均为间接带隙, 带隙值分别为1.695, 1.898, 1.965和1.984 eV. 对态密度分析可知在费米能级附近, 碱金属钒酸盐AV3O8 (A = Li, Na, K, Rb)导带底主要被V原子的最外层轨道所占据, 价带顶的主要贡献者是O-2p轨道, O原子的2p轨道还在费米能级附近表现出较强的局域性, 结合HOMO和LUMO以及布居分析说明在4种晶体中主要由V-3p轨道与O的2p轨道成键, V—O表现为强的共价键. 通过对晶体结构与光学性质关系的分析, 晶体较大的各向异性, 较高水平的响应电子分布各向异性指数, 阴离子基团的特殊排列和V-3d和O-2p轨道形成的d-p轨道杂化都是导致其大双折射率的主要原因, 经计算所得LiV3O8, NaV3O8, KV3O8和RbV3O8在1064 nm处的双折射率分别为0.28, 0.30, 0.28和0.27.Abstract: Birefringence, as a fundamental parameter of optical crystals, plays a vital role in numerous optical applications such as phase modulation, light splitting, and polarization, thereby making them key materials in laser science and technology. The significant birefringence of vanadate polyhedra provides a new approach for developing birefringent materials. In this study, first-principles calculations are used to investigate the band structures, density of states (DOS), electron localization functions (ELFs), and birefringence behaviors of four alkali metal vanadate crystals AV3O8 (A = Li, Na, K, Rb). The computational results show that all AV3O8 crystals have indirect band gaps, whose values are 1.695, 1.898, 1.965, and 1.984 eV for LiV3O8, NaV3O8, KV3O8, and RbV3O8, respectively. The DOS analysis reveals that near the Fermi level, the conduction band minima (CBM) in these vanadates are predominantly occupied by the outermost orbitals of V atoms, while the valence band maxima (VBM) are primarily contributed by O-2p orbitals. The O-2p orbitals also exhibit strong localization near the Fermi level. Combined with highest occupied molecular orbital-lowest unoccupied molecular orbital (HOMO-LUMO) analysis and population analysis, the bonding interactions in all four crystals mainly arise from the hybridization between V-3p and O-2p orbitals, indicating strong covalent bonding in V—O bonds. Through the analysis of structure-property relationships, the large birefringence is primarily attributed to the pronounced structural anisotropy, high anisotropy index of responsive electron distribution, unique arrangement of anionic groups, and d-p orbital hybridization between V-3d and O-2p orbitals. The calculated birefringence values at a wavelength of 1064 nm for LiV3O8, NaV3O8, KV3O8, and RbV3O8 are 0.28, 0.30, 0.28, and 0.27, respectively.
-
Key words:
- alkali metal vanadate /
- first principles /
- electronic structure /
- birefringence .
-
-
表 1 结构优化前后的AV3O8 (A = Li, Na, K, Rb)晶格参数
Table 1. Lattice parameters of AV3O8 (A = Li, Na, K, Rb) before and after geometry optimization.
Compounds a/nm b/nm c/nm β/(°) b/c Error V/nm3 LiV3O8 Before 0.668 0.360 1.203 107.830 0.299 1.67% 0.275 After 0.695 0.358 1.219 108.397 0.294 0.288 NaV3O8 Before 0.512 0.855 0.744 101.993 1.148 3.92% 0.319 After 0.531 0.877 0.795 113.173 1.103 0.350 KV3O8 Before 0.500 0.839 0.767 98.135 1.094 4.84% 0.319 After 0.516 0.865 0.831 103.812 1.041 0.361 RbV3O8 Before 0.501 0.842 0.791 96.943 1.064 2.16% 0.331 After 0.516 0.865 0.831 100.581 1.041 0.365 表 2 AV3O8 (A = Li, Na, K, Rb)的布居数 (Mulliken)
Table 2. Population of AV3O8 (A = Li, Na, K, Rb) (Mulliken).
Substance Species Atomic population Total Charge/e Bond Bond
populationLength/Å s p d LiV3O8 Li –0.20 0.00 0.00 –0.20 1.20 Li—O –0.06 2.16 V 0.20 0.31 3.20 3.71 1.29 O—V 0.01 2.86 O 1.90 4.68 0.00 6.59 –0.59 O—V 0.92 1.66 NaV3O8 Na 2.07 5.95 0.00 8.01 0.99 Na—O 0.02 2.38 V 0.15 0.33 3.18 3.66 1.34 O—V 0.92 1.64 O 1.89 4.72 0.00 6.62 –0.62 O—V 0.33 1.99 KV3O8 K 2.03 5.91 0.00 7.95 1.05 K—O 0.04 2.72 V 0.15 0.34 3.32 3.72 1.28 O—V 0.97 1.64 O 1.89 4.74 0.00 6.63 –0.63 O—V 0.33 1.98 RbV3O8 Rb 2.04 5.86 0.00 7.90 1.10 Rb—O 0.04 2.97 V 0.16 0.35 3.23 3.74 1.26 O—V 0.98 1.64 O 1.89 4.74 0.00 6.64 –0.64 O—V 0.34 1.98 表 3 AV3O8 (A = Li, Na, K, Rb)的响应电子分布各向异性指数(REDA)
Table 3. REDA of AV3O8 (A = Li, Na, K, Rb) of the electron distribution.
Compounds Groups δ Δn LiV3O8 VO4 0.011 0.28 NaV3O8 V3O8 0.019 0.30 KV3O8 V3O8 0.013 0.28 RbV3O8 V3O8 0.012 0.27 -
[1] Pedrotti F L, Pedrotti L M, Pedrotti L S 2018 Introduction to optics (England: Cambridge University Press) pp333–360 [2] Li X Z, Wang C, Chen X L, Li H, Jia L S, Wu L, Du Y X, Xu Y P 2004 Inorg. Chem. 43 8555 doi: 10.1021/ic049710m [3] Nomura H, Furutono Y 2008 Microelectron. Eng. 85 1671 doi: 10.1016/j.mee.2007.12.017 [4] Aoki K, Miyazaki H T, Hirayama H 2003 Nat. Mater. 2 117 doi: 10.1038/nmat802 [5] Lancry M, Desmarchelier R, Cook K, Poumellec B, Canning J 2014 Micromachines-Basel 5 825 doi: 10.3390/mi5040825 [6] Li R 2013 Z. Krist-Cryst. Mater. 228 526 doi: 10.1524/zkri.2013.1628 [7] Levy M, Jalali A A, Huang X 2009 J. Mater. Sci. Mater. El. 20 43 doi: 10.1007/s10854-007-9433-z [8] Zhang H, Zhang M, Pan S L, Yang Z H, Wang Z, Bian Q, Hou X L, Yu H W, Zhang F F, Wu K, Yang F, Peng Q J, Xu Z Y, Chang K B, Poeppelmeier K R 2015 Cryst. Growth Des. 15 523 doi: 10.1021/cg5016912 [9] Ghosh G 1999 Opt. Commun. 163 95 doi: 10.1016/S0030-4018(99)00091-7 [10] Luo H, Tkaczyk T, Sampson R, Dereniak E L 2006 Proc. SPIE 6119 136 doi: 10.1117/12.651149 [11] Guoqing Z, Jun X, Xingda C, Heyu Z, Siting W, Ke X, Fuxi G 1998 J. Cryst. Growth 191 517 doi: 10.1016/S0022-0248(98)00162-6 [12] Appel R, Dyer C D, Lockwood J N 2002 Appl. Opt. 41 2470 doi: 10.1364/AO.41.002470 [13] Cyranoski D 2009 Nature 457 953 doi: 10.1038/457953a [14] Krainer L, Paschotta R, Lecomte S, Moser M, Weingarten K J, Keller U 2002 IEEE J. Quantum Electron. 38 1331 doi: 10.1109/JQE.2002.802967 [15] Lisinetskii V A, Grabtchiko A S, Demidovich A A, Burakevich V N, Orlovich V A, Titov A N 2007 Appl. Phys. B 88 499 doi: 10.1007/s00340-007-2757-9 [16] Vodchits A I, Orlovich V A, Apanasevich P A 2012 J. Appl. Spectrosc. 78 918 doi: 10.1007/s10812-012-9554-x [17] Yu H, Liu J, Zhang H, Kaminskii A A, Wang Z, Wang J 2014 Laser Photonics Rev. 8 847 doi: 10.1002/lpor.201400022 [18] Lei B H, Kong Q, Yang Z H, Yang Y, Wang Y J, Pan S L 2016 J. Mater. Chem. C 4 6295 doi: 10.1039/C6TC01512C [19] Li K X, Zhang X Y, Chai B Q, Yu H W, Hu Z G, Wang J Y, Wu Y C, Wu H P 2025 Chem. Eur. J. 31 e202403515 doi: 10.1002/chem.202403515 [20] Li K X, Wu H P, Yu H W, Hu Z G, Wang J Y, Wu Y C 2024 Chem. Commun. 60 12734 doi: 10.1039/D4CC04530K [21] Lei B H, Yang Z H, Pan S L 2017 Chem. Commun. 53 2818 doi: 10.1039/C6CC09986F [22] Huang Y, Zhang X Y, Zhao S G, Mao J G, Yang B P 2024 J. Mater. Chem. C 12 7286 doi: 10.1039/D4TC01200C [23] Zhang S Z, Dong L F, Xu B H, Chen H G, Huo H, Liang F, Wu R, Gong P F, Lin Z S 2024 Inorg. Chem. Front. 11 5528 doi: 10.1039/D4QI01442A [24] Cheng J L, Xu D, Lu J, Zhang F F, Hou X L 2023 Inorg. Chem. 62 20340 doi: 10.1021/acs.inorgchem.3c03325 [25] Chen Z X, Xu F, Cao S N, Li Z F, Yang H X, Ai X P, Cao Y L 2017 Small 13 1603148 doi: 10.1002/smll.201603148 [26] Cao X Y, Yang Q, Zhu L M, Xie L L 2018 Ionics 24 943 doi: 10.1007/s11581-018-2477-7 [27] Yang H, Li J, Zhang X G, Jin Y L 2008 J. Mater. Process. Technol. 207 265 doi: 10.1016/j.jmatprotec.2007.12.122 [28] Zhu L M, Li W X, Xie L L, Yang Q, Cao X Y 2019 Chem. Eng. J. 372 1056 doi: 10.1016/j.cej.2019.05.009 [29] Feng L L, Zhang W, Xu L N, Li D Z, Zhang Y Y 2020 Solid State Sci. 103 106187 doi: 10.1016/j.solidstatesciences.2020.106187 [30] Kim H J, Jo J H, Choi J U, Voronina N, Myung S T 2020 J. Power Sources 478 229072 doi: 10.1016/j.jpowsour.2020.229072 [31] Shchelkanova M, Shekhtman G, Pershina S, Vovkotrub E 2021 Materials 14 6976 doi: 10.3390/ma14226976 [32] Wu W Z, Ding J, Peng H R, Li G C 2011 Mater. Lett. 65 2155 doi: 10.1016/j.matlet.2011.04.048 [33] Zhu J Z, Li X L, Chen S Y, Huang C M, Feng J J, Kuang Q, Fan Q H, Dong Y Z, Zhao Y M 2020 Electrochim. Acta 355 136799 doi: 10.1016/j.electacta.2020.136799 [34] Wadsley A D 1957 Acta Crystallogr. 10 261 doi: 10.1107/S0365110X57000821 [35] Bachmann H G, Barnes W H 1962 Can Mineral 7 219 [36] Baddour-Hadjean R, Boudaoud A, Bach S, Emery N, Pereira-Ramos J P 2014 Inorg. Chem. 53 1764 doi: 10.1021/ic402897d [37] Oka Y, Yao T, Yamamoto N 1997 Mater. Res. Bull. 32 1201 doi: 10.1016/S0025-5408(97)00096-2 [38] Segall M D, Lindan P J D, Probert M J 2002 J. Phys. : Condens. Matter 14 2717 doi: 10.1088/0953-8984/14/11/301 [39] Payne M C, Teter M P, Allan D C, Arias T A, Joannopoulos A J 1992 Rev. Mod. Phys. 64 1045 doi: 10.1103/RevModPhys.64.1045 [40] Srivastava G P, Weaire D 1987 Adv. Phys. 36 463 doi: 10.1080/00018738700101042 [41] Monkhorst H J, Pack J D 1976 Phys. Rev. B 13 5188 doi: 10.1103/PhysRevB.13.5188 [42] Perdew J P, Burke K, Ernzerhof M 1996 Phys. Rev. Lett. 77 3865 doi: 10.1103/PhysRevLett.77.3865 [43] Hamann D R, Schlüter M, Chiang C 1979 Phys. Rev. Lett. 43 1494 doi: 10.1103/PhysRevLett.43.1494 [44] Kong Q R, Yang Y, Liu L L, Bian Q, Lei B H, Li L P, Yang Z H, Su Z, Pan S L 2016 J. Mater. Res. 31 488 doi: 10.1557/jmr.2016.29 [45] 阎守胜 2011 固体物理基础(北京: 北京大学出版社) Yan S S 2011 Fundamentals of Solid-state Physics (Beijing: Peking University Press [46] Hiscocks J, Frisch M J 2009 Gaussian 09: IOps Reference 9 (USA: Gaussian [47] Riffet V, Contreras-Garcıa J, Carrasco J, Calatayud M 2016 J. Phys. Chem. C 120 4259 doi: 10.1021/acs.jpcc.5b11600 [48] Mulliken R S 1931 Chem. Rev. 9 347 doi: 10.1021/cr60034a001 [49] 张博, 王云杰, 齐亚杰, 和志豪, 丁家福, 苏欣 2024 人工晶体学报 53 999 doi: 10.16553/j.cnki.issn1000-985x.2024.06.007 Zhang B, Wang Y J, Qi Y J, He Z H, Ding J F, Su X 2024 J. Synth. Cryst. 53 999 doi: 10.16553/j.cnki.issn1000-985x.2024.06.007 [50] Su X, Wang Y J, Yang Z H, Huang X C, Pan S L, Li F, Lee M H 2013 J. Phys. Chem. C 117 14149 doi: 10.1021/jp4013448 [51] Tudi A, Han S, Yang Z H, Pan S L 2022 Coord. Chem. Rev. 459 214380 doi: 10.1016/j.ccr.2021.214380 [52] Wang X Y, Zhang B B, Yang D Q, Wang Y 2022 Dalton Trans. 51 14059 doi: 10.1039/D2DT02573F [53] Bai S, Yang D Q, Zhang B B, Li L, Wang Y 2022 Dalton Trans. 51 3421 doi: 10.1039/D2DT00183G [54] Chu Y, Wang H S, Chen Q, Su X, Chen Z X, Yang Z H, Li J J, Pan S L 2024 Adv. Funct. Mater. 34 2314933 doi: 10.1002/adfm.202314933 [55] Ding Y Y, Zhu M M, Wang J B, Li B, Qi H X, Liu L L, Chu Y Q 2024 Inorg. Chem. 63 20003 doi: 10.1021/acs.inorgchem.4c03593 [56] Lin L, Jiang X X, Wu C, Lin Z S, Huang Z P, Humphrey M G, Zhang C 2021 Dalton Trans. 50 7238 doi: 10.1039/D1DT00536G [57] Su X, Chu Y, Yang Z H, Lei B H, Cao C, Wang Y, Pan S L 2020 J. Phys. Chem. C 124 24949 doi: 10.1021/acs.jpcc.0c06991 [58] Bai S, Zhang X, Zhang B B, Li L, Wang Y J 2021 Inorg. Chem. 60 10006 doi: 10.1021/acs.inorgchem.1c01298 [59] Li S, Dou D, Chen C, Shi Q, Zhang B B, Wang Y J 2024 Inorg. Chem. 63 24076 doi: 10.1021/acs.inorgchem.4c04843 [60] Chu Y, Wang H S, Abutukadi T, Li Z, Mutailipu M, Su X, Yang Z H, Li J J, Pan S L 2023 Small 19 2305074 doi: 10.1002/smll.202305074 [61] Liu H J, Liang C W, Liang W I, Chen H J, Yang J C, Peng C Y, Chu Y H 2012 Phys. Rev. B Condens. Matter Mater. Phys. 85 014104 doi: 10.1103/PhysRevB.85.014104 [62] 王云杰, 和志豪, 丁家福, 苏欣 2025 人工晶体学报 54 85 doi: 10.16553/j.cnki.issn1000-985x.2024.0190 Wang Y J, He Z H, Ding J F, Su X 2025 J. Synth. Cryst. 54 85 doi: 10.16553/j.cnki.issn1000-985x.2024.0190 [63] 丁家福, 和志豪, 王云杰, 苏欣 2025 人工晶体学报 54 95 doi: 10.16553/j.cnki.issn1000-985x.2024.0193 Ding J F, He Z H, Wang Y J, Su X 2025 J. Synth. Cryst. 54 95 doi: 10.16553/j.cnki.issn1000-985x.2024.0193 [64] 储冬冬, 杨志华, 潘世烈 2024 人工晶体学报 53 1475 doi: 10.16553/j.cnki.issn1000-985x.20240731.005 Chu D D, Yang Z H, Pan S L 2024 J. Synth. Cryst. 53 1475 doi: 10.16553/j.cnki.issn1000-985x.20240731.005 -