绿色氟羟磷灰石热处理及颜色成因探究
岳素伟1,*, 剡晓旭2,*, 陈思敏2, 罗洁2
1.广州职业技术大学珠宝学院, 广东 广州 511483
2.广州城市理工学院珠宝学院, 广东 广州 510800
*通讯作者 e-mail: yuesw@gzpyp.edu.cn; yanxiaoxu@gcu.edu.cn

作者简介: 岳素伟, 1985年生,广州职业技术大学珠宝学院教授 e-mail: yuesw@gzpyp.edu.cn

摘要

磷灰石标准化学式为[Ca10(PO4)6 (F, Cl, OH)2], 是一种应用广泛的无机非金属材料和重要的地质指示矿物, 也是常见的宝石品种。 磷灰石中Ca2+可被二价过渡族金属阳离子(Fe2+、 Mn2+等)、 稀土元素离子(La3+、 Ce3+、 Gd3+等)或碱金属和碱土金属离子(K+、 Na+、 Sr2+等)类质同象替代, 从而呈现黄色、 绿色、 蓝色、 紫色及变色效应。 绿色氟羟磷灰石经650 ℃-700℃热处理转变为鲜艳的蓝色, 具有较高的观赏性和商业价值。 运用电子探针(EPMA)、 傅里叶变换中红外光谱(Mid FT-IR)、 紫外-可见光光谱(UV-Vis)和电子顺磁波谱(EPR), 测试与观察热处理前后氟羟磷灰石的化学组成以及不同温度下颜色变化。 结果显示, 样品化学成分主要由Ca2+和P5+及少量的N5+组成, N5+主要以—NH2形式少量替代[PO3], 表明样品形成与生物化学作用密切相关。 样品热处理前后的UV-Vis特征显示, 垂直于c轴方向的吸收带主要位于264~402和584~804 nm区域, 其颜色主要由Ce3+和Mn2+类质同象替代畸变八面体中的Ca2+形成蓝绿区透过窗所致。 热处理前绿色磷灰石样品中的缺陷主要为离子心F--O--F-, 以及由Ce3+类质同象替代CaⅠ位所致的空穴心O2--O--V和OH--O--V, 热处理后由于空穴心的恢复则主要表现为Mn2+替代CaⅠ和CaⅡ位畸变八面体引起的离子心缺陷。 结合EPR波谱表征可知, Mn2+在CaⅠ和CaⅡ位置均有分布, 是样品呈现蓝色直接原因; Ce3+在CaⅠ位置类质同象替代程度较高时, 样品颜色以绿色为主, 反之呈现蓝色。

关键词: 绿色磷灰石; 热处理; 颜色成因; EPR谱
中图分类号:P575.4 文献标志码:A
The Heat Treatment Process and Color Mechanism of Green Fluor-Hydroxyapatite
YUE Su-wei1,*, YAN Xiao-xu2,*, CHEN Si-min2, LUO Jie2
1. Jewelry Institute, Guangzhou Polytechnic University, Guangzhou 511483, China
2. School of Jewelry, Guangzhou City University of Technology, Guangzhou 510800, China
*Corresponding authors
Abstract

Apatite is widely used as an inorganic non-metallic material, a geological indicator mineral, and a common gemstone variety. The chemical formula for apatite is [Ca10(PO4)6(F, Cl, OH)2], Ca2+sites can be isomorphously substituted by divalent transition metal cations (Fe2+, Mn2+, etc. ), rare earth element ions (La3+, Ce3+, Gd3+, etc. ), or alkali and alkaline earth metal ions (K+, Na+, Sr2+, etc. ), resulting in various colors such as yellow, green, blue, purple, and photochromic effects. Green fluor-hydroxyapatite turns into a bright blue hue after heat treatment at 650~700 ℃, possessing high ornamental and commercial value. The results show that the samplesare mainly composed of Ca2+ and P5+ with a small amount of N5+, which is mainly substituted [PO3] in the form of —NH2, indicating biochemical behavior correlated in formation. UV-Vis results show that the sample's color is mainly caused by the isomorphic substitution of Ce3+ and Mn2+ in the distorted octahedron of Ca2+, forming a blue-green transmission window. The defects in the green apatite sample before heat treatment are mainly ionic centers F--O--F-, combining with O2--O--V and OH--O--V, caused by the isomorphic substitution of Ce3+ at the CaⅠ site. After heat treatment, the recovery of hole centers is mainly manifested as ionic center defects caused by Mn2+ substituting the distorted octahedron at CaⅠ and CaⅡ sites. Combining EPR spectra shows that Mn2+ can be distributed at both CaⅠ and CaⅡ positions, causing a blue hue in the samples; otherwise, when Ce3+ is highly substituted at the CaⅠ position, the samples turn green.

Keyword: Green apatite; Heat treatment; Color mechanism; EPR spectra
引言

磷灰石标准化学式[Ca10(PO4)6 (F, Cl, OH)2], 空间群为P63/m, 属六方晶系钙磷酸盐矿物, 晶胞参数为a=9.37 Å , c=6.88 Å [1]。 磷灰石中附加阴离子(Cl-、 F-和OH-)与Ca2+相连形成畸变八面体结构[图1(a)], 该结构可直接反映结晶和演化温压条件, 也会导致不同磷灰石间物理性质的差异。 根据磷灰石的成分差异, 将其分为氟磷灰石Ca10(PO4)6F2、 氯磷灰石Ca10(PO4)6Cl2及羟基磷灰石Ca10(PO4)6(OH)2亚种[2, 3, 4]。 [PO4]3-与四个Ca2+相连形成稳定的磷钙桥单元, 呈正四面体精细对称结构, 构成磷灰石基本晶体格架[图1(a, b)], 可部分被 SiO44-SO42-CO32-等络阴离子团替代[4, 5]。 Ca2+在晶格中存在CaⅠ 和CaⅡ 两种占位形式, 可被Fe2+、 Mn2+等二价过渡族金属离子, La3+、 Ce3+、 Gd3+等稀土元素离子或K+、 Na+、 Sr2+等碱/碱土金属离子类质同象替代[5, 6, 7]。 CaⅠ 被九个不同键长的氧原子包围, 六个氧原子(O1和O2)与较短的键相连呈畸变三角棱柱状; 另外三个氧原子(O3)则与较长的键相连, 点群对称性为c3[图1(c)]。 CaⅡ 位于氟离子链共面边缘, 由六个氧原子和一个氟或氯原子围绕形成畸变的八面体结构, 点群对称性为m[图1(d)]。

图1 氟磷灰石的六方晶格晶体结构与元素占位分布图(根据Piper and Prener, 1972; Chen et al., 2002a修改)[4, 11]Fig.1 The hexagonal lattice crystal structure of fluorapatite with elemental occupancies (Modified according to Piper and Prener, 1972; and Chen et al., 2002a)

磷灰石既是一种应用广泛的无机非金属材料, 也是重要的地质指示矿物, 常用于生物活性和相容性材料、 半导体材料、 发光材料、 矿床及围岩年代标定, 以及岩浆及热液活动的地球化学特征示踪等[3, 5, 8, 9]。 高透明度和颜色艳丽者是常见的宝石品种, 有无色、 黄色、 绿色、 蓝色、 紫色等, 可见变色效应, 主要产于缅甸、 斯里兰卡、 莫桑比克、 马达加斯加、 加拿大、 巴西和中国等[10, 11, 12]。 磷灰石含较高的稀土元素(REE), Ce3+和Nd3+在橙红色和紫外区有较强的吸收, 而使磷灰石呈现绿色、 蓝色或紫色[7, 11]。 Mn2+、 Fe2+和Cu2+等离子替代畸变八面体中的CaⅠ 和CaⅡ 而形成黄色、 绿色、 蓝色和紫色等[5, 13]

绿黄色磷灰石在氧化或还原气氛中加热至650~700 ℃转变为更受欢迎蓝色品种, 使其价值增大, 但其致色与变色机理存在争议: Gilinskaya和Mashkovtsev认为绿色和蓝色分别由复合缺陷 SiO3--(Ce3++K+)和电子心 SO3-所致[14]; Hughes和Rakovan认为蓝色与 MnO43-有关[15]; Zhang等认为绿色的形成与Ce3+引起的可见光区补色有关[11]

本文选取绿黄色磷灰石进行400~800 ℃热处理实验, 结合电子探针(EPMA)、 傅里叶变换中红外光谱(Mid FTIR)、 紫外-可见光光谱(UV-Vis)和电子顺磁波谱(EPR), 旨在探究其致色与热处理变色机理。

1 实验部分
1.1 常规宝石学特征

本文选取一整块天然绿色磷灰石单晶体, 按照平行于c轴方向切割为9块厚度约2 mm的样品, 依次编号为Ap-1至Ap-9(图2)。 天然绿色磷灰石样品具玻璃光泽、 参差状断口, 长短波紫外线下均呈荧光惰性。 样品相对密度为3.18~3.21, 折射率平均为1.634~1.637, 双折射率约为0.003, 绿黄至绿色多色性, 符合标准氟羟基磷灰石宝石学特征。

图2 热处理前后绿色磷灰石的颜色对比Fig.2 The color comparison of green apatite before and after heat treatments

1.2 热处理

采用马弗炉(产品型号为MF-1100C-S, 最高温度1 100 ℃, 工作电压220 V, 额定功率1 kW, 控温精度± 1 ℃, 加热元件: 电阻丝), 弧形刚玉坩埚(容量20 mL)和SiO2粉末进行磷灰石样品在空气气氛下的热处理实验。 首先将样品放入刚玉坩埚中, 充填SiO2粉末, 从室温加热至目标温度后, 保温时间1 h, 恒温结束后缓慢降温冷却。 进行还原气氛下的热处理实验时, 刚玉坩埚中填充碳粉。 从400 ℃开始, 直至温度达到800 ℃, 分别对第一块至第九块样品间隔50 ℃并恒温1 h进行热处理。

1.3 测试方法

采用爱色丽Ci760色差仪对热处理前后样品的颜色进行测量, 记录Lab色空间中的L* a* b* 值, 并计算颜色变化Δ E。 采用日本岛津EPMA-1600型电子探针进行热处理前后磷灰石样品的组分能谱测试, 实验条件: 束流30 nA, 束斑直径10 μm, 加速电压15 kV, 元素峰位测试时间为10~60 s, 于华南理工大学分析测试中心完成。 采用布鲁克TENSOR-27型傅里叶变换红外光谱仪对样品分别进行反射法和透射法测试, 扫描范围: 2 000~400和4 000~2 000 cm-1, 分辨率: 4 cm-1, 扫描次数: 32。 采用广州标旗GEM-3000型紫外-可见光分光光度计对样品进行透射法测试, 检测范围: 200~1 100 nm, 平滑宽度3, 电压220 V, 功率250 W。 红外光谱和紫外光谱测试, 均在广州城市理工学院珠宝学院实验室完成。 采用布鲁克EMXplus-10/12型电子顺磁共振波谱仪进行实验测试, 于广东工业大学分析检测中心完成, 测试条件为室温, 磁场范围(4 100± 2 000) G, 扫描时间为120 s, 工作频率为9 838 Mc· s-1

2 结果与讨论
2.1 热处理结果

随着目标温度的上升, 样品绿色逐渐变浅, 而蓝色逐渐加深[图2(a)]。 未经热处理样品的平均CIE Lab数值为: L* =45.2, a* =-15.1, b* =18.3, 根据CIE 1976公式计算得出Δ Eab=2.6± 1.3, 显示样品颜色在肉眼观察下为均匀绿色[图2(b)]。 当热处理温度达到600 ℃时, 样品开始呈现浅蓝色(Δ Eab=33.4); 升至650~750 ℃呈现鲜艳的蓝色(Δ Eab范围为35.7~38.1); 但超到750 ℃后样品裂隙明显, 并在达到800 ℃后开始褪色(Δ Eab=31.9)[图2(c)]。 因此, 综合样品的色差值与美观性认为, 650~700 ℃是绿色磷灰石热处理为蓝色的最佳温度区间。

2.2 化学成分分析

选取热处理效果较好的样品AP-6和AP-7进行热处理前后的EPMA测试(表1), 采用氧原子法(Atoms pfu=26)进行晶体化学式计算, 其中OH-化学计量系数根据Kectham提供的计算方法进行校正[1]。 结果显示磷灰石样品主要由Ca2+和P5+组成, 及少量的N5+、 Si4+、 F-和OH-(质量分数< 3%), 应定名为氟羟磷灰石。 N5+的出现表明样品为生物沉积型成因, 在磷灰石中N元素常以 NO3-离子和— NH2基团两种形式存在, 需要结合光谱学特征对其赋存状态进一步分析。 F-和OH-是样品的主要挥发性成分, 热处理前后二者化学分配系数比值(XF/XOH)变化范围分别为0.28~0.45和~0.37, 表现出逐渐均一化的现象, 推测与热处理导致样品结晶度提高和挥发分饱和有关。

表1 绿色磷灰石热处理前后的电子探针化学成分特征 Table 1 Chemical compositional characteristics of green apatite with EPMA before and after heat treatment
2.3 中红外反射和透过光谱(Mid FT-IR)特征

样品Ap-6和Ap-7中红外光谱反射和透过特征基本一致。 指纹区(400~1 500 cm-1)反射峰位于472、 518、 575、 605、 962、 1 057和1 103cm-1附近, 由 PO43-振动引起。 其中472 cm-1归属于[PO3]四面体结构的二级倍频弯曲振动γ 2(P— O— P); 518、 575和605 cm-1归属于[PO3]的四级倍频不对称弯曲振动γ 4(P— O— P); 962 cm-1归属于[PO3]对称伸缩振动ν s(P— O— P); 1 057 cm-1归属于[PO3] 不对称伸缩振动ν as(P— O— P); 以及1 103 cm-1处的 HPO42-三级倍频联合吸收[图3(a)], 未见 CO32-振动特征[11, 16, 17]。 官能团区(1 500~4 000 cm-1)透过峰位于1 645、 1 728、 1 759、 1 784、 1 859、 1 930、 2 384、 2 488、 2 513、 2 879、 2 956、 3 038、 3 097、 3 379和3 558 cm-1附近。 其中1 600~2 000 cm-1区间透过峰归属于— NH2胶原类有机物面内弯曲振动γ (N— H)[16, 17]; 2 384 cm-1归属于OH-的不对称伸缩振动ν as(O— H), 随着热处理温度的升高透过率增加; 2 400~3 400 cm-1归属于吸附水(H2O)键振动; 3 558 cm-1附近透过峰归属于OH-的对称伸缩振动ν s(O— H) [图3(b)][6, 8, 18]

图3 热处理绿色磷灰石的中红外反射和透射光谱特征Fig.3 The Mid FTIR characteristics of green and blue apatite before and after heat treatment

与初始样品相比, 热处理后样品的1 057 cm-1ν as(P— O— P)和1 103 cm-1 HPO42-三级倍频吸收增强, 表明与之相连的Ca2+位类质同象替代行为随着温度升高愈加频繁, 推测与样品的蓝色呈正相关联系[11]。 官能团区— NH2胶原类有机物γ (N— H)的存在导致样品结晶度降低以及吸附水含量的增加, 进一步证明其成矿机制为生物沉积型[17]。 2 400~3 400 cm-1范围的透过峰随温度升高而透过率增加, 说明热处理将导致吸附水首先挥发, 可能对样品色调及透明度产生影响。

2.4紫外-可见光吸收光谱(UV-Vis)特征

天然绿色磷灰石样品Ap-6和Ap-7垂直于c轴方向UV-Vis吸收光谱特征主要位于264~402和584~804 nm区域, 形成451~565 nm(以458 nm为中心)区域透过窗, 导致样品呈绿色[图4(a)]。

图4 热处理前后磷灰石样品的紫外-可见光吸收光谱特征对比Fig.4 The comparison of UV-Vis spectra between natural and heated apatite samples

333 nm弱吸收峰归属于Ce3+从4f— 5d轨道的平价允许跃迁, 402 nm强吸收峰归属于Ce3+和Ce4+之间的电荷转移[9]。 584和642 nm弱吸收带分别归属于Mn2+在3d轨道上的自旋禁戒跃迁4T1(4G)→ 6A1(6S)和自旋允许跃迁4T1(4G)→ 4T2(4G)[9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19]。 700~870 nm红外光区连续吸收窄带归属于REE离子电荷转移, 如Nd3+和Ce3+[11, 14]。 热处理后, 300~402 nm区域内出现部分弱吸收带缺失和蓝移, 推测与畸变八面体中Ce3+的解耦有关; 584~804 nm区域吸收带则向红区偏移约3 nm, 是由Mn2+类质同象浓度的增加引起的晶格收缩所致[6]。 二者联合作用形成442~563 nm(以503 nm为中心)透过窗, 导致样品呈现蓝色[图4(b)]。 新产生的322和366 nm弱吸收峰分别归属于Ce3+/Ce4+类质同象替代CaⅡ 位置所致的5d至4f轨道跃迁5D3/22F7/25D3/22 F52[9]

2.5 电子自旋共振波谱(EPR)特征

样品热处理前后EPR波谱特征集中于2 800~3 500 G段(图5)。 当晶轴垂直于磁场方向时(c⊥H), 样品Ap-6热处理前EPR波谱特征峰g值分别为2.228、 2.208、 2.184和2.163[图5(a)], 其中g=2.228处表征指示Mn2+间隔~90 G的六条等距超精细谱线, 具有ms: +1/2↔ -1/2和I=5/2的轴对称自旋系统, 处于CaⅠ 位置; g=2.208处单条窄线归属于离子心F--O--F-的折叠效应, 具S=1/2的轴对称性, 处于F位置; 相伴出现的单条窄线g=2.184和g=2.163均由Ce3+替代CaⅠ 所致的空位(V)引起, 分别归属于O2-O--V和OH--O-V[20, 21]。 样品经650 ℃热处理后EPR特征峰以g=2.240为主[图5(b)], 归属于Mn2+的等距六重超精细谱线。 当晶轴平行于磁场方向时(c∥H), 天然样品Ap-7显示出不同的EPR谱峰特征, g值分别为2.181和2.171[图5(c)], 其中g=2.181归属于离子心F--O--F-, 由三条相距约35 G的超精细结构谱线构成, 强度比为1∶ 2∶ 1, 分别为g=2.204、 g=2.183和g=2.157[21, 22]; g=2.171归属于F-缺失所致的空位OH--O-V[20]。 经700 ℃热处理后样品EPR特征峰分别由以g=2.157为中心六条等距超精细结构谱线(约96 G)和以g=2.153为中心五重双线结构谱线(强度比例约为1∶ 2)组成, 其中g=2.157归属于Mn2+类质同象替代CaⅡ 位畸变八面体所致的离子心, g=2.153则归属于与CaⅡ 位相连的离子心F--O--F-[9]

图5 热处理前后磷灰石样品的EPR波谱特征对比Fig.5 The comparison of EPR spectra between natural and heated apatite samples

本文中所测样品均为氟羟磷灰石, 热处理前后样品的XF/XOH变化范围分别为0.28~0.45和~0.37, 表现出逐渐均一化的现象, 推测与热处理导致样品结晶度提高和挥发分饱和有关。 进一步结合EPMA和Mid FT-IR结果可知, 热处理前后N5+主要以— NH2形式少量替代[PO3], 说明样品形成过程与生物化学作用密切相关。 同时, — NH2胶原类有机物的存在将导致样品结晶度降低以及吸附水含量的增加, 但更有利于过渡族和稀土金属离子的捕获, 为样品致色提供稳定的化学物质来源。

热处理前后UV-Vis和EPR谱学特征联合显示, Ce3+和Mn2+类质同象替代CaⅠ 和CaⅡ 位置分别产生264~402和584~804 nm组合强吸收宽带, 并形成蓝绿区透过窗。 热处理前绿色磷灰石样品中c⊥H方向的缺陷主要为离子心F--O--F-以及由Ce3+类质同象替代CaⅠ 位所致的空穴心O2-O--V和OH--O-V, 热处理后由于空穴心的恢复则主要表现为Mn2+替代CaⅠ 位畸变八面体引起的离子心缺陷; 在c∥H方向的缺陷主要表现为强度比1∶ 2∶ 1的延展离子心F--O--F-; 热处理后则变为Mn2+和F-频繁替代CaⅡ 位畸变八面体引起的离子心缺陷。 上述现象说明, 所测磷灰石样品中的Ce3+一般占据CaⅠ 位置, 通过空穴心的键长弥补, 在333和402 nm分别以平价允许跃迁以及电荷转移的形式产生紫区吸收。 随着Ce3+引起的空穴心减少, 畸变八面体中的Ce3+解耦导致300~402 nm区域内蓝移现象, 变为紫外区吸收为主。 而Mn2+在CaⅠ 和CaⅡ 位置均有分布, 通过类质同象替代的电子自旋禁戒跃迁和自旋允许跃迁, 无论是热处理前后在584和642 nm处均能形成稳定橙红区宽吸收。 二者相互叠加作用在蓝绿区形成透过窗, 并导致热处理前呈现绿色而热处理后呈现蓝色。

3 结论

(1)磷灰石样品主要由Ca2+和P5+组成, 及少量的N5+、 Si4+、 F-和OH-(wt%< 3%)应定名为氟羟基磷灰石。

(2)绿色磷灰石样品在还原性气氛中经热处理转变为鲜艳的蓝色, 其最佳热处理工艺条件为650~700 ℃且恒温1 h。

(3)红外光谱官能团区(1 600~2 000) cm-1范围透过峰归属于— NH2胶原类有机物面内弯曲振动γ (N— H), 同时, 在电子探针结果也检测出N5+的存在, 二者可联合作为判别生物沉积型成因磷灰石的重要依据。

(4)磷灰石样品的绿色和蓝色均由Ce3+和Mn2+类质同象替代CaⅠ 和CaⅡ 位置形成的蓝绿区透过窗所致。 Mn2+是样品呈现蓝色的直接原因, Ce3+在CaⅠ 位置类质同象替代的程度较高时, 样品颜色以绿色为主, 反之呈蓝色。

参考文献
[1] Ketcham R A. American Mineralogist, 2015, 100(7): 1620. [本文引用:2]
[2] XING Kai, SHU Qi-hai(邢凯, 舒启海). Mineral Deposits(矿床地质), 2021, 40(2): 189. [本文引用:1]
[3] Stock M J, Humphreys M C S, Smith V C, et al. Journal of Petrology, 2018, 59(12): 2463. [本文引用:2]
[4] Piper W W, Prener J S. Physical Review B, 1972, 6(7): 2547. [本文引用:2]
[5] Iskhakova K, Murzakhanov F, Mamin G, et al. IOP Conf. Series: Earth and Environmental Science, 2018, 155: 012002. [本文引用:4]
[6] Kamal H, Hezma A M. Physical Science International Journal, 2015, 7(3): 137. [本文引用:3]
[7] Pan Y M, Dong P, Chen N. Geochimica et Cosmochimica Acta, 2003, 67(10): 1889. [本文引用:2]
[8] Hammerli J, Hermann J, Tollan P, et al. Contributions to Mineralogy and Petrology, 2021, 176(12): 105. [本文引用:2]
[9] Zhou Q, Dolgov L, Srivastava A M, et al. Journal of Materials Chemistry C, 2018, 6(11): 2652. [本文引用:5]
[10] CHEN Chao-yang, HUANG Wei-zhi, SHAO Tian, et al(陈超洋, 黄伟志, 邵天, ). Spectroscopy and Spectral Analysis(光谱学与光谱分析), 2021, 41(5): 1483. [本文引用:2]
[11] Zhang Z Y, Xu B, Yuan P Y, et al. Crystals, 2022, 12(8): 1067. [本文引用:7]
[12] Chindudsadeegul P, Jamkratoke M. Spectrochimica Acta Part A: Molecular and Biomolecular Spectroscopy, 2018, 204: 276. [本文引用:2]
[13] Drouet C. Journal of Chemical Thermodynamics, 2019, 136: 182. [本文引用:2]
[14] Gilinskaya L G, Mashkovtsev R I. Journal of Structural Chemistry, 1995, 36(1): 76. [本文引用:3]
[15] Hughes J M, Rakovan J F. Elements, 2015, 11(3): 165. [本文引用:2]
[16] Tao J H. Methods in Enzymology, 2013, 532: 533. [本文引用:3]
[17] Grunenwald A, Keyser C, Sautereau A M, et al. Journal of Archaeological Science, 2014, 49: 134. [本文引用:4]
[18] TU Jiang-lei, GUO Fu-qiang, Chun-chun, et al(涂姜磊, 郭富强, 吕春春, ). Journal of Biomedical Engineering(生物医学工程学杂志), 2011, 28(1): 99. [本文引用:2]
[19] Zhang C T, Chen C Y, Li Z B, et al. Minerals, 2023, 13(9): 1139. [本文引用:1]
[20] Pan Y M, Chen N, Weil J A. American Mineralogist, 2002, 87, 10: 1333. [本文引用:2]
[21] CAI Xiu-cheng, FU Yu-de, TANG Rong-bing(蔡秀成, 富毓德, 唐荣炳). Acta Mineralogica Sinica(矿物学报), 1981, 1(2): 77. [本文引用:2]
[22] Mayer I, Cuisinier F J G, Popov I, et al. European Journal of Inorganic Chemistry, 2006, 7: 1460. [本文引用:1]