PHZS:Cu2+的EPR谱和局域结构研究

张华明

(1.南昌航空大学光电检测工程技术实验室 南昌 330063；2.南昌航空大学教育部无损检测重点实验室 南昌 330063)

晶体K2Zn(SO4)2·6H2O (PHZS)属于单斜晶系，空间群为(P121/a1)，每个晶胞含2个Zn原子，晶格常数[1]a≈ 0.903 4 nm，b≈ 1.218 4 nm，c≈0.614 8 nm，δ=104.8°，该晶体由一系列不规则的KO8多面体和Zn(H2O)6八面体通过氢键与SO4相连接而构成网格状结构。其中Zn2+与周围6个H2O分子构成[Zn(H2O)6]2+基团，该位置属正交(D2)点群对称。众所周知，晶体的光学、磁学等性能与掺杂离子的局域结构密切相关，而掺杂离子所处的局域环境往往不同于母体位置[2-3]。因此，研究晶体中掺杂离子的局域结构对理解掺杂离子影响材料性能的微观机理非常重要。电子顺磁共振(EPR)谱强烈依赖于顺磁离子所处局域环境，并可通过分析作为其实验结果的EPR参量定量地确定掺杂离子周围的局部结构[4]。文献[1]报道了PHZS:Cu2+的EPR参量(g因子gx,gy,gz和精细结构常数Ax,Ay,Az)的实验数据，但是上述EPR参量实验结果以及杂质离子局域结构信息至今尚未得到定量的理论分析。本文基于3d9离子EPR参量的高阶微扰公式计算PHZS: Cu2+的g因子和精细结构常数A因子，相关的正交晶场参量由重叠模型确定，并考虑了配体轨道和旋轨耦合作用以及基态波函数中两个态2A1g(θ)和2A1g(ε)之间混合对EPR参量的影响。

1 理论模型及计算公式

Cu2+掺入PHZS晶体后会替代母体Zn2+位置并与周围6个最近邻水分子形成[Cu(H2O)6]2+基团。根据文献[1]测得的g因子各向异性gz＞gx，gy＞ 2，可以判断杂质[Cu(H2O)6]2+基团属于正交伸长八面体。在此环境下，Cu2+(3d9)离子较低的轨道双重态(2Eg)将分裂成2A1g(θ)和2A1g′(ε)两个轨道单重态，较高的轨道三重态(2T2g)将分裂成2B1g(ζ)、2B2g(ξ)和2B3g(η)三个轨道单重态。由于正交对称下2A1g(θ)和2A1g′(ε)具有相同的不可约表示，二者将发生混合，因此基态应为2A1g(θ)和2A1g′(ε)的混合态，即：

式中，N为共价因子，表征体系共价性；α和β为态混合系数，并满足归一化关系：

由双旋−轨耦合模型和微扰理论，可得3d9离子在正交对称下EPR参量的高阶微扰公式为[5]：

式中，ζd0和ζp0为自由Cu2+和配体O2−的旋−轨耦合系数；Nγ(γ=t,e)和λ(或λs)表示归一化因子和轨道混合系数；A为积分R＜ns|∂/∂y|npy＞，其中R为金属−配体间距；κc和κc′分别为Cu2+的各向同性和各向异性芯区极化常数；P0(≈ 388 ×10−4cm−1)为自由Cu2+离子的偶极超精细结构参量[6-7]。由离子簇模型[8], 分子轨道系数Nt，Ne和λt，λe满足以下近似关系：

和归一化条件：

式(3)中能级差Ei(i=1,2,3,4)可由立方场参量Dq和正交场参量表示为：

以上立方和正交场参量可由重叠模型[9]得到：

以上公式中，指数律[10-13]系数为对应参考距离的内禀参量。为垂直于C2轴的杂质−配体键长，Rz为平行于C2轴的杂质−配体键。利用PHZS晶体中的金属−配体平均键长，基于Slater型自洽场波函数可以计算出群重叠积分归一化因子Nγ和轨道混合系数λγ可由式(4)～式(6)获得。利用自由Cu2+和配体O2−离子的旋轨耦合[7]系数151 cm−1，可得式(4)中的旋轨耦合系数ς和ς′和轨道缩小因子k和k′，所得结果如表1所示。对于八面体晶场下的3dn离子，文献[10,14]中广泛采用的关系式也用于此处。这样，掺杂离子的结构参数Rx,Ry和Rz(特别是正交畸变)便与晶场参量(特别是四角场参量)以及杂质局部结构相关联。因此，通过分析PHZS:Cu2+的EPR，可获得PHZS晶体中Cu2+杂质中心的局域结构信息。

表1 K2Zn(SO4)2·6H2O:Cu2+晶体的相关参数

根据已有[Cu(H2O)6]2+基团的光谱，式(8)中光谱参量Ā4(R)可以取为650 cm−1。对晶体中的Cu2+离子，各向同性芯区极化常数κc一般在0.2～0.3范围，取κc≈ 0.252。此外，考虑到杂质中心正交畸变引起的各向异性3d-3s(4s)轨道混合，各向异性芯区极化常数κc′取为0.033。

拟合EPR参量理论与实验值相符，得到PHZS:Cu2+的结构参数Rx，Ry，Rz和基态轨道波函数分别为：

对应的EPR参量计算结果(Cal.c)如表2所示。为了突出基态波函数混合和配体轨道以及旋轨耦合作用对EPR参量的贡献，忽略基态波函数混合(即取α= 1,β= 0)的计算结果(Cal.a)以及忽略配体轨道和旋轨耦合作用所得结果(Cal.b)也一并列于表2。表中，A因子单位为10−4/cm；Cal.a考虑配体轨道和旋轨耦合作用，但忽略基态波函数混合对EPR参量的贡献(即取α= 1,β= 0)的计算结果；Cal.b考虑基态波函数混合，但忽略配体轨道和旋轨耦合作用(即取ζ=ζ′ =Nζd0，k=k′ =N)的计算结果；Cal.c同时考虑配体轨道和旋轨耦合作用以及基态波函数混合的结果。

表2 K2Zn(SO4)2·6H2O:Cu2+的EPR参量理论计算及实验值

2 讨 论

由表1可以看出，基于高阶微扰式(3)并采用式(9)中的Cu2+离子结构参数所得的EPR参量理论值与实验符合很好。可见，本文所采用计算公式和通过拟合EPR参量所得基态波函数以及杂质中心的局域结构信息是合理的。

1)八面体晶场下的3d9离子随局部对称度降低，会导致立方基态两个分量2A1g(θ)和2A1g(ε)之间的混合，同时配体轨道和旋轨耦合作用也会对EPR参量产生贡献。从表1可以看出，考虑基态波函数之间混合以及配体轨旋轨耦合作用贡献的结果(Cal.c)与实验值符合很好，然而忽略基态波函数中能态2A1g(θ)和2A1g(ε)混合的计算(Cal.b)以及忽略配体轨道和旋轨耦合作用的结果(Cal.a)与实验值符合较差。通过拟合EPR参量所得基态波函数混合系数α≈ 0.978与前人通过分析EPR参量获得的Cu2+杂质在类似的Zn(C3H3O4)2(H2O)2晶体[5]中的数值α≈ 0.995接近，说明本文所得基态波函数是合理的。

2)掺杂离子局域结构通常难以由实验直接测出，而本文通过拟合EPR参量获得的PHZS:Cu2+的Cu2+-H2O键长(即Rx≈ 0.197 nm,Ry≈ 0.213 nm,Rz≈ 0.224 nm)，与通过X衍射实验获得的同系物K2Cu(SO4)2·6H2O晶体中Cu2+-H2O键长[15](即Rx≈ 0.194 3 nm，Ry≈0.206 9 nm，Rz≈0.227 8 nm)接近，进一步说明本文通过分析EPR参量获得PHZS: Cu2+的局域结构是合理的。

3)式(3)中各向同性芯区极化常数[16]κc可由式κc≈ −2χ/(3＜r−3＞)获得，其中＜r–3＞≈8.252 a.u为铜离子波函数负3次方期望值[17]；χ≈−3.12 a.u.为铜离子未成对自旋密度[16]。由此可得κc≈ 0.25，这与本文所采用的κc≈ 0.252较吻合，因而是合理的。

3 结束语

基于3d9(Cu2+)离子在正交伸长八面体对称下的高阶微扰公式，计算了PHZS: Cu2+的EPR参量。计算中考虑了基态波函数混合以及配体轨道和旋轨耦合作用的贡献，所得理论值与实验值符合很好。通过分析EPR参量，获得了[Cu(H2O)6]2+基团局域结构以及杂质Cu2+离子的基态波函数。