DHMS中掺杂Cu2+电子顺磁共振参量及局部结构

李超英，付金仙，董 明，黄 瑛，毛杰健

（1. 上饶师范学院物理与电子信息学院 江西 上饶 334000；2. 上饶师范学院图书馆 江西 上饶 334000）

·物理电子学·

DHMS中掺杂Cu2+电子顺磁共振参量及局部结构

李超英1，付金仙1，董 明2，黄 瑛1，毛杰健1

（1. 上饶师范学院物理与电子信息学院 江西 上饶 334000；2. 上饶师范学院图书馆 江西 上饶 334000）

基于晶体场模型，利用3d9离子斜方（C4v）伸长八面体对称EPR参量的高阶微扰公式，计算了（NH4）2Mg（SO4）2·6H2O：Cu2+的g因子gii和超精细结构常数Aii。公式中晶场参量由重叠模型确定，计算中考虑了d轨道基态波函数的混合。研究结果表明，晶体（NH4）2Mg（SO4）2·6H2O中络离子［Cu（H2O）6］2+的Cu2+-H2O键长Rx约为0.187 2 nm，Ry约为0.203 3 nm，Rz约为0.229 2 nm；中心金属离子基态波函数混合系数α和β分别为0.995和0.099 9。所得EPR参量理论计算与实验符合很好，并对上述结果进行了讨论。

晶体场理论； Cu2+离子； 缺陷结构； 晶体DHMS； 电子顺磁共振

晶体（NH4）2Mg（SO4）2·6H2O （DHMS）的空间结构为P21/a，每单个晶胞中有2个二价Mg2+离子（即Z= 2），晶格常数a = 0.931 6 nm，b = 1.258 0 nm，c = 0.620 2 nm，β = 107.09°［1］。晶体DHMS中金属Mg2+离子周围有6个最邻近H2O分子，形成的［Mg（H2O）6］2+基团属于斜方压缩八面体（C4v）点群对称结构［1-3］。研究成果均表明，当掺杂离子进入基质晶体后，掺杂离子与其所替换的母体离子所处的局域结构往往不同。掺杂离子的电子顺磁共振（electron paramagnetic resonance， EPR）谱对周围环境极为敏感。因此，对顺磁中心的EPR参量的研究可以用来确定掺杂过渡金属离子在晶体中所处的局部结构［4］。目前，对晶体DHMS中掺杂过渡金属离子的EPR实验及其局域结构的理论研究均有报道［5-6］。文献［7］报道了DHMS：Cu2+的EPR参量的实验数据（g因子gii和精细结构常数Aii，i = x， y， z），并采用简单二阶微扰公式计算DHMS： Cu2+的EPR参量，计算中未能将杂质中心的局域结构与晶场参量相互联系，因此未能获得杂质中心结构信息。材料的光学、磁学性能与掺杂离子在晶体中所处的局域结构密切相关，所以研究材料中掺杂离子周围的局域结构，对理解过渡金属离子与基质晶体的相互作用以及认识材料的微观机理都非常有意义。因此，对DHMS： Cu2+的EPR参量作进一步的合理解释并获得杂质中心的局域结构是必要的。本文采用斜方伸长八面体（C4v）对称下Cu2+离子的EPR参量高阶（三阶）微扰公式计算DHMS： Cu2+的g因子和精细结构常数A因子，公式中涉及的晶场参量由重叠模型并联系晶体中掺杂离子局域结构参量确定，并考虑了基态波函数中能态2A1g（θ）和2A1g（ε）之间的混合对EPR参量的贡献。

1 理论模型与计算公式

当Cu2+掺入晶体DHMS后，它将占据母体Mg2+离子位置，形成［Cu（H2O）6］2+基团并保留母体的斜方（C4v）对称。由文献［7］的EPR实验结果给出的gzz＞ gxx，gyy＞2可知，杂质Cu2+离子与周围的6个H2O分子构成［Cu（H2O）6］2+基团属于斜方伸长八面体（C4v）对称，不同于晶体DHMS中母体［Mg（H2O）6］2+基团的斜方压缩八面体对称。Cu2+离子在斜方伸长八面体晶场，能级较低的轨道双重态2Eg将分裂为两个轨道单重态2A1g（ ）和2A1g（ ），能级较高的轨道三重态2T2g将分裂成3个轨道单重态2B1g（ζ）、2B2g（ζ）和2B3g（η）［8-9］。斜方对称下能态2A1g（θ）和2A1g（ε）具有相同的不可约表示，此时，基态应为2A1g（θ）和2A1g（ε）之间的混合，基态波函数为［10-11］：

态混合系数α和β满足归一化关系［9］：

由微扰理论可得，斜方对称下3d9（Cu2+）离子的EPR参量（g因子gxx、gyy、gzz和精细结构常数Axx、Ayy、Azz）的高阶微扰公式分别为［11］：

式中，gs≈ 2.002 3为自由电子对应的g因子值；ζ = Nζd0为晶体中Cu2+离子的旋-轨耦合系数；ζd0≈ 829 cm-1［12］为自由Cu2+离子的旋-轨耦合系数；N为表征体系共价性的共价因子；κc为掺杂Cu2+离子的各项同性芯区极化常数，κc′为来自斜方畸变使中心离子3d-3s（4s）轨道混合引起的各向异性芯区极化常数；P0（P0≈ 388 ×10-4cm-1［13］）为自由Cu2+离子的偶极超精细结构常数。

式（3）中能级差Ei（i = 1～4）可由晶场参量Dq、Ds、Dt、Dξ和Dη表示为：

由重叠模型［14］联系晶体DHMS： Cu2+的局域结构，晶场参量Dq、Ds、Dt、Dξ、Dη可分别表示为：

式中，指数律系数t2≈ 3，t4≈ 5［3-4，8］；A2（R0）与4（R0）为对应参考距离R0的内禀参量，R0== （Rx+ Ry+ Rz）/3；Rx和Ry表示垂直于C4轴平面上的Cu-H2O键长；Rz表示平行于C4轴的Cu-H2O键长。研究表明，对3dn离子在许多晶体中处在八面体晶场的情况下，比率2（R0）/A4（R0）约为9～12之间［3，8-9，16］，则取A2（R0）/ A4（R0） = 12。这样，杂质中心［Cu（H2O）6］2+基团的局域结构参数Rx、Ry和Rz就与晶场参量相互关联，即与式（4）中能级Ei及与式（3）中的EPR参量相联系。

根据对类似晶体晶体（NH4）2Mg（SO4）2·6H2O中掺杂Cu2+离子形成的［Cu（H2O）6］2+基团的光谱研究，光谱参量A4（R0） ≈ 620 cm-1［3］。对Cu2+离子掺杂情况，其各向同性芯区极化常数κc处于0.2～0.5之间［7-9，17-18］，取κc≈ 0.26。由于Cu2+离子中心3d-3s（4s）轨道的混合对超精细结构常数各项异性的贡献，各向异性芯区极化常数κc′取为0.04。

拟合EPR参量理论和实验相符，DHMS晶体中Cu2+离子中心局域结构参数Rx、Ry及Rz分别为Rx≈0.187 2 nm，Ry≈ 0.203 3 nm，Rz≈ 0.229 2 nm。基态轨道波函数为：

相应的EPR参量计算结果（Cal.c）如表1所示。为了便于对比，忽略基态2A1g（θ）和2A1g（ε）之间的混合（即α = 1，β = 0）的计算结果（Cal.b）及利用简单二阶微扰公式所得结果（Cal.a）也列于表1中。Cal.a为文献［7］基于点电荷模型的简单微扰公式计算结果；Cal.b为基于式（3），但忽略了基态波函数中能态2A1g（θ）和2A1g（ε）的混合对EPR参量的贡献（即α = 1，β = 0）所得的计算结果；Cal.c为采用式（3），并考虑基态波函数中能态2A1g（θ）和2A1g（ε）的混合对EPR参量的贡献所得结果。

表1 （NH4）2Mg（SO4）2·6H2O： Cu2+的g因子和超精细结构常数理论和实验值

2 分析与讨论

由表1可知，计算所得EPR参量g因子和精细结构常数A因子理论值与实验值符合很好，其结果相比文献［7］也有很大的改进。因此，本文采用的斜方对称下Cu2+离子EPR参量计算式（1），及通过拟合EPR参量所得到的相关结构参数是合理的。

1） 在八面体晶场作用下，随着杂质中心对称度降低，基态之间的混合和配体离子的旋-轨耦合作用都能影响最后的计算结果。从表1可以发现，采用式（3）并考虑基态波函数中能态2A1g（θ）和2A1g（ε）的混合对EPR参量的贡献所得结果（Cal.c），相比忽略基态波函数中能态2A1g（θ）和2A1g（ε）的混合计算（Cal.b）及前人采用简单的二阶微扰公式所得结果（Cal.a），与实验值符合更好，而且，通过拟合EPR参量所得能态混合系数α与文献［10］通过分析晶体C6H7KO7：Cu2+的EPR参量所得混合系数α ≈ 0.995一致，进一步说明本文所得基态波函数是合理的。此外，由于中心Cu2+离子的旋轨耦合系数（ζd0≈ 829 cm-1 ［12］）比配体O2-离子）大很多，因此，以上计算忽略了来自配体旋-轨耦合作用的贡献。由于［Cu（H2O）6］2+基团的共价性很弱，这种近似算法是可靠的，并被很多研究者采用［4，9，11，18］。

2） 过渡金属离子掺杂局域结构往往很难被实验测定，本文通过分析晶体DHMS： Cu2+的EPR参量的方法，获得中杂质Cu2+离子周围配体的Cu2+-H2O键长分别为Rx≈ 0.187 2 nm， Ry≈ 0.203 3 nm， Rz≈0.229 2 nm，这与纯晶体（NH4）2Cu（SO4）2·6H2O中通过X衍射实验获得的Cu2+-H2O键长（即Rx≈0.196 6nm，Ry≈ 0.207 2nm，Rz≈ 0.223 0 nm［3］）非常接近，说明了本文所得晶体DHMS中［Cu（H2O）6］2+基团的局域结构是可靠的。

3） 各向同性芯区极化常数κc，由对Cu2+离子、中心金属离子未配对电子自旋密度χ ≈ -3.4 a.u.［13］，3d轨道径向波函数负3次方的期望值

［1］ MASLEN E N， RIDOUT S C， WATSON K J， et al. The structures of tutton's salts. I. diammonium hexaaquamagnesium-（II） sulfate［J］. Acta Cryst C， 1988，44（3）： 409-412.

［2］ PETRASHEN V E， YABLOKOV Y V， DAVIDOVICH R L. The lattice structure parameters and configuration of Cu2+Jahn-Teller centres in tutton salt crystals［J］. Phys Status Solidi （B）， 1980， 101（1）： 117-125.

［3］ ZHENG W C， ZHANG D T， SU P， et al. Studies of the optical band positions and EPR g factors for Cu（H2O）62+centers in tutton salt crystals［J］. Spectrochim Acta A， 2011，81（1）： 548-551.

［4］ 张华明， 易江林， 肖文波. 晶体 zinc 1-malate trihydrate：Cu2+的自旋哈密顿参量及局部结构的理论研究［J］. 人工晶体学报， 2014， 43（3）： 708-711. ZHANG Hua-ming， YI Jiang-lin， XIAO Wen-bo. Theoretical investigation of the local structures and spin Hamiltonian parameters for zinc 1-malate trihydrate： Cu2+［J］. J Synth Cryst， 2014， 43（3）： 708-711.

［5］ WEI Q. Investigations of the optical and EPR spectra for Cr3+ions in diammonium hexaaqua magnesium sulphate single crystal ［J］. Acta Phys Pol A， 2010， 118（4）： 670-672.

［6］ KRIPAL R， PANDEY S D， MISRA M G. EPR， optical absorption and superposition model studies of Fe3+-doped diammonium hexaaqua magnesium sulfate： A case of hyperfine structure［J］. Appl Magn Reson， 2013， 44（11）：1295-1310.

［7］ KRIPAL R， MISRA M G. EPR and optical absorption study of Cu2+-doped diammonium hexaaqua magnesium sulphate single crystals［J］. Appl Magn Reson， 2013， 44（3）： 411-434.

［8］ ZHANG H M， WU S Y， KUANG M Q， et al. Investigation of the EPR parameters and local structures for Cu2+in Bis（L-asparaginato） M（II） catalysts （M = Zn， Cd， Mg）［J］. J Phys Chem Solids， 2012， 73（7）： 846-850.

［9］ DONG H N， WU S Y， LI P. Theoretical explanation of EPR parameters for Cu2+ion in TiO2crystal［J］. Phys Status Solidi （B）， 2004， 241（8）： 1935-1938.

［10］ AYGUN Z. Variable temperature EPR studies of Cu2+and VO2+doped potassium dihydrogen citrate （C6H7KO7）［J］. Spectrochim Acta A， 2013， 104： 130-133.

［11］ HUANG Y P， WANG L J， FENG W L. Theoretical investigation of the local structure of the KHP2O4：Cu2+single crystal［J］. Radiat Eff Defect Solids， 2009， 164（3）：183-186.

［12］ WU S Y， ZHANG H M， XU P. et al. Studies on the defect structure for Cu2+in CdSe nanocrystals［J］. Spectrochim Acta A， 2010， 75（1）： 230-234.

［13］ MCGARVEY B R. The isotropic hyperfine interaction［J］. Phys Chem， 1967， 71（1）： 51-66.

［14］ NEWMAN D J， NG B. The superposition model of crystal fields［J］. Rep Prog Phys， 1989， 52（6）： 699-763.

［15］ KUANG M Q， WU， S Y， HU X F， et al. Theoretical studies of the defect structures and the spin Hamiltonian parameters for the Pd+and Pd3+centers in NH4Cl［J］. Comput Theor Chem， 2013， 1022： 45-49.

［16］ 邬劭轶， 林季资， 付强， 等. 锐钛矿中填隙V4+的EPR参量的理论研究［J］. 电子科技大学学报， 2008， 37（4）： 556-557. WU Shao-yi， LIN Ji-zi， FU Qiang， et al. Studies of the EPR parameters of the interstitial V4+in anatase［J］. Journal of University of Electronic Science and Technology of Chian，2008， 37（4）： 556-557.

［17］ KUANG M Q， WU S Y， ZHANG H M， et al. Theoretical studies on the local structure and spin Hamiltonian parameters for the orthorhombic Cu2+center in LiNbO3［J］. Optik， 2012， 123（18）， 1601-1604.

［18］ 李超英， 杨建荣. Ca（OD）2：Cu2+局部结构及自旋哈密顿参量研究［J］. 电子科技大学学报， 2012， 41（3）： 405-407. LI Chao-yin， YANG Jian-rong. Theoretical explanation of the local structures and spin Hamiltonian parameters in Ca（OD）2of the tetragonal Cu2+center［J］. Journal of University of Electronic Science and Technology of Chian，2012， 41（3）： 405-407.

［19］ ZHANG H M， WAN X， ZHANG Z M. Theoretical studies of the spin Hamiltonian parameters and local structures for the tetragonal Cu2+and Ni3+centers in Mg2TiO4［J］. J Alloy Compd， 2013， 549： 226-230.

［20］ WU S Y， ZHANG H M， XU P. Studies on the defect structure for Cu2+in CdSe nanocrystals［J］. Spectrochim Acta A， 2010， 75（1）： 230-234.

编 辑 黄 莘

TheoretiPcaarl aImnveetsetrisg aatniodn L oofc tahl eS E trluecctturorn es P foar r aDmHaMgnSe：t iCc uR2e+sonance

LI Chao-ying1， FU Jin-xian1， DONG Ming2， HUANG Ying1， and MAO Jie-jian1
（1.School of Physics and Electronic Information， Shangrao Normal University Shangrao Jiangxi 334000；2. Library of Shangrao Normal University Shangrao Jiangxi 334000）

Based on the crystal field theory， the anisotropic g factors giiand the hyperfine structure constants Aiifor （NH4）2Mg（SO4）2·6H2O：Cu2+are theoretically investigated by using the high-order perturbation formulas of these electron paramagnetic resonance （EPR） parameters for a 3d9ion in rhombically （C4v） elongated octahedra. In the calculated formulas， the crystal field parameters are set up from the superposition model and the admixture of d-orbitals in the ground state wave function is taken into account. According to the calculations， the local structure values ［Cu（H2O）6］2+clustering in （NH4）2Mg（SO4）2·6H2O crystal are： Rx≈ 0.1872 nm， Ry≈ 0.2033 nm， Rz≈ 0.2292 nm； and the mixing coefficients of the ground state wave function parameters α and β are 0.995 and 0.0999，respectively. The calculated results show good agreement with the experimental data， and the reasonableness of these results is discussed.

crystal field theory； Cu2+ion； defect structure； DHMS crystal； electron paramagnetic resonance

O737

A

10.3969/j.issn.1001-0548.2016.02.008

2014 - 07 - 24；

2014 - 12 - 19

国家自然科学基金（11365017， 11465015）；江西省教育厅项目（GJJ 151058）

李超英（1958 - ），女，教授，主要从事光电信息材料方面的研究.