1. 引言
近年来,核能应用领域不断拓展,辐射探测技术在科学研究和日常生活中发挥着越来越重要的作用[1]。核辐射探测技术应用广泛,在核医学、核电站安全监测,环境放射性探测和空间粒子物理学等诸多领域都展现出显著的应用价值[2]。根据核心材料和工作原理的不同,辐射探测器可以分为气体探测器、闪烁体探测器和半导体探测器。在这些类型中,闪烁体探测器因其独特的性能优势成为应用最为广泛的一类探测器,它可以将高能射线(如X射线和γ射线)转化为紫外至红外波长的可见光。这些光信号经光电倍增管(PMT)转换成电信号或脉冲信号,从而实现有效的射线检测。
闪烁体材料种类丰富多样,是闪烁体探测器的关键组成部分。其研发历史可以追溯到1903年William Crookes首次在实验中观察到高能α粒子吸收闪烁现象[3]。传统闪烁体大致可分为两类:无机闪烁体和有机闪烁体,二者区别主要在于其化学成分和结构[4] [5]。无机闪烁体由无机化合物构成,通常包含一种或多种金属离子。典型的无机闪烁体有NaI:Tl、CsI:Tl和Bi4Ge3O12(BGO) [6]-[8]。有机闪烁体是由有机分子构成的闪烁体,他们通常由含有芳香环和一些功能基团的化合物组成。典型的有机闪烁体包括有机晶体和塑料及液体闪烁体[9]。随着技术的发展,对闪烁体的要求也越来越高,包括大的有效原子序数、更高的光输出、更快的衰减和更高的能量分辨率。传统的闪烁体材料越来越难以满足所有要求,因此开发高性能的新型材料对于闪烁体探测技术的发展至关重要。
金属卤化物钙钛矿材料由于其大的有效原子质量、低成本的溶液制造工艺和优异的光电性能,因此在X射线和γ射线的直接和间接探测方面取得了一系列突破[10]-[12]。作为一种新发现的用于高能射线检测的闪烁体,金属卤化物钙钛矿在溶液加工制造、可见光谱可调发光和快速衰减时间显示出比传统闪烁体更强的优势。迄今为止,各类有机–无机和全无机金属卤化物钙钛矿已被广泛应用于闪烁体探测领域[13]。然而,由于其低激子结合能,大多数有机–无机钙钛矿材料在室温下表现出较差的稳定性和较低的光输出。相较之下,全无机钙钛矿被认为具有优异的环境稳定性,能够满足实际应用的需求。在全无机钙钛矿闪烁体中,铜基卤化物因其丰富性、低成本和无毒性,而受到广泛关注。CsCu2I3是铜基钙钛矿家族中一个相对稳定但研究较少的成员,其独特的物理和光学特性使其在高能射线探测领域展现出潜在的应用价值。其具备一维结构,能够实现单向和高效的电荷载流子传输,从而降低非辐射复合现象,这一特性使其在X射线成像领域具有独特的优势[14]。此外,CsCu2I3还表现出更快的光衰减特性,使其特别适合计算机断层扫描(CT)的X射线成像[15]。近年来,在X射线探测成像领域,CsCu2I3的应用取得了显著进展并展现出优异的性能。
作为高能辐射探测的闪烁体探测器,CsCu2I3钙钛矿闪烁体不可避免会面临辐照损伤的问题。在宏观上,材料在受到辐照后通常会发生体积膨胀、表面起泡和机械性能脆化等行为的转变。从微观的角度来看,在高能射线的辐照作用下,CsCu2I3闪烁体中会产生各种辐射缺陷,如体缺陷和点缺陷。体缺陷主要指晶界和错位,而点缺陷主要包括原子空位、原子间隙、原子反位替代等。这些缺陷对材料光电性质有着重要的影响。不幸的是,到目前为止,关于辐照缺陷是如何影响CsCu2I3闪烁体发光性能的理论计算研究却相对匮乏。为了揭示辐照缺陷影响闪烁体发光性能的深层物理机制,本文采用第一性原理计算方法,重点研究了Cs3Cu2I5晶体在不同空位缺陷条件下的缺陷形成能、电子结构和光学性质。本研究结果表明,CsCu2I3中I空位缺陷会在禁带中间引入深层缺陷能级,产生非辐射复合中心,降低闪烁体的发光性能。此外I空位缺陷还会提高闪烁体对可见光的吸收系数,增强了对发射光的自吸收能力,降低闪烁体探测器的探测效率。这些发现为优化CsCu2I3闪烁体提供了重要的理论依据。
2. 计算模型与方法
2.1. 模型构建
铜基卤化物钙钛矿闪烁体CsCu2I3属于正交晶系,所属空间群为Cmcm。在CsCu2I3的晶胞结构中,Cs离子占据八面体间隙位置,而[Cu2I3]−基团形成扭曲的八面体结构,整体呈现一维钙钛矿特征,如图1所示。本研究采用的是1 × 1 × 2的CsCu2I3的超胞模型,考虑了包括一个Cs位缺陷(VCs)、一个Cu位缺陷(VCu)以及两个I位缺陷(VI1、VI2)在内的不同元素的空位缺陷,如图2所示。
2.2. 计算方法
本研究的理论计算是在基于密度泛函理论框架的VASP (Vienna Ab initio Simulation Package)软件平台上完成[16] [17]。在计算参数的设置方面,采用投影增强赝势(PAW)方法处理核电子与价电子间的相互作用[18] [19],选取Perdew-Burke-Ernzerhof (PBE)形式的广义梯度近似(GGA)作为交换关联泛函[20]。为确保模拟计算的精确度,平面波基组的截断能设置为500 eV,经过测试验证该能量截断能够充分满足CsCu2I3闪烁体材料的计算要求。在结构优化过程中,我们采用全原子弛豫策略,设置Hellmann-Feynman力的收敛阈值在0.01 eV/Å,电子自洽迭代的能量收敛标准为1.0 × 10−5电子伏特。为准确描述布里渊区积分,在几何优化和电子自洽计算中均采用3 × 3 × 3的k点网格进行采样,确保计算结果的精确性和可靠性[21]。
Figure 1. The crystal structure diagram of CsCu2I3
图1. CsCu2I3的晶体结构示意图
Figure 2. The structure of defect system in CsCu2I3 scintillator: (a) VCs、(b) VCu、(c) VI1、(d) VI2, the red dashed circles indicate the positions of the vacancy defects
图2. CsCu2I3闪烁体缺陷体系的晶胞结构模型:(a) VCs、(b) VCu、(c) VI1、(d) VI2,其中红色虚线的圆圈表示空位缺陷的位置
3. 结果与讨论
3.1. 缺陷形成能计算
在缺陷形成能的计算中,本文采取了一种混合算法,这一方法通过特殊的k点取样方法计算体系种能量,再利用Γ点计算方法精确处理缺陷能级。缺陷α在电荷态q中的缺陷形成能Ef的计算公式如下:
(1)
式中Etotal和Eα分别代表不含缺陷和含α缺陷体系的总能量,ni表示原子数量的变化(ni > 0表示原子移除,ni < 0表示原子添加),μi是相应元素的化学势。EV表征价带顶能级的位置,EF为相对于价带顶位置的费米能级,其数值可在带隙范围内变动。
在闪烁体的生长过程中,不可避免地会产生各种各样的缺陷。这些缺陷可能会重新分布光生载流子改变其浓度分布,还可能会形成非辐射复合中心阻止电荷转移的过程。在热力学平衡状态下生长的Cs3Cu2I5晶体应该满足以下化学计量关系:
(2)
其中ΔμCs、ΔμCu和ΔμI被定义为从对应元素的相态转变到生长环境中所引起的化学势的变化。ΔHf(CsCu2I3)表示每化学式单位CsCu2I3的形成焓。为了避免Cs和Cu的生成以及I2的损失,应对Cs、Cu和I元素的化学势进行相应的限制。
(3)
为避免出现其它相互竞争的二元和三元化合物,如CsCu2I3、CsI、CsI3、CsI4、CuI和CuI4,必须仔细考虑这些化合物的形成条件。相应的化学势和焓变应满足以下条件,以确保目标化合物能够优先生长。
(4)
(5)
(6)
(7)
(8)
(9)
满足方程式(3)至方程式(9)的Cu和I元素的化学势显示在图3中的绿色线段,表明CsCu2I3的合成条件相当严格,需要对化学计量比进行精确的调控,以防止其他相的竞争形成。在这一可行的化学势范围内,我们选择了两个具有代表性的点A和C,以研究CsCu2I3的空位缺陷。点A代表富铜贫碘条件(ΔμCs = −2.87 eV, ΔμCu = −0 eV, ΔμI = −0.41 eV),而点C代表贫铜富碘的条件(ΔμCs = −3.12 eV, ΔμCu = −0.22 eV, ΔμI = −0.18 eV)。
Figure 3. The range of chemical potential limits corresponding to the equilibrium growth conditions of CsCu2I3
图3. CsCu2I3平衡生长条件所对应的化学势极限范围
每个缺陷电荷态的形成能与费米能级的函数关系如图4所示,其中线段的斜率表示缺陷的电荷态,而纽结则指示跃迁能级。费米能级在不同电荷态之间的结合位置,反映了带隙中热力学跃迁能级的分布。从图中可以观察到,所有的空位缺陷的形成能受生长条件和费米能级的影响。
通过对比分析,我们观察到某些缺陷类型,如VCs和VCu,这些缺陷在化学势的A点和C点均表现出较低的形成能,表明它们在这些条件下是占主导地位。在比较富铜贫碘与富碘贫铜两种合成环境下的缺陷形成能时,我们观察到显著的差异。在富碘贫铜的条件下,I空位缺陷的形成能显著升高,而Cs和Cu空位缺陷的形成能则相应降低。这一趋势表明,在富碘贫铜的条件下,I空位缺陷的形成受到抑制,而Cs和Cu空位缺陷的形成则相对促进。相反,在富铜贫碘的条件下,Cs和Cu空位缺陷的形成被抑制,而I空位缺陷的形成得到增强。因此,通过精细调控CsCu2I3闪烁体的合成环境,我们能够有效地调控缺陷的形成,这对于优化材料的性能具有重要意义。
Figure 4. The relationship between the formation of different charge states of vacancy defects in CsCu2I3 and the Fermi energy level, with the valence band maximum (VBM) and conduction band minimum (CBM) set at 0.00 and 1.83 eV, respectively. Only the charge states with the lowest formation energies for each defect are presented
图4. CsCu2I3中空位缺陷不同电荷态的形成与费米能级的关系,VBM和CBM分别被设置为0.00和1.83 eV,每个缺陷只给出了形成能最低的电荷态
3.2. 电子结构分析
固体能带理论为理解凝聚态物质的物理性质提供了关键的理论基础。固体是由大量相互作用的粒子构成的多粒子体系,其微观结构的变化难以通过直接求解原子核与电子简单体系的薛定谔方程计算。针对晶体的多体问题,能带理论采用单电子近似方法,将复杂的多体问题简化为容易处理的单电子问题。该理论通过构建周期性势场中的电子本征态,成功实现了对晶体电子结构的有效计算。这种简化的理论模型不仅可以降低计算的复杂程度,还能更好的解释材料的导电性、光学特性等物理性质。
为深入研究CsCu2I3中空位缺陷对发光性能的影响,本文基于缺陷形成能的计算结果,选取了形成能较低的Cs、Cu和I三种空位缺陷系统进行电子结构的计算分析。图5展示了初始体系和在这三种不同空位缺陷体系下CsCu2I3闪烁体的能带结构。晶体材料的物理性质主要取决于费米能级附近的电子态分布,因此本研究特别关注−2 eV~4 eV这一能量区间。理论计算结果表明,Cs3Cu2I5闪烁体的导带底(CBM)和价带顶(VBM)均处于布里渊区的Г点处,展现出典型的直接带隙特征。通过模拟计算得到的CsCu2I3的原胞带隙宽度为1.83 eV,与实验测得的值[22] (3.5 eV)相差较大,这主要源于PBE近似的影响,PBE泛函会低估了禁带宽度,在本研究中,我们主要定性分析空位缺陷对于本征CsCu2I3的影响,着重考虑相对变化而非绝对值,所以带隙的低估并不会影响我们的结论。计算采用PBE泛函,其结果与文献报道的数据[23] (约1.7 eV)表现出良好的一致性。
通过对空位缺陷体系的能带结构的分析,发现在VCs和VCu体系中,费米能级向下移动进入价带内部,带隙值分别提升至2.01 eV、2.05 eV。并且这两种缺陷会在VBM附近引入浅层缺陷能级,这些浅层能级能够拓宽热弛豫和辐射复合的路径,增强闪烁体的发光特性。在VI体系中,能带整体向低能量方向偏移,禁带宽度增至2.07 eV,同时在费米能级附近产生了新的能级,该能级距离VBM和CBM较远,通常被称为深层缺陷能级。与浅缺陷能级不同,由I空位产生的深缺陷能级可以“捕获”电子–空穴对,否则这些电子–空穴对将被发光中心有效捕获。这种消耗导致载流子从深缺陷能级逃逸的时间延长,从而影响发光中心的复合过程和发光寿命。此外,深缺陷能级可以作为非辐射复合中心,导致在这些中心大量捕获电子和空穴,从而降低载流子传输到辐射中心的效率,并最终导致闪烁体的发光性能和PLQY降低。浅缺陷能级和深缺陷能级还会分别增加载流子的散射和复合概率,显著改变载流子在闪烁体材料中的输运行为,进而延长闪烁体的响应时间。因此需要调控浅缺陷能级和深缺陷能级的分布,平衡载流子输运速度和发光性能。
Figure 5. The band structure of CsCu2I3: (a) initial system; (b) VCs; (c) VCu; (d) VI. The red lines indicate the positions of the defect levels
图5. CsCu2I3的能带结构:(a) 初始体系;(b) VCs;(c) VCu;(d) VI。红线表示缺陷能级位置
3.3. 态密度分析
通过投影态密度(PDOS)分析,我们能够更加深入研究Cs3Cu2I5闪烁体VBM和CBM的轨道组成特征。研究结果显示表明,VBM主要是由Cu的4d轨道和I的5p轨道构成的,而CBM的电子态主要来源于Cu的4s轨道和I的5p轨道的贡献。Cs元素的电子态在整个价带和导带边缘区域的分布密度均较低,仅在−9 eV价带内部存在峰值,这一现象表明Cs+在材料的本征发光过程中主要起结构支撑作用,而不直接参与辐射复合过程的电子跃迁。图6展示了三种空位缺陷体系的PDOS。在Cs和Cu空位缺陷体系的态密度分布中,态密度总体向右偏移,VBM穿过费米能级,因此被视为P型缺陷,主要由Cu的4d轨道及I的5p轨道贡献。在I空位缺陷体系中,态密度整体向低能量方向移动,并在费米能级附近引入深层缺陷态,因此被视为N型缺陷。这种深能级缺陷态的形成可以归因于空位缺陷附近I阴离子与Cu阳离子的重新杂化,主要来源于I的5p轨道和Cu的4s和3d轨道的相互作用,这种杂化模式与CBM的轨道组成具有高度的一致性。
Figure 6. The projected density of states of CsCu2I3: (a) initial system; (b) VCs; (c) VCu; (d) VI
图6. CsCu2I3的分波态密度:(a) 初始体系;(b) VCs;(c) VCu;(d) VI
3.4. 光学性质
此评估闪烁体材料优劣的关键参数包括较高的发光率、较低的光吸收系数和折射率等。为了建立固体电子结构与带间微观物理跃迁过程中的理论联系,我们引入了介电函数这一概念。
(10)
式中ε1(ω)为实部,ε2(ω)为虚部。其中ε2(ω)表示从占据态和未占据态之间的动量矩阵元素计算得到的虚部,其大小取决于电子在能带之间的跃迁:
(11)
式中,e代表电子的电荷,m代表自由电子的质量,ω代表入射光子的频率,
中m表示偶极矩阵,i和j分别代表初态和末态,fi表示波函数矢量k的第i个态的费米分布函数。
根据Kramers-Kronig色散关系及直接跃迁的概率,可以推导出ε1(ω)的表达式:
(12)
根据频率相关的介电函数,光学吸收系数α(ω)由下式计算:
(13)
图7展示了初始CsCu2I3晶体与三种空位缺陷体系(VCs、VCu和VI)的介电函数示意图。当入射光能量接近零时,介电函数实部纵坐标上的值反映了静介电常数。由图7(a)可知,初始体系与缺陷体系(VCs、VCu和VI)的静介电常数分别为4.55、8.55、10.87和4.69。所有缺陷体系的静介电常数均高于初始体系,这表明缺陷的存在增强了晶体的极化程度。其中,含Cu空位缺陷的CsCu2I3表现出最大的静介电常数,表明该体系具有最强的电荷束缚能力和最大的极化能力,从而使材料的光生电场变大,促进其内部光激发载流子的迁移,缩短探测响应时间。图7(b)为初始体系和缺陷体系的介电函数虚部,该参数直接反映了材料中电偶极子的能量损失过程,涉及电子的带间跃迁。计算结果表明,初始CsCu2I3的吸收边约在0.89 eV,首个吸收峰约在1.85 eV,这与CBM和VBM之间的电子跃迁相关。I空位缺陷体系在低能区显示出新的峰值,分别在1.23 eV和1.72 eV附近,这归因于碘空位缺陷内深层杂质电子跃迁至导带。Cs和Cu缺陷体系在可见光区域内未显现新的峰值,但在红外区(约0.1 eV)及1.5 eV附近出现了新的吸收峰,且可见光范围内吸收强度略有降低,这主要是Cs和Cu空位缺陷引起的浅层杂质能级间的跃迁所致。
Figure 7. The dielectric function of the CsCu2I3 primitive and defect system: (a) the real part and (b) the imaginary part
图7. CsCu2I3初始体系和缺陷体系的介电函数:(a) 实部和(b) 虚部
Figure 8. The variation of the absorption coefficient of CsCu2I3 scintillators in the primitive system and defect system with (a) photon energy and (b) wavelength
图8. CsCu2I3闪烁体在初始体系和缺陷体系中吸收系数随(a) 光子能量和(b) 波长的变化关系曲线
图8揭示了初始及缺陷CsCu2I3闪烁体的光学吸收系数。图8(a)数据表明,各缺陷体系的吸收带边缘均出现了红移现象,这与缺陷能级的形成和能带图中带隙增大的趋势相一致。由于光电倍增管的最佳工作波段位于可见光范围(380~740 nm),我们特别关注空位缺陷对CsCu2I3闪烁体在该波段内光吸收系数的影响。结果显示,含有Cs空位缺陷的CsCu2I3闪烁体在425~730 nm波段的光吸收能力低于本征CsCu2I3闪烁体。含有Cu空位缺陷的CsCu2I3闪烁体在410~730 nm波段的光吸收能力同样低于不含缺陷的CsCu2I3闪烁体。这表明Cs和Cu缺陷的存在显著降低了CsCu2I3闪烁体对可见光的光吸收。与之相反,含I空位缺陷的CsCu2I3闪烁体在可见光范围内的光吸收能力强于本征CsCu2I3闪烁体,这表明I缺陷提升了闪烁体材料对可见光的吸收能力。这一现象可能是由于Cu的4s、4p轨道及I的5p轨道之间的相互作用产生了深层缺陷能级,提高了电子跃迁效率并降低了所需能量,导致激发态电子数量增加,增强了对可见光的吸收能力。考虑到CsCu2I3闪烁体在高能射线激发下主要发射黄色光,发射波段在540~580 nm之间,自发形成的Cs和Cu空位缺陷会降低CsCu2I3闪烁体对发射光的自吸收,从而减少了向光电倍增管传输过程中的光损耗,这或许能够提高闪烁体探测器的探测效率。
4. 总结
本研究通过采用第一性原理计算方法,深入探讨了辐照损伤对CsCu2I3闪烁体探测性能的影响。从微观角度系统地分析了无铅铜基卤化物钙钛矿CsCu2I3中所有空位缺陷的性质,包括晶体结构、光电性质和缺陷容限,以及缺陷形成能。计算结果表明,Cs和Cu的空位缺陷是CsCu2I3中占主导地位的缺陷类型,且可以通过调控CsCu2I3闪烁体的合成环境,有效地控制相关缺陷的形成。I空位缺陷在带隙中引入了一个深能级缺陷态,这阻碍了载流子在发光中心的复合过程,降低了发光寿命。同时,这些缺陷产生了非辐射复合中心,导致电子和空穴通过非辐射方式(如热释放)复合,而非通过辐射(即发光)方式,从而降低了CsCu2I3闪烁体材料的发光效率和性能。与I空位缺陷不同,Cs和Cu空位缺陷引入的浅缺陷能级可以拓宽发光路径,并增加受热弛豫和辐射复合的速率,这对于提高闪烁体的探测效率至关重要。为避免吸收光谱与发射光谱重叠,并提高光电倍增管在可见光范围内的探测效率,降低闪烁体在可见光范围内的光学吸收系数是必要的。I缺陷的存在提高了CsCu2I3闪烁体在可见光范围内的光吸收系数,而Cs和Cu缺陷的存在则略微降低了对可见光的光学吸收。综上所述,通过抑制I缺陷的形成和适当促进Cs和Cu缺陷的形成,可以有效提高CsCu2I3闪烁体探测器的发光性能。这种策略还可以减少闪烁体内部发射光的自吸收,从而可能提高其探测效率。因此,碘化和退火处理以及金属元素掺杂有望改善材料的性能及其在光电器件中的应用表现。本研究不仅为理解CsCu2I3闪烁体在辐照环境下的行为提供了理论基础,也为设计和优化高性能闪烁体材料提供了重要的指导。
基金项目
本研究得到了国家重点研发计划项目(2022YFB1902700)、装备预研教育部联合基金项目(8091B042203)、国家自然科学基金项目(11875129)、强脉冲辐射环境模拟与效应国家重点实验室基金(SKLIPR1810)、等离子体物理全国重点实验室基金(6142A04240203)、核物理与核技术国家重点实验室开放课题(NPT2023KFY06)、中国铀业有限公司–东华理工大学核资源与环境国家重点实验室联合创新基金(2022NRE-LH-02)、中央高校基本科研业务费(2023JG001)等支持。
NOTES
*通讯作者。